Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular evolution”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Comparative modeling reveals the molecular determinants of aneuploidy fitness cost in a wild yeast model

Although implicated as deleterious in many organisms, aneuploidy can underlie rapid phenotypic evolution. However, aneuploidy will be maintained only if the benefit outweighs the cost, which remains incompletely understood. To quantify this cost and the molecular determinants behind it, we generated a panel of chromosome duplications in Saccharomyces cerevisiae and applied comparative modeling and molecular validation to understand aneuploidy toxicity. We show that 74%–94% of the variance in aneuploid strains’ growth rates is explained by the cumulative cost of genes on each chromosome, measured for single-gene duplications using a genomic library, along with the deleterious contribution of small nucleolar RNAs (snoRNAs) and beneficial effects of tRNAs. Machine learning to identify properties of detrimental gene duplicates provided no support for the balance hypothesis of aneuploidy toxicity and instead identified gene length as the best predictor of toxicity. Our results present a generalized framework for the cost of aneuploidy with implications for disease biology and evolution.

genic load↗
DOE OSTI2026

Transient Pulse-Response Time-of-Flight Mass Spectrometry for Complex, Deactivating Heterogeneous Catalytic Systems: Application to Ethane Dehydroaromatization

The study of complex, multistep bond-forming and -breaking reactions in heterogeneous catalytic systems often encounters challenges associated with the involvement of large numbers of intermediates among branching pathways. Kinetic information obtained from traditional steady-state measurements can be complemented with that from time-resolved methods to uncover details of the underlying chemistry. Herein, we describe an approach for tracking the complete time-resolved chemical composition (ca. 4–200 u) of a reactor effluent in response to a reactant pulse. We use a six-port rotary valve with a metered sampling loop to pulse reactants at ambient pressure into a flow reactor packed with a catalyst bed within the isothermal region of a heated furnace. The temporal evolution of effluent species is tracked using time-resolved molecular-beam time-of-flight mass spectrometry. We highlight the possibilities that this method has to offer by studying the complex bifunctional mechanism of ethane dehydroaromatization over an HZSM-5-supported platinum catalyst. We demonstrate that energy-tunable ionization sources, which facilitate isomer resolution, enable the measurement of the full mass spectral time-dependent system response. This includes the evolution of major products and mechanistically relevant reactive intermediates such as 1,3-butadiene and cyclopentadiene; these species have not previously been observed from this reaction. In additional studies, we also assess the role of platinum in the catalyst by examining temporal responses to ethane and ethylene feeds. Results show that two temporally distinct formation pathways exist for methane and benzene, and that their importance depends on both catalyst composition and reactant identity. Additionally, characteristics of catalyst deactivation are uniquely observable in the time-resolved mass spectral response, including the selective deactivation of a benzene formation pathway. The combination of time-of-flight mass spectrometry with tunable ionization enables the simultaneous observation of all effluents in a complex mixture of intermediates with isomer/isobar differentiation capabilities that can be applied to any complex reaction system.

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

Prediction of vacancy defect diffusion paths in high entropy alloys via machine learning on molecular dynamics data

Identifying the diffusion path of point defects is a critical step in understanding their evolution and the mechanisms of related phenomena. Defect diffusion occurs at small length and time scales, with impacts on material properties that may continue to evolve over ns to μs, ms, and the continuum scale (s, min, etc., and cm, m, etc.). The time scale accessible to molecular dynamics (MD) simulations is limited by small step sizes, typically in the fs range. Thus, surrogate models of MD simulations through machine learning (ML)-based algorithms are of great interest, especially for complex systems such as high entropy alloys (HEAs). In this work, dynamics governing vacancy migration in HEA were approximated with graph convolutional network (GCN) models as ansatzes for kinetic Monte Carlo (KMC) rate catalogs. Network design considered that diffusion in crystalline solids generally depends on interactions between defects and their immediate neighbor atoms. Graphs represented the vacancy surroundings, MD-generated trajectories provided training and comparison datasets, and unsupervised GCN models approximated interatomic dynamics governing vacancy migration in HEAs as ansatzes for KMC. A proof-of-concept model trained on MD data for the Fe, Ni, Cr, Co, and Cu HEA environment was used with two different neighbor interactions to assess the feasibility of training a GCN to predict vacancy defect transition rates in the HEA environment. The resulting setup rapidly generated MD-formatted synthetic trajectories based on dynamics learned from the MD training set, with a time acceleration of roughly two orders of magnitude and a similar diffusion coefficient to MD observations. Additionally, Nudged Elastic Band (NEB) calculations were performed on randomly generated FeNiCrCoCu HEA structures to determine vacancy migration barriers across nearest-neighbor sites. Transition probabilities for each jump, categorized by atomic type, were extracted from these calculations. NEB-based and GCN-based approaches led to similar outcomes.

Reimer, C↗
DOE OSTI2022

Analyzing the effect of pressure on the properties of point defects in $\gamma\mathrm{U–Mo}$ through atomistic simulations

Uranium–molybdenum (U–Mo) alloys in monolithic fuel foil are the primary candidate for the conversion of high-performance research reactors in the USA. Monolithic fuel is utilized in a plate-type design with a zirconium diffusion barrier and aluminum cladding. These fuel types are unique in that they contain no plenum for the release of fission gases, which, in conjunction with the aluminum cladding, can lead to large stress states within the fuel. The nature of how fundamental processes of radiation damage, including the evolution of point defects, under such stresses occur is unknown. In this work, we present molecular dynamics simulations of the formation energy of point defects under applied stress. We report this work will allow for the implementation of stress-dependent microstructural evolution models of nuclear fuels, including those for both fission gas bubble growth and creep, which are critical to ensure the stable and predictable behavior of research reactor fuels.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Understanding the Reaction of Nuclear Graphite with Molecular Oxygen

A thorough understanding of oxidation is important when considering the health and integrity of graphite components in graphite reactors. For the next generation of graphite reactors, HTGRs specifically, an unlikely air ingress has been deemed significant enough to have made its way into the licensing applications of many international licensing bodies. While a substantial body of literature exists on nuclear graphite oxidation in the presence of molecular oxygen and significant efforts have been made to characterize oxidation kinetics of various grades, the value of existing information is somewhat limited. Often, multiple competing processes, including reaction kinetics, mass transfer, and microstructural evolution, are lumped together into a single rate expression that limits the ability to translate this information to different conditions. This article reviews the reaction of graphite with molecular oxygen in terms of the reaction kinetics, gas transport, and microstructural evolution of graphite. It also presents the foundations of a model for the graphite-molecular oxygen reaction system that is kinetically independent of graphite grade, and is capable of describing both the bulk and local oxidation rates under a wide range of conditions applicable to air-ingress.

J. Kane, Joshua↗
DOE OSTI2026

pH-Driven Restructuring of Hydration Layers and Cation Ad-sorption at the Alumina-Water Interface

Oxide-water interfaces underpin ion separation, catalysis, and electrochemical energy technologies, where the electrical double layer (EDL) controls adsorption, transport, and reactivity. Yet, the molecular-scale link between pH-dependent surface protonation, hydration-layer structure, and counter-ion adsorption remains poorly defined. Here, we combine in situ crystal truncation rod (CTR) and resonant anomalous X-ray reflectivity (RAXR) with streaming potential measurements and ab initio molecular dynamics (AIMD) simulations to resolve the chemical and structural evolution of the EDL at the single-crystal alumina (012)-water interface in 10 mM Rb+ over pH 3-12. CTR measurements reveal two distinct adsorbed water layers at ~2.2 and ~3.5 Å above the surface that each shift toward the substrate at transition pHs near 6.5 and 10.6, respectively, directly reflecting changes in primary hydration layer structure in response to the deprotonation of bridging and terminal aluminol groups. RAXR shows a 10-fold increase in Rb+ coverage and a decrease in mean adsorption height from ~3.5 to ~2.7 Å with increasing pH, indicating enhanced counter-ion binding accompanied by Stern layer contraction. Streaming potential measurements demonstrate that the zeta potential, i.e., potential at the hydrodynamic shear plane, is positive at pH 3 and becomes negative at pH ≥3.5. This negative charge magnitude increases with increasing pH, consistent with progressive surface deprotonation at higher pH. AIMD identifies inner- and outer-sphere Rb+ complexes whose adsorption heights and coordination geometries depend sensitively on the protonation state of surface oxygens, providing atomistic support for the experimentally inferred trends. These measurements establish two discrete, site-specific pH transitions in hydration-layer structure that track aluminol (de)protonation and quantitatively link them to a pH-driven contraction of the Stern layer (increasing Rb+ coverage and decreasing adsorption height). This provides a direct structural basis for connecting surface acid-base chemistry to ion binding distances at an oxide-water interface.

Electrical double layer (EDL), Surface protonation↗
DOE OSTI2024

Femtosecond Core-Level Spectroscopy Reveals Involvement of Triplet States in the Gas-Phase Photodissociation of Fe(CO) 5

Excitation of iron pentacarbonyl [Fe(CO) 5 ], a prototypical photocatalyst, at 266 nm causes the sequential loss of two CO ligands in the gas phase, creating catalytically active, unsaturated iron carbonyls. Despite numerous studies, major aspects of its ultrafast photochemistry remain unresolved because the early excited-state dynamics have so far eluded spectroscopic observation. This has led to the long-held assumption that ultrafast dissociation of gas-phase Fe(CO) 5 proceeds exclusively on the singlet manifold. Herein, we present a combined experimental–theoretical study employing ultrafast extreme ultraviolet transient absorption spectroscopy near the Fe M 2,3 -edge, which features spectral evolution on 100 fs and 3 ps time scales, alongside high-level electronic structure theory, which enables characterization of the molecular geometries and electronic states involved in the ultrafast photodissociation of Fe(CO) 5 . We assign the 100 fs evolution to spectroscopic signatures associated with intertwined structural and electronic dynamics on the singlet metal-centered states during the first CO loss and the 3 ps evolution to the competing dissociation of Fe(CO) 4 along the lowest singlet and triplet surfaces to form Fe(CO) 3 . Calculations of transient spectra in both singlet and triplet states as well as spin–orbit coupling constants along key structural pathways provide evidence for intersystem crossing to the triplet ground state of Fe(CO) 4 . Finally, our work presents the first spectroscopic detection of transient excited states during ultrafast photodissociation of gas-phase Fe(CO) 5 and challenges the long-standing assumption that triplet states do not play a role in the ultrafast dynamics.

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

Data Generation for Machine Learning Interatomic Potentials and Beyond

The field of data-driven chemistry is undergoing an evolution, driven by innovations in machine learning models for predicting molecular properties and behavior. Recent strides in ML-based interatomic potentials have paved the way for accurate modeling of diverse chemical and structural properties at the atomic level. The key determinant defining MLIP reliability remains the quality of the training data. A paramount challenge lies in constructing training sets that capture specific domains in the vast chemical and structural space. This Review navigates the intricate landscape of essential components and integrity of training data that ensure the extensibility and transferability of the resulting models. We delve into the details of active learning, discussing its various facets and implementations. We outline different types of uncertainty quantification applied to atomistic data acquisition and the correlations between estimated uncertainty and true error. The role of atomistic data samplers in generating diverse and informative structures is highlighted. Furthermore, we discuss data acquisition via modified and surrogate potential energy surfaces as an innovative approach to diversify training data. The Review also provides a list of publicly available data sets that cover essential domains of chemical space.

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

Small-Molecule Models of Hydrogen-Evolving MX 2 (M = Mo, W; X = S, Se) Bulk Solids: Composition–Activity Relationships

Triangular metal chalcogenide clusters of the form [M 3 Q 7 L 3 ]An (M = Mo or W; Q = S or Se; L = i Bu 2 NCS 2 – , (CF 3 CH 2 ) 2 NCS 2 – , i Bu 2 NCSe 2 – , or i Bu 2 PS 2 – ; An = Cl – or I – ) have been investigated as molecular analogues of layered metal dichalcogenide (MX 2 ) H 2 -evolution catalysts. These clusters have been evaluated for their relative H 2 -evolving ability under a common photolysis protocol implementing [Ru(bpy) 3 ] 2+ as chromophore and Et 3 N as sacrificial electron donor. With M constant as Mo and with constant supporting ligand, clusters with an all-sulfide core enable greater H 2 -TON than clusters with an all-selenide core. A more active catalyst is produced by [Mo 3 S 7 (S 2 CN i Bu 2 ) 3 ] + I – than its W 3 analogue with the same core sulfide composition and supporting dithiocarbamate ligands. Dichalcogenocarbamate ligands provide more active catalysts than dialkyldithiophosphate ligated clusters, and within the dichalcogenocarbamate set, greater H 2 -turnovers correlate with more-electron-donating ligands (i.e., i Bu 2 NCS 2 – > (CF 3 CH 2 ) 2 NCS 2 – > i Bu 2 NCSe 2 – ). Cluster cations with Cl – as counteranion are very similar in activity H 2 -evolving levels to identical clusters with I – , ruling out any significant interfering effect by I – upon the electron transfer relay between Et 3 N and catalyst. In the aggregate, observations are consistent with a mechanism for H 2 evolution that involves reductive extrusion of H 2 from a metal hydride intermediate.

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

Sculpting 2D Crystals via Membrane Contractions before and during Solidification

When phospholipids crystallize within the otherwise fluid membranes of giant unilamellar vesicles, the resulting molecularly thin “2D” solids exhibit great variety in their morphology evolution. For example, within membranes containing moderate amounts of the crystallizing component, crystals grow with a fixed morphology depending on vesicle size. Conversely for membranes containing large amounts of the crystallizing species, we find small compact crystals on vesicles of all sizes. However, on large vesicles, growing crystals sprout flower petals that lengthen progressively. These behaviors result from two combined mechanisms: first, like other 2D solids, the shear rigidity of phospholipid crystals renders them intolerant to morphologies with nonzero Gaussian curvature. As a result and especially at elevated membrane tension, the cost of bending elasticity is reduced at the expense of line energy by the formation of flowers as opposed to compact crystals. Second, the composition-dependent tension rise during cooling relaxes via water permeation of the membrane with a time constant scaling as R2. The amount of crystal formed for a small decrease in temperature determines this composition-dependent increase in stress from thermal contractions versus solidification. Surface Evolver computations were motivated using the predicted tension evolution to develop a processing space that maps to experimental observations for initial and growing crystal morphology. Important variable groups are identified, including a scaled ratio of bending to line energy, a vesicle-size-independent group for membrane contractions, and a time constant for stress relaxation. Though processing stresses ultimately relax, the crystal morphology persists well beyond the processing window.

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

Vibrational Dynamics and Phase Transitions of Hydrazine to 50 GPa

The high pressure behavior of hydrazine, N 2 H 4 , has been investigated to 50 GPa at room temperature using infrared and Raman spectroscopy to explore pressure induced phase transitions and changes in hydrogen bonding. Three solid–solid phase transitions were detected at 11, 21, and 32 GPa on room temperature compression through dramatic changes in the lattice vibration and N–H stretching regions with increasing pressure in both measurement techniques. The transition to phase IV, which appears at 32 GPa, exhibits increased hydrogen bonding with significant hysteresis, persisting to 9 GPa on decompression. This work presents a detailed analysis of the pressure dependence of mode shifts and calculations of mode Grüneisen parameters as well as a determination of an approximate thermodynamic Grüneisen parameter. We compare these results to the behavior of other small molecular materials such as ammonia and water and explore the evolution of hydrogen bonding in hydrazine toward the symmetrically hydrogen bonded state, which has previously been suggested by theoretical computations.

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

pH-Driven Restructuring of Hydration Layers and Cation Adsorption at the Alumina–Water Interface

Oxide−water interfaces underpin ion separations, catalysis, and electrochemical energy technologies, where the electrical double layer (EDL) controls adsorption, transport, and reactivity. However, the molecular-scale links between pHdependent surface protonation, hydration-layer structure, and counterion adsorption remain poorly defined. Here, we combine in situ crystal truncation rod (CTR) and resonant anomalous X-ray reflectivity (RAXR) with streaming potential measurements and ab initio molecular dynamics (AIMD) simulations to resolve the chemical and structural evolution of the EDL at the single-crystal alumina (012)−water interface in 10 mM Rb + over pH 3−12. CTR measurements reveal two distinct adsorbed water layers at ∼2.2 and ∼3.5 Å above the surface. Each water layer shifts toward the substrate at transition pHs near 6.5 and 10.6, respectively, reflecting changes in primary hydration layer structure in response to the deprotonation of bridging and terminal aluminol groups. RAXR shows a 10-fold increase in Rb + coverage and a decrease in mean adsorption height from ∼3.5 to ∼2.7 Å with increasing pH, indicating enhanced counterion binding accompanied by Stern layer contraction. Streaming potential measurements demonstrate that the zeta potential, i.e., the potential at the hydrodynamic shear plane, is positive at pH 3 and becomes negative at pH ≥ 3.5. The negative charge magnitude increases with increasing pH, consistent with progressive surface deprotonation at higher pH. AIMD identifies inner- and outer-sphere Rb+ complexes whose adsorption heights and coordination geometries depend sensitively on the protonation state of surface oxygens, providing atomistic support for the experimentally inferred trends. These measurements establish two discrete, site-specific pH transitions in hydration-layer structure that track aluminol (de)protonation and quantitatively link them to a pH-driven contraction of the Stern layer (increasing Rb + coverage and decreasing adsorption height). This provides a direct structural basis for connecting surface acid−base chemistry to ion binding distances at an oxide−water interface.

Adsorption↗
DOE OSTI2023

Influence of near-surface oxide layers on TiFe hydrogenation: mechanistic insights and implications for hydrogen storage applications

The inevitable formation of passivating oxide films on the surface of the TiFe intermetallic compound limits its performance as a stationary hydrogen storage material. Extensive experimental efforts have been dedicated to the activation of TiFe, i.e. oxide layer removal prior to utilization for hydrogen storage. However, development of an efficient activation protocol necessitates a fundamental understanding of the composition and structure of the air-exposed surface and its interaction with hydrogen, which is currently absent. Therefore, in this study we explored the growth and nature of oxide films on the most exposed TiFe surface (110) in depth using static and dynamic first-principles methods. We identified the lowest energy structures for six oxygen coverages up to approximately 1.12 nm of thickness with a global optimization method and studied the temperature effects and structural evolution of the oxide phases in detail via ab initio molecular dynamics (AIMD). Based on structural similarity and coordination analysis, motifs for TiO 2 and TiFeO 3 as well as Ti(FeO 2 ) x (x = 2, 3 or 5) phases were identified. On evaluating the interaction of the oxidized surface with hydrogen, a minimal energy barrier of 0.172 eV was predicted for H 2 dissociation while H migration from the top of the oxidized surface to the bulk TiFe was limited by several high-lying energy barriers above 1.4 eV. Our mechanistic insights will prove themselves valuable for informed designs towards new activation methods of TiFe and related systems as hydrogen storage materials.

08 HYDROGEN↗
DOE OSTI2020

Genome resequencing reveals demographic history and genetic architecture of seed salinity tolerance in Populus euphratica

Populus euphratica is a dominant tree species in desert riparian forests and possesses extraordinary adaptation to salinity stress. Exploration of its genomic variation and molecular underpinning of salinity tolerance is important for elucidating population evolution and identifying stress-related genes. Here, we identify approximately 3.15 million single nucleotide polymorphisms using whole-genome resequencing. The natural populations of P. euphratica in northwest China are divided into four distinct clades that exhibit strong geographical distribution patterns. Pleistocene climatic fluctuations and tectonic deformation jointly shaped the extant genetic patterns. A seed germination rate-based salinity tolerance index was used to evaluate seed salinity tolerance of P. euphratica and a genome-wide association study was implemented. A total of 38 single nucleotide polymorphisms were associated with seed salinity tolerance and were located within or near 82 genes. Expression profiles showed that most of these genes were regulated under salt stress, revealing the genetic complexity of seed salinity tolerance. Furthermore, DEAD-box ATP-dependent RNA helicase 57 and one undescribed gene (CCG029559) were demonstrated to improve the seed salinity tolerance in transgenic Arabidopsis. These results provide new insights into the demographic history and genetic architecture of seed salinity tolerance in desert poplar.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Establishing defect-property relationships for 2D-nanomaterials (Final technical report)

Studies of new families of two-dimensional nanomaterials (2DNMs) have established that they possess unique properties that diverge from those of their bulk counterparts. This project aimed to determine how tolerant 2DNMs are to extreme photon and particle fluxes and to identify the mechanisms governing their radiation response. Early work indicated that graphene is not as representative of other 2DNMs as previously assumed; however, the origin of this difference remained unclear. To address these knowledge gaps, this project investigated the physical processes occurring at multiple length scales in transition metal dichalcogenides (TMDs) and quantified their structural stability and property evolution under far-from-equilibrium conditions. In-situ and ex-situ ion and electron irradiations were performed to directly control and monitor defect formation in TMDs. Complementary density functional theory (DFT) and molecular dynamics (MD) simulations were employed to elucidate the mechanisms of defect generation and evolution. The outcomes of this work established mechanistic understanding of irradiation-induced defect formation in 2DNMs, quantified the radiation tolerance of TMDs, and built a fundamental knowledge base for correlating defect structures with material properties. Collectively, these results provide new insights into the stability of low-dimensional materials in extreme environments and enable the predictive design of radiation-tolerant 2DNMs.

2D materials↗
DOE OSTI2022

Evolution of microscopic heterogeneity and dynamics in choline chloride-based deep eutectic solvents

Deep eutectic solvents (DESs) are an emerging class of non-aqueous solvents that are potentially scalable, easy to prepare and functionalize for many applications ranging from biomass processing to energy storage technologies. Predictive understanding of the fundamental correlations between local structure and macroscopic properties is needed to exploit the large design space and tunability of DESs for specific applications. Here, we employ a range of computational and experimental techniques that span length-scales from molecular to macroscopic and timescales from picoseconds to seconds to study the evolution of structure and dynamics in model DESs, namely Glyceline and Ethaline, starting from the parent compounds. We show that systematic addition of choline chloride leads to microscopic heterogeneities that alter the primary structural relaxation in glycerol and ethylene glycol and result in new dynamic modes that are strongly correlated to the macroscopic properties of the DES formed.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Discontinuous Galerkin discretization for quantum simulation of chemistry

All-electron electronic structure methods based on the linear combination of atomic orbitals method with Gaussian basis set discretization offer a well established, compact representation that forms much of the foundation of modern correlated quantum chemistry calculations—on both classical and quantum computers. Despite their ability to describe essential physics with relatively few basis functions, these representations can suffer from a quartic growth of the number of integrals. Recent results have shown that, for some quantum and classical algorithms, moving to representations with diagonal two-body operators can result in dramatically lower asymptotic costs, even if the number of functions required increases significantly. We introduce a way to interpolate between the two regimes in a systematic and controllable manner, such that the number of functions is minimized while maintaining a block-diagonal structure of the two-body operator and desirable properties of an original, primitive basis. Techniques are analyzed for leveraging the structure of this new representation on quantum computers. Empirical results for hydrogen chains suggest a scaling improvement from O(N 4.5 ) in molecular orbital representations to O(N 2.6 ) in our representation for quantum evolution in a fault-tolerant setting, and exhibit a constant factor crossover at 15 to 20 atoms. Moreover, we test these methods using modern density matrix renormalization group methods classically, and achieve excellent accuracy with respect to the complete basis set limit with a speedup of 1–2 orders of magnitude with respect to using the primitive or Gaussian basis sets alone. These results suggest our representation provides significant cost reductions while maintaining accuracy relative to molecular orbital or strictly diagonal approaches for modest-sized systems in both classical and quantum computation for correlated systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗