Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computer experiments”

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 1,585 records · Page 88
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Quantitatively accurate activity measurements with a dedicated cardiac SPECT camera: Physical phantom experiments

Purpose: Recently, there has been increased interest in dedicated cardiac single photon emission computed tomography (SPECT) scanners with pinhole collimation and improved detector technology due to their improved count sensitivity and resolution over traditional parallel-hole cameras. With traditional cameras, energy-based approaches are often used in the clinic for scatter compensation because they are fast and easily implemented. Some of the cardiac cameras use cadmium-zinc-telluride (CZT) detectors which can complicate the use of energy-based scatter correction (SC) due to the low-energy tail—an increased number of unscattered photons detected with reduced energy. Modified energy-based scatter correction methods can be implemented, but their level of accuracy is unclear. In this study, the authors validated by physical phantom experiments the quantitative accuracy and reproducibility of easily implemented correction techniques applied to {sup 99m}Tc myocardial imaging with a CZT-detector-based gamma camera with multiple heads, each with a single-pinhole collimator. Methods: Activity in the cardiac compartment of an Anthropomorphic Torso phantom (Data Spectrum Corporation) was measured through 15 {sup 99m}Tc-SPECT acquisitions. The ratio of activity concentrations in organ compartments resembled a clinical {sup 99m}Tc-sestamibi scan and was kept consistent across all experiments (1.2:1 heart to liver and 1.5:1 heart to lung). Two background activity levels were considered: no activity (cold) and an activity concentration 1/10th of the heart (hot). A plastic “lesion” was placed inside of the septal wall of the myocardial insert to simulate the presence of a region without tracer uptake and contrast in this lesion was calculated for all images. The true net activity in each compartment was measured with a dose calibrator (CRC-25R, Capintec, Inc.). A 10 min SPECT image was acquired using a dedicated cardiac camera with CZT detectors (Discovery NM530c, GE Healthcare), followed by a CT scan for attenuation correction (AC). For each experiment, separate images were created including reconstruction with no corrections (NC), with AC, with attenuation and dual-energy window (DEW) scatter correction (ACSC), with attenuation and partial volume correction (PVC) applied (ACPVC), and with attenuation, scatter, and PVC applied (ACSCPVC). The DEW SC method used was modified to account for the presence of the low-energy tail. Results: T-tests showed that the mean error in absolute activity measurement was reduced significantly for AC and ACSC compared to NC for both (hot and cold) datasets (p < 0.001) and that ACSC, ACPVC, and ACSCPVC show significant reductions in mean differences compared to AC (p ≤ 0.001) without increasing the uncertainty (p > 0.4). The effect of SC and PVC was significant in reducing errors over AC in both datasets (p < 0.001 and p < 0.01, respectively), resulting in a mean error of 5% ± 4%. Conclusions: Quantitative measurements of cardiac {sup 99m}Tc activity are achievable using attenuation and scatter corrections, with the authors’ dedicated cardiac SPECT camera. Partial volume corrections offer improvements in measurement accuracy in AC images and ACSC images with elevated background activity; however, these improvements are not significant in ACSC images with low background activity.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2022

Neutrinoless double-β decay: Combining quantum Monte Carlo and the nuclear shell model with the generalized contact formalism

Neutrinoless double beta decay searches can determine the Majorana nature of neutrinos, the absolute neutrino mass, and provide invaluable insights on the matter dominance of the universe. However, the uncertainty in the nuclear matrix elements that govern the decay limits the physics reach of these experiments. Here, we devise a novel framework based on the generalized contact formalism that combines the nuclear shell model and quantum Monte Carlo methods and compute the neutrinoless double-beta decay of nuclei used in the most advanced experiments, including 76 Ge, 130 Te, and 136 Xe. Our results cover all relevant terms, including the leading-order short-range operator recognized recently. We validate our method in light nuclei by comparing against accurate variational Monte Carlo results. On heavy systems we obtain reduced nuclear matrix elements compared with previous calculations due to additional correlations captured by quantum Monte Carlo and introduced within the generalized contact formalism, suggesting longer decay half-lives than previously considered. On the other hand, we find an enhancement of the nuclear matrix elements due to the new short-range operator.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

TREAT Neutronics Analysis and Design Support, Part II: Multi-SERTTA-CAL

Experiment vehicle design is necessary in preparation for Transient Reactor Test (TREAT) facility restart and the resumption of transient testing to support Accident Tolerant Fuel (ATF) characterization and other future fuels testing requirements. Currently the most mature vehicle design is the Multi-SERTTA (Static Environments Rodlet Transient Test Apparatuses), which can accommodate up to four concurrent rodlet-sized specimens under separate environmental conditions. Robust test vehicle design requires neutronics analyses to support design development, optimization of the power coupling factor (PCF) to efficiently maximize energy generation in the test fuel rodlets, and experiment safety analyses. In integral aspect of prior TREAT transient testing was the incorporation of calibration experiments to experimentally evaluate and validate test conditions in preparation of the actual fuel testing. The calibration experiment package established the test parameter conditions to support fine-tuning of the computational models to deliver the required energy deposition to the fuel samples. The calibration vehicle was designed to be as near neutronically equivalent to the experiment vehicle as possible to minimize errors between the calibration and final tests. The Multi-SERTTA-CAL vehicle was designed to serve as the calibration vehicle supporting Multi-SERTTA experimentation. Models of the Multi-SERTTA-CAL vehicle containing typical PWR-fuel rodlets were prepared and neutronics calculations were performed using MCNP6.1 with ENDF/B-VII.1 nuclear data libraries; these results were then compared against those performed for Multi-SERTTA to determine the similarity and possible design modification necessary prior to construction of these experiment vehicles. The estimated reactivity insertion worth into the TREAT core is very similar between the two vehicle designs, with the primary physical difference being a hollow Inconel tube running down the length of the calibration vehicle. Calculations of PCF indicate that on average there is a reduction of approximately 6.3 and 12.6%, respectively, for PWR fuel rodlets irradiated under wet and dry conditions. Changes to the primary or secondary vessel structure in the calibration vehicle can be performed to offset this discrepancy and maintain neutronic equivalency. Current possible modifications to the calibration vehicle include reduction of the primary vessel wall thickness, swapping Zircaloy-4 for stainless steel 316 in the secondary containment, or slight modification to the temperature and pressure of the water environment within the primary vessel. Removal of some of the instrumentation within the calibration vehicle can also serve to slightly increase the PCF. Future efforts include further modification and optimization of the Multi-SERTTA and Multi-SERTTA-CAL designs in preparation of actual TREAT transient testing. Experimental results from both test vehicles will be compared against calculational results and methods to provide validation and support additional neutronics analyses.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

TREAT neutronics analysis and design support, Part II: Multi-SERTTA-CAL

Experiment vehicle design is necessary in preparation for Transient Reactor Test (TREAT) facility restart and the resumption of transient testing to support Accident Tolerant Fuel (ATF) characterization and other future fuels testing requirements. Currently the most mature vehicle design is the Multi-SERTTA (Static Environments Rodlet Transient Test Apparatuses), which can accommodate up to four concurrent rodlet-sized specimens under separate environmental conditions. Robust test vehicle design requires neutronics analyses to support design development, optimization of the power coupling factor (PCF) to efficiently maximize energy generation in the test fuel rodlets, and experiment safety analyses. An integral aspect of prior TREAT transient testing was the incorporation of calibration experiments to experimentally evaluate and validate test conditions in preparation of the actual fuel testing. The calibration experiment package established the test parameter conditions to support fine-tuning of the computational models to deliver the required energy deposition to the fuel samples. The calibration vehicle was designed to be as near neutronically equivalent to the experiment vehicle as possible to minimize errors between the calibration and final tests. The Multi-SERTTA-CAL vehicle was designed to serve as the calibration vehicle supporting Multi-SERTTA experimentation. Models of the Multi-SERTTA-CAL vehicle containing typical PWR-fuel rodlets were prepared and neutronics calculations were performed using MCNP6.1 with ENDF/B-VII.1 nuclear data libraries; these results were then compared against those performed for Multi-SERTTA to determine the similarity and possible design modification necessary prior to construction of these experiment vehicles. The estimated reactivity insertion worth into the TREAT core is very similar between the two vehicle designs, with the primary physical difference being a hollow Inconel tube running down the length of the calibration vehicle. Calculations of PCF indicate that on average there is a reduction of approximately 6 % for PWR fuel rodlets irradiated under both wet and dry conditions. Changes to the primary or secondary vessel structure in the calibration vehicle can be performed to offset this discrepancy and maintain neutronic equivalency. Current possible modifications to the calibration vehicle include reduction of the primary vessel wall thickness, swapping Zircaloy-4 for stainless steel 316 in the secondary containment, or slight modification to the temperature and pressure of the water environment within the primary vessel. Removal of some of the instrumentation within the calibration vehicle can also serve to slightly increase the PCF. Future efforts include further modification and optimization of the Multi-SERTTA and Multi-SERTTA-CAL designs in preparation of actual TREAT transient testing. Experimental results from both test vehicles will be compared against calculational results and methods to provide validation and support additional neutronics analyses. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

Study of hypervelocity projectile impact on thick metal plates

Hypervelocity impacts generate extreme pressure and shock waves in impacted targets that undergo severe localized deformation within a few microseconds. These impact experiments pose unique challenges in terms of obtaining accurate measurements. Similarly, simulating these experiments is not straightforward. This paper proposed an approach to experimentally measure the velocity of the back surface of an A36 steel plate impacted by a projectile. All experiments used a combination of a two-stage light-gas gun and the photonic Doppler velocimetry (PDV) technique. The experimental data were used to benchmark and verify computational studies. Two different finite-element methods were used to simulate the experiments: Lagrangian-based smooth particle hydrodynamics (SPH) and Eulerian-based hydrocode. Both codes used the Johnson-Cook material model and the Mie-Grüneisen equation of state. Experiments and simulations were compared based on the physical damage area and the back surface velocity. Finally, the results of this study showed that the proposed simulation approaches could be used to reduce the need for expensive experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Visualizing the Nanoscale Oxygen and Cation Transport Mechanisms during the Early Stages of Oxidation of Fe–Cr–Ni Alloy Using In Situ Atom Probe Tomography

Abstract Understanding the early stages of interactions between oxygen and material surfaces—especially at very high spatial resolutions—is highly beneficial for fields ranging from materials degradation, corrosion, geological sciences, forensics, and catalysis. The ability of in situ atom probe tomography (APT) is demonstrated to track the diffusion of oxygen and metal ions at nanoscale spatial resolution during the early stages of oxidation of a model Fe–Cr–Ni alloy. Using 18 O isotope tracers in these in situ APT experiments and complementary ex situ multimodal microscopy, spectroscopy, and computational simulations allows to precisely analyze the kinetics of oxidation and determine that outward cation diffusion to oxide/air interface is the primary mechanism for intragranular oxide growth in this alloy at 300 °C. This unique in situ isotopic tracer APT approach and the insights gained can be highly beneficial for studying early stages of gas–surface reactions in a broad array of materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Performance Improvements on SNS and HFIR Instrument Data Reduction Workflows Using Mantid

Performance of data reduction workflows at the High Flux Isotope Reactor (HFIR) and the Spallation Neutron Source (SNS) at Oak Ridge National Laboratory (ORNL) is mainly determined by the time spent loading raw measurement events stored in large and sparse datasets. This paper describes: (1) our long-term view to leverage SNS and HFIR data management needs with our experience at ORNL’s world-class high performance computing (HPC) facilities, and (2) our short-term efforts to speed up current workflows using Mantid, a data analysis and reduction community framework used across several neutron scattering facilities. We show that minimally invasive short-term improvements in metadata management have a moderate impact in speeding up current production workflows. We propose a more disruptive domain-specific solution: the No Cost Input Output (NCIO) framework, we provide an overview, the risks and challenges in NCIO’s adoption by HFIR and SNS stakeholders.

Godoy, William↗
DOE OSTI2025

Robustness of Deep Learning Classification to Adversarial Input on GPUs: Asynchronous Parallel Accumulation Is a Source of Vulnerability

The ability of machine learning (ML) classification models to resist small, targeted input perturbations—known as adversarial attacks—is a key measure of their safety and reliability. We show that floating-point non associativity (FPNA) coupled with asynchronous parallel programming on GPUs is sufficient to result in misclassification, without any perturbation to the input. Additionally, we show that this misclassification is particularly significant for inputs close to the decision boundary and that standard adversarial robustness results may be overestimated up to 4.6 when not considering machine-level details. We first study a linear classifier, before focusing on standard Graph Neural Network (GNN) architectures and datasets used in robustness assessments. We develop a novel black-box attack using Bayesian optimization to discover external workloads that can change the instruction scheduling which bias the output of reductions on GPUs and reliably lead to misclassification. Motivated by these results, we present a new learnable permutation (LP) gradient-based approach to learning floating-point operation orderings that lead to misclassifications. The LP approach provides a worst-case estimate in a computationally efficient manner, avoiding the need to run identical experiments tens of thousands of times over a potentially large set of possible GPU states or architectures. Finally, using instrumentation-based testing, we investigate parallel reduction ordering across different GPU architectures under external background workloads, when utilizing multi-GPU virtualization, and when applying power capping. Our results demonstrate that parallel reduction ordering varies significantly across architectures under the first two conditions, substantially increasing the search space required to fully test the effects of this parallel scheduler-based vulnerability. These results and the methods developed here can help to include machine-level considerations into adversarial robustness assessments, which can make a difference in safety and mission critical applications.

Shanmugavelu, Sanjif [Maxeler Technologies, a Groq↗
DOE OSTI2025

Neutrino properties from muonium-antimuonium mixing

The nature of neutrino mass — whether neutrinos are Dirac or Majorana particles — remains one of the central open questions in particle physics. While observation of neutrinoless double beta decay would confirm the Majorana case, the absence of a signal offers no definitive insight. In light of this, we investigate muonium-antimuonium mixing — proposed for further study at the MACE experiment — as an alternative probe of neutrino properties. We compute the mixing amplitude in the Standard Model minimally extended to include massive Dirac or Majorana neutrinos, and correct previous calculations by properly treating the relevant infrared scales. As the GIM mechanism strongly suppresses the Dirac contribution, we explore whether relaxing unitarity of the PMNS matrix can enhance the mixing without obscuring neutrino properties. Surprisingly, the answer to this question is negative. We also examine the pseudo-Dirac case — predominantly Dirac neutrinos with small Majorana masses — and find that this scenario can significantly enhance the mixing compared to the pure Dirac case, especially for normal mass ordering.

Neutrino interactions↗
DOE OSTI2018

Study on the Impact Characteristics of Coherent Supersonic Jet and Conventional Supersonic Jet in EAF Steelmaking Process

Supersonic oxygen-supplying technologies, including the coherent supersonic jet and the conventional supersonic jet, are now widely applied in electric arc furnace steelmaking processes to increase the bath stirring, reaction rates, and energy efficiency. However, there has been limited research on the impact characteristics of the two supersonic jets. In the present study, by integrating theoretical modeling and numerical simulations, a hybrid model was developed and modified to calculate the penetration depth and impact zone volume of the coherent and conventional supersonic jets. The computational fluid dynamics results were validated against water model experiments. The results show that the lance height has significant influence on the jet penetration depth and jet impact zone volume. The penetration depth decreases with increasing lance height, whereas the jet impact zone volume initially increases and then decreases with increasing lance height. In addition, the penetration depth and impact zone volume of the coherent supersonic jet are larger than those of the conventional supersonic jet at the same lance height, which illustrates the advantages of the coherent supersonic jet in delivering great amounts of oxygen to liquid melt with a better stirring effect compared to the conventional supersonic jet. A newly defined parameter, the k value, reflects the velocity attenuation and the potential core length of the main supersonic jet. Finally, a hybrid model and its modifications can well predict the penetration depth and impact zone volume of the coherent and conventional supersonic jets.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Synergistic Effect of Nitrogen and Refractory Material on TiN Formation and Equiaxed Grain Structure of Ferritic Stainless Steel

The effect of nitrogen content on the formation of an equiaxed solidification structure of Fe-16Cr steel was investigated. Moreover, two different kinds of refractory materials, i.e., alumina and magnesia, were employed to control the type of oxide inclusion. The characteristics of TiN(-oxide) inclusions were quantitatively analyzed in both molten steel and solidified samples. When the melting was carried out in the alumina refractory, the grain size continuously decreased with increasing nitrogen content. However, a minimum grain size was observed at a specific nitrogen content (approx. 150 ppm) when the steel was melted in the magnesia refractory. Most of the single TiN particles had a cuboidal shape and fine irregularly shaped particles were located along the grain boundary due to the microsegregation of Ti at the grain boundary during solidification. The type of TiN-oxide hybrid inclusion was strongly affected by the refractory material where Al{sub 2}O{sub 3}-TiN and MgAl{sub 2}O{sub 4}-TiN hybrid-type inclusions were obtained in the alumina and magnesia refractory experiments, respectively. The formation of oxide inclusions was well predicted by thermochemical computations and it was commonly found that oxide particles were initially formed, followed by the nucleation and growth of TiN. When the nitrogen content increased, the number density of TiN linearly increased in the alumina refractory experiments. However, the number of TiN exhibits a maximum at about [N] = 150 ppm, at which a minimum grain size was obtained in the magnesia refractory experiments. Therefore, the larger the number density of TiN, the smaller the primary grain size after solidification. The number density of TiN in the steel melted in the magnesia refractory was greater than that in the steel melted in the alumina refractory at given Ti and N contents, which was due to the lower planar lattice disregistry of MgAl{sub 2}O{sub 4}-TiN interface rather than that of Al{sub 2}O{sub 3}-TiN interface. When ∆T{sub TiN} (= difference between the TiN precipitation temperature and the liquidus of the steel) was 20 K to 40 K, the number density of effective TiN was maximized and thus, the grain size was minimized after solidification. Finally, although most of the TiN particles were smaller than 1 μm in the molten steel samples irrespective of the nitrogen content, TiN particles larger than 10 μm were observed in the solidified samples when the nitrogen content was greater than 150 ppm. The growth of TiN particles during melting and solidification was well predicted by the combinatorial simulation of the ‘Ostwald ripening model’ based on the Lifshitz–Slyozov–Wagner theory in conjunction with the ‘Diffusion controlled model’ using Ohnaka’s microsegregation equation.

36 MATERIALS SCIENCE↗
DOE OSTI2018

A cooperative conjugate gradient method for linear systems permitting efficient multi-thread implementation

This paper revisits, in a multi-thread context, the so-called multi-parameter or block conjugate gradient (B-CG) methods, first proposed as sequential algorithms by O’Leary and Brezinski, for the solution of the linear system Ax=b, for an n-dimensional symmetric positive definite matrix A. Instead of the scalar parameters of the classical CG algorithm, which minimizes a scalar functional at each iteration, multiple descent and conjugate directions are updated simultaneously. Implementation involves the use of multiple threads and the algorithm is referred to as cooperative CG (CCG) to emphasize that each thread now uses information that comes from the other threads. It is shown that for a sufficiently large matrix dimension n, the use of an optimal number of threads results in a worst case flop count of O(n{sup 7/3}) in exact arithmetic. Numerical experiments on a multi-core, multi-thread computer, for synthetic and real matrices, illustrate the theoretical results.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Modeling low cycle fatigue (LCF) of additively manufactured Hastelloy X using An accelerated crystal plasticity fatigue damage model

This paper presents a microstructure-based model for low cycle fatigue (LCF) behavior and life of Nickel-based alloy Hastelloy X manufactured using laser-powder bed fusion (L-PBF) additive manufacturing (AM). AM Hastelloy X, a solution-strengthened alloy, is tested at elevated temperature under fully reversed LCF conditions at different strain levels. A generalized plane strain finite element model is generated from electron backscatter diffraction (EBSD) characterization. The constitutive behavior of the material under fatigue is modeled using crystal plasticity and calibrated with both monotonic tensile and cyclic stress–strain data. The fatigue micro-crack initiation and propagation in the microstructure is modeled using a modified Chaboche fatigue damage model. An embedded boundary condition with a homogenous medium is used to apply the cyclic deformation and prevent numerically introduced over-constraints during fatigue simulation. A ‘cycle-jump’ method is used to accelerate the fatigue simulation and reduce the computational cost. The simulation results are compared to LCF experiments, showing satisfactory matches in cyclic stress behavior and number of cycles to macro-crack initiation for all applied strain ranges. In addition, the model illustrates the potential for quantifying microscale fatigue life impacting factors such as microstructure and surface roughness, which is needed to accurately quantify the reliability of AM components in service.

36 MATERIALS SCIENCE↗
DOE OSTI2021

OC6 Phase Ib: Validation of the CFD predictions of difference-frequency wave excitation on a FOWT semisubmersible

During the previous OC5 11 project, state-of-the-art mid-fidelity engineering tools for floating wind systems were found to consistently underpredict the nonlinear, low-frequency responses of semisubmersible offshore wind platforms, leading to substantial errors in the structural loads because of low-frequency surge and pitch resonance. To examine this underprediction, a coordinated investigation with computational fluid dynamics (CFD) simulations and model-basin experiments was carried out. Both investigations involved a fixed and simplified OC5-DeepCwind semisubmersible in bichromatic waves. The wave excitations - especially the nonlinear, difference-frequency excitation - on the structure from the CFD simulations were compared to the experimental measurements for validation, with uncertainty analyses for both the experimental and the CFD results. Furthermore, the wave excitations on each column of the semisubmersible were measured separately in the experiment, allowing the validation of the CFD results to be done on a per-column basis. Overall, the CFD predictions of the difference-frequency excitations agree with the experimental measurements, suggesting the CFD solutions can be used as a reference for tuning and improving the engineering-level tools and can provide a means to better understand the underprediction at low frequencies.

16 TIDAL AND WAVE POWER↗
DOE OSTI2024

Ion-temperature- and density-gradient-driven instabilities and turbulence in Wendelstein 7-X close to the stability threshold

Electrostatic gyrokinetic instabilities and turbulence in the Wendelstein 7-X stellarator are studied. Particular attention is paid to the ion-temperature-gradient (ITG) instability and its character close to marginal stability [Floquet-type turbulence (Zocco et al., Phys. Rev. E, vol. 106, 2022, p. L013202) with no electron temperature gradient]. The flux tube version of the δƒ code stella (Barnes et al., J. Comput. Phys., vol. 391, 2019, pp. 365–380) is used to run linear and nonlinear gyrokinetic simulations with kinetic electrons. The nature of the dominant instability depends on the wavelength perpendicular to the magnetic field, and the results are conveniently displayed in stability diagrams that take this dependence into account. This approach highlights the presence of universal instabilities, which are less unstable but have longer wavelengths than other modes. A quasi-linear estimate of the heat flux suggests they are relevant for transport. Close to the stability threshold, the linear eigenmodes and turbulence form highly extended structures along the computational domain if the magnetic shear is small. Numerical experiments and diagnostics are undertaken to assess the resulting radial localisation of the turbulence, which affects the interaction of the latter with zonal flows. Increasing the amplitude of the magnetic shear (e.g.through current drive) has a stabilising effect on the turbulence and, thus, reduces the nonlinear energy transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Is solid copper oxalate a spin chain or a mixture of entangled spin pairs?

Macroscopic assemblies of interacting spins give rise to a broad spectrum of behaviors determined by the spatial arrangement of the magnetic sites and the electronic interactions between them. Compounds of copper (II), in which each copper carries spin 1/2, exhibit a vast variety of physical properties. For antiferromagnetically coupled spin sites, there are two limiting scenarios: spin chains in which the spins can exhibit a long-range order or a mixture of dimers in which the spins within each pair are entangled but do not communicate with the spins from other dimers. In principle, the two types can be distinguished on the basis of experimental observations and modeling using empirically parameterized effective Hamiltonians, but in practice, ambiguity may persist for decades, as is the case for copper oxalate. Here we use high-level ab initio calculations to establish the validity of the nearest-site Heisenberg model and to predict the interaction strength between the magnetic sites. Finally, the computed magnetic susceptibility provides an unambiguous interpretation of magnetic experiments performed throughout half a century, clearly supporting the infinite spin-chain behavior of solid copper oxalate.

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

Describing the Diverse Geometries of Gold from Nanoclusters to Bulk—A First-Principles-Based Hybrid Bond-Order Potential

Molecular dynamics simulations using empirical force fields (EFFs) are crucial for gaining fundamental insights into atomic structure and long time scale dynamics of Au nanoclusters with far-reaching applications in energy and devices. This approach is thwarted by the failure of currently available EFFs in describing the size-dependent dimensionality and diverse geometries exhibited by Au clusters (e.g., planar structures, hollow cages, tubes, pyramids, space-filled structures). Here, we mitigate this issue by introducing a new hybrid bond-order potential (HyBOP), which accounts for (a) short-range interactions via Tersoff-type BOP terms that accurately treat bond directionality and (b) long-range dispersion effects by a scaled Lennard–Jones term whose contribution depends on the local atomic density. We optimized the independent parameters for our HyBOP using a global optimization scheme driven by genetic algorithms. Moreover, to ensure good transferability of these parameters across different length scales, we used an extensive training data set that encompasses structural and energetic properties of 1000 13-atom Au clusters, surface energies, as well as bulk polymorphs, obtained from density functional theory (DFT) calculations. Our newly developed HyBOP has been found to accurately describe (a) global minimum energy configurations at different cluster sizes as well as order of stability of various cluster configurations at any size, (b) critical size of transition from planar to globular clusters, (c) evolution of structural motifs with cluster size, and (c) thermodynamics, structure, elastic properties, and energetic ordering of bulk condensed phases as well as surfaces, in excellent agreement with DFT calculations and spectroscopic experiments. This makes our newly developed HyBOP a valuable, computationally robust but inexpensive tool for investigating a wide range of materials phenomena occurring in Au at the atomistic level.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Folding and Assembly of Short α, β, γ -Hybrid Peptides: Minor Variations in Sequence and Drastic Differences in Higher-Level Structures

Multilevel protein structures typically involve polypeptides of sufficient lengths. Here we report the folding and assembly of seven short tetrapeptides sharing the same types of α-, β-, and aromatic γ-amino acid residues. These are two sets of hybrid peptides, with three members in one set and four in the other, having complementary hydrogen-bonding sequences that were hypothesized to pair into linear H-bonded duplexes. However, instead of undergoing the anticipated pairing, the initially examined three oligomers, 1 and 2a or 2b, differing only in their central αβ hybrid dipeptide sequence, do not associate with each other and exhibit distinctly different folding behavior. Experiments based on NMR and mass spectrometry, along with computational studies and systematic inference, reveal that oligomer 1 folds into an expanded β-turn containing an unusual hybrid α/β-amino acid sequence composed of glycine and β-alanine, two α- and β-amino acid residues that are conformationally most flexible, and peptides 2a and 2b adopt a noncanonical, extended helical conformation and dimerize into double helices undergoing rapid conformational exchange or helix inversion. The different central dipeptide sequences, αβ vs βα, result in drastically different intramolecular H-bonding patterns that are responsible for the observed folding behavior of 1 and 2. The revealed turn and double helix have few natural or synthetic counterparts, and provide novel and unique folding prototypes based on which chiral α- and β-amino acids are incorporated. The resultant derivatives 1a, 1b, 2c, and 2d follow the same folding and assembling behavior and demonstrate the generality of this system with the formation of expanded β-turns and double helices with enhanced folding stabilities, hampered helix inversion, as well as defined and dominant helical sense. This work has demonstrated the unique capability of synthetic foldamers in generating structures with fascinating folding and assembling behavior. As a result, the revealed systems offer ample opportunity for further structural optimization and applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗