Dancer–Yan Spikes in Critical Elliptic Problems

• Dancer–Yan spikes are spike solutions that emerge from the bifurcation of Aubin–Talenti bubbles when a compact direction alters the Euclidean concentration mechanism.
• They arise in energy-critical Lane–Emden equations on waveguide manifolds (ℝ^d × T) and are analyzed via variational methods with semivirial constraints and scaling arguments.
• In two-dimensional plasma free-boundary problems, these spikes exhibit a two-scale structure with refined inner–outer scaling and a distinctive Type I blow-up profile.

Dancer–Yan spikes are spike solutions associated with elliptic problems in which a compact or bounded direction alters the concentration mechanism available in the fully Euclidean setting. In the energy-critical Lane–Emden equation on the waveguide manifold $\mathbb R_x^d\times \mathbb T_y$, they are positive-frequency solutions that decay in the noncompact $x$-directions and are periodic in $y$, arising from the bifurcation of Aubin–Talenti bubbles after one direction is compactified (Luo, 1 Jun 2026). In a distinct two-dimensional singularly perturbed free-boundary problem from plasma physics, the same terminology designates a precise Type I blow-up profile characterized by a refined inner–outer scaling and an Emden–log limiting profile (Bartolucci et al., 28 Jul 2025).

1. Waveguide Lane–Emden setting

Let $D=d+1\ge 3$ and

$2^*=\frac{2D}{D-2}=2+\frac{4}{d-1}.$

On the waveguide

$$\mathcal X=\mathbb R_x^d\times \mathbb T_y,\qquad \mathbb T=\mathbb R/(2\pi\mathbb Z),$$

the relevant positive-frequency standing waves have the form $e^{i\sqrt\beta\,t}u(x,y)$ and satisfy

$$-\Delta_{x,y}u+\beta u=u^{2^*-1},\qquad u>0,\ \beta>0.$$

The variational formulation fixes the mass

$M(u)=\int_{\mathcal X}|u|^2,\qquad M(u)=c>0,$

and uses the energy

$$E(u)=\frac12\int_{\mathcal X}|\nabla_{x,y}u|^2-\frac{D-2}{2D}\int_{\mathcal X}|u|^{2^*}.$$

To isolate functions that “just fail to scatter” in the Euclidean directions, Luo introduces the semivirial functional

$x$0

and studies the constrained minimization problem

$x$1

A Lagrange-multiplier argument then shows that any optimizer solves the elliptic equation for some $x$2 (Luo, 1 Jun 2026).

This formulation is specific to the partially periodic geometry. The constraint $x$3 is not a generic Pohozaev identity imposed after the fact; it is built into the minimization problem to encode the Euclidean-direction criticality while keeping the periodic direction explicit.

2. Aubin–Talenti bubbles and the Dancer bifurcation picture

In the fully Euclidean space $x$4, the zero-frequency equation

$x$5

has the Aubin–Talenti family of optimizers. These satisfy

$x$6

where

$x$7

is the sharp Sobolev constant. They are the unique positive finite-energy solutions of the zero-frequency problem on $x$8, and on $x$9 they extend trivially by being constant in $y$0 (Luo, 1 Jun 2026).

Dancer’s bifurcation perspective begins precisely at this rigid Euclidean family. Once one breaks full Euclidean invariance by compactifying one direction to a torus, the bubble manifold bifurcates at small positive frequency $y$1 into genuinely $y$2-dependent solutions. The schematic Lyapunov–Schmidt expansion is

$y$3

with decay in $y$4 and $y$5-periodicity in $y$6. In Luo’s terminology, these are the Dancer–Yan spikes (Luo, 1 Jun 2026).

The bifurcation viewpoint is significant because it identifies the objects as perturbations of the sharp Sobolev optimizers. A plausible implication is that the waveguide problem should be analyzed relative to the Euclidean bubble threshold rather than by perturbative compactness alone.

3. Semivirial-vanishing geometry

Luo’s approach replaces the bifurcation construction with an energy-based variational method built around semivirial-vanishing geometry. The auxiliary functional is

$y$7

On the constraint set $y$8, one has $y$9. The geometry is organized by the $D=d+1\ge 3$0-scaling

$D=d+1\ge 3$1

for which

$D=d+1\ge 3$2

Two structural facts are central. First, if $D=d+1\ge 3$3, there is a unique $D=d+1\ge 3$4 such that $D=d+1\ge 3$5, and moreover $D=d+1\ge 3$6. Second,

$D=d+1\ge 3$7

These statements convert the problem into one where minimizing sequences can be projected onto the semivirial manifold without losing control of the energy (Luo, 1 Jun 2026).

The paper explicitly contrasts this mechanism with the usual Brezis–Nirenberg scenario. No lower-order focusing perturbation is available in the model. Instead, the energy drop is produced by the bounded periodic direction, while the focusing takes place through scaling in the Euclidean directions alone. This suggests that the compact geometry acts as an effective variational defect relative to the Euclidean critical problem.

4. Strict sub-bubbling and normalized ground states

The Euclidean bubble threshold is

$D=d+1\ge 3$8

The main variational input is the strict sub-bubbling estimate

$D=d+1\ge 3$9

To prove it, Luo truncates a rescaled Aubin–Talenti bubble to a single period slab in the periodic direction. The key asymptotics are that the lost gradient tail satisfies

$2^*=\frac{2D}{D-2}=2+\frac{4}{d-1}.$0

whereas the lost nonlinear tail is only $2^*=\frac{2D}{D-2}=2+\frac{4}{d-1}.$1. Correspondingly,

$2^*=\frac{2D}{D-2}=2+\frac{4}{d-1}.$2

for $2^*=\frac{2D}{D-2}=2+\frac{4}{d-1}.$3, and

$2^*=\frac{2D}{D-2}=2+\frac{4}{d-1}.$4

After projection back to $2^*=\frac{2D}{D-2}=2+\frac{4}{d-1}.$5 by the unique scaling $2^*=\frac{2D}{D-2}=2+\frac{4}{d-1}.$6, the energy changes only by $2^*=\frac{2D}{D-2}=2+\frac{4}{d-1}.$7, which produces a test function strictly below the bubble level (Luo, 1 Jun 2026).

The leading-order decomposition makes the mechanism more explicit. Writing

$2^*=\frac{2D}{D-2}=2+\frac{4}{d-1}.$8

one has

$2^*=\frac{2D}{D-2}=2+\frac{4}{d-1}.$9

with

$$\mathcal X=\mathbb R_x^d\times \mathbb T_y,\qquad \mathbb T=\mathbb R/(2\pi\mathbb Z),$$0

Thus the leading-order tail removal creates a negative contribution of order $$\mathcal X=\mathbb R_x^d\times \mathbb T_y,\qquad \mathbb T=\mathbb R/(2\pi\mathbb Z),$$1 in the gradient part, while the nonlinear loss is lower order. This geometric energy drop is what restores compactness below the sharp Sobolev threshold.

Once $$\mathcal X=\mathbb R_x^d\times \mathbb T_y,\qquad \mathbb T=\mathbb R/(2\pi\mathbb Z),$$2 is known, concentration of an Euclidean bubble can be excluded by the standard profile decomposition on $$\mathcal X=\mathbb R_x^d\times \mathbb T_y,\qquad \mathbb T=\mathbb R/(2\pi\mathbb Z),$$3. Weak vanishing is also ruled out, since $$\mathcal X=\mathbb R_x^d\times \mathbb T_y,\qquad \mathbb T=\mathbb R/(2\pi\mathbb Z),$$4 forces a uniform lower bound on the critical $$\mathcal X=\mathbb R_x^d\times \mathbb T_y,\qquad \mathbb T=\mathbb R/(2\pi\mathbb Z),$$5-norm. After compensating by at most an $$\mathcal X=\mathbb R_x^d\times \mathbb T_y,\qquad \mathbb T=\mathbb R/(2\pi\mathbb Z),$$6-translation, a minimizing sequence has a nonzero weak limit, and Brezis–Lieb splitting of mass, energy, and semivirial implies that the limit itself satisfies

$$\mathcal X=\mathbb R_x^d\times \mathbb T_y,\qquad \mathbb T=\mathbb R/(2\pi\mathbb Z),$$7

A secondary scaling argument yields $$\mathcal X=\mathbb R_x^d\times \mathbb T_y,\qquad \mathbb T=\mathbb R/(2\pi\mathbb Z),$$8, and positivity follows by the maximum principle. Therefore, for every $$\mathcal X=\mathbb R_x^d\times \mathbb T_y,\qquad \mathbb T=\mathbb R/(2\pi\mathbb Z),$$9, the infimum $e^{i\sqrt\beta\,t}u(x,y)$0 is attained by a positive solution

$e^{i\sqrt\beta\,t}u(x,y)$1

The same framework also records that there exists $e^{i\sqrt\beta\,t}u(x,y)$2 such that if $e^{i\sqrt\beta\,t}u(x,y)$3, any optimizer $e^{i\sqrt\beta\,t}u(x,y)$4 has $e^{i\sqrt\beta\,t}u(x,y)$5 (Luo, 1 Jun 2026).

5. Two-dimensional free-boundary formulation and the classical spike ansatz

A second appearance of Dancer–Yan spikes occurs in the two-dimensional free-boundary problem studied by Bartolucci, Jevnikar, Wei, and Wu. Let $e^{i\sqrt\beta\,t}u(x,y)$6 be a bounded $e^{i\sqrt\beta\,t}u(x,y)$7-domain, with $e^{i\sqrt\beta\,t}u(x,y)$8 and $e^{i\sqrt\beta\,t}u(x,y)$9. The plasma problem is

$$-\Delta_{x,y}u+\beta u=u^{2^*-1},\qquad u>0,\ \beta>0.$$0

The plasma region is $$-\Delta_{x,y}u+\beta u=u^{2^*-1},\qquad u>0,\ \beta>0.$$1, and $$-\Delta_{x,y}u+\beta u=u^{2^*-1},\qquad u>0,\ \beta>0.$$2 is the free boundary. After introducing $$-\Delta_{x,y}u+\beta u=u^{2^*-1},\qquad u>0,\ \beta>0.$$3, with $$-\Delta_{x,y}u+\beta u=u^{2^*-1},\qquad u>0,\ \beta>0.$$4, the paper rewrites the problem in the equivalent form

$$-\Delta_{x,y}u+\beta u=u^{2^*-1},\qquad u>0,\ \beta>0.$$5

In the singular limit $$-\Delta_{x,y}u+\beta u=u^{2^*-1},\qquad u>0,\ \beta>0.$$6, equivalently $$-\Delta_{x,y}u+\beta u=u^{2^*-1},\qquad u>0,\ \beta>0.$$7, one recovers

$$-\Delta_{x,y}u+\beta u=u^{2^*-1},\qquad u>0,\ \beta>0.$$8

with

$$-\Delta_{x,y}u+\beta u=u^{2^*-1},\qquad u>0,\ \beta>0.$$9

Unlike the higher-dimensional case, the usual rescalings $M(u)=\int_{\mathcal X}|u|^2,\qquad M(u)=c>0,$0 fail in dimension two because the planar limit

$M(u)=\int_{\mathcal X}|u|^2,\qquad M(u)=c>0,$1

is impossible (Bartolucci et al., 28 Jul 2025).

Dancer and Yan’s remedy is a two-scale construction. Let $M(u)=\int_{\mathcal X}|u|^2,\qquad M(u)=c>0,$2 solve

$M(u)=\int_{\mathcal X}|u|^2,\qquad M(u)=c>0,$3

radial and decreasing, and define

$M(u)=\int_{\mathcal X}|u|^2,\qquad M(u)=c>0,$4

With $M(u)=\int_{\mathcal X}|u|^2,\qquad M(u)=c>0,$5 so that $M(u)=\int_{\mathcal X}|u|^2,\qquad M(u)=c>0,$6, and with $M(u)=\int_{\mathcal X}|u|^2,\qquad M(u)=c>0,$7, $M(u)=\int_{\mathcal X}|u|^2,\qquad M(u)=c>0,$8 is chosen by

$M(u)=\int_{\mathcal X}|u|^2,\qquad M(u)=c>0,$9

The Dancer–Yan trial function in $$E(u)=\frac12\int_{\mathcal X}|\nabla_{x,y}u|^2-\frac{D-2}{2D}\int_{\mathcal X}|u|^{2^*}.$$0 is

$$E(u)=\frac12\int_{\mathcal X}|\nabla_{x,y}u|^2-\frac{D-2}{2D}\int_{\mathcal X}|u|^{2^*}.$$1

It satisfies

$$E(u)=\frac12\int_{\mathcal X}|\nabla_{x,y}u|^2-\frac{D-2}{2D}\int_{\mathcal X}|u|^{2^*}.$$2

In particular, when $$E(u)=\frac12\int_{\mathcal X}|\nabla_{x,y}u|^2-\frac{D-2}{2D}\int_{\mathcal X}|u|^{2^*}.$$3 and $$E(u)=\frac12\int_{\mathcal X}|\nabla_{x,y}u|^2-\frac{D-2}{2D}\int_{\mathcal X}|u|^{2^*}.$$4, it is an exact solution of

$$E(u)=\frac12\int_{\mathcal X}|\nabla_{x,y}u|^2-\frac{D-2}{2D}\int_{\mathcal X}|u|^{2^*}.$$5

The asymptotic features of this spike are also explicit. One has $$E(u)=\frac12\int_{\mathcal X}|\nabla_{x,y}u|^2-\frac{D-2}{2D}\int_{\mathcal X}|u|^{2^*}.$$6 and

$$E(u)=\frac12\int_{\mathcal X}|\nabla_{x,y}u|^2-\frac{D-2}{2D}\int_{\mathcal X}|u|^{2^*}.$$7

The inner core width is $$E(u)=\frac12\int_{\mathcal X}|\nabla_{x,y}u|^2-\frac{D-2}{2D}\int_{\mathcal X}|u|^{2^*}.$$8, the free-boundary layer is logarithmic in $$E(u)=\frac12\int_{\mathcal X}|\nabla_{x,y}u|^2-\frac{D-2}{2D}\int_{\mathcal X}|u|^{2^*}.$$9, and the peak height obeys

$x$00

so

$x$01

(Bartolucci et al., 28 Jul 2025).

6. Singular-limit classification, vanishing level, and multi-spike structure

The 2025 classification theorem concerns nonnegative solutions $x$02 of

$x$03

under the uniform bounds

$x$04

where $x$05 is defined by

$x$06

and $x$07. If $x$08 is a local maximizer with $x$09, then exactly one of four behaviors occurs (Bartolucci et al., 28 Jul 2025).

Case	Limiting description	Mass behavior
(A) Vanishing	$x$10 on compact subsets eventually	no local spike
(B–i) Type I	$x$11 on expanding balls	both $x$12- and $x$13-mass quantize
(B–ii) Type II	$x$14, $x$15	$x$16-mass quantizes, $x$17-mass vanishes in the spike
(B–iii) Fading spike	$x$18 bounded away from $x$19	both masses decay too fast for quantization

The refined scaling is

$x$20

followed by

$x$21

The classical Dancer–Yan spike is exactly case (B–i): $x$22, and on each ball $x$23 with $x$24,

$x$25

where $x$26 equals the Emden solution $x$27 for $x$28 and $x$29 outside. In this regime,

$x$30

and

$x$31

The “vanishing level” is simply the threshold $x$32, and the geometry near a Type I or II spike is encoded by

$x$33

Hence the free boundary lies at distance $x$34 from $x$35, while the amplitude above $x$36 is

$x$37

Equivalently,

$x$38

The global theory with Dirichlet boundary data $x$39 and non-vanishing total $x$40-mass shows that boundary regions contain no spikes, that only finitely many spike sequences occur, and that these are of Type I or II. The centers are separated by distances much larger than the local core scales, and away from the union of $x$41 one has $x$42 by Green’s representation. Matching inner expansions to outer Green expansions quantizes the $x$43-mass to $x$44, while a Pohozaev identity around each spike yields drift equations forcing the limiting centers to be critical points of the Kirchhoff–Routh Hamiltonian

$x$45

with

$x$46

A central clarification provided by the classification is that, in dimension $x$47, it is not true that any solution in the singular limit is a Dancer–Yan spike. The spiking structure is more rich. By contrast, for $x$48, simple $x$49-rescaling around a nondegenerate maximum leads to a unique radially decaying solution of

$x$50

and that classification forces a single spike type (Bartolucci et al., 28 Jul 2025).