Next Article in Journal
Scheduling of Multi-AGV Systems in Automated Electricity Meter Verification Workshops Based on an Improved Snake Optimization Algorithm
Next Article in Special Issue
A Necessary Optimality Condition on the Control of a Charged Particle
Previous Article in Journal
Time–Energy Uncertainty Relation for Neutrino Oscillations: Historical Development, Applications, and Future Prospects
Previous Article in Special Issue
Well-Posedness and Energy Decay Rates for a Timoshenko-Type System with Internal Time-Varying Delay in the Displacement
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Optimal Control Perspective on Classical and Quantum Physical Systems

by
Mauricio Contreras González
1,
Marcelo Villena
2 and
Roberto Ortiz Herrera
3,4,*
1
Departamento de Física, Facultad de Ciencias, Universidad Metropolitana de Ciencias de la Educación UMCE, Santiago 7760197, Chile
2
Departamento de Ingeniería Comercial, Universidad Técnica Federico Santa María, Santiago 7630000, Chile
3
Facultad de Ingeniería, Universidad Diego Portales, Santiago 8370191, Chile
4
Facultad de Ciencias Económicas y Administrativas FACEA, Universidad Católica de la Santísima Concepción, Concepción 4070129, Chile
*
Author to whom correspondence should be addressed.
Symmetry 2023, 15(11), 2033; https://doi.org/10.3390/sym15112033
Submission received: 22 July 2023 / Revised: 9 October 2023 / Accepted: 18 October 2023 / Published: 8 November 2023

Abstract

:
This paper analyzes classical and quantum physical systems from an optimal control perspective. Specifically, we explore whether their associated dynamics can correspond to an open- or closed-loop feedback evolution of a control problem. Firstly, for the classical regime, when it is viewed in terms of the theory of canonical transformations, we find that a closed-loop feedback problem can describe it. Secondly, for a quantum physical system, if one realizes that the Heisenberg commutation relations themselves can be considered constraints in a non-commutative space, then the momentum must depend on the position of any generic wave function. That implies the existence of a closed-loop strategy for the quantum case. Thus, closed-loop feedback is a natural phenomenon in the physical world. By way of completeness, we briefly review control theory and the classical mechanics of constrained systems and analyze some examples at the classical and quantum levels.

1. Introduction

Over the last decades, there has been increasing interest among physicists in applying their concepts and methods to economics, as seen in some classical texts [1,2,3,4,5,6,7,8] or exciting articles about quantum mechanics in finance [9,10,11,12,13,14,15,16]. In this paper, we want to do the exact opposite, that is, we want to apply ideas from the optimal control theory (usually used in finance and economics [17,18,19,20,21]) to interpret the classical and quantum physical world. Indeed, open- and closed-loop control problems are not in the toolbox of mathematical methods in physics; nevertheless, dynamic optimization with its optimal control theory is the cornerstone of modern economic analysis.
There is little in the literature on the possible relationships between optimal control theory and classical and quantum mechanics, despite the fact that they share certain similar mathematical structures. For example, the Pontryagin equations resemble the Hamilton equations, and the Bellman equation is analogous to the Hamilton–Jacobi Equation, which corresponds to the classical limit of quantum theory. Recently, some studies have appeared on how a control problem associated with an economic model can be interpreted as a second-class constrained physical system [22,23,24,25,26]. In [22], it is found that at the classical level, the constrained dynamics given by Dirac’s brackets are the same as those provided by the Pontryagin equations [27]. The quantization of this second-class constrained system (with a right-order operator scheme) ended with a Schrödinger equation, which is just the Hamilton–Jacobi–Bellman equation used in optimal control theory and economics [28]. It is surprising to find that after the quantization of the Pontryagin theorem, Bellman’s maximum principle is obtained.
Thus, if an economic feedback system (characterized by an optimal control theory model) can be seen as a physical system, one can naturally ask the inverse question: Can the physical systems be seen from an optimal control perspective? Or, more specifically, do there exist open-loop and closed-loop strategies in physical systems? If the answer is affirmative, then where are they? Moreover, how and why do these feedback problems appear? This paper gives some clues and answers to all of these questions.
From a control theory perspective, open-loop strategies for the Lagrangian multiplier characterize the Pontryagin theory, while the Hamilton–Jacobi–Bellman equation has intrinsically closed-loop strategies. Since Pontryagin’s approach is equivalent to a classical mechanical model and Bellman’s theory is a quantum mechanical one, that suggests that open-loop strategies can be related to classical mechanics and closed-loop feedbacks with quantum mechanics.
This article shows that classical mechanics can be represented as an open-loop feedback problem in terms of Hamilton’s equations and a closed-loop strategies problem when one uses the canonical transformations approach. Indeed, both schemes are equivalent due to the invariance of the classical Hamiltonian equations by canonical transformations.
In the same way, from the perspective of quantum physics, the Heisenberg canonical commutation relations can be seen as a constraint in the non-commutative phase-space; this necessarily implies the existence of a specific relation between the momentum and the position, which is equivalent to a closed-loop feedback. Thus, open- and closed-loop feedback are also natural phenomena in our physical world.
This work has two principal objectives: The first is to give an optimal control perspective of the classical and quantum systems and to understand how and why such feedback appears. The second one is to show the physicist and optimal control community how the optimal control theory and physical systems (with its classical and quantum mechanical description of the world) are related. For that, it is necessary to have a common language and give the principal ideas behind each theory. The second objective will be fulfilled with the lecture in Section 2 and Section 3 of this paper. After that, we progress to complete the first objective.
It should be noted that given the length of our article, we have not addressed interesting recent topics of stability in control theory (see, for example, [29,30]) and quantum calculus (see, for example, [31,32]). Thus, the main purpose of the paper is to provide a general discussion of the relationship between classical mechanics and quantum mechanics with optimal control theory.
The paper is structured as follows: Section 2 reviews, for the non-specialist reader, some central ideas in control theory, such as the Pontryagin/Bellman approaches and the concepts of closed- and open-loop strategies. Section 3 establishes the basic ideas of classical mechanics in Lagrangian and Hamiltonian forms. Also, the so-called constrained systems in physics are defined and given as an example of this class of system, the optimal control theory itself, at classical and quantum levels.
In Section 4, we start with the paper’s primary objective, so Section 4.1 gives the first clues to interpret classical mechanics as a closed-loop system in control theory. Section 4.2 briefly reviews canonical transformations, which is the standard way to obtain the Hamilton–Jacobi Equation in physics from a purely classical perspective (non-quantum). In Section 4.3, we consider quantum mechanics and discuss why quantum mechanical systems must be interpreted naturally as closed-loop ones. Section 5 gives several examples of the presence of closed-loop strategies in physics at the classical and quantum levels. Lastly, in Section 6, the conclusions of this work are presented.

2. Key Ideas in Control Theory

This section reviews some critical concepts in control theory, such as the Pontryagin approach with its open-loop strategies and the Bellman approach with its closed-loop strategies, that we will use later to interpret the physical systems from an optimal control point of view.

2.1. Dynamic Optimization: The Pontryagin Approach

A generic one-dimensional system will be considered to emphasize the fundamental ideas and keep the equations simple. Generalizations to higher dimensions are straightforward. Consider a standard economic optimal control problem [17,18] for which one must optimize functional
A [ x , u ] = t 0 t 1 F ( x , u , t ) d t ,
where x represents a state variable (for example, the production rate of a specific good), and u denotes the control variable (such as the cost associated with a firm marketing). The state variable x must obey the market dynamics defined by the differential equation
x ˙ = f ( x , u , t ) x ( t 0 ) = x 0 .
and the optimal control problem consists of determining the production path x = x ( t ) and the control path u = u ( t ) , which maximise or minimise the cost function. To achieve this goal, one must apply the Lagrange multiplier method, so one considers instead an improved functional A defined over the extended space ( x , u , λ ) , which is given by:
A [ x , u , λ ] = t 0 t 1 F ( x , u , t ) λ ( x ˙ f ( x , u , t ) ) d t .
To obtain the optimal solution, the integrand of (3) can be interpreted as the Lagrangian:
L ( x , u , λ , x ˙ , u ˙ , λ ˙ ) = F ( x , u , t ) λ ( x ˙ f ( x , u , t ) ) .
so the extremal curves then satisfy the Euler–Lagrange equations:
L λ d d t L λ ˙ = 0 , L x d d t L x ˙ = 0 , L u d d t L u ˙ = 0 .
These equations are also written in the form
x ˙ = H λ , λ ˙ = H x , H u = 0 ,
with the Hamiltonian H defined by
H = H ( x , u , λ , t ) = F ( x , u , t ) + λ f ( x , u , t ) .
Equation (6) are the Pontryagin equations, which can be obtained through the Pontryagin maximum principle, and these solve the optimization problem.

2.2. Open-Loop and Closed-Loop Strategies

Note that the action (3) can be written in a Hamiltonian form as
A [ x , u , λ ] = t 0 t 1 λ x ˙ + H ( x , u , λ , t ) d t .
Now, the Pontryagin equations can also be obtained by optimizing the action (8) for its three variables x, u and λ . If δ x , δ u , and δ λ denote the corresponding functional variations, by expanding the Hamiltonian in a Taylor series around the optimal solution x * ( t ) , u * ( t ) , λ * ( t ) as
x ( t ) = x * ( t ) + δ x ( t ) , u ( t ) = u * ( t ) + δ u ( t ) λ ( t ) = λ * ( t ) + δ λ ( t ) ,
one obtains, after standard manipulations (that is, (1) by keeping only the first-order terms, (2) making an integration by parts and (3) noting that the initial point is fixed, i.e., δ x ( t 0 ) = 0 ).
δ A = t 0 t 1 H λ x ˙ δ λ + H x + λ ˙ δ x + H u δ u d t λ ( t 1 ) δ x ( t 1 ) .
where quantities inside the square bracket in the integral are evaluated at the optimal solution x * , u * , λ * . To maximize or minimize the action δ A must vanish, that is, the action must be stationary under arbitrary functional variations
δ A = 0 .
Now, in optimal control theory, there are two classes of control strategies:
  • Open-loop strategies that depend only on time: u = u ( t ) ;
  • Closed-loop strategies that depend on the state variable x and time: u = u ( x , t ) [33].
In the case of an open-loop strategy, x ( t ) , u ( t ) , and λ ( t ) are independent variables, so δ x , δ u , and δ λ are linearly independent. Hence, Equation (11) implies that coefficients that multiplies the functional variations in (10) must be vanish, and one obtains then the Pontryagin Equation (6) and the transversality condition λ ( t 1 ) = 0 for the optimal solution
x = x * ( t ) , λ = λ * ( t ) , u = u * ( t ) .
What happens, however, with closed-loop strategies? In this case, due to the relations between u = u ( x , t ) and x, the variations are related by δ u = u x δ x . Substituting this into (10) yields
δ A = t 0 t 1 H λ x ˙ δ λ + H x + H u u x + λ ˙ δ x d t λ ( t 1 ) δ x ( t 1 ) .
If λ and x are independent, the optimization of the functional (11) implies
H λ x ˙ = 0 , H x + H u u x + λ ˙ = 0 ,
plus the transversality condition. Note that the optimization process does not give the equation associated to the optimal control in (6). Then, u = u ( x , t ) must be given explicitly in terms of x, in another case, one has three unknowns but only two equations of motion. Thus, the variational problem is not well-defined for an arbitrary closed-loop strategy u = u ( x , t ) because the equations of motion can not determine its structure. It must be given from the beginning.
Now, note that equation for the control u in (6) is not a differential equation, but is algebraic:
H u = F ( x , u , t ) u + λ f ( x , u , t ) u = 0 .
From the above equation u * can be solved (in principle) as
u = u * ( x , λ , t ) ,
Thus, the optimization problem implies that the optimal open-loop u * strategy (12) is a closed-loop one. How can one understand this contradiction? Could the optimization problem be inconsistent?
To give an answer, consider the following arbitrary (non-necessarily optimal) closed-loop strategy u = u ( x , λ , t ) , so δ u = u x δ x + u λ δ λ . After substituting this into (10)
δ A = t 0 t 1 H λ + H u u λ x ˙ δ λ + H x + H u u x + λ ˙ δ x d t λ ( t 1 ) δ x ( t 1 ) .
If x ( t ) and λ ( t ) are independent, then (11) gives the equations
H λ + H u u λ x ˙ = 0 , H x + H u u x + λ ˙ = 0 .
Clearly the above equations are not equivalent to the Pontryagin open-loop equations. but if one choose for u its optimal value u * defined by H u = 0 , then Equation (18) becomes
H ( x , u * ( x , λ , t ) , λ , t ) λ x ˙ = 0 , H ( x , u * ( x , λ , t ) , λ , t ) x + λ ˙ = 0 ,
which are the Pontryagin equations. The solutions of (19) give optimal paths for the state variable x = x * ( t ) and Lagrangian multipliers λ = λ * ( t ) respectively. Then, the optimal control open-loop strategy u * ( t ) in (12) is obtained from the optimal closed-loop strategy (16) as
u = u * ( t ) = u * ( x * ( t ) , λ * ( t ) , t )
In this way, the optimal closed-loop strategy u = u * ( x , λ , t ) evaluated over the optimal x * ( t ) and λ * ( t ) trajectories is equivalent to the optimal open-loop u = u * ( t ) one. So, from now on, we will not distinguish between them in the case of Pontryagin’s theory.
Lastly, consider the case in which λ and x are not independent but λ = λ ( x , t ) , then δ λ = λ x δ x . Replacing it into (17), choosing the optimal control u = u * , considering the transversality condition λ ( t 1 ) = 0 and the fact that λ ˙ = λ x x ˙ + λ t , one obtains
δ A = t 0 t 1 H λ λ x x ˙ λ x + H x + λ x x ˙ + λ t δ x d t .
From H = F ( x , u * , t ) + λ f ( x , u * , t ) , it follows that
H λ = f ( x , u * , t ) H x = F ( x , u * , t ) x + λ f ( x , u * , t ) x ,
and using these equations in (21) finally gives
δ A = t 0 t 1 f λ x + F x + λ f x + λ t δ x d t
or
δ A = t 0 t 1 d H * d x + λ t δ x d t .
where
H * = H * ( x , t ) = H ( x , u * ( x , λ ( x , t ) , t ) , λ ( x , t ) , t )
is the effective Hamiltonian in terms of state variable x. Thus, the optimization process (11) implies that λ = λ ( x , t ) satisfies the optimal consistency condition
d H * ( x , t ) d x + λ t = 0 ,
If λ * ( x , t ) satisfies (23), the optimal state variable x ( t ) is obtained from Equation (2) as
x ˙ = f ( x , u * ( x , λ * ( x , t ) , t ) , t ) .
Thus, there exist three types of strategies in optimal control theory: open-loop ( x = x ( t ) , λ = λ ( t ) , u = u * ( t ) ), inert closed-loop ( x = x ( t ) , λ = λ ( t ) , u = u * ( x ( t ) , λ ( t ) , t ) ) and λ closed-loop ( x = x ( t ) , λ = λ ( x ( t ) , t ) , u = u * ( x ( t ) , λ ( x ( t ) , t ) , t ) ) strategies. The first two are entirely equivalent because these give the same dynamical equations.

2.3. Dynamic Optimization: Bellman theory

A different approach to the optimal control problem is given by the dynamic programming theory developed by Richard Bellman [28] in the 1950s. In this case, the fundamental variable is not the state variable but is the optimal value of the action defined by
J ( x 0 , t 0 ) = max u t 0 t F ( x , u , t ) d t ,
where x satisfies (2). Bellman’s dynamic programming principle implies that J ( x , t ) satisfies a non-linear partial differential equation called the Hamilton–Jacobi–Bellman equation [17] given by
max u F ( x , u , t ) + J ( x , t ) x f ( x , u , t ) = J ( x , t ) t .
Note that the Lagrangian multiplier λ of the Pontryagin approach can be identified with J ( x , t ) x . Thus, from the Pontryagin perspective, Bellman’s theory is a model with a closed-loop λ strategy λ ( x , t ) = J ( x , t ) x . Maximizing the left side of (26) for the control variable gives
u * = u * ( x , t ) = u * ( x , J ( x , t ) x , t ) = u * ( x , λ ( x , t ) , t )
so the Hamilton–Jacobi–Bellman equation becomes
F ( x , u * , t ) + J ( x , t ) x f ( x , u * , t ) = J ( x , t ) t .
Differentiating (27) with respect to x, one obtains
F ( x , u * , t ) x + 2 J ( x , t ) x 2 f ( x , u * , t ) + J ( x , t ) x f ( x , u * , t ) x + ( F ( x , u * , t ) u * + J ( x , t ) x f ( x , u * , t ) u * ) d u * ( x , t ) d x = 2 J ( x , t ) x t .
Noting that
H u = F ( x , u * , t ) u * + J ( x , t ) x f ( x , u * , t ) u * = 0 ,
and replacing J ( x , t ) x by λ ( x , t ) , one arrives to
F ( x , u * , t ) x + λ ( x , t ) x f ( x , u * , t ) + λ f ( x , u * , t ) x = λ ( x , t ) t ,
or
d H * ( x , t ) d x + λ ( x , t ) t = 0 .
The last equation is the same equation (23). Thus, the previous analysis implies that the consistency condition (23) is just the derivative of the Hamilton–Jacobi–Bellman equation. Equation (23) is then according to (27) equal to
d d x F ( x , u * , t ) + J ( x , t ) x f ( x , u * , t ) + J ( x , t ) t = 0 .
By integrating the last equation one obtains
F ( x , u * , t ) + J ( x , t ) x f ( x , u * , t ) + J ( x , t ) t = g ( t ) ,
where g ( t ) is a time-dependent arbitrary function. Then, the Pontryagin approach gives for a closed-loop λ * ( x , t ) = J ( x , t ) x a non-homogeneous Hamilton–Jacobi–Bellman equation, whereas Bellman’s maximum principle gives instead a homogeneous Hamilton–Jacobi–Bellman equation.

3. Classical Mechanics

In this section we review the Lagrangian and Hamiltonian approaches to classical mechanics and the analysis of constrained systems in physics.

3.1. Lagrangian Mechanics

It is well-known in classical mechanics that for a conservative system, that is, the one for which the force is the gradient of a potential energy U ( r )
F = U ( r ) ,
the Newton equations
F = m a
can be obtained through an optimization process similar to an optimization dynamic economic problem. For simplicity, in the case of a one-dimensional non-relativistic system given by the x coordinate, the objective functional is the Action functional A [ x , x ˙ ] defined by
A [ x , x ˙ ] = t 0 t 1 L ( x , x ˙ ) d t
where L ( x , x ˙ ) = 1 2 m x ˙ 2 U ( x ) is the Lagrangian function. The action (31) optimization for x gives the well-known Euler–Lagrange equation
L x d d t L x ˙ = 0
that, of course, is the Newton equation for the x variable:
L x ˙ = m x ˙ , d d t L x ˙ = x ¨ = a , L x = U x = F x
so Equation (32) is just F x = m a x . Note that the Lagrangian approach gives a second-order differential equation for the coordinate x.

3.2. Hamiltonian Mechanics

Instead of using the space ( x , x ˙ ) associated with the Lagrangian formalism, one can take the phase-space ( x , p x ) , where the momentum variable p x is defined by
p x = L ( x , x ˙ ) x ˙ .
The last equation can be used to solve x ˙ in terms of p x as x ˙ = x ˙ ( p x ) . In this phase-space ( x , p x ) , the objective functional is the Hamiltonian Action defined by
A [ x , p x ] = t 0 t 1 p x x ˙ H ( x , p x ) d t ,
where the Hamiltonian function H ( x , p x ) is constructed from the Lagrangian L ( x , x ˙ ) as
H ( x , p x ) = p x x ˙ ( p x ) L ( x , x ˙ ( p x ) ) .
The equations of motion in the phase-space ( x , p x ) are obtained by the optimization of the Action functional (35) for x and p x as independent variables. In this case, by the same standard techniques used in (10) one has that
δ A = t 0 t 1 H p x x ˙ δ p x + H x + p x ˙ δ x d t + p x ( t 2 ) δ x ( t 2 ) p x ( t 1 ) δ x ( t 1 )
To obtain the optimal phase-space trajectory one need that
δ A = 0 ,
so all the terms in δ x and δ p x must vanish. For this, it is supposed that:
  • The end points of x ( t ) are fixed, that is, δ x ( t 1 ) = 0 and δ x ( t 0 ) = 0 ;
  • δ x and δ p x are linearly independent.
Therefore, the Hamiltonian equations of motion are obtained from (38) as:
x ˙ = H ( x , p x ) p x , p ˙ x = H ( x , p x ) x .
One must recall that Hamilton equations are first-order differential equations, one for each independent x and p x variable. Both first-order equations are equivalent to the single second-order Euler–Lagrange Equation (32).
Now, by introducing the Poisson brackets between phase-space functions A ( x , p x ) and B ( x , p x ) according to
{ A , B } = A x B p x B x A p x
the Hamilton equations of motion read as
x ˙ = { x , H } p ˙ x = { p x , H } .
Note that the dynamical evolution of any variable θ defined on phase-space will given by Hamiltonian-type equation
θ ˙ = { θ , H }

3.3. Constrained Systems in Hamiltonian Mechanics

In the usual cases, the velocities x ˙ can be solved in terms of p x from Equation (34). However, there are situations (the most interesting ones) where this cannot be carried out. For example, for the Lagrangian
L ( x , x ˙ ) = a x ˙ + b .
one has that
p x = L x ˙ = a .
The class of systems for which the velocities cannot be solved regarding the momentum are called constrained systems. In this case, the Lagrangian and Hamiltonian theories can be different. Paul Maurice Dirac, a British physicist and Nobel winner, developed a methodology for studying general constrained systems in the phase-space. This method is now called Dirac’s method [34,35,36,37]. One must recall that constrained systems are ubiquitous in physics because all important physical theories, such as Electromagnetism, Einstein’s General Relativity, and the Gauge theories of particle physics, are examples of this category.
In the following subsection, Dirac’s method will be used to understand control theory from a physicist’s point of view.

3.4. Optimal Control Theory as an Example of a Constrained System

Recently, studies have shown how a generic control problem can be reinterpreted as a second-class constrained physical system [22,23,24,25,26]. In [22], it is found that at the classical level, the optimal control problem (1) and (2) when writing in terms of the state variable x, control variable u and the Lagrangian multiplier λ as (see Equations (7) and (8))
A [ x , u , λ ] = t 0 t 1 λ x ˙ + H ( x , u , λ , t ) d t = t 0 t 1 λ x ˙ + F ( x , u , t ) + λ f ( x , u , t ) d t
can be mapped to the Hamiltonian Action of a fictitious particle. In fact, by comparing Equation (45) with (35), one notes that the Lagrangian multiplier λ can be identified with momentum p x , so the functional (45) can be written as
A [ x , u , p x ] = t 0 t 1 p x x ˙ F ( x , u , t ) + p x f ( x , u , t ) d t
The above action is just a Hamiltonian action modulo a minus sign, so that one can use the action
A [ x , u , p x ] = t 0 t 1 p x x ˙ F ( x , u , t ) + p x f ( x , u , t ) d t
The Lagrangian associated with Equation (47) is
L = p x x ˙ F ( x , u , t ) + p x f ( x , u , t )
One must note that L does not depend on u ˙ , so its momentum p u associated with the u coordinate becomes
p u = L u ˙ = 0
Equation (49) is then a constraint in the phase-space ( x , u , p x , p u ) (one cannot solve u ˙ in terms of p u from the above equation), and the theory of optimal control thus becomes a constrained system. We denote the constraint in (49) as
Φ 1 = p u = 0
These constraints are called primary constraints according to Dirac’s method. The corresponding Hamiltonian in the phase-space ( x , u , p x , p u ) is then
H x , u , p x , p u = p x x ˙ + p u u ˙ L = p x x ˙ + p u u ˙ + F ( x , u , t ) + p x f ( x , u , t ) p x x ˙
so
H x , u , p x , p u , t = Φ 1 u ˙ + F ( x , u , t ) + p x f ( x , u , t )
or
H x , u , p x , p u , t = Φ 1 u ˙ + H 0 x , u , p x , p u , t
with
H 0 x , u , p x , p u , t = F ( x , u , t ) + p x f ( x , u , t ) .
On the constraint surface Φ 1 = 0 , the Hamiltonian is then
H x , u , p x , p u , t = H 0 x , u , p x , p u , t
Thus, H 0 is the energy of the system. The constraint (49) would be incorporated in the model by the Lagrange multipliers method, so one should consider instead a modified Hamiltonian defined by
H ˜ x , u , p x , p u , t = H 0 + μ 1 Φ 1 = H 0 + μ 1 p u
where μ 1 is a Lagrange multiplier. Now, one can impose that the constraint Φ 1 be preserved in time with the modified Hamiltonian (56), so
Φ ˙ 1 = Φ 1 , H ˜ = 0 ,
where { A , B } denotes the Poison bracket in the x , u , p x , p u phase-space
{ A , B } = A x B p x B x A p x + A u B p u B u A p u
Equation (57) yields
Φ ˙ 1 = p u , H ˜ = H ˜ u = H 0 u = F ( x , u , t ) u + p x f ( x , u , t ) u = 0 .
From the last equation, one cannot solve any velocity x ˙ or u ˙ in terms of momentum p x or p u . Thus, (57) is a new secondary constraint and is just the optimal law for the control variable. Then, from a phase-space analysis, the third Pontryagin equation in (6) is just the secondary constraint
Φ 2 = H 0 u = F ( x , u , t ) u + p x f ( x , u , t ) u = 0 .
To add the effect of the secondary constraint, one again uses the Lagrange multipliers method and considers a new extended Hamiltonian
H ˜ 2 x , u , p x , p u = H 0 + μ 1 Φ 1 + μ 2 Φ 2 .
According to Dirac’s method, one must impose time preservation for both constraints Φ 1 and Φ 2 with the extended Hamiltonian H ˜ 2 :
Φ ˙ 1 = Φ 1 , H ˜ 2 = 0 , Φ ˙ 2 = Φ 2 , H ˜ 2 = 0 ,
that is
Φ 1 , H 0 + μ 2 Φ 1 , Φ 2 = 0 Φ 2 , H 0 + μ 1 Φ 2 , Φ 1 = 0 ,
The above system can be written as
Φ 1 , H 0 Φ 2 , H 0 + 0 Φ 1 , Φ 2 Φ 1 , Φ 2 0 μ 1 μ 2 = 0 0 .
The matrix
Δ = 0 Φ 1 , Φ 2 Φ 1 , Φ 2 0
is called Dirac matrix. If the Dirac matrix is invertible, the above equation system can be solved for μ 1 and μ 2 , and in this situation, no new constraint appears in the theory.
Now
Φ 1 , Φ 2 = p u , Φ 2 = Φ 2 u = 2 H 0 u 2
so
Δ = 0 2 H 0 u 2 2 H 0 u 2 0 .
If
2 H 0 u 2 0
then Δ matrix is invertible, and no new constraint appears. In this case, the constraint set Φ 1 , Φ 2 is termed a second-class one (see [34,35,36,37] for details). Note that for a large class of economic systems, such as the linear quadratic ones, the condition (63) is fulfilled.
In summary, from a phase-space optics, the optimal control theory corresponds to a second-class constrained system, defined by the Hamiltonian H 0 = F ( x , u , t ) + p x f ( x , u , t ) and the two second-class constraints Φ 1 = p u = 0 and Φ 2 = H 0 u = F ( x , u , t ) u + p x f ( x , u , t ) u = 0 .
Now, one can go further by trying to quantize this classical system by imposing canonical commutations relations for the phase-space variables ( x , u , p x , p u ) :
[ x ^ , p ^ x ] = i I [ u ^ , p ^ u ] = i I
and where any other commutator is zero. Thus, the positions x, u and the momentum variables p x , p u become operators, which can be represented in a function space (Hilbert space) as
x ^ = x p ^ x = i x u ^ = u p ^ u = i u
The time-dependent Schrödinger equation associated to the Hamiltonian H 0 is
H ^ 0 x ^ , u ^ , p ^ x , p ^ u Ψ ( x , u , t ) = i t Ψ ( x , u , t )
or
F ( x , u , t ) i f ( x , u , t ) x Ψ ( x , u , t ) = i Ψ ( x , u , t ) t ,
Note that for the quantization process, a right-side operator order was chosen for the momentum operator.
Following Dirac, the only physically admissible solutions Ψ P of the Schrödinger Equation (66) satisfy the constraints Φ 1 = p u = 0 and Φ 2 = H 0 u = F ( x , u , t ) u + p x f ( x , u , t ) u = 0 at quantum level, that is, as operator equations acting on the wave function Ψ P
Φ ^ 1 Ψ P = 0 , Φ ^ 2 Ψ P = 0 .
or explicitly
i u Ψ P = 0
F ( x , u , t ) u i f ( x , u , t ) u x Ψ P = 0
and also the Schrödinger equation
F ( x , u , t ) i f ( x , u , t ) x Ψ P = i t Ψ P .
Equation (67) implies that Ψ P does not depend on u, so Ψ P = Ψ P ( x , t ) . Now, if one write the wave function in the form
Ψ P ( x , t ) = e i S ( x , t )
the quantum Equations (68) and (69) becomes the following equations for S ( x , t )
F ( x , u , t ) u + f ( x , u , t ) u S ( x , t ) x = 0 ,
F ( x , u , t ) + f ( x , u , t ) S ( x , t ) x = S ( x , t ) t
Note that the above two equations are -independent, due to the linear character of the Hamiltonian operator H 0 and the quantum constraint (68), in the momentum operator p ^ x = i x . Also, Equations (70) and (71) are equivalent to a single Hamilton–Jacobi–Bellman Equation (26)
max u F ( x , u , t ) + f ( x , u , t ) S ( x , t ) x = S ( x , t ) t ,
if the S ( x , t ) is identified with the optimal value function J ( x , t ) [28]. Thus, the Bellman approach to the optimal control problem is the right quantization of the Pontryagin theory. For more details about this topic, the elimination of the second-class constraints at the classical level, different quantization schemes, and issues related to the quantization process, see [22,25].

4. Open/Closed-Loop Strategies and Physics

The past two sections deal with control theory and the physicist’s view of it as a second-class constrained system in the phase-space. So after this digression, we enter the principal objective of this paper, that is, to invert the optics of the problem: instead of seeing the control theory as a physical system, one can consider the classical and quantum physical systems from the point of view of control theory. We will start with an analysis of the classical theory.

4.1. Open/Closed-Loop Strategies and Classical Mechanics

From an economic or dynamic optimization point of view, the problem of optimizing the Action (35) is analogous to an optimal control problem but without the control variable u. As seen in the previous section, the solutions to the control problem are given by the Pontryagin Equation (6) plus the transversality condition. Note that the first two Pontryagin equations in (6) are precisely the Hamiltonian equations of motion (39) if one identifies the Lagrangian multiplier λ with the canonical momentum p x . Because Pontryagin theory has open-loop strategies ( x ( t ) and λ ( t ) ) are independent), the Hamiltonian theory can then be considered as an open-loop model, similar to Pontryagin’s theory.
One may ask: can closed-loop strategies occur in physical systems as in economic systems?
As is shown in this paper, the answer is affirmative, and they appear naturally in the context of canonical transformations. In order to give a first clue to the answer to the above question, suppose that one imposes a constraint over the phase-space of the form
Φ ( x , p x , t ) = 0 .
This constraint represents at each time a surface in the phase-space (a line in our bi-dimensional case) where the system can evolve (see Figure 1).
Actually, one can write p in terms of x from (72) by solving the constraint:
p x = p x ( x , t ) .
Using the analogy of the momentum p with the Lagrangian multiplier λ , Equation (73) corresponds to the closed-loop λ strategy from an economic point of view. Thus, x and p are not independent in this case, and their variations are related by δ p x = p x x δ x . By replacing in (37) and using the fact that the endpoints of x are fixed, one arrives at
δ A = t 0 t 1 p x t H x H p x p x x δ x d t .
By defining the reduced Hamiltonian H * ( x , t ) = H ( x , p x ( x , t ) , t ) , Equation (74) can be written as
δ A = t 0 t 1 p x t H * x δ x d t .
The optimization of the action then gives
p x t H * x = 0 .
This last equation is a consistency condition that the closed-loop strategy (73) must satisfy, to give an extremal of the action. That is, if p x = p x ( x , t ) satisfies (76), then p x = p x ( x , t ) is an optimal closed-loop strategy.
Now, if the closed-loop momentum strategy p x = p x ( x , t ) is just the derivative of some function S ( x , t ) , such as
p x ( x , t ) = S ( x , t ) x ,
condition (76) gives
x S ( x , t ) t H * ( x , t ) = 0 .
So, one obtains, by integration, that
S ( x , t ) t H ( x , S ( x , t ) x ) = g ( t )
for some function g of time. The above equation is just an inhomogeneous Hamilton–Jacobi Equation. Thus, the derivative of the Hamilton–Jacobi Equation can be seen as the consistency condition to give the Action an extremal in the closed-loop p x strategy case. Also, this small analysis implies that closed-loop momentum strategies are closely related to the Hamilton–Jacobi Equation.
Now, it is well known that the Hamilton–Jacobi Equation appears in classical mechanics in the context of the canonical transformations, but how and where do the closed-loop p x = p x ( x , t ) strategies appear there? In the next section, a short review of canonical transformations will be given, and it will elucidate how the closed-loop strategies appear in that context.

4.2. Canonical Transformations and Closed-Loop Strategies

The Hamilton–Jacobi Equation can be considered as a bridge between nature’s classical and quantum descriptions. The Hamilton–Jacobi Equation is obtained from the quantum sector, taking the limit 0 of a particular form of the Schrödinger equation called the quantum Hamilton–Jacobi Equation (see Section 5 for details). Instead, from the pure classical mechanical description, the Hamilton–Jacobi Equation can be obtained using the canonical transformation framework. Canonical transformations are coordinate transformations in the Hamiltonian mechanics’ phase-space ( x , p x ) . Thus, two observers are present in this framework, each associated with a coordinate system. The first one is the standard Cartesian system. But, the second one is special: it is constructed with the trajectories that are solutions of the same Hamiltonian equations, i.e., the solutions of the Newton equations. Thus, one hopes from the point of view of the second observer that the solution of the equation of motion is described by constant values.
To understand how the canonical transformations method works, consider a general phase-space coordinate transformation of the form
Q = Q ( x , p x ) P = P ( x , p x ) .
The transformation is called canonical if the Hamiltonian equations of motions are invariant under (80), that is, if
Q ˙ = H ˜ P , P ˙ = H ˜ Q ,
where
H ˜ ( P , Q , t ) = H ( x , p x , t ) + F ( x , Q , t ) t .
The function F is called the generator of the canonical transformation, and the coordinate transformation (80) can be reconstructed from F through Equations [38,39]
p x = F ( x , Q , t ) x , P = F ( x , Q , t ) Q .
One must note at this point that the Hamiltonian Equation (39) refers to a unique coordinate system (in this case, a Cartesian coordinate system). So, the Hamilton equations in (39) are “single observer” equations. Instead, the canonical transformation brings a new second observer into the problem, because one has two different coordinate systems: the initial Cartesian ( x , p x ) and the second one ( Q , P ) . So, the theory of canonical transformations is a “two observers” view of classical mechanics. This characteristic induces the closed-loop p x -strategies from a purely classical point of view (closed-loop p x -strategies can also be induced from quantum mechanics to the classical realm, as we shall see later).
The Hamilton–Jacobi theory relies on considerable freedom in choosing F ( x , Q , t ) . This theory does not work directly with F, but with its Legendre transformation S defined by
S ( x , P , t ) = F ( x , Q , t ) + P Q
In this case, the canonical transformation is reconstructed via the equations
p x = S ( x , P , t ) x
Q = S ( x , P , t ) P
and the respective Hamiltonians are related by
H ˜ ( P , Q , t ) = H ( x , p x , t ) + S ( x , P , t ) t
Note again that one has huge freedom in choosing S. The Hamilton–Jacobi theory corresponds to the choice of S that makes the second-observer Hamiltonian H ˜ ( Q , P ) equal zero:
H ˜ ( P , Q , t ) = 0 .
So, the equations of motion for the second observer are
Q ˙ = 0 , P ˙ = 0 .
Thus, for the second observer, the dynamical variable remains constant in time:
Q ( t ) = Q 0 , P ( t ) = P 0 .
Note that this is possible only if the second coordinate system ( P , Q ) corresponds to the trajectories of the solution of the Hamilton equations. Only in this way can the solutions of the equations of motion have constant values. However, what does the first observer see? First, due to Equation (90), the coordinate transformations (80) give
Q 0 = Q ( x , p , t ) , P 0 = P ( x , p , t ) ,
but each of these equations defines constant ( Q , P ) coordinate lines. These are constraints over the phase-space ( p x , x ) of the first observer, from which one can generate two different closed-loop p x -strategies according to Equation (72) (See Figure 2).
Thus, the “two observers” perspective of classical mechanics, through the method of canonical transformation, is responsible for the generation of the closed-loop p x -strategies. From (91), it is not clear if the p x closed-loops strategies thus generated satisfy the consistency condition (76) or if they satisfy the second condition (77) to obtain a Hamilton–Jacobi Equation as in (79) for S. Instead of using (91), one can see these constant coordinate lines in terms of Equations (85) and (86). These equations are equivalent to (91), because the canonical transformation can be reconstructed from (85) and (86).
A constant P 0 -line in (91) is equivalent (from (85)) to
p x = S ( x , P 0 , t ) x ,
thus, the constant P 0 -line in (91) satisfies (77). From (87) and (88), S ( x , t ) satisfies the Hamilton–Jacobi Equation as in (79) with g ( t ) = 0 . This implies that a closed-loop strategy generated by (92) satisfies the consistency relation (76), so (92) defines a true optimal closed-loop p x strategy.
For the Q 0 constant coordinate line in (91), one has, however, due to (86), that
Q 0 = S ( x , P , t ) P ,
but, from this equation, one cannot obtain p x in terms of x, so this line does not generate a p x closed-loop strategy at all. This is due to the structure of Equations (85) and (86). In order to reconstruct the canonical transformation (80), it is necessary to invert the system (85) and (86). Note that only from (85) can the momentum P be written in terms of the first-observer variables as P = P ( x , p x , t ) . But, from (86) alone, one cannot solve Q in terms of the x , p x . The other Equation (85) is needed to do that. Thus, a constant Q 0 -line alone can not generate a true closed-loop p x = p x ( x , t ) strategy.
In this way, closed-loop p x strategies appear in classical mechanics due to the two observers’ interpretation of the canonical transformation theory. These closed-loop strategies are inert, similarly to how the optimal u * closed-loop ones are inert in control theory. That is because both closed-loop approaches p x and u * give the same dynamical equation of the open-loop case. For the closed-loop p x case, the open-loop dynamics (analogous to those given by Pontryagin’s equations) are provided by the Hamiltonian equations of motion of the first observer in the ( x , p x ) phase-space.

4.3. Quantum Mechanics and Closed-Loop Strategies

In this section, the origins of the Hamilton–Jacobi–Bellman equation that appears in the limit 0 in the quantum phenomena will be explained as a consequence of the emergence of closed-loop p x -strategies in the quantum world.
Consider the Schrödinger equation for a non-relativistic particle of mass m:
2 2 m 2 Ψ ( x , t ) x 2 + U ( x ) Ψ ( x , t ) = i Ψ ( x , t ) t .
Writing the wave function in the form
Ψ ( x , t ) = e i S ( x , t ) ,
and by substituting (95) into the Schrödinger equation, the following equation, called the quantum Hamilton–Jacobi Equation for S ( x , t ) , is obtained:
1 2 m S ( x , t ) x 2 + U ( x ) i 2 m 2 S ( x , t ) x 2 = S ( x , t ) t .
Note that this equation is completely equivalent to Schrödinger’s equation, but here, the classical and quantum realms can be clearly identified. In fact, by taking the limit 0 in (96), one obtains
1 2 m S ( x , t ) x 2 + U ( x ) = S ( x , t ) t .
Equation (97) is just the classical Hamilton–Jacobi Equation
H ( x , S ( x , t ) x ) = S ( x , t ) t ,
associated with the classical Hamiltonian function of the non-relativistic particle
H ( x , p x ) = p x 2 2 m + U ( x ) ,
where one must identify p x with the derivative of S
p x = p x ( x , t ) = S ( x , t ) x .
to make contact with the classical Hamiltonian theory. And, it is precisely this identification which generates the closed-loop p x strategy through (100). Note that it is induced from the quantum realm to the classical world, in the limit 0 , of the quantum Hamilton–Jacobi Equation. The identification in (100) thus comes from a pure quantum description. Consider the momentum operator
p ^ x = i x .
This operator is characterized by its eigenfunctions and eigenvalues:
p ^ x Φ p x ( x ) = p x Φ p x ( x ) .
where the solution of this equation gives
Φ p x ( x ) = e i p x x .
In this context, p x and x are independent variables, and the eigenfunction (103) corresponds to states with well-defined values of the momentum.
Note now that if one applies the momentum operator to a generic wave function Ψ , which is a solution of the Schrödinger equation (written in the “momentum form” (95)), one obtains
p ^ x Ψ ( x , t ) = S ( x , t ) x Ψ ( x , t ) .
By looking at the wave function as a vector with a continuous index x, the above equation implies that (locally at each point x) the momentum operator is diagonal so that any wave function can be seen as an eigenstate of the momentum operator with momentum eigenvalue S ( x , t ) x . Thus, one must identify the momentum eigenvalue p x in this quantum state with the derivative of the S function through (100). This identification generates the closed-loop p x -strategies directly in the quantum world.
On the other hand, the same Heisenberg canonical commutation relations
[ x ^ , P ^ x ] = x ^ P ^ x P ^ x x ^ = i I ^ .
can be seen as a constraint in the non-commutative phase-space ( x ^ , P ^ x ) . Thus, from (105) one could “solve” the momentum operator P ^ x in terms of the x ^ operator. This necessarily implies the existence of a certain relation between P ^ x and x ^ or between their eigenvalues. The representation of the canonical operator as a differential operator acting on a function space or Hilbert space as
x ^ x , p ^ x i x
is equivalent to solving the constraint (105), because on any wave function Ψ ( x , t ) , Equation (105) is satisfied identically. The memory of the quantum constraint (105) is then transferred locally to the momentum eigenvalue, according to (104). In a sense, the representation of the wave function as Ψ ( x , t ) = e i S ( x , t ) locally diagonalizes the momentum operator p ^ x over any quantum state, and (96) is just the Schrödinger equation in this diagonal basis. Note that all of this is a kinematic effect created by the Heisenberg commutation relation (105); the dynamical effects appear when the explicit form of S ( x , t ) is needed, and for that, one must solve the full quantum Hamilton–Jacobi Equation (96) explicitly. Note that quantum mechanics, as in (105), can be viewed as a constrained system in a non-commutative space, so one would apply a generalization of Dirac’s method [34,35,36,37] to non-commutative spaces [40] to study quantum mechanical systems.
We can say that closed-loop p x -strategies correspond to a pure quantum phenomenon and are a consequence of Heisenberg’s uncertainty principle. In an arbitrary quantum state, momentum and position cannot be independent; they are related through the non-commutative character of the position and momentum operators defined by the canonical commutation relation (105). A less-defined position state would emerge in a more-defined momentum state. Thus, these two variables must depend on one another in some way. Relation (100) is tantamount to a conversation between them. Only in a pure-momentum state, as given in (103), the link disappears and position and momentum become independent variables.
In fact, in a pure-momentum state, p x ( x , t ) = p x 0 is constant, that is, all of the eigenvalues are the same, so Equation (100) gives
S ( x , t ) = p x 0 x + ϕ ( t )
as a solution, where ϕ ( t ) is some function of time. Thus, the wave function is
Ψ ( x , t ) = e i S ( x , t ) = e i ( p x 0 x + ϕ ( t ) ) = e i ϕ ( t ) Φ p x 0 ( x ) ,
which is the same momentum eigenstate amplified by a temporal arbitrary phase. Then, the linear character of S ( x , t ) in terms of x implies that p x and x are independent variables, and no closed-loop p x strategy exists in this case. The same can be said for a pure-position eigenstate.
Thus, closed-loop p x -strategies are an inherent part of the quantum mechanical world and permeate the classical world in the limit 0 through the Hamilton–Jacobi Equation.

5. Some Examples

In the following Section 5.1 and Section 5.2, we analyze some common textbook examples from closed-loop strategies’ point of view. Section 5.3 gives a quantum mechanical example, and Section 5.4 and Section 5.5 give some non-canonical quantum examples that illustrate the dependence of closed-loop strategy on the explicit form of the canonical commutations relations, that is, the dependence on the quantum constraint. Finally, in Section 5.5 we give a quantum control example.

5.1. The Stationary Case

The quantum and classical Hamilton–Jacobi Equations (96) and (97) are non-stationary equations, that is, they depend explicitly on time. In quantum mechanics, stationary states play a fundamental role. They are defined by
Ψ ( x , t ) = e i E t Φ ( x ) .
By substituting this into the time-dependent Schrödinger Equation (94), the time-independent or stationary Schrödinger equation is obtained:
2 2 m 2 Ψ ( x ) x 2 + U ( x ) Ψ ( x ) = E Ψ ( x ) .
Now, by writing
Ψ ( x ) = e i W ( x )
and substituting into (110), the stationary quantum Hamilton–Jacobi Equation holds:
1 2 m W ( x ) x 2 + U ( x ) i 2 m 2 W ( x ) x 2 = E .
Taking again the classical limit 0 in (112), the stationary classical Hamilton–Jacobi Equation appears:
1 2 m W ( x ) x 2 + U ( x ) = E .
Due to (109) and (111), we have Ψ ( x , t ) = e i S ( x , t ) = e i ( W ( x ) E t ) , which implies that
S ( x , t ) = W ( x ) E t
for the stationary case. In this case, the closed-loop p x -strategies are given by
p x = p x ( x , t ) = W ( x ) x .

5.2. The Non-Stationary Case

But what about the non-stationary closed-loop p x -strategies in the classical limit? In order to analyze this case, consider the example of a free particle, that is, U ( x ) = 0 . The non-stationary classical Hamilton–Jacobi Equation is now
1 2 m S ( x , t ) x 2 = S ( x , t ) t .
One can find a solution of the form S ( x , t ) = 1 2 a ( t ) x 2 , so by substituting in (115), one obtains a ( t ) = 1 ( P 0 t m ) , so S ( x , t ) = 1 2 x 2 ( P 0 t m ) , and the corresponding closed-loop p x strategy is
p x ( x , t ) = S ( x , t ) x = x ( P 0 t m ) .
One can evaluate x ( t ) using Equation (93)
Q 0 = S ( x , P 0 , t ) P 0 ,
where the integration constant P 0 must be identified with the constant momentum for the second observer in the coordinate system ( Q , P ) . Thus
Q 0 = 1 2 x 2 ( P 0 t m ) 2 ,
from which x ( t ) is computed as
x ( t ) = 2 P 0 t m 2 Q 0 = 2 Q 0 P 0 t m .
The associated open-loop p x ( t ) strategy is found by p x ( t ) = p x ( x ( t ) , t ) , similarly to Equation (20). Thus, by substituting (119) into (116):
p x ( t ) = x ( t ) ( P 0 t m ) = 2 Q 0 ( P 0 t m ) ( P 0 t m ) = 2 Q 0 = p 0 ,
so
x ( t ) = p 0 P 0 t m = p 0 t m + x 0 ,
where x 0 = P 0 p 0 . These last equations are the solutions for the motion of a free particle of course!
Now from the Hamiltonian equations, one obtains the open-loop dynamics for the free particle:
x ˙ = p x m , p ˙ x = U ( x ) x = 0 ,
so the open-loops dynamics are
x ( t ) = p 0 m t + x 0 , p x ( t ) = p 0 .
Then, the open-loop p x ( t ) strategy coming from the Hamiltonian equations of motion is equivalent to the non-stationary closed-loop p x ( x , t ) strategy coming from the non-stationary classical Hamilton–Jacobi Equation. This equivalence is valid because the classical Hamilton–Jacobi Equation approach corresponds to a “two-observer” point of view of classical mechanics. The S ( x , t ) function is just the generator of the canonical transformation, which leaves the Hamiltonian equations invariant. Thus, there are two schemes:
  • The Hamiltonian “one observer” approach with its open-loop p x ( t ) -strategies;
  • The Hamilton–Jacobi “two observer” approach with its closed-loop p x ( x , t ) -strategies.
Closed-loop strategies p x ( x , t ) coming from the Hamilton–Jacobi Equation are similar to the inert optimal closed-loop strategies u * ( x , t ) of Pontryagin’s approach (in optimal control theory) in the sense that they are equivalent to the open-loop ones.
Note that this equivalence is not generally valid for control theory when analyzed in the phase-space, as conducted in Section 3.4. As is shown in [37,41], for second-class systems, the description of a mechanical system in terms of canonical transformations (together with the Hamilton–Jacobi Equation for the S function) can be inconsistent. It is due that the second-class systems violate the Carathéodory’s integrability conditions. By using the methodology developed in [41], in [26], there is proof that there exists a class of optimal controls problems, called regular subclass, in which the functions F ( x , u , t ) and f ( x , u , t ) in (1) and (2) have the form
F ( x , u , t ) = F 0 ( x ) + F 1 u + 1 2 F 2 u 2 , f ( x , u , t ) = f 0 ( x ) + f 1 u ,
(where F 1 , F 2 and f 1 are constants and F 0 ( x ) and f 0 ( x ) are arbitrary continuous functions), such that for this subclass:
(i)
The Hamilton–Jacobi Equation is well-defined;
(ii)
The solution to the Hamilton–Jacobi Equation gives the same dynamics as the Pontryagin equations.
For more general forms of the F ( x , u , t ) and f ( x , u , t ) functions, it is not clear whether the problem is consistent or if the dynamics of the Hamilton–Jacobi–Bellman are equivalent to the Pontryagin equations in the phase-space.

5.3. The Pure Quantum Limit and Closed-Loop Strategies

In the previous section, the classical limit 0 was taken, and its characteristics were explored in terms of the closed-loop strategies. In this section, the inverse limit will be taken, that is, > > 1 , and the consequences of a higher non-commutative system of quantum variables will be explored.
Consider the non-stationary quantum Hamilton–Jacobi Equation (96). Taking the limit > > 1 and supposing that the time derivative of the S function has a higher value,
i 2 m 2 S ( x , t ) x 2 = S ( x , t ) t .
or, what is the same,
2 2 m 2 S ( x , t ) x 2 = i S ( x , t ) t .
But, this is the free-particle Schrödinger equation for S ( x , t ) . In this way, one can write
S ( x , t ) = e i T ( x , t )
where T ( x , t ) satisfies the quantum Hamilton–Jacobi Equation
1 2 m T ( x , t ) x 2 i 2 m 2 T ( x , t ) x 2 = T ( x , t ) t
and the wave function is given by
Ψ ( x , t ) = e i S ( x , t ) = e i e i T ( x , t ) .
The corresponding closed-loop p x strategy is
p x = S ( x , t ) x = i T ( x , t ) x S ( x , t )
Note that T ( x , t ) x can be interpreted as a closed-loop p T strategy for the quantum Hamilton–Jacobi Equation (127). Thus, denoting the closed-loop p x strategy for (124) by p S ( x , t ) , then both strategies are related by
p S ( x , t ) = i S ( x , t ) p T ( x , t )
Again, if the time derivative of T has a higher value, and as > > 1 , the quantum Hamilton–Jacobi for T ( x , t ) (127) is in this limit again a Schrödinger equation,
2 2 m 2 T ( x , t ) x 2 = i T ( x , t ) t .
Hence, we can write
T ( x , t ) = e i U ( x , t )
where U ( x , t ) satisfies
1 2 m U ( x , t ) x 2 i 2 m 2 U ( x , t ) x 2 = U ( x , t ) t .
and the wave function is
Ψ ( x , t ) = e i S ( x , t ) = e i e i e i U ( x , t )
Putting p U = U ( x , t ) x the closed-loop strategy associated to the quantum Hamilton–Jacobi Equation for U (133), the corresponding closed-loop p S strategy is then in this case
p S ( x , t ) = i p U ( x , t ) S ( x , t ) T ( x , t ) .
Since one can keep iterating this procedure to infinity, quantum mechanical systems can admit multistage closed-loop strategies, and they are connected in a strongly non-linear way, as in (134).

5.4. A Non-Canonical Example

It is well known that quantum gravity effects, modelled by string theory, loop quantum gravity, or black hole physics, predict the existence of a Generalized Uncertainty Principle (GUP), which can change the usual canonical commutation relations [42,43,44,45,46,47,48] and its implications for entanglement and the Hamilton–Jacobi Equation [49,50,51,52,53]. In this case, the canonical commutation relations
[ x ^ , p ^ x ] = i I
are replaced by a more general one of the type
[ x ^ , p ^ x ] = i I + F ( x ^ , p ^ x )
for some function F ( x , p x ) . The form of the quantum algebra (137) guarantees that the system has a classical limit when goes to zero. Now, if one sees the commutation relations as a constraint in a noncommutative space, then if one “solves” p ^ x in terms of x ^ from (137), then an explicit form of p ^ x in terms of x ^ would depend on the function F. Thus, the expression (100) for p x ( x , t ) would depend on F. To explicitly show that dependence, we consider two simple toy examples of GUP algebras below.

5.4.1. A GUP Algebra Depending on x ^

As an example this dependence on GUP, consider the following commutation relations
[ x ^ , p ^ x ] = i I + α x ^
for which F ( x , p ) = α x , and where α is a constant. An operator representation of the above commutation relations is
x ^ = x
p ^ x = i x + α x x = i ( 1 + α x ) x
The momentum eigenstates are, in this case,
p ^ x Φ p x ( x ) = p x Φ p x ( x )
i ( 1 + α x ) Φ p x x = p x Φ p x
so
Φ p x ( x ) = C ( 1 + α x ) i p x α = C e i p x ln [ ( 1 + α x ) 1 α ] .
Consider now the usual non-relativist classical Hamiltonian (99). The quantization of this Hamiltonian by the rule (138), implies the following Schrödinger equation
2 2 m ( 1 + α x ) x ( 1 + α x ) x Ψ ( x , t ) + V ( x ) Ψ ( x , t ) = i Ψ ( x , t ) t
or
2 2 m α ( 1 + α x ) Ψ ( x , t ) x + ( 1 + α x ) 2 2 Ψ ( x , t ) x 2 + V ( x ) Ψ ( x , t ) = i Ψ ( x , t ) t
The corresponding Hamilton–Jacobi Equation is obtained from (95), by using the fact that
Ψ ( x , t ) x = e i S S x i
2 Ψ ( x , t ) x 2 = e i S S x 2 i 2 + i e i S 2 S x 2
so (145) becomes
2 2 m α ( 1 + α x ) e i S S x i + ( 1 + α x ) 2 e i S i 2 2 S x 2 + x 2 i e i S 2 S x 2
+ V ( x ) e i S = i e i S i S t
or
( 1 + α x ) 2 2 m S x 2 + V ( x ) i 2 m ( 1 + α x ) S x + ( 1 + α x ) 2 2 S x 2 = S t
By taking the classical limit 0 , one arrives to the modified Hamilton–Jacobi equation
( 1 + α x ) 2 2 m S x 2 + V ( x ) = S t .
Now, one must identify p x ( x , t ) , not with the derivative of S as in (100), but instead with
p x ( x , t ) = ( 1 + α x ) S ( x , t ) x .
Thus, the closed-loop p x strategy, depends on the form of the GUP. In fact, the momentum operator (140) acting on (95) gives
p ^ x Ψ ( x , t ) = i ( 1 + α x ) x e i S ( x , t ) = ( 1 + α x ) S x Ψ ( x , t )
so from the above equation, one can identify again the local momentum p x ( x , t ) again with ( 1 + α x ) S x .
If, instead of the commuting relation (138), one can consider the generalization
[ x ^ , p ^ x ] = i I + α x ^ n
for n = 1 , 2 , 3 , the momentum operator has the representation
p ^ x = i ( 1 + α x n ) x
so when acting on the wave function (95) it gives
p ^ x Ψ ( x , t ) = i ( 1 + α x n ) x e i S ( x , t ) = ( 1 + α x n ) S x Ψ ( x , t )
then the corresponding closed-loop p x ( x , t ) strategy is in this case
p x ( x , t ) = ( 1 + α x n ) S ( x , t ) x .

5.4.2. A GUP Algebra Depending on the Momentum p ^ x

Consider now the following GUP commutation relations depending on the momentum
[ x ^ , p ^ x ] = i I + β p ^ x
These commutation relations can be represented by the differential operators
x ^ = x
p ^ x = 1 β e β p ^ 0 I
where p ^ 0 = i x is the standard momentum operator for β = 0 case.
The action of the momentum operator (160) on the wave function (95) is
p ^ x Ψ ( x , t ) = 1 β e β p ^ 0 I e i S ( x , t )
By expanding the exponential operator e β p ^ 0 , one can show that
1 β e β p ^ 0 I e i S ( x , t ) = 1 β e β i x I e i S ( x , t ) =
1 β e β ( S x ) 1 + i β e β i x I β i x S ( x , t ) + i Θ ( x , t ) e i S ( x , t )
where the Θ function contains products between powers of ( S x ) and high-order derivatives of S ( x , t ) . For example, the first terms of Θ are
Θ ( x , t ) = 1 2 β 2 S x 2 S x 2 + 1 4 β 3 S x 2 2 S x 2 +
1 8 i β 3 2 S x 2 2 + 1 6 i β 3 S x 3 S x 3 +
so
p ^ x Ψ ( x , t ) = p ( x , t ) Ψ ( x , t )
where the local quantum momentum p ( x , t ) is
p ( x , t ) = 1 β e β ( S x ) 1 + i β e β i x I β i x S ( x , t ) + i Θ ( x , t )
Note that p ( x , t ) is a complex number in this case. In the classical limit 0 , the local quantum momentum (167) becomes
p ( x , t ) = 1 β e β ( S x ) 1
because the operator
i β e β i x I β i x
goes to zero when 0 . The classical local momentum given by (168) could also be obtained from the corresponding Hamilton–Jacobi Equation. The Schrödinger equation associated with the Hamiltonian function (99) is in this case
1 2 m β 2 e β p ^ 0 I e β p ^ 0 I Ψ + V ( x ) Ψ = i Ψ t
or
1 2 m β 2 e 2 β p ^ 0 2 e β p ^ 0 + I e i S ( x , t ) + V ( x ) e i S ( x , t ) = S ( x , t ) t e i S ( x , t )
Due to the fact that in the classical limit 0
e β p ^ 0 e i S ( x , t ) = e β i x e i S ( x , t ) = e β ( S x ) e i S ( x , t )
Equation (171) becomes in this classical limit
1 2 m β 2 e 2 β ( S x ) 2 e β ( S x ) + 1 e i S ( x , t ) + V ( x ) e i S ( x , t ) = S ( x , t ) t e i S ( x , t )
that is
1 2 m β 2 e β ( S x ) 1 2 + V ( x ) = S ( x , t ) t
From the above equation, one can identify the corresponding classical closed-loop strategy p ( x , t ) as the same as that given in (168). Thus, we see that the quantum commutation relations (canonical or GUP) imply the existence of closed-loop strategies for the local momentum, because these quantum commutation relations are constraints in a non-commutative space.

5.5. A Quantum Control Example

At last, we consider a simple example of quantum control [54,55] defined by a single spin particle S = 1 2 ( σ x , σ y , σ z ) , being σ x , σ y , and σ z , which are the so-called Pauli matrices. The Schrödinger equation for this particle in a magnetic field B = ( B x , B y , B z ) is [54]
ψ t = 1 2 σ z u z + σ x u x ( t ) + σ y u y ( t ) ψ ,
where the components B x , B y are considered as time-dependent control variables u x ( t ) , u y ( t ) , and the B z component is a fixed parameter u z .
By separating the two component spinor
ψ = ψ R + i ψ I = x 1 x 2 + i x 3 x 4
in its real and imaginary parts and doing the same with the Pauli matrices, the Schrödinger (175) can be writing as a four vector equation of the form [54]
d x d t = 1 2 H x u x ( t ) + H y u y ( t ) + H z u z x .
where
H x = 0 0 0 1 0 0 1 0 0 1 0 0 1 0 0 0 H y = 0 1 0 0 1 0 0 0 0 0 0 1 0 0 1 0 H z = 0 0 1 0 0 0 0 1 1 0 0 0 0 1 0 0
and
x = x 1 x 2 x 3 x 4
The quantum control problems consist thus in to minimize the functional [54]
J u x , u y = Φ ( x ( T ) ) + k 0 T u x 2 ( t ) + u y 2 ( t ) d t
subject to the Schrödinger Equation (176).

5.5.1. Open-Loop Pontryagin Dynamics

Introducing a Lagrange multiplier for each scalar equation in the system (176), i.e., the vector
λ = λ 1 λ 2 λ 3 λ 4
the Pontryagin Hamiltonian h is
h = k u x 2 ( t ) + u y 2 ( t ) + 1 2 λ t H x u x ( t ) + H y u y ( t ) + H z u z x
The corresponding Pontryagin equations are for the state variables (spinor components)
x ˙ = 1 2 H x u x ( t ) + H y u y ( t ) + H z u z x = 1 2 H x
and for the Lagrange multipliers
λ ˙ = 1 2 H x u x ( t ) + H y u y ( t ) + H z u z λ = 1 2 H λ
where
H = 0 u y u z u x u y 0 u x u z u z u x 0 u y u x u z u y 0
The equation for the controls gives
u x = 1 4 k λ t H x x = λ 1 ( t ) x 4 ( t ) λ 2 ( t ) x 3 ( t ) + λ 3 ( t ) x 2 ( t ) + λ 4 ( t ) x 1 ( t ) 4 k
u y = 1 4 k λ t H y x = λ 1 ( t ) x 2 ( t ) + λ 2 ( t ) x 1 ( t ) λ 3 ( t ) x 4 ( t ) + λ 4 ( t ) x 3 ( t ) 4 k
By replacing (182) and (183) in (181) one has that
H [ x , λ ] = 0 λ 1 ( t ) x 2 ( t ) λ 2 ( t ) x 1 ( t ) + λ 3 ( t ) x 4 ( t ) λ 4 ( t ) x 3 ( t ) 4 k λ 1 ( t ) x 2 ( t ) + λ 2 ( t ) x 1 ( t ) λ 3 ( t ) x 4 ( t ) + λ 4 ( t ) x 3 ( t ) 4 k 0 u z λ 1 ( t ) x 4 ( t ) λ 2 ( t ) x 3 ( t ) + λ 3 ( t ) x 2 ( t ) + λ 4 ( t ) x 1 ( t ) 4 k λ 1 ( t ) x 4 ( t ) λ 2 ( t ) x 3 ( t ) + λ 3 ( t ) x 2 ( t ) + λ 4 ( t ) x 1 ( t ) 4 k u z u z λ 1 ( t ) x 4 ( t ) + λ 2 ( t ) x 3 ( t ) λ 3 ( t ) x 2 ( t ) λ 4 ( t ) x 1 ( t ) 4 k λ 1 ( t ) x 4 ( t ) + λ 2 ( t ) x 3 ( t ) λ 3 ( t ) x 2 ( t ) λ 4 ( t ) x 1 ( t ) 4 k u z 0 λ 1 ( t ) x 2 ( t ) λ 2 ( t ) x 1 ( t ) + λ 3 ( t ) x 4 ( t ) λ 4 ( t ) x 3 ( t ) 4 k λ 1 ( t ) x 2 ( t ) + λ 2 ( t ) x 1 ( t ) λ 3 ( t ) x 4 ( t ) + λ 4 ( t ) x 3 ( t ) 4 k 0
The Hamiltonian matrix (184) gives for equations (179) and (180) a Ricatti-type system of equations for both x and λ . One must integrate these equations with the initial conditions x ( 0 ) = x 0 for the state variables, and the final conditions (the transversality condition)
λ ( T ) = ( Φ x 1 , Φ x 2 , Φ x 3 , Φ x 4 )
for the Lagrangian multipliers. Note that due to the fact that the Hamiltonian matrix H in (184) is antisymmetric, one has that
d | | x | | 2 d t = d d t ( x t · x ) = 2 x t · d x d t = x t H x = i , j x i H i j x j = 0
and the same is true for λ
d | | λ | | 2 d t = d d t ( λ t · λ ) = 2 λ t · d λ d t = λ t H λ = i , j λ i H i j λ j = 0 .
The above two equations imply that the norm of both vectors x and λ remain constant during the dynamical evolution. Also,
d d t λ t · x = d λ t d t · x + λ t · d x d t = 1 2 ( H λ ) t x + 1 2 λ t H x = 1 2 λ t H t x + λ t H x = 1 2 λ t H x + λ t H x = 0
so
d d t | | λ | | | | x | | cos ( α ) = 0
and due to the fact that the norms of x and λ are constant, the last equation implies that
d d t cos ( α ) = 0
so, the angle α between x and λ remains constant. The vector system x ( t ) , λ ( t ) evolves in time as a rigid structure with no relative motion. In this way, one can always decompose λ ( t ) in the form
λ ( t ) = a x ( t ) + b x ( t )
where a and b are constant, and x is an orthogonal vector to x . Let x = ( a 1 , a 2 , a 3 , a 4 ) , due to the fact that
x · x = x 1 ( t ) a 1 ( t ) + x 2 ( t ) a 2 ( t ) + x 3 ( t ) a 3 ( t ) + x 4 ( t ) a 4 ( t ) = 0
we can choose x as
x = a 1 ( t ) a 2 ( t ) a 3 ( t ) x 1 ( t ) a 1 ( t ) + x 2 ( t ) a 2 ( t ) + x 3 ( t ) a 3 ( t ) x 4 ( t )
By replacing (190) in (179) and (180), one obtains
x ˙ = 1 2 H [ x , a x + b x ] x
x ˙ = 1 2 H [ x , a x + b x ] x .
In this case, the controls (182) and (183) read as
u x = 1 4 k ( a x + b x ) t H x x = b 4 k x t H x x
u y = 1 4 k ( a x + b x ) t H y x = b 4 k x t H y x
or explicitly
u x = b a 1 ( t ) x 1 2 ( t ) a 1 ( t ) x 4 2 ( t ) a 2 ( t ) x 1 ( t ) x 2 ( t ) a 2 ( t ) x 3 ( t ) x 4 ( t ) a 3 ( t ) x 1 ( t ) x 3 ( t ) + a 3 ( t ) x 2 ( t ) x 4 ( t ) 4 k x 4 ( t )
u y = b a 1 ( t ) x 1 ( t ) x 3 ( t ) a 1 ( t ) x 2 ( t ) x 4 ( t ) + a 2 ( t ) x 1 ( t ) x 4 ( t ) a 2 ( t ) x 2 ( t ) x 3 ( t ) a 3 ( t ) x 3 2 ( t ) a 3 ( t ) x 4 2 ( t ) 4 k x 4 ( t )
One can consider first the case b = 0 , so λ = a x . In this case, the controls (195) and (196) read as
u x = 1 4 k a x t H x x = 0
u y = 1 4 k a x t H y x = 0
so
H = 0 0 u z 0 0 0 0 u z u z 0 0 0 0 u z 0 0
and the equation of motion for x becomes
d x 1 d t d x 2 d t d x 3 d t d x 4 d t = 1 2 0 0 u z 0 0 0 0 u z u z 0 0 0 0 u z 0 0 x 1 x 2 x 3 x 4
or
d d t x 1 ( t ) d d t x 2 ( t ) d d t x 3 ( t ) d d t x 4 ( t ) = 1 2 u z x 3 ( t ) u z x 4 ( t ) u z x 1 ( t ) u z x 2 ( t )
so
x 1 ( t ) = A 1 sin ( u z 2 t + ϕ 1 ) x 2 ( t ) = A 2 sin ( u z 2 t + ϕ 2 ) x 3 ( t ) = A 1 cos ( u z 2 t + ϕ 1 ) x 4 ( t ) = A 2 cos ( u z 2 t + ϕ 2 )
which is the usual sinusoidal solution.
The systems for x are
d d t a 1 ( t ) d d t a 2 ( t ) d d t a 3 ( t ) a 1 ( t ) x 1 ( t ) a 2 ( t ) x 2 ( t ) a 3 ( t ) x 3 ( t ) d d t x 4 ( t ) x 4 2 ( t ) + a 1 ( t ) d d t x 1 ( t ) a 2 ( t ) d d t x 2 ( t ) a 3 ( t ) d d t x 3 ( t ) x 4 ( t ) + x 1 ( t ) d d t a 1 ( t ) x 2 ( t ) d d t a 2 ( t ) x 3 ( t ) d d t a 3 ( t ) x 4 ( t ) = 1 2 u z a 3 ( t ) u z a 1 ( t ) x 1 ( t ) a 2 ( t ) x 2 ( t ) a 3 ( t ) x 3 ( t ) x 4 ( t ) u z a 1 ( t ) u z a 2 ( t )
Note that there are three variables and four equations, so the system for x can be inconsistent. For the particular case b = 0 , if one replaces the derivatives of a 1 , a 2 , a 3 from the three first equations in the fourth one, it obtains the identity
u z a 2 ( t ) 2 = u z a 2 ( t ) 2
thus, there is no contradiction this case.
Consider now the case a = 0 , so λ = b x . The systems for x are
d d t x 1 ( t ) = + b 8 k a 1 ( t ) x 1 ( t ) + a 2 ( t ) x 2 ( t ) + a 3 ( t ) x 3 ( t ) x 1 ( t ) + a 1 ( t ) x 4 ( t ) a 2 ( t ) x 3 ( t ) + a 3 ( t ) x 2 ( t ) x 4 ( t ) + b 8 k a 1 ( t ) x 1 ( t ) + a 2 ( t ) x 2 ( t ) + a 3 ( t ) x 3 ( t ) x 3 ( t ) + a 1 ( t ) x 2 ( t ) + a 2 ( t ) x 1 ( t ) a 3 ( t ) x 4 ( t ) x 4 ( t ) x 2 ( t ) x 4 ( t ) 1 2 u z x 3 ( t )
d d t x 2 ( t ) = + b 8 k a 1 ( t ) x 1 ( t ) + a 2 ( t ) x 2 ( t ) + a 3 ( t ) x 3 ( t ) x 1 ( t ) + a 1 ( t ) x 4 ( t ) a 2 ( t ) x 3 ( t ) + a 3 ( t ) x 2 ( t ) x 4 ( t ) x 3 ( t ) x 4 ( t ) b 8 k b a 1 ( t ) x 1 ( t ) + a 2 ( t ) x 2 ( t ) + a 3 ( t ) x 3 ( t ) x 3 ( t ) + a 1 ( t ) x 2 ( t ) + a 2 ( t ) x 1 ( t ) a 3 ( t ) x 4 ( t ) x 4 ( t ) x 1 ( t ) x 4 ( t ) + 1 2 u z x 4 ( t )
d d t x 3 ( t ) = b 8 k a 1 ( t ) x 1 ( t ) + a 2 ( t ) x 2 ( t ) + a 3 ( t ) x 3 ( t ) x 1 ( t ) + a 1 ( t ) x 4 ( t ) a 2 ( t ) x 3 ( t ) + a 3 ( t ) x 2 ( t ) x 4 ( t ) x 2 ( t ) x 4 ( t ) + b 8 k a 1 ( t ) x 1 ( t ) + a 2 ( t ) x 2 ( t ) + a 3 ( t ) x 3 ( t ) x 3 ( t ) + a 1 ( t ) x 2 ( t ) + a 2 ( t ) x 1 ( t ) a 3 ( t ) x 4 ( t ) x 4 ( t ) + 1 2 u z x 1 ( t )
d d t x 4 ( t ) = b 8 k a 1 ( t ) x 1 ( t ) + a 2 ( t ) x 2 ( t ) + a 3 ( t ) x 3 ( t ) x 1 ( t ) + a 1 ( t ) x 4 ( t ) a 2 ( t ) x 3 ( t ) + a 3 ( t ) x 2 ( t ) x 4 ( t ) x 1 ( t ) x 4 ( t ) b 8 k a 1 ( t ) x 1 ( t ) + a 2 ( t ) x 2 ( t ) + a 3 ( t ) x 3 ( t ) x 3 ( t ) + a 1 ( t ) x 2 ( t ) + a 2 ( t ) x 1 ( t ) a 3 ( t ) x 4 ( t ) x 4 ( t ) x 3 ( t ) x 4 ( t ) 1 2 u z x 2 ( t )
whereas the systems for x are
d d t a 1 ( t ) = + b 8 k a 1 ( t ) x 1 ( t ) + a 2 ( t ) x 2 ( t ) + a 3 ( t ) x 3 ( t ) x 1 ( t ) x 4 2 ( t ) + b 8 k a 1 ( t ) x 4 ( t ) a 2 ( t ) x 3 ( t ) + a 3 ( t ) x 2 ( t ) x 4 ( t ) a 1 ( t ) x 1 ( t ) a 2 ( t ) x 2 ( t ) a 3 ( t ) x 3 ( t ) x 4 2 ( t ) + b 8 k a 1 ( t ) x 1 ( t ) + a 2 ( t ) x 2 ( t ) + a 3 ( t ) x 3 ( t ) x 3 ( t ) + a 1 ( t ) x 2 ( t ) + a 2 ( t ) x 1 ( t ) a 3 ( t ) x 4 ( t ) x 4 ( t ) a 2 ( t ) x 4 ( t ) 1 2 u z a 3 ( t )
d d t a 2 ( t ) = + b 8 k a 1 ( t ) x 1 ( t ) + a 2 ( t ) x 2 ( t ) + a 3 ( t ) x 3 ( t ) x 1 ( t ) + a 1 ( t ) x 4 ( t ) a 2 ( t ) x 3 ( t ) + a 3 ( t ) x 2 ( t ) x 4 ( t ) a 3 ( t ) x 4 ( t ) b 8 k a 1 ( t ) x 1 ( t ) + a 2 ( t ) x 2 ( t ) + a 3 ( t ) x 3 ( t ) x 3 ( t ) + a 1 ( t ) x 2 ( t ) + a 2 ( t ) x 1 ( t ) a 3 ( t ) x 4 ( t ) x 4 ( t ) a 1 ( t ) x 4 ( t ) + 1 2 u z a 1 ( t ) x 1 ( t ) a 2 ( t ) x 2 ( t ) a 3 ( t ) x 3 ( t ) x 4 ( t )
d d t a 3 ( t ) = b 8 k a 1 ( t ) x 1 ( t ) + a 2 ( t ) x 2 ( t ) + a 3 ( t ) x 3 ( t ) x 1 ( t ) + a 1 ( t ) x 4 ( t ) a 2 ( t ) x 3 ( t ) + a 3 ( t ) x 2 ( t ) x 4 ( t ) a 2 ( t ) x 4 ( t ) + b 8 k a 1 ( t ) x 1 ( t ) + a 2 ( t ) x 2 ( t ) + a 3 ( t ) x 3 ( t ) x 3 ( t ) x 4 2 ( t ) + b 8 k a 1 ( t ) x 2 ( t ) + a 2 ( t ) x 1 ( t ) a 3 ( t ) x 4 ( t ) x 4 ( t ) a 1 ( t ) x 1 ( t ) a 2 ( t ) x 2 ( t ) a 3 ( t ) x 3 ( t ) x 4 2 ( t ) + 1 2 u z a 1 ( t )
The last equation is
a 1 ( t ) x 1 ( t ) a 2 ( t ) x 2 ( t ) a 3 ( t ) x 3 ( t ) d d t x 4 ( t ) x 4 2 ( t ) + a 1 ( t ) d d t x 1 ( t ) a 2 ( t ) d d t x 2 ( t ) a 3 ( t ) d d t x 3 ( t ) x 1 ( t ) d d t a 1 ( t ) x 2 ( t ) d d t a 2 ( t ) x 3 ( t ) d d t a 3 ( t ) x 4 ( t ) = b 8 k a 1 ( t ) x 1 ( t ) + a 2 ( t ) x 2 ( t ) + a 3 ( t ) x 3 ( t ) x 1 ( t ) + a 1 ( t ) x 4 ( t ) a 2 ( t ) x 3 ( t ) + a 3 ( t ) x 2 ( t ) x 4 ( t ) a 1 ( t ) x 4 ( t ) b 8 k a 1 ( t ) x 1 ( t ) + a 2 ( t ) x 2 ( t ) + a 3 ( t ) x 3 ( t ) x 3 ( t ) + a 1 ( t ) x 2 ( t ) + a 2 ( t ) x 1 ( t ) a 3 ( t ) x 4 ( t ) x 4 ( t ) a 3 ( t ) x 4 ( t ) 1 2 u z a 2 ( t )
In this case, if one replaces the derivatives of a 1 , a 2 , a 3 in the fourth equation, one does not obtain an identity, but instead a contradiction, unless that b = 0 . The same can be said for the general case λ ( t ) = a x ( t ) + b x ( t ) .

5.5.2. Closed-Loop Bellman Dynamics

Consider now the closed-loop feedbacks case, so one must write the Hamilton–Jacobi–Bellman equation for the optimization problem (177)
max u x , u y k u x 2 + u y 2 + 1 2 i = 1 4 j = 1 4 S x i H i j ( u ) x j = S t
where H i j ( u ) are the matrix elements of the Hamiltonian matrix in (181). The maximization, with respect to the control, gives
2 k u x + 1 2 j S x i H x i j x j = 0 2 k u y + 1 2 j S x i H y i j x j = 0
or
u x = 1 4 k ( S ) t H x x = x 1 S x 4 + x 2 S x 3 x 3 S x 2 x 4 S x 1 4 k u y = 1 4 k ( S ) t H y x = x 1 S x 2 x 2 S x 1 + x 3 S x 4 x 4 S x 3 4 k
The HJB equation becomes
1 16 k ( 8 k u z x 1 S x 3 8 k u z x 2 S x 4 8 k u z x 3 S x 1 + 8 k u z x 4 S x 2 x 1 2 S x 2 2 x 1 2 S x 4 2 + 2 x 1 x 2 S x 1 S x 2 2 x 1 x 2 S x 3 S x 4 + 2 x 1 x 4 S x 1 S x 4 + 2 x 1 x 4 S x 2 S x 3 x 2 2 S x 1 2 x 2 2 S x 3 2 + 2 x 2 x 3 S x 1 S x 4 + 2 x 2 x 3 S x 2 S x 3 x 3 2 S x 2 2 x 3 2 S x 4 2 2 x 3 x 4 S x 1 S x 2 + 2 x 3 x 4 S x 3 S x 4 x 4 2 S x 1 2 x 4 2 S x 3 2 ) = S t
Because the Lagrangian multiplier is λ = S , then by Equation (190), one could try to write
S = a x + b x
or
S x 1 = a x 1 + b a 1 ( x ) S x 2 = a x 2 + b a 2 ( x ) S x 3 = a x 3 + b a 3 ( x ) S x 4 = a x 4 b x 1 a 1 ( x ) + x 2 a 2 ( x ) + x 3 a 3 ( x ) x 4
for some functions a 1 , a 2 , a 3 . In order to integrate S, one can choose a 1 ( x ) = a 2 ( x ) = a 3 ( x ) = F ( x 4 ) , so
S x 1 = a x 1 + b F ( x 4 ) S x 2 = a x 2 + b F ( x 4 ) S x 3 = a x 3 + b F ( x 4 ) S x 4 = a x 4 b F ( x 4 ) x 4 ( x 1 + x 2 + x 3 )
By integrating S from the three first equations, one obtains
S = 1 2 a ( x 1 2 + x 2 2 + x 3 2 + x 4 2 ) + b F ( x 4 ) ( x 1 + x 2 + x 3 ) .
By replacing S in the fourth equations, one has
a x 4 + b d F ( x 4 ) d x 4 ( x 1 + x 2 + x 3 ) = a x 4 b F ( x 4 ) x 4 ( x 1 + x 2 + x 3 )
or
d F ( x 4 ) d x 4 = F ( x 4 ) x 4
so
F ( x 4 ) = 1 x 4
Thus, the S function has the structure
S = 1 2 a ( x 1 2 + x 2 2 + x 3 2 + x 4 2 ) + b ( x 1 + x 2 + x 3 ) x 4
Now, one must verify that the S function above is the solution of the HJB equation. By replacing (216) in (209), one obtains
b 16 k x 4 4 ( b x 1 4 2 b x 1 3 x 2 2 b x 1 3 x 3 b x 1 2 x 2 2 2 b x 1 2 x 2 x 3 + 2 b x 1 2 x 2 x 4 2 b x 1 2 x 3 2 3 b x 1 2 x 4 2 + 2 b x 1 x 2 2 x 4 2 b x 1 x 2 x 3 2 2 b x 1 x 3 3 4 b x 1 x 3 x 4 2 + 2 b x 1 x 4 3 b x 2 2 x 3 2 2 b x 2 2 x 3 x 4 2 b x 2 2 x 4 2 2 b x 2 x 3 3 2 b x 2 x 3 2 x 4 b x 3 4 3 b x 3 2 x 4 2 2 b x 3 x 4 3 2 b x 4 4 + 4 k u z x 1 x 2 x 4 2 + 4 k u z x 1 x 4 3 + 4 k u z x 2 2 x 4 2 + 4 k u z x 2 x 3 x 4 2 4 k u z x 3 x 4 3 + 4 k u z x 4 4 ) = 0
so the only possibility is b = 0 , so there is no x component in the solution. Thus, finally,
S = 1 2 a ( x 1 2 + x 2 2 + x 3 2 + x 4 2 )
and the controls (208) give u x = 0 , u y = 0 , from which the equations of motion for the x vector is again given by the system (203). Thus, open- and closed-loop strategies give the same answer for the dynamics, because this example belongs to the regular subclass [26] of optimal control problems mentioned at the end of the Section 5.2.

6. Conclusions

In this article, we developed an optimal control perspective on the dynamical behaviour of classical and quantum physical systems. The most crucial element of this view is the presence of feedback characterized by open- or closed-loop strategies in the system.
Thus, in quantum theory, the closed-loop strategies appear naturally due to thinking that Heisenberg’s commutation relations are a constraint in a non-commutative phase-space. Hence, this implies a relation between any quantum state’s momentum and particle position.
By taking the classical limit 0 in the full quantum Hamilton–Jacobi Equation, one arrives at the closed-loop dynamics associated with the classical Hamilton–Jacobi theory. The non-commutative character of quantum theory (generated by quantum constraint) is transferred to the classical theory through the closed-loop p x = S ( x , t ) x strategy. Since S ( x , t ) satisfies the classical Hamilton–Jacobi Equation, the dynamics generated by S ( x , t ) (under the properties of canonical transformations, whose generator is just S ( x , t ) ) are completely equivalent to those open-loop dynamics dictated by the Hamiltonian equations of motion.
From a purely classical point of view, these closed-loop strategies can be explained by the canonical transformation theory’s “two observers” character. If the solutions of the equations of motion are constant for the second observer, then their solutions look like constraints for the first one. That necessarily relates the momentum of the particle with its position for the first observer, generating, in this way, the closed-loop p x = S ( x , t ) x strategy.

Author Contributions

Conceptualization, M.C.G.; methodology, M.C.G. and R.O.H.; formal analysis, M.C.G. and M.V.; investigation, M.C.G., M.V. and R.O.H.; writing—original draft preparation M.C.G. and M.V.; writing—review and editing, M.C.G. and R.O.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mantegna, R.N.; Stanley, H.E. An Introduction to Econophysics; Cambridge University Press: Cambridge, UK, 1999. [Google Scholar]
  2. Boucheaud, J.P.; Potters, M. Theory of Financial Risk and Derivative Pricing: From Statistical Physics to Risk Management; Cambridge University Press: Cambridge, UK, 2009. [Google Scholar]
  3. Baaquie, B.E. Quantum Finance: Path Integrals and Hamiltonians for Option and Interest Rates; Cambridge University Press: Cambridge, UK, 2007. [Google Scholar]
  4. Ilinski, K. Physics of Finance: Gauge Modelling in Non–Equilibrium Pricing; Willey: Hoboken, NJ, USA, 2001. [Google Scholar]
  5. Voit, J. The Statistical Mechanics of Financial Markets; Springer: Berlin/Heidelberg, Germany, 2005. [Google Scholar]
  6. Sinha, S.; Chatterjee, A.; Chakrabortia, A.; Chakrabarti, B.K. Econophysics: An Introduction; Willey–VCH: Weinheim, Germany, 2010. [Google Scholar]
  7. Johnson, N.F.; Jefferies, P.; Hui, P.M. Financial Market Complexity: What Physics Can Tell Us about Market Behaviour; Oxford University Press: Oxford, UK, 2003. [Google Scholar]
  8. Dash, J. Quantitative Finance and Risk Management: A Physicist’s Approach; World Scientific: Singapore, 2016. [Google Scholar]
  9. Haven, E. A Black–Scholes Schrödinger option price: ‘bit’ versus ‘qubit’. Phys. A 2003, 324, 201–206. [Google Scholar] [CrossRef]
  10. Haven, E. A discussion on embedding the Black–Scholes option price model in a quantum physics setting. Phys. A 2002, 304, 507–524. [Google Scholar] [CrossRef]
  11. Baaquie, B.E. Quantum field theory of treasury bonds. Phys. Rev. E 2001, 64, 016121. [Google Scholar] [CrossRef] [PubMed]
  12. Baaquie, B.E. Quantum finance Hamiltonian for coupon bond European and barrier options. Phys. Rev. E 2008, 77, 036106. [Google Scholar] [CrossRef] [PubMed]
  13. Baaquie, B.E. Interest rates in quantum finance: The Wilson expansion and Hamiltonian. Phys. Rev. E 2009, 80, 046119. [Google Scholar] [CrossRef]
  14. Contreras, M. Stochastic volatility models at ρ =±1 as a second-class constrained Hamiltonian systems. Phys. A 2015, 405, 289–302. [Google Scholar] [CrossRef]
  15. Contreras, M.; Bustamante, M. Multi-asset Black–Scholes model as a variable second-class constrained dynamical system. Phys. A 2016, 457, 540–572. [Google Scholar]
  16. Contreras, M.; Pellicer, R.; Villena, M.; Ruiz, A. A quantum model for option pricing: When Black–Scholes meets Schrödinger and its semi–classic limit. Phys. A 2010, 329, 5447–5459. [Google Scholar] [CrossRef]
  17. Kamien, M.I.; Schwartz, N.L. The Calculus of Variations and Optimal Control in Economics and Management; Elsevier Science: Amsterdam, The Netherlands, 1991. [Google Scholar]
  18. Sethi, S.P.; Thompson, G.L. Optimal Control Theory: Applications to Management Science and Economics, 2nd ed.; Springer Science + Business Media: New York, NY, USA, 2009. [Google Scholar]
  19. Caputo, M.R. Foundations of Dynamic Economic Analysis: Optimal Control Theory and Applications; Cambridge University Press: Cambridge, UK, 2005. [Google Scholar]
  20. Weitzman, M.L. Income, Wealth, and the Maximum Principle; Harvard University Press: Cambridge, MA, USA, 2007. [Google Scholar]
  21. Dockner, E.J.; Jorgensen, S.; Long, N.V. Differential Games in Economics and Management Science; Cambridge University Press: Cambridge, UK, 2001. [Google Scholar]
  22. Contreras, M.; Pellicer, R.; Villena, M. Dynamic optimization and its relation to classical and quantum constrained systems. Phys. A 2017, 479, 12–25. [Google Scholar] [CrossRef]
  23. Hojman, S.A. Optimal Control and Dirac’s Theory of Singular Hamiltonian Systems. Unpublished work.
  24. Itami, T. Quantum Mechanical Theory of Nonlinear Control (in IFAC Nonlinear Control Systems); IFAC Publications: St. Petersburg, Russia, 2001; p. 1411. [Google Scholar]
  25. Contreras, M.; Peña, J.P. The quantum dark side of the optimal control theory. Phys. A 2019, 515, 450–473. [Google Scholar]
  26. Contreras, M.; Peña, J.P.; Aros, R. Second class constraints and the consistency of optimal control theory in phase space. Phys. A 2021, 562, 125–335. [Google Scholar] [CrossRef]
  27. Pontryagin, L.S.; Boltyanskii, V.G.; Gamkrelidze, R.V.; Mishchenko, E.F. The Mathematical Theory of Optimal Processes; CRC Press: Boca Raton, FL, USA, 1987. [Google Scholar]
  28. Bellman, R. The theory of dynamic programming. Bull. Am. Math. Soc. 1954, 60, 503–516. [Google Scholar] [CrossRef]
  29. Yang, X.; Wu, L.; Zhang, H. A space-time spectral order sinc-collocation method for the fourth-order nonlocal heat model arising in viscoelasticity. Appl. Math. Comput. 2023, 457, 128192. [Google Scholar] [CrossRef]
  30. Yang, X.; Zhang, Q.; Yuan, G. On positivity preservation in nonlinear finite volume method for multi-term fractional subdiffusion equation on polygonal meshes. Nonlinear Dyn 2018, 92, 595–612. [Google Scholar] [CrossRef]
  31. Ali, M.A.; Budak, H.; Zhang, Z.; Yildirim, H. Some new Simpson’s type inequalities for coordinated convex functions in quantum calculus. Math. Methods Appl. Sci. 2020, 44, 4515–4540. [Google Scholar] [CrossRef]
  32. Budak, H.; Erden, S.; Ali, M.A. Simpson and Newton type inequalities for convex functions via newly defined quantum integrals. Math. Methods Appl. Sci. 2020, 44, 378–390. [Google Scholar] [CrossRef]
  33. Erickson, G.M. Differential game models of advertising competitions. J. Political Econ. 1973, 8, 637–654. [Google Scholar] [CrossRef]
  34. Dirac, P.A.M. Generalized Hamiltonian dynamics. Proc. R. Soc. Lond. A 1958, 246, 326–332. [Google Scholar] [CrossRef]
  35. Dirac, P.A.M. Lectures on Quantum Mechanics; Yeshiva University Press: New York, NY, USA, 1967. [Google Scholar]
  36. Teitelboim, C.; Henneaux, M. Quantization of Gauge Systems; Princeton University Press: Princeton, NJ, USA, 1994. [Google Scholar]
  37. Rothe, H.J.; Rothe, K.D. Classical and Quantum Dynamics of Constrained Hamiltonian Systems (World Scientific Lectures Notes in Physics v81); World Scientific: Singapore, 2010. [Google Scholar]
  38. Fetter, A.L.; Walecka, J.D. Theoretical Mechanics of Particles and Continua; Dover Publications: Mineola, NY, USA, 2003. [Google Scholar]
  39. Goldstein, H. Classical Mechanics, 3rd ed.; Pearson: London, UK, 2001. [Google Scholar]
  40. Contreras, G.M. Dirac’s Method in a Non-Commutative Phase Space; UMCE: Santiago, Chile, 2023; in preparation. [Google Scholar]
  41. Rothe, K.D.; Scholtz, F.G. On the Hamilton–Jacobi equation for second-class constrained systems. Ann. Phys. 2003, 308, 639–651. [Google Scholar] [CrossRef]
  42. Tawfik, A.N.; Diab, A.M. A review of the generalized uncertainty principle. Rep. Prog. Phys. 2015, 78, 126001. [Google Scholar] [CrossRef] [PubMed]
  43. Bruneton, J.P.; Larena, J. Quantum theory of the generalised uncertainty principle. Gen. Relativ. Gravit. 2017, 49, 56. [Google Scholar] [CrossRef]
  44. Pedram, P. A class of GUP solutions in deformed quantum mechanics. Int. J. Mod. Phys. D 2010, 19, 2003–2009. [Google Scholar] [CrossRef]
  45. Seifi, M.; Sefiedgar, A.S. The effects of the covariant generalized uncertainty principle on quantum mechanics. Can. J. Phys. 2022, 101, 242–247. [Google Scholar] [CrossRef]
  46. Luciano, G.G.; Petruzziello, L. Generalized uncertainty principle and its implications on geometric phases in quantum mechanics. Eur. Phys. J. Plus 2021, 136, 179. [Google Scholar] [CrossRef]
  47. Scardigli, F. The deformation parameter of the generalized uncertainty principle. J. Phys. Conf. Ser. 2019, 1275, 012004. [Google Scholar] [CrossRef]
  48. Casadio, R.; Scardigli, F. Generalized Uncertainty Principle, Classical Mechanics, and General Relativity. Phys. Lett. B 2020, 807, 135558. [Google Scholar] [CrossRef]
  49. Reginatto, M.; Hall, M.J.W. Entangling quantum fields via a classical gravitational interaction. J. Phys. Conf. Ser. 2019, 1275, 012039. [Google Scholar] [CrossRef]
  50. Övgün, A. Entangled Particles Tunneling From a Schwarzschild Black Hole immersed in an Electromagnetic Universe with GUP. Int. J. Theor. Phys. 2016, 55, 2919–2927. [Google Scholar] [CrossRef]
  51. Park, D. Quantum entanglement with generalized uncertainty principle. Nucl. Phys. B 2022, 977, 115736. [Google Scholar] [CrossRef]
  52. Guo, X.; Wang, P.; Yang, H. The classical limit of minimal length uncertainty relation: Revisit with the Hamilton-Jacobi method. J. Cosmol. Astropart. Phys. 2016, 2016, 62. [Google Scholar] [CrossRef]
  53. Reginatto, M. Exact Uncertainty Principle and Quantization: Implications for the Gravitational Field. Braz. J. Phys. 2005, 35, 476–480. [Google Scholar] [CrossRef]
  54. Dehaghani, N.B.; Pereira, F.L. Optimal Control of Quantum Systems by Pontryagin Maximum Principle. U.Porto J. Eng. 2022, 8, 194–201. [Google Scholar] [CrossRef]
  55. D’Alessandro, D.; Dahleh, M. Optimal control of two-level quantum systems. In Proceedings of the 2000 American Control Conference, ACC (IEEE Cat. No.00CH36334), Chicago, IL, USA, 28–30 June 2000; Volume 6, pp. 3893–3897. [Google Scholar] [CrossRef]
Figure 1. Constraint over the phase-space.
Figure 1. Constraint over the phase-space.
Symmetry 15 02033 g001
Figure 2. Q and P constraints over the phase-space.
Figure 2. Q and P constraints over the phase-space.
Symmetry 15 02033 g002
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Contreras González, M.; Villena, M.; Ortiz Herrera, R. An Optimal Control Perspective on Classical and Quantum Physical Systems. Symmetry 2023, 15, 2033. https://doi.org/10.3390/sym15112033

AMA Style

Contreras González M, Villena M, Ortiz Herrera R. An Optimal Control Perspective on Classical and Quantum Physical Systems. Symmetry. 2023; 15(11):2033. https://doi.org/10.3390/sym15112033

Chicago/Turabian Style

Contreras González, Mauricio, Marcelo Villena, and Roberto Ortiz Herrera. 2023. "An Optimal Control Perspective on Classical and Quantum Physical Systems" Symmetry 15, no. 11: 2033. https://doi.org/10.3390/sym15112033

APA Style

Contreras González, M., Villena, M., & Ortiz Herrera, R. (2023). An Optimal Control Perspective on Classical and Quantum Physical Systems. Symmetry, 15(11), 2033. https://doi.org/10.3390/sym15112033

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop