Search NASA⌕ Search

SEARCH · Search NASA

Results for “Problem structuring methods”

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 811 records · Page 45

tih_vqe [SWR-23-32]

This software supports the paper, "Exploring the scaling limitations of the variational quantum eigensolver with the bond dissociation of hydride diatomic molecules," published in the International Journal of Quantum Chemistry, whose abstract is as follows: Materials simulations involving strongly correlated electrons pose fundamental challenges to state-of-the-art electronic structure methods but are hypothesized to be the ideal use case for quantum computing. To date, no quantum computer has simulated a molecule of a size and complexity relevant to real-world applications, despite the fact that the variational quantum eigensolver (VQE) algorithm can predict chemically accurate total energies. Nevertheless, because of the many applications of moderately-sized, strongly correlated systems, such as molecular catalysts, the successful use of the VQE stands as an important waypoint in the advancement toward useful chemical modeling on near-term quantum processors. In this paper, we take a significant step in this direction. We lay out the steps, write, and run parallel code for an (emulated) quantum computer to compute the bond dissociation curves of the TiH, LiH, NaH, and KH diatomic hydride molecules using VQE. TiH was chosen as a relatively simple chemical system that incorporates d orbitals and strong electron correlation. Because current VQE implementations on existing quantum hardware are limited by qubit error rates, the number of qubits available, and the allowable gate depth, recent studies have focused on chemical systems involving s and p block elements. Through VQE + UCCSD calculations of TiH, we evaluate the near-term feasibility of modeling a molecule with d-orbitals on real quantum hardware. We demonstrate that the inclusion of d-orbitals and the use of the UCCSD ansatz, which are both necessary to capture the correct TiH physics, dramatically increase the cost of this problem. We estimate the approximate error rates necessary to model TiH on current quantum computing hardware using VQE+UCCSD and show them to likely be prohibitive until significant improvements in hardware and error correction algorithms are available.

Graf, Peter↗

Electric and magnetic dipole strength in 112,114,116,118,120,124 Sn

Background: There is renewed interest in electric dipole strength distributions for a variety of reasons including the extraction of the dipole polarizability related to properties of the symmetry energy and a measure for the neutron skin thickness, understanding the structure of low-energy $\textit{E}$1 strength in nuclei with neutron excess, and establishing the systematics of the isovector giant dipole resonance (IVGDR). Inelastic proton scattering at energies of a few hundred MeV and very forward angles including 0° has been established as a tool for the study of electric and magnetic dipole strength distributions in nuclei. Purpose: The present work aims at a systematic investigation of the electric and magnetic dipole strength distributions in the chain of stable even-mass tin isotopes. Methods: Inelastic proton scattering experiments were performed at the Research Center for Nuclear Physics, Osaka, with a 295-MeV beam covering laboratory angles 0°–6° and excitation energies 6–22 MeV. Cross sections due to $\textit{E}$1 and $\textit{M}$1 excitations were extracted with a multipole decomposition analysis (MDA) and then converted to reduced transition probabilities with the “virtual photon method” for $\textit{E}$1 and the “unit cross section method” for $\textit{M}$1 excitations, respectively. Including a theory-aided correction for the high-excitation-energy region not covered experimentally, the electric dipole polarizability was determined from the $\textit{E}$1 strength distributions. Results: In this work, total photoabsorption cross sections derived from the $\textit{E}$1 and $\textit{M}$1 strength distributions show significant differences compared to those from previous ($\textit{γ, xn}$) experiments in the energy region of the IVGDR. The widths of the IVGDR deduced from the present data with a Lorentz parametrization show an approximately constant value of about 4.5 MeV in contrast to the large variations between isotopes observed in previous work. The IVGDR centroid energies are in good correspondence to expectations from empirical systematics of their mass dependence. Furthermore, a study of the dependence of the IVGDR energies on bulk matter properties is presented. The $\textit{E}$1 strengths below neutron threshold show fair agreement in this work with results from ($\textit{γ, γ'}$) experiments on 112,116,120,124 Sn in the energy region between 6 and 7 MeV, where also isoscalar $\textit{E}$1 strength was found for 124 Sn. At higher excitation energies, large differences are observed, pointing to a different nature of the excited states with small ground-state branching ratios. The isovector spin-$\textit{M}$1 strengths exhibit a broad distribution between 6 and 12 MeV in all studied nuclei. Conclusions: The present results contribute to the solution of a variety of nuclear structure problems including the systematics of the energy and width of the IVGDR, the structure of low-energy $\textit{E}$1 strength in nuclei, new constraints to energy density functionals (EDFs) aiming at a systematic description of the dipole polarizability across the nuclear chart, from which properties of the symmetry energy can be derived, and the systematics of the isovector spin-$\textit{M}$1 strength in heavy nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Spectral Analysis of Regular Material Point Method and its Application to Study High Pressure Reverse Osmosis Membrane Compaction and Embossing

Material Point Method (MPM) is gaining widespread interest in applied continuum mechanics. The fact that all the continuum properties are stored on the particles (or material points) and the governing equations are solved on these material points makes MPM extremely suited to problems involving severe material deformations, such as crack propagation, soil movement, and fluid flows. Despite its popularity, only a few studies have focused on the numerical properties of MPM. This presentation introduces a global spectral analysis of the regular material point method. Contrary to previous studies, the analysis focuses on the numerical properties of the method in the spectral space. The amplification factor is derived as a function of the non- dimensional wave numbers. It provides insights into the stability and dissipative properties of the method for various CFL and Fourier numbers. The effect of the grid shape functions, number of particles per cell and their locations inside the grid cell are also analyzed. The EXAGOOP MPM solver (https://github.com/NREL/Exagoop.git) is developed at the National Renewable Energy Laboratory as a part of the NAWI UHPRO project and is based on the AMReX framework. A single-level, uniform cartesian grid is used as the background mesh, while the particle class in AMReX is used to manage the material point operations. Linear hat and B-splines are used as grid shape functions, while the time integration is performed using explicit Euler time integration. EXAGOOP is both CPU and GPU compatible and has been demonstrated to work well on multiple compute architectures. The performance of EXAGOOP on various computing architectures is presented along with its application to study compaction and embossing of high-pressure reverse osmosis membranes. The MPM solution accurately reproduces the membrane deformation. The deformed pore size and structure simulated using MPM also agree well with experimental SEM images.

material point method↗

Demonstration of CFD to support the coupled analysis of a reactor pressure vessel subjected to pressurized thermal shock

The structural components that comprise nuclear reactors and their supporting structures are subjected to harsh operating environments that can challenge their integrity, especially after exposure for extended duration or under accident condition. As one of the most significant components of a reactor, the Reactor Pressure Vessel (RPV) is exposed to an aggressive environment during the operation time (e.g. more than 40 years). Aging degradation mechanisms (e.g. thermo-fatigue) could grow initial defects up to a critical size, increasing the susceptibility to failure in the RPV. The conventional methods are mostly based on simple crack and structure geometries. Very limited studies consider the real conditions of the RPV subjected to a thermal shock due to a Loss of Coolant Accident (LOCA). During a LOCA event, the most severe conditions take place when the emergency core cooling (ECC) water is injected inside the cold legs filled initially with hotter water and/or steam. The rapid cooling of the down-comer and the internal RPV surface followed probably by re-pressurization of the RPV causes large temperature gradients and variation of pressure which induces thermal-mechanical stresses. In order to develop the model for integrity assessment of a reactor pressure vessel (RPV) subjected to pressurized thermal shock (PTS), a multi-physics simulation, which includes the thermo-hydraulic, thermo-mechanical and fracture mechanics analyses is necessary. The prediction of the temperature field is achieved by using computational fluid dynamics (CFD) simulation. In this report, a demonstration CFD standalone simulation is performed to support coupled analysis for Reactor Pressure Vessel (RPV) subjected to Pressurized Thermal Shock (PTS). The study use a simplified computational domain to represents a real RPV. The purpose of the study is to demonstrate the transient temperature response of RPV to ECC injection. The CFD model is built in a robust and efficient way for further coupled calculation. The next steps of this work, including the coupled thermal and tensor mechanics capabilities using Cardinal are expected to be complete by the end of FY21 for the demo problem. After this, into FY22, the capability will be demonstrated for a realistic RPV.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Large-Scale Materials Modeling at Quantum Accuracy: Ab Initio Simulations of Quasicrystals and Interacting Extended Defects in Metallic Alloys

Ab initio electronic-structure has remained dichotomous between achievable accuracy and length-scale. Quantum many-body (QMB) methods realize quantum accuracy but fail to scale. Density functional theory (DFT) scales favorably but remains far from quantum accuracy. We present a framework that breaks this dichotomy by use of three interconnected modules: (i) invDFT: a methodological advance in inverse DFT linking QMB methods to DFT; (ii) MLXC: a machine-learned density functional trained with invDFT data, commensurate with quantum accuracy; (iii) DFT-FE-MLXC: an adaptive higher-order spectral finite-element (FE) based DFT implementation that integrates MLXC with efficient solver strategies and HPC innovations in FE-specific dense linear algebra, mixed-precision algorithms, and asynchronous compute-communication. Furthermore, we demonstrate a paradigm shift in DFT that not only provides an accuracy commensurate with QMB methods in ground-state energies, but also attains an unprecedented performance of 659.7 PFLOPS (43.1% peak FP64 performance) on 619,124 electrons using 8,000 GPU nodes of Frontier supercomputer.

density functional theory↗

Melting of spatially modulated phases at domain wall/surface junctions in antiferrodistortive multiferroics

A physical understanding of the nature of spatially modulated phases (SMPs) in rare-earth-doped antiferrodistortive (AFD) multiferroics and how they behave close to surfaces and interfaces is lacking. Here the emergence of the antiferroelectric (AFE), ferroelectric (FE), or ferrielectric (AFE-FE) spatial modulation in the vicinity of the morphotropic phase transition in La x Bi 1 – x Fe O 3 ( x ~ 0.2 ) is explored on the atomic level using high-resolution scanning transmission electron microscopy (HRSTEM). The suppression, or “melting,” of the AFE-type SMP in the vicinity of the AFD twin wall/surface junction is revealed by HRSTEM in La 0.22 Bi 0.78 Fe O 3 films and explained by the hybrid approach combining Landau-Ginzburg-Devonshire (LGD) phenomenology and the semimicroscopic four-sublattice model (FSM). The LGD-FSM approach reduces the problem of AFE (or AFE-FE) SMP emergence and stability to the thermodynamic analysis of the free-energy functional with AFE, FE, and AFD long-range order parameters and two master parameters: the FE-AFE coupling strength between four neighboring A sites and the nonstoichiometry factor, which are proportional to the variations of La concentration in La x Bi 1 – x Fe O 3 films. Finally, we establish that the surface-induced melting of SMPs and the associated broadening of AFE AFD domain walls minimize the film free energy under certain conditions imposed on the master parameters and gradient energy below the critical value. The observed behavior provides insight into the origin of SMPs in AFD multiferroics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Enhanced MPM framework with multipatch isogeometric analysis for geotechnical applications

Achieving stable stress solutions at large strains using the Material Point Method (MPM) is challenging due to the accumulation of errors associated with geometry discretization, cell-crossing noise, and volumetric locking. Several simplified attempts exist in the literature to mitigate these errors, including higher-order frameworks. However, the stability of the MPM solution in such frameworks has been limited to simple geometries and the single-phase formulation (i.e., neglecting pore fluid). Although never explored, multipatch isogeometric analysis offers desirable qualities to simulate complex geometries while mitigating errors in the MPM. The degree of required high-order spatial integration has also never been investigated to infer a minimum limit for the stability of the stress solution in MPM. This paper presents a general-purpose numerical framework for simulating stable stresses in porous media, capturing both near incompressibility and multiphase interactions. First, the numerical framework is presented considering Non-Uniform Rational B-splines (NURBS) to perform isogeometric analysis (IGA) in MPM. Additionally, a volumetric strain smoothing algorithm is used to alleviate errors associated with volumetric locking. Second, the manifestation of cell-crossing errors is assessed via a series of problems with orders ranging from linear to cubic interpolation functions. Third, the use of NURBS is investigated and verified for problems with circular geometries. Finally, multipatch analysis is deployed to simulate plane strain and 3D penetration in soils, considering nearly incompressible elastoplastic (total stress) analysis and fully-coupled hydro-mechanical (effective stress) analysis. The stability of the solution is also analyzed for different constitutive models. From the results, it can be concluded that the framework using cubic interpolation functions with strain smoothing is the most convenient, presenting stable stress solutions for a broad range of multiphase geotechnical applications.

58 GEOSCIENCES↗

Direct Laser Writing of Multimetal Bifunctional Catalysts for Overall Water Splitting

Water electrolysis is of interest as a sustainable way to produce clean hydrogen and oxygen fuel and help mitigate the rising problems of climate change while meeting global energy demands. High-efficiency, stable, and earth-abundant bifunctional catalysts are needed to enable more effective electrochemical cells for the hydrogen evolution reaction (HER) and oxygen evolution reaction (OER). Here, we investigate the synthesis, composition, performance, and mechanism of multimetal catalysts serving dual functionality in both OER and HER of water electrolysis. Through a laser synthesis method, we synthesized heterogeneous catalysts of nanocrystalline multimetallic alloy pockets embedded within an amorphous oxide matrix. We evaluated the performance and composition of a range of mixed transition-metal oxide materials for both OER and HER, ultimately synthesizing a Cr 0.01 Fe 0.27 Co 0.34 Ni 0.38 O x /C y catalyst that has a stable, high-rate, and competitive overall water splitting performance of 1.76 V at 100 mA cm –2 in an alkaline medium. Using density functional theory to gain insight as the active site and mechanism, we propose that the inclusion of a minor amount of Cr increases the degeneracy of energetic states that lowers the cost of forming the O 2 p–d bond and H 1 s–d bond due to the hybridization of s, p, and d orbitals from Cr. Using a two-electrode water electrolysis cell with a constant potential of 1.636 V to mimic the setup for fuel production, we found the catalyst to be stable at 14–15 mA cm –2 for 40 h. Furthermore, this laser synthesis method allowing for facile and rapid synthesis of complex multimetal systems demonstrates how doping a Fe, Co, and Ni heterogeneous amorphous/nanocrystalline structure with small amounts of Cr is important for bifunctional catalytic behavior, particularly for increasing HER functionality in advancing our understanding for future electrocatalytic design.

08 HYDROGEN↗

Experimental methods for laboratory measurements of helium spectral line broadening in white dwarf photospheres

White Dwarf (WD) stars are the most common stellar remnant in the universe. WDs usually have a hydrogen or helium atmosphere, and helium WD (called DB) spectra can be used to solve outstanding problems in stellar and galactic evolution. DB origins, which are still a mystery, must be known to solve these problems. DB masses are crucial for discriminating between different proposed DB evolutionary hypotheses. Current DB mass determination methods deliver conflicting results. The spectroscopic mass determination method relies on line broadening models that have not been validated at DB atmosphere conditions. We performed helium benchmark experiments using the White Dwarf Photosphere Experiment (WDPE) platform at Sandia National Laboratories' Z-machine that aims to study He line broadening at DB conditions. Using hydrogen/helium mixture plasmas allows investigating the importance of He Stark and van der Waals broadening simultaneously. Accurate experimental data reduction methods are essential to test these line-broadening theories. In this paper, we present data calibration methods for these benchmark He line shape experiments. We give a detailed account of data processing, spectral power calibrations, and instrument broadening measurements. Uncertainties for each data calibration step are also derived. We demonstrate that our experiments meet all benchmark experiment accuracy requirements: WDPE wavelength uncertainties are <1 Å, spectral powers can be determined to within 15%, densities are accurate at the 20% level, and instrumental broadening can be measured with 20% accuracy. Fulfilling these stringent requirements enables WDPE experimental data to provide physically meaningful conclusions about line broadening at DB conditions.

79 ASTRONOMY AND ASTROPHYSICS↗

Space-Time Finite Element Tensor Network Approach for the Time-Dependent Convection–Diffusion–Reaction Equation with Variable Coefficients

In this paper, we present a new space-time Galerkin-like method, where we treat the discretization of spatial and temporal domains simultaneously. This method utilizes a mixed formulation of the tensor-train (TT) and quantized tensor-train (QTT) (please see Section Tensor-Train Decomposition), designed for the finite element discretization (Q1-FEM) of the time-dependent convection–diffusion–reaction (CDR) equation. We reformulate the assembly process of the finite element discretized CDR to enhance its compatibility with tensor operations and introduce a low-rank tensor structure for the finite element operators. Recognizing the banded structure inherent in the finite element framework’s discrete operators, we further exploit the QTT format of the CDR to achieve greater speed and compression. Additionally, we present a comprehensive approach for integrating variable coefficients of CDR into the global discrete operators within the TT/QTT framework. The effectiveness of the proposed method, in terms of memory efficiency and computational complexity, is demonstrated through a series of numerical experiments, including a semi-linear example.

convection–diffusion–reaction equation↗

Aberration corrected RF flipper for high resolution neutron spectroscopy

Project Summary Company: Adelphi Technology, Inc. Title: Aberration-corrected High Frequency RF Flipper for High-Resolution Neutron Spectroscopy PI: Dr. Jay Theodore Cremer Topic: C55-11 Enhancement of Scattering Instrumentation Technology Used at Pulsed and Continuous Sources Subtopic: d. Other Statement of the problem or situation that is being addressed. The quest to understand heterogeneous and hierarchical materials is gathering momentum, as described in a 2015 report by the Basic Energy Sciences Advisory Committee on Challenges at the Frontiers of Matter and Energy. For the past 40 years a technique called neutron spin echo (NSE) has been used to probe molecular motions in such non-crystalline materials over time scales from 10’s of picoseconds to 100’s of nanoseconds. The method has provided unique information about the dynamics of soft heterogeneous materials, including confirmation of the de Gennes model of polymer reptation and quantitative measurement of bending constants of biologically relevant lipid membranes. However, scientists continue to clamor for even higher resolution than NSE can provide. Biomaterials, polymers, glasses, and artificially nanostructured materials all manifest slow molecular motions because of weak or competing interactions between subunits and are amenable to study with neutrons, provided sufficiently long dynamical correlation times can be achieved. All these materials have important applications to advanced technologies so understanding them is key to technological progress. General statement of how this problem is being addressed. We will address the need for high-resolution neutron spectroscopy by using a technique called Neutron Resonance Spin Echo (NRSE). While similar to NSE in many respects, this method has the potential to exceed the NSE capabilities, if 2 technical hurdles can be overcome. The major impediments to successful high-resolution NRSE are the availability of two technologies: a very high frequency, efficient, radiofrequency (rf) flipper for neutrons and a method to correct certain magnetic aberrations. Based on previous STTR support and follow-on research we have developed a suitable rf flipper and we have invented a method to correct the magnetic aberrations. Both technologies need refinement to make them suitable for implementation at a neutron source such as the Oak Ridge National Laboratory nuclear reactor. In this proposal we seek to perfect the two technologies and to combine them into a single, operationally convenient device. Commercial Applications and Other Benefits In view of the increasing demand for the unique scientific information that high resolution neutron spectroscopy can provide, we expect several major instrumentation upgrades at both U.S. and foreign neutron centers will require make use of the NRSE method over the coming decade, creating a market for the devices we will design. These components will enhance scientists’ abilities to probe the time dependence of density fluctuations in a wide range of hierarchical and heterogeneous materials many of which are vital to existing and future technologies. Key Words – Polarized Neutrons, Neutron Spin Echo, Neutron Scattering, advanced materials. Summary for Members of Congress Neutron beams are a powerful materials-science probe that provide unique information about the structure of matter. The proposed devices will accelerate scientific discoveries required to achieve national goals for new technological materials.

36 MATERIALS SCIENCE↗

Multifidelity validation of digital surrogates using variable-density turbulent mixing models

In this study, ensembles of experimental data are presented and utilized to compare and validate two models used in the simulation of variable-density (Atwood = 0.22), compressible turbulent mixing. Though models of this kind (Reynolds averaged NavierStokes and large-eddy simulations) have been validated extensively with more canonical flows in previous studies, here the present approach offers novelty in the complexity of the geometry, the ensemble-based validation, and the uniformity of the computational framework on which the models are tested. Moreover, all experimental and computational tasks were completed by the authors which has led to a tightly coupled experimental configuration with its “digital twin.” The experimental divergent-shock-tube facility and its data acquisition methods are described and replicated in simulation space. A 2D Euler model which neglects the turbulent mixing at the interface is optimized to experimental data using a Gaussian process. This model then serves as the basis for both the 2D RANS and 3D LES studies that make comparisons to the mixing-layer data from the experiment. A relatively simple RANS model is shown to produce good agreement with experimental data only at late flow development times. The LES ensembles generally show good agreement with experimental data but display sensitivity to the characterization of initial conditions. Resolution-dependent behavior is also observed for certain higher-order statistics of interest. Overall, the LES model successfully captures the effects of divergent geometry, compressibility, and combined nonlinear instabilities inherent to the problem. The successful prediction of mixing width and its growth rate highlight the existence of three distinct regimes in the development of the instability, each with similarities to previously studied instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Reduced order models for thermal radiative transfer problems based on moment equations and data-driven approximations of the Eddington tensor

Here a new group of structure and asymptotic preserving reduced-order models (ROMs) for multidimensional nonlinear thermal radiative transfer (TRT) problems is presented. They are formulated by means of the nonlinear projective approach and data compression techniques. The nonlinear projection is applied to the Boltzmann transport equation (BTE) to derive a hierarchy of low-order moment equations. Approximation of the Eddington tensor that provides exact closure for the system of moment equations is found with projection-based data-driven methodologies. These include the (i) proper orthogonal decomposition (POD), (ii) dynamic mode decomposition (DMD) and (iii) a variant of the DMD. A parameterization is derived for this ROM for the temperature of radiation incoming to the problem domain (the radiation drive temperature). This parameterization is informed from results of a dimensionless study of the TRT problem. Analysis of the ROMs is performed on the classical Fleck-Cummings TRT multigroup test problem in 2D geometry with a radiation-driven Marshak wave. Numerical results are presented to demonstrate the performance of these ROMs for the simulation of evolving radiation and heat waves. Results show these models to be sufficiently accurate for practical computations with rather low-rank representations of the Eddington tensor. As the rank of the approximation is increased, the errors of solutions generated by the ROMs gradually decreases.

42 ENGINEERING↗

Optimal design of chemoepitaxial guideposts for the directed self-assembly of block copolymer systems using an inexact Newton algorithm

Directed self-assembly (DSA) of block copolymers (BCPs) is one of the most promising developments in the cost-effective production of nanoscale devices. The process makes use of the natural tendency for BCP melts to form nanoscale structures upon phase separation. The phase separation can be directed through the use of chemically patterned substrates to promote the formation of morphologies that are essential to the production of semiconductor devices. Moreover, the design of substrate pattern can be formulated as an optimization problem for which we seek optimal substrate designs that effectively produce given target morphologies. In this paper, we adopt a phase field model given by a nonlocal Cahn–Hilliard partial differential equation (PDE) based on the minimization of the Ohta–Kawasaki free energy, and present an efficient PDE-constrained optimization framework for the optimal design problem. The design variables are the locations of circular- or strip-shaped guiding posts that are used to model the substrate chemical pattern. To solve the ensuing optimization problem, we propose a variant of an inexact Newton conjugate gradient algorithm tailored to this problem. Additionally, we demonstrate the effectiveness of our computational strategy on numerical examples that span a range of target morphologies. Owing to our second-order optimizer and fast state solver, the numerical results demonstrate five orders of magnitude reduction in computational cost over previous work. The efficiency of our framework and the fast convergence of our optimization algorithm enable us to rapidly solve the optimal design problem in not only two, but also three spatial dimensions.

97 MATHEMATICS AND COMPUTING↗

The Uranium-Containing and Thorium-Containing Anions Studied by Photoelectron Spectroscopy

An in-depth knowledge of actinide chemistry is fundamental to many aspects of nuclear science and technology, including the synthesis and processing of materials and the remediation of waste disposal sites. Among the actinides, the chemical bonding behaviors of actinium and thorium resemble those of the transition metals; the 5f-electrons of protactinium, uranium, neptunium, and plutonium often play important roles in their bonding; and among the still heavier elements, their bonding tends to mimic the lanthanide elements in terms of electron shielding and their f-electron contributions. Bonding that involves 5f-electrons, however, is especially important, in part because of the significance of uranium and plutonium, but also because these elements are among the few where f-electron participation in bonding is relatively common. This work focused on studying uranium-containing and thorium-containing anions in the gas phase using negative ion photoelectron spectroscopy. Since this technique directly probed valence electrons, it was uniquely positioned to address open questions regarding molecular bonding and electron configurations. A particularly important issue concerned how bonding in actinide-containing molecules was affected by modifications to their actinide atoms’ environment, i.e., due to their interaction with ligands. A closely related question was how actinide atoms’ suborbitals were qualitatively reordered and their energies quantitatively shifted as a result of their ligated environments. These were especially relevant issues in regard to uranium due to it having multiple possible oxidation states (OS) and the potential for 5f electron participation in bonding. The effects of ligands on oxidation states and 5f-orbital energies in uranium bonding was expected to be pronounced. Both ligands and excess electrons were seen as probes of actinide atoms within actinide-containing molecules. Our strategy for advancing knowledge of chemical bonding in the actinide-containing species utilized the synergy between experiments and theory, where in some cases experimental results validated theory and where in others computational results assisted in interpreting experiments. Calculations on actinide systems are terrifically challenging due to large spin-orbit interactions, relativistic effects, and just the sheer number of electrons involved. Even in the simplest species, e.g., U and U2, the most sophisticated, modern calculations carried out by the most experienced theorists often only approximate experimentally-measured values, such as electron affinities. For theory to provide confident predictions that can be used to solve real problems it needed an iterative and ultimately corrective mechanism by which its methods can develop further. Experiments can be used to identify when theory has failed; whereupon the subsequent process of using the experiment-theory interplay can be used to find the cause of the failure. Upon fixing it in one case, different test species can be proposed and studied by the experiment-theory combination to determine whether the problem has been corrected. Thus, experiments not only measure the values of molecular properties, they also provide navigational 3 beacons that keep computations off the reefs in an otherwise dark sea with few reference points. Experimental measurements in the actinide field are not only important, they are in actuality essential to computational progress. While it was not always possible to compare the theoreticallydetermined quantity of interest directly with the same experimentally-measured observable, it was usually possible to compare consequential properties that are both calculable and measurable. In the work completed here electron affinities and electronic state spacings were often sensitive consequential parameters. Reasonable agreement between measured and computational values signaled that a calculation that was very likely to be on-track. We had established collaborative relationships with five computational groups, all of which have expertise in computational actinide chemistry. Their PI’s are L. Cheng, D. Dixon, L. Gagliardi, K. Peterson, and B. Vlaisavljevich. Our close interaction with our theory partners led to us suggesting systems to them and them to us. This reciprocal interaction between our experimental and their computational results was among the most important strengths of this work and was a thread woven throughout. Even though anion photoelectron spectroscopic studies are conducted on anions, much of the information that they provide, pertains to the electronic structure of the neutral counterparts of those anions; among these are electron affinities and electronically excited state spacings. Our experimental tools included several specialized ion sources for forming the anionic species of interest, a mass spectrometer for identifying and mass-selecting them, and an anion photoelectron spectrometer for determining their electron affinities (EA) and characterizing the electronic states of the selected anions’ neutral counterparts. Anion photoelectron spectroscopy is conducted by crossing a mass-selected beam of anions with a fixed-frequency laser beam and energy-analyzing the resultant photodetached electrons. The photodetachment process is governed by the energyconserving relationship: hν = EBE + EKE, where hν is the photon’s energy, EBE is the electron binding (photodetachment transition) energy, and EKE is the electron’s kinetic energy. In our apparatus mass-selection is accomplished via time-of-flight mass spectrometry (TOF-MS), electron energy analysis is achieved with either a magnetic bottle or by velocity mapped imaging. Photodetachment of electrons from anions is implemented via either Nd:YAG or excimer lasers. The photodetachment transition energy, i.e., the EBE, between the ground vibrational and electronic state of an anion and the ground vibrational and electronic state of that anion’s neutral counterpart is the adiabatic electron affinity (EA) of that neutral molecule. Likewise, photodetachment transitions between the ground vibrational and electronic state of an anion and the various electronically-excited states of that anion’s corresponding neutral map the electronic spectrum of that neutral species, i.e., the spectral spacings in the photoelectron spectrum are a mirror image of the neutral’s electronic spectrum. It was, of course, crucial to be able to form the anionic species of interest. There, we had a particularly broad field of anion sources from which to choose. These included several variants of pulsed laser vaporization (LV), laser photoemission, infrared desorption plus photoemission, pulsed arc discharge (PACIS), electrospray ionization (ESI), and Rydberg electron transfer (RET). Each of these anion sources were readily combined with, i.e., connected to, the anion photoelectron spectroscopic portion of our apparatus as described above. Among the sources that utilize lasers, visible light for LV sources as well as IR for desorption sources are provided by Nd:YAG lasers. Ultraviolet photons are provided by both Nd:YAG and excimer lasers, whereas the excitation wavelengths for RET experiments come from two Nd:YAG-pumped dye lasers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Benzenesulfonamide Derivatives as Complexants and Extractants for Addressing the Mercury Problem at the Savannah River Site

Mercury (Hg) is a major global pollutant arising from both natural and anthropogenic sources. Its widespread use in medicinal and industrial applications makes it a common chemical exposure and environmental pollutant. It can exist in several forms which include: Metallic mercury (Hg{sup 0}), mercurous (Hg{sub 2}{sup 2+}), mercuric salts (Hg{sup 2+}), and organic mercury (e.g. CH{sub 3}Hg{sup +}), with the latter being the most toxic of all the species. Due to mercury's high toxicity, new approaches towards its detection has received significant attention in the scientific community. Mercury exposure at the Savannah River Site (SRS) has been a recent concern especially with increasing amounts of organic mercury in the saltstone. It originates mainly from its use as an acidic dissolution catalyst of aluminum cladding from target fuels within the uranium and plutonium processing operations [1]. It is present to an amount of about 60 metric tons in the high-level waste (HLW) tanks. Organic Mercury species have been found in low activity waste (LAW) at the site that eventually ends up in the saltstone. Therefore, there is a need for: i) Converting organic mercury to other less toxic forms and ii) Complexation and removal of Hg prior of disposal of LAW in saltstone. Various methods have been developed for selective sensing of mercury in the presence of other toxic metals. These methods include using ligands that can form organo-soluble metal complexes with different optical and spectroscopic properties that can be used for toxic metal sensing. In 2005, our group pioneered an ion-exchange extraction method, in which o-phenylenediamine-derived disulfonamides were used to complex and selectively extract and sense Pb{sup 2+} from aqueous solutions into an organic phase [2,3]. Herein, a disulfonamide and a bis-dansylamide have been shown to extract, complex and sense Hg(II). Ligand 1: The crystal structure of the disulfonamide-Hg complex confirms the complexation of Hg(II) with the ligand. Complexation was corroborated by the {sup 1}H-NMR spectra obtained after contacting solutions of various concentrations of Hg{sup 2+} with 2 mM ligand in chloroform. Distinct resonances are observed at Hg/L ratio of 0.5 that are also observed for the isolated 1:2 complex. In the presence of excess mercury, new resonances, as well as the movement of Et{sub 3}N resonances indicate the formation of a different Hg-sulfonamide-triethylamine complex, presumably having 1:1 Hg:L stoichiometry. The electronic spectra of aqueous phases after extraction show that there was no free ligand absorption at 0.5 eq of Hg, indicating a complete complexation. Complexation was also confirmed by the UV-visible titrations with Hg{sup 2+} at constant ligand concentration. pH-dependent extraction carried out shows that extraction of Hg(II) by ligand 1 was over 90% for most alkaline pHs. Ligand 2: The crystal structure of the Ligand 2 complex formed with Hg(OAc){sub 2} shows a remarkable coordination pattern with 4:2 metal:ligand stoichiometry. The fluorescent bis-dansyl disulfonamide derivative was found to complex and sense HgCl{sub 2} and Hg(OAc){sub 2} by demonstrating fluorescence quenching upon Hg(II) addition in comparison with other metals (Zn(II), Cd(II), Pb(II)). No were observed for Cu(II), Ag(I) and Co(II). We have shown the complexation of Hg(II) by a disulfonamide and a bis-dansyl disulfonamide ligand using several spectroscopic methods. Ligand 1 was able to extract mercury into chloroform and form a complex in the presence of excess mercury by synergistic complexation with triethylamine acting as a co-ligand. X-ray and NMR both confirm a 1:2 HgL{sub 2} stoichiometry. Ligand 2 can be used for sensing of Hg(II) as fluorescence quenching was observed after addition of Hg(II)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Impact of Numerical Hydrodynamics in Turbulent Mixing Transition Simulations

Underresolved simulations are unavoidable in high Reynolds (Re) and Mach (Ma) number turbulent flow applications at scale. Implicit large-Eddy simulation (ILES) often becomes the effective strategy to capture the dominating effects of convectively driven flow instabilities. We evaluate the impact of three distinct numerical strategies in simulations of transition and turbulence decay with ILES: the Harten–Lax–van Leer (HLL) Riemann solver applying Strang splitting and a Lagrange-plus-Remap formalism to solve the directional sweep—denoted split; the Harten–Lax–Van Leer-Contact (HLLC) Riemann solver using a directionally unsplit strategy and parabolic reconstruction—denoted unsplit; and the HLLC Riemann solver using unsplit and a low-Ma correction (LMC)—denoted unsplit*. Three case studies are considered: (1) a shock tube problem prototyping shock-driven turbulent mixing, (2) the Taylor–Green Vortex (TGV) prototyping transition to turbulence, and, (3) an homogeneous isotropic turbulence (HIT) case, focusing on the impact of discretization on transition and decay from fixed well-characterized initial conditions. Significantly more accurate predictions are provided by the unsplit schemes, in particular, when augmented with the LMC. For given resolution, only the unsplit schemes predict the turbulent mixing transition after reshock observed in the shock tube experiments. Relevant comparisons of ILES based on Euler and Navier–Stokes equations addressing potential occurrence of low-Re regimes in the applications are presented. Unsplit schemes are instrumental in allowing to capture the spatial development of the TGV flow and its validation at prescribed Re with significantly less resolution. HIT analysis confirms higher simulated turbulence Re and increased small-scale content associated with the unsplit discretizations.

74 ATOMIC AND MOLECULAR PHYSICS↗

A coupled discontinuous Galerkin-Finite Volume framework for solving gas dynamics over embedded geometries

Herein, we present a computational framework for solving the equations of inviscid gas dynamics using structured grids with embedded geometries. The novelty of the proposed approach is the use of high-order discontinuous Galerkin (dG) schemes and a shock-capturing Finite Volume (FV) scheme coupled via an hp adaptive mesh refinement (hp-AMR) strategy that offers high-order accurate resolution of the embedded geometries. The hp-AMR strategy is based on a multi-level block-structured domain partition in which each level is represented by block-structured Cartesian grids and the embedded geometry is represented implicitly by a level set function. The intersection of the embedded geometry with the grids produces the implicitly-defined mesh that consists of a collection of regular rectangular cells plus a relatively small number of irregular curved elements in the vicinity of the embedded boundaries. High-order quadrature rules for implicitly-defined domains enable high-order accuracy resolution of the curved elements with a cell-merging strategy to address the small-cell problem. The hp-AMR algorithm treats the system with a second-order finite volume scheme at the finest level to dynamically track the evolution of solution discontinuities while using dG schemes at coarser levels to provide high-order accuracy in smooth regions of the flow. On the dG levels, the methodology supports different orders of basis functions on different levels. The space-discretized governing equations are then advanced explicitly in time using high-order Runge-Kutta algorithms. Numerical tests are presented for two-dimensional and three-dimensional problems involving an ideal gas. The results are compared with both analytical solutions and experimental observations and demonstrate that the framework provides high-order accuracy for smooth flows and accurately captures solution discontinuities.

97 MATHEMATICS AND COMPUTING↗