  By determining the type of all stationary points of the Gibbs free energy
functional for layered superconductors in parallel magnetic fields, we
establish the classification of all solutions to coupled static sine-Gordon
equations for the phase differences with respect to their stability. We prove
that the only minimizers of the free energy are the Meissner solution (the
"vacuum" state) and soliton vortex-plane solutions [S. V. Kuplevakhsky, Phys.
Rev. B vol. 60, 7496 (1999); ibid. vol. 63, 054508 (2001); cond-mat/0202293].
They are the actual equilibrium field configurations. We present a topological
classification of these solutions. In contrast, previously proposed non-soliton
configurations ("isolated fluxons", "triangular Josephson-vortex lattices",
etc.) are absolutely unstable and unobservable: They are nothing but saddle
points of the Gibbs free-energy functional and are not even stationary points
of the Helmholtz free-energy functional (obtained from the former by a Legendre
transformation). (Physically, non-soliton configurations violate conservation
laws for the current and the flux.) The obtained results allow us to explain
dynamic stability of vortex planes, noticed in numerical simulations, and to
provide a unified interpretation of the available experimental data. We hope
that the paper will stimulate interest in the subject of specialists in
different fields of physics and in applied mathematics.
