Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fock space”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Spin cutoff factor and level density for Ni 59 from an analysis of compound nuclear reactions

The spin cutoff parameter for 59 Ni has been studied from different types of experimental data, including neutron angular distributions from the 56 Fe(α,n) 59 Ni reaction, spin of discrete levels, the level density from the proton evaporation spectrum of the 54 Fe( 6 Li,p) 59 Ni reaction and neutron resonance spacing. Experimental data points were compared with calculations using models widely used in literature. It was found that the available empirical models overestimate data points in the energy region below the neutron separation energy, while microscopic calculations which take into account pairing correlations within Hartree-Fock +Bardeen-Cooper-Schrieffer approach are consistent with data. As a result, it confirmed earlier findings that pairing correlations play an important role and need to be taken into account when the spin cutoff parameter is calculated below the neutron separation energy.

(a,n) reaction↗
DOE OSTI2022

Comparing correlation components and approximations in Hartree–Fock and Kohn–Sham theories via an analytical test case study

The asymmetric Hubbard dimer is a model that allows for explicit expressions of the Hartree–Fock (HF) and Kohn–Sham (KS) states as analytical functions of the external potential, Δv, and of the interaction strength, U. We use this unique circumstance to establish a rigorous comparison between the individual contributions to the correlation energies stemming from the two theories in the {U, Δv} parameter space. Within this analysis of the Hubbard dimer, we observe a change in the sign of the HF kinetic correlation energy, compare the indirect repulsion energies, and derive an expression for the “traditional” correlation energy, i.e., the one that corrects the HF estimate, in a pure site-occupation function theory spirit [Eq. (45)]. Next, we test the performances of the Liu–Burke and the Seidl–Perdew–Levy functionals, which model the correlation energy based on its weak- and strong-interaction limit expansions and can be used for both the traditional and the KS correlation energies. Further, our results show that, in the Hubbard dimer setting, they typically work better for the HF reference, despite having been originally devised for KS. These conclusions are somewhat in line with prior assessments of these functionals on various chemical datasets. However, the Hubbard dimer model allows us to show the extent of the error that may occur in using the strong-interaction ingredient for the KS reference in place of the one for the HF reference, as has been carried out in most of the prior assessments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Building a quantum computing architecture using 3D superconducting cavities

Quantum computers promise advantages over classical machines for solving certain complex problems, but building processors that truly deliver this advantage remains a central challenge, particularly due to limited coherence times. Three-dimensional superconducting radio-frequency (SRF) cavities offer an attractive platform due to their exceptionally long lifetimes. However, since these harmonic systems require nonlinear elements, such as transmons, for control, additional losses are often introduced. In this talk, I will present a multimode quantum system based on an elliptical SRF cavity hosting two cavity modes weakly coupled to an ancillary transmon circuit. This architecture is carefully engineered to preserve coherence while enabling efficient control. By optimizing the design to mitigate transmon-induced decoherence, we realize single-photon lifetimes of 20.6 ms and 15.6 ms in the two modes, with pure dephasing times exceeding 40 ms. Using sideband interactions and error-resilient protocols, such as measurement-based correction and post-selection, we demonstrate high-fidelity state control, including preparation of Fock states up to N=20 with fidelities above 95% (to our knowledge, the highest reported to date), as well as high-fidelity two-mode entanglement. These results highlight 3D SRF cavities as a robust foundation for qudit-based quantum information processing, harnessing the large Hilbert space of cavity modes. I will conclude by outlining strategies to further enhance coherence in both cavities and ancilla qubits, and discuss pathways toward scaling this architecture into a larger quantum computing platform.

Roy, Tanay [Fermilab]↗
DOE OSTI2023

Use of FLOSIC for understanding anion-solvent interactions

An Achille’s heel of lower-rung density-functional approximations is that the highest-occupied-molecular-orbital energy levels of anions, known to be stable or metastable in nature, are often found to be positive in the worst case or above the lowest-unoccupied-molecular-orbital levels on neighboring complexes that are not expected to accept charge. A trianionic example, [Cr(C2O4)3]3−, is of interest for constraining models linking Cr isotope ratios in rock samples to oxygen levels in Earth’s atmosphere over geological timescales. Here we describe how crowd sourcing can be used to carry out self-consistent Fermi–Löwdin–Orbital-Self-Interaction corrected calculations (FLOSIC) on this trianion in solution. The calculations give a physically correct description of the electronic structure of the trianion and water. In contrast, uncorrected local density approximation (LDA) calculations result in approximately half of the anion charge being transferred to the water bath due to the effects of self-interaction error. Use of group-theory and the intrinsic sparsity of the theory enables calculations roughly 125 times faster than our initial implementation in the large N limit reached here. By integrating charge density densities and Coulomb potentials over regions of space and analyzing core-level shifts of the Cr and O atoms as a function of position and functional, we unambiguously show that FLOSIC, relative to LDA, reverses incorrect solute-solvent charge transfer in the trianion-water complex. In comparison to other functionals investigated herein, including Hartree–Fock and the local density approximation, the FLOSIC Cr 1s eigenvalues provide the best agreement with experimental core ionization energies.

Chemistry↗
DOE OSTI2023

An accurate and efficient fragmentation approach via the generalized many-body expansion for density matrices

With relevant chemical space growing larger and larger by the day, the ability to extend computational tractability over that larger space is of paramount importance in virtually all fields of science. The solution we aim to provide here for this issue is in the form of the generalized many-body expansion for building density matrices (GMBE-DM) based on the set-theoretical derivation with overlapping fragments, through which the energy can be obtained by a single Fock build. In combination with the purification scheme and the truncation at the one-body level, the DM-based GMBE(1)-DM-P approach shows both highly accurate absolute and relative energies for medium-to-large size water clusters with about an order of magnitude better than the corresponding energy-based GMBE(1) scheme. Simultaneously, GMBE(1)-DM-P is about an order of magnitude faster than the previously proposed MBE-DM scheme [F. Ballesteros and K. U. Lao, J. Chem. Theory Comput. 18, 179 (2022)] and is even faster than a supersystem calculation without significant parallelization to rescue the fragmentation method. For even more challenging systems including ion–water and ion–pair clusters, GMBE(1)-DM-P also performs about 3 and 30 times better than the energy-based GMBE(1) approach, respectively. In addition, this work provides the first overlapping fragmentation algorithm with a robust and effective binning scheme implemented internally in a popular quantum chemistry software package. Thus, GMBE(1)-DM-P opens a new door to accurately and efficiently describe noncovalent clusters using quantum mechanics.

Chemistry↗
DOE OSTI2021

Modeling quasielastic interactions of monoenergetic kaon decay-at-rest neutrinos

Monoenergetic muon neutrinos at 236 MeV are readily produced in intense medium-energy proton facilities (≳2–3GeV) when a positive kaon decays at rest (KDAR; K + →μ + ν μ ). These neutrinos provide a unique opportunity to both study the neutrino interaction and probe the nucleus with a monoenergetic weak-interaction-only tool. We present cross-section calculations for quasielastic scattering of these 236-MeV neutrinos off 12 C and 40 Ar, paying special attention to low-energy aspects of the scattering process. Our model takes the description of the nucleus in a mean-field approach as the starting point, where we solve Hartree-Fock equations using a Skyrme type nucleon-nucleon interaction. Thereby, we introduce long-range nuclear correlations by means of a continuum random-phase approximation (CRPA) framework where we solve the CRPA equations using a Green's function method. The model successfully describes (e,e') data on 12 C and 40 Ca in the kinematic region that overlaps with the KDAR ν μ phase space. In addition to these results, we present future prospects for precision KDAR cross-section measurements and applications of our calculations in current and future experiments that will utilize these neutrinos.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Cartesian equivariant representations for learning and understanding molecular orbitals

Qualitative and quantitative orbital properties such as bonding/antibonding character, localization, and orbital energies are critical to how chemists understand reactivity, catalysis, and excited-state behavior. Despite this, representations of orbitals in deep learning models have been very underdeveloped relative to representations of molecular geometries and Hamiltonians. Here, we apply state-of-the-art equivariant deep learning architectures to the task of assigning global labels to orbitals, namely energies characterizations, given the molecular coefficients from Hartree–Fock or density functional theory. The architecture we have developed, the Cartesian Equivariant Orbital Network (CEONET), shows how molecular orbital coefficients are readily featurized as equivariant node features common to all graph-based machine-learned potentials. We find that CEONET performs well at predicting difficult quantitative labels such as the orbital energy and orbital entropy. Furthermore, we find that the CEONET representation provides an intuitive latent space for differentiating orbital character for the qualitative assignment of e.g. bonding or antibonding character. In addition to providing a useful representation for further integrating deep learning with electronic structure theory, we expect CEONET to be useful for automatizing and interpreting the results of advanced electronic structure methods such as complete active space self-consistent field theory. In particular, the ability of CEONET to infer multireference character via the orbital entropy paves the way toward the machine-learned selection of active spaces.

chemical reactions↗
DOE OSTI2023

Weak-coupling theory of neutron scattering as a probe of altermagnetism

Inelastic neutron scattering provides a powerful probe of the magnetic excitations of quantum magnets. Altermagnets have recently emerged as a new class of magnets with vanishing net magnetization characteristic of antiferromagnets and with a spin-split electronic structure typical of ferromagnets. Here we introduce a minimal Hubbard model with two-sublattice orthorhombic anisotropy as a framework to study altermagnetism. Using unrestricted Hartree-Fock calculations, we find an altermagnetic state for this model that evolves from a metallic state to an insulating state with increasing Hubbard-U Coulomb repulsion. We then examine the inelastic neutron scattering response in these states using random-phase approximation calculations of the dynamic spin susceptibility χ"(q,ω). We find that the magnetic excitation spectrum depends on its chirality for q along certain directions in reciprocal space, an observation that may be used in inelastic neutron scattering experiments as a probe of altermagnetism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Utilizing Essential Symmetry Breaking in Auxiliary-Field Quantum Monte Carlo: Application to the Spin Gaps of the C 36 Fullerene and an Iron Porphyrin Model Complex

We present three distinct examples where phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) can be reliably performed with a single-determinant trial wave function with essential symmetry breaking. Essential symmetry breaking was first introduced by Lee and Head-Gordon [ Phys. Chem. Chem. Phys. 2019, 21, 4763-4778, 10.1039/C8CP07613H]. We utilized essential complex and time-reversal symmetry breaking with ph-AFQMC to compute the triplet–singlet energy gap in the TS12 set. We found statistically better performance of ph-AFQMC with complex-restricted orbitals than with spin-unrestricted orbitals. We then showed the utilization of essential spin symmetry breaking when computing the singlet–triplet gap of a known biradicaloid, C 36 . ph-AFQMC with spin-unrestricted Hartree–Fock (ph-AFQMC+UHF) fails catastrophically even with spin-projection and predicts no biradicaloid character. With approximate Brueckner orbitals obtained from regularized orbital-optimized second-order Møller–Plesset perturbation theory (κ-OOMP2), ph-AFQMC quantitatively captures strong biradicaloid character of C 36 . Lastly, we applied ph-AFQMC to the computation of the quintet–triplet gap in a model iron porphyrin complex where brute-force methods with a small active space fail to capture the triplet ground state. We show unambiguously that neither triplet nor quintet is strongly correlated using UHF, κ-OOMP2, and coupled-cluster with singles and doubles (CCSD) performed on UHF and κ-OOMP2 orbitals. There is no essential symmetry breaking in this problem. By virtue of this, we were able to perform UHF+ph-AFQMC reliably with a cc-pVTZ basis set and predicted a triplet ground state for this model geometry. The largest ph-AFQMC in this work correlated 186 electrons in 956 orbitals. Our work highlights the utility, scalability, and accuracy of ph-AFQMC with a single-determinant trial wave function with essential symmetry breaking for systems mainly dominated by dynamical correlation with little static correlation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

General, Rigorous Approach for the Treatment of Interfragment Covalent Bonds

Here, a generalized, projection-based transformation of the method-agnostic Fock operator in various ab initio fragment-based quantum chemistry methods has been developed for the treatment of interfragment covalent bonds. This transformation freezes the relevant localized molecular orbital associated with each interfragment bond, thereby restricting the variational subspace of the fragment wave functions, in order to maintain the proper physical characteristics of the involved covalent bonds. In addition, sets of orbitals that would lead to multiple occupancy of certain orbitals are explicitly removed from the variational space. The transformation is developed for the specific case of mutually orthonormal frozen and unfrozen orbitals within each fragment. The newly developed approach is then used to study model systems with two popular ab initio fragment-based methods, and the results of these calculations are compared to those obtained by existing methodologies. Analysis is focused on both quantitative and qualitative accuracy as well as computational scalability and stability. Other methods for which the developed formalisms are appropriate are outlined, and future extensions of the methods are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Valley-spin polarization at zero magnetic field induced by strong hole-hole interactions in monolayer WSe 2

Monolayer transition metal dichalcogenides have emerged as prominent candidates to explore the complex interplay between spin and valley degrees of freedom. Their strong spin-orbit interaction and broken inversion symmetry lead to the spin-valley locking effect, in which carriers occupying theKandK′ valleys of the reciprocal space must have opposite spins. This effect is particularly strong for holes due to a larger spin-orbit gap in the valence band. By reducing the dimensionality of a monolayer of WSe 2 to 1D via electrostatic confinement, we demonstrate that spin-valley locking and strong hole-hole interactions lead to a ferromagnetic state where hole transport is spin-valley polarized, even without an applied magnetic field. A massive Dirac fermion model in the Hartree-Fock approximation reveals that many-body hole-exchange interactions lead to this polarized ground-state. This observation opens the possibility of implementing a robust and stable valley-polarized system, essential in valleytronic applications.

Science & Technology - Other Topics↗
DOE OSTI2023

Phase diagram of twisted bilayer graphene at filling factor $v = ±3$

Here, we study the correlated insulating phases of twisted bilayer graphene (TBG) in the absence of lattice strain at integer filling $v = ±3$. Using the self-consistent Hartree-Fock method on a particle-hole symmetric model and allowing translation symmetry breaking terms, we obtain the phase diagram with respect to the ratio of AA interlayer hopping (w 0 ) and AB interlayer hopping (w 1 ). When the interlayer hopping ratio is close to the chiral limit (w 0 /w 1 ≲ 0.5), a quantum anomalous Hall state with Chern number v c = ±1 can be observed consistent with previous studies. Around the realistic value w 0 /w 1 ≈ 0.8, we find a spin and valley polarized, translation symmetry breaking, state with C 2⁢z T symmetry, a charge gap and a doubling of the moiré unit cell, dubbed theC 2⁢z T stripe phase. The real-space total charge distribution of this C 2⁢z T stripe phase in the flat band limit does not have modulation between different moiré unit cells, although the charge density in each layer is modulated, and the translation symmetry is strongly broken. Other symmetries, including C 2⁢z , C 2⁢x , particle-hole symmetry P, and the topology of the C 2⁢z T stripe phase, are also discussed in detail. We observed braiding and annihilation of the Dirac nodes by continuously turning on the order parameter to its fully self-consistent value, and provide a detailed explanation of the mechanism for the charge gap opening despite preserving C 2⁢z T and valley U⁡(1) symmetries. In the transition region between the quantum anomalous Hall phase and the C 2⁢z T stripe phase, we find an additional competing state with comparable energy corresponding to a phase with a tripling of the moiré unit cell.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

BCS-BEC crossover in a (t 2g ) 4 excitonic magnet

The condensation of spin-orbit-induced excitons in t 2 g 4 electronic systems is attracting considerable attention. At large Hubbard U , antiferromagnetism was proposed to emerge from the Bose-Einstein Condensation (BEC) of triplons ( J eff = 1 ). Furthermore, we show that even at intermediate U regimes, the spin-orbit exciton condensation is possible leading also to staggered magnetic order. The canonical electron-hole excitations (excitons) transform into local triplon excitations at large U , and this BEC strong coupling regime is smoothly connected to the intermediate U excitonic insulator region. We solved the degenerate three-orbital Hubbard model with spin-orbit coupling ( λ ) in one dimension using the density matrix renormalization group, while in two dimensions we use the Hartree-Fock approximation (HFA). Employing these techniques, we provide the full λ versus U phase diagrams for both one- and two-dimensional lattices. Our main result is that at intermediate Hubbard U , increasing λ at fixed U the system transitions from an incommensurate spin-density-wave metal to a Bardeen-Cooper-Schrieffer (BCS) excitonic insulator, with coherence length r coh of O ( a ) and O ( 10 a ) in 1 d and 2 d , respectively, with a being the lattice spacing. Further increasing λ , the system eventually crosses over to the BEC limit (with r coh ≪ a ).

1-dimensional systems↗
DOE OSTI2025

Nonunitary Variational Quantum Eigensolver with the Localized Active Space Method and Cost Mitigation

Accurately describing strongly correlated systems with affordable quantum resources remains a central challenge for quantum chemistry applications on near and intermediate term quantum computers. The localized active space self-consistent field (LASSCF) approximates the complete active space self-consistent field (CASSCF) by generating active space-based wave functions within specific fragments while treating interfragment correlation with mean-field approach, hence is computationally less expensive. Hardware-efficient ansatzes (HEA) offer affordable and shallower circuits, yet they often fail to capture the necessary correlation. Previously, Jastrow-factor-inspired nonunitary qubit operators were proposed to use with HEA for variational quantum eigensolver (VQE) calculations (so-called nuVQE), as they do not increase circuit depths and recover correlation beyond the mean-field level for Hartree–Fock initial states. Here, in this study, we explore running nuVQE with LASSCF as the initial state. The method, named LAS-nuVQE, is shown to recover interfragment correlations, reach chemical accuracy with a small number of gates (<70) in both H 4 and square cyclobutadiene (C 4 H 4 ), and produces more accurate energetics than its HEA counterparts at all circuit depths. To further address the inherent symmetry-breaking in HEA, we implemented spin-constrained LAS-nuVQE to extend the capabilities of HEA further and show spin-pure results for square cyclobutadiene. We also mitigate the increased measurement overhead of nuVQE via Pauli grouping and shot-frugal sampling, reducing measurement costs by up to 2 orders of magnitude compared to ungrouped operator, and show that one can achieve better accuracy with a small number of shots (10 3–4 ) per one expectation value calculation compared to noiseless simulations with one or two orders of magnitude more shots. Finally, wall clock time estimates show that, with our measurement mitigation protocols, nuVQE becomes a cheaper and more accurate alternative than vanilla VQE with HEA. Taken together, these developments illustrate a practical pathway toward performing multireference chemical simulations with accuracy and affordable resources on today’s quantum hardware, achieving both accuracy and affordability in challenging correlated systems.

Wang, Qiaohong [Univ. of Chicago, IL (United State↗
DOE OSTI2020

Constraining level densities through quantitative correlations with cross-section data

The adopted level densities (LD) for the nuclei produced through different reaction mechanisms significantly impact the accurate calculation of cross sections for the different reaction channels. Many common LD models make simplified assumptions regarding the overall behavior of the total LD and the intrinsic spin and parity distributions of the excited states. However, very few experimental constraints are taken into account in these models: LD at neutron separation energy coming from average spacings of s- and p-wave resonances ( D 0 and D 1 , respectively) whenever they have been previously measured, and the sometimes subjective extrapolation of discrete levels. These, however, constrain the LD only in very specific regions of excitation energy, and for specific spins and parities. This work aims to establish additional experimental constraints on LD through quantitative correlations between cross sections and LD. This allows for fitting and the determination of detailed structures in LD. For this we use the microscopic Hartree-Fock-Bogoliubov (HFB) LD as a starting point as the HFB LD provide a more realistic spin and parity distributions than phenomenological models such as Gilbert-Cameron (GC). We then associate variations predicted by the HFB model with the structure observed in double-differential cross sections at low outgoing neutron energy, a region that is dominated by the LD input. We also use (n, p) on 56 Fe , as an example case where angle-integrated cross sections are extremely sensitive to LD. For comparison purposes we also perform calculations with the GC model. With this approach we are able to perform fits of the LD based on actual experimental data, constraining the model and ensuring its consistency. This approach can be particularly useful in extrapolating the LD to nuclei for which high-excited discrete levels and/or values of D 0 or D 1 are unknown. Finally, it also predicts neutron-induced inelastic γ cross sections that in some cases can differ significantly from more standard phenomenological LD models such as GC.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Generative deep-learning reveals collective variables of Fermionic systems

Complex processes of fermionic systems ranging from protein folding to nuclear fission often follow a low-dimensional reaction path parametrized in terms of a few collective variables. In nuclear theory, variables related to the shape of the nuclear density in a mean-field picture are key to describing the large amplitude collective motion of the neutrons and protons. Exploring the adiabatic energy landscape spanned by these degrees of freedom reveals the possible reaction channels while simulating the dynamics in this reduced space yields their respective probabilities. Unfortunately, this theoretical framework breaks down whenever the systems encounters a quantum phase transition with respect to the collective variables. Here, in this study, we introduce a novel generative deep-learning algorithm designed to build reaction paths that ensure that the many-fermion wave function stays differentiable with respect to the collective variables. This approach is applicable to any fermionic system described by a coherent state. We use the case of potential energy curves in the 16 O nucleus within the Hartree-Fock theory to illustrate its main features.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Møller–Plesset and Density-Fixed Adiabatic Connections for a Model Diatomic System at Different Correlation Regimes

In recent years, adiabatic connection (AC) interpolations developed within density functional theory (DFT) have been found to provide good performances in the calculation of interaction energies when used with Hartree–Fock (HF) ingredients. The physical and mathematical reasons for such unanticipated performance have been clarified, to some extent, by studying the strong-interaction limit of the Møller–Plesset (MP) AC. In this work, we calculate both the MP and the DFT AC integrand for the asymmetric Hubbard dimer, which allows for a systematic investigation of different correlation regimes by varying two simple parameters in the Hamiltonian: the external potential, Δv, and the interaction strength, U. Notably, we find that, while the DFT AC integrand appears to be convex in the full parameter space, the MP integrand may change curvature twice. Furthermore, we discuss different aspects of the second-order expansion of the correlation energy in each AC, and we demonstrate why the derivative of the λ-dependent density in the MP AC at λ = 0 (i.e., at the HF density) is zero in the model. Concerning the strong-interaction limit of both ACs in the Hubbard dimer setting, we show that the asymptotic value of the MP AC, W ∞ HF , is lower than (or equal to) its DFT analogue, W ∞ KS , if the two are compared at a given density, just like in real space. Furthermore, we also show that this is not always the case if the two quantities are compared at a given external potential.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Hybrid gausslet/Gaussian basis sets

We introduce hybrid gausslet/Gaussian basis sets, where a standard Gaussian basis is added to a gausslet basis in order to increase accuracy near the nuclei while keeping the spacing of the grid of gausslets relatively large. The Gaussians are orthogonalized to the gausslets, which are already orthonormal, and approximations are introduced to maintain the diagonal property of the two electron part of the Hamiltonian so that it continues to scale as the second power of the number of basis functions rather than the fourth. We introduce several corrections to the Hamiltonian designed to enforce certain exact properties, such as the values of certain two-electron integrals. We also introduce a simple universal energy correction that compensates for the incompleteness of the basis stemming from the electron–electron cusps based on the measured double occupancy of each basis function. We perform a number of Hartree Fock and full configuration interaction (full-CI) test calculations on two electron systems and Hartree Fock on a ten-atom hydrogen chain to benchmark these techniques. The inclusion of the cusp correction allows us to obtain complete basis set full-CI results for the two electron cases at the level of several micro-Hartrees, and we see similar apparent accuracy for Hartree Fock on the ten-atom hydrogen chain.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗