Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “INTERMEDIATE COUPLING”

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 1,585 records · Page 88

Coupled geomechanical investigation of depletion-induced fault reactivation

Fault reactivation during subsurface fluid production pose significant challenges to safe and sustainable resource extraction. Here, this study presents a three-dimensional coupled geomechanical framework to investigate the processes driving fault reactivation, capturing the interactions between reservoir dynamics and geomechanical responses. Verification against theoretical estimations based on linear poroelasticity confirms the model's capacity in representing reservoir background stress responses. However, the study reveals that relying solely on background stress states can underestimate or overestimate fault reactivation potential, emphasizing the importance of including localized stress perturbations such as differential compaction and stress redistribution. Applied to a fault (M1) inspired by the geological characteristics of the Groningen field, the model shows slip initiation at 2965 m depth with 16.0 MPa depletion, aligning with field observations where seismicity occurred at approximately 3 km depth after 15.8 MPa depletion. Parametric studies reveal: (1) inelastic reservoir compaction delays fault reactivation and mitigates fault slip by reducing stress concentration, (2) higher intermediate in-situ stress magnitudes decrease the Coulomb Failure Stress (CFS) increase rate and reduce fault slip, (3) larger fault offsets amplify shear stress near the offset zone, promoting earlier reactivation and longer rupture propagation, and (4) fault permeability significantly influences pressure diffusion, with low-permeability faults leading to sharper stress changes and earlier fault destabilization. These insights highlight the critical role of geological and mechanical parameters in fault reactivation and provide a predictive framework for mitigating induced seismicity risks.

Coupled geomechanical modeling↗

Electrochemical Oxidation of Organic Molecules at Lower Overpotential: Accessing Broader Functional Group Compatibility with Electron-Proton Transfer Mediators

CONSPECTUS: Electrochemical organic oxidation reactions are highly appealing because protons are often effective terminal electron acceptors, thereby avoiding undesirable stoichiometric oxidants. These reactions are often plagued by high overpotentials, however, that greatly limit their utility. Single-electron transfer (SET) from organic molecules generates high energy radical-cations. Formation of such intermediates often requires electrode potentials far above the thermodynamic potentials of the reaction and frequently causes decomposition and/or side reactions of ancillary functional groups. In this Account, we show how electrocatalytic electron-proton transfer mediators (EPTMs) address this challenge. EPTMs bypass the formation of radical-cation intermediates by supporting mechanisms that operate at electrode potentials much lower (=1 V) than analogous direct electrolysis reactions. The stable aminoxyl radical TEMPO (2,2,6,6-tetramethylpiperidine N-oxyl) is an effective mediator for electrochemical alcohol oxidation, and we have employed such processes for applications ranging from pharmaceutical synthesis to biomass conversion. A complementary electrochemical alcohol oxidation method employs a cooperative Cu/TEMPO mediator system that operates at 0.5 V lower electrode potential than the TEMPO-only mediated process. This difference, which arises from a different catalytic mechanism, rationalizes the broad functional group tolerance of Cu/TEMPO-based aerobic alcohol oxidation catalysts. Aminoxyl mediators address long-standing challenges in the "Shono oxidation", an important method for a-C–H oxidation of tertiary amides and carbamates. Shono oxidations are initiated by a high-potential SET step that limits their utility. Aminoxyl-mediated Shono-type oxidations have been developed that operate at much lower potentials and tolerate diverse functional groups. Analogous reactivity underlies a-C–H cyanation of secondary cyclic amines, a new method that enables efficient diversification of piperidine-based pharmaceutical building blocks and preparation of non-natural amino acids. Electrochemical oxidations of benzylic C–H bonds are commonly initiated by SET to generate radical cations, but such methods are again plagued by large overpotentials. Mediated electrolysis methods that promote hydrogen-atom-transfer (HAT) from benzylic C–H bonds to Fe-oxo species and phthalimide N-oxyl (PINO) support C–H oxygenation, iodination, and oxidative-coupling reactions. A complementary method merges photochemistry with electrochemistry to achieve amidation of C(sp3)–H bonds. This unique process operates at much lower overpotentials compatible with diverse functional groups. These results have broad implications for organic electrochemistry, highlighting the importance of "overpotential" considerations and the prospects for expanding synthetic utility by using mediators to bypass high-energy outer-sphere electron-transfer mechanisms. Principles demonstrated here for oxidation are equally relevant to electrochemical reductions.

Wang, Fei↗

Probing surface-adsorbate interactions through active particle dynamics

Adsorbate molecules present in a reaction mixture may bind to and block catalytic sites. Measurement of the surface coverage of these molecules via adsorption isotherms is critical for modeling and design of catalytic reactions on surfaces. However, it is challenging to measure isotherms in solution in a way that is directly relevant to catalytic activity under reaction conditions, particularly since adsorbates may bind with an enormous range of surface affinity parameters. Here we used the motion of self-propelled catalytic Janus particles, which employ the decomposition of hydrogen peroxide fuel as a propulsion mechanism, to determine the effective surface coverage of thioglycerol, furfural, and ethanol on a platinum surface as a function of concentration in aqueous solution by measuring the decrease in active motion due to the blocking of active sites. For strongly adsorbing thioglycerol, this effective coverage was compared and contrasted to the total adsorbed amount measured using inductively-coupled plasma analysis. Demonstrating the broad applicability of this approach, the surface affinity of the three adsorbates spanned more than four orders of magnitude. For each species, the adsorbate-mediated attenuation of active motion occurred over a wide concentration range and was well-described by a Langmuir isotherm. The strongly interacting thioglycerol had the highest affinity towards the surface (K a = 15.5 ± 4.3 mM –1 ) and fully deactivated the active particle motion at surface saturation. Furfural had an intermediate affinity (K a = 0.42 ± 0.07 mM –1 ) but did not fully block H 2 O 2 access to the surface at apparent saturation, consistent with a maximum fractional surface coverage of θ max = 0.67. Ethanol exhibited even lower affinity (K a = 0.0025 ± 2x10 -4 mM –1 ) and its coverage saturated at only θ max = 0.38. Analysis of isotherms at elevated temperatures enabled direct extraction of the enthalpies of adsorption. The degree of surface coverage at adsorbate saturation appeared to correlate with the relative energies of adsorption for the different adsorbate species and was consistent with adsorbate saturation of one of multiple active site populations towards H 2 O 2 decomposition. Furthermore, computational investigations into solvent effects on furfural adsorption showed good quantitative agreement with the experimental results. This work leverages unique properties of active particles to explore fundamental catalysis questions and demonstrates a novel paradigm for significant and experimentally accessible multidisciplinary research.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Unraveling Electronic Trends in O* and OH* Surface Adsorption in the MO 2 Transition-Metal Oxide Series

Understanding the bond strength of O* and OH* intermediates to metal-oxide surfaces is key to predicting the catalytic activity in oxygen-based electrochemistry. Here, we uncover highly non-linear trends in O* and OH* adsorption energies across the 3d, 4d, and 5d series of MO 2 transition-metal (TM) oxide surfaces computed within Hubbard- U corrected density functional theory (DFT + U ). Investigating the electronic structure with crystal orbital Hamiltonian populations (COHP) of the relevant metal–oxygen bonds reveals that the spin-dependent coupling strength between metal-d and oxygen-2p atomic orbitals together with the extent of filling of bonding and anti-bonding orbitals are the primary contributors to the adsorption energy. Importantly, we show that the integrated COHP obtained purely from bulk calculations is a highly accurate descriptor for surface adsorption energetics that captures trends across the group 5–12 TM oxide series within 0.19–0.36 eV. Our results suggest a pathway to prediction of adsorption energies for an arbitrary metal–ligand catalyst system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Meta-Generalized Gradient Approximation for the Cavity-Dependent Exchange-Correlation Interaction in Strongly Coupled Light–Matter Systems

Strong light–matter coupling in optical cavities enables the manipulation of chemical and physical properties without altering molecular composition. Theoretical modeling of such phenomena requires exchange-correlation (XC) functionals that account for both electron–electron and electron–photon (ep) interactions within quantum electrodynamical density functional theory (QEDFT). In this work, we develop a meta-generalized gradient approximation (meta-GGA) specifically targeting the cavity-dependent XC interaction in strongly coupled light–matter systems. This novel approximation is built upon a new semilocal polarizability approximation, which draws from the jellium-with-a-gap model, and can be extended to a “global hybrid” variant that goes beyond the isotropic model from previous approximations. The polarizability model yields significantly improved dispersion coefficients and benchmark calculations with the cavity-dependent XC functional demonstrate improved agreement with QED Hartree–Fock (QED-HF) reference energies. Application to the regioselectivity of brominated nitrobenzene intermediates reveals the functional’s capacity to capture cavity-induced energetic shifts. In conclusion, our results advance the Jacob’s ladder of functionals for QEDFT and provide a practical tool for modeling polaritonic chemistry.

Approximation↗

Direct Observation of Shock‐Induced Disordering of Enstatite Below the Melting Temperature

Abstract We report in situ structural measurements of shock‐compressed single crystal orthoenstatite up to 337 ± 55 GPa on the Hugoniot, obtained by coupling ultrafast X‐ray diffraction to laser‐driven shock compression. Shock compression induces a disordering of the crystalline structure evidenced by the appearance of a diffuse X‐ray diffraction signal at nanosecond timescales at 80 ± 13 GPa on the Hugoniot, well below the equilibrium melting pressure (>170 GPa). The formation of bridgmanite and post‐perovskite have been indirectly reported in microsecond‐scale plate‐impact experiments. Therefore, we interpret the high‐pressure disordered state we observed at nanosecond scale as an intermediate structure from which bridgmanite and post‐perovskite crystallize at longer timescales. This evidence of a disordered structure of MgSiO 3 on the Hugoniot indicates that the degree of polymerization of silicates is a key parameter to constrain the actual thermodynamics of shocks in natural environments.

Hernandez, J. ‐A.↗

Multiphoton ionization and multiphoton resonances in the tunneling regime

The rate of ionization of an atom of helium, argon, or hydrogen exposed to an intense monochromatic laser field and the quasienergy spectrum of their dressed states are studied for values of the Keldysh parameter between 1 and 0.6 and wavelengths between 390 and 1300 nm. The calculations are carried out within the non-Hermitian Floquet theory. Resonances with intermediate excited states significantly affect ionization from the dressed ground state at all the intensities and all the wavelengths considered. The dressed excited states responsible for these structures are large-{alpha}{sub 0} states akin to the Kramers-Henneberger states of the high-frequency Floquet theory. Within the single-active-electron approximation, these large-{alpha}{sub 0} states become species independent at sufficiently high intensity or sufficiently long wavelength. Apart for the resonance structures arising from multiphoton coupling with excited states, the ab initio Floquet ionization rate is in excellent agreement with the predictions of two different calculations in the strong field approximation, one based on a length-gauge formulation of this approximation and one based on a velocity-gauge formulation. The calculations also confirm the validity of the {omega}{sup 2} expansion as an alternative to the strong field approximation for taking into account the nonadiabaticity of the ionization process in intense low-frequency laser fields.

74 ATOMIC AND MOLECULAR PHYSICS↗

Real-Time Monitoring of Fracture Dynamics with a Contrast Agent-Assisted Electromagnetic Method

In collaboration with the Advanced Energy Consortium, our team has previously demonstrated that the placement of electrically active proppants (EAPs) in a hydraulic fracture surveyed by electromagnetic (EM) methods can enhance the imaging of the stimulated reservoir volumes during hydraulic fracturing. That work culminated in constructing a well-characterized EAP-filled fracture anomaly at the Devine field pilot site (DFPS). In subsequent laboratory studies, we observed that the electrical conductivity of our EAP correlates with changes in pressure, salinity, and flow. Thus, we postulated that the EAP could be used as an in-situ sensor for the remote monitoring of these changes in previously EAP-filled fractures. This paper presents our latest field data from the DFPS to demonstrate such correlations at an intermediate pilot scale. We conducted surface-based EM surveys during freshwater (200 ppm) and saltwater (2,500 ppm) slug injections while running surfaced-based EM surveys. Simultaneously, we measured the following: 1) bottomhole pressure and salinity in five monitoring wells; 2) injection rate using high-precision data loggers; 3) distributed acoustic sensors in four monitoring wells; and 4) tiltmeter data on the survey area. We demonstrated that injections into an EAP-filled fracture could be successfully coupled with real-time electric field measurements on the surface, leading to remote monitoring of dynamic changes within the EAP-filled fracture. Furthermore, by comparing the electrical field traces with the bottomhole pressure, flow rate, and salinity, we concluded that the observed electric field in our study is influenced by fracture dilation and flow rate. Salinity effect was observed when saltwater was injected. EM simulations solely based on assumptions of fracture conductivity changes during injection did not reproduce all of the measured electric field magnitudes. Preliminary estimates showed that including streaming potential in our geophysical model may be needed to reduce the simulation mismatch. The methods developed and demonstrated during this study will lead to a better understanding of the extent of fracture networks, formation stress states, fluid leakoff and invasion, characterizations of engineered fracture systems, and other applications where monitoring subsurface flow tracking is deemed important.

02 PETROLEUM↗

An infrared search for extraterrestrial laser signals

The focus of project SETI is on microwave frequencies, where receivers fundamentally have the best sensitivity for the detection of narrow band signals. Such receivers, when coupled to existing radio telescopes, form an optimum system for broad area searches over the sky. Detection of narrow band infrared signals is best done with a laser heterodyne reciever similar in function to a microwave spectral line receiver. A receiver was built for astrophysical observations at 30 THz (10 microns) and the spectrometer is being adapted for SETI work. The receiver uses a small CO2 laser as the local oscillator, a HgCdTe diode as the photomixer, and a multichannel intermediate frequency (IF) filterbank. An advanced multichannel IF processor is now being built to detect infrared line radiation in 1000 spectral channels each 1 MHz wide. When completed this processor will be used with a ground based telescope next year for a survey of several hundred selected stars for narrow band CO2 laser signals at 30 THz.

Betz, A.↗

Mesoporous Colloidal Superparticles of Platinum-Group Nanocrystals with Surfactant-Free Surfaces and Enhanced Heterogeneous Catalysis

Synthesis of colloidal superparticles (CSPs) of nanocrystals, a class of assembled nanocrystals in the form of colloidal particles, has been emerging as a new frontier in the field of nanotechnology because of their potential novel properties originated from coupling of individual nanocrystals in CSPs. Here, a facile approach is reported for the controlled synthesis of mesoporous CSPs made of various platinum-group nanocrystals that exhibit high colloidal stability and ligand-free surfaces to significantly benefit their applications in solution-phase heterogeneous catalysis. The synthesis relies on self-limiting growth of composite particles through coprecipitation of both Pt-group nanocrystals (or their precursor compounds) and silver halides on sacrificial substrates of colloidal silver particles. The intermediate silver halides in the composite particles play the critical role in limiting the continuous growth (and/or coalescence) of individual Pt-group nanocrystals and they can be selectively dissolved to create nanoscale pores in the resulting CSPs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Heavy photon search experiment at JLAB

The Heavy Photon Search (HPS) experiment in Hall-B at Jefferson Lab will search for new heavy vector boson(s), aka 'heavy photons', in the mass range of 20 MeV/c{sup 2} to 1000 MeV/c{sup 2} using the scattering of high energy, high intensity electron beams off a high Z target. The proposed measurements will cover the region of parameter space favored by the muon g-2 anomaly, and will explore a significant region of parameter space, not only at large couplings (α′/α > 10{sup −7}), but also in the regions of small couplings, down to α′/α∼10{sup −10}. The excellent vertexing capability of the Si-tracker uniquely enables HPS to cover the small coupling region. Also, HPS will search for heavy photons in an alternative to the e{sup +}e{sup −} decay mode, in the heavy photon's decay to μ{sup +}μ−.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Quantum Spin Dynamics Due to Strong Kitaev Interactions in the Triangular-Lattice Antiferromagnet CsCeSe 2

The extraordinary properties of the Kitaev model have motivated an intense search for new physics in materials that combine geometrical and bond frustration. In this Letter, we employ inelastic neutron scattering, spin wave theory, and exact diagonalization to study the spin dynamics in the perfect triangular-lattice antiferromagnet (TLAF) CsCeSe 2 . This material orders into a stripe phase, which is demonstrated to arise as a consequence of the off-diagonal bond-dependent terms in the spin Hamiltonian. By studying the spin dynamics at intermediate fields, we identify an interaction between the single-magnon state and the two-magnon continuum that causes decay of coherent magnon excitations, level repulsion, and transfer of spectral weight to the continuum that are controlled by the strength of the magnetic field. Our results provide a microscopic mechanism for the stabilization of the stripe phase in TLAF and show how complex many-body physics can be present in the spin dynamics in a magnet with strong Kitaev coupling even in an ordered ground state. Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Identification of the molecular-weight growth reaction network in counterflow flames of the C 3 H 4 isomers allene and propyne

The reaction networks responsible for aromatics formation in counterflow flames of the C 3 H 4 isomers allene and propyne are identified through a combined experimental and modeling study. Mole fraction profiles of near-atmospheric pressure (933 mbar) diffusion flames fueled by the C 3 H 4 isomers are analyzed by means of a newly assembled, chemically detailed kinetic mechanism. The experiment consists of a counterflow burner system that is coupled to a high-resolution time-of-flight molecular-beam mass spectrometer with single-photon ionization via synchrotron-generated vacuum-ultraviolet photons. Flame-sampled, mass-specific photoionization efficiency curves are used to identify the presence of aliphatically substituted aromatic species in addition to the commonly considered pericondensed ring structures. The new mechanism describes the formation and growth of aromatics through repetitive sequences of radical–radical and radical–molecule reactions that include C 1 6 intermediates. Higher concentrations of aromatic species are observed in the allene flame and the new mechanism captures the observed experimental trends very accurately. The results indicate the importance of the aliphatically substituted aromatics and of ring-enlargement reactions for the growth reactions. According to the model simulations, radical+radical recombination and PAH-radical+molecule reactions play an important role in PAH growth.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The monzonorite-anorthosite connection: The petrogenesis of terrestrial KREEP

There is a suite of rocks typically associated with Proterozoic massif anorthosites that bear some interesting similarities to lunar KREEP. In many cases these rocks are plutonic and have traditionally been referred to as the jotunite-mangerite-+/-charnockite-+/-syenite suite. However, in the Rogaland district of southwestern Norway, where they are referred to as 'monzonorites', these rocks are also present as fine-grained dikes and as the chill margin of a layered intrusion, and thus approximate magmatic liquid compositions are readily obtained by chemical analysis. Monzonorites are typically enriched in incompatible lithophile elements such as K (alkali feldspar is present), the rare earths (REE), and P. They have intermediate to low Mg', low-Ca pyroxene, and more evolved types have low Ti/Sm ratios. Much debate has developed over attempts to explain the link between monzonorites and massif anorthosites. One feature seems clear: monzonorites and associated anorthosites have different initial isotopic ratios, so a simple relation is not possible. However, there is apparently a continuum in major elements between the monzonorites and gabbros believed to represent magmas parental to the anorthosites. This continuum suggests a link via high pressure fractionation coupled with assimilation. Although more complicated, this scenario is similar to that evoked for the early Moon: following the formation of ferroan anorthosites, continued fractional crystallization of the residual liquids at the base of the crust led to the formation of KREEP. An attempt is made here to establish a link between monzonorites and high-Al gabbros which are nearly always found as ancillary intrusions associated with anorthosites, and which may record processes in lower crustal magma chambers.

Longhi, J.↗

Photoinduced Electron Transfer Informs on Pathway Coupling in Flavin-Based Electron Bifurcation

Flavin-based electron bifurcation (FBEB) is an enzymatic mechanism that generates extremely high-energy electrons to drive unfavorable chemical reactions. It is utilized by the NADHdependent ferredoxin:NADP + -oxidoreductase (Nfn) enzyme in hyperthermophile Pyrococcus furiosus to bifurcate electrons from NADPH into the coupled low-potential (endergonic) and highpotential (exergonic) pathways. This process enables P. furiosus to live in harsh and uninviting environments. Despite its biological importance, the mechanisms used by Nfn to facilitate exceptional directional control over short-lived, high-energy electrons and to prevent undesired transfer, particularly along the low-potential pathway, are still not well understood. To elucidate how the protein environment contributes to electronic control in the lowpotential pathway, new techniques must be utilized to probe these unstable intermediates. In this study, we have adapted lowtemperature photoexcitation combined with electron paramagnetic resonance (EPR) to accumulate the short-lived intermediate and place it in the context of the other cofactors involved in the low-potential pathway of Nfn. We observed coincident growth of both the radical intermediate and its nearby [4Fe-4S] cluster over 4.5 h of illumination with NADPH at cryogenic temperatures. The photogenerated paramagnetic species were stable in LN 2 storage indefinitely and recombined when warmed to higher temperatures. The results provide insights into the electron transfer steps and cofactor interactions along the low potential pathway, facilitating a more robust mechanistic understanding of the high-energy events of electron bifurcation. Furthermore, through comparison of cryogenic and room temperature experiments, a potential gating step involving the movement of key residues important for the reversibility of electron flow along this pathway is suggested.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

From amorphous to microcrystalline: Phase transition in rapid synthesis of hydrogenated silicon thin film in low frequency inductively coupled plasmas

Hydrogenated silicon (Si:H) thin films were fabricated on glass substrates by low frequency inductively coupled plasma-assisted chemical vapor deposition using a silane precursor with low hydrogen dilution at room temperature. The crystallinity and microstructure properties of the Si:H thin films deposited at different inductive radio-frequency (rf) power density were systematically studied by Raman spectroscopy, x-ray diffraction, and scanning electron microscopy. We found that at a low rf power density of 16.7 to 20.8 mW/cm{sup 3}, the structure of silicon thin films evolves from a completely amorphous phase to an intermediate phase containing both amorphous and microcrystalline silicon. As the power density is increased to a moderate value of 25 mW/cm{sup 3}, a highly crystallized (111)-preferred hydrogenated microcrystalline silicon ({mu}c-Si:H) film featuring a vertically aligned cone-shaped structure, is emerging. Both the crystallinity and deposition rate exhibit a monotonic increase with the increase in the rf power density, reaching a maximum value of 85% and 1.07 nm/s, respectively, at a power density of 41.7 mW/cm{sup 3}. Scanning electron microscopy reveals that continuous and dense {mu}c-Si:H films with grain size of tens to hundreds nanometers can be achieved deterministically without the formation of amorphous incubation layer, and this is of great importance for synthesis of multilayer structures in p-i-n solar cells. The formation mechanism of the {mu}c-Si:H films and the elimination of the amorphous incubation layer are explained in terms of the high electron density and the plasma-surface interactions.

36 MATERIALS SCIENCE↗

Symmetry and lattice mismatch induced strain accommodation near and away from correlated perovskite interfaces

Distinct MnO{sub 6} octahedral distortions near and away from the La{sub 0.67}Sr{sub 0.33}MnO{sub 3}/SrTiO{sub 3}(001) (LSMO/STO) interface are quantified using synchrotron x-ray diffraction and dynamical x-ray diffraction simulations. Three structural regions of stress accommodation throughout the film thickness were resolved: near the LSMO/STO interface, intermediate region farther from the interface, and the main layer away from the interface. The results show that within the first two unit cells stress is accommodated by the suppression of octahedral rotations in the film, leading to the expansion of the c-axis lattice parameter. Farther from the interface film structure acquires octahedral tilts similar to thicker perovskite films under tensile stress, leading to a reduced c-axis parameter. We demonstrate that these regions are related to two different strain coupling mechanisms: symmetry mismatch at the interface and lattice mismatch in the rest of the film. The findings suggest new routes for strain engineering in correlated perovskite heterostructures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

GPU acceleration of rank-reduced coupled-cluster singles and doubles

Here, we have developed a graphical processing unit (GPU) accelerated implementation of our recently introduced rank-reduced coupled-cluster singles and doubles (RR-CCSD) method. RR-CCSD introduces a low-rank approximation of the doubles amplitudes. This is combined with a low-rank approximation of the electron repulsion integrals via Cholesky decomposition. The result of these two low-rank approximations is the replacement of the usual fourth-order CCSD tensors with products of second- and third-order tensors. In our implementation, only a single fourth-order tensor must be constructed as an intermediate during the solution of the amplitude equations. Owing in large part to the compression of the doubles amplitudes, the GPU-accelerated implementation shows excellent parallel efficiency (95% on eight GPUs). Our implementation can solve the RR-CCSD equations for up to 400 electrons and 1550 basis functions—roughly 50% larger than the largest canonical CCSD computations that have been performed on any hardware. In addition to increased scalability, the RR-CCSD computations are faster than the corresponding CCSD computations for all but the smallest molecules. We test the accuracy of RR-CCSD for a variety of chemical systems including up to 1000 basis functions and determine that accuracy to better than 0.1% error in the correlation energy can be achieved with roughly 95% compression of the ov space for the largest systems considered. We also demonstrate that conformational energies can be predicted to be within 0.1 kcal mol -1 with efficient compression applied to the wavefunction. Finally, we find that low-rank approximations of the CCSD doubles amplitudes used in the similarity transformation of the Hamiltonian prior to a conventional equation-of-motion CCSD computation will not introduce significant errors (on the order of a few hundredths of an electronvolt) into the resulting excitation energies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗