Search

A summary of the technical advances that are incorporated in the fourth major release of the Q-Chem quantum chemistry program is provided, covering approximately the last seven years. These include developments in density functional theory methods and algorithms, nuclear magnetic resonance (NMR) property evaluation, coupled cluster and perturbation theories, methods for electronically excited and open-shell species, tools for treating extended environments, algorithms for walking on potential surfaces, analysis tools, energy and electron transfer modelling, parallel computing capabilities, and graphical user interfaces. In addition, a selection of example case studies that illustrate these capabilities is given. These include extensive benchmarks of the comparative accuracy of modern density functionals for bonded and non-bonded interactions, tests of attenuated second order Møller-Plesset (MP2) methods for intermolecular interactions, a variety of parallel performance benchmarks, and tests of the accuracy of implicit solvation models. Some specific chemical examples include calculations on the strongly correlated Cr2 dimer, exploring zeolite-catalysed ethane dehydrogenation, energy decomposition analysis of a charged ter-molecular complex arising from glycerol photoionisation, and natural transition orbitals for a Frenkel exciton state in a nine-unit model of a self-assembling nanotube.

We demonstrate reactivity between a β-diketiminate-supported niobium(III) imido complex and alkyl azides to form niobatetrazene complexes (BDI)Nb(N t Bu)(RNNNNR) (BDI = N , N -bis(2,6-diisopropylphenyl)-3,5-dimethyl-β-diketiminate; R = cyclohexyl ( 1 ), benzyl ( 2 )). Intriguingly, niobatetrazene complexes 1 and 2 can be interconverted via addition of an appropriate alkyl azide, likely through a series of concerted [3 + 2] cycloaddition and retrocycloaddition reactions in which π-loaded bis(imido) intermediates are formed. The bis(imido) intermediates were trapped upon addition of alkyl isocyanides to yield five-coordinate bis(imido) complexes (BDI)Nb(N t Bu)(NCy)(CNR) (R = tert -butyl ( 4a ), cyclohexyl ( 4b )). Two computational methods─density functional theory and density functional tight binding (DFTB)─were employed to calculate the lowest energy pathway across the potential energy surface for this multistep transformation. Reaction path calculations for individual cycloaddition or retrocycloaddition processes along the multistep reaction pathway showed that these transformations occur via a concerted, yet highly asynchronous mechanism, in which the two bond-breaking or -making events do not occur simultaneously. The use of the DFTB method in this work highlights its advantages and utility for studying transition metal systems.

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.

We investigate the interaction between a molecule and a pore mouth-a critical step in adsorption processes-by characterizing the conformation of a macrocyclic calix[4]arene-Ti IV complex, which is grafted on the external surface of a zeotype (*-SVY). X-ray absorption and 13 C{ 1 H} CPMAS NMR spectroscopies independently detect a unique conformation of this complex when it is grafted at crystallographically equivalent locations that lie at the interface of 7 Å hemispherical microporous cavities and the external surface. Electronic structure calculations support the presence of this unique conformation, and suggest that it is brought about by a specific orientation of the macrocycle that maximizes non-covalent interactions between calix[4]arene upper-rim tert-butyl substituents and the microporous-cavity walls. Our comparative study provides a rare "snapshot" of a molecule partially confined at a pore mouth, an essential intermediate for adsorption into micropores, and demonstrates how surrounding environment controls this confinement in a sensitive fashion., (© 2021 Wiley-VCH GmbH.)

Electron paramagnetic resonance (EPR) studies of the rhenium(II) complex Re(η 5 -Cp)(BDI) ( 1 ; BDI = N , N '-bis(2,6-diisopropylphenyl)-3,5-dimethyl-β-diketiminate) have revealed that this species reversibly binds N 2 in solution: flash frozen toluene solutions of 1 disclose entirely different EPR spectra at 10 K when prepared under N 2 versus Ar atmospheres. This observation was additionally verified by the synthesis of stable CO and 2,6-xylylisocyanide (XylNC) adducts of 1 , which display EPR features akin to those observed in the putative N 2 complex. While we found that 1 displays an extremely large g max value of 3.99, the binding of an additional ligand leads to substantial decreases in this value, displaying g max values of ca. 2.4. Following the generation of isotopically enriched 15 N 2 , HYSCORE experiments allowed for the measurement of the corresponding hyperfine couplings associated with spin delocalization onto the electron-accepting ligands in these species, which proved to be small. A cumulative assessment of the EPR data, when combined with insights provided by near-infrared (NIR) spectroscopy and time-dependent density functional theory (TDDFT) calculations, indicated that while the binding of electron acceptors to 13 CO adducts of 1 , HYSCORE experiments allowed for the measurement of the corresponding hyperfine couplings associated with spin delocalization onto the electron-accepting ligands in these species, which proved to be small. A cumulative assessment of the EPR data, when combined with insights provided by near-infrared (NIR) spectroscopy and time-dependent density functional theory (TDDFT) calculations, indicated that while the binding of electron acceptors to 1 does lead to decreases in g max in relative accord with the field strength (i.e., π-acidity) of the variable ligand, the magnitude of these decreases is primarily due to the changes in electronic structure at the Re center.

The focus of this work is OS-CCSD-SPT(2), which is a second-order similarity transformed perturbation theory correction to opposite spin coupled cluster singles doubles, where in the latter the same-spin amplitudes are removed and the opposite-spin ones are solved self-consistently. OS-CCSD-SPT(2) is free of empirical parameters, has an instrinsic scaling of O ( N 6 ), and makes no use of triples. We demonstrate that, for non-multireference molecules, OS-CCSD-SPT(2) produces relative energies whose accuracy is significantly higher than what is generally expected of a triples-free model. For example, using PBE0 orbitals in the reference, OS-CCSD-SPT(2) exhibits a mean absolute deviation (MAD) of 1.13 kcal/mol with respect to CCSD(2 F ) benchmark values for the non-multireference subset of W4-08 atomization energies (cf. a MAD > 6.5 kcal/mol for CCSD) and a MAD of 0.68 kcal/mol for the energies of reactions generated from the W4-08 molecules. These MADs are reduced to 0.61 and 0.63 kcal/mol, respectively, by a simple one-parameter spin-component scaling of the OS-CCSD-SPT(2) same-spin correlation energy. OS-CCSD is also naturally amenable to higher order corrections: the associated third-order correction, OS-CCSD-SPT(3), which does involve connected triples and quadruples, exhibits a MAD of 0.44 kcal/mol for the same atomization-energy benchmark.

Machine learning (ML) approximations to density functional theory (DFT) potential energy surfaces (PESs) are showing great promise for reducing the computational cost of accurate molecular simulations, but at present, they are not applicable to varying electronic states, and in particular, they are not well suited for molecular systems in which the local electronic structure is sensitive to the medium to long-range electronic environment. With this issue as the focal point, we present a new machine learning approach called "BpopNN" for obtaining efficient approximations to DFT PESs. Conceptually, the methodology is based on approaching the true DFT energy as a function of electron populations on atoms; in practice, this is realized with available density functionals and constrained DFT (CDFT). The new approach creates approximations to this function with neural networks. These approximations thereby incorporate electronic information naturally into a ML approach, and optimizing the model energy with respect to populations allows the electronic terms to self-consistently adapt to the environment, as in DFT. We confirm the effectiveness of this approach with a variety of calculations on Li n H n clusters.

Anti-Markovnikov additions to alkenes have been a longstanding goal of catalysis, and anti-Markovnikov addition of arenes to alkenes would produce alkylarenes that are distinct from those formed by acid-catalysed processes. Existing hydroarylations are either directed or occur with low reactivity and low regioselectivity for the n-alkylarene. Herein, we report the first undirected hydroarylation of unactivated alkenes with unactivated arenes that occurs with high regioselectivity for the anti-Markovnikov product. The reaction occurs with a nickel catalyst ligated by a highly sterically hindered N-heterocyclic carbene. Catalytically relevant arene- and alkene-bound nickel complexes have been characterized, and the rate-limiting step was shown to be reductive elimination to form the C-C bond. Density functional theory calculations, combined with second-generation absolutely localized molecular orbital energy decomposition analysis, suggest that the difference in activity between catalysts containing large and small carbenes results more from stabilizing intramolecular non-covalent interactions in the secondary coordination sphere than from steric hindrance.

A dicopper complex featuring a symmetrically bridging nitrile ligand and supported by a binucleating naphthyridine-based ligand, [Cu 2 ( μ-η 1 : η 1 -MeCN)DPFN](NTf 2 ) 2 , was treated with phosphaalkynes (RC[triple bond, length as m-dash]P, isoelectronic analogues of nitriles) to yield dicopper complexes that exhibit phosphaalkynes in rare μ-η 2 : η Bu (up). With intermediate steric bulk (R = Ad), the energy difference between the two geometries is small enough that both are observed in solution, and NMR spectroscopy and computations indicate that the solid-state structure corresponds to the minor isomer observed in solution. Meanwhile, treatment of [Cu 2 binding coordination modes. X-ray crystallography revealed that these unusual "tilted" structures exist in two isomeric forms (R "up" vs. R "sideways"), depending on the steric profile of the phosphaalkyne's alkyl group (R = Me, Ad, or t Bu). Only one isomer is observed in both solution and the solid state for R = Me (sideways) and t Bu (up). With intermediate steric bulk (R = Ad), the energy difference between the two geometries is small enough that both are observed in solution, and NMR spectroscopy and computations indicate that the solid-state structure corresponds to the minor isomer observed in solution. Meanwhile, treatment of [Cu 2 ( μ-η 1 : η 1 -MeCN)DPFN](NTf 2 ) 2 with 2-butyne affords [Cu 2 ( μ-η 2 : η 2 -(MeC[triple bond, length as m-dash]CMe))DPFN](NTf 2 ) 2 : its similar ligand geometry demonstrates that the tilted μ-η 2 : η 2 binding mode is not limited to phosphaalkynes but reflects a more general trend, which can be rationalized via an NBO analysis showing maximization of π-backbonding., (This journal is © The Royal Society of Chemistry 2020.)

In this work, we revisited the idea of using the coupled-cluster (CC) ground state formalism to target excited states. Our main focus was targeting doubly excited states and double core hole states. Typical equation-of-motion (EOM) approaches for obtaining these states struggle without higher-order excitations than doubles. We showed that by using a non-Aufbau determinant optimized via the maximum overlap method, the CC ground state solver can target higher energy states. Furthermore, just with singles and doubles (i.e., CCSD), we demonstrated that the accuracy of ΔCCSD and ΔCCSD(T) (triples) far surpasses that of EOM-CCSD for doubly excited states. The accuracy of ΔCCSD(T) is nearly exact for doubly excited states considered in this work. For double core hole states, we used an improved ansatz for greater numerical stability by freezing core hole orbitals. The improved methods, core valence separation (CVS)-ΔCCSD and CVS-ΔCCSD(T), were applied to the calculation of the double ionization potential of small molecules. Even without relativistic corrections, we observed qualitatively accurate results with CVS-ΔCCSD and CVS-ΔCCSD(T). Remaining challenges in ΔCC include the description of open-shell singlet excited states with the single-reference CC ground state formalism as well as excited states with genuine multireference character. The tools and intuition developed in this work may serve as a stepping stone toward directly targeting arbitrary excited states using ground state CC methods.

We show that using complex, spin-restricted orbitals in Kohn-Sham (KS) density functional theory allows one to access a new class of densities that is not accessible by either spin-restricted (RKS) or spin-unrestricted (UKS) orbitals. We further show that the real part of a complex RKS (CRKS) density matrix can be nonidempotent when the imaginary part of the density matrix is not zero. Using CRKS orbitals shows significant improvements in the triplet-singlet gaps of a benchmark set, called TS12, for well-established, widely used density functionals. Moreover, it was shown that RKS and UKS yield qualitatively wrong charge densities and spin densities, respectively, leading to worse energetics. We demonstrate that representative modern density functionals show surprisingly no improvement even with a qualitatively more accurate density from CRKS orbitals. To this end, our work not only provides a way to escape the symmetry dilemma whenever there exists a CRKS solution, but also suggests a new route to design better approximate density functionals.

The effect of dynamic reorganization and confinement of isolated Ti IV catalytic centers supported on silicates is investigated for olefin epoxidation. Active sites consist of grafted single-site calix[4]arene-Ti IV centers or their calcined counterparts. Their location is synthetically controlled to be either unconfined at terminal T-atom positions (denoted as type-(i)) or within confining 12-MR pockets (denoted as type-(ii); diameter ∼7 Å, volume ∼185 Å 3 ) composed of hemispherical cavities on the external surface of zeotypes with *-SVY topology. Electronic structure calculations (density functional theory) indicate that active sites consist of cooperative assemblies of Ti IV centers and silanols. When active sites are located at unconfined type-(i) environments, the rate constants for cyclohexene epoxidation (323 K, 0.05 mM Ti IV , 160 mM cyclohexene, 24 mM tert-butyl hydroperoxide) are 9 ± 2 M -2 s -1 ; whereas within confining type-(ii) 12-MR pockets, there is a ∼5-fold enhancement to 48 ± 8 M -2 s -1 . When a mixture of both environments is initially present in the catalyst resting state, the rate constants reflect confining environments exclusively (40 ± 11 M -2 s -1 ), indicating that dynamic reorganization processes lead to the preferential location of active sites within 12-MR pockets. While activation enthalpies are Δ H ‡ app = 43 ± 1 kJ mol -1 irrespective of active site location, confining environments exhibit diminished entropic barriers (Δ S ‡ app = -68 J mol -1 K -1 for unconfined type-(i) vs -56 J mol -1 K -1 for confining type-(ii)), indicating that confinement leads to more facile association of reactants at active sites to form transition state structures (volume ∼ 225 Å 3 ). These results open new opportunities for controlling reactivity on surfaces through partial confinement on shallow external-surface pockets, which are accessible to molecules that are too bulky to benefit from traditional confinement within micropores.

We report the failure of coupled-cluster valence-bond (CCVB) theory with two-pair configurations [D. W. Small and M. Head-Gordon, J. Chem. Phys. 130 , 084103 (2009)] for open-shell (OS) spin-frustrated systems where including three-pair configurations is necessary to properly describe strong spin-correlations. We extend OS-CCVB by augmenting the model with three-pair configurations within the independent amplitude approximation. The resulting new electronic structure model, OS-CCVB+i3, involves only a quadratic number of independent wavefunction parameters. It includes the recently reported closed-shell CCVB+i3 as a special case. Its cost is dominated by integral transformations, and it is capable of breaking multiple bonds exactly for all systems examined so far. The strength of OS-CCVB+i3 is highlighted in realistic systems including the [CaMn 3 O 4 ] cubane subunit of the oxygen-evolving complex and a molecular magnet with the [Cr 9 ] core unit as well as model systems such as N 3 , V 3 O 3 , and P 5 . We show that OS-CCVB+i3 is only slightly dependent on the underlying perfect-pairing reference, while OS-CCVB shows a stronger dependence. We also emphasize the compactness of the OS-CCVB+i3 wavefunction compared to the heat-bath configuration interaction wavefunction, a recently introduced soft exponential-scaling approach.

Coupled cluster valence bond (CCVB) is a simple electronic structure method based on a perfect pairing (PP) reference with 2-pair recouplings for strong electron correlation problems. CCVB is spin-pure, size-consistent, and can exactly (in its active space) separate any molecule into atoms for which unrestricted Hartree-Fock (UHF) at dissociation is the sum of the ground state UHF energies of the atoms. However CCVB is far from a complete description of strong correlations. Its first failure to exactly describe spin-recouplings arises at the level of 3 electron pairs, such as the recoupling of 3 triplet oxygen atoms in the dissociation of singlet ozone. Such situations are often associated with spin frustration. To address this limitation, an extension of CCVB, termed CCVB+i3, is reported here that includes an independent (i) amplitude approximation to the 3-pair recouplings. CCVB+i3 thereby has the same basic computational requirements as those of CCVB, which has previously been shown to be an efficient method. CCVB+i3 correctly separates molecules that CCVB cannot. As a by-product, an independent 2-pair amplitude approximation to CCVB, called PP+i2, is also defined. Remarkably, PP+i2 can also correctly separate systems that CCVB cannot. CCVB+i3 is validated on the symmetric dissociation of D 3 h ozone. CCVB+i3 is then used to explore the role of 3-pair recouplings in an [Fe 4 S 4 (SCH 3 ) 4 ] 2- cluster that has been used to model the iron-sulfur core of [Fe 4 S 4 ] ferredoxins. Using localized PP orbitals, such recouplings are demonstrated to be large in some low-lying singlet excited states of the cluster. Significant 3 pair recoupling amplitudes include the usual triangular motif associated with spin frustration and other geometric arrangements of the 3 entangled pairs across the 4 iron centers.

The Coupled Cluster Valence Bond (CCVB) method, previously presented for closed-shell (CS) systems, is extended to open-shell (OS) systems. The theoretical development is based on embedding the basic OS CCVB wavefunction in a fictitious singlet super-system. This approach reveals that the OS CCVB amplitude equations are quite similar to those of CS CCVB, and thus that OS CCVB requires the same level of computational effort as CS CCVB, which is an inexpensive method. We present qualitatively correct CCVB potential energy curves for all low-lying spin states of P 2 and Mn 2 + . CCVB is successfully applied to the low-lying spin states of some model linear polycarbenes, systems that appear to be a hindrance to standard density functionals. We examine an octa-carbene dimer in a side-by-side orientation, which, in the monomer dissociation limit, exhibits maximal strong correlation over the length of the polycarbene.

We demonstrate a block-tensor based implementation of coupled-cluster valence-bond singles and doubles (CCVB-SD) [Small, D. W.; Head-Gordon M. J. Chem. Phys. 2012, 137, 114103] which is a simple modification to restricted CCSD (RCCSD) that provides a qualitatively correct description of valence correlations even in strongly correlated systems. We derive the Λ-equation of CCVB-SD and the corresponding unrelaxed density matrices. The resulting production-level implementation is applied to oligoacenes, correlating up to 318 electrons in 318 orbitals. CCVB-SD shows a qualitative agreement with exact methods for short acenes and reaches the bulk limit of oligoacenes in terms of natural orbital occupation numbers, whereas RCCSD shows nonvariational behavior even for relatively short acenes. A significant reduction in polyradicaloid character is found when correlating all valence electrons instead of only the π-electrons.

We introduce a necessary and sufficient condition for an arbitrary wavefunction to be collinear, i.e., its spin is quantized along some axis. It may be used to obtain a cheap and simple computational procedure to test for collinearity in electronic structure theory calculations. We adapt the procedure for Generalized Hartree Fock (GHF), and use it to study two dissociation pathways in CO2. For these dissociation processes, the GHF wave functions transform from low-spin Unrestricted Hartree Fock (UHF) type states to noncollinear GHF states and on to high-spin UHF type states, phenomena that are succinctly illustrated by the constituents of the collinearity test. This complements earlier GHF work on this molecule.

Restricted Hartree Fock using complex-valued orbitals (cRHF) is studied. We introduce an orbital pairing theorem, with which we obtain a concise connection between cRHF and real-valued RHF, and use it to uncover the close relationship between cRHF, unrestricted Hartree Fock, and generalized valence bond perfect pairing. This enables an intuition for cRHF, contrasting with the generally unintuitive nature of complex orbitals. We also describe an efficient computer implementation of cRHF and its corresponding stability analysis. By applying cRHF to the Be + H2 insertion reaction, a Woodward-Hoffmann violating reaction, and a symmetry-driven conical intersection, we demonstrate in genuine molecular systems that cRHF is capable of removing certain potential energy surface singularities that plague real-valued RHF and related methods. This complements earlier work that showed this capability in a model system. We also describe how cRHF is the preferred RHF method for certain radicaloid systems like singlet oxygen and antiaromatic molecules. For singlet O2, we show that standard methods fail even at the equilibrium geometry. An implication of this work is that, regardless of their individual efficacies, cRHF solutions to the HF equations are fairly commonplace.

We describe an efficient implementation of the coupled cluster valence bond (CCVB) model. CCVB captures a certain essential part of the description of molecules with strong correlations (SC), which allows it to achieve correct energy profiles when covalent bonds are broken, while maintaining proper spin symmetry and size extensivity. To illustrate treatment of SC in bond breaking, we examine the symmetric dissociation of the sulfur allotropes S6 and S8 into triplet S atoms. To show applicability to larger systems and to explore whether CCVB can capture aspects of SC that arise in extended π systems, we report results for a series of acenes up to 12 fused benzene rings, with active spaces of up to 228 correlated electrons. The lowest-energy CCVB solutions found for two of the largest acenes show signatures consistent with multi-electron SC and partial delocalization.

Closed-shell coupled cluster singles and doubles (CCSD) is among the most important of electronic-structure methods. However, it fails qualitatively when applied to molecular systems with more than two strongly correlated electrons, such as those with stretched or broken covalent bonds. We show that it is possible to modify the doubles amplitudes to obtain a closed-shell CCSD method that retains the computational cost and desirable features of standard closed-shell CCSD, e.g., correct spin symmetry, size extensivity, orbital invariance, etc., but produces greatly improved energies upon bond dissociation of multiple electron pairs; indeed, under certain conditions the dissociation energies are exact.

We use time-resolved interferometry, fluorescence assays and computational fluid dynamics (CFD) simulations to examine the viability of confluent adherent cell monolayers to selection via laser microbeam release of photoresist polymer micropallets. We demonstrate the importance of laser microbeam pulse energy and focal volume position relative to the glass-pallet interface in governing the threshold energies for pallet release as well as the pallet release dynamics. Measurements using time-resolved interferometry show that increases in laser pulse energy result in increasing pallet release velocities that can approach 10 m s(-1) through aqueous media. CFD simulations reveal that the pallet motion results in cellular exposure to transient hydrodynamic shear stress amplitudes that can exceed 100 kPa on microsecond timescales, and which produces reduced cell viability. Moreover, CFD simulation results show that the maximum shear stress on the pallet surface varies spatially, with the largest shear stresses occurring on the pallet periphery. Cell viability of confluent cell monolayers on the pallet surface confirms that the use of larger pulse energies results in increased rates of necrosis for those cells situated away from the pallet centre, while cells situated at the pallet centre remain viable. Nevertheless, experiments that examine the viability of these cell monolayers following pallet release show that proper choices for laser microbeam pulse energy and focal volume position lead to the routine achievement of cell viability in excess of 90 per cent. These laser microbeam parameters result in maximum pallet release velocities below 6 m s(-1) and cellular exposure of transient hydrodynamic shear stresses below 20 kPa. Collectively, these results provide a mechanistic understanding that relates pallet release dynamics and associated transient shear stresses with subsequent cellular viability. This provides a quantitative, mechanistic basis for determining optimal operating conditions for laser microbeam-based pallet release systems for the isolation and selection of adherent cells.

We give a pedagogical overview of our recently introduced electronic-structure method, Coupled Cluster Valence Bond (CCVB). We show that CCVB can be viewed as an approximation to the accurate, yet very expensive, Spin Coupled Valence Bond model (SCVB). Both of these models are intended for use on strongly correlated molecular systems, especially when the strong correlations are due to electron spin coupling. Using familiar ideas from electronic-structure theory, we provide definitions for these strong-correlation concepts. We show that CCVB and SCVB generally produce similar results, with more substantial discrepancies occurring for systems displaying electronic resonance. We conclude that CCVB is a useful, inexpensive alternative to SCVB., (This journal is © the Owner Societies 2011)

Spin unrestriction is typically defined as a free variation of the molecular orbitals of different spins in order to lower the molecular energy. When applied to approximate active space electron correlation methods, such as generalized valence bond coupled cluster (GVB-CC) models, this approach often leads to undesirable artifacts. We therefore present an alternative unrestricted in active pairs (UAP) procedure for spin polarization in GVB-CC methods, which resembles the corresponding orbitals of unrestricted Hartree-Fock theory. The UAP procedure permits spin polarization only within the two orbital subspaces describing each electron pair (consisting of one nominally occupied and one virtual orbital). This great reduction to just a linear number of degrees of freedom associated with spin polarization eliminates many of the undesirable artifacts associated with unconstrained spin unrestriction. The UAP procedure is tested on a variety of potential curves for bond breaking and the properties of small radicals as well as larger polyenyl radicals, allyl, pentadienyl, and heptatrienyl.

For a given number of electrons, total spin, and matching spin z-component, we construct a set that spans the many-electron spin subspace associated with these spin values. Each vector in the set is tensorially related to spin-pure vectors of six electrons or less. We show that in the limit of separated atoms coupled to any allowed overall spin, the corresponding spin vector has a simple form relative to the introduced sets. From this, we set up a model that is computationally simple, spin pure, size consistent, and able to properly treat molecules as they dissociate into atoms or fragments.

We develop a correction for the coupled cluster version of the perfect pairing (PP) model. The correction is based on finding modified values of the PP amplitudes such that the second coupled cluster central moment defined in the space of all valence single and double substitutions vanishes and, subject to this constraint, minimizing the deviation between the modified and unmodified PP amplitudes with respect to a chosen metric. We discuss how this correction can be generalized to other constrained doubles models, such as local correlation and active-space models. While the correction is not strictly size consistent and retains some of the deficiencies of the PP model, numerical results indicate that much of the missing active-space coupled cluster singles and doubles correlation energy is recovered.

Molecular dynamics simulations are presented for condensed-phase electron transfer (ET) systems where the electronic polarizability of both the solvent and the solute is incorporated. The solute polarizability is allowed to change with electronic transition. The results display notable deviation from the standard free energy parabolas of traditional ET theories. A new three-parameter ET model is applied, and the theory is shown to accurately model the free energy surfaces. This paper presents conclusive evidence that the traditional theory for the free energy barrier of ET reactions requires modification.

Bransford and Franks (1971) showed that people retain the ideas expressed in linguistic input, rather than retaining the individual sentences. The extent to which this phenomenon reflects organizational' processes based on co-occurrence information was investigated by using Bransford and Franks' procedure but substituting arbitrarily chosen groups of words for sentences. The subjects gave higher recognition ratings to large new groups (analogous to complex sentences) than to small new groups (analogous to simpler ones), so long as the size of the group was in the range normally presented. They also showed the same ordering for old groups; but they rated old groups higher than new groups within each group size. Hierarchical recognition response structure can be established without an ideational basis, and it does not require complete absence from memory of original input items.