Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Complex chemical systems”

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

Synergistic Solvent-Surface Interactions Enable Alkyne Semihydrogenation at Palladium

Enabling higher yield and better selectivity for fine-chemical synthesis through heterogeneous catalysis is intricately linked to the interplay of active sites, reaction conditions, and mass transfer influence provided by the catalyst. Alkyne semihydrogenation is ubiquitous in the production of bulk chemicals in the pharmaceutical, polymer, or fine-chemical industries, but product selectivity remains a major challenge. Here, in this study, we demonstrate that the design of catalysts encompassing nickel (Ni) foams as contiguous monolith supports, decorated with ultralow loading of Pd/PdO x nanoparticles on a carbonized polydopamine interface and tuned with a thin layer of Al 2 O 3 , in conjunction with an optimized reaction environment leads to highly selective alkyne semihydrogenation. The reactions demonstrate good functional group tolerance and applicability to flow reactor systems. Combined computational and experimental studies are presented to describe the synergistic effect between the solvent-surface interaction and the degree of Pd surface reduction that are necessary to promote this selectivity. The system highlights the opportunity for catalyst-solvent codesign as a benign alternative to more complex reactants featuring extrinsic poisons or less-favored dopants.

atomic layer deposition↗

Crosslinked polymeric ionic liquid stationary phases incorporating a combination of silver(I) and copper(II) ions for olefin separations by gas chromatography

The development of efficient and selective olefin separation technologies is a critical challenge within the chemical industry given the high energy demands of conventional cryogenic distillation methods. To address this, the study explores the effectiveness of polymeric ionic liquid (PIL) stationary phases in gas chromatography (GC) that incorporate both silver(I) ([Ag + ]) and copper(II) ([Cu 2+ ]) ions for enhanced olefin separation via π-complexation interactions. Further, the effects of both non-crosslinked and crosslinked PIL stationary phases containing varying ratios of [Ag⁺] and [Cu²⁺] ions were systematically investigated to evaluate their impact on the separation of olefins from paraffins and structurally similar olefins. Retention measurements using a broad range of olefinic probes revealed that [Ag⁺] ions provided stronger and more consistent π-complexation in PIL-based stationary phases, while [Cu²⁺] ions preferentially enhance olefin retention in crosslinked systems. Furthermore, thermal stability studies showed that [Ag + ]-containing PIL stationary phases maintain stable olefin separation efficiency over prolonged exposure under elevated temperatures, while crosslinked PIL stationary phases containing [Cu 2+ ] ions exhibited stronger retention due to increased surface adsorption. These findings provide insight into the tunability of metal–olefin interactions using PIL-based stationary phases and suggest design strategies for achieving thermally stable and highly selective chromatographic separations of unsaturated hydrocarbons.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Supervised learning of a chemistry functional with damped dispersion

Abstract Kohn–Sham density functional theory is widely used in chemistry, but no functional can accurately predict the whole range of chemical properties, although recent progress by some doubly hybrid functionals comes close. Here, we optimized a singly hybrid functional called CF22D with higher across-the-board accuracy for chemistry than most of the existing non-doubly hybrid functionals by using a flexible functional form that combines a global hybrid meta-nonseparable gradient approximation that depends on density and occupied orbitals with a damped dispersion term that depends on geometry. We optimized this energy functional by using a large database and performance-triggered iterative supervised training. We combined several databases to create a very large, combined database whose use demonstrated the good performance of CF22D on barrier heights, isomerization energies, thermochemistry, noncovalent interactions, radical and nonradical chemistry, small and large systems, simple and complex systems and transition-metal chemistry.

Liu, Yiwei (ORCID:0000000288126163)↗

Decoupling Complex Multi–Length–Scale Morphology in Non–Fullerene Photovoltaics with Nitrogen K–Edge Resonant Soft X–ray Scattering

Complex morphology in organic photovoltaics (OPVs) and other functional soft materials commonly dictates performance. Such complexity in OPVs originates from the mesoscale kinetically trapped non-equilibrium state, which governs device charge generation and transport. Resonant soft X-ray scattering (RSoXS) has been revolutionary in the exploration of OPV morphology in the past decade due to its chemical and orientation sensitivity. However, for non-fullerene OPVs, RSoXS analysis near the carbon K-edge is challenging, due to the chemical similarity of the materials used in active layers. Here, an innovative approach is provided by nitrogen K-edge RSoXS (NK-RSoXS), utilizing the spatial and orientational contrasts from the cyano groups in the acceptor materials, which allows for determination of phase separation. NK-RSoXS clearly visualizes the combined feature sizes in PM6:Y6 blends from crystallization and liquid-liquid demixing, while PM6:Y6:Y6-BO ternary blends with reduced phase-separation size and enhanced material crystallization can lead to current amplification in devices. Nitrogen is common in organic semiconductors and other soft materials, and the strong and directional N 1s → π* resonances make NK-RSoXS a powerful tool to uncover the mesoscale complexity and open opportunities to understand heterogeneous systems.

36 MATERIALS SCIENCE↗

DECOVALEX-2019: An international collaboration for advancing the understanding and modeling of coupled thermo-hydro-mechanical-chemical (THMC) processes in geological systems

The DECOVALEX Project is an international research collaboration for advancing the understanding and modeling of coupled thermo-hydro-mechanical-chemical (THMC) processes in geological systems. DECOVALEX stands for “DEvelopment of COupled Models and VALidation against EXperiments”. The creation of this international initiative, now running for almost 30 years, was motivated by the recognition that prediction of these coupled effects is an essential part of the performance and safety assessment of geologic disposal systems for radioactive waste and spent nuclear fuel, and also for a range of other sub-surface engineering activities. DECOVALEX emphasizes joint analysis and comparative modeling of state-of-the-art field and laboratory experiments, across a range of host rock options and repository designs. Participating research teams are from radioactive waste management organizations, national research institutes, regulatory agencies, universities, and consulting groups, providing a wide range of perspectives and solutions to these complex problems. The most recent phase of the DECOVALEX Project, here referred to as DECOVALEX-2019, started in 2016 and ended in 2019. Modeling teams from 13 international partner organizations participated in the comparative evaluation of seven modeling tasks involving complex field and/or laboratory experiments. Furthermore, this Virtual Special Issue on DECOVALEX-2019 provides an in-depth overview of these collaborative research efforts and how these have advanced the state-of-the-art of understanding and modeling coupled THMC processes.

58 GEOSCIENCES↗

Using a Chromatographic Pseudophase Model To Elucidate the Mechanism of Olefin Separation by Silver(I) Ions in Ionic Liquids

Silver(I) ions undergo selective olefin complexation and have been utilized in various olefin/paraffin separation techniques such as argentation chromatography and facilitated transport membranes. Ionic liquids (ILs) are solvents known for their low vapor pressure, high thermal stability, low melting points, and ability to promote a favorable solvation environment for silver(I) ion–olefin interactions. To develop highly selective separation systems, a fundamental understanding of analyte partitioning to the stationary phase and the thermodynamic driving forces behind solvation is required. In this study, a chromatographic model treating silver(I) ions as a pseudophase is constructed and employed for the first time to investigate the olefin separation mechanism in silver(I) salt/IL mixtures. Stationary phases containing varying amounts of noncoordinated silver(I) salt ([Ag + ][NTf 2 – ]) dissolved in the 1-decyl-3-methylimidazolium bis[(trifluoromethyl)sulfonyl]imide ([C 10 MIM + ][NTf 2 – ]) IL are utilized to determine the partition coefficients of various analytes including alkanes, alkenes, alkynes, aromatics, aldehyde, esters, and ketones. As ligand coordination to silver(I) ions is known to lower its olefin complexation capability, this study also examines two different types of coordinated silver(I) ion pseudophases, namely, monocoordinated silver(I) salt ([Ag + (1-decyl-2-methylimidazole, DMIM)][NTf 2 – ]) and dicoordinated silver(I) salt ([Ag + (1-methylimidazole, MIM)(DMIM)][NTf 2 – ]). The extent of olefin partitioning to the coordinated silver(I) ion pseudophases over the carrier gas and IL decreased by up to two orders of magnitude. Values for enthalpy, entropy, and free energy of solvation were determined for the three silver(I) ion-containing systems. Olefin retention was observed to be enthalpically dominated, while ligand coordination to the silver(I) ion pseudophase resulted in variations for both enthalpic and entropic contributions to the free energy of solvation. Furthermore, the developed model can be used to study chemical changes that occur in silver(I) ions over time as well as identify optimal silver(I) salt/IL mixtures that yield high olefin selectivity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Low Yield Nuclear Monitoring Physics Experiment 1 – Integrated Data Acquisition System Design and Initial Observations

The report documents the design of the Integrated Data AcQuisition (IDAQ) system and observations recorded during the first in a series of underground chemical explosions conducted on the Nevada National Security Site (NNSS) in southern Nevada. Experiments are funded as part of Low Yield Nuclear Monitoring (LYNM) research and development within the United States National Nuclear Security Administration NA-22 nuclear non-proliferation program. The series is part of the broader Physical Experiment 1 (PE1) being conducted in and around the P-tunnel facility on the NNSS. Each explosive experiment utilizes several tons of comp-B to generate signals recorded by a broad suite of instrumentation. The IDAQ serves as the backbone for all subsurface instrumentation providing precise time synchronization, remote control, data exfiltration and backup, along with recording several sensing modalities throughout the underground complex that includes ground motion, environmental conditions, and electromagnetic signals.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Rheology guiding the design and printability of aqueous colloidal composites for additive manufacturing

Vat photopolymerization (VP) and direct ink write (DIW) additive manufacturing (AM) provide complex geometries with precise spatial control employing a vast array of photo-reactive polymeric systems. Although VP is recognized for superior resolution and surface finish, DIW provides versatility for higher viscosity systems. However, each AM platform presents specific rheological requirements that are essential for successful 3D printing. First, viscosity requirements constrain VP polymeric materials to viscosities below 10 Pa s. Thus, this requirement presents a challenging paradox that must be overcome to attain the physical performance of high molecular weight polymers while maintaining suitable viscosities for VP polymeric materials. Second, the necessary rheological complexity that is required for DIW pastes requires additional rheological measurements to ensure desirable thixotropic behavior. Further, this manuscript describes the importance of rheological measurements when designing polymeric latexes for AM. Latexes effectively decouple the dependency of viscosity on molecular weight, thus enabling high molecular weight polymers with low viscosities. Photo-crosslinking of water-soluble monomers and telechelic oligomeric diacrylates in the presence of the latex enables the fabrication of a scaffold, which is restricted to the continuous aqueous phase and effectively surrounds the latex nanoparticles enabling the printing of otherwise inaccessible high molecular weight polymers. Rheological testing, including both steady and oscillatory shear experiments, provides insights into system properties and provides predictability for successful printing. This perspective article aims to provide an understanding of both chemical functionality (photo- and thermal-reactivity) and rheological response and their importance for the successful design and evaluation of VP and DIW processable latex formulations.

36 MATERIALS SCIENCE↗

Ultrafast Electron Transfer in Au-Cyanobacteria Hybrid for Solar to Chemical Production

The rise of inorganic-biohybrid organisms for solar-to-chemical production has spurred mechanistic investigations into the dynamics of the biotic-abiotic interface to drive the development of next-generation hybrid systems. The model system, cyanobacteria-gold nanoparticle hybrids, combines a light harvester with a photosynthetic bacterium to drive the reduction of CO2 to glycerol with improved efficiency and increased glycerol production by 14.6%, in comparison to cyanobacteria only. In this work, we report insights into this unique photochemical behavior and propose a charge-transfer pathway from Au nanoparticle to cyanobacteria. Transient absorption (TA) spectroscopy revealed that photoexcited electron transfer rates are on the order of a few ps to the potential electron acceptor in photosystem II. This work represents a promising platform to utilize a conventional spectroscopic methodology to extract insights from more complex biotic-abiotic hybrid systems.

charge transfer↗

Tank Waste LDR Organics Data Summary for Sample-and-Send (Rev.1A)

The presence of organic chemicals regulated under the Resource Conservation and Recovery Act (RCRA) Land Disposal Restrictions (LDR) adds complexity to treating and disposing of the low activity fraction of Hanford tank waste if a low temperature treatment method such as grouting is used (SRNL-STI-2020-00228). The complexity arises from the fact that the baseline vitrification method is considered by the Washington State Department of Ecology (Ecology) as providing adequate thermal treatment for organics; a status not automatically extended to a lowtemperature process, such as solidifying the waste in a cementitious waste form. In addition, the Environmental Protection Agency (EPA) LDR program is intended to ensure that wastes are properly treated prior to disposal. Proper treatment makes hazardous waste less harmful to groundwater by reducing the mobility and/or toxicity of the hazardous constituents in the waste. EPA guidance indicates that stabilization/solidification of waste for organics could be considered impermissible dilution under the LDR dilution prohibition. In addition, waste storage activities at Hanford have required transferring and blending waste within the tank system and these activities have potentially altered the concentrations of the hazardous constituents. The LDR dilution prohibition found in 40 Code of Federal Regulations (CFR) 268.3 states that “… no generator, transporter, handler, or owner or operator of a treatment, storage, or disposal facility shall in any way dilute a restricted waste or the residual from treatment of a restricted waste as a substitute for adequate treatment …”. Hence, if LAW is to be treated using low temperature stabilization (such as cementation), then it is important to demonstrate both how past storage activities have contributed to the removal (by vacuum evaporation), or destruction (by in situ decomposition) of the LDR organics and how future retrieval and waste feed preparation will contribute to their removal (by filtration and ion exchange). Demonstrating these processes helps validate that cementation without additional organic treatment does not necessarily represent impermissible dilution. To aid in implementing cementitious solidification of Low Activity Waste (LAW), WRPS has been developing a regulatory and processing LDR treatment variance strategy termed “Sampleand-Send” that relies, in part, on demonstrating that in situ decomposition reactions along with historic evaporation of tank waste has destroyed or removed most of the LDR organics possibly associated with Hanford Tank Waste (SRNL-STI-2020-00582, SRNL-STI-2021-00453, SRNL-STI-2022-00391). Under the Sample-and-Send concept, Hanford tank waste would be retrieved, processed through a Tank-Side Cesium Removal-like system, and staged as a candidate feed that would then be sampled to confirm the waste acceptance criteria is met for solidification in an LAW cementitious treatment facility. If it can be shown that LDR organics are at concentrations below the waste acceptance criteria (WAC) for cementitious stabilization and have been sufficiently removed (by historic evaporation or by filtration and ion exchange during Cs removal), destroyed (by historic in situ decomposition), or are not soluble in LAW above the WAC then additional organic treatment is not needed prior to creating a cementitious final waste form and the concept of Sample-and-Send would be proposed to establish a non-rulemaking site-specific treatment variance using the specified method of treatment “STABL” to remove sampling requirements of the waste form after treatment. Waste not meeting the WAC could either be routed to the Hanford Waste Treatment and Immobilization Plant for LAW vitrification, or further processed by evaporation or chemical oxidation before solidifying in a cementitious waste form. A key component in implementing the Sample-and-Send strategy is identifying which of the 207 LDR organic compounds associated with the RCRA Part A permit application waste codes for the Double Shell Tanks (DSTs) and Single Shell Tanks (SSTs) and any applicable Underlying Hazardous Constituents (UHCs) from 40 CFR 268.48 should be considered as potentially present and thus subject to regulation. In addition, it is also necessary to understand the solubility volatility, and reactivity of these compounds in LAW to identify which of the potentially present LDR organic compounds are not soluble above regulatory levels or are likely to have been removed by historic evaporation or destroyed by in situ decomposition reactions. If there are potentially present LDR organic compounds that have not been removed or destroyed and are soluble above regulatorily significant concentrations then a treatability variance may be needed for these species to eliminate any concerns pertaining to impermissible dilution. The spreadsheet accompanying this calculation report contains the data and logic computations needed to screen the list of 207 LDR organics associated with Hanford tank waste to identify those potentially present and to indicate which compounds may need to be included in a treatability variance.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Vacancy-mediated complex phase selection in high entropy alloys

Phase selection in Ti-Zr-Hf-Al high-entropy alloys was investigated by in-situ high-energy x-ray diffraction, single-crystal x-ray diffraction, and KKR-CPA electronic-structure methods that address disorder and vacancies, predicting formation enthalpy and chemical short-range order (SRO). Samples with varying Al content were synthesized, characterized, and computationally assessed to ascertain the composition-dependent phase selection, as increased Al content often acts as a stabilizer of a body-centered-cubic structure. Equiatomic TiZrHfAl was especially interesting due to its observed bcc superstructure – a variant of γ-brass with 4 vacancies per cell (not 2 as in γ-brass). We highlight how vacancy ordering mediates selection of this variant of γ-brass, which is driven by vacancy-atom SRO that dramatically suppress all atomic SRO. As vacancies are inherent in processing refractory systems, we expect that similar discoveries await in other high entropy alloys or in revisiting older experimental data.

36 MATERIALS SCIENCE↗

Multi-Scale Modeling Framework for Mercury Biogeochemistry

Multi-Scale modeling of mercury (Hg) geochemical speciation and reactions has been performed by integrating atomistic quantum chemical calculations with continuum scale speciation models. Major progress has been made in the improvement of quantum chemical models to calculate critical thermodynamic data for Hg complexes in aquatic environments. Rapid and reliable quantum chemical approaches have been developed for calculating acid dissociation constants (pK a ) and stability constants (log K), with calculated mean unsigned errors of 0.5 and 1.5 log units, respectively for ligand molecules and Hg complexes. At the continuum scale, systematic analysis of uncertainty propagation in mercury (Hg) speciation modeling has been conducted and was used to identify environmental conditions under which thermodynamic constant uncertainties are significant and recommended to be accounted for. The integrated framework for multi-scale modeling of mercury geochemistry is open to the research community through the web-based multiscale modeling aqueous speciation resource, AQUA-MER. The improved quantum chemical approaches for thermodynamic constant calculations are accessible through AQUA-MER and can be used to provide the missing constants in the continuum scale speciation calculations. In addition to low molecular mass Hg complex speciation, modeling natural aquatic environments also involve the transport of high molecular weight dissolved organic matter (DOM) in reactive flows simultaneously with equilibrium and kinetic reactions. To this end, atomistic MD simulations were performed to capture the details of aggregation, mechanisms and distribution of functional groups in DOM at the molecular level. The elemental composition and calculated bulk properties of the DOM models are in close agreement with experimental measurements. A travel-time based reactive transport model in the hyporheic zone of stream corridors was established for the multicomponent Hg-DOM-S system and implemented through PFLOTRAN.

54 ENVIRONMENTAL SCIENCES↗

Building a predictive model for polycyclic aromatic hydrocarbon dosimetry in organotypically cultured human bronchial epithelial cells using benzo[ a ]pyrene

The airway epithelium is a primary route of exposure for inhaled toxicants, and organotypic culture models represent an important advancement for toxicity testing compared to simple in vitro models that may lack metabolic capability and multicellular structure/communication associated with the bronchial epithelium in vivo. A quantitative understanding of chemical dosimetry is key for interpreting and extrapolating study results; however, dosimetry is understudied in organotypic models limiting ability to predict toxicity. We developed a dosimetry model for primary human bronchial epithelial cells (HBECs) cultured at the air-liquid interface (ALI) using benzo[a]pyrene (BAP), a representative polycyclic aromatic hydrocarbon. Dose and time course evaluation of metabolite formation and enzyme activity and expression were utilized to parameterize a cellular dosimetry model to improve the utility of ALI-HBECs for assessing chemical risk. Dosimetry analysis demonstrated absorption of BAP into cells and an increase in Phase 1 and 2 metabolites over time that correlated with regulation of metabolizing enzymes. BAP was cleared from cells by 48 h after exposure, and the primary metabolites generated in ALI-HBECs were BAP-3-phenol, BAP-4,5-dihydrodiol, BAP-7,8-dihydrodiol, BAP-9,10-dihydrodiol, BAP-7,8,9,10-tetrol, BAP-3-phenol-glucuronide, BAP-4,5-dihydrodiol-glucuronide, and BAP-9,10-dihydrodiol-glucuronide. The resulting dosimetry model described BAP and 7,8-dihydrodiol toxicokinetics in ALI-HBECs and suggested active excretion of 7,8-dihydrodiol. Overall, this study demonstrates metabolic competency of ALI-HBECs for BAP metabolism, demonstrates the usefulness of complex in vitro systems for human-relevant toxicity data, and exhibits how in silico models can be utilized for understanding the dosimetry of test compounds to aid in in vitro to human extrapolation of toxicity data for risk assessments.

Benzo[a]pyrene↗

A Geometric Measure Theory Approach to Identify Complex Structural Features on Soft Matter Surfaces

The structural features that protrude above or below a soft matter interface are well-known to be related to interfacially mediated chemical reactivity and transport processes. It is a challenge to develop a robust algorithm for identifying these organized surface structures, as the morphology can be highly varied and they may exist on top of an interface containing significant interfacial roughness. A new algorithm that employs concepts from geometric measure theory, algebraic topology, and optimization is developed to identify candidate structures at a soft matter surface, and then, using a probabilistic approach, to rank their likelihood of being a complex structural feature. The algorithm is tested for a surfactant laden water/oil interface, where it is robust to identifying protrusions responsible for water transport against a set identified by visual inspection. To our knowledge, this is the first example of applying geometric measure theory to analyze the properties of a chemical/materials science system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Reactive Transport Modeling of Mineral Precipitation and Carbon Trapping in Discrete Fracture Networks

In this study we use numerical experiments to analyze reactive flow and transport behavior in discrete intersecting fracture networks, focusing on (a) how reaction-induced changes in physical and chemical properties affect flow connectivity and (b) how fracture networks developed in the Earth's critical zone contribute to carbon sequestration via mineral weathering reactions. In the first part of the study, we used two-dimensional reactive flow and transport simulations to analyze the impacts of mixing in a natural discrete fracture network. We concluded that reaction-induced changes can substantially alter the flow connectivity, especially at fracture intersections. The second set of simulations considered the problem of natural weathering of fractured mafic and ultramafic rocks in the partially saturated Earth's critical zone as a function of infiltration rates, fracture permeability, and partially saturated flow parameters. As a model system, we considered an incongruent reaction network with dissolution of forsterite and precipitation of magnesite. The behavior is complex in terms of the rate-controlling processes because of the multicomponent nature of the system as shown by the grid Peclet number: the CO 2 behavior is gas diffusion-controlled in the partially saturated zone, while the rate of water flow via the Damkӧhler number controls Mg 2+ transport through the fracture network. The amounts of carbon that can be trapped are modest, but the naturally fractured domain considered here provides a useful “base case” against which various engineered solutions can be compared.

58 GEOSCIENCES↗

Niobium coaxial cavities with internal quality factors exceeding 1.4 × 10 9 for circuit quantum electrodynamics

Group-V materials such as niobium and tantalum have become popular choices for extending the performance of circuit quantum electrodynamics (cQED) platforms, allowing for quantum processors and memories with reduced error rates and more modes. The complex surface chemistry of niobium, however, makes identifying the main modes of decoherence difficult at millikelvin temperatures and single-photon powers. We use niobium coaxial quarter-wave cavities to study the impact of etch chemistry, prolonged atmospheric exposure, and the significance of cavity conditions prior to and during cooldown—in particular, niobium hydride evolution—on single-photon coherence. We demonstrate cavities with quality factors Q int ≳ 1.4 × 10 9 in the single-photon regime, a 15-fold improvement over aluminum cavities of the same geometry. We rigorously quantify the sensitivity of our fabrication process to various loss mechanisms and demonstrate a two- to fourfold reduction in the two-level system loss tangent and a three- to fivefold improvement in the residual resistance over traditional buffered chemical polishing etching techniques. Finally, we demonstrate transmon integration and coherent cavity control while maintaining a cavity coherence of 11.3 ms. The accessibility of our method, which can be easily replicated in academic laboratory settings, together with the demonstration of its performance, mark an advancement in three-dimensional cQED.

Oriani, Andrew E. [Chicago U., EFI; Chicago U.] (O↗

Study of Self-Interaction Errors in Density Functional Calculations of Magnetic Exchange Coupling Constants Using Three Self-Interaction Correction Methods

We examine the role of self-interaction errors (SIE) removal on the evaluation of magnetic exchange coupling constants. In particular we analyze the effect of scaling down the self-interaction-correction (SIC) for three non-empirical density functional approximations (DFAs) namely, the local spin density approximation, the Perdew-Burke-Ernzerhof generalized gradient approximation, and recent SCAN family of meta-GGA functionals. To this end, we employ three one-electron SIC methods: Perdew-Zunger [Perdew, J. P.; Zunger, A. Phys. Rev. B, 1981, 23, 5048] SIC, the orbitalwise scaled SIC method [Vydrov, O. A. et al., J. Chem. Phys. 2006, 124, 094108], and the recent local scaling method [Zope, R. R. et al., J. Chem. Phys. 2019, 151, 214108]. We compute the magnetic exchange coupling constants using the spin projection and non projection approaches for sets of molecules composed of dinuclear and polynuclear H-He models, organic radical molecules, and chlorocuprate, and compare these results against accurate theories and experiment. Our results show that for the systems that mainly consist of single electron regions, PZSIC performs well but for more complex organic systems and the chlorcuprates, an overcorrecting tendency of PZSIC combined with the DFAs utilized in this work is more pronounced, and in such cases LSIC with kinetic energy density ratio performs better than PZSIC. Furthermore, analysis of the results in terms of SIC corrections to the density and to the total energy shows that both density and energy correction are required to obtain an improved prediction of magnetic exchange couplings.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Many-body reactive force field development for carbon condensation in C/O systems under extreme conditions

In this paper, we describe the development of a reactive force field for C/O systems under extreme temperatures and pressures, based on the many-body Chebyshev Interaction Model for Efficient Simulation (ChIMES). The resulting model, which targets carbon condensation under thermodynamic conditions of 6500 K and 2.5 g cm –3 , affords a balance between model accuracy, complexity, and training set generation expense. We show that the model recovers much of the accuracy of density functional theory for the prediction of structure, dynamics, and chemistry when applied to dissociative condensed phase systems at 1:1 and 1:2 C:O ratios, as well as molten carbon. Our C/O modeling approach exhibits a 104 increase in efficiency for the same system size (i.e., 128 atoms) and a linear system size scalability over standard quantum molecular dynamics methods, allowing the simulation of significantly larger systems than previously possible. We find that the model captures the condensed-phase reaction-coupled formation of carbon clusters implied by recent experiments, and that this process is susceptible to strong finite size effects. Overall, we find the present ChIMES model to be well suited for studying chemical processes and cluster formation at pressures and temperatures typical of shock waves. We expect that the present C/O modeling paradigm can serve as a template for the development of a broader high pressure–high temperature force-field for condensed phase chemistry in organic materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗