Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Doubling method”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Reduced-Cost Four-Component Relativistic Double Ionization Potential Equation-of-Motion Coupled-Cluster Approaches with 4-Hole–2-Particle Excitations and Three-Body Clusters

The double ionization potential (DIP) equation-ofmotion (EOM) coupled-cluster (CC) method with 4-hole−2- particle (4h-2p) excitations on top of the CC with singles, doubles, and triples calculation, abbreviated as DIP-EOMCCSDT(4h-2p), along with its perturbative DIP-EOMCCSD(T)(a)(4h-2p) approximation, are extended to a relativistic four-component (4c) framework. In addition, we introduce and test a new computationally practical DIP-EOMCC approach, which we call DIPEOMCCSD( T)(ã)(4h-2p), that approximates the treatment of 4h- 2p correlations within the DIP-EOMCCSD(T)(a)(4h-2p) method and reduces the $\mathcal{N}$ 8 scaling characterizing DIP-EOMCCSDT(4h- 2p) and DIP-EOMCCSD(T)(a)(4h-2p) to $\mathcal{N}$ 7 with the system size $\mathcal{N}$. Further improvements in computational efficiency are obtained using the frozen natural spinor (FNS) approximation to reduce the numbers of unoccupied spinors entering the correlated steps of the DIP-EOMCC calculations according to a well-defined occupation-number-based threshold. The resulting 4c-FNS-DIPEOMCC approaches are used to compute DIPs for the series of inert gas atoms from argon to radon as well as the vertical DIPs in Cl 2 , Br 2 , HBr, and HI, which have been experimentally examined in the past. We demonstrate that, when using complete basis set extrapolations and FNS truncation threshold of 10 −4.5 , the 4c-FNS-DIP-EOMCCSD(T)(ã)(4h-2p) calculations are capable of predicting DIPs in agreement with experimental data, improving upon their nonrelativistic and spin-free scalar-relativistic counterparts, particularly when examining DIPs characterized by stronger spin−orbit coupling effects.

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

Adhesively-bonded Metal-CFRTP Bi-materials: Enhanced Crack Growth Resistance via Plasma and Quantified Fracture via Size Effect Method

This work investigated air plasma effect on the Mode I crack growth resistance of adhesively-bonded metal-CFRTP dissimilar joints by using aluminum alloy (AA6061) and short-carbon-fiber-reinforced polyamide 66 (CFRPA66) as an example, and the Double Cantilever Beam (DCB) fracture testing as an evaluation method. The results show that air plasma treatment can significantly improve the fracture resistance of adhesively-bonded AA6061-CFRPA66 dissimilar joints with about 140% enhanced Mode I fracture energy in maximum compared to non-treated ones. The quantified fracture energies from size effect method for both non-treated and plasma-treated cases are geometry-independent, whereas this is not true for modified beam theory causing geometry-dependent results. The foregoing improvement was confirmed from the failure surface morphology of plasma-treated specimens, showing fibers peeling off from CFRPA66 surface due to plasma-enhanced bonding between CFRPA66 and adhesive by the formation of covalent bonds at their interface. This study demonstrated air plasma as an efficient surface modification method to enhance the fracture resistance of adhesively-bonded metal-CFRTP dissimilar joints, and size effect method as a characterization method to properly quantify their fracture properties. These aspects are valuable in the area of multi-materials joining.

Plasma, Adhesive Bonding, Slow Crack Growth, fract↗
DOE OSTI2021

Texture Memory in Hexagonal Metals and Its Mechanism

Texture memory is a phenomenon in which retention of initial textures occurs after a complete cycle of forward and backward transformations, and it occurs in various phase-transforming materials including cubic and hexagonal metals such as steels and Ti and Zr alloys. Texture memory is known to be caused by the phenomena called variant selection, in which some of the allowed child orientations in an orientation relationship between the parent and child phases are preferentially selected. Without such variant selection, the phase transformations would randomize preferred orientations. In this article, the methods of prediction of texture memory and mechanisms of variant selections in hexagonal metals are explored. The prediction method using harmonic expansion of orientation distribution functions with the variant selection in which the Burgers orientation relationship, {110}β//{0001} α-hex <1$\bar{1}$1>β//$\langle2\bar{11}0\rangle$ α-hex , is held with two or more adjacent parent grains at the same time, called “double Burgers orientation relation (DBOR)”, is introduced. This method is shown to be a powerful tool by which to analyze texture memory and ultimately provide predictive capabilities for texture changes during phase transformations. Variation in nucleation and growth rates on special boundaries and an extensive growth of selected variants are also described. Analysis of textures of commercially pure Ti observed in situ by pulsed neutron diffraction reveals that the texture memory in CP-Ti is indeed quite well predicted by consideration of the mechanism of DBOR. The analysis also suggests that the nucleation and growth rates on the special boundary of 90° rotation about $\langle2\bar{11}0\rangle$ α-hex should be about three times larger than those of the other special boundaries, and the selected variants should grow extensively into not only one parent grain but also other grains in α-hex(hexagonal)→β(bcc) transformation. The model calculations of texture development during two consecutive cycles of α-hex→β→α-hex transformation in CP-Ti and Zr are also shown.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Ensemble methods for quantification of potassium oxide in ChemCam Mars and laboratory spectra

In this paper we test new approaches for predicting the amount of element oxides in rock samples from the ChemCam instrument suite onboard the NASA Curiosity rover by focusing on K 2 O. Using the expanded dataset compiled by Gasda et al. (2021) with and without the Earth to Mars (E2M and NoE2M) transformation discussed in Clegg et al. (2017) we trained blended submodels using the “double blending” technique and compared these to ensemble methods (Random Forest, ExtraTrees, and Gradient Boosting Regression). We found that ensemble methods performed similar to blended submodels when looking at RMSE-P on the laboratory spectra and provided significant advantages when looking at spectra coming from Mars. For the full model, blended submodels achieved an RMSE-P of 0.62 and 0.60 (E2M and NoE2M respectively) while Gradient Boosting Regression resulted in a slightly improved RMSE-P of 0.59 and 0.60. More importantly, by employing a local RMSE-P estimation technique where model performance is evaluated based on nearby test samples we found that using ensemble methods can lower the quantification limit for K 2 O from the current value of ≈0.6 wt% to ≈0.08 wt% using Extra Trees and Random Forest. This would allow for a much larger range of K 2 O values to be quantified on Mars with greater certainty given that most targets seen on Mars tend to have <1 wt% K2O. Finally, we used both Mean Decrease in Impurity (MDI) and permutation importance techniques to investigate the wavelengths used by the ensemble methods and found that they correspond to known potassium emission lines. This suggests that ensemble methods can provide an easier to train and improved alternative to blended submodels for predicting potassium compositions from Laser Induced Breakdown Spectroscopy (LIBS) data.

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

A bi-layer barrier design for 122-type iron-based superconducting wires and tapes

Iron-based superconducting wires and tapes hold great promise for high-field magnet applications. A promising design for 122-type wires and tapes based on the powder-in-tube method is using silver and copper double-layer sheaths. For this design a heat treatment temperature below ~ 779 °C is required to prevent Ag-Cu liquid formation. However, this may be below the optimal heat treatment temperature for the critical current density, and still cannot prevent Ag-Cu interdiffusion occurring in the solid state. In this work we propose adding a niobium or tantalum or vanadium (or their alloys) barrier layer between the Ag and Cu to solve the Ag-Cu interdiffusion issue, given that the group-VB metals (vanadium, niobium, tantalum) are relatively inert to both Ag and Cu. To investigate the effectiveness of this design, BaFe 1.84 Co 0.16 As 2 wires and tapes with Ag/Cu and Ag/Ta/Cu sheaths, as well as Ba 0.6 K 0.4 Fe 2 As 2 wires and tapes with Ag/Cu and Ag/Nb/Cu sheaths, were fabricated. It was found that both the Ta and Nb layers kept integral after wire drawing, but after a large flat-rolling reduction the Ta layer broke while the Nb layer kept integral. In the tapes with Ag/Cu sheaths (without the Ta or Nb layer) Cu diffused through the Ag layer and into the powder cores during 740 °C heat treatment, while in the tapes with Ag/Nb/Cu sheaths the Nb layer effectively blocked Ag-Cu interdiffusion even at 900 °C. Further, this work demonstrates that Ta is a suitable barrier material for 122-type wires, while Nb is suitable for both wires and tapes. In this design using Ag/Nb (or Ta)/Cu sheaths, we can regard the outer Cu as the conductor matrix while the Ag and Nb (or Ta) serve as two layers of barriers that suppress reactions between the components. Thus, we call this design a “bi-layer barrier” design for 122-type wires and tapes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

An interactive machine learning platform for analyzing multi-particle coincidence data from cold target recoil ion momentum spectroscopy

We present SCULPT (Supervised Clustering and Uncovering Latent Patterns with Training), a comprehensive software platform for analyzing tabulated high-dimensional multi-particle coincidence data from Cold Target Recoil Ion Momentum Spectroscopy (COLTRIMS) experiments. The software addresses critical challenges in modern momentum spectroscopy by integrating advanced machine learning techniques with physics-informed analysis in an interactive web-based environment. SCULPT implements uniform manifold approximation and projection for non-linear dimensionality reduction to reveal correlations in high-dimensional data. We also discuss potential extensions to deep autoencoders for feature learning and genetic programming for automated discovery of physically meaningful observables. A novel adaptive confidence scoring system provides quantitative reliability assessments by evaluating user-selected clustering quality metrics with predefined weights that reflect each metric’s robustness. The platform features configurable molecular profiles for different experimental systems, interactive visualization with selection tools, and comprehensive data filtering capabilities. Utilizing a subset of SCULPT’s capabilities, we analyze photo-double-ionization data measured using the COLTRIMS method for three-body dissociation of the D 2 O molecule, revealing distinct fragmentation channels and their correlations with physics parameters. The software’s modular architecture and web-based implementation make it accessible to the broader atomic and molecular physics community, significantly reducing the time required for complex multi-dimensional analyses. This opens the door to finding and isolating rare events exhibiting non-linear correlations on the fly during experimental measurements, which can help steer exploration and improve the efficiency of experiments.

Artificial neural networks↗
DOE OSTI2020

Effects of Testosterone Supplementation on Ghrelin and Appetite During and After Severe Energy Deficit in Healthy Men

Abstract Background Severe energy deficits cause interrelated reductions in testosterone and fat free mass. Testosterone supplementation may mitigate those decrements, but could also reduce circulating concentrations of the orexigenic hormone ghrelin, thereby exacerbating energy deficit by suppressing appetite. Objective To determine whether testosterone supplementation during severe energy deficit influences fasting and postprandial ghrelin concentrations and appetite. Design and methods Secondary analysis of a randomized, double-blind trial that determined the effects of testosterone supplementation on body composition changes during and following severe energy deficit in nonobese, eugonadal men. Phase 1 (PRE-ED): 14-day run-in; phase 2: 28 days, 55% energy deficit with 200 mg testosterone enanthate weekly (TEST; n = 24) or placebo (PLA; n = 26); phase 3: free-living until body mass recovered (end-of-study; EOS). Fasting and postprandial acyl ghrelin and des-acyl ghrelin concentrations and appetite were secondary outcomes measured during the final week of each phase. Results Fasting acyl ghrelin concentrations, and postprandial acyl and des-acyl ghrelin concentrations increased in PLA during energy deficit then returned to PRE-ED values by EOS, but did not change in TEST (phase-by-group, P < 0.05). Correlations between changes in free testosterone and changes in fasting acyl ghrelin concentrations during energy deficit (ρ = -0.42, P = 0.003) and body mass recovery (ρ = -0.38; P = 0.01) were not mediated by changes in body mass or body composition. Transient increases in appetite during energy deficit were not affected by testosterone treatment. Conclusions Testosterone supplementation during short-term, severe energy deficit in healthy men prevents deficit-induced increases in circulating ghrelin without blunting concomitant increases in appetite. Clinical Trials Registration www.clinicaltrials.gov NCT02734238 (registered 12 April 2016).

Karl, J. Philip↗
DOE OSTI2023

Construction and Commissioning of Gas Electron Multiplier (GEM) Detectors in Advanced Assembly Design for Low-Energy Applications at High Rates and Analysis of GEM Data From the Muse Experiment at PSI

The search beyond the Standard Model explores dark matter and a potential fifth force. DarkLight@ARIEL investigates the dark photon as a mediator between ordinary and dark matter. The experiment aims to measure the process e?Ta ? e?TaX ? e?Ta(e?e+), by detecting a charged lepton pair in the final state. The spectrometers will be instrumented with Gas Electron Multiplier (GEM) detectors with minimal material budget for tracking. A novel GEM construction technique is employed for fabricating 25 cm x 40 cm GEM detectors, where all layers are mechanically stretched and assembled within a double frame. The dissertation outlines physics motivation, methods, experimental setup, and the role of GEM detectors in DarkLight. A comprehensive discussion on GEM detectors, focusing on the "NS2" technique, including design, parts, assembly, testing, and performance evaluation will be discussed. The proton has garnered attention due to discrepancies in measuring its charge radius using muonic hydrogen and electron-based methods. The 2010 muonic hydrogen measurement of the proton charge radius Rp = 0.84184(67) fm, showed a significant 7? discrepancy compared to the previously known value of Rp = 0.8775(51) fm, giving rise to the proton radius puzzle. Through simultaneous measurements of e ? p and ? ? p elastic scattering, MUSE facilitates a precise and direct comparison of the proton radius. The physics background of the proton radius puzzle, measurement techniques, the approach adopted by the MUSE experiment, its significance and the experimental setup will be discussed. Furthermore, the role of GEM detectors, along with an in-depth analysis of their data and efficiency will be discussed and evaluated.

Mohammed Prem Nazeer, Sahara↗
DOE OSTI2023

Defect-Promoted Ni-Based Layer Double Hydroxides with Enhanced Deprotonation Capability for Efficient Biomass Electrooxidation

Ni-based hydroxides are promising electrocatalysts for biomass oxidation reactions, supplanting the oxygen evolution reaction (OER) due to lower overpotentials while producing value-added chemicals. The identification and subsequent engineering of their catalytically active sites are essential to facilitate these anodic reactions. Herein, the proportional relationship between catalysts’ deprotonation propensity and Faradic efficiency of 5-hydroxymethylfurfural (5-HMF)-to-2,5 furandicarboxylic acid (FDCA, FEFDCA) is revealed by thorough density functional theory (DFT) simulations and atomic-scale characterizations, including in situ synchrotron diffraction and spectroscopy methods. The deprotonation capability of ultrathin layer-double hydroxides (UT-LDHs) is regulated by tuning the covalency of metal (M)-oxygen (O) motifs through defect site engineering and selection of M 3+ co-chemistry. NiMn UT-LDHs show an ultrahigh FE FDCA of 99% at 1.37 V versus reversible hydrogen electrode (RHE) and retain a high FE FDCA of 92.7% in the OER-operating window at 1.52 V, about 2× that of NiFe UT-LDHs (49.5%) at 1.52 V. Ni–O and Mn–O motifs function as dual active sites for HMF electrooxidation, where the continuous deprotonation of Mn–OH sites plays a dominant role in achieving high selectivity while suppressing OER at high potentials. The results showcase a universal concept of modulating competing anodic reactions in aqueous biomass electrolysis by electronically engineering the deprotonation behavior of metal hydroxides, anticipated to be translatable across various biomass substrates.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Improving Chirped Fiber Bragg Grating Resolution for Position-Sensitive Sensors in Shock- and Detonation-Driven Experiments

Chirped fiber Bragg gratings (CFBGs) are robust diagnostic sensors that are widely used to track detonation-driven and shock wave propagation. CFBGs are inscribed with a linearly chirped periodic index of refraction changes that alter the Bragg wavelength along the length of the probe. The light return of each individual Bragg element is captured by a detector at a unique time to map the full reflected spectrum. The CFBG spectrum is measured with a dispersive Fourier transform of the reflected light that temporally stretches the spectrum to increase spatial resolution and make a one-to-one map of the wavelength on a time axis. Here, we propose an improvement of CFBG temporal resolution by incorporating two co-linear laser pulses with orthogonal polarization states and a 5 ns time offset. The two separate signals were split and tracked by two separate detectors. An oscilloscope captured good separation in the signals, and two separate spectrograms were generated and interleaved in the post-processing of the data. This novel technique doubled the CFBG temporal resolution and led to a doubled location resolution. As a proof-of-concept of this technique, the resolution improvement was compared between standard CFBG measurements and the two polarization states method on a position-sensitive CFBG sensor. CFBG resolution doubling will advance sensor capabilities and will have a direct impact on improving capture and analysis in dynamic, high-explosive experiments.

42 ENGINEERING↗
DOE OSTI2020

Determination of background doping polarity of unintentionally doped semiconductor layers

Here, we present a method of determining the background doping type in semiconductors using capacitance–voltage measurements on overetched double mesa p–i–n or n–i–p structures. Unlike Hall measurements, this method is not limited by the conductivity of the substrate. By measuring the capacitance of devices with varying top and bottom mesa sizes, we were able to conclusively determine which mesa contained the p–n junction, revealing the polarity of the intrinsic layer. This method, when demonstrated on GaSb p–i–n and n–i–p structures, concluded that the material is residually doped p-type, which is well established by other sources. The method was then applied to a 10 monolayer InAs/10 monolayer AlSb superlattice, for which the doping polarity was unknown, and indicated that this material is also p-type.

42 ENGINEERING↗
DOE OSTI2021

Coupling 1D xRAGE simulations with machine learning for graded inner shell design optimization in double shell capsules

Advances in machine learning provide the ability to leverage data from expensive simulations of high-energy-density experiments to significantly cut down on computational time and costs associated with the search for optimal target designs. This study presents an application of cutting-edge Bayesian optimization methods to the one-dimensional (1D) design optimization of double shell graded layer targets for inertial confinement fusion experiments. This investigation attempts to reduce hydrodynamic instabilities while retaining high yields for future NIF experiments. Machine learning methods can use predictive physics simulations to identify graded layer designs from within the vast design space that demonstrate high predicted performance, including novel designs with high uncertainty in performance that may hold unexpected promise. By applying machine learning tools to the simulation design, we map the trade-off between 1D yield and instability, specifically isolating parameter ranges, which maintain high performance while showing significantly improved Rayleigh–Taylor stability over the point design. Furthermore, the groundwork laid in this study will be a useful design tool for future NIF experiments with graded layer targets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Improving the representation of shallow cumulus convection with the simplified-higher-order-closure–mass-flux (SHOC+MF v1.0) approach

Abstract. Parameterized boundary layer turbulence and moist convection remain some of the largest sources of uncertainty in general circulation models. High-resolution climate modeling aims to reduce that uncertainty by explicitly attempting to resolve deep moist convective motions. An example of such a model is the Simple Cloud-Resolving E3SM Atmosphere Model (SCREAM) with a target global resolution of 3.25 km, allowing for a more accurate representation of complex mesoscale deep convective dynamics. Yet, small-scale planetary boundary layer turbulence and shallow convection still need to be parameterized, which in SCREAM is accomplished through the turbulent-kinetic-energy-based (TKE-based) simplified higher-order closure (SHOC) – a simplified version of the assumed-double-Gaussian-PDF (probability density function) higher-order-closure method. In this paper, we implement a stochastic-multiplume mass-flux (MF) parameterization of dry and shallow convection in SCREAM to go beyond the limitations of double-Gaussian-PDF closures and couple it to SHOC (SHOC+MF). The new parameterization implemented in a single-column model type version of SCREAM produces results for two shallow cumulus convection cases (marine and continental shallow convection) that agree well with the reference data from large-eddy simulations, thus improving the general representation of the thermodynamic quantities and their turbulent fluxes as well as cloud macrophysics in the model. Furthermore, SHOC+MF parameterization shows weak sensitivity to the vertical grid resolution and model time step.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

In-Situ Visualization of a Growing Brittle Crack in Aluminum Oxynitride Using Synchrotron X-Rays and the Double-Cleavage Drilled Compression Geometry

Brittle fracture is difficult to study in situ due to the speed of a growing crack and the often-catastrophic nature of failure in brittle materials. As a result, the influence of microstructural considerations, such as orientation, grain boundary locations, and strain field, on the crack path remains poorly understood. Presented in this study is a method addressing this knowledge gap, which utilizes the double-cleavage drilled compression geometry to achieve quasi-stable fracture in aluminum oxynitride (AlON). Synchrotron X-ray micro-computed tomography is used to characterize the crack shape and length, while high-energy diffraction microscopy provides information on the strains, orientations, and shapes of grains in the microstructure surrounding the crack tip. During testing, the crack grew in discrete and irregular jumps while the fracture toughness falls within reported ranges. The crack in AlON is found to have no greater tendency to crack intergranularly as compared to transgranularly, and grains which are cracked transgranularly do not display a trend in orientation or stress when compared to those around which the crack followed a grain boundary. The high resolution of the crack path and microstructural data provides a path forward for modeling and understanding 3D brittle fracture.

Gorske, Sara F.↗
DOE OSTI2025

Impact of solvation on the electronic resonances in uracil

Interactions of low-energy electrons with the DNA and RNA nucleobases are known to form metastable states, known as electronic resonances. In this work, we study electron attachment to solvated uracil, an RNA nucleobase, using the orbital stabilization method at the Equation of Motion-Coupled Cluster for Electron Affinities with Singles and Doubles (EOM-EA-CCSD) level of theory with the Effective Fragment Potential (EFP) solvation method. We benchmarked the approach using multireference methods, as well as by comparing EFP and full quantum calculations. The impact of solvation on the first one particle (1p) shape resonance, formed by electron attachment to the π* LUMO orbital, as well as the first two particle one hole (2p1h) resonance, formed by electron attachment to neutral uracil's π–π* excited state, was investigated. We used molecular dynamics simulations for solvent configurations and applied charge stabilization technique-based biased sampling to procure configurations adequate to cover the entire range of the electron attachment energy distribution. The electron attachment energy in solution is found to be distributed over a wide range of energies, between 4.6 eV to 6.8 eV for the 2p1h resonance, and between −0.1 eV to 2 eV for the 1p resonance. The solvent effects were similar for the two resonances, indicating that the exact electron density of the state is not as important as the solvent configurations. Multireference calculations extended the findings showing that solvation effects are similar for the lowest four resonances, further indicating that the specific solute electron density is not as important, but rather the water configurations play the most important role in solvation effects. Lastly, by comparing bulk solvation to clusters of uracil with a few water molecules around it, we find that the impact of microsolvation is very different from that of bulk solvation.

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

Self‐Assembly of Helical Nanofibrous Chiral Covalent Organic Frameworks

Abstract Despite significant progress on the design and synthesis of covalent organic frameworks (COFs), precise control over microstructures of such materials remains challenging. Herein, two chiral COFs with well‐defined one‐handed double‐helical nanofibrous morphologies were constructed via an unprecedented template‐free method, capitalizing on the diastereoselective formation of aminal linkages. Detailed time‐dependent experiments reveal the spontaneous transformation of initial rod‐like aggregates into the double‐helical microstructures. We have further demonstrated that the helical chirality and circular dichroism signal can be facilely inversed by simply adjusting the amount of acetic acid during synthesis. Moreover, by transferring chirality to achiral fluorescent molecular adsorbents, the helical COF nanostructures can effectively induce circularly polarized luminescence with the highest luminescent asymmetric factor ( g lum ) up to ≈0.01.

Tang, Xihao↗
DOE OSTI2020

The Generality of the GUGA MRCI Approach in COLUMBUS for Treating Complex Quantum Chemistry

The core part of the program system COLUMBUS allows highly efficient calculations using variational multireference (MR) methods in the framework of configuration interaction with single and double excitations (MR-CISD) and averaged quadratic coupled-cluster calculations (MR-AQCC), based on uncontracted sets of configurations and the graphical unitary group approach (GUGA). The availability of analytic MR-CISD and MR-AQCC energy gradients and analytic nonadiabatic couplings for MR-CISD enables exciting applications including, e.g., investigations of π-conjugated biradicaloid compounds, calculations of multitudes of excited states, development of diabatization procedures, and furnishing the electronic structure information for on-the-fly surface nonadiabatic dynamics. With fully variational uncontracted spin-orbit MRCI, COLUMBUS provides a unique possibility of performing high-level calculations on compounds containing heavy atoms up to lanthanides and actinides. Crucial for carrying out all of these calculations effectively is the availability of an efficient parallel code for the CI step. Configuration spaces of several billion in size now can be treated quite routinely on standard parallel computer clusters. Emerging developments in COLUMBUS, including the all configuration mean energy multiconfiguration self-consistent field method and the graphically contracted function method, promise to allow practically unlimited configuration space dimensions. Spin density based on the GUGA approach, analytic spin-orbit energy gradients, possibilities for local electron correlation MR calculations, development of general interfaces for nonadiabatic dynamics, and MRCI linear vibronic coupling models conclude this overview.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗