Here, we present a general framework for computing the infrared anharmonic vibrational spectra of polyatomic molecules using the Born-Oppenheimer molecular dynamics (BOMD) with the software PyRAMD. To account for the nuclear quantum effects, we suggest using the simplified Wigner sampling (SWS) approach [\textit{Phys. Chem. Chem. Phys.}, 2023, \textbf{25}, 18406-18423] simultaneously coupled with the Andersen and Berendsen thermostats. We propose a new criterion for selecting the parameter of the SWS based on the molecules' harmonic vibrational frequencies and usage of the large-time-step blue shift correction, allowing for a decrease in computational expenses. For the Fourier transform of the dipole moment autocorrelation function, we propose to use the regularized least-squares analysis, which allows us to obtain higher frequency resolution than the direct application of fast Fourier transform. Finally, we suggest the usage of the pre-parameterized scaling factors for the IR spectra from BOMD, also providing the scaling factors for the spectra at the BLYP-D3(BJ)/6-31G, PBE-D3(BJ)/6-31G, and PBEh-3c levels of theory.

Ultrafast pump-probe spectroscopic studies allow for deep insights into the mechanisms and timescales of photophysical and photochemical processes. Extracting valuable information from these studies, such as reactive intermediates' lifetimes and coherent oscillation frequencies, is an example of the inverse problems of chemical kinetics. This article describes a consistent approach for solving this inverse problem that avoids the common obstacles of simple least-squares fitting that can lead to unreliable results. The presented approach is based on the regularized Markov Chain Monte Carlo sampling for the strongly nonlinear parameters, allowing for a straightforward solution of the ill-posed nonlinear inverse problem. The software to implement the described fitting routine is introduced and the numerical examples of its application are given. We will also touch on critical experimental parameters, such as the temporal overlap of pulses and cross-correlation time and their connection to the minimal reachable time resolution.