Next Article in Journal
Postural Complexity during Listening in Young and Middle-Aged Adults
Next Article in Special Issue
How to Secure Valid Quantizations
Previous Article in Journal
On the Choice of the Item Response Model for Scaling PISA Data: Model Selection Based on Information Criteria and Quantifying Model Uncertainty
Previous Article in Special Issue
Entangled Qubit States and Linear Entropy in the Probability Representation of Quantum Mechanics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Symplectic Radon Transform and the Metaplectic Representation

by
Maurice A. de Gosson
Faculty of Mathematics (NuHAG), University of Vienna, Oskar-Morgenstern-Platz 1, 1090 Vienna, Austria
Entropy 2022, 24(6), 761; https://doi.org/10.3390/e24060761
Submission received: 19 April 2022 / Revised: 13 May 2022 / Accepted: 14 May 2022 / Published: 28 May 2022
(This article belongs to the Special Issue Quantum Mechanics and Its Foundations II)

Abstract

:
We study the symplectic Radon transform from the point of view of the metaplectic representation of the symplectic group and its action on the Lagrangian Grassmannian. We give rigorous proofs in the general setting of multi-dimensional quantum systems. We interpret the Radon transform of a quantum state as a generalized marginal distribution for its Wigner transform; the inverse Radon transform thus appears as a “demarginalization process” for the Wigner distribution.

1. Introduction

The idea of using what is today called the “Radon transform” to reconstruct a function from partial data goes back to the 1917 work [1] by the Austrian mathematician Johann Radon. While Radon originally only considered two or three-dimensional systems (in which case it is called the “X-ray transform”), the theory has since then been generalized to arbitrary Euclidean spaces. The applications of the Radon transform to quantum mechanics and optics have been developing rapidly these last years (see the paper [2] by Vogel and Risken). Using an approach from Mancini, Man’ko, and Tombesi [3], several authors [4,5,6,7,8,9,10] study what they call the “symplectic Radon transform” of a mixed quantum state; in a recent paper [11], some of these results are extended to the framework of C algebras. The aim of the present paper is to give a simple rigorous approach to the theory of the symplectic Radon transform in several degrees of freedom. For this, we will use systematically the theory of the metaplectic group as developed in our previous work [12], together with the elementary theory of Lagrangian subspaces of the standard symplectic space. Our main observation is the following: integration of the Wigner transform W ψ ( x , p ) (of a function, or of a state) along the x and p coordinate planes give the correct probability distributions | ψ ( x ) | 2 and | ψ ^ ( p ) | 2 in position and momentum space). However, these are not sufficient to reconstruct the state ψ (this is an aspect of the “Pauli problem”, see Section 4). However, applying a metaplectic transform U ^ associated with a symplectic rotation U to the state ψ transforms the Wigner transform following the rule W U ^ ψ ( x , p ) = W ψ ( U 1 ( x , p ) ) (this is the “symplectic covariance property” of the Wigner transform, which is well-known in harmonic analysis, especially in the Weyl–Wigner–Moyal approach to quantum mechanics). This use of metaplectic transforms leads, by calculating the marginals of W U ^ ψ ( x , p ) along the x and p coordinate planes to infinitely many probability distributions, and these allow the reconstruction of the state. More precisely, we will redefine the Radon transform as being given by the formula
R ψ ( X , A , B ) = det Λ 1 | U ^ A , B ψ ( Λ 1 X ) | 2
where X R n and U ^ A , B is the metaplectic operator associated with the symplectic rotation
U A B = Λ 1 A Λ 1 B Λ 1 B Λ 1 A
where A , B are square matrices with rank ( A , B ) = n and Λ = ( A T A + B T B ) 1 / 2 . We thereafter prove (under suitable conditions on ψ ) the inversion formula
W ψ ( x , p ) = ( 2 π ħ ) 2 n 2 R ψ ( X , A , B ) e i ħ ( X A x B p ) d X d A d B
which reduces for n = 1 to the usual inversion formula found in the literature (the case n = 1 is discussed in Section 2 to motivate the results in the general case). As an illustration, we apply our constructions to the generalized Gaussian states
ψ V , W ( x ) = 1 π ħ n / 4 ( det V ) 1 / 4 e 1 2 ħ ( V + i W ) x 2 .
in Section 3, which gives us the opportunity to shortly discuss the Pauli problem for Gaussians, thus generalizing results in [13].
This paper is structured as follows:
  • In Section 2, we study the “easy case” of systems with one degree of freedom; this allows us to present the main ideas of this paper without using too much new terminology (for instance, any line through the origin of the phase plane is a Lagrangian subspace).
  • In Section 3, we deal with the multi-dimensional case and give a working definition of the symplectic Radon transform from the point of view of operator theory, with a particular emphasis on the metaplectic representation of the symplectic group. We prove an inversion formula allowing the reconstruction of a state from the knowledge of its Radon transform; in particular, the case of Gaussian states is discussed, which allows to give new insight in the old “Pauli problem”.
  • For the reader’s convenience, we have included two appendices where we shortly review the main properties of the metaplectic group (Appendix A) and of the Lagrangian Grassmannian (Appendix B).

2. The Case n = 1

Let ρ ^ be a mixed quantum state with one degree of freedom: ρ ^ is a positive-semidefinite trace class operator on L 2 ( R ) with trace Tr ( ρ ^ ) = 1 . In view of the spectral theorem, there exists a sequence ( ψ j ) with ψ j L 2 ( R ) and a sequence of non-negative numbers ( λ j ) with j λ j = 1 such that ρ ^ = j λ j | ψ j ψ j | where | ψ j ψ j | is the orthogonal projection of L 2 ( R ) onto the ray C ψ j . By definition [14,15], the Wigner distribution of ρ ^ is the convex sum
ρ = j λ j W ψ j
where W ψ j is the usual Wigner transform of ψ j , defined for ψ L 2 ( R ) by
W ψ ( x , p ) = 1 2 π ħ e i ħ p y ψ ( x + 1 2 y ) ψ ( x 1 2 y ) d y .
Recall [14,15,16] that the marginal properties
W ψ ( x , p ) d p = | ψ ( x ) | 2 , W ψ ( x , p ) d x = | ψ ^ ( p ) | 2
make sense provided that ψ and its Fourier transform are, in addition to being square integrable, absolutely integrable: ψ , ψ ^ L 1 ( R ) .
In most texts studying the tomographic picture of quantum mechanics, the symplectic Radon transform of the quantum state ρ ^ is defined by the integral
R ρ ^ ( X , a , b ) = ρ ( x , p ) δ ( X a x b p ) d p d x
where a and b are real numbers, and it is claimed that the following essential reconstruction formula holds
ρ ψ ( x , p ) = 1 2 π ħ R ρ ^ ( X , a , b ) e i ħ ( X a x b p ) d X d a d b .
The following result is simultaneously a rigorous restatement and a justification of these formulas. It will be extended to the case of quantum states with an arbitrary number n of freedom in the forthcoming sections. Among other things, we see that the inverse Radon transform can be viewed as a “demarginalization process” [17] for the Wigner distribution.
We will use the following notation: we set U a , b = a / λ b / λ b / λ a / λ where λ = a 2 + b 2 ; clearly, U a , b is a rotation in the x , p plane.
Theorem 1.
Let ρ ^ be a pure quantum state: ρ = 2 π ħ W ψ for some ψ L 2 ( R ) . We assume that in addition ψ , ψ ^ L 1 ( R ) , which ensures that the marginal properties are satisfied.
(i) 
The Radon transform R ρ ^ ( X , a , b ) is given by the formula
R ρ ^ ( X , a , b ) = λ 1 | U ^ a , b ψ ( λ 1 X ) | 2
where U ^ a , b Mp ( n ) is anyone of the two metaplectic operators covering the rotation U a , b .
(ii) 
The inverse Radon transform is given by the formula:
W ψ ( x , p ) = 1 2 π ħ R ρ ^ ( X , a , b ) e i ħ ( X a x b p ) d X d a d b .
(iii) 
The Radon transform of ψ is given by the line integral
R ρ ^ ( X , a , b ) = W ψ ( z ( t ) ) | z ˙ ( t ) | d t
where t z ( t ) is a parametrization of the straight line a , b X in R 2 with equation a x + b p = X .
Proof. 
It is sufficient to assume that ρ ^ is a pure state, that is, ρ = 2 π ħ W ψ for some ψ . (i) Let us make the change of variables
u v = a / λ b / λ b / λ a / λ x p
in the integral (5). This leads to the expression
R ρ ^ ( X , a , b ) = W ψ ( U a , b 1 ( u , v ) ) δ ( X λ u ) d u d v .
Since δ ( X λ u ) = λ 1 δ ( λ 1 X u ) , this can be rewritten
R ρ ^ ( X , a , b ) = λ 1 W ψ ( λ 1 ( a u b v , b u + a v ) ) δ ( λ 1 X u ) d u d v .
In view of the symplectic covariance property [12,16,18] of the Wigner transform, we have
W ψ ( U a , b 1 ( u , v ) ) = W U ^ a , b ψ ( u , v )
where U ^ a , b is any one of the two metaplectic operators (see the Appendix A) covering U and hence (12) yields
R ρ ^ ( X , a , b ) = λ 1 W U ^ a , b ψ ( λ 1 X , v ) δ ( λ 1 X u ) d u d v = λ 1 W U ^ a , b ψ ( λ 1 X , v ) d v ,
hence Formula (7) using the marginal properties (4). (ii) Let us denote A the right-hand side of the equality (8). Using the first marginal property (4), we have
A = λ 1 1 2 π ħ R 3 | U ^ a , b ψ ( λ 1 X ) | 2 e i ħ ( X a x b p ) d X d a d b = λ 1 1 2 π ħ R 4 W U ^ a , b ψ ( λ 1 X , P ) e i ħ ( X a x b p ) d X d P d a d b .
Replacing X with λ X and using the symplectic covariance property (13), we get
A = 1 2 π ħ R 4 W W U ^ a , b ψ ( X , P ) e i ħ ( λ X a x b p ) d X d P d a d b = 1 2 π ħ R 4 W ψ ( U a , b 1 ( X , P ) ) e i ħ ( λ X a x b p ) d X d P d a d b = 1 2 π ħ R 4 W ψ ( ( a / λ ) X ( b / λ ) P , ( b / λ ) X + ( a / λ ) P ) ) e i ħ ( λ X a x b p ) d X d P d a d b .
Setting Y = ( a / λ ) X ( b / λ ) P and Z = ( b / λ ) X + ( a / λ ) P (and hence λ X = a Y + b Z ), we have d X d P = d Y d Z so that
A = 1 2 π ħ R 4 W ψ ( Y , Z ) ) e i ħ ( a ( Y x ) + b ( Z p ) ) d Y d Z d a d b .
In view of the Fourier inversion formula, written formally as
e i ħ ( a ( Y x ) + b ( Z p ) ) d a d b = 2 π ħ δ ( Y x , Z p )
we thus have
A = W ψ ( x , p ) δ ( Y x , Z p ) d Y d Z = W ψ ( x , p )
which was to be proven. (iii) It is sufficient to show that (9) holds for one parametrization. In view of Formula (11), we have
R ρ ^ ( X , a , b ) = λ 1 W ψ ( U a , b 1 ( u , v ) ) δ ( λ 1 X u ) d u d v
= λ 1 W ψ ( U a , b 1 ( u , v ) ) δ ( λ 1 X u ) d u d v
= W ψ ( U a , b 1 ( λ 1 X , v ) ) d v
that is, since U a , b 1 = a / λ b / λ b / λ a / λ , and replacing v with t,
R ρ ^ ( X , a , b ) = W ψ ( a λ 2 X v λ 1 t , v λ 2 X + a λ 1 t ) d t .
Set now x ( t ) = a λ 2 X v λ 1 t and p ( t ) = v λ 2 X + a λ 1 t . Then, a x ( t ) + b p ( t ) = X and x ˙ ( t ) 2 + p ˙ ( t ) 2 = 1 ; hence, (17) implies (9). □
Remark 1.
The first part of the theorem above uses the physicist’s definition (5) and can thus be taken as a mathematically correct redefinition of the Radon transform. We will exploit this fact in the next section.

3. The Multivariate Case

For n 1 , we consider R 2 n T R n equipped with its standard symplectic structure, which is defined by
σ ( z , z ) = J z · z , J = 0 I I 0 .
We denote by Sp ( n ) the symplectic group of ( R 2 n , σ ) and by M p ( n ) its unitary representation of its double cover (the metaplectic group; see Appendix A).

3.1. Definitions

Let A , B be two real square n × n matrices with such that A T B = B A T and rank ( A , B ) = n . Setting
M A B = A B B A , Λ = ( A T A + B T B ) 1 / 2
and noting that A T A + B T B is invertible, we have the factorization
M A B = Λ 0 0 Λ Λ 1 A Λ 1 B Λ 1 B Λ 1 A
where
U A , B = Λ 1 A Λ 1 B Λ 1 B Λ 1 A U ( n )
is a symplectic rotation. Note that its inverse is
U A , B 1 = A T Λ 1 B T Λ 1 B T Λ 1 A T Λ 1 .
Definition 1.
The symplectic Radon transform of ψ is the transformation
R ψ ( · , A , B ) : L 2 ( R n ) L 1 ( R n )
defined, for ψ L 2 ( R n ) by
R ψ ( X , A , B ) = det Λ 1 | U ^ A , B ψ ( Λ 1 X ) | 2
where Λ = ( A T A + B T B ) 1 / 2 and the unitary operator
U ^ A , B : L 2 ( R n ) L 2 ( R n )
is any of the two elements of Mp ( n ) covering the symplectic rotation U A , B .
That R ψ ( · , A , B ) maps L 2 ( R n ) into L 1 ( R n ) is clear, since U ^ A , B is unitary on L 2 ( R n ) .
When det B 0 , the metaplectic operator U ^ A , B is defined by
U ^ A , B ψ ( x ) = 1 2 π ħ n / 2 i m n / 2 | det B 1 | R n e i ħ A ( x , x ) ψ ( x ) d x
A ( x , x ) = 1 2 A B 1 x · x B 1 x · x + 1 2 B 1 A x · x
after one has made a choice of the integer m modulo 4 (see Appendix A).
Recall that in the multi-dimensional case, the Wigner transform W ψ of ψ L 2 ( R n ) is given by the integral
W ψ ( x , p ) = 1 2 π ħ n R n e i ħ p y ψ ( x + 1 2 y ) ψ ( x 1 2 y ) d y
and that the marginal properties generalizing (4) hold
R n W ψ ( x , p ) d p = | ψ ( x ) | 2 , R n W ψ ( x , p ) d x = | ψ ^ ( p ) | 2
for ψ , ψ ^ L 1 ( R n ) ; the Fourier transform ψ ^ = F ψ is here given by
F ψ ( p ) = 1 2 π ħ n / 2 R n e i ħ p · x ψ ( x ) d x .

3.2. The Radon Inversion Formula

We have the following straightforward generalization of the inversion result (ii) in Theorem 1:
Theorem 2.
Viewing A and B as elements of R n 2 , we have
W ψ ( x , p ) = ( 2 π ħ ) 2 n 2 R n ( 2 n + 1 ) R ψ ( X , A , B ) e i ħ ( X A x B p ) d X d A d B
where d A = i , j = 1 n d a i j and d B = i , j n d b i j if A = ( a i j ) 1 i , j n , B = ( b i j ) 1 i , j n .
Proof. 
It goes exactly as the proof of Theorem 1 (ii). Let us denote by W the integral in the right-hand side of (26). We have, using the first marginal condition (24),
W = det Λ 1 | U ^ A , B ψ ( Λ 1 X ) | 2 e i ħ ( X A x B p ) d X d A d B = det Λ 1 W ( U ^ A , B ψ ) ( Λ 1 X , P ) e i ħ ( X A x B p ) d P d X d A d B = W ( ψ ) ( U A , B 1 ( X , P ) ) e i ħ ( Λ X A x B p ) d P d X d A d B = W ( ψ ) ( U A , B 1 ( X , P ) ) e i ħ ( Λ X A x B p ) d P d X d A d B ,
that is, by Formula (19):
W = W ( ψ ) ( A T Λ 1 X B T Λ 1 P , B T Λ 1 X + A T Λ 1 P ) ) × e i ħ ( Λ X A x B p ) d P d X d A d B .
Setting Y = A T Λ 1 X B T Λ 1 P and Z = B T Λ 1 X + A T Λ 1 P ) , we have d Y d Z = d P d X and hence
W = W ψ ( Y , Z ) e i ħ ( A ( Y x ) + B ( Z p ) ) d Y d Z d A d B .
Integration of the exponential with respect to the variables A and B yields
R 2 n 2 e i ħ ( A ( Y x ) + B ( Z p ) ) d A d B = ( 2 π ħ ) 2 n 2
and hence
W = ( 2 π ħ ) 2 n 2 W ψ ( Y , Z ) δ ( Y x , Z p ) | d Y d Z = ( 2 π ħ ) 2 n 2 W ψ ( x , p )
which was to be proven. □

3.3. Interpretation as Generalized Marginals

Let us return to the definition (20) of the Radon transform:
R ψ ( X , A , B ) = det Λ 1 | U ^ A , B ψ ( Λ 1 X ) | 2 .
If we choose A = I and B = 0 , this reduces to the formula
R ψ ( X , I , 0 ) = | ψ ( X ) | 2 = R n W ψ ( X , P ) d P ;
similarly, if A = 0 and B = I , we get
R ψ ( X , A , B ) = | F ψ ( X ) | 2 = R n W ψ ( P , X ) d X .
definition (20) reduces to the formula
R ψ ( X , A , B ) = | U ^ A , B ψ ( X ) | 2
showing that the Radon transform is essentially a margin property for the “rotated” Wigner transform of ψ . In fact, we can view R ψ ( X , A , B ) as the surface integral of the Wigner transform on the (affine) Lagrangian plane (see Appendix B)
A , B X = { ( x , p ) R 2 n : A x + B p = X } .
In view of the marginal properties and the symplectic covariance of the Wigner transform, followed by the change of variables P Λ 1 P , we have
| U ^ A , B ψ ( Λ 1 X ) | 2 = R n W ( U ^ A , B ψ ) ( Λ 1 X , P ) d P = R n W ψ ( U A , B 1 ( Λ 1 X , P ) ) d P = det Λ R n W ψ ( U A , B 1 ( Λ 1 X , Λ 1 P ) ) d P ,
that is, explicitly,
| U ^ A , B ψ ( Λ 1 X ) | 2 = det Λ × R n W ψ ( A T Λ 1 X B T Λ 1 P , B T Λ 1 X + A T Λ 1 P ) d P .
Using the multi-parametrization
X ( P ) = A T Λ 1 X B T Λ 1 P P ( P ) = B T Λ 1 X + A T Λ 1 P
we have
R ψ ( X , A , B ) = R n W ψ ( X ( P ) , P ( P ) ) d P
we have A X ( P ) + B P ( P ) = X , so we can interpret the formula above as a surface integral
| U ^ A , B ψ ( Λ 1 X ) | 2 = A , B X W ψ ( Z ) d μ ( Z )
where d μ ( Z ) is the Lebesgue measure on A , B X .

4. Radon Transform of Generalized Gaussians

4.1. Generalized Gaussians

By generalized (centered) Gaussian, we mean a function
ψ V , W ( x ) = 1 π ħ n / 4 ( det V ) 1 / 4 e 1 2 ħ ( V + i W ) x 2 .
where V and W are real symmetric n × n matrices with V > 0 (i.e., positive definite). Such centered Gaussians are generalizations of the usual squeezed coherent states appearing in the physical literature; see [12,18]. The function ψ V , W is normalized to unity: | | ψ V , W | | L 2 = 1 and its Wigner transform is given by
W ψ V , W ( z ) = 1 π ħ n e 1 ħ G z · z
where G is the symmetric and symplectic matrix
G = V + W V 1 W W V 1 V 1 W V 1 .
That G Sp ( n ) easily follows from the observation that G = S T S where
S = V 1 / 2 0 V 1 / 2 W V 1 / 2
clearly is symplectic. Let Gauss 0 ( n ) be the set of all centered Gaussians (28); a central (and often implicitly used) result is that the metaplectic group Mp ( n ) acts transitively on Gauss 0 ( n ) : we have an action
Mp ( n ) × Gauss 0 ( n ) Gauss 0 ( n )
which is totally described by the action of S ^ Mp ( n ) on the fiducial coherent state [18]
ϕ 0 ( x ) = 1 π ħ n / 4 e 1 2 ħ | x | 2
whose Wigner transform is
W ϕ 0 ( z ) = 1 π ħ n e 1 ħ | z | 2 .
Using the symplectic covariance of the Wigner transform, it is not difficult to prove that if S ^ covers the symplectic matrix S = A B C D , then
S ^ ϕ 0 ( x ) = e i γ ψ V , W ( x )
where γ R and
U = ( A A T + B B T ) 1 V = ( C A T + D B T ) ( A A T + B B T ) 1 .

4.2. The Radon Transform of ψ V , W

Let us calculate the Radon transform
R ψ U , V ( X , A , B ) = ( det Λ 1 ) | U ^ A , B ψ V , W ( Λ 1 X ) | 2
of ψ V , W using Formula (20). For this, we have to determine U ^ A , B ψ V , W , where U ^ A , B Mp ( n ) covers the symplectic rotation
U A , B = Λ 1 A Λ 1 B Λ 1 B Λ 1 A .
The most natural (and easiest) way to determine U ^ A , B ψ V , W is to use the symplectic covariance formula
W ( U ^ A , B ψ V , W ) ( z ) = W ψ V , W ( U A , B 1 z )
of the Wigner transform; it yields, taking (29) into account and using the relation U A , B 1 = U A , B T
W ( U ^ A , B ψ V , W ) ( z ) = 1 π ħ n e 1 ħ U A , B G U A B T z · z .
To explicitly determine G = U A , B G U A B T is a rather lengthy (though straightforward) calculation; however, since we have an action (32), we know (by (30)) that there will exist V and W such that
G = V + W V 1 W W V 1 V 1 W V 1
corresponding to the Gaussian
U ^ A , B ψ V , W = ψ V , W ( x ) = 1 π ħ n / 4 ( det V ) 1 / 4 e 1 2 ħ ( V + i W ) x 2 .
Now, since it is only the (squared) modulus of U ^ A , B ψ V , W which appears in (35), we will have
| U ^ A , B ψ V , W | 2 = 1 π ħ n / 2 ( det V ) 1 / 2 e 1 2 ħ V x 2
so that it suffices to determine V , whose inverse is the lower right block of G . This is easily done using the relation G = U A , B G U A B T and one finds, after a few calculations and simplifications,
V = Λ B V B T + ( A B W ) V 1 ( A B W ) T 1 Λ .
Summarizing, we have, after insertion in (35),
R ψ V , W ( X , A , B ) = C × exp 1 2 ħ B V B T + ( A B W ) V 1 ( A B W ) T 1 X · X
C = 1 π ħ n / 2 ( det V ) 1 / 2 det B V B T + ( A B W ) V 1 ( A B W ) T 1 .
Suppose for instance ψ V , W is a “squeezed coherent state”; then, W = 0 and ψ V = ψ V , 0 is
ψ V ( x ) = 1 π ħ n / 4 ( det V ) 1 / 4 e 1 2 ħ V x 2 .
Its Radon transform is then
R ψ V ( X , A , B ) = C exp 1 2 ħ ( B V B T + A V 1 A T ) 1 X · X
C = 1 π ħ n / 2 ( det V ) 1 / 2 det ( B V B T + A V 1 A T ) 1
V = Λ B V B T + A V 1 A T 1 Λ

4.3. Application: The Pauli Problem

When n = 1 , the Formula (38) takes the simple form (replacing A , B , V , W with scalars a , b , v , w )
R ψ v , w ( X , a , b ) = C exp 1 2 ħ b 2 v + ( a b w ) 2 v 1 1 X · X
where
ψ v , w ( x ) = 1 π ħ 1 / 4 v 1 / 4 e 1 2 ħ ( v + i w ) x 2 .
Let us apply this formula to the “Pauli problem”. This problem goes back to the question Pauli asked in [19], whether the probability densities | ψ ( x ) | 2 and | ψ ^ ( p ) | 2 uniquely determine the wavefunction ψ ( x ) . The answer is “no”: consider in fact the Gaussian wavepacket
ψ ( x ) = 1 2 π σ x x 1 / 4 e x 2 4 σ x x e i σ x p 2 ħ σ x x x 2
whose Fourier transform of ψ is given by
ψ ^ ( p ) = e i γ 1 2 π σ p p 1 / 4 e p 2 4 σ p p e i σ x p 2 ħ σ p p p 2
where γ is an unessential constant real phase. Thus,
| ψ ( x ) | 2 = 1 2 π σ x x 1 / 2 e x 2 2 σ x x , | ψ ^ ( p ) | 2 = 1 2 π σ p p 1 / 2 e p 2 2 σ p p
and these relations imply the knowledge of σ x x and of σ p p but not of the covariance σ x p (the latter can actually be determined up to a sign using the fact that ψ saturates the Robertson–Schrödinger uncertainty principle: we have σ x x σ p p σ x p 2 = 1 4 ħ 2 ). Let us calculate the Radon transform of ψ using Formula (43). We have here v = ħ / 2 σ x x , w = σ x p / σ x x ; hence, Formula (43) becomes
R ψ v , w ( X , a , b ) = C exp 1 2 ħ b 2 ħ 2 σ x x + a + b σ x p σ x x 2 2 σ x x ħ 1 X 2 .
Notice that, as expected,
R ψ v , w ( X , 1 , 0 ) = 1 2 π σ x x 1 / 2 e X 2 4 σ x x = | ψ ( x ) | 2
and, using the relation σ x x σ p p σ x p 2 = 1 4 ħ 2 ,
R ψ v , w ( X , 0 , 1 ) = 1 2 π σ p p 1 / 2 e X 2 2 σ p p = | ψ ^ ( X ) | 2 .
These relations show why we cannot recover the state ψ using the two radon transforms R ψ v , w ( X , 1 , 0 ) and R ψ v , w ( X , 0 , 1 ) : none of them allow us to determine the covariance σ x p . However, it suffices with one clever choice of the parameters a and b in Formula (47). Suppose indeed that we have measured, for some values of the parameters a and b, the positive quantity
K = b 2 ħ 2 σ x x + a + b σ x p σ x x 2 2 σ x x ħ
appearing in the exponent of (47). Assuming that the variances σ x x and σ p p are known, we can find the covariance σ x p as follows: viewing (48) as a quadratic equation in the unknown σ x p , we demand that this equation has exactly one real root. This imposes the relation
K = b 2 ħ 2 σ x x
and reduces (48) to the equation
a + b σ x p σ x x = 0
from which σ x p is unambiguously determined.

5. Concluding Remarks

In this paper, we outlined a novel approach to the symplectic Radon transform, which we believe is conceptually very simple once one has realized the fundamental role played in quantum mechanics by the metaplectic representation of the symplectic group. As we have discussed elsewhere [20] a few years ago, this is the shortest bridge between classical (Hamiltonian) mechanics, and its refinement, quantum mechanics. This being said, our approach is somewhat sketchy, since we have not characterized the classes of functions (or states) to which we can apply the radon transform, limiting ourselves, for simplicity, to the square integrable case. It is however well-known (at least by people belonging to the harmonic analysis community) that there are function spaces invariant under metaplectic transformations which are larger than the space of square integrable functions. We are, among other possibilities, thinking about Feichtinger’s modulation spaces [21], which are a very flexible tool for creating quantum states, and which can be used together with Shubin’s pseudodifferential calculus [22] to extend the theory (one could, for instance, envisage the reconstruction of general observables along these lines).
We will definitely come back to these topics in a near future.

Funding

This work has been financed by the Grant P 33447 N of the Austrian Research Foundation FWF.

Data Availability Statement

Not applicable.

Acknowledgments

The author wishes to thank the referees for useful remarks and suggestions, and acknowledge Open Access Funding by the Austrian Science Fund (FWF).

Conflicts of Interest

The author declares no conflict of interest.

Appendix A. The Metaplectic Group Mp(n)

For a detailed study of the metaplectic group Mp ( n ) see [12,23]. For a rather “soft” (but still rigorous) approach, see [18].
Let S = A B C D be a real 2 n × 2 n matrix, where the “ blocks” A , B , C , D are n × n matrices. Let J = 0 I I 0 the standard symplectic matrix. We have S Sp ( n ) (the symplectic group) if and only S J S T = S T J S = J . These relations are equivalent to any of the two sets of conditions
A T C , B T D are symmetric , and A T D C T B = I
A B T , C D T are symmetric , and A D T B C T = I .
One says that S is a free symplectic matrix if B is invertible, i.e., det B 0 . To a free symplectic matrix is associated a generating function: it is the quadratic form
A ( x , x ) = 1 2 D B 1 x · x B 1 x · x + 1 2 B 1 A x · x .
The terminology comes from the fact that the knowledge of A ( x , x ) uniquely determines the free symplectic matrix S: we have
x p = A B C D x p p = x A ( x , x ) p = x A ( x , x )
as can be verified by a direct calculation (the quadratic form A is thus a generating function of the free symplectic matrix S).
Now, to every free symplectic matrix S A , we associate two operators S ^ A , m by the formula
S ^ A , m ψ ( x ) = 1 2 π ħ n / 2 i m n / 2 | det B 1 | e i ħ A ( x , x ) ψ ( x ) d n x
where m (“Maslov index” [12]) corresponds to a choice of argument for det B 1 : m = 0   mod 2 if det B 1 > 0 and m = 1 mod 2 if det B 1 < 0 (m is defined modulo 4). It is not difficult to prove that the generalized Fourier transforms S ^ A , m are unitary operators on L 2 ( R n ) . These operators generate a group, the metaplectic group Mp ( n ) . One shows that every S ^ Mp ( n ) can be written (non uniquely) as a product S ^ A , m S ^ A , m . This group is a double covering of Sp ( n ) , the covering projection being defined by
π Mp : Mp ( n ) Sp ( n ) , π Mp ( S ^ A , m ) = S A .

Appendix B. The Lagrangian Grassmannian

A linear subspace of the symplectic space ( R 2 n , σ ) is called a Lagrangian subspace (or plane) if it is maximally isotropic for the skew orthogonality relation σ ( z , z ) = 0 . It must thus have dimension dim = n and σ vanishes identically on . The set of all Lagrangian planes in R 2 n is called the Lagrangian Grassmannian and is denoted by Lag ( n ) . The symplectic group Sp ( n ) acts transitively on Lag ( n ) ; in fact, the action
Sp ( n ) × Lag ( n ) Lag ( n )
thus defined induces a transitive action
U ( n ) × Lag ( n ) Lag ( n )
where U ( n ) is the image in Sp ( n ) of the unitary group U ( n , C ) by the monomorphism ι : u U defined by U ( x , p ) = ( x , p ) if u ( x + i p ) = x + i p . In matrix notation
ι ( A + i B ) = A B B A .
The conditions (A1) and (A2) are here equivalent to
A T B is symmetric , and A T A + B T B = I
A B T is symmetric , and A A T + B B T = I .
The elements of U ( n ) are the symplectic rotations of ( R 2 n , ω ) :
U ( n ) = Sp ( n ) O ( 2 n , R )
and the transitivity of the action (A7) follows: let ( , ) Lag ( n ) × Lag ( n ) and choose two bases { e 1 , , e n } and { e 1 , , e n } of and , respectively. Then, { e 1 , , e n ; J e 1 , , J e n } and { e 1 , , e n ;   J e 1 , , J e n } are both symplectic and orthogonal bases of ( R 2 n , σ ) . The automorphism U of R 2 n taking the first basis to the second is thus in Sp ( n ) O ( 2 n , R ) and we have = U .
Lemma A1.
Let Lag ( n ) . There exist real n × n matrices A , B with rank ( A , B ) = n and A T B = B T A and A B T = B A T such that
= { ( x , p ) R 2 n : A x + B p = 0 } .
Proof. 
It is a straightforward consequence of the transitivity of the action (A7) of U ( n ) on Lag ( n ) . □
Note that in particular, the equation p = M x defines a Lagrangian plane such that ( 0 × R n ) = 0 if and only if M is symmetric and det M 0 .

References

  1. Radon, J. On the Determination of Functions from their Integrals along Certain Manifolds. Ber. Sachsige Akad. Wiss. 1917, 69, 262–277. (In German) [Google Scholar] [CrossRef] [PubMed]
  2. Vogel, K.; Risken, H. Determination of quasiprobability distributions in terms of probability distributions for the rotated quadrature phase. Phys. Rev. A 1989, 40, 2847. [Google Scholar] [CrossRef] [PubMed]
  3. Mancini, S.; Man’ko, V.I.; Tombesi, P. Different realizations of the tomographic principle in quantum state measurement. J. Mod. Opt. 1997, 44, 2281–2292. [Google Scholar] [CrossRef]
  4. Asorey, M.; Facchi, P.; Man’ko, V.I.; Marmo, G.; Pascazio, S.; Sudarshan, E.C.G. Generalized quantum tomographic maps. Phys. Scr. 2012, 85, 85065001. [Google Scholar] [CrossRef] [Green Version]
  5. Facchi, P.; Ligabo, M.; Pascazio, S. On the inversion of the Radon transform: Standard versus M2 approach. J. Mod. Opt. 2010, 57, 239–243. [Google Scholar] [CrossRef]
  6. Chernega, V.; Man’ko, O.V.; Man’ko, V.I. Entangled Qubit States and Linear Entropy in the Probability Representation of Quantum Mechanics. Entropy 2022, 24, 527. [Google Scholar] [CrossRef] [PubMed]
  7. Ibort, A.; Man’ko, V.I.; Marmo, G.; Simoni, A.; Ventriglia, F. An introduction to the tomographic picture of quantum mechanics. Phys. Scr. 2009, 79, 065013. [Google Scholar] [CrossRef] [Green Version]
  8. Ibort, A.; Man’ko, V.I.; Marmo, G.; Simoni, A.; Ventriglia, F. On the tomographic picture of quantum mechanics. Phys. Lett. A 2010, 374, 2614–2617. [Google Scholar] [CrossRef] [Green Version]
  9. López-Saldívar, J.A.; Man’ko, M.A.; Man’ko, V.I. Measurement of the Temperature Using the Tomographic Representation of Thermal States for Quadratic Hamiltonians. Entropy 2021, 23, 1445. [Google Scholar] [CrossRef] [PubMed]
  10. Man’ko, O.V.; Man’ko, V.I. Probability Representation of Quantum States. Entropy 2021, 23, 549. [Google Scholar] [CrossRef] [PubMed]
  11. Ibort, A.; López-Yela, A. Quantum tomography and the quantum Radon transform. Inverse Probl. Imaging 2021, 15, 893–928. [Google Scholar] [CrossRef]
  12. de Gosson, M. Symplectic Geometry and Quantum Mechanics; Series “Operator Theory: Advances and Applications”; Birkhäuser: Basel, Switzerland, 2006; Volume 166. [Google Scholar]
  13. de Gosson, C.; Gosson, M.d. The Pauli Problem for Gaussian Quantum States: Geometric Interpretation. Mathematics 2021, 9, 2578. [Google Scholar] [CrossRef]
  14. de Gosson, M. Quantum Harmonic Analysis of the Density Matrix. Quanta 2018, 7, 74–110. [Google Scholar] [CrossRef] [Green Version]
  15. de Gosson, M. Quantum Harmonic Analysis: An Introduction; Volume 4 in the Series Advances in Analysis and Geometry; De Gruyter: Berlin, Germany, 2021. [Google Scholar]
  16. de Gosson, M. The Wigner Transform, World Scientific; Series Advanced Texts in Mathematics; World Scientific: Singapore, 2017. [Google Scholar]
  17. Park, J.; Lu, Y.; Lee, J.; Shen, Y.; Zhang, K.; Zhang, S.; Zubairy, M.S.; Kim, K.; Nha, H. Revealing nonclassicality beyond Gaussian states via a single marginal distribution. Proc. Natl. Acad. Sci. USA 2017, 114, 891–896. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Littlejohn, R.G. The semiclassical evolution of wave packets. Phys. Rep. 1986, 138, 193–291. [Google Scholar] [CrossRef] [Green Version]
  19. Pauli, W. General Principles of Quantum Mechanics; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2012; [original title: Prinzipien der Quantentheorie; Handbuch der Physik: V.5.1, 1958]. [Google Scholar]
  20. de Gosson, M.; Hiley, B. Imprints of the Quantum World in Classical Mechanics. Found. Phys. 2011, 41, 1415–1436. [Google Scholar] [CrossRef] [Green Version]
  21. Gröchenig, K. Foundation of Time-Frequency Analysis; Birkhäuser: Boston, MA, USA, 2001. [Google Scholar]
  22. Shubin, M.A. Pseudodifferential Operators and Spectral Theory; Springer: Berlin/Heidelberg, Germany, 1987. [Google Scholar]
  23. Leray, J. Lagrangian Analysis and Quantum Mechanics, a Mathematical Structure Related to Asymptotic Expansions and the Maslov Index; MIT Press: Cambridge, UK, 1982. [Google Scholar]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

de Gosson, M.A. Symplectic Radon Transform and the Metaplectic Representation. Entropy 2022, 24, 761. https://doi.org/10.3390/e24060761

AMA Style

de Gosson MA. Symplectic Radon Transform and the Metaplectic Representation. Entropy. 2022; 24(6):761. https://doi.org/10.3390/e24060761

Chicago/Turabian Style

de Gosson, Maurice A. 2022. "Symplectic Radon Transform and the Metaplectic Representation" Entropy 24, no. 6: 761. https://doi.org/10.3390/e24060761

APA Style

de Gosson, M. A. (2022). Symplectic Radon Transform and the Metaplectic Representation. Entropy, 24(6), 761. https://doi.org/10.3390/e24060761

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