Initial conditions for surface hopping trajectories from the VSCF-Wigner distribution

In trajectory surface hopping (TSH) simulations, initial conditions are typically generated using harmonic Wigner distributions, which assume independent harmonic normal modes. While this assumption fails for anharmonic systems, it remains unclear under which conditions harmonic Wigner sampling becomes unreliable in photochemical simulations and whether anharmonicity alone is a sufficient criterion for reconsidering the use of harmonic Wigner sampling. In the present study, we introduce a sampling strategy based on vibrational self-consistent field (VSCF) theory to construct a VSCF Wigner quasiprobability distribution that incorporates anharmonic effects while retaining mode separability. Analytical expressions are derived in both harmonic and distributed Gaussian bases, enabling the implementation in TSH simulations. The method is applied to malonaldehyde and methyl hydroperoxide, which exhibit moderate and strong anharmonicity, respectively. For malonaldehyde, VSCF-based Wigner sampling and harmonic Wigner sampling yield similar results, indicating that harmonic Wigner sampling remains reliable despite the anharmonicity. In methyl hydroperoxide, where torsional motion significantly influences the excited-state character, VSCF Wigner sampling yields results comparable to those of the quantum thermostat approach while offering a computationally efficient and systematically improvable route to initial-condition sampling.