Note on origin-shift invariants and phase constraints for momentum space wave functions

The concept of phase invariance on shifting the origin formulated for the structure factors of crystallography is transferred to momentum space quantum chemistry. The resulting mathematical constraints in the form of inequalities are given for tractable one-electron molecules. © 1996 John Wiley & Sons, Inc.     

Confined helium atom low-lying S states analyzed through correlated Hylleraas wave functions and the Kohn-Sham model

Calculation including the electron correlation effects is reported for the ground 1 1S and lowest triplet 1 3S state energies of the confined helium atom placed at the center of an impenetrable spherical box. While the adopted wave-functional treatment involves optimization of three nonlinear parameters and 10, 20, and 40 linear coefficients contained in wave functions expressed in a generalized Hylleraas basis set that explicitly incorporates the interelectronic distance r12, via a Slater-type exponent and through polynomial terms entering the expansion, the Kohn-Sham model employed here uses the Perdew and Wang exchange-correlation functional in its spin-polarized version within the local-density approximation (LDA) with and without the self-interaction correction. All these calculations predict a systematic increase in the singlet-triplet energy splitting toward the high confinement regime, i.e., when the box radius is reduced. By using the variational results as benchmark, it is found that the LDA underestimates the singlet-triplet energy splitting, whereas the self-interaction correction overestimates such a quantity.     

Calculations of D‐wave bound states and resonance states of the screened helium atom using correlated exponential wave functions

We have investigated the effects of screened Coulomb (Yukawa) potentials on the bound ^1,3D states and the doubly excited ^1,3 D^e resonance states of helium atom using highly correlated exponential basis functions. The Density of resonance states are calculated using stabilization method. Highly correlated exponential basis functions are used to consider the correlation effect between the charged particles. A total of 18 resonances (nine each for ¹ D^e and ³ D^e states) below the n = 2 He ⁺ threshold has been calculated. For each spin states, this includes four members in the 2pnp series, three members in the 2snd series, and two members in 2pnf series. The resonance energies and widths for various screening parameters ranging from infinity to a small value for these ^1,3 D^e resonance states are reported along with the bound‐excited 1s3d ^1,3 D state energies. Overall behavior of the spectral profile of 1s3d ¹D state of helium atom due to electron‐electron and electron‐nucleus screening are also presented. Accurate resonance energies and widths are also reported for He in vacuum. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Analytic functions for the three-body potential of the helium trimer

The three-body potential for the ground state of the helium trimer is determined by an extended geminal model. The basis set for the calculation is an uncontracted (19s,7p,6d,5f,4g,2h) set of Gaussian-type functions. Three different types of configurations were considered: (i) equilateral triangles, (ii) linear configurations with R12=R23, and (iii) a set of pseudorandom configurations. The interatomic distances were selected within the interval [3.0,9.0] bohrs. The computed points have been fitted to global potential functions. The fit is characterized by a maximum absolute error equal to 0.69 microEh and a mean error equal to -0.018 microEh.     

Explicitly correlated wave functions: summary and perspective

We summarize explicitly correlated electronic structure theory in perspective of future work in the field. Earlier stages of approaches with different Ansätze in physics and chemistry are described. We then discuss recent advances focusing on explicitly correlated wave functions using cusp conditions. Removal of Coulomb singularities in terms of the rational generator is brought out from the viewpoint of many-body perturbation theory. On the basis of decomposition schemes for many-electron integrals in R12 and F12 methods, we further discuss the possibility of increasing the accuracy of molecular numerical integration and massively parallel calculations of explicitly correlated methods.     

Optimization of augmentation functions for correlated calculations of spin-spin coupling constants and related properties

A new hierarchy of augmented basis sets optimized for the calculation of molecular properties such as indirect spin-spin coupling constants is presented. Based on the Dunning hierarchy of cc-pVXZ (X = D, T, Q, and 5) basis sets augmentation functions with tight exponents have been optimized for coupled-cluster calculations of indirect spin-spin coupling constants. The optimal exponents for these tight functions have been obtained by optimizing the sum of the absolute values of all contributions to the coupling constant. On the basis of a series of test cases (CO, HF, N(2), F(2), H(2)O, NH(3), and CH(4)) we propose a set of tight s, p, and d functions to be added to the uncontracted Dunning basis sets, and, subsequently, to recontract. The resulting ccJ-pVXZ (X = D, T, Q, and 5) basis sets demonstrate excellent cost efficiency in benchmark calculations. These new basis sets should generally be applicable for the calculation of spin-spin coupling constants and other properties that have a strong dependence on powers of 1r or even contain a delta distribution for correlated ab initio methods.     

Three stages of optimization and simple correlated wave functions

It is advocated to carry out an optimization procedure, which is based upon the variational method, in such a way that the optimum values of the variational parameters are expressed as functions of physical constants, such as the atomic number, Z. The three stages involved in this treatment are illustrated by the optimization of nine correlated wave functions, which describe the ground states of atomic two-electron systems. An analysis of the Z-expansions of the total energies associated with these functions leads to the concept of a class of variational functions. The performances of functions belonging to the same class differ only marginally, especially at larger values of Z. Consequently, the concept of class may be used to bring some order in the plethora of variational functions. © 1994 John Wiley & Sons, Inc.     

Product of group-function expansions for correlated wave functions

A theorem is proved which demonstrates the relationship between a product of group functions describing the correlated motion of a particular group of electrons in an N-electron system and a wave function obtained from the exact wave function which describes the correlation of the same group of electrons. By considering such products of group functions as elements in a variational wave function, an expansion for correlated wave functions is suggested, which emphasizes the correlated motion of groups of electrons in the whole system.     

Diatomic interactions in momentum space. The effect of polarization and floating functions

The method of momentum density for interatomic interactions is used to investigate the pictures and roles of the polarization and floating functions in momentum (p-) space. Referring to the previous results from the minimal LCAO (Finkelstein-Horowitz) momentum density, we quantitatively discuss the effect of these functions for the bonding process in the ground state of H ₂ ⁺ system. The essence of the polarization and floating effects is found to be a modulation of the oscillation in the two-center part of the momentum density. The polarization function introduces a term with a phase and the floating function enlarges the period of the oscillation. An increased migration of the density from the one-center to the two-center part is also important. As a result, both the functions contribute to emphasize the contraction and expansion of momentum density observed previously. However, the floating function disturbs the density distribution in high momentum region, reflecting the destruction of cusps in position (r-) space. We point out an error in the pioneer work of Duncanson.     

Correlated electron-pair properties of the Li atom in momentum space

On the basis of multiconfiguration Hartree–Fock calculations, correlated electron-pair intracule (relative motion) and extracule (center-of-mass motion) properties are reported for the Li atom in momentum space. The present results are more accurate and consistent than those in the literature. Received: 10 September 2001 / Accepted: 11 December 2001 / Published online: 22 March 2002     

Localized molecular orbital studies in momentum space. II. Correspondence between coordinate and momentum space properties of various electron pairs

The momentum space properties of the ten-electron systems Ne, HF, H₂O, NH₃ and CH₄ as well as those of CH₃CH₃, CH₃NH₂, CH₃OH and FCH₂OH were investigated using localized molecular orbitals (LMO) obtained from ab initio self-consistent-field (SCF) wavefunctions constructed from double zeta quality gaussain basis sets.Compton profiles of various LMO electron pairs (CC, CN, CO, CF; CH, NH, OH, FH bond pairs and C, N, O, F lone pairs) are tabulated. In order to understand the correspondence between the momentum and the coordinate space properties of those electron pairs, the concept of the size and the shape of an LMO electron pair charge distribution has been utilized. The use of the intermediate expectation values of pⁿ is introduced for the purpose of interpreting the momentum space properties.The dependence of molecular property partitioning on different localization schemes and on different basis sets is also studied by using the H₂O profile as an example.     

On calculations of correlated wave functions with heavy configurational mixing

We discuss aspects of the theory and computation of wave functions and energies of discrete states of polyelectronic atoms that are represented in zero order by configurations with holes in subshells below the valence subshell. Both in zero order and in the remaining correlation components, such wave functions have particularities stemming from the state‐specific self‐consistent field and the heavy configurational mixing associated with near‐degeneracies and hole‐filling correlations. By referring to a variety of examples from small‐ and large‐scale calculations, it is noted that appropriate penetration into the many‐body problem can provide, in an economic and physically transparent way, reliable interpretations and semi‐ and fully quantitative understanding of issues related to states with inner holes and to cases of near‐degeneracies that result in strongly correlated wave functions. Whenever hole‐filling correlations are allowed, multiple correlations (i.e., beyond single‐ and double‐orbital substitutions in the single reference configuration) acquire increased importance relative to that in ordinary electronic structures. This is demonstrated via large‐scale multiconfigurational Hartree–Fock (MCHF) plus configuration interaction (CI) calculations on the Cl KL3s3p^{6 2}S discrete state, which is the lowest of its symmetry. The calculations incorporated correlations up to selected sextuple orbital excitations from the M shell. MCHF plus CI calculations at the level of quadruple orbital substitutions were also carried out for the Cl KL3s²3p^{5 2}P^o ground state and the excitation energy at this level of calculation was found to be 85,364 cm^?1, in excellent agreement with the experimental value of the fine‐structure‐weighted average, 85,385 cm^?1 (10.59 eV). Within the approximations of the calculation, the hole‐filling triple and quadruple orbital correlations, which, of course, are absent from the ²P^o state, contribute about 1 eV, which is significant. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Analytic structure of ground-state energies and wave functions for the inhomogeneous electron liquid in non-relativistic He-like atomic ions with nuclear charge Ze

The ground state of the He atom for fixed nucleus remains intractable so far as regards exact analytic solutions. However, some important results already exist pertaining to its ground-state wave function Ψ and corresponding electron density n(r). Here, we extend the existing studies by focussing attention on the non-relativistic series of He-like atomic ions with nuclear charge Z. We then find it instructive to start from the energy E(Z) of such a two-electron spin-compensated problem. This is known to have non-analytic behaviour at a critical Z, say Z _c , equal to 0.911028. A form of Darboux transformation going back at very least to Brändas and Goscinski [E. Brändas and O. Goscinski, Int. J. Quantum Chem. 6, 59 (1972)] is refined somewhat here, and compared with a more intuitive approach of Callan [E. Callan, Int. J. Quantum Chem. 6, 431 (1972)]. The important 1/Z expansion of E(Z) is also invoked. The electron density n(r) and the ground-state wave function Ψ are then treated in turn, in a related manner; especially their asymptotic behaviour far from the nucleus. Finally, two exact wave functions for analytically solvable two-Fermion models are shown to sum the infinite series proposed by Fock [V. Fock, Izv. Akad. Nauk SSSR, Ser. Fiz. 18, 161 (1954)].     

Accurate explicitly correlated wave functions for two electrons in a square

An explicitly correlated linear-r(12) variational method is developed for a system of two electrons confined to a two-dimensional square well with infinite walls. The wave function is written as an expansion in products of non-negative integer powers of the relative and center-of-mass electronic coordinates and powers of r(12) restricted to 0 and 1. This form indirectly includes higher powers of the interelectronic distance and exhibits a much faster convergence than a similar expansion without r(12)-dependent terms. The method is implemented using high-precision floating-point arithmetic. Ground-state total energies are reported with at least 12 accurate significant figures for squares with sides from 1 to 50 bohrs. The method can be used "as is" for excited states and for two-dimensional rectangular wells. We also show that wave functions for two electrons in a square and in a rectangle have a higher symmetry than can be accounted for by the point group of the system.     

Excited states of boron isoelectronic series from explicitly correlated wave functions

The ground state and some low-lying excited states arising from the 1s2 2s2p2 configuration of the boron isoelectronic series are studied starting from explicitly correlated multideterminant wave functions. One- and two-body densities in position space have been calculated and different expectation values such as , , , , , and , where r, r12, and R stand for the electron-nucleus, interelectronic, and two electron center of mass coordinates, respectively, have been obtained. The energetic ordering of the excited states and the fulfillment of the Hund's rules is analyzed systematically along the isoelectronic series in terms of the electron-electron and electron-nucleus potential energies. The effects of electronic correlations have been systematically studied by comparing the correlated results with the corresponding noncorrelated ones. All the calculations have been done by using the variational Monte Carlo method.     

Accurate calculations of intermolecular interaction energies using explicitly correlated wave functions

Explicitly correlated second-order M?ller-Plesset (MP2-F12) calculations of intermolecular interaction energies for the S22 benchmark set of Jurecka, Sponer, Cerny, and Hobza (Chem. Phys. Phys. Chem. 2006, 8, 1985) are presented and compared with standard MP2 results. The MP2 complete basis set limits are estimated using basis set extrapolation and augmented quadruple-zeta and quintuple-zeta basis sets. Already with augmented double-zeta basis sets the MP2-F12 interaction energies are found to be closer to the complete basis set limits than standard MP2 calculations with augmented quintuple-zeta basis sets. Various possible approximations in the MP2-F12 method are systematically tested. Best results are obtained with localized orbitals and the diagonal MP2-F12/C(D) ansatz. Hybrid approximations, in which some contributions of the auxiliary basis set are neglected and which considerably reduce the computational cost, have a negligible effect on the interaction energies. Also the orbital-invariant fixed-amplitude approximation of Ten-no leads to only slightly less accurate results. Preliminary results for the neon and benzene dimers, obtained with the recently proposed CCSD(T)-F12a approximation, indicate that the CCSD(T) basis set limits can also be very closely approached using augmented triple-zeta basis sets.     

Analytic evaluation of infrared intensities and polarizabilities by two-configuration self-consistent field wave functions

For more than a few molecular electronic states, the simplest qualitatively correct picture of the electronic structure is provided by the two-configuration self-consistent-field (TCSCF) method. Here, analytic methods are reported for the evaluation of TCSCF infrared intensities and polarizabilities. These new methods have been implemented and applied to the molecules CH₂(¹A₁), CF₂, CH ₃ ^– , NH₃, HF and O₃. Nine different basis sets, ranging from 3–21G to triple zeta plus double polarization (TZ + 2P), have been used. In several cases one finds qualitative differences between the analogous SCF and TCSCF predictions.     

Evaluation of the second-order reduced density matrix for correlated electronic wave functions

Formulas are presented for the evaluation of the second-order reduced density matrix for correlated wave functions, including wave functions constructed from pairwise nonorthogonal orbitals. Numerical results are provided for the correlation holes and conditional nuclear spin densities for the water molecule.     

Two-electron integrations in the quantum theory of atoms in molecules with correlated wave functions

A recent method proposed to compute two-electron integrals over arbitrary regions of space [Martin Pendas, A. et al., J Chem Phys 2004, 120, 4581] is extended to deal with correlated wave functions. To that end, we use a monadic factorization of the second-order reduced density matrix originally proposed by E. R. Davidson [Chem Phys Lett 1995, 246, 209] that achieves a full separation of the interelectronic components into one-electron terms. The final computational effort is equivalent to that found in the integration of a one determinant wave function with as many orbitals as occupied functions in the correlated expansion. Similar strategies to extract the exchange and self-interaction contributions from the two-electron repulsion are also discussed, and several numerical results obtained in a few test systems are summarized.