Communication: multireference equation of motion coupled cluster: a transform and diagonalize approach to electronic structure.

The novel multireference equation-of-motion coupled-cluster (MREOM-CC) approaches provide versatile and accurate access to a large number of electronic states. The methods proceed by a sequence of many-body similarity transformations and a subsequent diagonalization of the transformed Hamiltonian over a compact subspace. The transformed Hamiltonian is a connected entity and preserves spin- and spatial symmetry properties of the original Hamiltonian, but is no longer Hermitean. The final diagonalization spaces are defined in terms of a complete active space (CAS) and limited excitations (1h, 1p, 2h, …) out of the CAS. The methods are invariant to rotations of orbitals within their respective subspaces (inactive, active, external). Applications to first row transition metal atoms (Cr, Mn, and Fe) are presented yielding results for up to 524 electronic states (for Cr) with an rms error compared to experiment of about 0.05 eV. The accuracy of the MREOM family of methods is closely related to its favorable extensivity properties as illustrated by calculations on the O2-O2 dimer. The computational costs of the transformation steps in MREOM are comparable to those of closed-shell Coupled Cluster Singles and Doubles (CCSD) approach.

External coupled-cluster perturbation theory: description and application to weakly interaction dimers. Corrections to the random phase approximation.

The formalism for developing perturbation theory by using an arbitrary fixed (external) set of amplitudes as an initial approximation is presented in a compact form: external coupled-cluster perturbation theory (xCCPT). Nonperturbative approaches also fit into the formalism. As an illustration, the weakly interacting dimers Ne(2) and Ar(2) have been studied in the various ring-coupled-cluster doubles (CCD) approximations; ring, direct-ring, antisymmetrized ring, and antisymmetrized direct ring, and a second-order correction in the xCCPT approach is added. The direct approaches include the summation of just Coulomb terms with the intention of selectively summing the largest terms in the perturbation first. "Coulomb attenuation" is effected by taking the random phase approximation to define such amplitudes, whose results are then improved upon using perturbation theory. Interaction energies at the ring-CCD level are poor but the xCCPT correction employed predicts binding energies which are only a few percent from the coupled-cluster single double (triple) values for the direct ring-CCD variants. Using the MP2 amplitudes which neglect exchange, the initial Coulomb-only term, leads to very accurate Ne(2) and Ar(2) potentials. However, to accurately compute the Na(2) potential required a different initial wavefunction, and hence perturbation. The potential energy surfaces of Ne(2) and Ar(2) are much too shallow using linear coupled-cluster doubles. Using xCCPT(2) with these amplitudes as the initial wavefunction led to slightly worse results. These observations suggest that an optimal external set of amplitudes exists which minimizes perturbational effects and hence improve the predictability of methods.

Parallel implementation of electronic structure energy, gradient, and Hessian calculations.

ACES III is a newly written program in which the computationally demanding components of the computational chemistry code ACES II [J. F. Stanton et al., Int. J. Quantum Chem. 526, 879 (1992); [ACES II program system, University of Florida, 1994] have been redesigned and implemented in parallel. The high-level algorithms include Hartree-Fock (HF) self-consistent field (SCF), second-order many-body perturbation theory [MBPT(2)] energy, gradient, and Hessian, and coupled cluster singles, doubles, and perturbative triples [CCSD(T)] energy and gradient. For SCF, MBPT(2), and CCSD(T), both restricted HF and unrestricted HF reference wave functions are available. For MBPT(2) gradients and Hessians, a restricted open-shell HF reference is also supported. The methods are programed in a special language designed for the parallelization project. The language is called super instruction assembly language (SIAL). The design uses an extreme form of object-oriented programing. All compute intensive operations, such as tensor contractions and diagonalizations, all communication operations, and all input-output operations are handled by a parallel program written in C and FORTRAN 77. This parallel program, called the super instruction processor (SIP), interprets and executes the SIAL program. By separating the algorithmic complexity (in SIAL) from the complexities of execution on computer hardware (in SIP), a software system is created that allows for very effective optimization and tuning on different hardware architectures with quite manageable effort.

Efficient electronic integrals and their generalized derivatives for object oriented implementations of electronic structure calculations.

For the new parallel implementation of electronic structure methods in ACES III (Lotrich et al., in preparation) the present state-of-the-art algorithms for the evaluation of electronic integrals and their generalized derivatives were implemented in new object oriented codes with attention paid to efficient execution on modern processors with a deep hierarchy of data storage including multiple caches and memory banks. Particular attention has been paid to define proper integral blocks as basic building objects. These objects are stand-alone units and are no longer tied to any specific software. They can hence be used by any quantum chemistry code without modification. The integral blocks can be called at any time and in any sequence during the execution of an electronic structure program. Evaluation efficiency of these integral objects has been carefully tested and it compares well with other fast integral programs in the community. Correctness of the objects has been demonstrated by several application runs on real systems using the ACES III program.

A measure of response to perturbation used to assess structural change in some polluted and unpolluted stream fish communities.

A new method for measuring structural change in sets of species which have been subjected to natural or experimental perturbation is developed and is shown to be superior to static diversity and evenness measures for this purpose. Three parameters, HΔ', JΔ', and XΔ;} are shown to provide necessary and sufficient information on the severity of a perturbation as well as the uniformity of its effect on all species in the set. When positive and negative changes in species abundance are considered separately, the method is sensitive to compensatory changes which are not detected by static measures.The parameters are then calculated for some data sets on polluted and unpolluted fish communities in second and third order streams from the Clemons Fork watershed in eastern Kentucky. Results indicate that HΔ', the diversity of change over two sampling seasons, is high for perturbed and unperturbed systems, but JΔ' the eveness of change is lower for the communities which were polluted in the second sampling season. Severe pollution results in the suppression of most major fish species, whereas more moderate pollution results in a large number of compensatory changes. The biological basis for such an outcome is discussed, and the notion of these three parameters as the "vital signs" of a healthy ecosystem is presented.