1. The Problem
Unlike the Roman Deity Janus enshrined with two equal faces, in the calculus of variations most of the attention is given to the direct venue which emerges naturally when a scientific problem is ab initio introduced via the variational principle. However, in many scientific endeavors, in particular when the underlying problem is partially phenomenological, models are derived on the level of differential equations, thus raising a natural quest for an underlying Lagrangian structure which apart from “scientific beauty” would carry many advantages, with conservation laws being one of them. The ‘upstream’ path from Equations to the underlying variational structure is known as an inverse problem of calculus of variations. Unfortunately, our ability to tackle the inverse problem and thus deduce the underlying Lagrangian structure from a given dynamical system is far more modest in spite of the long-standing attempts, c.f., [
1,
2,
3,
4,
5], which date back to Helmholtz, as those efforts, rather than provide us with a definite algorithm to construct the Lagrangian, provide us the conditions for its existence or statements of equivalence rather than with a specific ‘how’ (notably, there is a far richer literature on Hamiltonian structures than on Lagrangians, though a more balanced approach may be found in texts on dispersive waves. c.f., [
6,
7,
8]).
In the present communication, we unfold the algorithm to determine the underlying Lagrangian of a particular dynamical system as it relates to the
equations
and a variety of their extensions. Equation (
1) were introduced by us some time ago, refs. [
9,
10] and for
beget the compaction, solitary wave with a finite span. In the particular case
, its compactons take a simple trigonometric form. For instance, when
with
H being the Heaviside function. However, besides compactons, the
equations had many other fascinating features; at the time the model was introduced, Lagrangian structure did not seem to be one of them. The apparent lack of a Lagrangian structure was the main critique directed at the
models by Cooper et al. [
11], who introduced an alternative, compacton-supporting model derived from an ab initio assumed Lagrangian. In what follows, for the
and
cases we shall construct the seemingly ’missing’ Lagrangian and thus, at least in part, refute their critique, but first let us note that to account for other physically relevant cases [
12], the
model was extended to the
setup [
10]
which conserves
where
. For
, Equation (
3) reduces to the
setup and with model in [
11] being a particular
case.
As shown in [
10], for
, i.e.,
, the
equations follow from a conventional Lagrangian which in terms of
reads
adding to the conservation roster the conservation of energy. However, since in the
equations
, in this setup only the semi-linear, Korteveg–deVries-like sub-cases are endowed with the Lagrangian structure (5).
Yet, insofar as the
and
cases are concerned, i.e., for the
and
cases, the claim in [
11] was too sweeping as it
refers only to the conventional Lagrangian structure (5). We have recently demonstrated [
13], that both setups
admit a Lagrangian structure which, however, is ‘hidden from sight’ and, as we shall see shortly, quite different from the conventional one (5). Consequently, the two conservation laws in (4) were appended with an additional, non-local, conservation law.
As we shall see shortly, in both and setups, the derivation of the non-standard Lagrangian was based on non-obvious steps/tricks. The present communication aims to re-derive the underlying Lagrangian of the setup via a more algorithmic approach, which exploits a certain formal affinity between the equations and a dynamical system describing the motion of a Newtonian mass-spring chain which is not only easily deduced from a Lagrangian setup, but is often ab initio formulated via its Lagrangian. Although tying two very different setups may seem at first to be yet another ad hoc trick, as we shall see the idea may be easily extended to other dynamical systems, and thus has potentially a much wider scope. To avoid misunderstanding, we note that we assume the considered problems to be stated on the whole line with solutions vanishing at infinity and thus satisfying the natural boundary conditions implied by the variational derivation.
Let us first briefly summarize the derivation in [
13]; see also [
10]. To this end, the
equations are rewritten as a Hamiltonian system
This step is deductive. The maneuver which leads to the breakthrough is to introduce a new variable
v via
with
understood in a pseudo-differential sense. In terms of
, we have
The next, final step consists of introducing a ‘potential’ variable
[
7,
8,
9], where
that casts Equation (
8) into a form which, as per ref. [
4], is akin to the Lagrangian
of what, prior to their map, were
equations.
Two remarks are in order:
(1) As noted, in the preceding discussion we have tacitly assumed that our system is given on a line with all data vanishing at infinity, thus the natural boundary conditions assumed tacitly in the variational derivation hold trivially. If non-trivial boundary conditions are at play, the issue may become more evolved. However, recalling the very recent work of Olver [
14] who further developed the old observation [
1,
2,
3,
4,
5] that one may modify the variational problem without altering the corresponding Euler–Lagrange equations by adding a null Lagrangian to the integrand, we note that since this modification changes the associated natural boundary conditions, it enables one to enlarge the range of boundary value problems akin to variational techniques. Thus, the core difficulty remains to find the underlying Lagrangian, whereas its extension to handle non-natural boundary conditions may be delegated to the unfolding of a proper null Lagrangian.
(2) Note that the last step of introducing a potential in (5) is well known and apart from its more recent use in dispersive systems [
6,
7,
8], it could be found in the classical textbooks in electrodynamics., c.f., [
15], or continuum [
16,
17].
With a Lagrangian and its invariant properties at hand, in addition to the conservation of mass and energy, we now also have conservation of the momentum
, which in terms of the original variables, turns into the conservation of
It is easily seen that though the presented approach was applied to a very particular problem it may be quite naturally extended. Thus, for instance, one may consider
where
is any ‘reasonable’ linear operator and
any nice function [
4]. For instance, in Ref. [
18] addressing synchronisation of oscillators
was assumed. Using the previous notation, in terms of
v we have
In terms of
,
, the underlying Lagrangian reads,
As another extension, consider the
second order in time system (c.f., Equation (
27))
In terms of
,
, its Lagrangian reads,
and
.
In the next section, we take a detour needed to tie our problem with the obvious Lagrangian structure of a Newtonian mass-spring lattice which will render a more natural path toward the deduction of a Lagrangian structure of the
equations. Yet another take on the lattice–Lagrangian relations is briefly discussed in the
Appendix A.
2. The Newtonian Chain
We start with the Hamiltonian of an unrelated mass-particle chain [
19,
20],
Let
,
,
,
, and
, then expansion of the potential
P yields
To this order of expansion and the standard discrete-continuum association, the ab initio available Lagrangian density becomes
with its equations of motion
Thus, if
and
, then after one differentiation the resulting PDE reads
However, we shall also need an alternative description. Let
and
v a new variable defined
via
then in terms of
v, the Hamiltonian (13) reads
which yields
Acting on (20) with
and
we obtain
Actually, without invoking
v, we may use y ab initio and rewrite the Hamiltonian (13) in terms of
y
with its Lagrangian density
which yields
Of course, since
, acting with
, turns Equation (
24) into (21).
Note that if, say,
and
, obtained after carrying the expansion of
to fourth order, Equation (
21) turns into a second order in time variant of the
equation
With its Lagrangian structure given via (
14). Both Equations (
18) and (
27) and their underlying Lagrangians provide a quasi-continuum rendition of the same lattice, with expansions centering around different nodal locations.
Bridging between the Two Systems
The Newtonian system is second order in time. To turn it into a structurally analogous, but of first order in time system, we modify the kinetic part in (23)
with the resulting Lagrangian becoming
Identifying
P with
G and
with
y, begets a Lagrangian form identical to (12), with a related dynamical system which is first order in time. Setting
yields Equation (
11). In particular,
yields the
setup with the resulting sought-after Lagrangian structure.
Finally, consider again the direct expansion of the lattice and its underlying standard Lagrangian (14). Assuming again that
,
and repeating the
association, we recover module normalizable coefficients, the Lagrangian (5) and the underlying Equation (
3);
, with
.
Thus, as with Equations (21) and (25), with each of the expansions being based on a different nodal location on the lattice, we obtain a different equation of motion and a different corresponding Lagrangian: one which is conventional and readily available, and the other hidden and subject of the present exploration.