Search NASA⌕ Search

SEARCH · Search NASA

Results for “CONSISTENCY”

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 163 records · Page 9

Transfer learning-based soybean LAI estimations by integrating PROSAIL, UAV, and PlanetScope imagery

Accurate Leaf Area Index (LAI) estimations at the soybean plot scale is achievable using high-resolution Unmanned Aerial Vehicle (UAV) imagery and field measurement samples. However, the limited coverage of UAV flights restricts large-scale remote sensing monitoring in expansive soybean fields. This study leverages the broad coverage and 3-m resolution of PlanetScope satellite imagery to extend LAI prediction from UAV to satellite scales through transfer learning, using UAV-scale LAI estimates as a benchmark to validate cross-scale consistency. To address this challenge, this study proposed the LAI-TransNet, a two-stage transfer learning framework designed for precise and scalable soybean LAI prediction across large areas, demonstrating its effectiveness in cross-scale monitoring. In Stage 1, a UAV-scale benchmark is established using PROSAIL-simulated UAV reflectance data (UAV-Sim) and field-measured soybean LAI. Traditional machine learning, deep learning, and transfer learning models are trained on a hybrid UAV-Sim and field-measured dataset (UAV-Sim_Measured), with the transfer learning model CNN-TL, fine-tuned using pre-trained weights derived from UAV-Sim, achieving the highest accuracy (R 2 = 0.81, RMSE = 0.64 m 2 /m 2 , rRMSE = 11.5 %). In Stage 2, LAI-TransNet is developed by fine-tuning the CNN-TL model on PlanetScope simulated data (PS-Sim), preprocessed via cross-domain mapping to align UAV and satellite spectral features. Real PlanetScope imagery is corrected for reflectance consistency with reference to UAV imagery spectral profiles. LAI-TransNet outperforms other deep learning models trained directly on PS-Sim (R 2 = 0.69 vs. 0.60–0.63), ensuring robust cross-scale consistency. In conclusion, by bridging UAV and satellite scales, LAI-TransNet enables large-scale soybean LAI monitoring, enhancing precision agriculture management through improved monitoring with the PlanetScope imagery.

Leaf area index (LAI)↗

Sub-melt nanosecond pulsed-laser induced densification and strain-field relaxation in single-crystal diamond

Dislocations and polishing-induced defect networks in synthetic diamond introduce local strain fields and broaden Raman features, limiting performance in optical, thermal, and electronic applications. Laser annealing is emerging as a promising approach to repair surface and near-surface defects in diamond without entering the melt regime, yet surface densification, defect-state modification, and associated structural changes have not been well quantified. In this work, we show that sub-melt nanosecond pulsed-laser annealing (PLA) induces near-surface densification and defect-mediated strain relaxation in single-crystal Chemical Vapor Deposition (CVD) diamond. Single- and two-pulse PLA were applied, and structural evolution was quantified using co-registered ISO 25,178 white-light interferometry, depth-resolved Raman spectroscopy, and cross-sectional STEM with geometric phase analysis (GPA). Across a 5 × 6 grid (n = 30), responsive regions exhibit large reductions in local slope (Sdq 45–65%), developed area (Sdr 60–90%), height spread (Sp, Sz 30–65%), void volume (Vv 57–60%), and roughness amplitude (Sa, Sq 48–57%), consistent with densification of ∼4–6.5 nm. Raman profiling shows narrowing of the diamond line and improved spectral uniformity to depths of ∼2–3 μm. Given that the Raman probing depth significantly exceeds the densified layer thickness, this response is interpreted as consistent with long-range strain-field redistribution originating from the near-surface region. STEM-GPA strain maps further support this interpretation, showing smoother strain fields, suppressed hotspots, and redistribution of localized strain concentrations following PLA. These results are consistent with defect-mediated strain relaxation and densification-driven modification of the near-surface energy state. The approach provides a scalable pathway for improving near-surface structural quality in diamond relevant to electronic, photonic, and quantum applications.

Areal surface metrology (ISO 25,178)↗

Kinetic and structural evidence for specific DMSO interference with reversible binding of uncharged bis-oximes to hAChE and their reactivation kinetics of OP-hAChE

The structural basis of inhibitory effect of organic solvent dimethyl sulfoxide (DMSO) on human acetylcholinesterase (EC 3.1.1.7; hAChE) was inferred from the effect of DMSO on kinetics of reversible inhibition of uncharged, heterocyclic bis-oximes to hAChE, from DMSO effect on rates of reactivation of inactive organophosphate (OP)-hAChE conjugates by bis-oximes and by X-ray structures of bis-oxime and DMSO binding to hAChE. Here, the reversible inhibition constant of DMSO for hAChE in 0.1 M phosphate buffer pH 7.4 at 22 °C, was K i = (0.32 ± 0.04) % (or 45 ± 5 mM). The K i of the bis-oxime LG-703 for hAChE was 3.2-fold larger in 1 % DMSO, consistent with direct competition between LG-703 and DMSO. The X-ray structure of the LG-703∗hAChE complex (PDB ID: 6U3P) shows DMSO and LG-703 bound to individual hAChE monomers, LG-703 in the chain A and DMSO in the chain B. In the co-crystallization both small molecules were present at a similar excess over their corresponding K i values for hAChE (7.8-fold for DMSO and 6.5-fold for LG-703) and formation of two different complexes (DMSO∗hAChE and LG-703∗hAChE), in the same crystal, appears consistent with inhibition kinetics. Furthermore, rates of reactivation of paraoxon-inhibited hAChE (POX-hAChE) and of VX-hAChE by LG-703 and by a novel heterocyclic bis-oxime LG-1922 were reduced 2 – 3-fold in DMSO, consistent with observation of the active-center-bound DMSO molecules in the newly solved structure of the LG-1922∗POX-hAChE complex presented here and in our POX-hAChE structure (PDB ID: 8DT2) showing obstruction of the reactivator access to the conjugated P atom.

Acetylcholinesterase inhibition↗

Amantadine interactions with phase separated lipid membranes

Amantadine, a small amphilphic organic compound that consists of an adamantane backbone and an amino group, was first recognized as an antiviral in 1963 and received approval for prophylaxis against the type A influenza virus in 1976. Since then, it has also been used to treat Parkinson’s disease-related dyskinesia and is being considered as a treatment for corona viruses. Since amantadine usually targets membrane-bound proteins, its interactions with the membrane are also thought to be important. Biological membranes are now widely understood to be laterally heterogeneous and certain proteins are known to preferentially co-localize within specific lipid domains. Does amantadine, therefore, preferentially localize in certain lipid composition domains? To address this question, here, we studied amantadine’s interactions with phase separating membranes composed of cholesterol, DSPC (1,2-distearoyl-sn-glycero-3-phosphocholine), POPC (1-palmitoyl-2-oleoyl-glycero-3-phosphocholine), and DOPC (1,2-dioleoyl-sn-glycero-3-phosphocholine), as well as single-phase DPhPC (1,2-diphytanoyl-sn-glycero-3-phos-phocholine) membranes. From Langmuir trough and differential scanning calorimetry (DSC) measurements, we determined, respectively, that amantadine preferentially binds to disordered lipids, such as POPC, and lowers the phase transition temperature of POPC/DSPC/cholesterol mixtures, implying that amantadine increases membrane disorder. Further, using droplet interface bilayers (DIBs), we observed that amantadine disrupts DPhPC membranes, consistent with its disordering properties. Finally, we carried out molecular dynamics (MD) simulations on POPC/DSPC/cholesterol membranes with varying amounts of amantadine. Consistent with experiment, MD simulations showed that amantadine prefers to associate with disordered POPC-rich domains, domain boundaries, and lipid glycerol backbones. Since different proteins co-localize with different lipid domains, our results have possible implications as to which classes of proteins may be better targets for amantadine.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Experimental and theoretical investigation into the high pressure deflagration products of 2,6-diamino-3,5-dinitropyrazine-1-oxide (LLM-105)

Diamond anvil cell (DAC) laser ignition experiments and reactive ab initio molecular dynamics (AIMD) simulations were performed on the high explosive (HE) LLM-105 to investigate its high pressure (HP) deflagration chemistry. Raman and optical spectroscopy measurements reveal LLM-105 reacts into an opaque carbonaceous product at 4–25 GPa. At pressures >~ 27 GPa, the reaction product consists of an amorphous optically transparent solid and nitrogen (N 2 ) in the solid phase. While not a one-to-one comparison due to the small time and length scales, the HP AIMD simulations show that some of the product is molecular N 2 , in qualitative agreement with experiment, while above 20 GPa most of the product consists of large amorphous C x H y N z O k clusters. Clustering is enhanced with pressure and reduces with temperature. In the experiments with initial sample pressure >~ 25 GPa, the pressure within the DAC decreases with minimal change in DAC cavity area. At initial sample pressures of 43.9 GPa, when quenched to 0 K, simulations predict a product experiencing a lower pressure consistent with the experimental measurement at lower load pressures. In conclusion, the results are important for understanding the HP deflagration chemistry of LLM-105.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Two-dimensional burn velocity and analysis of burn products of nitromethane at high pressure

We have studied the laser-initiated deflagration of nitromethane (CH 3 NO 2 ) under high static pressure in the diamond anvil cell. Time-resolved images of the deflagrations have been obtained using intensified CCDs (ICCD). In contrast with previous work, we rely on spontaneous emission from the reaction, rather than changes in the speckle pattern produced by artificial illumination of the sample. Furthermore, as opposed to the 1D records obtained previously with streak cameras, ICCDs permit imaging of the burn in two dimensions providing the ability to directly observe anisotropic deflagration behavior. Here, we report several examples of this behavior and discuss its possible origins. We have also investigated the products of the reaction using Raman spectroscopy. At pressures below 25 GPa, the burn product is observed to be opaque and has a Raman spectrum consistent with a carbonaceous soot. At pressures above 25 GPa, the burn product is observed to be transparent (consistent with earlier reports) and its Raman spectrum reveals the presence of molecular N 2 and a nitrogen-carbonate species. The latter species can be recovered to ambient pressure and has a Raman spectrum consistent with ammonium carbonate or ammonium bicarbonate. To aid interpretation of the Raman spectra, measurements have also been made on the isotopologues CH 3 15 NO 2 and 13 CH 3 NO 2 . This work establishes the importance of advanced time-resolved imaging to reveal the details of deflagration under high pressure and further advances our understanding of the corresponding chemistry.

Deflagration↗

Generative learning for slow manifolds and bifurcation diagrams

In dynamical systems characterized by separation of time scales, the approximation of so called “slow manifolds”, on which the long term dynamics lie, is a useful step for model reduction. Initializing on such slow manifolds is a useful step in modeling, since it circumvents fast transients, and is crucial in multiscale algorithms (like the equation-free approach) alternating between fine scale (fast) and coarser scale (slow) simulations. In a similar spirit, when one studies the infinite time dynamics of systems depending on parameters, the system attractors (e.g., its steady states) lie on bifurcation diagrams (curves for one-parameter continuation, and more generally, on manifolds in state parameter space. Sampling these manifolds gives us representative attractors (here, steady states of ODEs or PDEs) at different parameter values. Algorithms for the systematic construction of these manifolds (slow manifolds, bifurcation diagrams) are required parts of the “traditional” numerical nonlinear dynamics toolkit. In more recent years, as the field of Machine Learning develops, conditional score-based generative models (cSGMs) have been demonstrated to exhibit remarkable capabilities in generating plausible data from target distributions that are conditioned on some given label. It is tempting to exploit such generative models to produce samples of data distributions (points on a slow manifold, steady states on a bifurcation surface) conditioned on (consistent with) some quantity of interest (QoI, observable). In this work, we present a framework for using cSGMs to quickly (a) initialize on a low-dimensional (reduced-order) slow manifold of a multi-time-scale system consistent with desired value(s) of a QoI (a “label”) on the manifold, and (b) approximate steady states in a bifurcation diagram consistent with a (new, out-of-sample) parameter value. This conditional sampling can help uncover the geometry of the reduced slow-manifold and/or approximately “fill in” missing segments of steady states in a bifurcation diagram. Finally, the quantity of interest, which determines how the sampling is conditioned, is either known a priori or identified using manifold learning-based dimensionality reduction techniques applied to the training data.

Dynamical systems↗

Comparison of interlaminar damage modeling strategies for hybrid composite/aluminum laminates subjected to low-velocity impact

Low-velocity impact of hybrid metal-composite structures was investigated experimentally and computationally. Composite laminates consisting of 2D woven glass fiber reinforced polymer (GFRP) and carbon fiber reinforced polymer (CFRP) were joined with a 6061-T6 aluminum plate using an epoxy adhesive. Two variations of the structure were studied; one consisting of all plies oriented at 0° and one consisting of all plies oriented at 45°. A drop tower was used to impact structures at a range of energies, including energies above and below the threshold at which the aluminum layer was perforated. Numerical simulations were implemented using Sierra/SM, an in-house transient dynamics finite element code developed at Sandia National Laboratories. A Hosford plasticity model was used to describe the response of the aluminum layer. A newly implemented orthotropic continuum damage mechanics (CDM) constitutive model was used to represent the composite laminate. This 3D-CDM model was compared to a cohesive zone model (2D-CDM/CZM) to investigate efficacy of aluminum perforation energy prediction, delamination prediction, and computational cost. Accuracy of each model was evaluated using the experimental results. Each showed good agreement with the tests for both the force and velocity histories, as well as the observed damage mechanisms. The 2D-CDM/CZM model was marginally more accurate in capturing both the composite and aluminum behavior — this model averaged error percentages of -11.2% and 10.8% for residual velocity and peak force, respectively. Meanwhile, the 3D-CDM model predictions yielded average error percentages of -35.5% (velocity) and 22.6% (force). However, the 3D-CDM model generally resulted in a decreased computational cost; the average run time was 14% shorter than the 2D-CDM/CZM model and 3x as many timesteps per hour were computed using the same computational resources. In conclusion, new experimental data on the impact and perforation resistance of metal-composite laminates is presented in addition to numerical predictions of the impact behavior.

Carbon fiber↗

Framework to select robust energy retrofit measures for residential communities

Residential building energy retrofits are essential for enhancing environmental sustainability and reducing energy costs. The selection of retrofit measures is influenced by factors such as building systems, occupant behavior, government policy, weather variability, and climate change, all of which can significantly impact energy performance. Compared to retrofitting individual homes, evaluating and selecting optimal retrofit solutions for an entire community is challenging due to diverse residential compositions and variability present. Therefore, engineering robustness is crucial for ensuring consistent energy performance and resilience across different conditions. In this context, robustness refers to the ability of a retrofit measure to maintain its functionality and remain an optimal choice despite external disturbances or changes in inputs and conditions. This study presents a framework for evaluating the robustness of multiple retrofit measures across various building systems, occupant behaviors, and environmental scenarios at the community level. The framework comprises five key steps: scenario model development, integration of the National Residential Efficiency Measures database, energy performance simulation, cost-benefit aggregation, and retrofit solution selection. Each step enhances the framework’s robustness by incorporating the diversity of building characteristics, occupant behaviors, environmental conditions, retrofit options, and evaluation criteria. The framework’s effectiveness is demonstrated through a case study in southern Michigan in the United States, which includes 63 one-story single-family houses, 121 two-story single-family houses, and 8 townhouses. The study identifies furnace retrofits as the most robust solution for the entire community, consistently achieving source energy reductions of 4.7 %–8.0 % and payback period of 10–20 years across various scenarios. These findings are consistent with previous research, indicating the framework’s potential for broader applications in optimizing community-scale residential energy retrofits.

Shu, Lei↗

Physics vs structure: A systematic benchmark of learning strategies for multi-zone building thermal dynamics

Recent advances in physics-informed and data-driven machine learning promise improved thermal models for advanced building control, yet there is limited quantitative evidence on when added physics structure and architectural complexity are beneficial. Here, this work presents a systematic benchmark of five representative system identification methods for modeling multi-zone building thermal dynamics: linear state-space models, multi-layer perceptrons, neural state-space models, neural ordinary differential equations, and physically-consistent neural networks. The methods are evaluated across multiple data regimes and zone coupling strategies. Using a high-fidelity multi-zone commercial building emulator, we examine short-term and long-term prediction accuracy, computational efficiency, and ease of development. Our results reveal critical trade-offs between prediction performance, model complexity, and physical consistency. We demonstrate that decoupled, nonlinear black-box models consistently outperform coupled physics-constrained architectures in both predictive accuracy and out-of-distribution robustness in majority of the test cases for the building type considered in the study. Our findings quantify the cost of complexity in building thermal modeling and provide concrete, actionable, scenario-based guidelines for selecting model classes for control-oriented applications.

Building thermal modeling↗

Hourglass control in staggered-grid hydrodynamics using virtual element stabilization techniques

Numerical simulations using the staggered-grid hydrodynamics (SGH) discretization suffer from hourglass instabilities. In this work, we develop a stabilization method to suppress the hourglass instabilities using techniques from the virtual element method (VEM). The stiffness matrix of the VEM consists of two terms: the consistency matrix which is rank deficient and the stability matrix. Here, we first show that in two dimensions and on general polygons, the stiffness matrix of the SGH is identical to the consistency matrix of the linear VEM for both the diffusion equation and the linear elasticity equation. These analyses explain the origin of the hourglass instabilities of the SGH discretization method, and establish a theoretical foundation for our proposed stabilization method by augmenting the stiffness matrix of the SGH discretization using the VEM stability matrix. Then, we present numerical examples using Lagrangian SGH simulations. The numerical experiments demonstrate that the proposed VEM stabilization method is effective at eliminating hourglass modes in the SGH discretization.

97 MATHEMATICS AND COMPUTING↗

Nondestructive ultrasonic characterization of a triple-weld-bead wire arc additively manufactured ER70S-6 S-curved wall

Understanding build-scale microstructural variation in wire arc additive manufacturing (WAAM) of low-carbon steels is essential for ensuring consistent structure–property relationships throughout large components. Conventional destructive characterization techniques, such as scanning electron microscopy (SEM) and electron backscatter diffraction (EBSD), are time-intensive and limited to localized regions, making comprehensive evaluation of large WAAM structures challenging. In this study, a nondestructive ultrasonic approach was employed to characterize a 252 mm tall ER70S-6 S-curved WAAM wall produced using a triple-bead deposition strategy. Optical and SEM analysis revealed a repeating dual-region microstructure consisting of uniform polygonal ferrite at melt pool centers and heterogeneous ferrite with coarse and fine grains near melt pool boundaries, attributed to cyclic thermal conditions. Longitudinal ultrasonic backscatter imaging was used to evaluate the continuity of this periodicity along the full build height. The ultrasonic response exhibited a consistent repeating pattern that correlated with the observed layer-wise microstructural variation. X-ray computed tomography confirmed the absence of detectable porosity, indicating that ultrasonic contrast is primarily governed by grain morphology. Overall, the results indicate that longitudinal backscatter ultrasound is a promising nondestructive characterization technique for validating microstructural variations along the s-curved WAAM wall, with significant potential for microstructure optimization and process control.

36 MATERIALS SCIENCE↗

Hydrogenation of calcite and change in chemical bonding at high pressure: Diamond formation above 100GPa

Synchrotron X-ray diffraction (XRD) and Raman spectroscopy in laser heated diamond anvil cells and first principles molecular dynamics (FPMD) calculations have been used to investigate the reactivity of calcite and molecular hydrogen (H 2 ) at high pressures up to 120 GPa. We find that hydrogen reacts with calcite starting below 0.5 GPa at room temperature forming chemical bonds with carbon and oxygen. This results in the unit cell volume expansion; the hydrogenation level is much higher for powdered samples. Single-crystal XRD measurements at 8–24 GPa reveal the presence of previously reported III, IIIb, and VI calcite phases; some crystallites show up to 4% expansion, which is consistent with the incorporation of ≤ 1 hydrogen atom per formula unit. At 40–102 GPa XRD patterns of hydrogenated calcite demonstrate broadened features consistent with the calcite VI structure with incorporated hydrogen atoms. Above 80 GPa, the C–O stretching mode of calcite splits suggesting a change in the coordination of C–O bonds. Laser heating at 110 GPa results in the formation of C–C bonds manifested in the crystallization of diamond recorded by in situ XRD at 300 K and 110 GPa and by Raman spectroscopy on recovered samples commenced with C 13 calcite. We explored several theoretical models, which show that incorporation of atomic hydrogen results in local distortions of CO 3 groups, formation of corner-shared C–O polyhedra, and chemical bonding of H to C and O, which leads to the lattice expansion and vibrational features consistent with the experiments. In conclusion, the experimental and theoretical results support recent reports on tetrahedral C coordination in high-pressure carbonate glasses and suggest a possible source of the origin of ultradeep diamonds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Study of WH production through vector boson scattering and extraction of the relative sign of the W and Z couplings to the Higgs boson in proton-proton collisions at $\sqrt{s} = 13$ TeV

A search for the production of a W boson and a Higgs boson through vector boson scattering (VBS) is presented, using CMS data from proton-proton collisions at $\sqrt{s} = 13$ TeV collected from 2016 to 2018. The integrated luminosity of the data sample is 138 fb −1 . Selected events must be consistent with the presence of two jets originating from VBS, the leptonic decay of the W boson to an electron or muon, possibly also through an intermediate τ lepton, and a Higgs boson decaying into a pair of b quarks, reconstructed as either a single merged jet or two resolved jets. A measurement of the process as predicted by the standard model (SM) is performed alongside a study of beyond-the-SM (BSM) scenarios. The SM analysis sets an observed (expected) 95% co­fidence level upper limit of 14.3 (9.9) on the ratio of the measured VBS WH cross section to that expected by the SM. The BSM analysis, conducted within the so-called 𝜅 framework, excludes all scenarios with 𝜆 WZ < 0 that are consistent with current measurements, where 𝜆 WZ = 𝜅 W ∕𝜅 Z and 𝜅 W and 𝜅 Z are the HWW and HZZ coupling mod­fiers, respectively. The significance of the exclusion is beyond 5 standard deviations, and it is consistent with the SM expectation of 𝜆 WZ = 1.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Stochastic modal velocity field in rough-wall turbulence

Stochastically generated instantaneous velocity profiles are used to reproduce the outer region of rough-wall turbulent boundary layers in a range of Reynolds numbers extending from the wind tunnel to field conditions. Each profile consists in a sequence of steps, defined by the modal velocities and representing uniform momentum zones (UMZs), separated by velocity jumps representing the internal shear layers. Height-dependent UMZ is described by a minimal set of attributes: thickness, mid-height elevation, and streamwise (modal) and vertical velocities. These are informed by experimental observations and reproducing the statistical behaviour of rough-wall turbulence and attached eddy scaling, consistent with the corresponding experimental datasets. Sets of independently generated profiles are reorganized in the streamwise direction to form a spatially consistent modal velocity field, starting from any randomly selected profile. The operation allows one to stretch or compress the velocity field in space, increases the size of the domain and adjusts the size of the largest emerging structures to the Reynolds number of the simulated flow. By imposing the autocorrelation function of the modal velocity field to be anchored on the experimental measurements, we obtain a physically based spatial resolution, which is employed in the computation of the velocity spectrum, and second-order structure functions. The results reproduce the Kolmogorov inertial range extending from the UMZ and their attached-eddy vertical organization to the very-large-scale motions (VLSMs) introduced with the reordering process. The dynamic role of VLSM is confirmed in the –u'w' co-spectra and in their vertical derivative, representing a scale-dependent pressure gradient contribution.

42 ENGINEERING↗

Rhenium-Sulfido and -Dithiolato Corroles: Reflections on Chalcophilicity

The high-temperature (~180 °C) reaction between free-base meso-triarylcorroles and Re 2 (CO) 10 , followed by exposure to PCl 3 and thiols (or elemental sulfur), affords rhenium-sulfido (ReS) corroles in 67-76% yields. The use of shorter reaction times, lower temperatures (~130 °C), and a dithiol (e.g., ethane-1,2-dithiol) also allows the isolation of rhenium-dithiolato corroles, presumptive intermediates on the path to ReS corroles. The ReS corroles exhibit high thermal stability and two reversible oxidations and reductions in their cyclic voltammograms, with redox potentials nearly identical to those observed for analogous ReO corroles. The electrochemical HOMO-LUMO gaps of the complexes, at 2.2 eV, are consistent with ligand-centered oxidation and reduction. The UV-vis spectra of ReS corroles, on the other hand, differ significantly from those of their ReO counterparts. Scalar-relativistic DFT calculations suggest that this difference reflects low-energy LUMO+2 and LUMO+3 levels, consisting of Re-S π-antibonding interactions; the ReO corroles, in contrast, exhibit a larger LUMO+1/LUMO+2 gap, as expected for a relatively classical Gouterman-type metalloporphyrin analogue. The high stability of ReS corroles is consistent with geochemists' view of rhenium as a moderately chalcophilic element (i.e., one that partitions into sulfide melts) as well as with a recent quantitative analysis of thiophilicity, which indicates that rhenium's oxophilicity and thiophilicity are essentially evenly balanced.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Strategies to Obtain Reliable Energy Landscapes from Embedded Multireference Correlated Wavefunction Methods for Surface Reactions

Embedded correlated wavefunction (ECW) theory is a powerful tool for studying ground- and excited-state reaction mechanisms and associated energetics in heterogeneous catalysis. Several factors are important to obtaining reliable ECW energies, critically the construction of consistent active spaces (ASs) along reaction pathways when using a multireference correlated wavefunction (CW) method that relies on a subset of orbital spaces in the configuration interaction expansion to account for static electron correlation, e.g., complete AS self-consistent field theory, in addition to the adequate partitioning of the system into a cluster and environment, as well as the choice of a suitable basis set and number of states included in excited-state simulations. Here, in this work, we conducted a series of systematic studies to develop best-practice guidelines for ground- and excited-state ECW theory simulations, utilizing the decomposition of NH 3 on Pd(111) as an example. We determine that ECW theory results are relatively insensitive to cluster size, the aug-cc-pVDZ basis set provides an adequate compromise between computational complexity and accuracy, and that a fixed-clean-surface approximation holds well for the derivation of the embedding potential. Additionally, we demonstrate that a merging approach, which involves generating ASs from the molecular fragments at each configuration, is preferable to a creeping approach, which utilizes ASs from adjacent structures as an initial guess, for the generation of consistent potential energy curves involving open-d-shell metal surfaces, and, finally, we show that it is essential to include bands of excited states in their entirety when simulating excited-state reaction pathways.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Prediction of Redox Potentials for U, Np, Pu, and Am in Aqueous Solution

The redox properties of the actinides in aqueous solution are important for fuel production/reprocessing and understanding the environmental impact of nuclear waste. The redox potentials for U, Np, Pu, and Am in oxidation states from 0 up to VII (as appropriate) in aqueous solutions have been predicted at the density functional theory level with the B3LYP functional, Stuttgart small core pseudopotential basis sets for the actinides, and explicit (30H 2 O molecules)/implicit treatment of the aqueous solvent using the self-consistent reaction field COSMO and SMD approaches for the implicit solvation. The predictions of the structural parameters of clusters incorporating first and second solvation shells are consistent with the available experimental data., Our results are typically within 0.2 V of the available experimental data using two explicit solvation shells with an implicit solvent model. The use of the PW91 functional substantially improved the prediction of the Pu(VI/V) redox couple. The redox couples for An(VI/IV) and An(V/IV) which involve the addition of protons and removal of the actinyl oxygens led to slightly larger differences from experiment. Here, the An(IV/0) and An(III/0) couples were reliably predicted with our approach. Predictions of the unknown An(II/I) redox potentials were negative, consistent with expectations, and predictions for unknown An(VII/VI), An(III/II), and An(II/0) redox couples improve prior estimates.

Actinides↗