Quantum dynamics of the CH3 fragment: A curvilinear coordinate system and kinetic energy operators

QuTree: A tree tensor network package

We present QuTree, a C++ library for tree tensor network approaches. QuTree provides class structures for tensors, tensor trees, and related linear algebra functions that facilitate the fast development of tree tensor network approaches such as the multilayer multiconfigurational time-dependent Hartree approach or the density matrix renormalization group approach and its various extensions. We investigate the efficiency of relevant tensor and tensor network operations and show that the overhead for managing the network structure is negligible, even in cases with a million leaves and small tensors. QuTree focuses on providing simple, high-level routines while retaining easy access to the backend to facilitate novel developments. We demonstrate the capabilities of the package by computing the eigenstates of coupled harmonic oscillator Hamiltonians and performing random circuit simulations on a virtual quantum computer.

Eigenstate calculation in the state-averaged (multi-layer) multi-configurational time-dependent Hartree approach

A new approach for the calculation of eigenstates with the state-averaged (multi-layer) multi-configurational time-dependent Hartree (MCTDH) approach is presented. The approach is inspired by the recent work of Larsson [J. Chem. Phys. 151, 204102 (2019)]. It employs local optimization of the basis sets at each node of the multi-layer MCTDH tree and successive downward and upward sweeps to obtain a globally converged result. At the top node, the Hamiltonian represented in the basis of the single-particle functions (SPFs) of the first layer is diagonalized. Here p wavefunctions corresponding to the p lowest eigenvalues are computed by a block Lanczos approach. At all other nodes, a non-linear operator consisting of the respective mean-field Hamiltonian matrix and a projector onto the space spanned by the respective SPFs is considered. Here, the eigenstate corresponding to the lowest eigenvalue is computed using a short iterative Lanczos scheme. Two different examples are studied to illustrate the new approach: the calculation of the vibrational states of methyl and acetonitrile. The calculations for methyl employ the single-layer MCTDH approach, a general potential energy surface, and the correlation discrete variable representation. A five-layer MCTDH representation and a sum of product-type Hamiltonian are used in the acetonitrile calculations. Very fast convergence and order of magnitude reductions in the numerical effort compared to the previously used block relaxation scheme are found. Furthermore, a detailed comparison with the results of Avila and Carrington [J. Chem. Phys. 134, 054126 (2011)] for acetonitrile highlights the potential problems of convergence tests for high-dimensional systems.

A non-hierarchical correlation discrete variable representation

The correlation discrete variable representation (CDVR) facilitates (multi-layer) multi-configurational time-dependent Hartree (MCTDH) calculations with general potentials. It employs a layered grid representation to efficiently evaluate all potential matrix elements appearing in the MCTDH equations of motion. The original CDVR approach and its multi-layer extension show a hierarchical structure: the size of the grids employed at the different layers increases when moving from an upper layer to a lower one. In this work, a non-hierarchical CDVR approach, which uses identically structured quadratures at all layers of the MCTDH wavefunction representation, is introduced. The non-hierarchical CDVR approach crucially reduces the number of grid points required, compared to the hierarchical CDVR, shows superior scaling properties, and yields identical results for all three representations showing the same topology. Numerical tests studying the photodissociation of NOCl and the vibrational states of CH₃ demonstrate the accuracy of the non-hierarchical CDVR approach.

Vibrational control of the reaction pathway in the H + CHD3 → H2 + CD3 reaction

An accurate full-dimensional quantum state-to-state simulation of the six-atom title reaction based on first-principles theory is reported. Counterintuitive effects are found: Increasing the energy in the reactant's CD₃ umbrella vibration reduces the energy in the corresponding product vibration. An in-depth analysis reveals the crucial role of the effective dynamical transition state: Its geometry is controlled by the vibrational states of the reactants and subsequently controls the quantum state distribution of the products. This finding enables generalizing the concept of transition state control of chemical reactions to the quantum state-specific level.

Simulation of the photodetachment spectra of the nitrate anion (NO3-) in the B 2E' energy range and non-adiabatic electronic population dynamics of NO3

The photodetachment spectrum of the nitrate anion (NO3-) in the energy range of the NO3 second excited state is simulated from first principles using quantum wave packet dynamics. The prediction...

Accurate quantum dynamics simulation of the photodetachment spectrum of the nitrate anion (NO3 -) based on an artificial neural network diabatic potential model

The photodetachment spectrum of the nitrate anion (NO₃ ^-) is simulated from first principles using wavepacket quantum dynamics propagation and a newly developed accurate full-dimensional fully coupled five state diabatic potential model. This model utilizes the recently proposed complete nuclear permutation inversion invariant artificial neural network diabatization technique [D. M. G. Williams and W. Eisfeld, J. Phys. Chem. A 124, 7608 (2020)]. The quantum dynamics simulations are designed such that temperature effects and the impact of near threshold detachment are taken into account. Thus, the two available experiments at high temperature and at cryogenic temperature using the slow electron velocity-map imaging technique can be reproduced in very good agreement. These results clearly show the relevance of hot bands and vibronic coupling between the X̃ ²A₂ ^' ground state and the B̃ ²E' excited state of the neutral radical. This together with the recent experiment at low temperature gives further support for the proper assignment of the ν₃ fundamental, which has been debated for many years. An assignment of a not yet discussed hot band line is also proposed.

Full-dimensional quantum mechanical calculations of the reaction probability of the H + CH4 reaction based on a mixed Jacobi and Radau description

A full-dimensional time-dependent wave packet study using mixed polyspherical Jacobi and Radau coordinates for the title reaction has been reported. The non-reactive moiety CH₃ has been described using three Radau vectors, whereas two Jacobi vectors have been used for the bond breaking/formation process. A potential-optimized discrete variable representation basis has been employed to describe the vibrational coordinates of the reagent CH₄. About one hundred billion basis functions have been necessary to achieve converged results. The reaction probabilities for some initial vibrational states are given. A comparison between the present approach and other methods, including reduced and full-dimensional ones, is also presented.

Efficient anharmonic vibrational spectroscopy for large molecules using local-mode coordinates

This article presents a general computational approach for efficient simulations of anharmonic vibrational spectra in chemical systems. An automated local-mode vibrational approach is presented, which borrows techniques from localized molecular orbitals in electronic structure theory. This approach generates spatially localized vibrational modes, in contrast to the delocalization exhibited by canonical normal modes. The method is rigorously tested across a series of chemical systems, ranging from small molecules to large water clusters and a protonated dipeptide. It is interfaced with exact, grid-based approaches, as well as vibrational self-consistent field methods. Most significantly, this new set of reference coordinates exhibits a well-behaved spatial decay of mode couplings, which allows for a systematic, a priori truncation of mode couplings and increased computational efficiency. Convergence can typically be reached by including modes within only about 4 Å. The local nature of this truncation suggests particular promise for the ab initio simulation of anharmonic vibrational motion in large systems, where connection to experimental spectra is currently most challenging.

Resonances in the entrance channel of the elementary chemical reaction of fluorine and methane

Extending the fully quantum-state-resolved description of elementary chemical reactions beyond three or four atom systems is a crucial issue in fundamental chemical research. Reactions of methane with F, Cl, H or O are key examples that have been studied prominently. In particular, reactive resonances and nonintuitive mode-selective chemistry have been reported in experimental studies for the F+CH4 →HF+CH3 reaction. By investigating this reaction using transition-state spectroscopy, this joint theoretical and experimental study provides a clear picture of resonances in the F+CH4 system. This picture is deduced from high-resolution slow electron velocity-map imaging (SEVI) spectra and accurate full-dimensional (12D) quantum dynamics simulations in the picosecond regime.

Layered discrete variable representations and their application within the multiconfigurational time-dependent Hartree approach

The multiconfigurational time-dependent Hartree (MCTDH) approach facilitates multidimensional quantum dynamics calculations by employing a layered representation: in the upper layer of the representation the multidimensional wave function is expanded in a basis of time-dependent single-particle functions, in the lower layer then these single-particle functions are represented in a larger time-independent basis. This two layer approach can be generalized to a multilayer MCTDH representation to further increase the numerical efficiency. In recent work [U. Manthe, J. Chem. Phys. 128, 164116 (2008)], a quadrature scheme has been introduced which facilitates multilayer MCTDH calculations for Hamiltonians with general potential energy functions. Based on time-dependent discrete variable representation (DVR) grids corresponding to the single-particle function bases the scheme employs the correlation DVR (CDVR) approach. The present work shows that serious problems arise when this original multilayer CDVR approach is applied for accurate calculations of vibrational states. A solution to these problems is presented which uses a revised scheme to define the time-dependent DVR grids entering into the CDVR scheme. Applications studying the vibrational states of CH(3) demonstrate that the resulting revised multilayer MCTDH/CDVR approach yields smoothly convergent and accurate results.