Trygve Helgaker

BASIS-SET CONVERGENCE OF CORRELATED CALCULATIONS ON WATER

The basis-set convergence of the electronic correlation energy in the water molecule is investigated at the second-order Mo/ller–Plesset level and at the coupled-cluster singles-and-doubles level with and without perturbative triples corrections applied. The basis-set limits of the correlation energy are established to within 2 mEh by means of (1) extrapolations from sequences of calculations using correlation-consistent basis sets and (2) from explicitly correlated calculations employing terms linear in the interelectronic distances rij. For the extrapolations to the basis-set limit of the correlation energies, fits of the form a+bX−3 (where X is two for double-zeta sets, three for triple-zeta sets, etc.) are found to be useful. CCSD(T) calculations involving as many as 492 atomic orbitals are reported.

BASIS-SET CONVERGENCE IN CORRELATED CALCULATIONS ON NE, N2, AND H2O

Valence and all-electron correlation energies of Ne, N2, and H2O at fixed experimental geometries are computed at the levels of second-order perturbation theory (MP2) and coupled cluster theory with singles and doubles excitations (CCSD), and singles and doubles excitations with a perturbative triples correction (CCSD(T)). Correlation-consistent polarized valence and core-valence basis sets up to sextuple zeta quality are employed. Guided by basis-set limits established by rij-dependent methods, a number of extrapolation schemes for use with the correlation-consistent basis sets are investigated. Among the schemes considered here, a linear least-squares procedure applied to the quintuple and sextuple zeta results yields the most accurate extrapolations.

The Dalton quantum chemistry program system

Dalton is a powerful general‐purpose program system for the study of molecular electronic structure at the Hartree–Fock, Kohn–Sham, multiconfigurational self‐consistent‐field, Møller–Plesset, configuration‐interaction, and coupled‐cluster levels of theory. Apart from the total energy, a wide variety of molecular properties may be calculated using these electronic‐structure models. Molecular gradients and Hessians are available for geometry optimizations, molecular dynamics, and vibrational studies, whereas magnetic resonance and optical activity can be studied in a gauge‐origin‐invariant manner. Frequency‐dependent molecular properties can be calculated using linear, quadratic, and cubic response theory. A large number of singlet and triplet perturbation operators are available for the study of one‐, two‐, and three‐photon processes. Environmental effects may be included using various dielectric‐medium and quantum‐mechanics/molecular‐mechanics models. Large molecules may be studied using linear‐scaling and massively parallel algorithms. Dalton is distributed at no cost from http://www.daltonprogram.org for a number of UNIX platforms.

Excitation energies from the coupled cluster singles and doubles linear response function (CCSDLR). Applications to Be, CH+, CO, and H2O

The linear response function for a coupled cluster singles and doubles wave function is used to calculate vertical electronic energies for the closed shell system Be, CH+, CO, and H2O. It is shown that excitations of single electron replacement character can be described accurately in such an approach. Improved convergence is obtained using a preconditioned form of the coupled cluster linear response matrix.

Basis-set convergence of the energy in molecular Hartree–Fock calculations

Abstract The basis-set convergence towards the numerical limit of the Hartree–Fock total energy and binding energy is investigated for the correlation-consistent cc-pVXZ basis sets. For both energies, solid improvements are obtained with each increment in X . The basis-set errors for the total energy (Δ E ) fit an exponential form better than a power form, and the total energy is better fitted than the binding energy. It is difficult to find generally reliable extrapolation schemes for the total energy. In most cases, the most successful scheme gives results extrapolated beyond a given X that are comparable to the cc-pV(X+1)Z results, but occasionally it fails dramatically for large X . Indeed, explicit calculation of the energy in a larger basis set, especially the cc-pV6Z set for which Δ E ⩽0.1 mE h , gives the most reliable estimate of the basis-set limit.

Recent advances in wave function-based methods of molecular-property calculations.

Recent Advances in Wave Function-Based Methods of Molecular-Property Calculations Trygve Helgaker,* Sonia Coriani, Poul Jørgensen, Kasper Kristensen, Jeppe Olsen, and Kenneth Ruud Centre for Theoretical and Computational Chemistry, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway Dipartimento di Scienze Chimiche e Farmaceutiche, Universit a degli Studi di Trieste, Via Giorgieri 1, I-34127 Trieste, Italy Lundbeck Center for Theoretical Chemistry, Department of Chemistry, Aarhus University, 8000 Aarhus C, Denmark Centre for Theoretical and Computational Chemistry, Department of Chemistry, University of Tromsø, N-9037 Tromsø, Norway

The CC3 model: An iterative coupled cluster approach including connected triples

An alternative derivation of many-body perturbation theory (MBPT) has been given, where a coupled cluster parametrization is used for the wave function and the method of undetermined Lagrange multipliers is applied to set up a variational coupled cluster energy expression. In this variational formulation, the nth-order amplitudes determine the energy to order 2n+1 and the nth-order multipliers determine the energy to order 2n+2. We have developed an iterative approximate coupled cluster singles, doubles, and triples model CC3, where the triples amplitudes are correct through second order and the singles amplitudes are treated without approximations due to the unique role of singles as approximate orbital relaxation parameters. The compact energy expressions obtained from the variational formulation exhibit in a simple way the relationship between CC3, CCSDT-1a [Lee et al., J. Chem. Phys. 81, 5906 (1984)] CCSDT-1b models [Urban et al., J. Chem. Phys. 83, 4041 (1985)], and the CCSD(T) model [Raghavachari et...

Density-functional theory of linear and nonlinear time-dependent molecular properties

We present density-functional theory for linear and nonlinear response functions using an explicit exponential parametrization of the density operator. The response functions are derived using two alternative variation principles, namely, the Ehrenfest principle and the quasienergy principle, giving different but numerically equivalent formulas. We present, for the first time, calculations of dynamical hyperpolarizabilities for hybrid functionals including exchange-correlation functionals at the general gradient-approximation level and fractional exact Hartree–Fock exchange. Sample calculations are presented of the first hyperpolarizability of the para-nitroaniline molecule and of a porphyrin derived push–pull molecule, showing good agreement with available experimental data.

Assessment of a Coulomb-attenuated exchange–correlation energy functional

The recently proposed CAM-B3LYP exchange-correlation energy functional, based on a partitioning of the r operator in the exchange interaction into long- and short-range components, is assessed for the determination of molecular thermochemistry, structures, and second order response properties. Rydberg and charge transfer excitation energies and static electronic polarisabilities are notably improved over the standard B3LYP functional; classical reaction barriers also improve. Ionisation potentials, bond lengths, NMR shielding constants and indirect spin-spin coupling constants are comparable with the two functionals. CAM-B3LYP atomisation energies and diatomic harmonic vibrational wavenumbers are less accurate than those of B3LYP. Future research directions are outlined.

Basis set convergence of the interaction energy of hydrogen-bonded complexes

The Hartree-Fock and correlation contributions to the interaction energy of the hydrogen-bonded complexes (HF)2, (HCl)2, H2OHF, HCNHF, and (H2O)2 are computed in conventional calculations employing the aug-cc-pVXZ series of basis sets at the levels of Hartree-Fock theory, second-order perturbation theory, and coupled-cluster theory with single and double excitations augmented by a perturbative triples correction. The basis set convergence of the interaction energy is examined by comparison with results obtained with an explicitly correlated wave function model. The counterpoise-corrected and uncorrected Hartree-Fock interaction energies both converge very unsystematically. The convergence of the uncorrected correlation contribution is also very unsystematic because the basis set superposition error and the error from the incomplete description of the electronic Coulomb cusp both are present. Once the former has been effectively removed by the counterpoise correction, the cusp dominates and the convergence...