Abstract

A procedure is described that can be used to reconstruct the quantum state of a molecular ensemble from time-dependent internuclear probability density functions determined by time-resolved electron diffraction. The procedure makes use of established techniques for evaluating the density matrix and the phase-space joint probability density, that is, the Wigner function. A novel expression for describing electron diffraction intensities in terms of the Wigner function is presented. An approximate variant of the method, neglecting the off-diagonal elements of the density matrix, was tested by analyzing gas electron diffraction data for N2 in a Boltzmann distribution and TRED data obtained from the 193 nm photodissociation of CS2 to carbon monosulfide, CS, at 20, 40, and 120 ns after irradiation. The coherent changes in the nuclear subsystem by time-resolved electron diffraction method determine the fundamental transition from the standard kinetics to the dynamics of the phase trajectory of the molecule and the tomography of molecular quantum state.

1. Introduction

In accordance with basic quantum principles, the state of an individual molecule cannot be determined experimentally [1]. However, for an ensemble of similarly prepared systems, it is possible to determine their state operator, the so-called density matrix. Knowing the state of a system means having the maximum possible information about all physical quantities of interest available [2]. The density matrix and the joint phase-space probability density, or Wigner function, [3, 4] have a one-to-one correspondence [5] that describes the maximal statistical information available. Thus, in the following text, when the term molecular quantum state is used, we mean the quantum state of an ensemble of similarly prepared molecular species.

In 1933 it was demonstrated by Freenberg [6] (see also [1, page 71]) that, in principle, a pure quantum state can be reconstructed from the time-dependent coordinate probability density and its derivative . It was shown by Weigert [7] that a pure quantum state may also be reconstructed by measuring the distribution at time and monitoring its evolution after a short time interval ; that is, .

For mixed quantum states, the method of optical homodyne tomography to measure the Wigner function (and density matrix) was first demonstrated by Smithey et al. [8, 9] for both vacuum and squeezed vacuum states of a single spatial-temporal mode in an applied electromagnetic field. (A number of investigations into the preparation and measurement of the quantum state of light may be found in [10].)

In experiments of Kurtsiefer et al. [11, 12], it is demonstrated that Wigner’s function can be reconstructed from double-slit experiments involving a coherent beam of helium atoms. Their result [11, 12] shows that the joint phase-space probability density , where is the momentum, exhibits negative areas indicating purely nonclassical behavior, such as what might be expected of a particle in an atomic interferometer when it is in a linear superposition of states at two separate locations (see also comments in [13]).

In studies of molecular states, measurements of time-dependent fluorescence spectra make it possible to reconstruct Wigner’s function for single molecular vibrational modes excited by short optical pulses. This technique is referred to as Molecular Emission Tomography (MET) [14] and was used by Dunn et al. in studies of electronically excited Na2 wave packets, moving in a suggested harmonic potential.

In the current paper, it is described how time-resolved electron diffraction (TRED) may be used in tomographic studies of the molecular state. TRED [1517] provides experimental observations having a direct correspondence to the molecular time-dependent probability density function . In TRED experiments, a pulsed electron source is used to probe an ensemble of similarly prepared species throughout a sequence in time. The diagnostic electron pulses and the concomitant exciting laser pump pulses are synchronized in such way that stroboscopic diffraction pictures are obtained. In this way, time is introduced as a new variable into the otherwise conventional diffraction technique [1517]. Thus, reacting laser-excited species, transient states of chemical reactions, or molecular wave packets can be investigated [1520]. The dynamics of wave packets, as the coherent superposition of quantum states created by short, phase-controlled optical pulses [2124], can also be probed. Coherent nuclear dynamics can be followed in this way [15, 2527] because the elastic scattering of fast electrons proceeds on an ultrashort (ultimately attosecond) time scale and represents a nondestructive probing technique [15, 25, 26].

2. Theory

2.1. General Aspects

TRED experiments yield time-dependent molecular scattering intensity functions, [15, 25, 26], where , is the momentum vector of the scattered electrons, and is the scattering angle: In (1), is the incident electron beam profile, is the scattering distance, are the scattering functions [28], and is the dot product of the scattering vector and internuclear distance vector . Integration in (1) and in the following equations is over the infinite space if not indicated otherwise.

In classical mechanics, there are no equations that describe the evolution of the probability densities (or , where is momentum); only joint probability densities, , can be expressed by the Liouville equation. Therefore, corresponding quantum equations for , as used in (1), and cannot exist. However, the well known Wigner-Liouville equation [4] can be used to describe the evolution of the Wigner function, . In order to derive an equation that expresses the time-dependent molecular intensities in terms of the Wigner function, the marginal properties of [3, 4] must be considered. That is, Hence, (1) can be expressed as follows: Equation (4) is the most general representation of TRED molecular intensities in terms of the Wigner function. In this representation, can be interpreted as the filtered projection of the Wigner function, with the scattering operator as the filter, modified by the scattering functions . For a number of problems, can be evaluated analytically (see, e.g., the review in [4]) or by solving the Wigner-Liouville equation numerically with appropriate molecular potential functions.

In (1) and (2), . The wave function can be expanded in terms of the orthonormal basis functions in the following way [1, 2]: where is the quantum number of the state energy , is the oscillator frequency, and is the amplitude. Then, (1). can be written as where the are the density matrix elements.

From (6) it is seen that UED molecular intensities depend explicitly on the quantum state of a molecular system. Likewise, the obtained from TRED data, that is, the Fourier transform of (1), depends on both the position and the time and contains implicitly all information on the quantum state of the system, being a projection or “shadow” [29] of the Wigner function. A simple sequence of TRED measurements of the position-dependent probability density provides the necessary information for tomographic reconstruction of the Wigner function. For this purpose, the inverse Radon transformation can be applied [29, 30].

Another possibility for molecular quantum state reconstruction is the evaluation of the density matrix elements, , which have a one-to-one correspondence [5] to Wigner’s function: where is the density matrix [2]:

2.2. Approximate Procedures

When the TRED measurements do not span a complete cycle and the data are incomplete, only the diagonal density matrix elements, , can be evaluated from the probability density averaged over the time interval , with , as shown by Richter and Wünsche [3134], Leonhardt [29, 35], Leonhardt and Raymer [36], and others [3739]. This method uses the superposition principle combined with the structure of Schrödinger’s equation to reveal the quantum state of a molecule from the time-averaged coordinate-dependent probability density function. In general, and for a complete data set, where is the Fourier transform of the time-dependent probability density [29, equation 5.59], the are the normalized regular solutions of the time-independent Schrödinger equation, and the are the linearly independent nonnormalizable irregular second solutions [33, 36, 40].

In the particular case of incomplete data, one can express the diagonal density matrix elements as Time-dependent probability density function averaged over electron pulse profile of duration provides evolution in time of the diagonal density matrix elements in (11).

The terms state, density matrix, or Wigner function are, for most practical purposes, synonymous. Thus, when the available data are incomplete and only the diagonal density matrix elements are used, one has to keep in mind that it is not possible to determine the complete Wigner function. However, it is possible to use the diagonal density matrix elements to calculate an approximate or truncated Wigner function, which contains useful information when the off-diagonal elements are small. The truncated Wigner functions can be used to calculate approximate von Neumann entropies and purity parameters which, nonetheless, may contain meaningful information on the quantum state of a system.

3. Experiment

Instrumentation for TRED developed at the University of Arkansas, spanning the time domain from microseconds to picoseconds, has been described before [41]. A UV laser-driven tantalum photocathode in a Pierce-type [42] electron gun assembly was used for generation of the probe electron pulses. The mechanics of the photoemission process suggest that the electron pulse should follow the temporal profile of the incident optical pulse [43, 44]. In the experiments used for this study, the output from an excimer laser (Questek Model 2260, 193 nm (ArF), pulse duration ~15 ns) was used to generate electron pulses (duration ~15 ns FWHM, ~1010 electrons per pulse [41]) which were accelerated through a potential of ~40 keV. For the current study, carbon disulfide, CS2, target molecules were irradiated with similar 193 nm laser pulses from a second excimer at 50–130 mJ per pulse. Delay times between excitation and electron scattering varied from 20 to 120 ns and the diffraction intensities were recorded online using a phosphor screen-photodiode array detector [41]. The target gases were admitted into the vacuum chamber through a nozzle orifice (I.D. 200 μm) at a pushing pressure of 30–50 Torr. The electron scattering occurred in a volume immediately adjacent to the nozzle exit in which the gas pressure is estimated to be ~10 Torr. Such conditions preclude the formation of clusters, and the number of collisions per particle in the target gas is estimated to be in the range of in 15 ns. Also, time averaging is determined by the electron pulse profile and averaging time  ns. The experimental TRED data sets [45] and sampling functions are presented in Figure 1.

To further illustrate the described procedures, N2 in an equilibrium Boltzmann distribution was chosen as a test case because its potential energy function is well known from high-resolution spectroscopic measurements. Data were recorded throughout the temperature range from 298 K to 773 K, with a continuous electron beam [15, 41]. Experimental data and sampling functions are presented in Figure 2.

4. Results and Discussion

The photochemistry of CS2 has been thoroughly studied in the absorption region from 180 to 210 nm [4547]. At 193 nm, and in single-photon processes, the photodissociation proceeds according to the following two-channel scheme: The structure of excited state of CS2 is believed to be quasilinear, but with different geometrical parameters compared to the ground state; further, predissociation occurs on a time scale of ~1 ps.

In supersonic molecular beams of 193 nm irradiated CS2, it was found [46, 47] that the nascent vibrational distribution of CS is strongly inverted for both the and dissociation channels, exhibiting a bimodal [46] or broad [47] distribution with low occupancy at and , as summarized in Table 1.

To apply the quantum sampling procedure to N2 as a test equilibrium case, position-probability densities were first derived from the diffraction intensities by using the stochastic approach to data analysis [45]. This approach [45] allows one to refine TRED data of laser-excited species without imposing any model restrictions on intramolecular motion. It is an advantage of the procedure that the requisite molecular parameters may be refined from the experimental data without any information on the potential energy surface. The probability density functions obtained in this way were then sampled with the functions [34]: where and are, respectively, the normalized regular and nonnormalizable irregular wave functions referred to above. The sampling functions for the first 3 levels are shown in Figure 2. Calculated density matrix elements for the Boltzmann distribution of N2 at are compared with the experimentally determined values in Table 1.

The truncated Wigner function, together with the associated marginal position- and momentum-dependent probability density functions ((2) and (3)), is shown in Figure 3. The results represent, for the nearly harmonic one-dimensional oscillator at equilibrium, a Gaussian-like function and Gaussian-like marginal position- and momentum-probability densities.

For CS, we derived time-dependent probability density functions at several delay times from the experimental time-dependent molecular intensity functions, using the stochastic approach [15, 45]. To calculate the sampling functions of (14), the Dunham vibrational potential of CS, as determined spectroscopically up to sixth order [48], was used. The regular wave functions were evaluated by solving Schrödinger’s equation for the same Dunham potential [48] using the variational technique of [49].

The values of the diagonal density matrix elements, , obtained by applying (11) to the experimental at different delay times, are represented in Table 1. Truncated Wigner functions for , 40 ns, and 120 ns are shown in Figure 4 and compared with those for the nascent state distribution of [47] (see also Table 1). For and 120 ns the truncated Wigner functions and the corresponding marginal density functions, and ((2) and (3)), are similar. In contrast, these functions exhibit significant changes between , 20, and 40 ns. We interpret these changes in the following way.

Changes in the first 20 ns of the evolution of the system are indicative of rapid equilibration, predominantly by noncollisional intermolecular vibrational energy transfer involving, for example, dipole-dipole interactions. During this time span bimodal distributions, both in position and momentum space, evolve to monomodal and narrower distributions. Between 20 and 40 ns, collision-induced energy transfer becomes more important than collisionless energy transfer processes. Thus, we believe that the changes between 20 and 40 ns indicate electronic to vibrational energy transfer induced by collisions between atoms and molecules. Between 40 and 120 ns, this process shows features of saturation.

It can be inferred from Table 1 that, keeping the foregoing approximations in mind, there is significant depopulation of the first three vibrational levels of CS between 20 and 120 ns, very likely indicating the transfer of vibrational energy to rotational and translational degrees of freedom. This interpretation is supported by the fact that the rotational temperature of CS in this process increases from at 20 ns to 5000 ± 600 K at 40 ns and to 5200 ± 800 K at 120 ns [45]. The truncated Wigner functions at 20, 40, and 120 ns are positive throughout the entire phase space, in agreement with the hypothesis that the molecular quantum state of CS after 20 ns can be characterized as a classical-like statistical mixture, rather than a coherently prepared superposition of the lower vibrational states. This confirms the increasing importance of collisions after 20 ns and is in contrast to the nascent molecular distribution, for which the complete Wigner function exhibits significant negative areas (not visible in Figure 4). Furthermore in this case, the bimodality of the marginal distributions, and , reflects the fact that the 193 nm photodissociation of CS2 proceeds via two different channels, producing and . The symmetric nature of the bimodality shown in Figure 4 is an artifact due to the truncation of the Wigner function.

We have also investigated the processes described above in terms of the von Neumann entropy [29]: where is the trace of the operator and is the density matrix. The von Neumann entropy vanishes for pure states, exceeds zero for mixed states, and is an extensive quantity for nonentangled subsystems; thus, it is a fundamental measure of the preparation impurity of a quantum state [29]. It is also informative to consider the so-called purity parameter, , which represents a lower bound for the von Neumann entropy; that is, : In the current investigation, the diagonal density matrix elements were used to derive approximate estimates for and . In spite of their approximate nature, it is seen from Table 2 that the estimates obtained are in agreement with previous descriptions [45] of the 193 nm photodissociation of CS2. That is, at 20 ns both the approximated von Neumann entropy and purity parameters are smaller than at (nascent CS), indicating the transition from bimodal to monomodal distributions. At times later than 20 ns, both of these parameters increase due to increasing populations of vibrational states above .

5. Conclusion

The material presented above shows that time-resolved electron diffraction measurements offer a new avenue for quantum state reconstruction, yielding the Wigner function from a tomographically complete set of time-dependent position-probability densities. When the TRED data are incomplete, a truncated Wigner function can be determined from the diagonal density matrix elements derived from the time-averaged experimental probability density . Even though the diagonal terms do not describe the full quantum state of a system, they still contain useful information.

Deriving the truncated Wigner function for N2 in an equilibrium Boltzmann distribution, the expected Gaussian-like features were found. In the case of the 193 nm photodissociation of CS2, the analysis indicates collisionless vibrational energy transfer mechanisms in the photogenerated ensemble of nascent CS and a collision-induced electronic to vibrational energy transfer process, from to at later times. Apart from quantum state reconstruction, (4) represents a novel expression for analyzing diffraction intensities directly in terms of the Wigner function.

Molecular quantum state measurements by TRED are complementary to quantum state preparation and control of matter by laser fields. Coherent nuclear dynamics can be followed in this way [15, 2527] because the elastic scattering of fast electrons proceeds on an ultrashort (ultimately attosecond) time scale and represents a nondestructive probing technique. Adequately, new avenues in physical chemistry are in development based on studies of coherent nuclear dynamics [15, 16, 50].

Since the 1980, the scientific world made intensive efforts in order to register a movie about the coherent nuclei dynamics in the molecules, the fast dynamic processes in biological tissues and cells, and the structure dynamics of the solid in nanovolumes in time. The observed coherent changes in the nuclear subsystem by time-resolved electron diffraction method determine the fundamental transition from the standard kinetics to the dynamics of the phase trajectory of the molecule and the tomography of molecular quantum state.

Acknowledgments

Academician of the RAS, A. L. Buchachenko, Professor G. V. Fetisov (Chemistry Department of Moscow Lomonosov State University), Professor L. Schafer (University of Arkansas at Fayetteville, USA), Professor E. A. Ryabov (Institute of Spectroscopy RAS), and Professor V. N. Bagratashvili (IPLIT RAS) are greatly thanked for their helpful discussions. This work was supported by RFBR Grants 12-02-00840-a and 12-02-12048-ofi_m.