Multidimensional time-dependent discrete variable representations in multiconfiguration Hartree calculations

In the multiconfiguration time-dependent Hartree (MCTDH) approach, the wave function is expanded in time-dependent basis functions, called single-particle functions, to increase the efficiency of the wave-packet propagation. The correlation discrete variable representation (CDVR) approach, which is based on a time-dependent discrete variable representation (DVR), can be employed to evaluate matrix elements of the potential energy. The efficiency of the MCTDH method can be further enhanced by using multidimensional single-particle functions. However, up to now the CDVR approach could not be used in MCTDH calculations employing multidimensional single-particle functions, since this would require a general multidimensional non-direct-product DVR scheme. Recently, Dawes and Carrington presented a practical scheme to implement general non-direct-product multidimensional DVRs [R. Dawes and T. Carrington, Jr., J. Chem. Phys. 121, 726 (2004)]. The present work utilizes their scheme in the MCTDH/CDVR approach. The accuracy is tested using the photodissociation of NOCl as example. The results show that the CDVR scheme based on multidimensional time-dependent DVRs allows for an accurate evaluation of the potential in MCTDH calculations with multidimensional single-particle functions.     

A multilayer multiconfigurational time-dependent Hartree approach for quantum dynamics on general potential energy surfaces

The multiconfigurational time-dependent Hartree (MCTDH) approach facilitates multidimensional quantum dynamics calculations by representing the wavepacket in an optimal set of time-dependent basis functions, called single-particle functions. Choosing these single-particle functions to be themselves multidimensional wavefunctions which are represented using a MCTDH representation, a multilayer MCTDH scheme has been constructed and used for quantum dynamics calculations treating up to 1000 degrees of freedom rigorously [Wang and Thoss, J. Chem. Phys. 199, 1289 (2003)]. The present work gives a practical scheme which facilitates the application of the multilayer MCTDH approach, which previously has only been employed to study systems described by model-type Hamiltonians, to molecular systems described by more complicated Hamiltonians and general potential energy surfaces. A multilayer extension of the correlation discrete variable representation (CDVR) scheme employed in MCTDH calculations studying quantum dynamics on general potential energy surfaces is developed and tested in a simple numerical application. The resulting multilayer MCTDH/CDVR approach might offer a perspective to rigorously describe the quantum dynamics of larger polyatomic systems.     

Rotational excitations in para-H2+para-H2 collisions: full- and reduced-dimensional quantum wave packet studies comparing different potential energy surfaces

We study the process of rotational excitation in the collisions of para-H(2) with para-H(2) by propagating wave packets with the multiconfiguration time-dependent Hartree (MCTDH) algorithm. Transition probabilities are then calculated by the method of Tannor and Weeks based on time-correlation functions. Calculations were carried out up to a total angular momentum of J=70 to compute integral cross sections up to 1.2 eV in collision energy and thermal rate coefficients from 100 to 3000 K. The process is studied on the full-dimensional potential energy surface of Boothroyd-Martin-Keogh-Peterson (BMKP) as well as on the rigid rotor surface of Diep and Johnson. We test the validity of the rigid rotor approximation by also considering two rigid rotor restrictions of the BMKP potential energy surface (PES). Additionally, we investigate a variant of the BMKP PES suggested by Pogrebnya and Clary [Chem. Phys. Lett. 363, 523 (2002)] with reduced anisotropy. We compare our results with previous theoretical data for the cross sections and with experimental data for the rate coefficients at low temperatures.     

Quantum dynamics of dissociative chemisorption of CH(4) on Ni(111): Influence of the bending vibration

Two-dimensional, three-dimensional, and four-dimensional quantum dynamic calculations are performed on the dissociative chemisorption of CH(4) on Ni(111) using the multiconfiguration time-dependent Hartree (MCTDH) method. The potential energy surface used for these calculations is 15-dimensional (15D) and was obtained with density functional theory for points which are concentrated in the region that is dynamically relevant to reaction. Many reduced dimensionality calculations were already performed on this system, but the molecule was generally treated as pseudodiatomic. The main improvement of our model is that we try to describe CH(4) as a polyatomic molecule by including a degree of freedom describing a bending vibration in our three-dimensional and four-dimensional models. Using a polyspherical coordinate system, a general expression for the 15D kinetic energy operator is derived, which discards all the singularities in the operator and includes rotational and Coriolis coupling. We use seven rigid constraints to fix the CH(3) umbrella of the molecule to its gas phase equilibrium geometry and to derive two-dimensional, three-dimensional, and four-dimensional Hamiltonians, which were used in the MCTDH method. Only four degrees of freedom evolve strongly along the 15D minimum energy path: the distance of the center of mass of the molecule to the surface, the dissociative C[Single Bond]H bond distance, the polar orientation of the molecule, and the bending angle between the dissociative C[Single Bond]H bond and the umbrella. A selection of these coordinates is included in each of our models. The polar rotation is found to be important in determining the mode selective behavior of the reaction. Furthermore, our calculations are in good agreement with the finding of Xiang et al. [J. Chem. Phys. 117, 7698 (2002)] in their reduced dimensional calculation that the helicopter motion of the umbrella symmetry axis is less efficient than its cartwheel motion for promoting the reaction. The effect of pre-exciting the bend modes is qualitatively incorrect at higher energies, suggesting the necessity of including additional rotational and vibrational degrees of freedom in the model.     

Complexes of type C6H7(+)·L (L = N2 and CO2) studied by explicitly correlated coupled cluster theory

Complexes of the benzenium ion (C(6)H(7)(+)) with N(2) or CO(2) have been studied by explicitly correlated coupled cluster theory at the CCSD(T)-F12x (x = a, b) level [T. B. Adler et al., J. Chem. Phys. 127, 221106 (2007)] and the double-hybrid density functional B2PLYP-D [T. Schwabe and S. Grimme, Phys. Chem. Chem. Phys. 9, 3397 (2007)]. Improved harmonic vibrational wavenumbers for C(6)H(7)(+) have been obtained by CCSD(T?)-F12a calculations with the VTZ-F12 basis set. Combining them with previous B2PLYP-D anharmonic contributions we arrive at anharmonic wavenumbers which are in excellent agreement with recent experimental data from p-H(2) matrix isolation IR spectroscopy [M. Bahou et al., J. Chem. Phys. 136, 154304 (2012)]. The energetically most favourable conformer of C(6)H(7)(+)·N(2) shows a π-bonded structure similar to C(6)H(7)(+)·Rg (Rg = Ne, Ar) [P. Botschwina and R. Oswald, J. Phys. Chem. A 115, 13664 (2011)] with D(e) ≈ 870 cm(-1). For C(6)H(7)(+)·CO(2), a slightly lower energy is calculated for a conformer with the CO(2) ligand lying in the ring-plane of the C(6)H(7)(+) moiety (D(e) ≈ 1508 cm(-1)). It may be discriminated from other conformers through a strong band predicted at 1218 cm(-1), red-shifted by 21 cm(-1) from the corresponding band of free C(6)H(7)(+).     

Full dimensional (15-dimensional) quantum-dynamical simulation of the protonated water dimer. II. Infrared spectrum and vibrational dynamics

The infrared absorption spectrum of the protonated water dimer (H5O2+) is simulated in full dimensionality (15 dimensional) in the spectral range of 0-4000 cm(-1). The calculations are performed using the multiconfiguration time-dependent Hartree (MCTDH) method for propagation of wavepackets. All the fundamentals and several overtones of the vibrational motion are computed. The spectrum of H5O2+ is shaped to a large extent by couplings of the proton-transfer motion to large amplitude fluxional motions of the water molecules, water bending and water-water stretch motions. These couplings are identified and discussed, and the corresponding spectral lines are assigned. The large couplings featured by H5O2+ do not hinder, however, to describe the coupled vibrational motion by well defined simple types of vibration (stretching, bending; etc.) based on well defined modes of vibration, in terms of which the spectral lines are assigned. Comparison of our results to recent experiments and calculations on the system is given. The reported MCTDH IR spectrum is in very good agreement to the recently measured spectrum by Hammer et al. [J. Chem. Phys. 122, 244301 (2005)].     

Quantum Monte Carlo calculations of the dissociation energy of the water dimer

We report diffusion quantum Monte Carlo (DMC) calculations of the equilibrium dissociation energy D(e) of the water dimer. The dissociation energy measured experimentally, D(0), can be estimated from D(e) by adding a correction for vibrational effects. Using the measured dissociation energy and the modern value of the vibrational energy Mas et al., [J. Chem. Phys. 113, 6687 (2000)] leads to D(e)=5.00+/-0.7 kcal mol(-1), although the result Curtiss et al., [J. Chem. Phys. 71, 2703 (1979)] D(e)=5.44+/-0.7 kcal mol(-1), which uses an earlier estimate of the vibrational energy, has been widely quoted. High-level coupled cluster calculations Klopper et al., [Phys. Chem. Chem. Phys. 2, 2227 (2000)] have yielded D(e)=5.02+/-0.05 kcal mol(-1). In an attempt to shed new light on this old problem, we have performed all-electron DMC calculations on the water monomer and dimer using Slater-Jastrow wave functions with both Hartree-Fock approximation (HF) and B3LYP density functional theory single-particle orbitals. We obtain equilibrium dissociation energies for the dimer of 5.02+/-0.18 kcal mol(-1) (HF orbitals) and 5.21+/-0.18 kcal mol(-1) (B3LYP orbitals), in good agreement with the coupled cluster results.     

Quantum-classical dynamics of scattering processes in adiabatic and diabatic representations

We demonstrate the workability of a TDDVR based [J. Chem. Phys. 118, 5302 (2003)], novel quantum-classical approach, for simulating scattering processes on a quasi-Jahn-Teller model [J. Chem. Phys. 105, 9141 (1996)] surface. The formulation introduces a set of DVR grid points defined by the Hermite part of the basis set in each dimension and allows the movement of grid points around the central trajectory. With enough trajectories (grid points), the method converges to the exact quantum formulation whereas with only one grid point, we recover the conventional molecular dynamics approach. The time-dependent Schrodinger equation and classical equations of motion are solved self-consistently and electronic transitions are allowed anywhere in the configuration space among any number of coupled states. Quantum-classical calculations are performed on diabatic surfaces (two and three) to reveal the effects of symmetry on inelastic and reactive state-to-state transition probabilities, along with calculations on an adiabatic surface with ordinary Born-Oppenheimer approximation. Excellent agreement between TDDVR and DVR results is obtained in both the representations.     

Full-dimensional quantum dynamics of vibrationally highly excited NHD2

We report on full-dimensional vibrational quantum dynamics of the highly excited ammonia isotopologue NHD(2) using a newly developed potential energy surface and the MCTDH program package. The calculations allow to realistically simulate an infrared laser induced stereomutation reaction at the pyramidal nitrogen atom in the femtosecond time domain. Our results allow for a thorough qualitative and quantitative understanding of infrared photoinduced stereomutation kinetics, the underlying quantum dynamics, and the reaction mechanisms. Comparison is made with a previous, reduced dimensionality study of the same reaction [R. Marquardt, M. Quack, I. Thanopulos, and D. Luckhaus, J. Chem. Phys. 118, 643 (2003)], and it is shown that slight variances of reduced spaces lead to significantly different kinetics. Because the quantum dynamics depends subtly on variances of reduced spaces, reduced dimensionality treatments are not reliable even for qualitative predictions of the stereomutation kinetics. The first direct comparison between the Multiconfigurational Time Dependent Hartree [M. H. Beck, A. Ja?ckle, G. A. Worth et al., Phys. Rep. 324, 1 (2000)] and Unimolecular Reactions Induced by Monochromatic Infrared Radiation [M. Quack and E. Sutcliffe, QCPE Bulletin 6, 98 (1986)] program packages on a specific, four dimensional quantum dynamical problem allows for their full validation in the present work.     

Universal state-selective corrections to multi-reference coupled-cluster theories with single and double excitations

The recently proposed universal state-selective (USS) corrections [K. Kowalski, J. Chem. Phys. 134, 194107 (2011)] to approximate multi-reference coupled-cluster (MRCC) energies can be commonly applied to any type of MRCC theory based on the Jeziorski-Monkhorst [B. Jeziorski and H. J. Monkhorst, Phys. Rev. A 24, 1668 (1981)] exponential ansatz. In this paper we report on the performance of a simple USS correction to the Brillouin-Wigner and Mukherjee's MRCC approaches employing single and double excitations (USS-BW-MRCCSD and USS-Mk-MRCCSD). It is shown that the USS-BW-MRCCSD correction, which employs the manifold of single and double excitations, can be related to a posteriori corrections utilized in routine BW-MRCCSD calculations. In several benchmark calculations we compare the USS-BW-MRCCSD and USS-Mk-MRCCSD results with the results obtained with the full configuration interaction method.     

State-to-state rotational transitions in H2+H2 collisions at low temperatures

We present quantum mechanical close-coupling calculations of collisions between two hydrogen molecules over a wide range of energies, extending from the ultracold limit to the superthermal region. The two most recently published potential energy surfaces for the H(2)-H(2) complex, the so-called Diep-Johnson (DJ) [J. Chem. Phys. 112, 4465 (2000); 113, 3480 (2000)] and Boothroyd-Martin-Keogh-Peterson (BMKP) [J. Chem. Phys. 116, 666 (2002)] surfaces, are quantitatively evaluated and compared through the investigation of rotational transitions in H(2)+H(2) collisions within rigid rotor approximation. The BMKP surface is expected to be an improvement, approaching chemical accuracy, over all conformations of the potential energy surface compared to previous calculations of H(2)-H(2) interaction. We found significant differences in rotational excitation/deexcitation cross sections computed on the two surfaces in collisions between two para-H(2) molecules. The discrepancy persists over a large range of energies from the ultracold regime to thermal energies and occurs for several low-lying initial rotational levels. Good agreement is found with experiment B. Mate et al., [J. Chem. Phys. 122, 064313 (2005)] for the lowest rotational excitation process, but only with the use of the DJ potential. Rate coefficients computed with the BMKP potential are an order of magnitude smaller.     

Effect of the overall rotation on the cis-trans isomerization of HONO induced by an external field

Rovibrational eigenenergies of HONO are computed and compared to experimental energies available in the literature. For their computation, we use a previously developed potential energy surface (PES) and a newly derived exact kinetic energy operator (KEO) including the overall rotation for a tetra-atomic molecule in non-orthogonal coordinates. In addition, we use the Heidelberg Multi-Configuration Time-Dependent Hartree (MCTDH) package. We compare the experimental rovibrational eigenvalues of HONO available in the literature with those obtained with MCTDH and a previously developed potential energy surface (PES) [F. Richter et al., J. Chem. Phys., 2004, 120, 1306.] for the cis geometry. The effect of the overall rotation on the process studied in our previous work on HONO [F. Richter et al., J. Chem. Phys., 2007, 127, 164315.] leading to the cis→trans isomerization of HONO is investigated. This effect on this process is found to be weak.     

Parallel algorithm for the computation of the Hartree-Fock exchange matrix: gas phase and periodic parallel ONX

In this paper we present an efficient parallelization of the ONX algorithm for linear computation of the Hartree-Fock exchange matrix [J. Chem. Phys. 106, 9708 (1997)]. The method used is based on the equal time (ET) partitioning recently introduced [J. Chem. Phys. 118, 9128 (2003)] and [J. Chem. Phys. 121, 6608 (2004)]. ET exploits the slow variation of the density matrix between self-consistent-field iterations to achieve load balance. The method is presented and some benchmark calculations are discussed for gas phase and periodic systems with up to 128 processors. The current parallel ONX code is able to deliver up to 77% overall efficiency for a cluster of 50 water molecules on 128 processors (2.56 processors per heavy atom) and up to 87% for a box of 64 water molecules (two processors per heavy atom) with periodic boundary conditions.     

Theoretical investigation of highly excited vibrational states in DFCO: calculation of the out-of-plane bending states and simulation of the intramolecular vibrational energy redistribution

A previously developed modified Davidson scheme [C. Iung and F. Ribeiro, J. Chem. Phys. 121, 174105 (2005)] is applied to compute and analyze highly excited (nu2,nu6) eigenstates in DFCO. The present paper is also devoted to the simulations of the intramolecular vibrational energy redistribution (IVR) initiated by an excitation of the out-of-plane bending vibration (nnu6, n=2,4,6, . . . ,18, and 20). The multiconfiguration time-dependent Hartree method is exploited to propagate the corresponding six-dimensional wave packets. A comprehensive comparison with experimental data as well as with previous simulations of IVR in HFCO [G. Pasin et al. J. Chem. Phys. 124, 194304 (2006)] is presented.     

A density-functional approach to polarizable models: a Kim-Gordon response density interaction potential for molecular simulations

A combined linear-response-frozen electron-density model has been implemented in a molecular-dynamics scheme derived from an extended Lagrangian formalism. This approach is based on a partition of the electronic charge distribution into a frozen region described by Kim-Gordon theory [J. Chem. Phys. 56, 3122 (1972); J. Chem. Phys. 60, 1842 (1974)] and a response contribution determined by the instantaneous ionic configuration of the system. The method is free from empirical pair potentials and the parametrization protocol involves only calculations on properly chosen subsystems. We apply this method to a series of alkali halides in different physical phases and are able to reproduce experimental structural and thermodynamic properties with an accuracy comparable to Kohn-Sham density-functional calculations.     

Electronic structure and spectrum of UO2 2+ and UO2Cl4 2-

A theoretical study is presented of the electronic spectra of the UO(2) (2+) and UO(2)Cl(4) (2-) ions, based on multiconfigurational perturbation theory (CASSCF/CASPT2), combined with a recently developed method to treat spin-orbit coupling [P.-A. Malmqvist et al., Chem. Phys. Lett. 357, 230 (2002); B. O. Roos and P.-A. Malmqvist, Phys. Chem. Chem. Phys. 6, 2919 (2004)]. The results are compared to the experimental spectroscopic data obtained for uranyl ions in Cs(2)UO(2)Cl(4) crystals from Denning [Struct. Bonding (Berlin) 79, 215 (1992)] and to previous theoretical calculations performed using a combined configuration-interaction spin-orbit treatment [Z. Zhang and R. M. Pitzer, J. Phys. Chem. A 103, 6880 (1999); S. Matsika and R. M. Pitzer, J. Phys. Chem. A. 105, 637 (2001)]. As opposed to the latter results, the calculations performed in this work point to a significant effect of the weakly bound equatorial chlorine ligands on the excitation energies.     

Theoretical investigation of intramolecular vibrational energy redistribution in highly excited HFCO

The present paper is devoted to the simulations of the intramolecular vibrational energy redistribution (IVR) in HFCO initiated by an excitation of the out-of-plane bending vibration [nnu(6)=2,4,6,...,18,20]. Using a full six-dimensional ab initio potential energy, the multiconfiguration time-dependent Hartree (MCTDH) method was exploited to propagate the corresponding six-dimensional wave packets. This study emphasizes the stability of highly excited states of the out-of-plane bending mode which exist even above the dissociation threshold. More strikingly, the structure of the IVR during the first step of the dynamics is very stable for initial excitations ranging from 2nu(6) to 20nu(6). This latter result is consistent with the analysis of the eigenstates obtained, up to 10nu(6), with the aid of the Davidson algorithm in a foregoing paper [Iung and Ribeiro, J. Chem. Phys. 121, 174105 (2005)]. The present study can be considered as complementary to this previous investigation. This paper also shows how MCTDH can be used to predict the dynamical behavior of a strongly excited system and to determine the energies of the corresponding highly excited states.     

Polarizable and flexible model for ethanol

A polarizable, flexible model for ethanol is obtained based on an extensive series of B3LYP/6-311++G(d,p) calculations and molecular dynamics simulations. The ethanol model includes electric-field dependence in both the atomic charges and the intramolecular degrees of freedom. Field-dependent intramolecular potentials have been attempted only once previously, for OH and HH stretches in water [P. Cicu et al., J. Chem. Phys. 112, 8267 (2000)]. The torsional potential involving the hydrogen-bonding hydrogen in ethanol is found to be particularly field sensitive. The methodology for developing field-dependent potentials can be readily generalized to other molecules and is discussed in detail. Molecular dynamics simulations of bulk ethanol are performed and the results are assessed based on comparisons with the self-diffusion coefficient [N. Karger et al., J. Chem. Phys. 93, 3437 (1990)], dielectric constant [J. T. Kindt and C. A. Schmuttenmaer, J. Phys. Chem. 100, 10373 (1996)], enthalpy of vaporization [R. C. Wilhoit and B. J. Zwolinski, J. Phys. Chem. Ref. Data, Suppl. 2, 2 (1973)], and experimental interatomic distributions [C. J. Benmore and Y. L. Loh, J. Chem. Phys. 112, 5877 (2000)]. The simultaneous variation of the atomic charges and the intramolecular potentials requires modified equations of motion and a multiple time step algorithm has been implemented to solve these equations. The article concludes with a discussion of the bulk structure and properties with an emphasis on the hydrogen bonding network.