Non-collinear spin structure for DFT calculations of γ- and ε-Pu [Poster]
Optimized non-collinear spin structures for γ- and ε-Pu improve the reliability of DFT calculations.
SEARCH · Search NASA
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.
Optimized non-collinear spin structures for γ- and ε-Pu improve the reliability of DFT calculations.
Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the cause of coronavirus disease 2019 (COVID-19), has spread globally and is being surveilled with an international genome sequencing effort. Surveillance consists of sample acquisition, library preparation, and whole genome sequencing. This has necessitated a classification scheme detailing Variants of Concern (VOC) and Variants of Interest (VOI), and the rapid expansion of bioinformatics tools for sequence analysis. These bioinformatic tools are means for major actionable results: maintaining quality assurance and checks, defining population structure, performing genomic epidemiology, and inferring lineage to allow reliable and actionable identification and classification. Additionally, the pandemic has required public health laboratories to reach high throughput proficiency in sequencing library preparation and downstream data analysis rapidly. However, both processes can be limited by a lack of a standardized sequence dataset. We identified six SARS-CoV-2 sequence datasets from recent publications, public databases and internal resources. In addition, we created a method to mine public databases to identify representative genomes for these datasets. Using this novel method, we identified several genomes as either VOI/VOC representatives or non-VOI/VOC representatives. To describe each dataset, we utilized a previously published datasets format, which describes accession information and whole dataset information. Additionally, a script from the same publication has been enhanced to download and verify all data from this study.
The lifetime of a power electronics package is, to a large extent, determined by the reliability of its bonded interfaces under major loading conditions in an operational environment. Based on the application-level requirements of the package, bonded interface material is selected based on the results obtained from accelerated tests such as thermal cycling and power cycling. In addition to evaluating the reliability through accelerated tests, it is important to consider the impact of other component layers on the thermomechanical performance of the interface material, both from a material and geometric perspective. The co-efficient of thermal expansion (CTE) mismatch introduced by the use of different materials within a package and its structural design plays a critical role in determining the interface material reliability. In this paper, we present the results of a parametric modeling study of a power electronics package under thermal cycling in which the materials and geometric design of the different component layers were varied with respect to a baseline design to understand their impact on the reliability of the interface material. We chose the volume-averaged strain energy density per cycle computed at the corner region of the interface material as the metric for the reliability comparisons. Our results indicate that in addition to the CTE mismatch, the stiffness of individual component layers has a major impact on reliability. Among the different baseplates that we studied, aluminum silicon-carbide baseplates offered superior reliability over their copper and aluminum counterparts. We also found that the magnitude of the impact of stiffness variation—amongst the adjacent layers—on the reliability of the solder joint is dependent on the inherent CTE mismatch between the layers.
The lifetime of a power electronics package is, to a large extent, determined by the reliability of its bonded interfaces under major loading conditions in an operational environment. Based on the application-level requirements of the package, bonded interface material is selected based on the results obtained from accelerated tests such as thermal cycling and power cycling. In addition to evaluating the reliability through accelerated tests, it is important to consider the impact of other component layers on the thermomechanical performance of the interface material, both from a material and geometric perspective. The co-efficient of thermal expansion (CTE) mismatch introduced by the use of different materials within a package and its structural design plays a critical role in determining the interface material reliability. In this paper, we present the results of a parametric modeling study of a power electronics package under thermal cycling in which the materials and geometric design of the different component layers were varied with respect to a baseline design to understand their impact on the reliability of the interface material. We chose the volume-averaged strain energy density per cycle computed at the corner region of the interface material as the metric for the reliability comparisons. Our results indicate that in addition to the CTE mismatch, the stiffness of individual component layers has a major impact on reliability. Among the different baseplates that we studied, aluminum silicon-carbide baseplates offered superior reliability over their copper and aluminum counterparts. We also found that the magnitude of the impact of stiffness variation—amongst the adjacent layers—on the reliability of the solder joint is dependent on the inherent CTE mismatch between the layers.
Optimizations of atomic positions belong to the most commonly performed tasks in electronic structure calculations. Many simulations like global minimum searches or characterizations of chemical reactions require performing hundreds or thousands of minimizations or saddle computations. To automatize these tasks, optimization algorithms must not only be efficient but also very reliable. Unfortunately, computational noise in forces and energies is inherent to electronic structure codes. This computational noise poses a severe problem to the stability of efficient optimization methods like the limited-memory Broyden–Fletcher–Goldfarb–Shanno algorithm. We here present a technique that allows obtaining significant curvature information of noisy potential energy surfaces. We use this technique to construct both, a stabilized quasi-Newton minimization method and a stabilized quasi-Newton saddle finding approach. We demonstrate with the help of benchmarks that both the minimizer and the saddle finding approach are superior to comparable existing methods.
Since March 2010 the preliminary design of the ITER ICRF Antennas have been developed by CYCLE, a consortium consisting of IPP (Garching), CCFE (Culham), CEA (Cadarache), Politecnico di Torino (Torino) and LPPERM/KMS (Brussels). This paper describes the steps taken to develop the present geometry of the triplet pair Strap Housing Module from a thermal and structural perspective, and shows the critical areas of the structure. Key issues are the manufacturability, (achieved by HIPing - Hot Isostatic Pressing), the ability to handle the radiating plasma thermal flux of 0.35 MW/m{sup 2}, the RF losses and the neutronic radiation. HIPing is necessary to achieve the complicated system of cooling channels inside the structure, which divides the coolant equally in order to supply each strap in the triplet with 1 l/s of water. The components have also to withstand the strong mechanical forces generated by plasma disruptions affecting all internal structures and the elevated design cooling water pressure of 5MPa. In order to maximise reliability, joints between different materials in the cooling water system have been kept to a minimum. Therefore, in the interests of fabricability and availability, the whole structure is manufactured out of stainless steel (316L(N)IG). The low conductivity of 316L(N)IG demands small wall thicknesses to avoid hot spots; however this reduces the mechanical strength. Consequently an in depth FEM analysis is presented, which was used to find and to improve the critical aspects of this important component and was the best means of finding the optimum between thermal and mechanical performance.
Stockpile stewardship requires the description of weapons performance without resorting to underground nuclear testing. In the earlier tests, selected isotopes were used as detectors, and recovered after irradiation. Aspects of nuclear device performance were inferred by comparing the measured isotopic ratios to those predicted from simulations. The reaction flows that produce the final isotopic distributions proceed through regions of the nuclear chart that include unstable nuclei. Presently, improved nuclear data input is required to reanalyze prior tests and to certify the stockpile's reliability and safety. Many important cross sections are unknown, as is shown in the example of the Yttrium reaction network (Figure 1). The relevant reactions include (n,2n), (n,n'), (n,gamma), (n,p) and other charged-particle emitting reactions. The cross sections have to be calculated or inferred from indirect measurements. In both cases, reliable optical models that are valid a few nucleons away from stability are needed. The UNEDF Nuclear Reaction activities address this need by combining nuclear-structure input from UNEDF structure calculations with modern reaction theory and large-scale computational capabilities to develop microscopic nucleon-nucleus optical potentials that can be extrapolated to unstable nuclei. In addition, the reaction calculation tools and optical models developed in this context are proving valuable for planning and interpreting indirect (surrogate) measurements of the required cross sections.
In the present work, the near edge X-ray absorption spectroscopy (NEXAFS) spectra at both C and N K-edges of pyridine, 2-fluoropyridine, and 2,6-difluoropyridine have been studied both experimentally and theoretically. From an electronic point of view, both transition potential density functional theory and time-dependent density functional theory approaches lead to reliable results provided that suitable basis sets and density functionals are employed. In this connection, the global hybrid B3LYP functional in conjunction with the EPR-III basis set appears particularly suitable after constant scaling of the band positions. For the N K-edge, vertical energies obtained at these levels and broadened by symmetric Gaussian distributions provide spectra in reasonable agreement with the experiment. Vibronic contributions further modulate the band-shapes leading to a better agreement with the experimental results, but are not strictly necessary for semi-quantitative investigations. On the other hand, vibronic contributions are responsible for strong intensity redistribution in the NEXAFS C K-edge spectra, and their inclusion is thus mandatory for a proper description of experiments. In this connection, the simple vertical gradient model is particularly appealing in view of its sufficient reliability and low computational cost. For more quantitative results, the more refined vertical Hessian approach can be employed, and its effectiveness has been improved thanks to a new least-squares fitting approach.
The efforts in this project were focused in three primary objectives, provided below: • To develop the needed microstructural processing and performance relationships and associated material models for specific constituents in fabricated weldments (such as the parent material, heat affected zone regions and weld metal), • Apply these metallurgical relationships through modeling of a welded pressure bearing power plant component subjected to cyclic operational conditions under both mechanical and thermal loading, and • Validate the model through novel structural feature and component tests. An attempt was made to optimize the heat treatment schedule to reduce or eliminate the presence of boron nitride (BN) in the matrix of the Grade 92 material. BN has been shown to be a primary nucleation source for creep damage. Proper evaluation of material behavior, and as assessed in this research, included the influence of geometric notches (such as those from poor fabrication or design) and metallurgical notches (such as weldments or poorly processed parent material). Better practice methodologies were implemented in the theoretical design of a high-temperature, critical component for cyclic operation. It was demonstrated that vessel testing can be reliably performed and provides a critical link between standard testing, structural assessment and validation. While the goal of an improved Grade 92 steel was not achieved in this project, there is considerable value in applying the documented testing and assessment methodologies to assess material and/or component behavior. This is especially true for materials and/or designs which are expected to be subjected to flexible operating modes. Design of a cycling state of the art thick wall steam header for long term operation with a state of the art ferritic steel was a task in this project. Implementation of that task showed how many factors including plant operating mode, creep strength, creep ductility, fatigue strength, and resistance to crack growth should be considered and evaluated to properly design high temperature components for long term service (200,000 operating hours or more) in a modern Fossil energy power plant. The current ASME non-nuclear boiler design code does not mandate this level of complex analysis nor does it provide all the metallurgical and mechanical properties data needed to perform this kind of design analysis.
The shifts of core-level binding energies can provide powerful information about the electronic structure of a material. Understanding the physical origin of these shifts for catalytically relevant oxides may provide important insight into their properties. This requires reliable theoretical methods which are able to relate the binding energy shifts to the electronic structure. In order to establish such a methodology, the CaO(100) surface to bulk core-level binding energy shifts have been studied with Hartree-Fock and density-functional theory methods using both cluster and periodic slab models. The shifts obtained from the different theoretical methods are compared with each other and with data from synchrotron-based x-ray photoelectron spectroscopy (XPS) measurements. With a common approximation for the slab model treatment of XPS, the predicted binding energy shifts are seriously in error. In this work, the origin of the error is identified as arising from a flawed treatment of the surface atom binding energies, and a method for correcting the failure is presented.
While the formalism of multiresolution analysis (MRA), based on wavelets and adaptive integral representations of operators, is actively progressing in electronic structure theory (mostly on the independent-particle level and, recently, second-order perturbation theory), the concepts of multiresolution and adaptivity can also be utilized within the traditional formulation of correlated (many-particle) theory which is based on second quantization and the corresponding (generally nonorthogonal) tensor algebra. In this paper, we present a formalism called scale-adaptive tensor algebra (SATA) which exploits an adaptive representation of tensors of many-body operators via the local adjustment of the basis set quality. Given a series of locally supported fragment bases of a progressively lower quality, we formulate the explicit rules for tensor algebra operations dealing with adaptively resolved tensor operands. The formalism suggested is expected to enhance the applicability and reliability of local correlated many-body methods of electronic structure theory, especially those directly based on atomic orbitals (or any other localized basis functions).
Knowledge of the electron density distribution in the solar corona put constraints on the magnetic field configurations for coronal modeling and on initial conditions for solar wind modeling. We work with polarized SOHO/LASCO-C2 images from the last two recent minima of solar activity (1996–1997 and 2008–2010), devoid of coronal mass ejections. The goals are to derive the 4D electron density distributions in the corona by applying a newly developed time-dependent tomographic reconstruction method and to compare the results between the two solar minima and with two magnetohydrodynamic models. First, we confirm that the values of the density distribution in thermodynamic models are more realistic than in polytropic ones. The tomography provides more accurate distributions in the polar regions, and we find that the density in tomographic and thermodynamic solutions varies with the solar cycle in both polar and equatorial regions. Second, we find that the highest-density structures do not always correspond to the predicted large-scale heliospheric current sheet or its helmet streamer but can follow the locations of pseudo-streamers. We deduce that tomography offers reliable density distributions in the corona, reproducing the slow time evolution of coronal structures, without prior knowledge of the coronal magnetic field over a full rotation. Finally, we suggest that the highest-density structures show a differential rotation well above the surface depending on how they are magnetically connected to the surface. Such valuable information on the rotation of large-scale structures could help to connect the sources of the solar wind to their in situ counterparts in future missions such as Solar Orbiter and Solar Probe Plus.
The mercury dihalides show a remarkable diversity in the structural preferences in their minimum energy structure types, spanning molecular to strongly bound ionic solids. A challenge in the development of density functional methods for extended systems is to arrive at strategies that serve equally well such a broad range of bonding modes or structural preferences. The chemical bonding and the stabilities of mercury dihalides and the general utility and reliability of the van der Waals density functional with C09 exchange (vdW-DF-C09) in predicting or describing the energetics and structural preferences in these metal dihalides is examined. Here, we show that, in contrast with the uncorrected generalized gradient approximation of the Perdew-Burke-Erzenhoff (PBE) exchange-correlation functional, qualitative and quantitative patterns in the bonding of the mercury dihalide solids are well reproduced with vdW-DF-C09 for the full series of HgX 2 systems for X = F, Cl, Br, and I. The possible existence of a low-temperature cotunnite polymorph for HgF 2 and PbF 2 is posited.
Dynamic structural systems either operational or under test require complex spectral analysis in order to characterize their modal responses. In some applications constant vigilance in terms of analysis or potential failures demand an accurate methodology to estimate both modal frequencies as well as mode shapes. Typically, classical periodogram spectra are windowed and averaged applying the well-known Welch periodogram methodology (WPM) for spectral estimation [1]. In this report we discuss an alternative technique—the Multitaper Method (MTM) that can be applied to solve this challenging problem especially in the case of noisy, uncertain accelerometer measurements [2]. This approach is based on the development of a set of orthogonal tapers or windows enabling a superior statistical performance alleviating the need to section/overlap the measured data. Therefore, the MTM is capable of producing reliable and accurate spectral estimates enabling the extraction of the structural modal vibration frequencies.
Moderate-to-large magnitude earthquakes represent a significant risk around the world and a major issue across the Department of Energy (DOE) mission interests, ranging from DOE's one-of-a-kind mission critical facilities to national energy systems (nuclear power plants, electric/gas distribution systems, and renewable energy production facilities). As a result, for the DOE there is a strong motivation in promoting and developing a risk-informed and performance-based approach for both the assessment of existing facilities and the design of new structures. DOE National Labs have historically made multidisciplinary contributions to the solution of challenging seismic hazard and risk problems and have been at the forefront for the development of design standards for major facilities. The current approach for estimating the ground motion hazard for critical facilities is based upon a Probabilistic Seismic Hazard Analysis (PSHA), whereby measured ground motion data are collected and homogenized into a probabilistic framework to estimate the hazard for future events. However, recent studies have highlighted the site-specific nature of earthquakes and shown the limitations of a purely empirical approach. Moreover, the design of infrastructure systems has been moving toward a performance-based approach, in which specific performance goals are specified at different levels of excitation and for different limit states. Such an approach requires the availability of robust and reliable computational tools for the nonlinear dynamic analysis of complex structural systems. Ongoing research at Lawrence Berkeley National Laboratory (LBNL) and Lawrence Livermore National Laboratory (LLNL) aims at developing an advanced simulation environment at the exascale for a coupled assessment of earthquake hazard (ground motions) and earthquake risk (structural demands)
Corrections of order {alpha}{sup 5} and {alpha}{sup 6} to the hyperfine structure of S- and P-wave energy levels of the muonic-helium ion are calculated. Electron-vacuum-polarization effects, corrections for the nuclear structure, and recoil effects are taken into account. The numerical values obtained for respective hyperfine splitting, -1334.73 meV (1S), -166.64 meV (2S), -58 712.90 {mu}eV (2P{sub 1/2}), and -24 290.69 {mu}eV (2P{sub 3/2}), can be viewed as a reliable estimate for a comparison with experimental data, and the hyperfine-structure interval of {Delta}{sub 12} = 8{Delta}E{sup hfs}(2S) - {Delta}E{sup hfs}(1S) = 1.59 meV can be used to test QED predictions.
The reliability of photovoltaics is commonly studied at the module level. Many reliability problems originate from module attributes, such as metal interconnections to cells, junction boxes. However, significant work in reliability can also be done prior to module design. Testing for reliability earlier in the research cycle increases the probability of avoiding common module reliability problems before cell changes are implemented on a large scale. Cell-level reliability studies can thus lower the rates of module failures in the field and provide confidence to investors that new technologies will perform as advertised. This report summarizes how we investigated three reliability concerns in Cu(In,Ga)Se 2 (CIGS) photovoltaics at the cell level: metastability, shading-induced damage, and potential-induced degradation (PID). We find that examining these concerns required developing robust measurement protocols including the fabrication of novel testing structures. This information will allow readers to incorporate sound metrics for investigating reliability phenomena and aid their studies of cell and module reliability improvements.
Purpose: Precise contour delineation of tumor targets and critical structures from CT simulations is essential for accurate radiotherapy (RT) treatment planning. However, manual and automatic delineation processes can be error prone due to limitations in imaging techniques and individual anatomic variability. Tedious and laborious manual verification is hence needed. This study develops a general framework for automatically assessing RT contours for head-neck cancer patients using geometric attribute distribution models (GADMs). Methods: Geometric attributes (centroid and volume) were computed from physician-approved RT contours of 29 head-neck patients. Considering anatomical correlation between neighboring structures, the GADM for each attribute was trained to characterize intra- and interpatient structure variations using principal component analysis. Each trained GADM was scalable and deformable, but constrained by the principal attribute variations of the training contours. A new hierarchical model adaptation algorithm was utilized to assess the RT contour correctness for a given patient. Receiver operating characteristic (ROC) curves were employed to evaluate and tune system parameters for the training models. Results: Experiments utilizing training and non-training data sets with simulated contouring errors were conducted to validate the framework performance. Promising assessment results of contour normality/abnormality for the training contour-based data were achieved with excellent accuracy (0.99), precision (0.99), recall (0.83), and F-score (0.97), while corresponding values of 0.84, 0.96, 0.83, and 0.9 were achieved for the non-training data. Furthermore, the areas under the ROC curves were above 0.9, validating the accuracy of this test. Conclusion: The proposed framework can reliably identify contour normality/abnormality based upon intra- and inter-structure constraints derived from clinically-approved contours. It also allows physicians to analytically determine the system parameters to fit various clinic requirements (e.g. as-low-as-possible false positives). It has great potential for improving RT work flow. More geometric attributes and training sets will be investigated to improve framework performance in the future.