Spatially modulated vacua in a Lorentz-invariant scalar field theory

Spatial modulation has been studied for a long time in condensed matter, nuclear matter and quark matter, where the manifest Lorentz invariance is lost due to the finite density/temperature effects and so on. In this paper, spatially modulated vacua at zero temperature and zero density are studied in Lorentz invariant field theories. We first propose an adaptation of the Nambu–Goldstone theorem to higher derivative theories under the assumption of the absence of ghosts: when a global symmetry is spontaneously broken due to vacuum expectation values of space-time derivatives of fields, a Nambu–Goldstone (NG) boson appears without a canonical kinetic (quadratic derivative) term with a quartic derivative term in the modulated direction while a Higgs boson appears with a canonical kinetic term. We demonstrate this in a simple model allowing (meta)stable modulated vacuum of a phase modulation (Fulde–Ferrell state), where an NG mode associated with spontaneously broken translational and U(1) symmetries appears.


Introduction
Spatially modulated ground states were theoretically proposed in superconductors a half century ago [1][2][3], and such states are now called Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) states. More precisely, Fulde-Ferrell (FF) and Larkin-Ovchinnikov (LO) states denote modulations of a phase and amplitude of a condensation, respectively. The LO states were shown to be ground states in the presence of a magnetic field inducing the spin imbalance for a Cooper pair of a superconductor [4]. In the last couple of years, there a e-mail: nitta@phys-h.keio.ac.jp b e-mail: shin-s@kitasato-u.ac.jp c e-mail: ryokokur@keio.jp have been several claims of its observation (see Ref. [5] for a review). Recently, ultracold atomic Fermi gases have renewed interest in FFLO states (see Refs. [6,7] for a review). The spin polarized superfluid state was observed in Ref. [8] and it was claimed that the FFLO state has been achieved in this experiment. FFLO states in a ring were also proposed in cold Fermi gases [9] and in superconductors [10].
FFLO states, called twisted kink crystals, were also studied in the chiral Gross-Neveu model in 1+1 dimensions [11][12][13] (see [14] for application to a superconductor). Spatially modulated chiral condensations, such as FF states (called dual chiral density wave or chiral spiral) [15,16] and LO states (called real kink crystal) [17,18], have been proposed to appear in a certain region of the phase diagram of QCD in 3+1 dimensions (see Ref. [19] as a review). Although the Cooper pair is usually refers to the particle-particle condensates, the chiral condensation is related to the particleantiparticle (or hole) pairing. They were also proposed in diquark condensations exhibiting color superconductivity in high density QCD (see [20,21] as a review) and were also discussed in the context of the AdS/CFT correspondence [22][23][24][25].
These spatial modulations were originally proposed in condensations of fermions forming Cooper pairs. In terms of the Ginzburg-Landau effective theory, which is a scalar field theory, these states are realized as ground states of the theory due to the presence of a wrong sign of a gradient term and positive higher derivative terms. In general, these kinds of inhomogeneous states spontaneously break translational as well as rotational symmetries. Nambu-Goldstone (NG) modes associated with these broken symmetries in such backgrounds were studied in Refs. [26,27]. After all, inhomogeneous states in condensed matter, nuclear matter and quark matter studied so far are all realized in theories where the Lorentz invariance is explicitly broken due to the finite density/temperature effects and so on.
In this paper, we study spatially modulated vacua at zero temperature and zero density (but not ground states in finite density and/or temperature) in manifestly Lorentz invariant field theories, with a particular attention to spontaneous symmetry breaking and NG bosons. From a viewpoint of lowenergy effective theories, field theories generically receive higher derivative corrections.
We assume that there is no ghost in the theory implying the absence of more than one derivative on one field which can not be eliminated by partial integration. For example, the term ∂ 2 ϕ = ∂ m ∂ m ϕ with space-time index m, generally causes the so-called Ostrogradski instability [28]. This is a crucial difference with non-relativistic cases. Then, all higher derivative terms come in a way that only a single space-time derivative acting on one field, ∂ m ϕ. Thus, the effective theory is in general a function of ∂ m ϕ (complemented by a potential term). In this set-up we study an adaptation of the NG theorem to higher derivative theories , stating that when a global symmetry is spontaneously broken due to vacuum expectation values of space-time derivatives of fields, an NG boson appears without canonical kinetic (quadratic derivative) terms with a quartic derivative term in the modulated direction, while a Higgs boson appears with a non-zero canonical kinetic term.
After giving general discussion of the stability of general higher derivative models, we give a simple model illustrating this. Our model admits (meta)stable modulated vacuum of a phase modulation (Fulde-Ferrell state), where an NG mode associated with spontaneously broken translational and U (1) symmetries appears.

Adaptation of the Nambu-Goldstone theorem to higher derivative theories
In this section, we apply the NG theorem to the case that global symmetries of a Lagrangian are spontaneously broken due to vacuum expectation values (VEVs) of space-time derivatives of fields. Here, we consider the case that there is no Ostrogradski instability [28] assuming that there is no more than one space-time derivative on a field. We show that an analogue of the NG boson appears without canonical kinetic term with a quartic derivative term. In addition, we will show that a Higgs boson, which is defined by a mode that is orthogonal to the above mentioned NG mode, appears with a non-zero canonical kinetic term in the vacuum. In the following, we consider d-dimensional relativistic field theories where the Lorentz invariant Lagrangian L is given by a functional of ∂ m ϕ a . Here m = 0, 1, . . . , d − 1 is the space-time index and ϕ a (a = 1, . . . , N ) are complex scalar fields. The energy functional E of the theories depends only on the first time derivative of fields which we denote I = ∂ m ϕ a , † I = ∂ mφa . The index I = 1, . . . , d N labels the fields and directions of the space-time derivatives. Vacua |0 in the theories are defined such a way that they provide extrema of the energy E with respect to I , † I : In these extrema, we assume that the fields I , † I develop VEVs: Here some of these VEVs are not zero and they need not to be constants in general. Indeed, as we will see later, they are spatially varying functions for modulated vacua. Since I are in fact given by the space-time derivative of the fields (therefore they are Lorentz vectors), the non-zero VEVs (2) generically break the translational and rotational symmetries.
Hereafter, we assume that VEVs are spacelike vectors. 1 Now we introduce the dynamical fields˜ I = ∂ mφa ,˜ † I = ∂ mφ † a as fluctuations around a vacuum determined by the condition (1). We shift the fields around the VEVs, I → v I +˜ I , and the energy is expanded as where we have defined the following Hermitian matrix: Here the symbol * | v stands for the values evaluated in the vacuum. We note that the matrix M determined by the second derivatives of E is just the curvature of the energy density and it is in general a function of x i (i = 1, 2, 3). In order that the extrema defined by (1) become local minima of the energy, M should be a positive semi-definite matrix for all the regions in x. These conditions of vacua do not guarantee that they are global minima but meta-stable local minima are allowed in general. From the expression (3), one observes that the eigenvalues of M correspond to coefficients of the quadratic kinetic terms of the dynamical fieldsφ a ,φ † a . Since M is a positive semi-definite matrix, there are no fluctuation modes whose kinetic terms have the wrong sign (negative sign in the energy functional). However we stress that there are possible zero eigenvalues for a general M. When M has zero eigenvalues, the quadratic terms of the corresponding modes vanish.
In order to see the meaning of the zero eigenvalues of M, we elucidate the relation between the matrix M and the spontaneous symmetry breaking. The fields I , † I transform according to symmetries of theories: Here Q A are generators of the symmetry groups and the Hermitian matrices T A are an irreducible representation of Q A . In a vacuum |0 , some of the fields I develop non-zero VEVs and we have We now define the following vector, For generators that satisfy TÂ v = 0, the corresponding symmetry is preserved in the vacuum while for T A v = 0 the symmetry is spontaneously broken. The energy functional E is invariant under the following transformation, where ε A are infinitesimal real parameters. Then we have By differentiating the above relation with respect to I , † I and evaluate the result in a vacuum, we find Therefore T A v are eigenvectors associated with the zero eigenvalues of M. The relation (10) indicates that when some of the fields I , † I develop VEVs that spontaneously break symmetries, then the canonical quadratic kinetic terms for the modes that correspond to the zero eigenvalues of M vanish. We call these Nambu-Goldstone (NG) modes. On the other hand, the modes that are orthogonal to the NG modes appear with quadratic kinetic terms in the energy functional. We call these Higgs modes. We note that since the vector T A v generically depends on x i , there is a possibility that T A v(x) vanishes at some specific points x i = x i 0 in a general setup. At these points, the broken symmetries are recovered locally and one expects that a non-zero quadratic term associated with the NG mode recovers. We do not exclude this possibility but it is nevertheless not always the case.
Indeed, as we will show in an explicit example of the spatially dependent VEV in the next section, the vector T A v never vanishes at special points and the theorem discussed in this section completely works well in all regions in spacetime.

A model for spatially modulated vacua
In order to understand the discussion in the previous section concretely, we introduce a Lorentz invariant scalar field model where, in addition to the canonical quadratic kinetic terms, higher derivative corrections are involved. We begin with the observation that the global stability of modulation is guaranteed when the highest power of the derivative terms |∂ m ϕ| 2 is odd and an appropriate sign of the terms are chosen. We propose a scalar field model of the simplest example where a spatially modulated state is allowed as a (meta-)stable vacuum. We then apply the Nambu-Goldstone theorem discussed in the previous section to the model and show that there are modes where the quadratic kinetic terms vanish (NG modes). We demonstrate that there are always associated modes with the non-zero quadratic kinetic terms (Higgs modes).

Global stability of modulation
Let us consider a complex scalar field ϕ. The general Lorentz invariant Lagrangian containing finite powers of |∂ϕ| 2 = ∂ m ϕ∂ mφ is where n ∈ Z is the highest power of the derivative terms and · · · implies lower orders. The space-time index m is contracted by η mn = diag.(−1, 1, 1, 1). The dot inφ stands for the derivative of the field with respect to x 0 and ∇ is spatial derivatives. The canonical conjugate momentum is Then, the Hamiltonian associated with the Lagrangian (11) reads Let us discuss the stability of a vacuum in the model. First, by looking at the second term in (13), we see that the energy is bounded from below only when one chooses the upper sign in (11). Otherwise, the spatial gradient of the field causes an instability as |∇ϕ| 2 → ∞. Second, the first term in (13) implies that the energy is bounded from below only when the highest order n is odd. For even n, an instability in the temporal direction grows |φ| 2 → ∞. Therefore, the simplest Lagrangian is of the third order in |∂ϕ| 2 containing six derivatives. In the next subsection, we consider an example of the third order Lagrangian allowing a modulated vacuum.

A model and vacua
We propose a four-dimensional Lorentz invariant complex scalar field model whose Lagrangian is given by Here k > 0, λ > 0, α > 0 are real constants. 2 The Lagrangian (14) contains the ordinary kinetic term of the complex scalar field ϕ and the higher derivative corrections. The Lagrangian (14) is invariant under a global U (1) transformation ϕ → e iθ ϕ with a constant θ , in addition to the Poincaré symmetry including the SO(3, 1) Lorentz and translational symmetries. The Lagrangian can contain a potential term of ϕ too. In this paper, for simplicity we do not consider a potential term. In this case, the Lagrangian possesses a shift symmetry where c is a constant. The energy functional associated with the Lagrangian (14) is calculated to be where i, j = 1, 2, 3, π ϕ , πφ are the canonical momenta for ϕ,φ andφ = ∂ϕ ∂ x 0 ,φ = ∂φ ∂ x 0 . The spatially modulated vacua (ground states) are expected to appear at extrema of the energy functional (16).
We now employ an ansatz ϕ = ϕ(x 1 ) for one-dimensional spatial modulation along the x 1 -direction: 0|∂ 1 ϕ|0 = 0. We also assume static configurations. Then the energy functional (16) becomes The function E(x) is interpreted as a potential for x = |∂ 1 ϕ| 2 ≥ 0. It is obvious that E(x)(x ≥ 0) has a local minimum at x = 0 in which the vacuum energy is E(0) = 0 and the scalar field ϕ has a constant VEV. Whether E(x) has another minimum or not crucially depends on the parame- has another vacuum. In this case, the function E(x) has extrema at Since k > 0, x = x − corresponds to a local maximum while x = x + = 0 is a minimum which is a candidate of a modulated vacuum. Note that λ should be positive in order that x + > 0. The condition α > 0 is necessary in order that the potential is bounded from below. We find that the vacuum energy is classified according to the discriminant condition of the function αx 2 − λx + k. We have three distinct types of vacua. When the parameters k, λ, α satisfy the condition λ 2 − 4αk < 0, then the function E(x) becomes positive definite. The local vacuum energy at x = x + is positive E(x + ) > 0. This is a meta-stable vacuum which decays to the global minimum (true vacuum) x = 0 within a finite time. See Fig. 1a for the potential profile. When the parameters satisfy the condition λ 2 − 4αk = 0, the potential function E(x) looks like Fig. 1b. In addition to the global vacuum x = 0, we have a local vacuum x = x + in which E(x + ) = 0. They are actually global vacua and are degenerated. Finally, when the condition λ 2 − 4αk > 0 is satisfied, the function E(x) looks like Fig. 1c. Now the vacuum x = 0 becomes meta-stable and the vacuum at x = x + is energetically favoured. Then the true vacuum is located In each vacuum we have |∂ 1 ϕ| 2 = x + . The general solution that satisfies this relation is where c is a constant and F(s) is a real function. We are interested in a spatially modulated vacuum state. The most conservative choice is a linear function F(s) = ps where p is a constant. As we will see below, this vacuum preserves the highest symmetry in the theory. The vacuum is then given by where the constants p, ϕ 0 satisfy p 2 |ϕ 0 | 2 = x + . This is the ground state where the spatial modulation along the x 1direction occurs. The period of the modulation is given by 2π/ p. It is easy to confirm that the modulated vacuum (20) satisfies the equation of motion as follows. The equation of motion forφ in the model (14) is given by For a configuration ϕ = ϕ(x 1 ), we have In the modulated vacuum (20), we have |∂ 1 ϕ| 2 = x + = const. Therefore the Eq. (22) becomes Since 3αx 2 + − 2λx + + k = 0, we have shown that the modulated vacuum (20) satisfies the equation of motion.
We have a spatially modulated vacuum (20) along the x 1 -direction in which 0|∂ 1 ϕ|0 = v = i pϕ 0 e i px 1 = 0. It is obvious that the translational symmetry along the x 1direction and the rotational symmetries in the (x 1 , x 2 ) and (x 1 , x 3 ) planes are spontaneously broken in the modulated vacuum. Then the four-dimensional Poincaré symmetry is broken down to that in three dimensions.
As a matter of fact, due to the U (1) symmetry ϕ → e iθ ϕ, the simultaneous operation of the translation x 1 → x 1 + a and the U (1) transformation ϕ → e −i pa ϕ is preserved in the modulated vacuum (20). Here a is a constant. Meanwhile, the combination of the translation and the inverse rotation x 1 → x 1 + a, ϕ → e +i pa ϕ is broken. Therefore the global symmetry including the translational operation along the x 1 -direction is broken in such a way that Here P 1 represents translation along the x 1 -direction and "sim" means the simultaneous operation. As we have remarked above, this symmetry breaking pattern is a consequence of the simplest choice of the modulated vacuum (20). No other choices of F(s) results in this [P 1 ×U (1)] sim symmetry. Note that the translation P 1 and the rotations in the (x 1 , x 2 ) and (x 1 , x 3 ) planes are not independent each other [32]. Therefore we expect that there is one NG mode associated with the spontaneous symmetry breaking P 1 × U (1) → [P 1 × U (1)] sim . 3 We will clarify this issue in the followings subsections.

Linear analysis: Nambu-Goldstone and Higgs bosons
In this subsection, we show the stability analysis in our model, but the analysis employed here is completely general at the linear level for any model admitting local vacuum exhibiting a spatial modulation. We now shift the field from the modulated vacuum (20) and introduce the fluctuationφ as a dynamical field: where ϕ = ϕ 0 e i px 1 is the modulating VEV. In the following, we will show that there are no fluctuation modes that cause instabilities of the vacuum (20). The quadratic terms of the dynamical scalar fieldφ are extracted from the energy: Here the indexm = 0, 2, 3 is contracted by δmn. The vector ϕ and the matrix M are defined by Each block diagonal sector is found to be We have separated the quadratic terms to the SO(2, 1) Lorentz invariant sector (transverse direction) and the direction of the modulation. Since M is an Hermitian matrix, it is diagonalized by a unitary matrix: Here U 1 , U 2 are 2×2 unitary matrices. The eigenvalues s 1 , s 2 of M 1 are found to be It is easy to show that s 1 > 0 and s 2 = 0 by the definition of x + . The eigenvalues t 1 , t 2 of M 2 are calculated as Again we find t 1 > 0 and t 2 = 0. There are positive and zero eigenvalues as anticipated. This implies that our assumptioṅ ϕ = ∂ 2 ϕ = ∂ 3 ϕ = 0 guarantees the minimization condition of the energy. Therefore, the eigenvalues and eigenvectors of M are given by The unitary matrices are We are faced with the fact that there are modes whose quadratic kinetic terms disappear in the transverse (s 2 = 0) and the modulation (t 2 = 0) sectors. In order to understand the nature of the zero eigenvalues of M, we analyze the broken generator of the symmetry in the vacuum. The vacuum vector v is non-zero only in the modulation sector. Namely we have ∂mϕ = 0 and The action of the translation P 1 and the U (1) transformation generators T P 1 , T U (1) on the VEVs are given by The generator associated with the unbroken symmetry is defined by T ub = T P 1 − T U (1) . Indeed, we find that the action of T ub on v gives the vanishing result T ub v = 0. On the other hand, the generator associated with the broken symmetry is given by T b = T P 1 + T U (1) . Then we find This is exactly the eigenvector for the zero eigenvalue t 2 = 0 in the modulation sector. We note that the other zero eigenvalue s 2 = 0 corresponds to the flat direction in the SO(2, 1) invariant sector which does not accompany with the spontaneous symmetry breaking. By using the unitary matrix U 2 in (32), the matrix M 2 is diagonalized and we derive the Higgs and the NG modes associated with t 1 > 0, t 2 = 0 as These are the linear combinations of the fields that contain the derivative along the x 1 -direction. It is natural to define the modes (36) as the derivative of the Higgs and NG fields ∂ 1 H (x) and ∂ 1 N (x). As we have clarified, the quadratic term of the NG field N (x) vanishes in the energy functional (Lagrangian) while the Higgs field H (x) appears without a mass gap. Although, the NG and the Higgs modes emerge as an derivative modes in the modulation direction, we are interested in the modes that propagate in the transverse directions to the modulation. In order to examine this, we perform the linear transformation of the upper half part of ϕ and diagonalize the left upper half of M. Using the unitary matrix U 1 in (32), the matrix M 1 is diagonalized and we find that the Lagrangian in the quadratic order of fields is rewritten as Here and hereafter, the indexm is contracted by ηmn = diag. (−1, 1, 1), and we have defined Hereφm ,c is the mode associated with s 1 > 0 and therefore it has the quadratic canonical kinetic term. On the other hand, ϕm ,0 is the one for s 2 = 0 whose quadratic kinetic term vanishes. Note that they are distinguished from the Higgs and the NG modes in the modulation direction. Again we define the fields ∂m A, ∂m B that correspond to the modes (38). Then, the linear transformations (38) are interpreted as the following field re-definition: 4 Then, the NG and the Higgs modes in the modulation directions are represented in terms of A, B: We obtain the Lagrangian for the dynamical fields A, B in the modulated vacuum at the quadratic order as One observes that the field A does not propagate in the modulation direction while B does not propagate in the transverse direction. Only the gradient of B in the modulation direction contributes to the energy. This is a reflection of the fact that the term ∂ 1 A is included in the NG modeφ NG and it never appears in the Lagrangian at the quadratic order. A similar analysis has been done in [35] for the dispersion relation of NG and Higgs modes in a plane-wave type ground state in a Lorentz non-invariant theory. We did not consider a potential term for ϕ, and consequently the system has the shift symmetry in Eq. (15). What we have identified as a "Higgs boson" here is actually an NG boson associated with the spontaneous breaking of the shift symmetry. If we add a potential term in the original Lagrangian, the Higgs boson obtains a mass. Therefore, the gapless property is originated from the shift symmetry, but what we have found here is the existence of the quadratic kinetic term of the Higgs boson, in contrast to its absence for the NG boson.

Higher order terms
Here, we study higher order expansion, and show that the cubic order of the expansion of the Lagrangian contains no term consisting of only the NG bosonφ 3 NG , while in the quartic order there exists a term consisting of only the NG bosoñ ϕ 4 NG . In general, we cannot exclude a possibility ofφ 3 NG a priori, since the cubic derivative terms (∂ 1φ ) 3 exist after translational symmetry along x 1 -direction and the rotational symmetries in the (x 1 , x 2 ) and (x 1 , x 3 ) are broken. To see this, we calculate the cubic derivative terms of fluctuation ∂φ in the Lagrangian (14) with the introduction of the fluctuation ϕ → ϕ +φ in (24). The explicit calculation leads to the cubic derivative terms of the Lagrangian L cub. as By using Eq. (36), we find where the ellipsis · · · means that the terms with ∂mφ. Thus, we find that the pure cubic NG termφ 3 NG vanishes. Now we see that the NG mode appears with a quartic derivative term. The quartic derivative terms in the Lagrangian (14) L quart. is similarly obtained as The quartic derivative term containing purely NG modeφ NG is found by using (36) as where the ellipsis · · · expresses the terms other thanφ 4 NG . Therefore, we conclude that a term consisting of only the NG mode appears with the quartic derivative term.

Conclusion and discussions
In this paper, we have studied the spatially modulated vacua in a Lorentz invariant field theory where no finite density/temperature effects are included. The NG theorem, for a global symmetry spontaneously broken due to vacuum expectation values of space-time derivatives of fields, states that there appears an NG boson without a canonical quadratic kinetic term but with a quadratic derivative term in the modulated direction and a Higgs boson. We demonstrated this in the simple model whose energy functional can be written by the derivative terms of the scalar fields. The potential for the derivative terms allows a local vacuum as the modulated vacuum, where translational symmetry along one direction which we choose x 1 and the rotational symmetries in the (x 1 , x 2 ), (x 1 , x 3 ) planes are spontaneously broken together with the U (1) symmetry. A simultaneous transformation of P 1 and U (1) is preserved in the modulated vacuum. This modulated vacuum is completely consistent with the equation of motion. The vacuum energy depends on the parameters. In this paper, we have focused on the vacuum where the modulation takes place only in the x 1 -direction. Since the transverse directions xm, (m = 0, 2, 3) preserve the SO(2, 1) Lorentz symmetry, there are no mixing terms between ∂ 1 ϕ and ∂mϕ in the Lagrangian in the quadratic order of fields. This results in the complete block diagonal form of the matrix M in (26). We have explicitly shown that the "mass eigenstates" in the modulation and the transverse directions are different. Therefore we are not able to perform the diagonalization in these directions simultaneously. We have employed the linear combination of the fields (39) and represented the NG and Higgs modes in terms of A and B. The A-mode propagates in the transverse directions while the B-mode only oscillates in the modulation direction. Finally, we have demonstrated that a term containing only NG modes appears in the quartic derivative order. The Higgs mode, which is defined as the orthogonal mode to the NG mode in our discussion, is indeed an NG mode associated with the spontaneously broken shift symmetry. We note that this Higgs mode has its non-zero quadratic kinetic term. This is a consequence of the application of the NG theorem to higher derivative field theories.
Although we have illustrated the stability of the modulated vacuum in our simplest model, we would like to emphasize that the stability analysis employed in this paper is general at the linear level for any model admitting local vacuum exhibiting a spatial modulation.
Among general solutions in Eq. (19) which are energetically degenerated, we have focused on the FF state, which has the highest unbroken symmetry. Which vacuum is chosen among energetically degenerated vacua with different unbroken symmetry is known as a vacuum alignment problem first discussed in the context of technicolor models [36,37]. In such the cases, quantum corrections pick up the vacuum with the highest unbroken symmetry, and therefore we expect the same happens in our case too. We note that the structure of the vacuum modulation crucially depends on models we consider. For example, inhomogeneous chiral condensates in dense QCD appears to be a FFLO instead of a FF one.
We have studied the modulated vacuum in the Ginzburg-Landau type effective theory in the Lorentz invariant framework, without assuming any underlying microscopic theory. However, there is an argument about a no-go theorem for modulated vacua by fermion condensates in relativistic QCD-like theories [38]. It is an interesting open question whether our model can be obtained as the low-energy theory of a fermion condensation in relativistic theories or not. Possible future directions are in order. Beyond the semiclassical level in this paper, we need a more rigorous proof for the generalized NG theorem in full quantum level. The Higgs mechanism in a U (1)-gauged model, and spatial modulations along two or more directions [39] as well as a temporal modulation [40] are interesting directions. Applying our discussion to more general higher derivative theories such as higher-order Skyrme model [41] is also one of future directions. We also would like to embed our model to supersymmetric theories based on the formalism in Refs. [42][43][44], and supersymmetry breaking in modulated vacua will be reported elsewhere [45].