An Implementation Method for the Geodesics with Constraints on Heisenberg Manifolds ()
1. Introduction
Geodesics plays an important role in many applications, especially in nuclear physics, image processing, ··· Ovidiu Calin and Vittorio Mangione [1] considered the Heisenberg manifold structure to provide a qualitative characterization for geodesics under nonholonomic constraints. This method offers an excellent description or the solution of Euler—Lagrange equation associated to lagrangians with linear and quadratic constraints. Therefore it is highly desirable to consider the implementation of this interesting mathematical method. In this paper we investigate the implementation of the method presented in reference [1]. Due to the fact that the mathematical work in reference [1] lacks some details that are necessary for implementation, we address this issue by including such required details with complete proofs, after presenting the method described elsewhere [1] in an appropriate manner we implement and simulate mathematical. This implementation approach can successfully well illustrate the variation of some parameters such as Christoffel coefficients, ··· and ··· that are required in the determination of tensors. Our approach is also of great importance especially in the case where the determination of the geodesics is carried out by minimizing a performance index. Therefore our implementation approach can be considered as an attractive complement for the work of [1].
The working hypotheses:
We take the examples studied in [1] with the following:
1) Expressions of the vector fields
(1)
2) The Heisenberg Laplacian operator
(2)
3) The sub-Riemannian geometry can be defined on
by:
(3)
4) If
is the trajectory of a particle of mass m = 1 its energy is given by
(4)
where
, and
is a Riemannian metric which will be specified later.
5) We consider successively the two expressions:
(5)
![](https://www.scirp.org/html/33-7400972\78d8a491-f3b7-45d7-92e9-a714d6c32848.jpg)
where w is the 1-form such
![](https://www.scirp.org/html/33-7400972\ac18eef6-9236-434c-836f-b5419c223fdd.jpg)
And
is a constant, which has the physical significance of a charge.
is a 1-form which will be defined later.
6) We consider as in [2], the Heisenberg group
![](https://www.scirp.org/html/33-7400972\1e636619-660f-437d-ae9b-3fd9136f3e9f.jpg)
This group is non commutative and the law of the group is polynomial and can be written in ![](https://www.scirp.org/html/33-7400972\58cad18c-84ff-4d7f-be9e-d4e70c8a50f4.jpg)
![](https://www.scirp.org/html/33-7400972\51aaf3bc-cf11-43d4-b8a8-50213681193e.jpg)
The Lie algebra of H is spanned by the matrices
![](https://www.scirp.org/html/33-7400972\90076930-dde1-489d-9288-e5f0bfa6d6ad.jpg)
for which the following relations hold
![](https://www.scirp.org/html/33-7400972\c57cdd42-3476-4700-87a8-636438a11e39.jpg)
Proposition:
If A and B are in H,
![](https://www.scirp.org/html/33-7400972\ca27432d-227b-4936-b137-5fe01d6448c4.jpg)
This relation is not so senseless even if it can be very easily proved with a little computation.
It is indeed coming from the Baker-Campb-Hausdorff formula which expresses the product of the exponential of two matrices as the exponential of some quantity.
To be more precise, for two matrices M and N:
![](https://www.scirp.org/html/33-7400972\8ea970da-e139-410c-b632-04c828843ad4.jpg)
where
is a Lie series which depends on the iterated brackets of M and N:
![](https://www.scirp.org/html/33-7400972\03efeec4-2c7c-43d1-a814-a9c6eea76f3b.jpg)
In this case of the Heisenberg group whose Lie algebra is nilpotent of order 2, this serie stops after the first bracket term. We prefer to work with the exponential coordinates are the coordinates in the Lie algebra.
We identify
then with the triple
such that:
![](https://www.scirp.org/html/33-7400972\5a99ee21-413e-4604-97e4-dceefd2152f2.jpg)
The group law in these coordinates becomes:
![](https://www.scirp.org/html/33-7400972\a448bf28-004c-4928-a931-f28918632fcc.jpg)
And the inverse element is:
![](https://www.scirp.org/html/33-7400972\55670637-3fa9-4af6-aa6b-13e07ab99bc6.jpg)
The expressions of the left invariant vector fields in these exponential coordinates are then
(6)
Whereas the right-invariant vector fields is written:
![](https://www.scirp.org/html/33-7400972\7caff551-335c-4f2e-a372-6a6dc8c04550.jpg)
Reference [3] was the first to check easily that:
(7)
And
![](https://www.scirp.org/html/33-7400972\3ba81655-8bbe-479e-bf24-400af1730b3f.jpg)
The left invariant frame
,
is a basis for the horizontal fibration
![](https://www.scirp.org/html/33-7400972\a58bd4c4-47b6-424f-b3b4-93455a70b115.jpg)
Note that
(8)
Is a contact form in
i.e.
(
never vanishes); since
, it follows that K is not involutive.The distribution K will be called the horizontal distribution.
A more detailed about sub-Riemannian structures see in [4].
A curve
is called horizontal if
![](https://www.scirp.org/html/33-7400972\683ce41f-ab99-4fa2-a4c1-b6d69f631e2c.jpg)
Or
(9)
The vector fields (6) defines a unique Riemannian metric h such as
![](https://www.scirp.org/html/33-7400972\8f543edc-75f4-482f-8441-c5986d9c81ee.jpg)
And
.
For the coefficients calculus, we have used a little Mathematica program:
![](https://www.scirp.org/html/33-7400972\f967a60f-4064-42e3-ae8e-28a7509ee07f.jpg)
![](https://www.scirp.org/html/33-7400972\ec670fb0-955e-4b0e-99ac-f6e9e609846c.jpg)
![](https://www.scirp.org/html/33-7400972\196876fc-809f-4073-a46f-62a89e87b975.jpg)
![](https://www.scirp.org/html/33-7400972\4004174b-8dab-4816-8b81-9c3e86a6e7b5.jpg)
![](https://www.scirp.org/html/33-7400972\5574a739-697a-46c2-b46f-29009a56d14a.jpg)
![](https://www.scirp.org/html/33-7400972\3a4f9fef-11cd-4cd3-a65b-41051d74fd3e.jpg)
![](https://www.scirp.org/html/33-7400972\603dea0a-4bfa-4d6e-98fa-903df4107621.jpg)
![](https://www.scirp.org/html/33-7400972\372d044d-ae3a-4cea-9987-9a251ed90f06.jpg)
![](https://www.scirp.org/html/33-7400972\c4f982d8-fcd4-4fb2-a775-4b75b817d6ff.jpg)
![](https://www.scirp.org/html/33-7400972\4b83e7a7-c56f-4a91-8beb-db1a3048db22.jpg)
![](https://www.scirp.org/html/33-7400972\903782d6-aa0e-4f35-9abe-be557734f85c.jpg)
![](https://www.scirp.org/html/33-7400972\4591a880-7876-41be-a4cf-a4c9e29ac5b6.jpg)
![](https://www.scirp.org/html/33-7400972\4c193659-4223-42de-9067-de6085d91486.jpg)
![](https://www.scirp.org/html/33-7400972\d02fca9c-d89b-4ceb-9b40-201b7cad88f4.jpg)
![](https://www.scirp.org/html/33-7400972\da210648-feaa-4ecf-a26b-2b0d71b46b96.jpg)
![](https://www.scirp.org/html/33-7400972\89877ba1-7849-4056-8a86-11c00d73c80a.jpg)
So, we have :
,
,
,
,![](https://www.scirp.org/html/33-7400972\33971a8c-19a7-4800-a166-15c21ad85eeb.jpg)
.
As the coefficients
are symmetric in i and j the matrix of coefficients is
(10)
For a detailed study of sub-Riemannian geodesics on Heisenbeg group see in [5].
2. Main Results [6]
Heisenberg group case:
We shall construct the Euler-Lagrange equation for the Lagrangian (5) in the Levi-Civita connection form.
Lemma 1
If
are the components of the Ricci tensor with respect to the metric
on
then
(11)
Proof
We will calculate the coefficients of the Ricci tensor from the following relation:
![](https://www.scirp.org/html/33-7400972\e20c9902-83b7-4e0a-a7e1-7760a8e15fe6.jpg)
After, we will compare these results with those of theorem.
For the first calculation, we use the mathematica program:
Clear [coord, metric, inversemetric, affine, riemann, ricci,scalar, x, y, t]
![](https://www.scirp.org/html/33-7400972\bd85d65c-91b6-4f83-9097-ba027535a534.jpg)
![](https://www.scirp.org/html/33-7400972\b05ff1a0-c522-4763-ac98-96b555f9ec28.jpg)
![](https://www.scirp.org/html/33-7400972\d283159e-5eb7-470b-b9a6-f606685b619a.jpg)
![](https://www.scirp.org/html/33-7400972\96e683e5-1b97-49f5-a60e-f6f9d440b8e1.jpg)
Metric // MatrixForm
![](https://www.scirp.org/html/33-7400972\4586e3bf-eda8-4d44-8ce4-618713caf6ac.jpg)
Imversemetric = Simplify [Inverse [metric]]
![](https://www.scirp.org/html/33-7400972\2acc50aa-32e6-4081-92e9-2190b876a6da.jpg)
inversemetric // MatrixForm
![](https://www.scirp.org/html/33-7400972\cb142d4a-c03a-4b01-9489-710ad344a7a2.jpg)
![](https://www.scirp.org/html/33-7400972\8be4169e-9792-4858-9ac3-ead8503f8690.jpg)
![](https://www.scirp.org/html/33-7400972\05e5f594-b7a8-4d0a-b3be-153b751f1401.jpg)
![](https://www.scirp.org/html/33-7400972\5ac182fe-150a-45c2-9cb7-1f27b03e6c3c.jpg)
![](https://www.scirp.org/html/33-7400972\7a3b9562-c55f-4bfd-b15f-b4434970791b.jpg)
![](https://www.scirp.org/html/33-7400972\ea13fe0a-5b0e-4676-af8f-d36cc526ff6b.jpg)
![](https://www.scirp.org/html/33-7400972\fd69348c-930f-4d75-8efa-9d6558bc76d8.jpg)
![](https://www.scirp.org/html/33-7400972\f4b75737-5f87-4fde-aa27-1ea21df393ca.jpg)
![](https://www.scirp.org/html/33-7400972\f5fcebeb-70d1-4dc8-b7f2-ba7776dcaad2.jpg)
TableForm
,
![](https://www.scirp.org/html/33-7400972\b6502e54-07c4-4fb9-a01f-d54950e057a9.jpg)
![](https://www.scirp.org/html/33-7400972\5d3891ae-f7ec-43c2-adb5-d252815174f2.jpg)
![](https://www.scirp.org/html/33-7400972\1c19393f-c4f8-44f6-8d16-d00600331421.jpg)
![](https://www.scirp.org/html/33-7400972\90111b53-2256-4977-986d-ac369eec424a.jpg)
![](https://www.scirp.org/html/33-7400972\b275308e-d8a4-4df4-8012-c87e25597c4a.jpg)
![](https://www.scirp.org/html/33-7400972\a70612f6-ca61-4d2f-bad9-cd748093634c.jpg)
![](https://www.scirp.org/html/33-7400972\539e249e-27e2-48d7-90bc-ff9d137d164c.jpg)
LikewiseUsing (10) where we swap
with
, we get the Equation (11) on components:
![](https://www.scirp.org/html/33-7400972\a281dddd-6c22-4951-b049-e6ea29e9c266.jpg)
![](https://www.scirp.org/html/33-7400972\26c28c48-9a8c-4918-b45f-0c03a699b268.jpg)
![](https://www.scirp.org/html/33-7400972\ad0176f6-1bb3-47e2-93a8-dc66e70d41f1.jpg)
![](https://www.scirp.org/html/33-7400972\5e22f6da-fe08-47f7-8542-f2946a34956b.jpg)
Remarque:
To show the values of the Riemann tensor, add the following command :
![](https://www.scirp.org/html/33-7400972\97804a09-f02a-444f-960e-62717b91a2e2.jpg)
Corollary 1
If
is an horizontal curve ,then
(12)
where
does not depend on ![](https://www.scirp.org/html/33-7400972\a6a05076-d933-4fdb-b851-d5e1f9855c18.jpg)
Proof If
is horizontal,
. In this case,
![](https://www.scirp.org/html/33-7400972\00ec36db-54e2-40ce-bd69-7e93539c05c1.jpg)
so we have
![](https://www.scirp.org/html/33-7400972\ab3fb464-74e7-4f55-90bc-fe0978789ccb.jpg)
We see that
does not depend on
.
Using lemma1:
![](https://www.scirp.org/html/33-7400972\7ecba27c-39a1-4136-898a-e325598ac9cf.jpg)
The next proposition is a generalization of the previous corollary to any vector field.
Proposition 1
For any vector field V we have
(13)
Proof Using the metric
, we can write
(14)
where
![](https://www.scirp.org/html/33-7400972\3d046a34-1b81-4e77-8161-2e0fc33fc16b.jpg)
we have
![](https://www.scirp.org/html/33-7400972\d5fa3024-7cb0-41c7-abd6-3dcbdca9b3e4.jpg)
And (14) becomes
![](https://www.scirp.org/html/33-7400972\769b62e4-b351-4422-982c-048a142c0d90.jpg)
Using lemma1:
![](https://www.scirp.org/html/33-7400972\6dba01a3-6a0c-484a-a7b7-786926bbb7ff.jpg)
Corollary 2
The actions
and
both with respect to the metric
reach the extrema for the same functions
.
In particular, the extrema will be geodesics in the metric with coefficients
.
It is interesting that, even if the Lagrangian (5) has a non-holonomic constraint, the minimizers still behave as geodesics in a certain metric. This is given in the following result.
Theorem 1
The Euler-Lagrange equation for the Lagrangian (5) is
(15)
where
, with the coefficients given by
(16)
3. A More General Case [6]
Heisenberg manifold case
We have investigated the case when the vector fields are given by the formula (3). We shall deal in this section with the more general case of vector fields.
And ![](https://www.scirp.org/html/33-7400972\f88655e0-a3c1-45ce-bdb9-3bf590e4cc36.jpg)
With A(y), B(x) are smooth functions. The 1-form w in this case is
(17)
One may check that
(18)
Another important 1-form is
(19)
A computation shows the vector fields
,
and
are orthonormal in the Riemannian metric.
For computing the coefficients, we have use a little Mathematica program.
![](https://www.scirp.org/html/33-7400972\6c7498ff-ba86-4077-8ad5-2e0023879b90.jpg)
![](https://www.scirp.org/html/33-7400972\63204561-5ebf-4468-921c-777c051eacd5.jpg)
![](https://www.scirp.org/html/33-7400972\41fe781b-1604-4163-8c65-2e6f8d4890ca.jpg)
![](https://www.scirp.org/html/33-7400972\40ac341e-4de8-47dc-b755-db6907f2ea42.jpg)
![](https://www.scirp.org/html/33-7400972\f4b2cf77-be1f-4e4f-bfd3-24a0d3f205e7.jpg)
![](https://www.scirp.org/html/33-7400972\250140a9-2445-460e-9b9f-92cec71caa55.jpg)
![](https://www.scirp.org/html/33-7400972\7e4e89bc-5141-4d96-aa7c-2aa69f040582.jpg)
![](https://www.scirp.org/html/33-7400972\4a7ff357-93df-4388-9acf-9281562c55b3.jpg)
NSolve
![](https://www.scirp.org/html/33-7400972\d7dfa1c7-e28f-4bf8-bbe1-21e11e738dba.jpg)
![](https://www.scirp.org/html/33-7400972\92d054e7-871c-4e7f-baa3-6beaf7649ada.jpg)
![](https://www.scirp.org/html/33-7400972\b4d95e0f-88e4-4949-8780-971f4cde91c0.jpg)
![](https://www.scirp.org/html/33-7400972\8f5663a4-e785-47e1-8981-c216f97064a4.jpg)
![](https://www.scirp.org/html/33-7400972\f10a8bcf-4ee4-44d6-8dbe-4b209423f35d.jpg)
![](https://www.scirp.org/html/33-7400972\47bf16ae-cafd-4d05-8345-28a3ac6ac349.jpg)
![](https://www.scirp.org/html/33-7400972\c7426cb2-6f29-473b-b94c-6e0b66c1f79e.jpg)
![](https://www.scirp.org/html/33-7400972\4070c805-fafc-471d-878a-9dd863915d32.jpg)
![](https://www.scirp.org/html/33-7400972\89d52f90-eaa9-4285-af7d-1caaa737719e.jpg)
![](https://www.scirp.org/html/33-7400972\80aa2039-efb3-4c79-b941-1fa12267bc38.jpg)
,
![](https://www.scirp.org/html/33-7400972\e2c58440-eea9-4862-b4fc-e08da91ace02.jpg)
So, we have:
, ![](https://www.scirp.org/html/33-7400972\500cfafe-f57f-434d-b67b-7d28c6711b80.jpg)
,
![](https://www.scirp.org/html/33-7400972\22967725-920c-4350-857d-66f1f32e4d7a.jpg)
,
.
As the Riemannian metric is symmetric, we obtain
(20)
We shall consider the following Lagrangian with a quadratic potential constraint
(21)
When
and
we get the Lagrangian in (5). In this case
.
The following result is a generalization of lemma1.We shall denote by
,
(22)
Lemma 2
If
are the components of the Ricci tensor with respect to the metric given in (20), then
(23)
where
is obtained by flipping the sign in(20)
![](https://www.scirp.org/html/33-7400972\2b4f0fc3-860d-432e-bb30-be5c9d15d03c.jpg)
And R is the Ricci scalar
Proof
We will just replace the coefficients of
and
in the relation (24), and next we will use the following relation to compare the two results :
![](https://www.scirp.org/html/33-7400972\c6cd9bbd-a0fa-4666-a1dc-7ba949d8b321.jpg)
For this, we use the mathematica program.
1) We have
, then
![](https://www.scirp.org/html/33-7400972\71970641-eba1-4ca4-8989-ef6d7cdccd24.jpg)
![](https://www.scirp.org/html/33-7400972\701bcc2f-0823-4d93-ae44-d3fdbb2e881e.jpg)
![](https://www.scirp.org/html/33-7400972\612d73bd-0ebe-4260-bd71-31171ec1be67.jpg)
![](https://www.scirp.org/html/33-7400972\acbce8db-3fac-4489-b7a0-d0f05a9ba4e6.jpg)
![](https://www.scirp.org/html/33-7400972\f85e15e7-c808-4a74-8e8b-2efac9a4a7a9.jpg)
And the Ricci scalar is :
(24)
2) we put :
, ![](https://www.scirp.org/html/33-7400972\17272488-c928-4e3e-8972-d1732cce82cd.jpg)
Clear [coord, metric, inversemetric, affine, riemann, ricci, scalar, x, y, t, A, B]
![](https://www.scirp.org/html/33-7400972\91dbd96e-45fa-4831-ab35-6982b01b18a9.jpg)
![](https://www.scirp.org/html/33-7400972\01c19eb5-fbb1-45e4-9290-2403c657eaba.jpg)
![](https://www.scirp.org/html/33-7400972\b07c2145-19a3-4aa7-bca7-cea177affb32.jpg)
![](https://www.scirp.org/html/33-7400972\ecde77e5-86a6-44fe-8265-60721f034ddc.jpg)
![](https://www.scirp.org/html/33-7400972\d0a150a9-21ff-4f5a-b0a4-4c507cfe497c.jpg)
![](https://www.scirp.org/html/33-7400972\fb7d9145-f5d1-4105-b9ff-e8a6de60576b.jpg)
![](https://www.scirp.org/html/33-7400972\48d7bb33-2b52-45e9-9e76-df449f7aaf12.jpg)
![](https://www.scirp.org/html/33-7400972\79bb940f-1fcf-47d4-8fb2-f90bef10c342.jpg)
![](https://www.scirp.org/html/33-7400972\f4e316d4-4722-443c-a608-6dc718b8579a.jpg)
![](https://www.scirp.org/html/33-7400972\54cfd8c9-8fe0-43b5-a44a-643bedffa7fb.jpg)
![](https://www.scirp.org/html/33-7400972\d083c67e-8543-492b-884c-7aaaeabe6cd7.jpg)
![](https://www.scirp.org/html/33-7400972\acf18904-dba1-4ee4-980d-49d26f6b7731.jpg)
![](https://www.scirp.org/html/33-7400972\e3dcbaf0-5ca5-426f-ad95-b8f5a0e30d9d.jpg)
![](https://www.scirp.org/html/33-7400972\aa283a5f-9e6b-4ad4-91ba-a419f05bc22c.jpg)
![](https://www.scirp.org/html/33-7400972\8e46d7bf-9e8b-4fee-89a6-f3dfa6eb9895.jpg)
![](https://www.scirp.org/html/33-7400972\d23a176a-664e-4481-bb33-257529e2404f.jpg)
![](https://www.scirp.org/html/33-7400972\a4868840-980c-4be5-9592-bcb2f30074f9.jpg)
![](https://www.scirp.org/html/33-7400972\58cdcea3-5fe8-43c7-9c6b-090082e429c5.jpg)
![](https://www.scirp.org/html/33-7400972\6568a9dd-f152-4ace-bb03-0b4890ecae34.jpg)
![](https://www.scirp.org/html/33-7400972\34e54340-65b8-48fa-ab1f-b0d4b49d8278.jpg)
![](https://www.scirp.org/html/33-7400972\e691377a-f867-4cd7-a1e5-999cf12a00ab.jpg)
![](https://www.scirp.org/html/33-7400972\1d8a42af-694b-4412-9cc9-0f19f686bed0.jpg)
![](https://www.scirp.org/html/33-7400972\0ca0a161-198c-4693-9056-40a9fbb8cc4e.jpg)
![](https://www.scirp.org/html/33-7400972\f575f762-4726-41a8-b355-115267ba4dcc.jpg)
![](https://www.scirp.org/html/33-7400972\ab77f45f-0e71-4624-8e69-57d3b7d90ad3.jpg)
![](https://www.scirp.org/html/33-7400972\b14712eb-df17-4acf-93a4-2bb3469663c6.jpg)
![](https://www.scirp.org/html/33-7400972\8a3f2a30-e36d-4ca4-94c8-aac67d6080a1.jpg)
![](https://www.scirp.org/html/33-7400972\3efc3330-7317-441e-b9a8-767b738fbabe.jpg)
![](https://www.scirp.org/html/33-7400972\c9a4377b-a4aa-4900-88a0-865b8f2003da.jpg)
![](https://www.scirp.org/html/33-7400972\cc591664-ec6e-40f6-83aa-dbe78a2514ce.jpg)
![](https://www.scirp.org/html/33-7400972\538612c4-8d2d-4cba-b602-fca985470315.jpg)
![](https://www.scirp.org/html/33-7400972\c48f483a-08a1-463a-8a10-284c96658e6a.jpg)
,
![](https://www.scirp.org/html/33-7400972\65b4e1e1-0035-4454-bd09-e3a5b9fa1421.jpg)
![](https://www.scirp.org/html/33-7400972\0c8d1cca-2434-4864-b636-ab48b1640ca4.jpg)
![](https://www.scirp.org/html/33-7400972\f7c670c0-8142-4de2-9d42-304d8b12e6aa.jpg)
![](https://www.scirp.org/html/33-7400972\96ce5ca1-acaa-4e02-b68d-eab2f7bc2fed.jpg)
![](https://www.scirp.org/html/33-7400972\d339c78d-0de5-4766-9930-407f7b582fb3.jpg)
![](https://www.scirp.org/html/33-7400972\57316303-bb39-45fb-9be4-ef2be2a65be8.jpg)
![](https://www.scirp.org/html/33-7400972\ebc3e18d-af90-42b9-9f02-55f70e68e6c5.jpg)
![](https://www.scirp.org/html/33-7400972\02875e20-5c36-42ef-88f1-41dc85247e99.jpg)
![](https://www.scirp.org/html/33-7400972\067105eb-68b9-4735-abee-c0c205a6d201.jpg)
![](https://www.scirp.org/html/33-7400972\ca38898a-eb4d-4903-8196-a56488e642b0.jpg)
![](https://www.scirp.org/html/33-7400972\7ac615c7-12bf-4581-97fe-2d5a4d1e753a.jpg)
![](https://www.scirp.org/html/33-7400972\feab6026-2eea-49e8-80b7-d8a6c4a67a60.jpg)
![](https://www.scirp.org/html/33-7400972\96597d67-8db6-448f-9110-f61d9d01aac3.jpg)
![](https://www.scirp.org/html/33-7400972\a2121b95-f503-44eb-ad7c-444cadf05aa3.jpg)
Remarque:
To show the values of the Riemann tensor ,add the following command :
![](https://www.scirp.org/html/33-7400972\f4274346-8b72-43cd-b5c3-0ff37ef2d29e.jpg)
![](https://www.scirp.org/html/33-7400972\b9d446fd-13bb-4d4c-900e-33ace7ee6299.jpg)
Lemma 3
Consider the curve
(25)
Then
![](https://www.scirp.org/html/33-7400972\3f869042-d73c-4cb9-a55b-a724204205bf.jpg)
Proof
Let
be the tangent vector field Then
![](https://www.scirp.org/html/33-7400972\292542e8-b0c6-4d5f-b733-e7c72fea61ea.jpg)
We note that
and w measure the departure from the Heisenberg structure and horizontality respectively. When the structure is Heisenberg,
and when V is horizontal,
.
In the following we give a global interpretation for the difference
in terms of the horizontal distribution
.
Lemma 4
(26)
With
![](https://www.scirp.org/html/33-7400972\38644dc9-1368-4747-90b0-335a625d7147.jpg)
Furthermore
![](https://www.scirp.org/html/33-7400972\1809579e-da07-4fd4-be1c-a133592b9133.jpg)
Proof
For the relation (2.28):
![](https://www.scirp.org/html/33-7400972\c3876843-fa55-4c4d-b376-e56b976ece81.jpg)
For the second part we have:
![](https://www.scirp.org/html/33-7400972\1d86ef76-495f-42f9-ad24-6ee763c9f4b7.jpg)
Denote
, ![](https://www.scirp.org/html/33-7400972\84adb410-6776-4fb5-8494-694b5f094785.jpg)
Denote also
.
Theorem 2
and
reach the extrema for the same functions. In particular, the extrema will be geodesics in the metric
and obey the equation
where
is the Levi-Civita type connection defined by ![](https://www.scirp.org/html/33-7400972\91b1758d-3e3b-49fd-8e6e-6d7d2e5d4162.jpg)
Proof
From Lemma 2 and Lemma 3 we have
![](https://www.scirp.org/html/33-7400972\507fd2b1-d9c5-426f-9fdb-6a0634e425bd.jpg)
Using Lemma 4 we get
![](https://www.scirp.org/html/33-7400972\49fe7bfa-4984-4e4f-bc2e-e51761ef541a.jpg)
From (24)
![](https://www.scirp.org/html/33-7400972\4758d704-efdc-4209-9355-2000dded6590.jpg)
Hence
![](https://www.scirp.org/html/33-7400972\8b33b583-f4e6-49f3-9c46-eaa7ead30baa.jpg)
The last equation can be written also as
![](https://www.scirp.org/html/33-7400972\0ea6c3b1-640a-4ab6-8b1a-c3cc861b168c.jpg)
which completes the proof.
4. Natural Levi-Civita Connection on Heisenberg Group [6]
We start with the properties of the Levi-Civita connection on
. For each metric
one has a natural Levi-Civita connection
defined by
(27)
where the Christoffel symbols are defined by the metric (1).They depend linearly on
and are given by using Mathematica program.
![](https://www.scirp.org/html/33-7400972\990b580f-231d-4502-8fba-3d6e6f29426a.jpg)
![](https://www.scirp.org/html/33-7400972\cd8db00b-54cd-41da-866b-a9134eb447c8.jpg)
![](https://www.scirp.org/html/33-7400972\83424c1d-9006-463b-8ef0-a5a9551dd408.jpg)
![](https://www.scirp.org/html/33-7400972\35637e7e-9956-4146-bf12-6c56450ddd49.jpg)
![](https://www.scirp.org/html/33-7400972\060f9078-200f-4781-9c33-bdb640e57d95.jpg)
![](https://www.scirp.org/html/33-7400972\a7547089-f3be-42ac-8170-8e6ec2ed59c5.jpg)
![](https://www.scirp.org/html/33-7400972\68359e0b-878a-4ed8-8c44-5bf36569b0a0.jpg)
![](https://www.scirp.org/html/33-7400972\e1ace610-374a-43b4-a2ba-0ad3d8a09046.jpg)
![](https://www.scirp.org/html/33-7400972\ac8f67f7-9ecd-4f98-8a95-daccd17090db.jpg)
![](https://www.scirp.org/html/33-7400972\61245063-5938-4eee-9834-9f83b3ea492e.jpg)
![](https://www.scirp.org/html/33-7400972\40b69550-be67-4be6-b023-26643f21e8b6.jpg)
![](https://www.scirp.org/html/33-7400972\048d8480-863f-4c1b-9e54-04f517c62762.jpg)
![](https://www.scirp.org/html/33-7400972\10353431-3cda-44d5-97c8-606d0c2e5094.jpg)
![](https://www.scirp.org/html/33-7400972\db3a3e54-494c-48a1-9302-ba4de0bf3052.jpg)
![](https://www.scirp.org/html/33-7400972\bad61f2d-c7c9-460c-9e26-266b52d74b87.jpg)
![](https://www.scirp.org/html/33-7400972\88e230f9-0988-460d-85d6-5b3bd6e7eea4.jpg)
![](https://www.scirp.org/html/33-7400972\3ee5b482-d05f-4dd7-b109-2dbefb7f18d2.jpg)
![](https://www.scirp.org/html/33-7400972\0d5f940f-3418-49eb-83d2-df9806785a5a.jpg)
![](https://www.scirp.org/html/33-7400972\ee2fc686-7e0f-4b22-9d58-cb2cd5ed6fbe.jpg)
![](https://www.scirp.org/html/33-7400972\40cc1aed-ccf3-444e-8a81-4c59525c4d03.jpg)
![](https://www.scirp.org/html/33-7400972\afa337b8-7e52-4def-b53d-4e7d583a390d.jpg)
![](https://www.scirp.org/html/33-7400972\492e8389-5ed5-4d60-ad5d-9bfe12bb7aca.jpg)
![](https://www.scirp.org/html/33-7400972\5564e832-4ce2-448d-902c-becc66a1ae0e.jpg)
![](https://www.scirp.org/html/33-7400972\ad6cae85-b048-4cdd-8ff7-ee8dd03fb10d.jpg)
![](https://www.scirp.org/html/33-7400972\069ba6c2-66d0-4945-8ba1-0453f61ffd77.jpg)
![](https://www.scirp.org/html/33-7400972\093c8921-2095-4a32-bcaf-7baab1410ba2.jpg)
![](https://www.scirp.org/html/33-7400972\a044d43f-3ac7-44d5-9d56-e75b8d9d97e4.jpg)
The following result states tha the Levi-Civita connection with respect to
,
,
is always a linear combination of
and
and hence belongs to the distribution generated by this vectors
Lemma 5
For every
we have
,
![](https://www.scirp.org/html/33-7400972\6a359056-3dc4-45fe-ad15-bf7f6efc4a71.jpg)
![](https://www.scirp.org/html/33-7400972\3161fac8-36ba-4ea5-9d76-f2921a6e599a.jpg)
![](https://www.scirp.org/html/33-7400972\55a4a292-bddf-401a-a64e-411ec94862bd.jpg)
![](https://www.scirp.org/html/33-7400972\582940e5-785c-41b9-9e5c-0439b6676ac8.jpg)
![](https://www.scirp.org/html/33-7400972\be3fa6ca-22c8-458d-b54b-8583f6037e16.jpg)
Proof
For the demonstration, see in [6]
The following proposition shows that the vector fields are geodesic vector fields.
Proposition 2
For every ![](https://www.scirp.org/html/33-7400972\c1221632-1f95-46c5-9df5-e280013ec451.jpg)
![](https://www.scirp.org/html/33-7400972\9f588f77-8d15-4b59-b72f-f5a722ae44a6.jpg)
![](https://www.scirp.org/html/33-7400972\9165bce0-b7d4-4089-b2e5-ef26d512ee8e.jpg)
Proof
See in [6]
5. Conclusions
In this paper, we have proposed an approach for the implementation of the Geodesic method with constraints on Heisenberg Manifolds by including more details with complete proofs that are required for the implementation.
The method has been implemented on Mathematic.8. This implementation extends the rang of applications of the method with more flexibility and rapidity. It has been shown that the Ricci tensor can easily be determined using our implementation.