Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ADIABATIC PROCESSES”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Microstructural characteristics of adiabatic shear localization in a metastable beta titanium alloy deformed at high strain rate and elevated temperatures

The microstructural evolution and grain refinement within adiabatic shear bands in the Ti6554 alloy deformed at high strain rates and elevated temperatures have been characterized using transmission electron microscopy. No stress drops were observed in the corresponding stress–strain curve, indicating that the initiation of adiabatic shear bands does not lead to the loss of load capacity for the Ti6554 alloy. The outer region of the shear bands mainly consists of cell structures bounded by dislocation clusters. Equiaxed subgrains in the core area of the shear band can be evolved from the subdivision of cell structures or reconstruction and transverse segmentation of dislocation clusters. It is proposed that dislocation activity dominates the grain refinement process. The rotational recrystallization mechanism may operate as the kinetic requirements for it are fulfilled. The coexistence of different substructures across the shear bands implies that the microstructural evolution inside the shear bands is not homogeneous and different grain refinement mechanisms may operate simultaneously to refine the structure. - Graphical abstract: Display Omitted - Highlights: • The microstructure within the adiabatic shear band was characterized by TEM. • No stress drops were observed in the corresponding stress–strain curve. • Dislocation activity dominated the grain refinement process. • The kinetic requirements for rotational recrystallization mechanism were fulfilled. • Different grain refinement mechanisms operated simultaneously to refine the structure.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Shear localization in metallic materials at high strain rates

Three factors govern adiabatic shear localization: strain hardening (or softening), strain-rate hardening, and thermal softening. It is typically associated with large shear strains (>1), high strain rates (10 3 –10 7 s –1 ), and high temperatures (40–100% of melting point), all of which happen within narrow regions with widths of about 1–200 μm. It is often an undesirable phenomenon, leading to failure, but there are situations where it is desirable, e. g., the generation of machining chips. Here, we review the development of both theoretical and experimental achievements, from the initiation of shear bands to their propagation with emphasis on three aspects: novel experimental techniques, novel materials, and nano/microstructural effects. The principal characteristics of adiabatic shear bands in metallic materials at the nano- and micro-scale are described. Bands that were formerly identified as transformed actually consist of nanocrystalline/ultrafine grains. These grains result from the breakup of the microstructure by a rotational recrystallization process. The evolution of the microstructure inside shear bands and their interactions for hcp, bcc, and fcc alloys, high-entropy alloys, nanocrystalline alloys, and metallic glasses are analyzed mechanistically. In conclusion, the gaps in the field and opportunities for future research are identified. Modern experimental characterization and computational techniques enable a more profound and predictive understanding of adiabatic shear localization and its avoidance in advanced materials.

36 MATERIALS SCIENCE↗
DOE OSTI2013

On the structure of plasma liners for plasma jet induced magnetoinertial fusion

The internal structure and self-collapse properties of plasma liners, formed by the merger of argon plasma jets, have been studied via 3-dimensional numerical simulations using the FronTier code. We have shown that the jets merger process is accomplished through a cascade of oblique shock waves that heat the liner and reduce its Mach number. Oblique shock waves and the adiabatic compression heating have led to the 10 times reduction of the self-collapse pressure of a 3-dimensional argon liner compared to a spherically symmetric liner with the same pressure and density profiles at the merging radius. We have also observed a factor of 10 variations of pressure and density in the leading edge of the liner along spherical surfaces close to the interaction with potential plasma targets. Such a non-uniformity of imploding plasma liners presents problems for the stability of targets during compression.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Fitting coupled potential energy surfaces for large systems: Method and construction of a 3-state representation for phenol photodissociation in the full 33 internal degrees of freedom using multireference configuration interaction determined data

A recently reported algorithm for representing adiabatic states coupled by conical intersections using a quasi-diabatic state Hamiltonian in four and five atom systems is extended to treat nonadiabatic processes in considerably larger molecules. The method treats all internal degrees of freedom and uses electronic structure data from ab initio multireference configuration interaction wave functions with nuclear configuration selection based on quasi-classical surface hopping trajectories. The method is shown here to be able to treat ∼30 internal degrees of freedom including dissociative and large amplitude internal motion. Two procedures are introduced which are essential to the algorithm, a null space projector which removes basis functions from the fitting process until they are needed and a partial diagonalization technique which allows for automated, but accurate, treatment of the vicinity of extended seams of conical intersections of two or more states. These procedures are described in detail. The method is illustrated using the photodissociaton of phenol, C{sub 6}H{sub 5}OH(X{sup ~1}A{sup ′}) + hv → C{sub 6}H{sub 5}OH(A{sup ~1}A{sup ′}, B{sup ~1}A{sup ′′}) → C{sub 6}H{sub 5}O(X{sup ~2}B{sub 1}, A{sup ~2}B{sub 2}) + H as a test case. Ab initio electronic structure data for the 1,2,3{sup 1}A states of phenol, which are coupled by conical intersections, are obtained from multireference first order configuration interaction wave functions. The design of bases to simultaneously treat large amplitude motion and dissociation is described, as is the ability of the fitting procedure to smooth the irregularities in the electronic energies attributable to the orbital changes that are inherent to nonadiabatic processes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Investigation of plasmon relaxation mechanisms using nonadiabatic molecular dynamics

Hot carriers generated from the decay of plasmon excitation can be harvested to drive a wide range of physical or chemical processes. However, their generation efficiency is limited by the concomitant phonon-induced relaxation processes by which the energy in excited carriers is transformed into heat. However, simulations of dynamics of nanoscale clusters are challenging due to the computational complexity involved. Here, in this paper, we adopt our newly developed Trajectory Surface Hopping (TSH) nonadiabatic molecular dynamics algorithm to simulate plasmon relaxation in Au 20 clusters, taking the atomistic details into account. The electronic properties are treated within the Linear Response Time-Dependent Tight-binding Density Functional Theory (LR-TDDFTB) framework. The relaxation of plasmon due to coupling to phonon modes in Au 20 beyond the Born–Oppenheimer approximation is described by the TSH algorithm. The numerically efficient LR-TDDFTB method allows us to address a dense manifold of excited states to ensure the inclusion of plasmon excitation. Starting from the photoexcited plasmon states in Au 20 cluster, we find that the time constant for relaxation from plasmon excited states to the lowest excited states is about 2.7 ps, mainly resulting from a stepwise decay process caused by low-frequency phonons of the Au 20 cluster. Furthermore, our simulations show that the lifetime of the phonon-induced plasmon dephasing process is ~10.4 fs and that such a swift process can be attributed to the strong nonadiabatic effect in small clusters. Our simulations demonstrate a detailed description of the dynamic processes in nanoclusters, including plasmon excitation, hot carrier generation from plasmon excitation dephasing, and the subsequent phonon-induced relaxation process.

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

Theoretical Models of Optical Transients. I. A Broad Exploration of the Duration–Luminosity Phase Space

The duration–luminosity phase space (DLPS) of optical transients is used, mostly heuristically, to compare various classes of transient events, to explore the origin of new transients, and to influence optical survey observing strategies. For example, several observational searches have been guided by intriguing voids and gaps in this phase space. However, we should ask, do we expect to find transients in these voids given our understanding of the various heating sources operating in astrophysical transients? In this work, we explore a broad range of theoretical models and empirical relations to generate optical light curves and to populate the DLPS. We explore transients powered by adiabatic expansion, radioactive decay, magnetar spin-down, and circumstellar interaction. For each heating source, we provide a concise summary of the basic physical processes, a physically motivated choice of model parameter ranges, an overall summary of the resulting light curves and their occupied range in the DLPS, and how the various model input parameters affect the light curves. We specifically explore the key voids discussed in the literature: the intermediate-luminosity gap between classical novae and supernovae, and short-duration transients (≲10 days). We find that few physical models lead to transients that occupy these voids. Moreover, we find that only relativistic expansion can produce fast and luminous transients, while for all other heating sources events with durations ≲10 days are dim (M{sub R}≳−15 mag). Finally, we explore the detection potential of optical surveys (e.g., Large Synoptic Survey Telescope) in the DLPS and quantify the notion that short-duration and dim transients are exponentially more difficult to discover in untargeted surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1973

Condensation in CO2 free jet expansions. I - Dimer formation.

A theoretical and experimental investigation of the kinetics of CO2 dimer formation in the adiabatically expanding core of a free jet expansion has been carried out. The theoretical analysis is based on dimer formation by an efficient termolecular process and dimer destruction by bimolecular collisions. Account is taken of the expected variation of the specific heat ratio during the expansion and of the effective increase in the termolecular collision rate at low temperatures due to the presence of loosely bound orbiting pairs. Experimental dimer concentrations were measured mass spectrometrically by forming a molecular beam from the expanded gas mixture after the onset of molecular flow. Observed dimer concentrations, for a range of pre-expansion conditions of 400 to 900 torr pressure and room temperature, are in good agreement with the theoretical analysis, provided that account is taken of the effect of loosely bound orbiting pairs at low temperatures.

Dorfeld, W. G.↗
NASA NTRS1982

Structure of the magnetopause rotational discontinuity

A model is developed for the rotational discontinuities associated with the reconnection configuration on the dayside magnetopause and the energy transfer process in the tail magnetopause. The model, in which the ion dynamics is described by the fluid equations and electrons are assumed to move adiabatically along the magnetic field lines, is valid for a rotational discontinuity with a thickness greater than a few ion gyroradii. It is shown, by the inclusion of self-consistency for trapped electrons, that (1) the trapped electron density profile is uniquely related to the rate of angular rotation of the magnetic field, and (2) the sense of magnetic field rotation is determined by the ratio of the normal and tangential components of the magnetic field. It is found that the electron polarization of magnetopause rotational discontinuities should be in agreement with satellite observations.

Lee, L. C.↗
NASA NTRS2019

High Temperature Superconductor Lead Assemblies for XRISM

The RESOLVE instrument on the upcoming XRISM mission is a rebuild of the SXS instrument that successfully demonstrated high resolution X-ray spectroscopy on the Hitomi mission using a microcalorimeter array at 50 mK. As with SXS, the RESOLVE instrument requires High Temperature Superconductor (HTS) lead assemblies to carry high current to the Adiabatic Demagnetization Refrigerator (ADR) while meeting a very stringent heat conductance budget. For the RESOLVE HTS lead assemblies, we have used new materials and processes to improve the stability and consistency of solder joints to the HTS material. We have also designed a new thermal strap that eliminates stress concentrations on the HTS tapes. We present the results of testing on the engineering model and flight model units, as well as the results of long term testing on prototype solder joints.

Canavan, Edgar R.↗
DOE OSTI2019

Constitutive Behavior and Hot Workability of a Hot Isostatic Pressed Ti-22Al-25Nb Alloy during Hot Compression

A Ti-22Al-25Nb alloy was fabricated from prealloyed powders by hot isostatic pressing for 2 h at a temperature of 1050 °C and pressure of 100 MPa. The hot deformation behavior of the Ti-22Al-25Nb alloy was characterized by isothermal compression testing at deformation temperatures between 900 and 1060 °C and strain rates between 0.001 and 1 s{sup −1}. Based on the true stress–strain curves, a constitutive equation was constructed to describe the flow stress as a function of the strain rate and deformation temperature. Three-dimensional (3D) processing maps were developed based on dynamic material model theory using the stress flow data to identify the instability and optimization regions of the hot processing parameters. The results show that the flow stress decreases with increasing deformation temperature and decreasing strain rate, and the softening mechanisms are different under different deformation conditions. The apparent activation energies in the (α{sub 2} + β/B2 + O) and (α{sub 2} + B2) phase regions are calculated to be 865.177 and 590.661 kJ/mol, respectively, which suggests that the hot isostatic pressed (HIPed) Ti-22Al-25Nb alloy requires a high hot deformation activation energy. Combined with microstructural observations, the optimal processing domains are determined to be a temperature range of 1030-1060 °C and strain rate range of 0.01-0.1 s{sup −1} in the α{sub 2} + B2 phase region. Moreover, the results indicate that adiabatic shear bands and severe inhomogeneous deformation cause flow instability at lower temperatures (900-1005 °C) and higher strain rates (> 0.1 s{sup −1}).

36 MATERIALS SCIENCE↗
DOE OSTI2019

The effects of ultra-fine-grained structure and cryogenic temperature on adiabatic shear localization in titanium

The deformation at low temperatures (173 K and 77 K) in ultrafine-grained (100 and 500 nm) titanium is investigated and its effect on adiabatic shear localization is established. In comparison with coarse grained titanium, the strength of ultrafine-grained titanium is higher due to the classic Hall-Petch effect while the strain hardening approaches zero. Our results show that shear localization in dynamic deformation is also altered. The width of the shear band of coarse-grained titanium decreases from 30 to 18 mm (by 40%) with decreasing the initial deformation temperature to 77 K. In contrast, for 100 nm titanium, the width of shear band decreases more significantly, from 4 mmat room temperature to 1 mm (a 75% decrease) at 77 K. This difference is attributed to the combined effects of the decrease in the thermal conductivity and specific heat capacity, and the increase in thermal softening rate. These changes in the width are consistent with the predictions of the Grady and Bai-Dodd theories. Ultrafine- and nanorecrystallized grains are observed inside the bands which are dependent on initial grain size and initial deformation temperature. Here, the dislocation evolution is evaluated for the different conditions using a Kocks-Mecking formulation; the rotational dynamic recrystallization mechanism responsible for forming the ultrafine/nanosized grains (40-250 nm) is successfully modeled incorporating the differences in initial temperature and grain size. Our results and analysis are important in enhancing the understanding of the structural evolution processes under high strain-rates and cryogenic temperatures.

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

Single-Bubble Flow Boiling Phenomena - Interface Tracking Simulation

Boiling, as one of the most efficient heat transfer mechanisms, is widely used in various engineering systems. Better understanding and modeling of this process remains a major challenge in multiphase flow research. The distribution of vapor in boiling system affects the heat transfer rate and may cause unfavorable conditions, such as heater burn-out. A number of flow regimes have been observed experimentally. In general, there are three regimes under boiling conditions: nucleate boiling, transition boiling, and film boiling. The nucleate boiling regime is categorized into partial boiling regime and fully developed nucleate boiling, according to behavior of bubble dynamics and heat transfer characteristics. In the first region, as the wall temperature increases, the fraction of the wall surface subject to nucleate boiling increases until bubble formation occupies the entire heated surface. Heat transfer in this region is a complex mixture of single-phase forced convection and nucleate boiling. As bubble density increases, bubbles turn to coalesce and form insulating vapor patches, which impedes liquid flow back to the surface. The heat transfer rate under this condition reduces dramatically. This point is called Critical Heat Flux (CHF), which leads to boiling crisis. The behavior of boiling crisis is dependent on local fluid conditions due to bubble growth and departure. Experiments are the ultimate tool for boiling model development and verification, but they still lack some capabilities needed for better understanding of the local mechanism of this phenomena. It's difficult to capture many important details at the pressure/temperature conditions of interest using experimental techniques, such as full three dimensional measurements of velocity and temperature distribution and estimates of heat flux partitioning during the boiling process. The interface tracking simulation (ITS) is one of the promising approaches to describe heat transfer of boiling phenomena and their underlying mechanisms. The advances of high performance computing (HPC) in recent years made it possible to apply DNS to a wide variety of adiabatic flows. Some work has been done to implement evaporation / condensation models as well. Lee and Nydahl performed a numerical simulation of bubble growth and departure assuming that the bubble has the shape of a truncated sphere and a micro-layer exists during bubble growth. A 2-D cylindrical bubbles simulation has been done by Son with a simplified model to predict the heat flux in a thin liquid micro-layer. A direct simulation of 2D film boiling is conducted by Juric and Tryggvason. The film boiling, below critical pressure with a modified Level Set Method, is studied by Son and Dhir. Ose and Kunugi have developed a boiling and condensation model on subcooled boiling phenomena. Fuchs performed numerical simulations of bubble growth and departure considering wall conduction effects. Fully transient heat and fluid flow is modeled with a free surface of a rising bubble and a periodic calculation of the whole cycling of a growing, detaching, and rising bubble. Mukherjee and Kandikar have done a numerical simulation of an evaporating meniscus on a moving heated surface with level set method. The heat transfer associated with the advancing and receding contact angle has been studied. To develop the capability of physics-based boiling process, in the presented work the energy and momentum equations are coupled to simulate single-bubble growth scenarios. For a single bubble simulation in infinite liquid, the analytical solution exists and can be compared with the numerical results for verification purposes. The presented capabilities will be utilized in the large scale DNS of boiling flows in the future. The finite-element based code, PHASTA is used for the simulation. It is a parallel, hierarchic, higher-order accurate, adaptive, stabilized, transient analysis flow solver, which has been shown to be an effective tool for a wide range of single and two-phase flow problems. It has been shown to be highly scalable on top supercomputers ({approx}85% scaling on up to 3 million mesh partitions using a 92 billion element mesh . PHASTA has been shown to reliably predict the details of adiabatic bubbly flows as well as single-phase flows with heat transfer. To develop the capability of physics-based boiling process, in the presented work the energy and momentum equations are coupled to simulate a single-bubble growth scenario. For a single bubble growth in superheated water, the analytical solution exists and can be compared with the numerical results for verification purposes. The second-order accurate spatial and time discretization is used for the flow solver in the presented work. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Dynamics of collision-induced energy transfer

Collisions are of fundamental importance in many relevant physico-chemical phenomena in the gas phase and at the gas-surface interface. Collision-induced energy transfer not only offers a sensitive probe of interatomic interactions, but also impacts an array of processes ranging from transport to reactivity. This perspective highlights some recent advances on collision-induced energy transfer dynamics. In the gas phase, cold collisions revealed remarkable quantum effects and stereodynamics of bimolecular inelastic scattering, which motivated development of new state-to-state quantum scattering methods. For surface processes, experiments of hydrogen atom scattering from various surfaces combined with first-principles molecular dynamics simulations shed valuable light on the relative importance of adiabatic energy transfer to surface phonons and nonadiabatic energy transfer to surface electrons. We outline further challenges in these fields and expect more fruitful interplay between experiment and theory.

Cold collision↗
DOE OSTI2020

Effect of magnetic field ripple on parallel electron transport during microwave plasma heating in the Proto-MPEX linear plasma device

The Prototype Material Plasma Exposure eXperiment (Proto-MPEX) is a linear plasma device located at Oak Ridge National Laboratory to develop the plasma source concept for the Materials Plasma Exposure eXperiment (MPEX). Recent experiments have demonstrated the heating of electrons in a high-density deuterium plasma (n e ~ 5 × 10 19 m -3 ) via electron Bernstein waves and source plasma production via helicon waves in this linear configuration. Moreover, experimental observations suggest that the magnetic ripple adversely affects the parallel transport of the heated electrons toward the target where material samples are to be exposed to the plasma. In order to understand the transport process during microwave application, a test-particle Monte-Carlo (MC) code has been developed that incorporates the effects of Coulomb collisions, magnetic mirror adiabatic trapping and electron cyclotron interaction via a quasilinear RF heating operator. MC calculations indicate that the absence of 2nd turning points along the trajectory of cyclotron heated electrons significantly reduces adiabatic trapping and maximizes power transport to the Target surface. Thermalization of cyclotron heating electrons is analyzed with the MC code and it is found that for conditions relevant to ProtoMPEX, a significant fraction of the absorbed heating power is coupled to the Target surface via ballistic fast electrons; however, under certain conditions, up to 60% of the absorbed RF power can be dissipated via thermalization of the fast electrons on the background plasma. Results are compared with experimental measurements. Finally, the implications of the results are discussed in the context of the upcoming MPEX device.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Dinuclear systems in complete fusion reactions

Formation and evolution of dinuclear systems in reactions of complete fusion are considered. Based on the dinuclear system concept, the process of compound nucleus formation is studied. Arguments confirming the validity of this concept are given. The main problems of describing the complete fusion in adiabatic approximation are listed. Calculations of evaporation residue cross sections in complete fusion reactions leading to formation of superheavy nuclei are shown. Isotopic trends of the cross sections of heavy nuclei formation in complete fusion reactions are considered.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1978

The role of O-type neutral lines in magnetic merging during substorms and solar flares

An analytical model of O-type configurations was derived, for visualizing the geometry and for numerical treatment of plasma flows. It was found that at a certain distance from the neutral line the mean particle motion became decoupled from that of magnetic field lines (which obey the MHD condition). The decoupling distance depended on initial conditions in momentum space, suggesting that the MHD approximation which averages out such conditions may not suffice for describing plasma dynamics near the neutral line. It was also found that after inflowing particles are decoupled from the field line motion, they go over to a mode of runaway acceleration along the neutral line. It is concluded that if merging occurs at an X-O pair, two particle populations may be expected - low energy particles accelerated adiabatically by earthward convection past the X-type line and high-energy particles convected towards the O-type line. The second acceleration process depends critically on the rapidity of merging and is therefore expected to vary considerably from event to event.

Stern, D. P.↗
DOE OSTI2016

Drizzle formation in stratocumulus clouds: Effects of turbulent mixing

The mechanism of drizzle formation in shallow stratocumulus clouds and the effect of turbulent mixing on this process are investigated. A Lagrangian–Eularian model of the cloud-topped boundary layer is used to simulate the cloud measured during flight RF07 of the DYCOMS-II field experiment. The model contains ~ 2000 air parcels that are advected in a turbulence-like velocity field. In the model all microphysical processes are described for each Lagrangian air volume, and turbulent mixing between the parcels is also taken into account. It was found that the first large drops form in air volumes that are closest to adiabatic and characterized by high humidity, extended residence near cloud top, and maximum values of liquid water content, allowing the formation of drops as a result of efficient collisions. The first large drops form near cloud top and initiate drizzle formation in the cloud. Drizzle is developed only when turbulent mixing of parcels is included in the model. Without mixing, the cloud structure is extremely inhomogeneous and the few large drops that do form in the cloud evaporate during their sedimentation. Lastly, it was found that turbulent mixing can delay the process of drizzle initiation but is essential for the further development of drizzle in the cloud.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

Symmetrical windowing for quantum states in quasi-classical trajectory simulations: Application to electron transfer

It has recently been shown [S. J. Cotton and W. H. Miller, J. Chem. Phys. 139, 234112 (2013)] that a symmetrical windowing quasi-classical (SQC) approach [S. J. Cotton and W. H. Miller, J. Phys. Chem. A 117, 7190 (2013)] applied to the Meyer-Miller model [H.-D. Meyer and W. H. Miller, J. Chem. Phys. 70, 3214 (1979)] for the electronic degrees of freedom in electronically non-adiabatic dynamics is capable of quantitatively reproducing quantum mechanical results for a variety of test applications, including cases where “quantum” coherence effects are significant. Here we apply this same SQC methodology, within a flux-side correlation function framework, to calculate thermal rate constants corresponding to several proposed models of electron transfer processes [P. Huo, T. F. Miller III, and D. F. Coker, J. Chem. Phys. 139, 151103 (2013); A. R. Menzeleev, N. Ananth, and T. F. Miller III, J. Chem. Phys. 135, 074106 (2011)]. Good quantitative agreement with Marcus Theory is obtained over several orders of magnitude variation in non-adiabatic coupling. Moreover, the “inverted regime” in thermal rate constants (with increasing bias) known from Marcus Theory is also reproduced with good accuracy by this very simple classical approach. The SQC treatment is also applied to a recent model of photoinduced proton coupled electron transfer [C. Venkataraman, A. V. Soudackov, and S. Hammes-Schiffer, J. Chem. Phys. 131, 154502 (2009)] and population decay of the photoexcited donor state is found to be in reasonable agreement with results calculated via reduced density matrix theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗