The ionization and dissociation of molecules in intense laser fields represent a central topic in attosecond physics. In particular, the quantum entanglement between the emitted photoelectron and the residual molecular ion has garnered significant attention for its potential in ultrafast quantum control. Quantifying this correlation requires the experimental retrieval of the reduced density matrix of the molecular ion's internal degrees of freedom, specifically its vibrational states. However, the complexity of coupled electron-nuclear dynamics makes such measurements challenging.
In this work, we present a theoretical framework for reconstructing the vibrational reduced density matrix of H₂+ using a two-color extreme ultraviolet (XUV)-ultraviolet (UV) pump-probe scheme. Based on the Intense-Field Many-body S-matrix Theory (IMST), we derive the transition amplitude for dissociative ionization, where an XUV pump pulse single-photon ionizes H₂ into H₂+ (1sσ_g), followed by a time-delayed UV probe pulse that excites the ion to a repulsive state, leading to fragmentation. We demonstrate that the Kinetic Energy Release (KER) spectrum of the ionic fragments, as a function of the pump-probe delay, contains sufficient information to retrieve the vibrational density matrix.
Numerical simulations reveal that the scheme successfully reconstructs the diagonal elements (populations) of the low-energy vibrational states. However, significant deviations are observed in the off-diagonal elements (coherences). Systematic analysis identifies two primary sources of error: (1) An inherent algorithmic approximation where the nuclear wave packet evolution is fixed at the pump pulse center, ignoring the temporal distribution of ionization events across the pulse duration; and (2) Numerical grid mismatch, where the uniform frequency grid from Fourier transforming discrete delay times fails to match the non-uniform spacing of molecular vibrational levels, leading to ill-conditioned matrices and amplified noise in off-diagonal terms.
Furthermore, we investigate the dependence of reconstruction fidelity on laser parameters. Reducing the XUV pump wavelength, increasing its intensity, or shortening the UV probe wavelength significantly improves the completeness and accuracy of the retrieved matrix. Conversely, variations in the probe intensity have negligible effects.
This study validates the feasibility of characterizing molecular quantum states via pump-probe spectroscopy and provides practical guidelines for experimental parameter optimization in ultrafast quantum dynamics and quantum information science.