The authors present scalar-relativistic energy-consistent Hartree-Fock pseudopotentials for the main-group elements. The pseudopotentials do not exhibit a singularity at the nucleus and are therefore suitable for quantum Monte Carlo (QMC) calculations. They demonstrate their transferability through extensive benchmark calculations of atomic excitation spectra as well as molecular properties. In particular, they compute the vibrational frequencies and binding energies of 26 first- and second-row diatomic molecules using post-Hartree-Fock methods, finding excellent agreement with the corresponding all-electron values. They also show their pseudopotentials give superior accuracy than other existing pseudopotentials constructed specifically for QMC. Finally, valence basis sets of different sizes (VnZ with n=D,T,Q,5 for first and second rows, and n=D,T for third to fifth rows) optimized for our pseudopotentials are also presented.

Second order Møller-Plesset perturbation theory at the complete basis set limit and diffusion quantum Monte Carlo are used to examine several low energy isomers of the water hexamer. Both approaches predict the so-called prism to be the lowest energy isomer, followed by cage, book, and cyclic isomers. The energies of the four isomers are very similar, all being within 10-15 meV/ H 2 O. These reference data are then used to evaluate the performance of several density-functional theory exchange-correlation ͑xc͒ functionals. A subset of the xc functionals tested for smaller water clusters ͓I. Santra et al., J. Chem. Phys. 127, 184104 ͑2007͔͒ has been considered. While certain functionals do a reasonable job at predicting the absolute dissociation energies of the various isomers ͑coming within 10-20 meV/ H 2 O͒, none predict the correct energetic ordering of the four isomers nor does any predict the correct low total energy isomer. All xc functionals tested either predict the book or cyclic isomers to have the largest dissociation energies. A many-body decomposition of the total interaction energies within the hexamers leads to the conclusion that the failure lies in the poor description of van der Waals ͑dispersion͒ forces in the xc functionals considered. It is shown that the addition of an empirical pairwise ͑attractive͒ C 6 R −6 correction to certain functionals allows for an improved energetic ordering of the hexamers. The relevance of these results to density-functional simulations of liquid water is also briefly discussed.

No abstract

We present a simple, robust, and highly efficient method for optimizing all parameters of many-body wave functions in quantum Monte Carlo calculations, applicable to continuum systems and lattice models. Based on a strong zero-variance principle, diagonalization of the Hamiltonian matrix in the space spanned by the wave function and its derivatives determines the optimal parameters. It systematically reduces the fixed-node error, as demonstrated by the calculation of the binding energy of the small but challenging C(2) molecule to the experimental accuracy of 0.02 eV.

We present a mixed time-dependent density-functional theory (TDDFT)/classical trajectory surface hopping (SH) study of the photochemical ring opening in oxirane. Previous preparatory work limited to the symmetric CC ring-opening pathways of oxirane concluded that the Tamm-Dancoff approximation (TDA) is important for improving the performance of TDDFT away from the equilibrium geometry. This observation is supported by the present TDDFT TDA/SH calculations which successfully confirm the main experimentally derived Gomer-Noyes mechanism for the photochemical CO ring opening of oxirane and, in addition, provide important state-specific information not easily accessible from experiments. In particular, we find that, while one of the lowest two excited states is photochemically relatively inert, excitation into the other excited state leads predominantly to rapid ring opening, cyclic-C(2)H(4)O-->(*)CH(2)CH(2)O(*). This is followed by hopping to the electronic ground state where hot (4000 K) dynamics leads to further reactions, namely, (*)CH(2)CH(2)O()-->CH(3)CHO-->(*)CH(3)+(*)CHO and CH(4)+CO. We note that, in the dynamics, we are not limited to following minimum energy pathways and several surface hops may actually be needed before products are finally reached. The performance of different functionals is then assessed by comparison of TDDFT and diffusion Monte Carlo potential energy curves along a typical TDDFT TDA/SH reaction path. Finally, although true (S(0),S(1)) conical intersections are expected to be absent in adiabatic TDDFT, we show that the TDDFT TDA is able to approximate a conical intersection in this system.

The simplest cyanine dye series [H2N(CH)nNH2](+) with n = 1, 3, 5, 7, and 9 appears to be a challenge for all theoretical excited-state methods since the experimental spectra are difficult to predict and the observed deviations cannot be easily explained with standard arguments. We compute here the lowest vertical excitation energies of these dyes using a variety of approaches, namely, complete active space second-order perturbation theory (CASPT2), quantum Monte Carlo methods (QMC), coupled cluster linear response up to third approximate order (CC3), and various flavors of time-dependent density functional theory (TDDFT), including the recently proposed perturbative correction scheme (B2PLYP). In our calculations, all parameters such as basis set, active space, and geometry dependence are carefully analyzed. We find that all wave function methods give reasonably close excitation energies, with CASPT2 yielding the lowest values, and that the B2PLYP scheme gives excitations in satisfactory agreement with CC3 and DMC, significantly improving on the generalized gradient and hybrid approximations. Finally, to resolve the remaining discrepancy between predicted excitation energies and experimental absorption spectra, we also investigate the effect of excited-state relaxation. Our results indicate that a direct comparison of the experimental absorption maxima and the theoretical vertical excitations is not possible due to the presence of nonvertical transitions. The apparent agreement of earlier CASPT2 calculations with experiments was an artifact of the choice of active space and the use of an older definition of the zero-order Hamiltonian.

We perform a thorough comparative investigation of the excitation energies of the anionic and neutral forms of the green fluorescent protein (GFP) chromophore in the gas phase using a variety of first-principle theoretical approaches commonly used to access excited state properties of photoactive molecules. These include time-dependent density functional theory (TDDFT), complete-active-space second-order perturbation theory (CASPT2), equation-of-motion coupled cluster (EOM-CC), and quantum Monte Carlo (QMC) methods. We find that all approaches give roughly the same vertical excitation for the anionic form, while TDDFT predicts an excitation for the neutral chromophore significantly lower than the highly correlated methods. Our findings support the picture emerging from the extrapolation of the Kamlet-Taft fit of absorption experimental data in solution and indicate that the protein gives rise to a considerable bathochromic shift with respect to vacuum. These results also open some questions on the interpretation of photodestruction spectroscopy experiments in the gas phase as well as on the accuracy of previous theoretical calculations in the more complex protein environment.

We present a systematic investigation of the structural relaxation in the excited state of model retinal chromophores in the gas phase using the complete-active-space self-consistent theory (CASSCF), multiconfigurational second-order perturbation theory (CASPT2), quantum Monte Carlo (QMC), and coupled cluster (CC) methods. In contrast to the CASSCF photoisomerization mechanism of bond inversion followed by torsion around formal double bonds, we find that the other approaches predict an initial skeletal relaxation which does not lead to bond inversion but to a rather flexible retinal chromophore with longer bonds and with the bondlength pattern of the ground state being partly preserved. The relaxation proceeds then preferentially via partial torsion around formal single bonds and does not reach a conical intersection region. Our findings are compatible with solution experiments which point to the existence of multiple minima and relaxation pathways, some of which are nonreactive, do not lead to photoproducts via conical intersection, and are dominant in solution. Our results also demonstrate the importance of a balanced description of dynamical and static correlation in the excited-state gradients and raise serious concerns on the common use of the CASSCF method to investigate structural properties of photoexcited retinal systems.