Excited solutions in a Skyrme--Chern-Simons model in $2+1$ dimensions

This paper investigates excited solutions in a 2+1 dimensional Skyrme-Chern-Simons model, demonstrating that a Lagrange multiplier method is essential to resolve discontinuities in finding these states and revealing that while global charges exhibit non-standard parameter dependence, the fundamental solutions ($p=0$) consistently possess the minimal energy.

The Gist

Imagine the universe as a giant, invisible fabric. In physics, we often try to understand how this fabric can twist, knot, or form stable shapes called solitons (think of them like stable whirlpools in a river that don't just disappear).

This paper is about a specific type of knot in a 2D world (a flat sheet) that involves some very complex math. Here is the story of what the authors discovered, explained without the heavy jargon.

1. The Setup: Two Types of Knots

The authors are studying a model called the Skyrme-Chern-Simons (SCS) model.

• The "Skyrme" part: Imagine a rubber sheet that can be twisted into a ball. Usually, this ball has a "ground state" (the simplest, lowest-energy shape) and "excited states" (twisted, higher-energy shapes).
• The "Chern-Simons" part: This is like adding a magical, invisible magnetic wind that blows through the sheet. In previous studies, this wind caused some weird effects, like the energy of the knot changing in strange ways depending on how fast it spins or how much electric charge it holds.

The authors wanted to see: Does this magical wind change the rules for the "excited" (twisted) knots, or do they behave normally?

2. The Problem: A Broken Map

To find these excited knots, the physicists usually use a "map" (a mathematical formula) to describe the shape of the knot.

• The Trap: When they tried to use the standard map (called a "constraint compliant parametrization") to find the excited knots, the map suddenly broke. It was like trying to drive a car across a bridge that disappears right in the middle. The math showed a "discontinuity"—a jump or a glitch—specifically in a variable they called $g(r)$.
• The Consequence: Because the map broke, they couldn't find the excited knots using the old method. It was as if the excited knots were hiding in a dimension the old map couldn't see.

3. The Solution: The "Lagrange Multiplier" Trick

To fix the broken map, the authors used a clever mathematical tool called the Lagrange Multiplier method.

• The Analogy: Imagine you are trying to tie a shoelace, but you keep getting the knot wrong because you are trying to force the lace into a specific shape. Instead, you put a "guide rail" (the Lagrange multiplier) next to the lace. You don't force the lace into the shape; you just let the guide rail gently push it into the right place while you tie the knot.
• The Result: By using this "guide rail," they bypassed the glitch. They successfully found the excited knots that were previously hidden.

4. What They Found: The "Excited" vs. "Fundamental" Knots

They discovered that these knots are labeled by a number, $p$.

• $p = 0$ (Fundamental): This is the calm, simple knot. It has the lowest energy. It's like a smooth, round ball.
• $p \neq 0$ (Excited): These are the twisted, complex knots.
  • The "Node" Count: The number $p$ tells you how many times the knot "wiggles" or crosses zero as you move from the center to the edge. Think of it like the number of bumps on a rollercoaster track.
  • The "Odd" vs. "Even" Surprise: They found that knots with an odd number of wiggles ($p=1, 3, 5$) behave very strangely. They can split into two different "branches" of solutions. It's like a fork in the road where you can go left or right, and both paths lead to valid knots, but they look very different.
  • The "Hidden" Knots: Some of these excited knots are actually "imposters." They look like they are using the full 5-dimensional magic, but they are actually just simpler 3-dimensional knots hiding inside. The math revealed that these "imposter" knots have a glitch in their shape (the $g(r)$ variable jumps suddenly), which is why the old map couldn't find them.

5. The Big Conclusion: The Rules Didn't Change

The most important question was: Does the magical wind (the Chern-Simons term) change the energy hierarchy?

• The Expectation: Maybe the wind would make the twisted knots ($p=1, 2...$) cheaper (lower energy) than the simple knot ($p=0$).
• The Reality: No. Even with the magical wind, the simple, calm knot ($p=0$) is still the cheapest and most stable. The excited knots always cost more energy.
• The Takeaway: The magical wind adds some interesting flavor and weird behaviors (like the split branches for odd numbers), but it doesn't flip the basic order of the universe. The "ground floor" is still the ground floor.

Summary in One Sentence

The authors used a clever mathematical "guide rail" to find hidden, twisted knots in a 2D universe that were previously invisible due to a mathematical glitch, and they discovered that while these knots have some quirky behaviors, the simplest knot remains the most stable and lowest-energy option.

Technical Summary

1. Problem Statement and Motivation

The paper investigates excited soliton solutions within a Skyrme–Chern-Simons (SCS) model in 2+1 dimensions.

• Context: Previous work (Ref. [12]) studied a "contracted" version of this model where an $O(5)$ Skyrme scalar was reduced to an effective $O(3)$ scalar. That study found that the SCS term induces unusual behaviors in the energy ($E$), angular momentum ($J$), and electric charge ($Q_e$) relationships, such as non-monotonic slopes. However, that study was limited to "fundamental" solutions ($p=0$) where the scalar field configuration was static in a specific angular variable.
• The Gap: The authors aim to study excited solutions ($p \neq 0$) in the full $O(5)$ Skyrme model gauged by an $SO(2) \times SO(2)$ connection.
• The Technical Obstacle: A major difficulty arises when using a "constraint-compliant parametrization" (where the scalar field magnitude is fixed to 1 via trigonometric functions). In this parametrization, the function $g(r)$ (an angular variable of the scalar field) exhibits a discontinuity in its asymptotic boundary conditions when the asymptotic values of the gauge potentials satisfy $|c_\infty| = |d_\infty|$. This discontinuity prevents standard numerical integration methods from finding excited solutions, particularly those that are "embedded" $O(3)$ solutions within the larger $O(5)$ space.

2. Methodology

To overcome the numerical challenges and explore the full solution space, the authors employ the following methods:

• Model Definition:

  • The system consists of two $SO(2)$ gauge fields ($A_\mu, B_\mu$) and an $O(5)$ Skyrme scalar field ($\theta^a, a=1..5$).
  • The Lagrangian includes Maxwell terms, quadratic and quartic Skyrme kinetic terms, a "pion mass" potential, and the Skyrme-Chern-Simons (SCS) density term ($\Omega_{SCS}$).
  • The SCS term is derived from a descent of Skyrme-Chern-Pontryagin densities and is gauge-invariant in the equations of motion, though not in the action itself.

• Ansatz and Symmetry:

  • The authors assume azimuthal symmetry in the static limit.
  • The gauge potentials are defined by radial functions $a(r), c(r)$ and $b(r), d(r)$ with winding numbers $n$ and $m$.
  • The scalar field is parametrized using functions $f(r), g(r), \psi=n\phi, \chi=m\phi$.
  • Excited Solutions Condition: Excited solutions require $|n| \neq |m|$ and $|n|, |m| \geq 2$.

• Numerical Strategy (The Key Innovation):

  • Problem with Constraint Parametrization: Using the standard parametrization ($\theta^a = (\sin f \sin g, \dots)$) leads to a discontinuity in $g(r)$ at infinity for certain parameter regimes, making numerical integration impossible for excited states.
  • Lagrange Multiplier Method: The authors switch to a non-constraint compliant parametrization using Cartesian-like functions $R(r), S(r), T(r)$ subject to the algebraic constraint $R^2 + S^2 + T^2 = 1$.
  • They introduce a Lagrange multiplier $\beta(r)$ to enforce this constraint.
  • This transforms the problem into a system of 7 ODEs plus 1 algebraic constraint, which is solved using the COLDAE software package (a collocation method for differential-algebraic equations).
  • A compactified radial coordinate $x = r/(1+r)$ is used to handle the infinite domain.

3. Key Contributions

1. Resolution of the Discontinuity: The paper demonstrates that the Lagrange multiplier method is strictly necessary to obtain excited solutions in this model. It bypasses the discontinuity in the function $g(r)$ that plagues the constraint-compliant parametrization.
2. Characterization of Excited States: The authors define excited solutions by an integer $p$, which corresponds to the number of nodes of the function $S(r)$ (or "steps" in $g(r)$) in the radial domain.
   • $p=0$: Fundamental solutions (effectively $O(3)$).
   • $p \neq 0$: Excited solutions (fully $O(5)$ or embedded $O(3)$).
3. Discovery of Branching Solutions: For odd values of $p$ (e.g., $p=1$), the numerical results reveal two distinct solution branches for the same set of asymptotic parameters ($c_\infty, d_\infty$):
   • Lower Branch: Fully $O(5)$ Skyrme scalar solutions.
   • Upper Branch: Embedded $O(3)$ solutions. These are "excited" because $p \neq 0$, yet they effectively reduce to $O(3)$ symmetry. Crucially, these embedded solutions cannot be computed using the constraint-compliant parametrization due to the jump in $g(r)$.

4. Key Results

• Energy Hierarchy: The energy $E$ generally increases with the excitation number $p$. The fundamental solution ($p=0$) always possesses the minimal energy for a given set of parameters.
• Non-Standard Charge Behavior:
  • The global charges (Energy $E$, Angular Momentum $J$, Electric Charges $Q_e, Q_g$) exhibit complex dependencies on the asymptotic gauge potentials ($c_\infty, d_\infty$).
  • Discontinuity in Individual Charges: The individual electric charges $Q_e$ and $Q_g$ are discontinuous at the critical line $|c_\infty| = |d_\infty|$.
  • Continuity of Total Charge: The total electric charge $Q_t = nQ_e + mQ_g$ remains continuous across this boundary, resolving the anomaly.
• Energy Landscape:
  • For $p > 0$, the energy surface $E(J, Q_t)$ is multivalued and non-monotonic.
  • Local energy minima occur at non-zero angular momentum and electric charge (rotating, charged solutions), unlike the $p=0$ case where the minimum is at $J=0, Q_t=0$.
• Role of the SCS Term: The presence of the SCS term does not alter the fundamental pattern of energy levels. The excited states ($p > 0$) remain higher in energy than the fundamental state ($p=0$). The SCS term introduces exotic features (like negative slopes in $E$ vs $J$ plots) but does not invert the stability hierarchy.

5. Significance

• Theoretical Advancement: This work resolves a technical impasse in studying Skyrme-Chern-Simons models, proving that excited states exist and are accessible only through specific numerical techniques (Lagrange multipliers) that avoid coordinate singularities.
• Prototype for Higher Dimensions: The 2+1 dimensional model serves as a prototype for understanding anomaly densities in 3+1 dimensions. The authors note that constructing anomaly densities for gauged Skyrme fields in 3+1 dimensions is an active area of research, and the insights gained here (regarding the necessity of full scalar multiplets and the behavior of excited states) are crucial for future 3+1 dimensional studies.
• Physical Implications: The discovery of excited embedded $O(3)$ solutions that are inaccessible via standard parametrization suggests that the solution space of gauged Skyrme models is richer and more complex than previously thought, containing states that are topologically distinct yet dynamically embedded in lower-dimensional subspaces.

In summary, the paper successfully constructs and analyzes excited solitons in a 2+1 dimensional SCS model, overcoming significant numerical hurdles to reveal a complex energy landscape where the SCS term induces non-standard charge behaviors but preserves the fundamental energy hierarchy.