Spatial Periodic Forcing Can Displace Patterns It Is Intended to Control

Spatial periodic forcing of pattern-forming systems is an important, but lightly studied, method of controlling patterns. It can be used to control the amplitude and wave number of one-dimensional periodic patterns, to stabilize unstable patterns, and to induce them below instability onset. We show that, although in one spatial dimension the forcing acts to reinforce the patterns, in two dimensions it acts to destabilize or displace them by inducing two-dimensional rectangular and oblique patterns. Pattern formation phenomena are found in a wide variety of physical, chemical, and biological contexts. Examples include embryonic pattern formation [1], cardiac arrhythmias [2], bacterial colonies [3], nanoparticle assemblies [4], two-phase mixtures [5], thermal convection [6], nonlinear optics [7], chemical [8] and electrochemical reactions [9], and environmental pattern formation [10]. In some contexts, pattern formation is essential for the functioning of the system. This is the case with embryonic pattern formation or with vegetation patterning—a mechanism by which vegetation copes with water stress. In other contexts, pattern formation is an undesired outcome. This is the case with spiral waves in the heart muscle [11], dewetting of liquid films [12], or spatial patterning in the transverse directions of a laser beam [13]. In order to eliminate, modify, or induce patterns, means of controlling and manipulating them are needed. These means may consist of basic parameter tuning or may involve external intervention such as feedback control [14] or periodic forcing in time [15] and space [16]. Control of periodic patterns by spatial periodic forcing is achieved by locking the pattern's wave number, k, to a rational fraction of the forcing wave number, k f. The locking typically occurs over a limited wave-number range that increases with the forcing amplitude—the ''resonance tongue.'' Within this range, the wave number of the locked or resonant pattern is controllable by tuning the forcing wave number. Recent studies have shown that parametric one-dimensional (1d) spatial forcing of systems supporting stationary stripe patterns can also induce resonant two-dimensional (2d) patterns of rectangular and oblique forms [16,17]. Such patterns exist over a wide range of forcing wave numbers that includes resonance tongues of stripe patterns, in particular, the basic 1:1 (k ¼ k f) resonance, which would generally be the first choice for control. In this Letter, we address the question of how the 2d rectangular and oblique patterns interfere with the control of stripe patterns. Focusing on universal aspects of spatially forced pattern-forming systems, …

Pattern formation phenomena are found in a wide variety of physical, chemical, and biological contexts.Examples include embryonic pattern formation [1], cardiac arrhythmias [2], bacterial colonies [3], nanoparticle assemblies [4], two-phase mixtures [5], thermal convection [6], nonlinear optics [7], chemical [8] and electrochemical reactions [9], and environmental pattern formation [10].In some contexts, pattern formation is essential for the functioning of the system.This is the case with embryonic pattern formation or with vegetation patterning-a mechanism by which vegetation copes with water stress.In other contexts, pattern formation is an undesired outcome.This is the case with spiral waves in the heart muscle [11], dewetting of liquid films [12], or spatial patterning in the transverse directions of a laser beam [13].In order to eliminate, modify, or induce patterns, means of controlling and manipulating them are needed.These means may consist of basic parameter tuning or may involve external intervention such as feedback control [14] or periodic forcing in time [15] and space [16].
Control of periodic patterns by spatial periodic forcing is achieved by locking the pattern's wave number, k, to a rational fraction of the forcing wave number, k f .The locking typically occurs over a limited wave-number range that increases with the forcing amplitude-the ''resonance tongue.''Within this range, the wave number of the locked or resonant pattern is controllable by tuning the forcing wave number.
Recent studies have shown that parametric onedimensional (1d) spatial forcing of systems supporting stationary stripe patterns can also induce resonant twodimensional (2d) patterns of rectangular and oblique forms [16,17].Such patterns exist over a wide range of forcing wave numbers that includes resonance tongues of stripe patterns, in particular, the basic 1:1 (k ¼ k f ) resonance, which would generally be the first choice for control.In this Letter, we address the question of how the 2d rectangular and oblique patterns interfere with the control of stripe patterns.
Focusing on universal aspects of spatially forced pattern-forming systems, we study the Swift-Hohenberg (SH) equation as a minimal model that captures the relevant mathematical construct of a stationary nonuniform instability of a stationary uniform state [18].Adding parametric forcing, the equation reads where " is the distance from the pattern-forming instability of the uniform stationary state, u ¼ 0, of the unforced system; k 0 $ Oð1Þ is the wave number of the mode that begins to grow at the instability point; k f is the forcing wave number; and > 0 is the forcing amplitude [19].
Resonant stripe patterns of Eq. ( 1) exist in tongue-shaped domains in the forcing parameter plane k f =k 0 À.We begin by identifying these domains for patterns near the instability point j"j ( 1.We consider Eq. ( 1) in 1d and approximate the resonant stripe solutions as where the amplitude A is small, jAj $ Oð ffiffiffiffiffiffi j"j p Þ, and slowly varying in space and time, jA x j $ Oðj"jÞ, jA t j $ Oðj"j 3=2 Þ, and c.c. stands for the complex conjugate.The parameter n ¼ 1; 2; . . . is an integer representing the type of resonance k f :k ¼ n:1.Within and in the vicinities of the resonant tongues, the detunings n ¼ k 0 À k f =n from the exact resonances n:1 are small.Assuming n $ Oð ffiffiffiffiffiffi j"j p Þ, $ Oð ffiffiffiffiffiffi j"j p Þ and using multiple-scale analysis to order j"j 3=2 , we find the amplitude equation where and n;1 is the Kronecker delta.Note that d AE diverges for k f ¼ Ç2k 0 , that is, for the exact 2:1 resonance.This resonance, which requires a different scaling of the forcing, $ Oðj"jÞ, was studied earlier [17].
Constant solutions of Eq. ( 3) represent n:1 wave-numberlocked, or resonant, stationary stripe patterns.For n Þ 2, they have the form The phase is constant and equal to zero for n ¼ 1 but undetermined for higher resonances for the order j"j 3=2 of our calculation.The resonant stripe solutions exist for > n , where Wave-number-locked patterns corresponding to the 2:1 resonance have the form [17] with ¼ 0; .These solutions exist for > 2 , where Figure 1 shows the tongue-shaped existence ranges of n:1 resonant stripe patterns with n ¼ 1; . . .; 4, for parameters above, " > 0, and below, " < 0, the pattern-forming instability.The solid lines in the figure are the results of the analysis from Eqs. ( 5) and (7), and the shaded regions are numerical results from solving for stationary solutions of the SH equation (1) using a continuation method [20].
Of all resonances shown in Fig. 1, the 2:1 resonance region stands out in its robustness.It is wider and, for " < 0, i.e., below the pattern-forming instability, it appears at lower forcing amplitude .This is because the forcing is parametric, involving the linear term in the SH equation.Parametric forcing of higher-order terms will single out higher resonances.
A stronger expression of the special role the 2:1 resonance plays appears in two-space dimensions.In that case, a purely 1d forcing, k f ¼ k f x, can induce stable 2d patterns [17]-oblique patterns for < " and rectangular patterns for > ".These are resonant patterns that respond to the spatial forcing by locking the wave vector components in the forcing direction in 2:1 resonance, k x ¼ k f =2, and creating a wave vector component in the orthogonal direction, k y , to compensate for the unfavorable forcing wave number, so that The range of existence of these new 2d patterns is very wide in the forcing wave number k f but bounded from above by the 2:1 resonance of stripe patterns, i.e., 0<k f <2k 0 , since at k f ¼ 2k 0 the component k x attains its maximal possible value, k 0 .Figure 2 shows the existence domains of resonant rectangular and oblique patterns superimposed onto the tongue diagram of resonant 1:1 and 2:1 stripe patterns.
We now address the overlap domains of rectangular and oblique patterns with the 1:1 resonance tongue of stripe patterns.In order to study the interaction of the patterns, we approximate solutions to Eq. (1) as a superposition of a stripe mode and two oblique modes u % Ae ik f x þ ae iðk x xþk y yÞ þ be iðk x xÀk y yÞ þ c:c:; (9) where . Using multiple-scale analysis, we find the amplitude equations [21] FIG. 1. Existence domains of resonant stripe solutions of Eq. (1) (a) above the pattern-forming instability (" > 0) and (b) below it (" < 0).The shaded regions indicate the range of resonant solutions computed from stationary solutions of Eq. ( 1), and the solid curves show the approximate boundary of resonant solutions based on the amplitude equation approach.The agreement for the lower resonances is very good for sufficiently small values, and for the higher resonances it remains surprisingly good, even for large values.Parameters: (a) " ¼ 0:001, (b) " ¼ À0:001.
where (10) of the form ðA; 0; 0Þ represent 1:1 stripe patterns, while solutions of the form ð0; a; bÞ represent 2d patterns.Stationary stripe solutions are given by A ¼ 1 , where 1 is given by Eq. ( 5) with n ¼ 1. Stationary rectangular solutions are given by a ¼ R e i , b ¼ R e Ài , where is an arbitrary phase and while stationary oblique solutions are given by a where To see how the 2d resonant patterns affect the stability of the 1:1 resonant stripes, we analyze the stability of the solution ðA; a; bÞ ¼ ð 1 ; 0; 0Þ.The eigenvalue analysis of the two oblique modes shows that the stability region of the stripe solution is smaller than the existence region and has two distinct shapes, depending on the value of ".For larger values, " > " c ¼ 9=143, there is a continuous range in which stripe solutions are stable, while, for " < " c , the stability range is split into two ranges.Figure 3 shows the shapes of the stable solution ranges for two values of ", one above and one below the critical value " c .The significance is that for " < " c there is an intermediate range of forcing amplitude where the forcing destabilizes the stripe patterns even at exact resonance k f ¼ k 0 [22].This is in contrast to the behavior of 1d systems for which the forcing always acts to stabilize the stripe patterns.
The stability ranges of both the 1:1 and 2:1 stripe patterns are actually bistability ranges of the stripes and 2d patterns [23].This raises the question of which pattern is dominant; that is, which pattern invades the other in these ranges.To study this question, we calculated the energy (Lyapunov) functional of Eq. ( 1), using the analytic forms of the approximate stripe, rectangular, and oblique solutions; the pattern that has lower energy is dominant.Figure 4 shows the energies of rectangular patterns, resonant 1:1 stripe patterns, and resonant 2:1 stripe patterns, in their existence range along the k f =k 0 axis.The energy of the 1:1 stripe pattern is higher than that of the rectangular pattern, implying that the latter is dominant.This is supported by numerical solutions of Eq. ( 1), according to which the rectangular patterns invade the 1:1 stripe, as the snapshots in Fig. 4(a) show.This result holds even at exact resonance (k f ¼ k 0 ).The situation is different within the resonance range of 2:1 stripes; the energies of the stripe and rectangular patterns cross one another and split the range into a low-k f part where the rectangular patterns are dominant and a high-k f part where the stripe pattern is Two intriguing results stand out in the analysis described above: (1) the most obvious control practice, i.e., 1:1 periodic forcing, can destabilize the stripe pattern it is intended to control and (2) even when the forcing leaves the stripes linearly stable, the stripes are displaced by the 2d patterns that the forcing induces.These are counterintuitive results because, naively, we would expect the forcing to reinforce the stripe patterns.The forcing indeed reinforces the stripe patterns by increasing their amplitudes [see Eq. ( 5)] and widening their wave-number range, but it also induces the growth of the oblique modes which either destabilize or displace the stripe patterns.
These outcomes may have important implications for practical applications of spatial forcing in various fields of science, including nonlinear optics [24,25] and restoration ecology.An interesting example of the latter field is rehabilitation of banded vegetation on hill slopes by water harvesting [26].Water-harvesting methods often involve parallel contour ditches that accumulate runoff and along which the vegetation is planted [27].Our results suggest that the system may not respond as expected in a 1:1 resonance but rather form 2d patterns that involve longlasting processes of mortality and regrowth.Because the analysis is based on universal amplitude equations near an instability point, we expect it to hold for a wide variety of spatially forced systems. The

FIG. 2 (
FIG.2(color online).The existence domains of resonant rectangular and oblique patterns are very wide and overlap with the resonance tongues of 1:1 and 2:1 stripe patterns.The domain size for all four patterns is 30 Â 30.Parameters: " ¼ 0:1.
support of the United States-Israel Binational Science Foundation (Grant No. 2008241) is gratefully acknowledged.Part of this work was funded by the Department of Energy at Los Alamos National Laboratory under Contract No. DE-AC52-06NA25396 and by the DOE Office of Science Advanced Computing Research (ASCR) program in Applied Mathematical Sciences.