Hamiltonian matrix

In mathematics, a Hamiltonian matrix is a {{math|2n}}-by-{{math|2n}} matrix {{mvar|A}} such that {{math|JA}} is symmetric, where {{mvar|J}} is the skew-symmetric matrix

: $J =
\begin{bmatrix}
0_n & I_n \\
-I_n & 0_n \\
\end{bmatrix}$

and {{math|I_n}} is the {{mvar|n}}-by-{{mvar|n}} identity matrix. In other words, {{mvar|A}} is Hamiltonian if and only if {{math|1 = (JA)^T = JA}} where {{math|()^T}} denotes the transpose.{{citation | first = Khakim D. | last = Ikramov | title = Hamiltonian square roots of skew-Hamiltonian matrices revisited | journal = Linear Algebra and its Applications | volume = 325 | year = 2001 | pages = 101–107 | doi = 10.1016/S0024-3795(00)00304-9 | doi-access = free }}.

(Not to be confused with Hamiltonian (quantum mechanics))

Properties

Suppose that the {{math|2n}}-by-{{math|2n}} matrix {{mvar|A}} is written as the block matrix

: $A = \begin{bmatrix} a & b \\ c & d \end{bmatrix}$

where {{mvar|a}}, {{mvar|b}}, {{mvar|c}}, and {{mvar|d}} are {{mvar|n}}-by-{{mvar|n}} matrices. Then the condition that {{math|A}} be Hamiltonian is equivalent to requiring that the matrices {{math|b}} and {{math|c}} are symmetric, and that {{math|1 = a + d^T = 0}}.{{citation | first1 = K. R. | last1 = Meyer | first2 = G. R. | last2 = Hall | title = Introduction to Hamiltonian dynamical systems and the {{math|N}}-body problem | publisher = Springer | year = 1991 | isbn = 0-387-97637-X }}. Another equivalent condition is that {{math|A}} is of the form {{math|1 = A = JS}} with {{math|S}} symmetric.{{rp|34}}

It follows easily from the definition that the transpose of a Hamiltonian matrix is Hamiltonian. Furthermore, the sum (and any linear combination) of two Hamiltonian matrices is again Hamiltonian, as is their commutator. It follows that the space of all Hamiltonian matrices is a Lie algebra, denoted {{math|sp(2n)}}. The dimension of {{math|sp(2n)}} is {{math|2n² + n}}. The corresponding Lie group is the symplectic group {{math|Sp(2n)}}. This group consists of the symplectic matrices, those matrices {{mvar|A}} which satisfy {{math|1 = A^TJA = J}}. Thus, the matrix exponential of a Hamiltonian matrix is symplectic. However the logarithm of a symplectic matrix is not necessarily Hamiltonian because the exponential map from the Lie algebra to the group is not surjective.{{rp|34–36}}{{citation | first = Alex J. | last = Dragt | doi = 10.1196/annals.1350.025 | title = The symplectic group and classical mechanics | journal = Annals of the New York Academy of Sciences | year = 2005 | volume = 1045 | issue = 1 | pages=291–307 | pmid=15980319}}.

The characteristic polynomial of a real Hamiltonian matrix is even. Thus, if a Hamiltonian matrix has {{math|λ}} as an eigenvalue, then {{math|−λ}}, {{math|λ^*}} and {{math|−λ^*}} are also eigenvalues.{{rp|45}} It follows that the trace of a Hamiltonian matrix is zero.

The square of a Hamiltonian matrix is skew-Hamiltonian (a matrix {{mvar|A}} is skew-Hamiltonian if {{math|1 = (JA)^T = −JA}}). Conversely, every skew-Hamiltonian matrix arises as the square of a Hamiltonian matrix.{{citation | first = William C. | last = Waterhouse | authorlink = William C. Waterhouse | doi = 10.1016/j.laa.2004.10.003 | title = The structure of alternating-Hamiltonian matrices | journal = Linear Algebra and its Applications | volume = 396 | year = 2005 | pages = 385–390 | doi-access = free }}.

Extension to complex matrices

As for symplectic matrices, the definition for Hamiltonian matrices can be extended to complex matrices in two ways. One possibility is to say that a matrix {{mvar|A}} is Hamiltonian if {{math|1 = (JA)^T = JA}}, as above. Another possibility is to use the condition {{math|1 = (JA)^* = JA}} where the superscript asterisk ({{math|(⋅)^*}}) denotes the conjugate transpose.{{citation | first1 = Chris | last1 = Paige | first2 = Charles | last2 = Van Loan | authorlink2 = Charles F. Van Loan | title = A Schur decomposition for Hamiltonian matrices | journal = Linear Algebra and its Applications | volume = 41 | year = 1981 | pages = 11–32 | doi = 10.1016/0024-3795(81)90086-0 | doi-access = free }}.

Hamiltonian operators

Let {{mvar|V}} be a vector space, equipped with a symplectic form {{math|Ω}}. A linear map $A : \; V \mapsto V$ is called a Hamiltonian operator with respect to {{math|Ω}} if the form $x, y \mapsto \Omega(A(x), y)$ is symmetric. Equivalently, it should satisfy

: $\Omega(A(x), y) = -\Omega(x, A(y))$

Choose a basis {{math|e₁, …, e_2n}} in {{mvar|V}}, such that {{math|Ω}} is written as $\sum_i e_i \wedge e_{n+i}$ . A linear operator is Hamiltonian with respect to {{math|Ω}} if and only if its matrix in this basis is Hamiltonian.

References

Category:Matrices (mathematics)