It is shown that using an appropriate localized molecular orbital (LMO) basis, one is able to calculate coupled-cluster singles and doubles (CCSD) wave functions and energies for very large systems by performing full CCSD calculations on small subunits only. This leads to a natural linear scaling coupled-cluster method (NLSCC), in which total correlation energies of extended systems are evaluated as the sum of correlation energy contributions from individual small subunits within that system. This is achieved by defining local occupied orbital correlation energies. These are quantities, which in the LMO basis become transferable between similar molecular fragments. Conventional small scale existing molecular CCSD codes are all that is needed, the local correlation effect being simply transmitted via the appropriate LMO basis. Linear scaling of electronic correlation energy calculations is thus naturally achieved using the NLSCC approach, which in principle can treat nonperiodic extended systems of infinite basis set size. Results are shown for alkanes and several polyglycine molecules and the latter compared to recent results obtained via an explicit large scale LCCSD calculation. (c) 2004 American Institute of Physics.

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.

The universe is permeated by magnetic fields, with strengths ranging from a femtogauss in the voids between the filaments of galaxy clusters to several teragauss in black holes and neutron stars. The standard model behind cosmological magnetic fields is the nonlinear amplification of seed fields via turbulent dynamo to the values observed. We have conceived experiments that aim to demonstrate and study the turbulent dynamo mechanism in the laboratory. Here, we describe the design of these experiments through simulation campaigns using FLASH, a highly capable radiation magnetohydrodynamics code that we have developed, and large-scale three-dimensional simulations on the Mira supercomputer at the Argonne National Laboratory. The simulation results indicate that the experimental platform may be capable of reaching a turbulent plasma state and determining the dynamo amplification. We validate and compare our numerical results with a small subset of experimental data using synthetic diagnostics

The natural linear-scaled coupled-cluster (NLSCC) method ( Flocke, N.; Bartlett, R. J. J. Chem. Phys. 2004, 121, 10935 ) is extended to include approximate triple excitations via a coupled-cluster with single, double, and triple excitation method (CCSDT-3). The triples contribution can potentially be embedded in a larger singles and doubles region. NLSCC exploits the extensivity of the CC wave function to represent it in terms of transferable natural localized molecular orbitals (NLMOs) or functional groups thereof that are obtained from small quantum mechanical (QM) regions. Both occupied and virtual NLMOs are local because they derive from the single-particle density matrix. Noncanonical triples amplitudes are avoided by applying the unitary localization matrix to the canonical CC wave function for a QM region. A generalized NLMO code interfaced to the ACES II quantum chemistry software package provides NLMOs for the relevant number of atoms in a given functional group. Applications include linear polyglycine and the pentapeptide met-enkephalin, which was chosen as a more realistic three-dimensional system with nontrivial side chains. The results show that the triples contributions are quite large for aromatic bonds suggesting an interesting active space method for triples in which different bonds require different excitation levels. The NLSCC approach recovers a very large percentage (>99%) of the CCSD or CCSDT-3 correlation energy.