DFTpy : An efficient and object-oriented platform for orbital-free DFT simulations
Not provided.
Chemistry↗
SEARCH · Engineering Papers
Results for “complex systems theory”
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.
A multiscale anisotropic polymer network model coupled with phase field fracture
The study of polymers has continued to gain substantial attention due to their expanding range of applications, spanning essential engineering fields to emerging domains like stretchable electronics, soft robotics, and implantable sensors. These materials exhibit remarkable properties, primarily stemming from their intricate polymer chain network, which, in turn, increases the complexity of precisely modeling their behavior. Especially for modeling elastomers and their fracture behavior, accurately accounting for the deformations of the polymer chains is vital for predicting the rupture in highly stretched chains. Despite the importance, many robust multiscale continuum frameworks for modeling elastomer fracture tend to simplify network deformations by assuming uniform behavior among chains in all directions. Recognizing this limitation, our study proposes a multiscale fracture model that accounts for the anisotropic nature of elastomer network responses. At the microscale, damage in the chains is assumed to be driven by both the chain's entropy and the internal energy due to molecular bond distortions. In order to bridge the stretching in the chains to the macroscale deformation, we employ the maximal advance path constraint network model, inherently accommodating anisotropic network responses. As a result, chains oriented differently can be predicted to exhibit varying stretch and, consequently, different damage levels. To drive macroscale fracture based on damages in these chains, we utilize the micromorphic regularization theory, which involves the introduction of dual local-global damage variables at the macroscale. The macroscale local damage variable is obtained through the homogenization of the chain damage values, resulting in the prediction of an isotropic material response. The macroscale global damage variable is subjected to nonlocal effects and boundary conditions in a thermodynamically consistent phase field continuum formulation. Moreover, the total dissipation in the system is considered to be mainly due to the breaking of the molecular bonds at the microscale. To validate our model, we employ the double-edge notched tensile test as a benchmark, comparing simulation predictions with existing experimental data. Additionally, to enhance our understanding of the fracturing process, we conduct uniaxial tensile experiments on a square film made up of polydimethylsiloxane (PDMS) rubber embedded with a hole and notches and then compare our simulation predictions with the experimental observations. Furthermore, we visualize the evolution of stretch and damage values in chains oriented along different directions to assess the predictive capacity of the model. In conclusion, the results are also compared with another existing model to evaluate the utility of our model in accurately simulating the fracture behavior of rubber-like materials.
42 ENGINEERING↗
Advanced thermal/environmental barrier coatings of high-entropy rare earth disilicates tuned by strong anharmonicity of Eu 2 Si 2 O 7
Advancing thermal/environmental barrier coating (TEBC) materials with integrated thermal-mechanical functions is paramount for safeguarding SiC-based ceramic matrix composites (CMCs) in high-efficiency gas turbines. Herein, we employ a synergistic approach, combining density functional theory (DFT) methods and combinatorial chemistry techniques, to design high-performance and low-cost RE 2 Si 2 O 7 (RE = rare earth elements) TEBC materials tailored for enhanced compatibility with SiC-based CMCs. Expanding on phase stability of alloying pure RE 2 Si 2 O 7 , the investigation extends to the mechanical and thermal properties of solid solution systems, including Er 1/2 Y 3/4 Yb 3/4 Si 2 O 7 , Gd 1/4 Er 1/4 Y 3/4 Yb 3/4 Si 2 O 7 , and Eu 1/4 Er 1/4 Y 3/4 Yb 3/4 Si 2 O 7 . The solid solution systems exhibit a major reduction in lattice thermal conductivity relative to their pure counterparts, achieving ultralow values of 0.25 to 0.39 W m −1 K −1 at 1500 K. Furthermore, the coefficients of thermal expansion (CTE) of these solid solutions are precisely tuned within the desired range for SiC (4.4 to 5.5 × 10 −6 K −1 ), while maintaining good mechanical properties. Here, in particular, the addition of Eu 2 Si 2 O 7 demonstrates to be an important variable to the tuning of CTE and lattice thermal conductivity by leveraging its strong anharmonicity, presenting a pioneering avenue for fine-tuning material properties. In summary, this research not only identifies promising TEBC materials with superior thermal properties, but also introduces a valuable computational material design methodology for the rapid discovery of complex materials for harsh environments.
36 MATERIALS SCIENCE↗
Density functional modeling of the binding energies between aluminosilicate oligomers and different metal cations
Interactions between negatively charged aluminosilicate species and positively charged metal cations are critical to many important engineering processes and applications, including sustainable cements and aluminosilicate glasses. In an effort to probe these interactions, here we have calculated the pair-wise interaction energies (i.e., binding energies) between aluminosilicate dimer/trimer and 17 different metal cations M n+ (M n+ = Li + , Na + , K + , Cu + , Cu 2+ , Co 2+ , Zn 2+ , Ni 2+ , Mg 2+ , Ca 2+ , Ti 2+ , Fe 2+ , Fe 3+ , Co 3+ , Cr 3+ , Ti 4+ and Cr 6+ ) using a density functional theory (DFT) approach. Analysis of the DFT-optimized structural representations for the clusters (dimer/trimer + M n+ ) shows that their structural attributes (e.g., interatomic distances) are generally consistent with literature observations on aluminosilicate glasses. The DFT-derived binding energies are seen to vary considerably depending on the type of cations (i.e., charge and ionic radii) and aluminosilicate species (i.e., dimer or trimer). A survey of the literature reveals that the difference in the calculated binding energies between different M n+ can be used to explain many literature observations associated with the impact of metal cations on materials properties (e.g., glass corrosion, mineral dissolution, and ionic transport). Analysis of all the DFT-derived binding energies reveals that the correlation between these energy values and the ionic potential and field strength of the metal cations are well captured by 2nd order polynomial functions ( R 2 values of 0.99–1.00 are achieved for regressions). Given that the ionic potential and field strength of a given metal cation can be readily estimated using well-tabulated ionic radii available in the literature, these simple polynomial functions would enable rapid estimation of the binding energies of a much wider range of cations with the aluminosilicate dimer/trimer, providing guidance on the design and optimization of sustainable cements and aluminosilicate glasses and their associated applications. Finally, the limitations associated with using these simple model systems to model complex interactions are also discussed.
Gong, Kai↗
Design principles and physical properties of two-dimensional heterostructured borides
Principles of design to create dynamically stable transition metal-, lanthanide-, and actinide-based low-dimensional borides are presented. A charge transfer analysis of donor metal atoms to electron deficient honeycombed B lattices allows to predict complex covalent heterostructures hosting Dirac states. The applicable guidelines are supported with the analysis of phonon spectra computed with first-principles calculations to demonstrate the physical stability of nanometer-thick heterostructures. Similar or dissimilar layered borides can be stacked on top of each other in a layer-by-layer fashion creating compounds that can be fundamentally different from the individual layers, opening a rich playground to explore novel physical properties and new materials. Functionalities such as multiple Dirac states, highly dispersive electronic bands, and decoupled acoustic-optical phonon branches are studied. The combination of appealing electronic properties and physical realization makes of predicted layered borides promising materials to integrate a new generation of two-dimensional materials.
2-dimensional systems↗
Systematic Improvement of Quantum Monte Carlo Calculations in Transition Metal Oxides: sCI-Driven Wavefunction Optimization for Reliable Band Gap Prediction
Accurate determination of the electronic properties of correlated oxides remains a significant challenge for computational theory. Traditional Hubbard-corrected density functional theory (DFT+U) frequently encounters limitations in precisely capturing electron correlation, particularly in predicting band gaps. We introduce a systematic methodology to enhance the accuracy of diffusion Monte Carlo (DMC) simulations for both ground and excited states, focusing on LiCoO 2 as a case study. By employing a selected configuration interaction (sCI) approach, we demonstrate the capability to optimize wavefunctions beyond the constraints of single-reference DFT+U trial wavefunctions. Here, we show that the sCI framework enables accurate prediction of band gaps in LiCoO 2 , closely aligning with experimental values and substantially improving traditional computational methods. The study uncovers a nuanced mixed state of t 2g and e g orbitals at the band edges that is not captured by conventional single-reference methods, further elucidating the limitations of PBE+U in describing d-d excitations. Our findings advocate for the adoption of beyond-DFT methodologies, such as sCI, to capture the essential physics of excited-state wavefunctions in strongly correlated materials. The improved accuracy in band gap predictions and the ability to generate more reliable trial wavefunctions for DMC calculations underscore the potential of this approach for broader applications in the study of correlated oxides. This work not only provides a pathway for more accurate simulations of electronic structures in complex materials but also suggests a framework for future investigations of the excited states of other challenging systems.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Protection and Restoration Solutions to Reliable and Resilient Integration of Grid-connected PV Installations and Distributed Energy Resources: Design, Testbed, Proof of Work and Impact Studies (Final Report)
To gain a better understanding of the complex transients in a utility grid with a large number of solar PV installations coupled by PWM inverters during the protection and restoration period of the power grid, in this project a kW-level experimental system with dominating inverter-based resources and a hardware-in-loop simulation system with transmission and distribution models, as well as several SEL protective relays have been built and used. The major research findings of the project are summarized below in two perspectives: Experimental research: A self-organized ultra-high frequency solitary waveform is discovered and demonstrated in the hardware experimental system. Despite the familiarities, such a solitary waveform is distinct from any harmonics, transient, or resonance waves in that, it is 1) non-dispersive over time and space, 2) not associated with any active source or the linear superposition of sources, 3) half-cycle asymmetric, 4) not responding to filters or change of system characteristic resonance frequency, 5) ubiquitous as it occurs simultaneously everywhere in the system from DC supply, power lines, and the utility grid, and 6) explosive through tripping the protection or damaging susceptible devices or circuits. Analytical study: The nature of such a new waveform and an analytical explanation of its formation are studied based on related physics and non-linear science principles, with the following highlights: 1) such a new waveform follows the solutions to the Non-linear Schrödinger equation, so the classic linear perturbation theory is unable to explain or predict such a unique waveform as confirmed with the research team of RTDS; 2) the critical condition of occurrence of such a waveform is derived, which indicates that the breakings of solitary wave is a system synchronous issue between the utility grid in the 60Hz phasor domain and duty-ratio modulation of DC sources in the inverter switching frequency domain; 3) such a waveform carries energy mass so it could be detrimental; 4) such a waveform is deceiving as it is not readily detectable in the energy propagation direction by the primary protection equipment, while it is detrimental in the perpendicular direction of energy propagation, i.e. voltage direction. Thus, it has more likely challenges to the second primary equipment and devices, particularly at the weakest link and point such as aging insulation and inappropriate setting of susceptible devices. Therefore, such a waveform can be easily ignored in the aftermath investigation. The intellectual merits: The experimental discovery reveals certain unfamiliar transient phenomena that could challenge the integration of large-scale grid-connected solar PV installation during the group ride-through period of solar PV installations, commissioning of large-solar farms, or dramatic change of solar radiation conditions. Particularly, the research findings suggest that the occurrence of the unfamiliar transient phenomena is uniquely associated with the inverter-based solar PV installation, which is less likely to happen for rotary energy systems. The physics and non-linear science-based research work laid down an analytical path to the challenging transient stability problems including those that have been observed currently and future calls. Both experimental and analytical research results explain the limitation of many current research effects including some DoE research undertakings as well as possible solutions. The broader impacts: The research outcome of the project advances the understanding of the possible transient problems for the integration of inverter-based solar PV installation. It also highlights the theoretical and technical barriers for applying the conventional theory and technology to design and implement countermeasures against their adverse impacts of the large-scale solar PV installation and operation on grid reliability and security. More broadly, the research outcomes have shed some light on the open, fundamental challenges in the integration of large-scale solar photovoltaic energy into current and next-generation power grids across the country, and worldwide. The advanced physics and non-linear science-based analysis open up a path to harmonization of the renewables for societal energy needs via building a resilient and sustainable electric power infrastructure.
14 SOLAR ENERGY↗
Electronic structure analysis and DFT benchmarking of Rydberg-type alkali-metal-crown ether, -cryptand, and -adamanzane complexes
Density functional theory (DFT) and electron propagator theory (EPT) calculations were performed to study ground and excited electronic structures of alkali-metal (M) coordinated 9-crown-3, 24-crown-8, [2.1.1]cryptand, o-Me 2 -1.1.1, and 3 6 Adamanzane complexes. Each complex bears an expanded electron in the periphery and occupies diffuse 1p-, 1d-, 1f-type molecular orbitals (or superatomic 1P, 1D, 1F orbitals) in excited electronic states. The calculated superatomic shell model of the M(9-crown-3) 2 is 1S, 1P, 1D, 1F, 2S, 2P, 2D, 1G and it is held by all other complexes up to the studied 1F level. Due to the highly diffuse nature of the electron, the ionization energies of these complexes are significantly lower (1.6–2.0 eV) and hence these complexes belong to the superalkali category. Further, the ab initio EPT ionization energy and the excitation energies of the Li(9-crown-3) 2 were used to evaluate DFT errors associated with a series of exchange correlation functionals that span multiple rungs of Jacob's ladder (i.e., GGA, meta-GGA, global GGA hybrid, meta-GGA hybrid, range-separated hybrid, double-hybrid). Among these, the best performing functional is the range-separated hybrid CAM-B3LYP and the errors are within 6% of high-level ab initio EPT results. The accuracy of CAM-B3LYP is indeed transferable to similar complexes and hence the findings are expected to accelerate the progression of studies of Rydberg-type systems.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Using interfaces to: create strongly-coupled magnetic-ferroelectrics
Starting at the level of electrons and atoms our long-term goal is to rationally design complex oxide heterostructures and interface-materials with targeted emergent behaviors. This is not a matter of simply optimizing material parameters, but rather begins with understanding a mechanism to control the interplay between the diverse microscopic degrees of freedom prevalent in complex oxides in order to create targeted macroscopic phenomena and ends with the design of new material realizations. These realizations are in turn created with atomic-layer precision, structurally assessed to see that they are the intended realization, and finally their relevant properties are measured. During this program we have developed the scientific ideas necessary to apply this materials-by-design paradigm to the creation of materials offering electrical control of magnetism. At the heart of our approach lies a partnership between theory, synthesis, and characterization. These are the areas of expertise of the three co-PIs, who have worked closely together since their arrival at Cornell. During this project bold ideas on where to put the atoms to provide electrical control of magnetism were met with the latest methods of synthesis science and electron microscopy in an attempt create targeted new multiferroics. Our work applying this materials-by-design methodology to multiferroics has enjoyed multiple successes. The Co-PIs have invented new mechanisms for multiferroics, which provide strong coupling between polarization and magnetism, and reduced them to practice. Specifically, we pioneered multiferroics based on spin-phonon coupling, rotation-driven multiferroicity, and novel geometric multiferroics. Importantly, the latter two of these mechanisms enable the deterministic switching of magnetism by an applied electric field, which is key to the application of multiferroics to future device technologies. Most recently, we performed detailed theoretical analysis on a system involving the combination of magnetoelectric and colossal magnetoresistance behaviors and predict that electrical control of a metal-insulator transition is possible in appropriately strained SmBaMn2O6. The research enabled by this grant led to 63 publications in leading refereed journals including Nature, Nature Materials, Physical Review Letters, and Advanced Materials demonstrating a materials-by-design paradigm to the creation of materials offering electrical control of magnetism.
36 MATERIALS SCIENCE↗
Quantifying Errors in Effective Cluster Interactions of Lattice Gas Cluster Expansions
The promise of lattice gas (LG) cluster expansions (CEs) is that they can describe a given system property to any level of accuracy since the orthogonal “cluster basis functions” span the complete space of available configurations. Such an approach can be constructed to an arbitrarily large surface of a finite number of distinct adsorption sites. Unfortunately, this is only true for the case of an ideal, fixed lattice decorated with components at precise lattice points (the lattice “sites”) with no distortions or relaxations subsequently allowed. Since most systems, and surfaces specifically, do not conform to such an ideal set of constraints, errors in LG CEs must be expected or CE convergence severely hampered. Beyond this, numerical errors in the provided data can complicate the proper construction of a truly predictive and/or physically significant CE. In this work, we show here how reliance on typical statistical tools like confidence intervals cannot be expected to provide an accurate representation of the uncertainty of effective cluster interactions (ECIs) in the CE due to the nature of the target ab initio data and the nature of CEs themselves. We develop a method for estimating these errors that does not rely on statistical assumptions about the model or data. We then use these ECI errors to quantify fundamental consequences on the uncertainty of ECIs in CEs built from O/Fe(100) data whose surface and adsorbates have been allowed to relax in a typical manner and from O/Fe(100) data whose surface and adsorbates are fixed in ideal lattice positions. We also quantify the effect of using a different density functional theory exchange–correlation functional, using these ECI errors to assess the significance in any deviations. In both cases, our method is shown to have remarkable utility in the quantification of errors in the ECIs of CEs. While we stick to the lattice gas convention in this work, the method is equally applicable to the Ising convention or, in principle, any linear model of sufficient complexity.
08 HYDROGEN↗
Symmetry-specific orientational order parameters for complex structures
A comprehensive framework of characterizing complex self-assembled structures with a set of orientational order parameters is presented. It is especially relevant in the context of using anisotropic building blocks with various symmetries. Two classes of tensor order parameters are associated with polyhedral nematic and bond orientational order. For the latter, a variation of classical bond order parameters (BOPs) is introduced, which takes advantage of the symmetry of constituent particles and/or expected crystalline phases. These symmetrized BOPs can be averaged over an entire system or assigned locally to an individual bond. By combining that with the bond percolation procedure, one is able to identify coherent domains within a self-assembled structure. As a demonstration of the proposed framework, we apply it to a simulated hybrid system that combines isotropic and patchy particles with octahedral symmetry. In conclusion, not only does the methodology allow one to identify individual crystalline domains but also it detects coherent clusters of a peculiar compact amorphous structure that is not space-filling and lacks any long-range order.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Investigating trends in actinide covalency and magnetism with 35/37 Cl SSNMR spectroscopy and first-principles calculations
Solid-state NMR (SSNMR) spectroscopy is a powerful technique for studying actinide chemistry but has been significantly limited due to the complex paramagnetism, and radiological hazards presented by these materials. Lanthanide and actinide salts often feature magnetic ordering and can be paramagnetic, ferromagnetic, or antiferromagnetic depending on temperature and electronic structure. Paramagnetic interactions can manifest in SSNMR both as secular spectral shifts and/or couplings as well as contributions from non-secular relaxation. Both effects can be directly measured with NMR and used to extrapolate rich chemical information such as coordination environments, bonding characteristics, local molecular dynamics, and correlation times. Typically, these studies are carried out on high-γ and highly abundant NMR-active isotopes (e.g., 1 H, 6/7 Li, 19 F, 23 Na, etc.) or on enriched rare isotopes (e.g., 2 H and 17 O), which can be expensive. Herein, we present a facile methodology to measure the 35/37 Cl electric-field gradient (EFG) and paramagnetic shift anisotropy (SA) tensor components using static wideline SSNMR measurements of LaCl 3 , NdCl 3 , UCl 3 , and UCl 4 . The static powder spectra were measured with both 35 Cl and 37 Cl SSNMR to increase the fidelity of the extracted tensor parameters. Variable temperature NMR of a select case confirms the Curie–Weiss paramagnetism. Relaxation measurements of both nuclei further corroborate observations owing to the paramagnetic relaxation enhancement and reveal simultaneous quadrupolar relaxation mechanisms. Density functional theory (DFT) calculations using Hubbard U corrections to the uranium valence orbitals show excellent agreement with experimental EFG tensor parameters and help describe the bonding characteristics in these lanthanide and actinide systems.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Measurement of the angular dependence of the spontaneous Raman scattering in anisotropic crystalline materials using spherical samples: Potassium dihydrogen phosphate as a case example
A specialized experimental configuration was created to allow for more-accurate characterization of the spontaneous Raman scattering properties in anisotropic materials. This need stems from the challenges, arising from the complexity of light propagation, in obtaining accurate measurements of the angular dependence of the Raman scattering cross section in birefringent materials. The nonlinear optical material KH 2 PO 4 (KDP) is used as a model medium. Our report is motivated by the need to improve our understanding and management of transverse stimulated Raman scattering in KDP crystals and its deuterated analog, DKDP, typically used for frequency conversion and polarization control in large-aperture laser systems. Key to this experimental platform is the use of high-quality spherical samples that enable one to measure the Raman scattering cross section in a wide range of geometries using only a single sample. The effect of polarization rotation of both the pump light and the collected Raman signal must be carefully considered in data analysis and can give rise to artifacts, which can, in part, be mitigated by reducing the input and collection cone angles.
47 OTHER INSTRUMENTATION↗
Role of orbital selectivity on crystal structures and electronic states in BiMnO 3 and LaMnO 3 perovskites
Correlated oxides, such as BiMnO 3 and LaMnO 3 , show a complex interplay of electronic correlations and crystal structure exhibiting multiple first-order phase transitions, some without a clear order parameter. The quantitative theoretical description of this temperature-dependent electronic-structural interplay in the vicinity of a Mott transition is still a challenge. Here we address this issue by simultaneously considering both structural and electronic degrees of freedom, within a self-consistent density functional theory with embedded dynamical mean field theory. Our results show the existence of novel electronic states characterized by the coexistence of insulating, semimetallic, and metallic orbitals. This state is in agreement with resonant x-ray scattering. Here we also show that electronic entropy plays a decisive role in both electronic and structural phase transitions. By self-consistent determination of both the electronic state and the corresponding crystal structure, we show that the temperature evolution of these phases can be quantitatively explained from first principles, thus demonstrating the predictive power of the theoretical method for both the structural and the electronic properties.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Unified understanding of the impact of semiflexibility, concentration, and molecular weight on macromolecular-scale ring diffusion
Conformationally fluctuating, globally compact macromolecules such as polymeric rings, single-chain nanoparticles, microgels, and many-arm stars display complex dynamic behaviors due to their rich topological structure and intermolecular organization. Synthetic rings are hybrid objects with conformations that display both ideal random walk and compact globular features, which can serve as models of genomic DNA. To date, emphasis has been placed on the effect of ring molecular weight on their unusual behaviors. Here, we combine simulations and a microscopic force-level theory to build a unified understanding for how key aspects of ring dynamics depend on different tunable molecular properties including backbone rigidity, monomer concentration, degree of traditional entanglement, and molecular weight. Our large-scale molecular dynamics simulations of ring melts with very different backbone stiffnesses reveal unanticipated behaviors which agree well with our generalized theory. This includes a universal master curve for center-of-mass diffusion constants as a function of molecular weight scaled by a chemistry and thermodynamic state-dependent critical molecular weight that generalizes the concept of an entanglement cross-over for linear chains. The key physics is how backbone rigidity and monomer concentration induced changes of the entanglement length, interring packing, degree of interpenetration, and liquid compressibility slow down space-time dynamic-force correlations on macromolecular scales. A power law decay of the center-of-mass diffusion constant with inverse molecular weight squared is the first consequence, followed by an ultraslow activated hopping transport regime. Our results set the stage to address slow dynamics and kinetic arrest in different families of compact synthetic and biological polymeric systems.
Science & Technology - Other Topics↗
Tuning magnetic antiskyrmion stability in tetragonal inverse Heusler alloys
The identification of materials supporting complex, tunable magnetic order at ambient temperatures is foundational to the development of new magnetic device architectures. In this work, we report the design of Mn 2 XY tetragonal inverse Heusler alloys that are capable of hosting magnetic antiskyrmions whose stability is sensitive to elastic strain. Here We first construct a universal magnetic Hamiltonian capturing the short- and long-range magnetic order which can be expected in these materials. This model reveals critical combinations of magnetic interactions that are necessary to approach a magnetic phase boundary, where the magnetic structure is highly susceptible to small perturbations such as elastic strain. We then computationally search for quaternary Mn 2 (X 1 , X 2 ) Y alloys where these critical interactions may be realized and which are likely to be synthesizable in the inverse Heusler structure. We identify the Mn 2 Pt 1–z X z Ga family of materials with X = Au, Ir, Ni as an ideal system for accessing all possible magnetic phases, with several critical compositions where magnetic phase transitions may be actuated mechanically.
36 MATERIALS SCIENCE↗
235-F GoldSim Fate and Transport Model: Uncertainty Quantification
Building 235-F was configured with two missions in mind: Actinide Billet Line (ABL) and the fabrication of Pu-238 oxide for space program applications. ABL produced Np-237 billets for use in SRS reactors, whereas the design process, fabrication, and examination of Pu-238 oxide powder occurred within the following areas, respectively: Plutonium Experimental Facility (PEF), Plutonium Fuel Form (PuFF), and Old Metallography Lab (OML). By 1990 production ceased and by 2006 de-inventory occurred; however, assays have shown significant holdup remains within ABL and PuFF. As a result, 235-F is a Category 2 nuclear facility, with plans to undergo deactivation and decommission (D and D) via In-Situ Disposal (ISD). The purpose of this project is to ensure United States Environmental Protection Agency (USEPA) groundwater radiation maximum contaminant level (MCL) and dosage standards are met during the D and D of 235-F by quantifying uncertainty through probabilistic modeling and evaluation of various ISD alternatives. GoldSim is a dynamic modeling software package with a graphical, object-oriented interface capable of capturing the influence of complex system input variability on probabilistic system outcomes. A GoldSim stochastic fate and transport model for 235-F was developed and matched with a PORFLOW deterministic model to simulate probabilistic release and flow of radionuclides from ABL and PuFF into the vadose zone, the Upper Three Runs (UTR) Aquifer, and UTR Creek. The GoldSim model was used to probabilistically evaluate four ISD scenarios against USEPA groundwater radiation MCLs and dosage standards. The deterministic 235-F GoldSim fate and transport model continues to be refined to match the results of the PORFLOW deterministic model to ensure the model accurately represents radionuclide movement through the groundwater system. The stochastic variables that are utilized within the GoldSim model are founded on the most current data; a conservative perspective is taken where needed. Alignment with the PORFLOW deterministic model, coupled with input stochastic variability, allows the probabilistic 235-F GoldSim model to capture the conservative breadth of possible outcomes for radionuclide fate within this particular system. This ensures that the USEPA MCLs and dosage limits hold even in the worst case scenarios.
12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
Theoretical studies of chemical reactions related to the formation and growth of polycyclic aromatic hydrocarbons (PAH) and molecular properties of their key intermediates (Final Progress Report)
The formation mechanisms of polycyclic aromatic hydrocarbons, (PAHs) – organic molecules carrying fused benzene rings – are of great interest to scientists and engineers due to their importance in combustion chemistry and astrochemistry. On Earth, PAHs are largely produced in incomplete combustion of fossil fuel and are considered as critical precursors to unwanted soot particles leading to combustion inefficiency and causing air pollution along with detrimental health effects. Simple PAH molecules initially formed in the gas phase, are further involved in a build-up process in combustion flames leading to larger PAH, bowl-shaped nanostructures, fullerenes, and solid-phase species including carbonaceous dust, graphene particles, and soot. In deep space, PAH and their derivatives are potential key intermediates and nucleation sites leading eventually to carbonaceous nanoparticles (“interstellar grains”). Therefore, the understanding of the key processes in the synthesis of PAHs along with their precursors and their degradation mechanisms in combustion systems and in interstellar, circumstellar, and planetary atmospheric environments will provide critical insights into how complex aromatic structures, carbonaceous nanoparticles, and fullerenes are formed and destroyed. Achieving this understanding is an important step in the development of the efficient combustion processes and of the ecofriendly devices with reduced environmental pollution as well as technological strategies for the production of hydrogen and solid carbon through thermal or plasma-assisted pyrolysis of natural gas and biomass. Also, the understanding of the key processes of PAH and soot growth will help in our comprehension of chemical evolution in the universe. Detailed information on the mechanisms and reliable rate constants of the key elementary chemical reactions involved in PAH formation and destruction processes and in inception of soot particles is often missing, with the main deficiencies being the absence of temperature- and pressure-dependent rate constants for the broad range of conditions occurring in various terrestrial and interstellar processes and the lack of data on the reaction products and their branching ratios. Complementary to experimental studies, these gaps in knowledge can be filled by using quantum chemical calculations of reaction potential energy surfaces providing us with accurate energies of reaction products, intermediates, and transition states, revealing the reaction mechanism, and giving the molecular properties required to compute rate constants for relevant reaction steps and product branching ratios using the RRKM-Master Equation (ME) method. Molecular dynamics (MD) simulations can be used in cases when a reaction rate cannot be properly described by statistical theories. During the terminal renewal project period we employed these ab initio/RRKM-ME and MD approaches to complete our studies on several key reactions relevant to the formation/growth of PAH and inception of soot particles including (1) the reaction mechanism and kinetics of the resonance stabilized fulvenallenyl radical with propargyl and C 3 H 4 isomers; (2) the reaction mechanism and kinetics for the C + indene and C 2 + styrene reactions producing naphthyl or azulenyl radicals in low-temperature environments; (3) the MD study of non-equilibrium dimerization of acepyrene and coronene and its radical. The information derived from our theoretical calculations contributed to a better fundamental understanding of the reaction mechanisms and provide missing critical kinetic data to improve combustion models of hydrocarbon fuels and astrochemical models of the growth of carbonaceous molecules and particles in cold molecular clouds, circumstellar envelopes, and planetary atmospheres.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗