## Mass mixing between QCD axions

Hai-Jun Li 1, ∗ and Yu-Feng Zhou 1, 2, 3, 4, †

1 Key Laboratory of Theoretical Physics, Institute of Theoretical Physics, Chinese Academy of Sciences, Beijing 100190, China 2 School of Physical Sciences, University of Chinese Academy of Sciences, Beijing 100049, China 3 School of Fundamental Physics and Mathematical Sciences, Hangzhou Institute for Advanced Study, UCAS, Hangzhou 310024, China

4 International Centre for Theoretical Physics Asia-Pacific, Beijing/Hangzhou, China (Dated: August 2, 2024)

We introduce a novel level crossing in the mass mixing between two QCD axions, one canonical QCD axion and one Z N QCD axion. The level crossing can take place at the QCD phase transition critical temperature or slightly before it, depending on the ratio of the axion decay constants ∼ 1 69. . The cosmological evolution of the mass eigenvalues in these two cases is similar, however, the transition of axion energy density is completely different. Finally, we estimate the relic density of the QCD axion dark matter. This level crossing may also have some cosmological implications.

Introduction. -The QCD axions are attractive candidates for the cold dark matter (DM). The canonical QCD axion was predicted by the Peccei-Quinn (PQ) mechanism [1, 2] to solve the strong CP problem in the Standard Model (SM) [3-8]. It obtains a tiny mass from the QCD non-perturbative effects [9, 10]. As the DM candidate, it can be non-thermally produced in the early Universe through the misalignment mechanism [11-13]. At high cosmic temperatures, the QCD axion is massless and obtains a non-zero mass during the QCD phase transition, the axion oscillation when its mass is comparable to the Hubble parameter explains the observed DM abundance. See e g . . Ref. [14] for a recent review. The reduced-mass QCD axion, Z N axion [15], can also solve the strong CP problem with N ⩾ 3 [16] and account for the DM through the trapped plus kinetic misalignment mechanism [17, 18]. Here N is an odd number, the N mirror worlds that are nonlinearly realized by the axion field under a Z N symmetry can coexist in Nature. Due to the suppressed non-perturbative effects on axion potential from the N degenerate QCD groups, the Z N axion mass is exponentially suppressed at the QCD phase transition. To avoid confusion, the following 'QCD axion' only stands for the canonical QCD axion.

The mass mixing in the multiple axions model [19] has attracted extensive attention in recent years. 1 Considering the non-zero mass mixing between these axions, the cosmological evolution process of the mass eigenvalues called the level crossing can take place and induce the adiabatic transition of the axion energy density, which is similar to the MSW effect [23-25] in the neutrino oscillations. The level crossing has some interesting cosmological implications, such as the modification of the axion relic density and isocurvature perturbations [26-28],

1 For instance, the mass mixing between the QCD axion and axionlike particle (ALP) [20] or sterile axion [21], and the mixing between the Z N axion and ALP [22], etc.

the domain walls formation [29], the gravitational waves emission and primordial black holes formation [30-32], and also the dark energy composition [33, 34], etc.

In this letter, we investigate the mass mixing between two QCD axions, one QCD axion and one Z N axion. We find that the level crossing can take place at the QCD phase transition critical temperature or slightly before it, depending on the relation of the axion decay constants. The conditions for level crossing to occur in these two cases are discussed in detail. The cosmological evolution of the mass eigenvalues in these two cases is similar, while the transition of energy density between the axions is different. Finally, we briefly estimate the relic density of the QCD axion and Z N axion DM through the misalignment mechanism. To our knowledge, this is the first time to investigate the effect of level crossing in the mass mixing with a multiple QCD axions model.

The model. -Here we consider a minimal multiple QCD axions model, one QCD axion ϕ and one Z N axion φ , with the low-energy effective Lagrangian

$$\mathbb { m } ^ { \infty } _ { \mathbb { m } ^ {$$

where m T a ( ) and m N ( T ) are the temperature-dependent QCD axion and Z N axion masses, respectively, f a and F a are the QCD axion and Z N axion decay constants. The QCD axion mass m T a ( ) can be described by

$$\intertext { e \vee e \in e \in e \in e \in e \$$

where m π and f π are the mass and decay constant of the pion, respectively, z ≡ m /m u d ≃ 0 48 is the ratio of . the up to down quark masses, T QCD ≃ 150 MeV is the

QCD phase transition critical temperature, and b ≃ 4 08 . is an index taken from the dilute instanton gas approximation. The first term in Eq. (2) corresponds to the zero-temperature QCD axion mass m a, 0 . While the Z N axion mass m N ( T ) can be described by

$$\begin{cases} \frac { m _ { \pi } f _ { \pi } } {$$

where γ ∈ (0 1) is a temperature parameter. , The first term in Eq. (3) also corresponds to the zero-temperature Z N axion mass m N , 0 , which is suddenly exponentially suppressed at T QCD due to the Z N symmetry. Then the equations of motion of ϕ and φ are given by

$$\begin{array} {$$

$$\ddot { \varphi } + 3 H \dot { \varphi } + \frac {$$

where H T ( ) is the Hubble parameter, and the dots represent derivatives with respect to the physical time t . Diagonalizing the mass mixing matrix

$$\mathbf M ^ { 2 } = \left ( \begin{array} { c c c c } \mu _ { a } ^ { 2 } ( T ) + \frac { m _ { \mathcal { N } } ^ { 2 } ( T ) F _ { a } ^ { 2 } } { \mathcal { N } ^ { 2 } f _ { a } ^ { 2 } } & \frac { m _ { \mathcal { N } } ^ { 2 } ( T ) F _ { a } } { \mathcal { N } f _ { a } } \\ \frac { m _ { \mathcal { N } } ^ { 2 } ( T ) F _ { a } } { \mathcal { N } f _ { a } } & m _ { \mathcal { N } } ^ { 2 } ( T ) \end{array} \right ) \,, \quad ( 6 ) \quad \text{as th} \, \\ \mu _ { a } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { a } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \nu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \mu _ { b } \, \nu _ { b } \, \mu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \mu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \mu _ { b } \, \nu _ { b } \, \mu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \mu _ { b } \, \nu _ { b } \, \nu _ { b } \, \mu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \mu _ { b } \, \mu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \mu _ { b } \, \mu _ { b } \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b. } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { b } \, \nu _ { c } \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \,\nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolimits \, \nolim$$

we derive the heavy ( a h ) and light ( a l ) mass eigenstates

$$\left ( \begin{array} { c } a _ { h } \\ a _ { l } \end{array} \right ) = \left ( \begin{array} { c c } \cos \alpha & \sin \alpha \\ - \sin \alpha & \cos \alpha \end{array} \right ) \left ( \begin{array} { c } \phi \\ \varphi \end{array} \right ) \,, \quad ( 7 )$$

with the corresponding mass eigenvalues m h,l ( T )

$$\text{with the corresponding mass eigenvalues } m _ { h, l } ( T ) \\ m _ { h, l } ^ { 2 } ( T ) & = \frac { 1 } { 2 } \left [ m _ { a } ^ { 2 } ( T ) + m _ { \mathcal { N } } ^ { 2 } ( T ) + \frac { m _ { \mathcal { N } } ^ { 2 } ( T ) F _ { a } ^ { 2 } } { \mathcal { N } ^ { 2 } f _ { a } ^ { 2 } } \right ] \quad \text{where} \quad \text{in or} \quad \text{zero-} \\ & \pm \frac { 1 } { 2 \mathcal { N } ^ { 2 } f _ { a } ^ { 2 } } \left [ - 4 \mathcal { N } ^ { 4 } m _ { a } ^ { 2 } ( T ) m _ { \mathcal { N } } ^ { 2 } ( T ) f _ { a } ^ { 4 } \quad \text{the z} \\ & + \left ( \left ( m _ { a } ^ { 2 } ( T ) + m _ { \mathcal { N } } ^ { 2 } ( T ) \right ) \mathcal { N } ^ { 2 } f _ { a } ^ { 2 } \\ & + m _ { \mathcal { N } } ^ { 2 } ( T ) F _ { a } ^ { 2 } \right ) ^ { 2 } \right ] ^ { 1 / 2 }, \\ \text{where} \, \alpha \, \text{is the mass minor and} \quad \text{do}$$

where α is the mass mixing angle

$$\cos ^ { 2 } \alpha = \frac { 1 } { 2 } \left ( 1 + \frac { m _ { \mathcal { N } } ^ { 2 } ( T ) - m _ { a } ^ { 2 } ( T ) - \frac { m _ { \mathcal { N } } ^ { 2 } ( T ) F _ { a } ^ { 2 } } { \mathcal { N } ^ { 2 } f _ { a } ^ { 2 } } } { m _ { h } ^ { 2 } ( T ) - m _ { l } ^ { 2 } ( T ) } \right ) \,. \ ( 9 ) \quad \text{which} \quad \text{ an illu}$$

FIG. 1. The illustration of level crossing in the Case I. The red and blue lines represent the mass eigenvalues m h ( T ) and m T l ( ), respectively. The black dashed line represents the Hubble parameter H T ( ). The gray dashed line represents the temperature T QCD . We set f a = 10 12 GeV, F a = 10 11 5 . GeV ( ⇒ ζ ≃ 0 32), . N = 17, and γ = 0 25. .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00267v1-figure-8b1ee2ee-a0a0-437b-95b4-a1425a549dc4.jpg)

The so-called axion level crossing is that when considering the interaction between two axion fields, there is a nontrivial mass mixing between their temperaturedependent mass eigenvalues m h,l ( T ) that as functions of the cosmic temperature T . It can take place when the difference of m 2 h ( T ) -m T 2 l ( ) gets a minimum value.

Case I: Level crossing at T QCD . -We first consider a case that the zero-temperature QCD axion mass m a, 0 is smaller than the second term in Eq. (3) (we can define it as the mass m N ,π )

$$\frac { m _ { \pi } f _ { \pi } } { f _ { a } } \frac { \sqrt { z } } { 1 + z } < \frac { m _ { \pi } f _ { \pi } } { F _ { a } } \sqrt { \frac { z } { 1 - z ^ { 2 } } } \,, \quad \quad ( 1 0 )$$

it can be characterized as

$$\zeta \equiv \frac { F _ { a } } { f _ { a } } < \sqrt { \frac { 1 + z } { 1 - z } } \simeq 1. 6 9 \,, \quad \quad \quad ( 1 1 )$$

where we have defined the ratio ζ . On the other hand, in order for the level crossing to occur in this case, the zero-temperature QCD axion mass should be larger than the zero-temperature Z N axion mass m N , 0 , we have

$$\frac { m _ { \pi } f _ { \pi } } { f _ { a } } \frac { \sqrt { z } } { 1 + z } > \frac { m _ { \pi } f _ { \pi } } { \sqrt { \pi } F _ { a } } \sqrt { \frac { 1 - z } { 1 + z } } \mathcal { N } ^ { 3 / 4 } z ^ { \mathcal { N } / 2 } \,, \quad ( 1 2 )$$

it can also be characterized as

$$\zeta & > \frac { 1 } { \sqrt { \pi } } \, \sqrt { \frac { \left ( 1 - z \right ) \left ( 1 + z \right ) ^ { 3 } } { z ^ { 2 } } } \mathcal { N } ^ { 3 / 4 } z ^ { \mathcal { N } / 2 } \\ & \simeq 1. 2 4 \times 0. 4 8 ^ { \mathcal { N } / 2 } \mathcal { N } ^ { 3 / 4 } \,,$$

which depends on the value of N . Then we show in Fig. 1 an illustration of level crossing in this case. We have set

FIG. 2. The conditions for level crossing to occur in the {N , ζ } plane. The red and blue shadow regions correspond to the Cases I and II, respectively. The gray line represents ζ = 1 69. . Note that N = 2 k +1, k ∈ N + .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00267v1-figure-b760dff9-154c-449d-a596-0bb76c0aa1c2.jpg)

f a = 10 12 GeV, F a = 10 11 5 . GeV, N = 17, and γ = 0 25. . The red and blue lines correspond to the temperaturedependent mass eigenvalues m h ( T ) and m T l ( ), respectively. Note that the level crossing occurs at the QCD phase transition critical temperature

$$T _ { \times } = T _ { \text{QCD} } \,. \quad \quad \ \ ( 1 4 ) \quad \text{ob}$$

Let us discuss the cosmological evolution of the mass eigenvalues in this case. At high temperatures, the heavy mass eigenvalue m h ( T ) corresponds to the Z N axion, while the light one m T l ( ) corresponds to the QCD axion. They will approach to each other at T QCD and then move away from each other. When T &lt; T QCD , the m h ( T ) represents the QCD axion, and the m T l ( ) represents the Z N axion. For a more intuitive description, they can be expressed as

$$m _ { h } ( T ) \Rightarrow \begin{cases} m _ { a } ( T ) \,, & T < T _ { \times } \\ m _ { \mathcal { N } } ( T ) \,, & T > T _ { \times } \end{cases} \quad \text{(15)}$$

$$m _ { l } ( T ) \Rightarrow \begin{cases} m _ { \mathcal { N } } ( T ) \,, & T < T _ { \times } \\ m _ { a } ( T ) \,. & T > T _ { \times } \end{cases} \quad ( 1 6 ) \quad \text{where} \quad \text{where}$$

The above discussion is focused on the Case I that level crossing occurs at T QCD .

To compare with another case discussed below, we show the condition for level crossing to occur in the Case I by using Eqs. (11) and (13). See Fig. 2 with the red shadow region in the {N , ζ } plane. The blank regions represent that no level crossing occurs.

Case II: Level crossing slightly before T QCD . -Next we consider another case that the zero-temperature QCD axion mass is larger than m N ,π , we have

$$\zeta > \sqrt { \frac { 1 + z } { 1 - z } } \simeq 1. 6 9 \,.$$

FIG. 3. Same as Fig. 1 but for level crossing in the Case II. The gray dashed line represents the level crossing temperature T × ≃ 176 3 MeV. . We set f a = 10 12 GeV, F a = 10 12 5 . GeV ( ⇒ ζ ≃ 3 16), . N = 13, and γ = 0 25. .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00267v1-figure-7e3d8a67-91bd-413f-9d0e-9ddfcffcd421.jpg)

We show an illustration of this case in Fig. 3. We have set f a = 10 12 GeV, F a = 10 12 5 . GeV, N = 13, and γ = 0 25. . Different from the Case I, here the level crossing will take place slightly before the temperature T QCD . By solving d m ( 2 h ( T ) -m T 2 l ( )) /dT = 0 at T QCD &lt; T &lt; T QCD /γ , we obtain the level crossing temperature

$$\ s s _ { \substack { r y \\ \text{n.} } } & & T _ { \times } = T _ { \text{QCD} } \left ( \sqrt { \frac { 1 - z } { 1 + z } } \zeta \right ) ^ { 1 / b } \left ( 1 - \frac { \zeta ^ { 2 } } { \mathcal { N } ^ { 2 } } \right ) ^ { - 1 / ( 2 b ) } \,, \ \ ( 1 8 )$$

which is slightly larger than T QCD . The QCD axion mass at T × is given by

$$m _ { a } ( T _ { \times } ) \simeq m _ { \mathcal { N }, \pi } \sqrt { 1 - \frac { \zeta ^ { 2 } } { \mathcal { N } ^ { 2 } } } \,. \quad \quad ( 1 9 )$$

Note that ζ &lt; N should be satisfied in Eq. (18). This level crossing will last for a parametric duration

$$\Delta T _ { \times } = \left | \frac { 1 } { \cos \alpha ( T _ { \times } ) } \frac { d \cos \alpha ( T _ { \times } ) } { d T } \right | ^ { - 1 }, \quad \ \ ( 2 0 )$$

where α T ( ) is the mass mixing angle given by Eq. (9). This is the Case II. The temperature-dependent behavior of the mass eigenvalues is similar to the Case I. However, the transition of energy density between them is completely different, which will be discussed later. In fact, the Case I can transform into the Case II, depending on the value of the ratio ζ ∼ 1 69. .

Then we show in Fig. 2 (the blue shadow region) the condition for level crossing to occur in the Case II by using Eq. (17) and ζ &lt; N . Additionally, we also note that at T QCD /γ the QCD axion mass should be smaller than the Z N axion mass

$$\mu _ { 0 } \, \mu _ { 0 } \, \nu _ { 0 } \, \mu _ { 0 } \, \mu _ { 0 } \, \nu _ { 0 } \, \mu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \mu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \mu _ { 0 } \, \nu _ { 0 } \, \mu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \mu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \mu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \mu _ { 0 } \, \mu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \mu _ { 0 } \, \mu _ { 0 } \, \nu _ { 0 } \, \mu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \mu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \mu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \mu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \mu _ { 0 } \, \nu _ { 0 } \, \mu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \mu _ { 0 } \, \nu _ { 0 } \, \mu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \mu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \mu _ { 0 } \, \nu _ { 0 } \, \mu _ { 0 } \, \mu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \mu _ { 0 } \, \nu _ { 0 } \, \mu _ { 0 } \, \mu _ { 0 } \, \nu _ { 0 } \, \mu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \mu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \mu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nolimits \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nolimits \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 1 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 2 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { - 1 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _ { 0 } \, \nu _$$

it can be characterized as

$$\zeta < \sqrt { \frac { 1 + z } { 1 - z } } \gamma ^ { - b } \simeq 1. 6 9 \, \gamma ^ { - b } \,, \quad \quad ( 2 2 )$$

which depends on the value of γ . Since for the small γ , the constraint of ζ &lt; N is more stringent than Eq. (22), we do not show it in Fig. 2.

Axion relic density. -Here we briefly estimate the relic density of the axion DM in the Cases I and II through the misalignment mechanism. The pre-inflationary scenario is considered, in which the PQ symmetry is spontaneously broken during inflation.

Wefirst discuss the Case I. As an example, we focus our attention on the Z N axion DM. We begin with the QCD axion field, which is frozen at high temperatures with an arbitrary initial misalignment angle θ 1 ,a , and starts to oscillate at T 1 ,a . The oscillation temperature T 1 is given by m T ( ) = 3 H T ( ) with the Hubble parameter

$$H ( T ) = \sqrt { \frac { \pi ^ { 2 } g _ { * } ( T ) } { 9 0 } } \frac { T ^ { 2 } } { m _ { \text{Pl} } } \,, \quad \quad ( 2 3 ) \quad \text{whi} \\ \text{erg}.$$

where m Pl is the reduced Planck mass, and g ∗ is the number of effective degrees of freedom of the energy density. In the following, the subscript ' x ' stands for the physical quantity at T x or corresponding to x . The initial energy density in the QCD axion field is

$$\rho _ { a, 1 } = \frac { 1 } { 2 } m _ { a, 1 } ^ { 2 } f _ { a } ^ { 2 } \theta _ { 1, a } ^ { 2 } \,.$$

At T QCD &lt; T &lt; T 1 ,a , the QCD axion energy density is adiabatic invariant with the comoving number N a ≡ ρ a /m a 3 a , where a is the scale factor. Then we have the QCD axion energy density just before T QCD as

$$\rho _ { a, \text{QCD} } = \frac { 1 } { 2 } m _ { a, 0 } m _ { a, 1 } f _ { a } ^ { 2 } \theta _ { 1, a } ^ { 2 } \left ( \frac { a _ { 1, a } } { a _ { \text{QCD} } } \right ) ^ { 3 } \,. \quad ( 2 5 ) \quad \text{N} \epsilon$$

In this case, the QCD axion is trapped around θ tr ≡ θ a ( T QCD ) ∼ π until T QCD with the initial axion velocity ˙ θ tr ≡ ˙ θ a ( T QCD ). At T = T QCD , the light mass eigenvalue m T l ( ) will comprise the Z N axion, the axion mass is suddenly suppressed and the true minimum ∼ 0 develops. Then we have the mean velocity of the Z N axion

$$\sqrt { \langle \dot { \theta } _ { \text{tr} } ^ { 2 } \rangle } = \frac { 1 } { \sqrt { 2 } \zeta } \sqrt { m _ { a, 0 } m _ { a, 1 } } \theta _ { 1, a } \left ( \frac { a _ { 1, a } } { a _ { \text{QCD} } } \right ) ^ { 3 / 2 } \,. \quad ( 2 6 )$$

Note that the Z N axion energy density at T QCD is nonadiabatic, which just after T QCD is given by

$$\rho _ { \mathcal { N }, \text{tr} } = \frac { 1 } { 2 } F _ { a } ^ { 2 } \dot { \theta } _ { \text{tr} } ^ { 2 } + 2 \frac { m _ { \mathcal { N } } ^ { 2 } F _ { a } ^ { 2 } } { \mathcal { N } ^ { 2 } } \,, \quad \quad ( 2 7 ) \quad \stackrel { \dots \dots } { \rho _ { a, \times } } } { \text{at $1$} }$$

where we consider a case that the axion has a large mean velocity. Then the Z N axion will start to oscillate at T 2 ,

FIG. 4. The three terms in Eq. (32) as functions of the temperature T . The black, red, and blue lines represent ∆ t × , 2 π/m l ( T ), and 2 π/ m ( h ( T ) -m T l ( )), respectively. The gray dashed line represents T × ≃ 176 3 MeV. .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00267v1-figure-5e448084-8dc2-4355-be90-b464c4b8b784.jpg)

which is defined as the temperature when the kinetic energy is equal to the barrier height. At T 2 &lt; T &lt; T QCD , the Z N axion energy density is conserved with the comoving PQ charge q kin = ˙ θ N a 3 , and we have the scale factor at T 2 as a 2 = ( N ˙ θ tr / (2 m N )) 1 / 3 a QCD . Then at T 0 &lt; T &lt; T 2 , the adiabatic approximation is valid, and we have the Z N axion energy density at present T 0 as

$$\rho _ { \mathcal { N }, 0 } = \frac { m _ { \mathcal { N }, 0 } F _ { a } ^ { 2 } \dot { \theta } _ { \text{tr} } } { \mathcal { N } } \left ( \frac { a _ { \text{QCD} } } { a _ { 0 } } \right ) ^ { 3 } \,. \quad \quad ( 2 8 )$$

Substituting the mean velocity, it reads

$$\stackrel { \ m o w _ { J } } { N _ { a } } \equiv \\ \text{e the} \quad \rho _ { \mathcal { N }, 0 } = C \frac { m _ { \mathcal { N }, 0 } \sqrt { m _ { a, 0 } m _ { a, 1 } } \theta _ { 1, a } F _ { a } ^ { 2 } } { \sqrt { 2 } \mathcal { N } \zeta } \left ( \frac { \sqrt { a _ { 1, a } a _ { Q C D } } } { a _ { 0 } } \right ) ^ { 3 } \,, \, ( 2 9 )$$

where C ≃ 2 is a constant.

Next we discuss the Case II, focusing on the QCD axion DM. We begin with the Z N axion field, it starts to oscillate at T 1 , N with the initial misalignment angle θ 1 , N , and the initial energy density is given by

$$\rho _ { \mathcal { N }, 1 } = \frac { 1 } { 2 } m _ { \mathcal { N }, 1 } ^ { 2 } F _ { a } ^ { 2 } \theta _ { 1, \mathcal { N } } ^ { 2 } \,.$$

At T × &lt; T &lt; T 1 , N , the Z N axion energy density is adiabatic invariant, and at T × we have

$$\rho _ { \mathcal { N }, \times } = \frac { 1 } { 2 } m _ { \mathcal { N }, \pi } m _ { \mathcal { N }, 1 } F _ { a } ^ { 2 } \theta _ { 1, \mathcal { N } } ^ { 2 } \left ( \frac { a _ { 1, \mathcal { N } } } { a _ { \times } } \right ) ^ { 3 } \,. \quad ( 3 1 )$$

At T = T × , the level crossing occurs, the heavy mass eigenvalue m h ( T ) will comprise the QCD axion and the energy density ρ N × , is transferred to the QCD axion ρ a, × . To ensure this axion energy transition is adiabatic at T × , we need

$$\Gamma _ { \ast } \gg \max \left [ \frac { 2 \pi } { m _ { l } ( T _ { \times } ) }, \frac { 2 \pi } { m _ { h } ( T _ { \times } ) - m _ { l } ( T _ { \times } ) } \right ] \,, \quad ( 3 2 )$$

where the time ∆ t × corresponding to the temperature in Eq. (20). We also check it in Fig. 4 with the parameter set shown in Fig. 3. Then at T 0 &lt; T &lt; T × , the adiabatic approximation is valid again, and we have the present QCD axion energy density

$$\rho _ { a, 0 } = \frac { 1 } { 2 } m _ { a, 0 } m _ { \mathcal { N }, 1 } F _ { a } ^ { 2 } \theta _ { 1, \mathcal { N } } ^ { 2 } \left ( \frac { a _ { 1, \mathcal { N } } } { a _ { 0 } } \right ) ^ { 3 } \,. \quad \left ( 3 3 \right ) \quad \begin{pmatrix} 3 \\ 3 \\ 5 \end{pmatrix} \,.$$

Note that we only show two examples above to estimate the relic density of the Z N axion and QCD axion DM in the Cases I and II, respectively, one can also calculate others in a similar way.

Conclusion. -In summary, we have investigated the mass mixing between the QCD axion and Z N axion, and introduced a novel level crossing. In the mixing, we find that the level crossing can take place at the QCD phase transition critical temperature T QCD (Case I) or slightly before T QCD (Case II), depending on the ratio of the axion decay constants ζ ∼ 1 69. . The conditions for level crossing to occur in the Cases I and II are discussed in detail. In the Case I, the zero-temperature QCD axion mass m a, 0 should be smaller than the mass m N ,π and larger than the zero-temperature Z N axion mass m N , 0 , which can be characterized as 1 24 . × 0 48 . N / 2 N 3 / 4 &lt; ζ &lt; 1 69. . While in the Case II, the zero-temperature QCD axion mass m a, 0 should be larger than m N ,π and the ratio should be ζ &lt; N , which can also be characterized as 1 69 . &lt; ζ &lt; N . Note that in this case we also have a relatively weak constraint ζ &lt; 1 69 . γ -b , which will be stringent for the large value of γ .

In these two cases, the cosmological evolution of the mass eigenvalues is similar, however, the transition of energy density between the axions is different. In the Case I, the axion energy transition at the level crossing is non-adiabatic ( T QCD ), while in the Case II it is considered to be adiabatic ( T × ). Finally, we briefly estimate the relic density of the axion DM through the misalignment mechanism, focusing on the Z N axion and QCD axion DM in the Cases I and II, respectively. The level crossing in this scenario may also have some other interesting cosmological implications.

Acknowledgments. -This work was supported by the CAS Project for Young Scientists in Basic Research YSBR-006, the National Key R&amp;D Program of China (Grant No. 2017YFA0402204), and the National Natural Science Foundation of China (Grants No. 11821505, No. 11825506, and No. 12047503).

∗ lihaijun@itp.ac.cn

† yfzhou@itp.ac.cn

[1] R.D. Peccei and Helen R. Quinn, 'CP Conservation in the Presence of Instantons,' Phys. Rev. Lett. 38 , 14401443 (1977).

- [2] R.D. Peccei and Helen R. Quinn, 'Constraints Imposed by CP Conservation in the Presence of Instantons,' Phys. Rev. D 16 , 1791-1797 (1977).
- [3] Steven Weinberg, 'A New Light Boson?' Phys. Rev. Lett. 40 , 223-226 (1978).
- [4] Frank Wilczek, 'Problem of Strong P and T Invariance in the Presence of Instantons,' Phys. Rev. Lett. 40 , 279282 (1978).
- [5] Jihn E. Kim, 'Weak Interaction Singlet and Strong CP Invariance,' Phys. Rev. Lett. 43 , 103 (1979).
- [6] Mikhail A. Shifman, A. I. Vainshtein, and Valentin I. Zakharov, 'Can Confinement Ensure Natural CP Invariance of Strong Interactions?' Nucl. Phys. B 166 , 493-506 (1980).
- [7] Michael Dine, Willy Fischler, and Mark Srednicki, 'A Simple Solution to the Strong CP Problem with a Harmless Axion,' Phys. Lett. B 104 , 199-202 (1981).
- [8] A. R. Zhitnitsky, 'On Possible Suppression of the Axion Hadron Interactions. (In Russian),' Sov. J. Nucl. Phys. 31 , 260 (1980).
- [9] Gerard 't Hooft, 'Symmetry Breaking Through BellJackiw Anomalies,' Phys. Rev. Lett. 37 , 8-11 (1976).
- [10] Gerard 't Hooft, 'Computation of the Quantum Effects Due to a Four-Dimensional Pseudoparticle,' Phys. Rev. D 14 , 3432-3450 (1976).
- [11] John Preskill, Mark B. Wise, and Frank Wilczek, 'Cosmology of the Invisible Axion,' Phys. Lett. B 120 , 127132 (1983).
- [12] L.F. Abbott and P. Sikivie, 'A Cosmological Bound on the Invisible Axion,' Phys. Lett. B 120 , 133-136 (1983).
- [13] Michael Dine and Willy Fischler, 'The Not So Harmless Axion,' Phys. Lett. B 120 , 137-141 (1983).
- [14] Ciaran A. J. O'Hare, 'Cosmology of axion dark matter,' PoS COSMICWISPers , 040 (2024), arXiv:2403.17697 [hep-ph].
- [15] Anson Hook, 'Solving the Hierarchy Problem Discretely,' Phys. Rev. Lett. 120 , 261802 (2018), arXiv:1802.10093 [hep-ph].
- [16] Luca Di Luzio, Belen Gavela, Pablo Quilez, and Andreas Ringwald, 'An even lighter QCD axion,' JHEP 05 , 184 (2021), arXiv:2102.00012 [hep-ph].
- [17] Luca Di Luzio, Belen Gavela, Pablo Quilez, and Andreas Ringwald, 'Dark matter from an even lighter QCD axion: trapped misalignment,' JCAP 10 , 001 (2021), arXiv:2102.01082 [hep-ph].
- [18] Raymond T. Co, Lawrence J. Hall, and Keisuke Harigaya, 'Axion Kinetic Misalignment Mechanism,' Phys. Rev. Lett. 124 , 251802 (2020), arXiv:1910.14152 [hepph].
- [19] Christopher T. Hill and Graham G. Ross, 'Models and New Phenomenological Implications of a Class of Pseudogoldstone Bosons,' Nucl. Phys. B 311 , 253-297 (1988).
- [20] Ryuji Daido, Naoya Kitajima, and Fuminobu Takahashi, 'Level crossing between the QCD axion and an axionlike particle,' Phys. Rev. D 93 , 075027 (2016), arXiv:1510.06675 [hep-ph].
- [21] David Cyncynates and Jedidiah O. Thompson, 'Heavy QCD axion dark matter from avoided level crossing,' Phys. Rev. D 108 , L091703 (2023), arXiv:2306.04678 [hep-ph].
- [22] Hai-Jun Li, Ying-Quan Peng, Wei Chao, and YuFeng Zhou, 'Light QCD axion dark matter from double level crossings,' Phys. Lett. B 849 , 138444 (2024), arXiv:2310.02126 [hep-ph].

- [23] L. Wolfenstein, 'Neutrino Oscillations in Matter,' Phys. Rev. D 17 , 2369-2374 (1978).
- [24] S. P. Mikheyev and A. Yu. Smirnov, 'Resonance Amplification of Oscillations in Matter and Spectroscopy of Solar Neutrinos,' Sov. J. Nucl. Phys. 42 , 913-917 (1985).
- [25] S. P. Mikheev and A. Yu. Smirnov, 'Resonant amplification of neutrino oscillations in matter and solar neutrino spectroscopy,' Nuovo Cim. C 9 , 17-26 (1986).
- [26] Naoya Kitajima and Fuminobu Takahashi, 'Resonant conversions of QCD axions into hidden axions and suppressed isocurvature perturbations,' JCAP 01 , 032 (2015), arXiv:1411.2011 [hep-ph].
- [27] Shu-Yu Ho, Ken'ichi Saikawa, and Fuminobu Takahashi, 'Enhanced photon coupling of ALP dark matter adiabatically converted from the QCD axion,' JCAP 10 , 042 (2018), arXiv:1806.09551 [hep-ph].
- [28] David Cyncynates, Tudor Giurgica-Tiron, Olivier Simon, and Jedidiah O. Thompson, 'Resonant nonlinear pairs in the axiverse and their late-time direct and astrophysical signatures,' Phys. Rev. D 105 , 055005 (2022), arXiv:2109.09755 [hep-ph].
- [29] Ryuji Daido, Naoya Kitajima, and Fuminobu Takahashi, 'Domain Wall Formation from Level Crossing
- in the Axiverse,' Phys. Rev. D 92 , 063512 (2015), arXiv:1505.07670 [hep-ph].
- [30] David Cyncynates, Olivier Simon, Jedidiah O. Thompson, and Zachary J. Weiner, 'Nonperturbative structure in coupled axion sectors and implications for direct detection,' Phys. Rev. D 106 , 083503 (2022), arXiv:2208.05501 [hep-ph].
- [31] Zhe Chen, Archil Kobakhidze, Ciaran A. J. O'Hare, Zachary S. C. Picker, and Giovanni Pierobon, 'Cosmology of the companion-axion model: dark matter, gravitational waves, and primordial black holes,' (2021), arXiv:2110.11014 [hep-ph].
- [32] Hai-Jun Li and Yu-Feng Zhou, 'Gravitational waves from axion domain walls in double level crossings,' (2024), arXiv:2401.09138 [hep-ph].
- [33] Kristjan M¨¨rsepp, Enrico Nardi, uu and Clemente Smarra, 'Can the QCD axion feed a dark energy component?' (2024), arXiv:2405.00090 [hep-ph].
- [34] Kristjan M¨¨rsepp, 'On the viability of the QCD axion uu - dark energy transition,' (2024) arXiv:2405.20478 [hepph].
Hai-Jun Li Yu-Feng Zhou
Mass mixing between QCD axions
We introduce a novel level crossing in the mass mixing between two QCD
axions, one canonical QCD axion and one $Z_{\mathcal N}$ QCD axion. The level
crossing can take place at the QCD phase transition critical temperature or
slightly before it, depending on the ratio of the axion decay constants
$\sim1.69$. The cosmological evolution of the mass eigenvalues in these two
cases is similar, however, the transition of axion energy density is completely
different. Finally, we estimate the relic density of the QCD axion dark matter.
This level crossing may also have some cosmological implications.