No Arabic abstract
Theoretical calculations are contributing a significantly higher proportion of data to contemporary spectroscopic databases, which have traditionally relied on experimental observations and semi-empirical models. It is now a common procedure to extend calculated line lists to include ro-vibrational transitions between all bound states of the ground electronic state up to the dissociation limit. Advanced ab initio methods are utilized to calculate the potential energy and dipole moment surfaces (PESs and DMSs), and semi-empirical PESs are then obtained by combining ab initio and experimental data. The objective is to reach high accuracy in the calculated transition intensities for all parts of spectrum, i.e. to increase the predictive power of the model. We show that in order to perform this task, one needs, in addition to the standard improvements of the PES and DMS in the spectroscopically accessible regions, to extend the ab initio calculations of the PES towards the united-atom limit along the stretching coordinates. The argument is based on the correlation between the intensities of high-overtone transitions and the repulsive potential wall that has previously been theoretically established for diatomic molecules and is empirically extended here to linear and nonlinear triatomic molecules. We generate partial line lists for water and ozone, and together with an already available line list for carbon dioxide, we derive the normal intensity distribution, which is a direct consequence of this correlation. The normal distribution is not an instrument to compute highly accurate intensities, rather it is a means to analyze the intensities computed by the traditional methods.
We formulate a microcanonical distribution for an arbitrary one-electron triatomic molecule. This distribution can be used to describe the initial state in strongly-driven two-electron triatomic molecules. Namely, in many semiclassical models that describe ionization of two-electron molecules driven by intense infrared laser fields in the tunneling regime initially one electron tunnels while the other electron is bound. The microcanonical distribution presented in this work can be used to describe the initial state of this bound electron.
Using a semi-classical model, we investigate frustrated double ionization (FDI) in $mathrm{D_3^+}$, a two-electron triatomic molecule, when driven by an intense, linearly polarized, near-infrared (800 nm) laser field. We compute the kinetic energy release of the nuclei and find a good agreement between experiment and our model. We explore the two pathways of FDI and show that, with increasing field strength, over-the-barrier ionization overtakes tunnel ionization as the underlying mechanism of FDI. Moreover, we compute the angular distribution of the ion fragments for FDI and identify a feature that can potentially be observed experimentally and is a signature of only one of the two pathways of FDI.
Universal properties of the Coulomb interaction energy apply to all many-electron systems. Bounds on the exchange-correlation energy, inparticular, are important for the construction of improved density functionals. Here we investigate one such universal property -- the Lieb-Oxford lower bound -- for ionic and molecular systems. In recent work [J. Chem. Phys. 127, 054106 (2007)], we observed that for atoms and electron liquids this bound may be substantially tightened. Calculations for a few ions and molecules suggested the same tendency, but were not conclusive due to the small number of systems considered. Here we extend that analysis to many different families of ions and molecules, and find that for these, too, the bound can be empirically tightened by a similar margin as for atoms and electron liquids. Tightening the Lieb-Oxford bound will have consequences for the performance of various approximate exchange-correlation functionals.
Vibrational spectra and wavefunctions of polyatomic molecules can be calculated at low memory cost using low-rank sum-of-product (SOP) decompositions to represent basis functions generated using an iterative eigensolver. Using a SOP tensor format does not determine the iterative eigensolver. The choice of the interative eigensolver is limited by the need to restrict the rank of the SOP basis functions at every stage of the calculation. We have adapted, implemented and compared different reduced-rank algorithms based on standard iterative methods (block-Davidson algorithm, Chebyshev iteration) to calculate vibrational energy levels and wavefunctions of the 12-dimensional acetonitrile molecule. The effect of using low-rank SOP basis functions on the different methods is analyzed and the numerical results are compared with those obtained with the reduced rank block power method introduced in J. Chem. Phys. 140, 174111 (2014). Relative merits of the different algorithms are presented, showing that the advantage of using a more sophisticated method, although mitigated by the use of reduced-rank sum-of-product functions, is noticeable in terms of CPU time.
Including quantum mechanical effects on the dynamics of nuclei in the condensed phase is challenging, because the complexity of exact methods grows exponentially with the number of quantum degrees of freedom. Efforts to circumvent these limitations can be traced down to two approaches: methods that treat a small subset of the degrees of freedom with rigorous quantum mechanics, considering the rest of the system as a static or classical environment, and methods that treat the whole system quantum mechanically, but using approximate dynamics. Here we perform a systematic comparison between these two philosophies for the description of quantum effects in vibrational spectroscopy, taking the Embedded Local Monomer (LMon) model and a mixed quantum-classical (MQC) model as representatives of the first family of methods, and centroid molecular dynamics (CMD) and thermostatted ring polymer molecular dynamics (TRPMD) as examples of the latter. We use as benchmarks D$_2$O doped with HOD and pure H$_2$O at three distinct thermodynamic state points (ice Ih at 150K, and the liquid at 300K and 600K), modeled with the simple q-TIP4P/F potential energy and dipole moment surfaces. With few exceptions the different techniques yield IR absorption frequencies that are consistent with one another within a few tens of cm$^{-1}$. Comparison with classical molecular dynamics demonstrates the importance of nuclear quantum effects up to the highest temperature, and a detailed discussion of the discrepancies between the various methods let us draw some (circumstantial) conclusions about the impact of the very different approximations that underlie them. Such cross validation between radically different approaches could indicate a way forward to further improve the state of the art in simulations of condensed-phase quantum dynamics.