Vibrationally correlated calculations in polyspherical coordinates: Taylor expansion-based kinetic energy operators.

The efficiency of quantum chemical simulations of nuclear motion can in many cases greatly benefit from the application of curvilinear coordinate systems. This is rooted in the fact that a set of smartly selected curvilinear coordinates may represent the motion naturally well, thus decreasing the couplings between motions in these coordinates. In this study, we assess the validity of different Taylor expansion-based approximations of kinetic energy operators in a (curvilinear) polyspherical parametrization. To this end, we investigate the accuracy as well as the numerical performance of the approximations in time-independent vibrational coupled cluster and full vibrational interaction calculations for several test cases ranging from tri- to penta-atomic molecules. We find that several of the proposed schemes reproduce the vibrational ground state and excitation energies to a decent accuracy, justifying their application in future investigations. Furthermore, due to the restricted mode coupling and their inherent sum-of-products form, the new approximations open up the possibility of treating large molecular systems with efficient vibrational coupled cluster schemes in general coordinates.

Anharmonic Vibrational Treatment Exclusively in Curvilinear Valence Coordinates: The Case of Formamide.

A highly correlated approach using curvilinear valence coordinates is applied to calculate the vibrational fundamentals and some combination modes of the formamide molecule with high accuracy. A series of potential energy surfaces (PESs) has been generated by AGAPES, a program for adaptive generation of adiabatic PESs, at various electronic structure qualities until excellent nonaccidental agreement with the experimentally assigned fundamental transitions was reached at the CCSDT(T)-F12a/aug-cc-pVTZ level of theory using the improved relaxation method of the Heidelberg multiconfiguration time-dependent Hartree (MCTDH) package in connection with an exact expression for the kinetic energy in valence coordinates generated by the TANA program. By comparison of the overtone series ν1-3ν1 to experiment, we demonstrate that the known problems concerning the floppy ν1 wagging motion are solved within this approach. The potential energy coupling as well as the vibrational coupling in curvilinear coordinates is discussed together with the efficiency of this approach.

Explicitly correlated treatment of H2NSi and H2SiN radicals: electronic structure calculations and rovibrational spectra.

Various ab initio methods are used to compute the six dimensional potential energy surfaces (6D-PESs) of the ground states of the H(2)NSi and H(2)SiN radicals. They include standard coupled cluster (RCCSD(T)) techniques and the newly developed explicitly correlated RCCSD(T)-F12 methods. For H(2)NSi, the explicitly correlated techniques are viewed to provide data as accurate as the standard coupled cluster techniques, whereas small differences are noticed for H(2)SiN. These PESs are found to be very flat along the out-of-plane and some in-plane bending coordinates. Then, the analytic representations of these PESs are used to solve the nuclear motions by standard perturbation theory and variational calculations. For both isomers, a set of accurate spectroscopic parameters and the vibrational spectrum up to 4000 cm(-1) are predicted. In particular, the analysis of our results shows the occurrence of anharmonic resonances for H(2)SiN even at low energies.

Theoretical spectroscopy of trans-HNNH(+) and isotopomers.

The six-dimensional potential energy surface of the electronic ground state of trans-HNNH(+) (X (2)A(g)) is mapped at the RCCSD(T)/aug-cc-pV5Z level of theory. This potential energy surface is incorporated later into perturbative and variational treatments to solve the nuclear motion and to derive a set of spectroscopic data for trans-HNNH(+), trans-HNND(+), and trans-DNND(+). Our vibrational spectra are compared with those deduced from the earlier photoelectron spectra by Frost et al. [J. Chem. Phys. 64, 4719 (1976)], for which a good agreement between the theoretical and experimental results is found. Our calculations reveal the presence of strong anharmonic resonances between the vibrational levels of these cations even at low energies, thus complicating even more their assignment by vibrational quantum numbers. These resonances should participate in the transfer of intensities between the active modes during the direct photoionization of the neutral molecule and the combination modes and overtones of the inactive modes belonging to the totally symmetric irreducible representation.

Vibrational computing: simulation of a full adder by optimal control.

Within the context of vibrational molecular quantum computing, we investigate the implementation of a full addition of two binary digits and a carry that provides the sum and the carry out. Four qubits are necessary and they are encoded into four different normal vibrational modes of a molecule. We choose the bromoacetyl chloride molecule because it possesses four bright infrared active modes. The ground and first excited states of each mode form the one-qubit computational basis set. Two approaches are proposed for the realization of the full addition. In the first one, we optimize a pulse that implements directly the entire addition by a single unitary transformation. In the second one, we decompose the full addition in elementary quantum gates, following a scheme proposed by Vedral et al. [Phys. Rev. A 54, 147 (1996)]. Four elementary quantum gates are necessary, two two-qubit CNOT gates (controlled NOT) and two three-qubit TOFFOLI gates (controlled-controlled NOT). All the logic operations consist in one-qubit flip. The logic implementation is therefore quasiclassical and the readout is based on a population analysis of the vibrational modes that does not take the phases into account. The fields are optimized by the multitarget extension of the optimal control theory involving all the transformations among the 2(4) qubit states. A single cycle of addition without considering the preparation or the measure or copy of the result can be carried out in a very competitive time, on a picosecond time scale.

Optimal control simulation of the Deutsch-Jozsa algorithm in a two-dimensional double well coupled to an environment.

The quantum Deutsch-Jozsa algorithm is implemented by using vibrational modes of a two-dimensional double well. The laser fields realizing the different gates (NOT, CNOT, and HADAMARD) on the two-qubit space are computed by the multitarget optimal control theory. The stability of the performance index is checked by coupling the system to an environment. Firstly, the two-dimensional subspace is coupled to a small number Nb of oscillators in order to simulate intramolecular vibrational energy redistribution. The complete (2+Nb)D problem is solved by the coupled harmonic adiabatic channel method which allows including coupled modes up to Nb=5. Secondly, the computational subspace is coupled to a continuous bath of oscillators in order to simulate a confined environment expected to be favorable to achieve molecular computing, for instance, molecules confined in matrices or in a fullerene. The spectral density of the bath is approximated by an Ohmic law with a cutoff for some hundreds of cm(-1). The time scale of the bath dynamics (of the order of 10 fs) is then smaller than the relaxation time and the controlled dynamics (2 ps) so that Markovian dissipative dynamics is used.