QMCPACK v3.16.0
QMCPACK is an open-source production level many-body ab initio Quantum Monte Carlo code for computing the electronic structure of atoms, molecules, and solids with full performance portable GPU support.
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.
QMCPACK is an open-source production level many-body ab initio Quantum Monte Carlo code for computing the electronic structure of atoms, molecules, and solids with full performance portable GPU support.
QMCPACK is an open-source production level many-body ab initio Quantum Monte Carlo code for computing the electronic structure of atoms, molecules, and solids with full performance portable GPU support.
QMCPACK is an open-source production level many-body ab initio Quantum Monte Carlo code for computing the electronic structure of atoms, molecules, and solids with full performance portable GPU support.
QMCPACK is an open-source production level many-body ab initio Quantum Monte Carlo code for computing the electronic structure of atoms, molecules, and solids with full performance portable GPU support.
The goal of the DOE OE Energy Storage System Safety Roadmap is to foster confidence in the safety and reliability of energy storage systems. There are three interrelated objectives to support the realization of that goal: research, codes and standards(C/S)and communication/coordination. The objective focused on C/S is "To apply research and development to support efforts that refocused on ensuring that codes and standards are available to enable the safe implementation of energy storage systems in a comprehensive, non-discriminatory and science-based manner." The following activities support that objective and realization of the goal: (a.) Review and assess C/S which affect the design, installation, and operation of energy storage systems(ESS). (b.) Identify gaps in knowledge that require research and analysis that can serve as a basis for criteria in those C/S. (c.) Identify areas in C/S that are potentially in need of revision or enhancement and can benefit from activities conducted under research and development. (d.) Develop input for new or revisions to existing C/S through individual stakeholders, facilitated task forces, or through laboratory staff supporting these efforts. The purpose of this Codes and Standards Update is to support the above activities by providing information on efforts being conducted by U.S. standards developing organizations(SDOs)and other entities that are focused on ESS safety. The information is organized relative to the scope of each document in relation to ESS from the "macro to the micro"(e.g., from overarching covering considerable scope ,to installation specific, to ESS and then ESS components). Note that more macro documents are also likely to adopt by reference more micro documents.
Optimized operation of fusion devices demands detailed understanding of plasma transport, a problem that must be addressed with advances in both measurement and data analysis techniques. In this work, we adopt Bayesian inference methods to determine experimental particle transport, leveraging opportunities from high-resolution He-like ion spectra in a tokamak plasma. The Bayesian spectral fitting code is used to analyze resonance (w), forbidden (z), intercombination (x, y), and satellite (k, j) lines of He-like Ca following laser blow-off injections on Alcator C-Mod. This offers powerful transport constraints since these lines depend differently on electron temperature and density, but also differ in their relation to Li-like, He-like, and H-like ion densities, often the dominant Ca charge states over most of the C-Mod plasma radius. Using synthetic diagnostics based on the AURORA package, we demonstrate improved effectiveness of impurity transport inferences when spectroscopic data from a progressively larger number of lines are included.
Optimized operation of fusion devices demands detailed understanding of plasma transport, a problem that must be addressed with advances in both measurement and data analysis techniques. In this work, we adopt Bayesian inference methods to determine experimental particle transport, leveraging opportunities from high-resolution He-like ion spectra in a tokamak plasma. The Bayesian spectral fitting code is used to analyze resonance (w), forbidden (z), intercombination (x, y), and satellite (k, j) lines of He-like Ca following laser blow-off injections on Alcator C-Mod. This offers powerful transport constraints since these lines depend differently on electron temperature and density, but also differ in their relation to Li-like, He-like, and H-like ion densities, often the dominant Ca charge states over most of the C-Mod plasma radius. Using synthetic diagnostics based on the AURORA package, we demonstrate improved effectiveness of impurity transport inferences when spectroscopic data from a progressively larger number of lines are included.
Abstract Increased interest in compact heat exchangers (CHXs) to serve as intermediate heat exchangers of very high temperature reactors resulted in significant research and development on their design, analysis, and construction. Printed circuit heat exchangers are a type of CHXs with high thermal efficiency and compactness achieved through diffusion bonding a stack of etched plates with millimeter scaled channels. The diffusion bonding process changes the microstructural and mechanical properties of the wrought metal plates. The current nonnuclear design code ASME section VIII, division 1 captures the material property change through a “joint efficiency factor.” However, the current nuclear design code ASME section III, division 5 does not address or support the diffusion bonded material properties. Hence, there is a need to develop allowable stresses, isochronous curves, and fatigue life curves for various diffusion bonded alloys. In this study, Alloy 800H material was selected to establish the diffusion bonded material properties under tension, creep, fatigue, and creep-fatigue loads at elevated temperatures in the range 550–760 °C. A set of tests on diffusion bonded Alloy 800H (DB 800H) were performed and the acquired data are used in developing allowable stresses Sy, Su, Sr, Sm, St, Smt, So, isochronous curves and fatigue life curves according to the ASME section III, division 5 requirements. This paper also presents detailed procedures used in developing the ASME code section III division 5 design provisions for diffusion bonded Alloy 800H.
Here, we have begun a preliminary investigation into the opacity of Sn at low temperatures (< 50 eV). The emissivity and opacity of Sn is a crucial factor in determining the utility of Sn in EUV lithography, with numerous industrial implications. To this end, we have been exploring the accuracy of some approximations used in opacity models for the relevant ion stages of Sn (neutral through ~ 18 times ionized). We also find that the use of intermediate-coupling, as compared to full configuration-interaction, is not adequate to obtain accurate line positions of the important bound-bound transitions in Sn. One requires full configuration-interaction to properly describe the strong mixing between the various n=4 sub-shells that give rise to the Δn= 0 transitions that dominate the opacity spectrum at low temperatures. Furthermore, since calculations that include full configuration-interaction for large numbers of configurations quickly become computationally prohibitive, we have explored hybrid calculations, in which full configuration-interaction is retained for the most important transitions, while intermediate-coupling is employed for all other transitions. After extensive exploration of the atomic structure properties, local-thermodynamic-equilibrium (LTE) opacities are generated using the ATOMIC code at selected temperatures and densities and compared to experiment.
Magnetic confinement of plasmas is widely studied for restricting a plasma's extent to a given region. Presented is one scheme for what is referred to as magnetic plasma expulsion, which would prevent the plasma from entering a given region. Magnetic plasma expulsion techniques may allow instrumentation access to the interiors of plasmas. Magnetic plasma expulsion techniques may also make alternative approaches for magnetic plasma confinement possible. Particle-in-cell simulation with the Warp code is used to study magnetic plasma expulsion. Some conditions for achieving magnetic plasma expulsion are reported.
Using temporally shaped laser pulses, multiple shocks can be launched in direct-drive inertial confinement fusion implosion experiments to set the shell on a desired isentrope or adiabat. The velocity of the first shock and the times at which subsequent shocks catch up to it are measured through the VISAR diagnostic on OMEGA. Simulations reproduce these velocity and shock-merger time measurements when using laser pulses designed for setting mid-adiabat (α ~ 3) implosions, but agreement degrades for lower-adiabat (α ~ 1) designs. Several possibilities for this difference are studied: errors in placing the target at the center of irradiation (target offset), variations in energy between the different incident beams (power imbalance), and errors in modeling the laser energy coupled into the capsule. Simulation results indicate that shock timing is most sensitive to details of the density and temperature profiles in the coronal plasma, which influences the laser energy coupled into the target, and only marginally sensitive to target offset and beam power imbalance. A new technique under development to infer coronal profiles using x-ray self-emission imaging can be applied to the pulse shapes used in shock-timing experiments. In conclusion, this will help identify improved physics models to implement in codes and consequently enhance shock-timing predictive capability for low-adiabat pulses.
An experimental scan of the electron cyclotron waves (EC) toroidal injection angle in plasma breakdown is performed at the DIII-D tokamak. The second harmonic, extraordinary mode EC is used for the study. The dependence of n e and T e on the EC injection angle cannot be conclusively drawn from this study due to the large error bars in the n e and T e measurements. On the other hand, high T e data points are observed in some discharges which can be explained by nonlinear heating. The D α emission measurement shows a clear relation between the breakdown time and the injection angle. An experimental investigation of the cause of the dependence of breakdown delay on the EC injection angle suggests that when the injection angle is oblique, the EC heating after the reflection at the inboard wall may become ineffective and cause the breakdown delay even when the EC heating before and directly upon injection remains effective. A preliminary run of the heat and transport balance code DYON indicates that the obtained dataset is suitable for a quantitative validation of EC absorption models.
This report assesses the ASME Section III, Division 5, Subsection HB, Subpart B rules covering the design and construction of high temperature Class A nuclear reactor components for their robustness against creep notch effects. The creep notch effect combines the effect of multiaxial stresses on material creep deformation and damage. This report considers both effects, including the potential for creep mechanism shifts affecting the extrapolation of design data from high stress, short-time experimental data to low stress, long time operating component conditions. The report summarizes the history of and literature on multiaxial creep and surveys the current ASME design rules dealing with multiaxial effects. The report then describes the results of two dedicated numerical studies, one assessing the robustness of the ASME rules against uncertainty in extrapolating from uniaxial creep test data to multiaxial component conditions and a second study examining the potential effects of a mechanism shift from high stress, dislocation-mediated creep to diffusion-dominated creep at lower stresses. The final conclusion of the report is that the ASME rules adequately guard against multiaxial creep failure, though there are several aspects of the Code that could be optimized to provide less over-conservative design predictions or to provide a more consistent design margin as a function of temperature and stress. In a few areas, particularly the potential for mechanism shifts outside the currently-available experimental data, the Code rules should be revaluated as additional experimental data and new modeling and simulation results become available
The purpose of this project is to develop an improved Navy-wide approach for estimating risk and calculating cleanup goals for radiologically impacted environmental media (e.g., soils). This comparative analysis focused on the new features added (or improved) to the PRG/DCC Calculators since the 2015 comparison study (Yu et al. 2015a). RESRAD-ONSITE code Version 7.2 and the PRG/DCC Calculators available on the USEPA website between January and June of 2024 were used.
Version 2.0 of the RESRAD-RDD&IND computer code is designed to support the implementation of protective action guides (PAGs) after a nuclear emergency incident including a radiological dispersal device (RDD) and/or an improvised nuclear device (IND) incident (EPA 2017). Eight different group types, addressing various decisions, are available for selection. The RESRAD-RDD&IND code calculates radiological doses, stay times, etc., for the selected group that the user wishes to focus on. (That is, the results for all the groups are not calculated simultaneously, and the input for those other groups do not matter, although some parameter values are shared between groups.) Version 2.0 has a user-friendly interface so that the RESRAD-RDD&IND code can be used with minimal training. For example, the user can select the major characteristics of the problem-event type, source term, and decision type from the left side of the interface and then calculate the results with the default assumptions for the exposure scenarios. More in-depth analysis would include specifying site-specific exposure scenario characteristics in the right side of the interface. The procedures for data entry and results viewing are self-explanatory. This is because common window maneuvering features and text instructions were incorporated in the interface design. General and context-specific help are available to aid users entering parameter values, as well. The RESRAD-RDD&IND computer code gives the user the option to select either an RDD or IND incident for analysis. For an RDD event analysis, 11 radionuclides (Am-241, Cf-252, Cm-244, Co-60, Cs-137, Ir-192, Po-210, Pu-238, Pu-239, Ra-226, and Sr-90) are included. These 11 radionuclides are the radionuclides most likely used for an RDD. More than 90 radionuclides can be selected for an IND event analysis. Initial default concentrations are provided for 44 radionuclides for a uranium-fueled IND event. These 44 radionuclides are those that would contribute significantly to the radiation dose associated with a uranium-fueled bomb detonation. The radionuclides generated from ingrowth of these 44 initial radionuclides are also automatically included in the analysis. Pu-239, Cs-134m, Ru-105, and Rb-89 and their progeny can be selected for analysis if they are detected and their concentrations are determined. This user’s guide, which is Volume 2 of the User’s Manual for RESRAD-RDD&IND Code Version 2, provides instructions to users on how to install the RESRAD-RDD&IND code, navigate the interface, and use the various features, including those discussed above, to set up an analysis and view/print the results in text outputs. Volume 1 of the User’s Manual for RESRAD-RDD&IND Code Version 2 (Yu et al. 2026), which contains descriptions of the methodology and theoretical basis for dose modeling and the mathematical equations implemented in the code, can be accessed and viewed through the Help menu in the code or can be downloaded from the RESRAD website (https://resrad.evs.anl.gov).
This dataset contains solid-state 13C NMR data and atomistic molecular dynamics simulation files supporting the study of nanoscale secondary cell wall architecture across 13 genetically diverse Populus trichocarpa genotypes grown under uniform greenhouse conditions in 13C-enriched CO2 atmospheres (~89% 13C enrichment).The dataset contains two collections of solid-state 13C NMR data. (1) 200 MHz data (Bruker Avance III HD, 4 mm HX probe, 10 kHz MAS): raw Bruker TopSpin experiment folders and DMFIT-exported ascii spectra for selective and non-selective 1D 13C-13C spin diffusion experiments (3000 ms mixing) used to quantify inter-polymer spatial proximities, and short-mixing (1 ms) reference spectra used for polymeric abundance quantification by spectral deconvolution. (2) 600 MHz data (Bruker Avance III, 1.6 mm PhoenixNMR HXY probe, 30 kHz MAS): raw Bruker TopSpin experiment folders containing 2D CORD, 2D CP-INADEQUATE, and 13C/1H relaxation (T1, T1rho) experiments for all 13 genotypes, with processed Excel workbooks per experiment type. Molecular dynamics simulation code, coordinate files, and analysis scripts (NAMD/CHARMM/Python) for six atomistic cell wall models are included. Summarized ssNMR data are compiled into a single excel file and subjected to statistical analysis. Multivariate analysis code (PCA, Pearson correlation) and summary data are provided as excel worksheets and Jupyter notebooks (Python 3).
A key goal of exoplanet spectroscopy is to measure atmospheric properties, such as abundances of chemical species, in order to connect them to our understanding of atmospheric physics and planet formation. In this new era of high-quality JWST data, it is paramount that these measurement methods are robust. When comparing atmospheric models to observations, multiple candidate models may produce reasonable fits to the data. Typically, conclusions are reached by selecting the best-performing model according to some metric. This ignores model uncertainty in favor of specific model assumptions, potentially leading to measured atmospheric properties that are overconfident and/or incorrect. In this paper, we compare three ensemble methods for addressing model uncertainty by combining posterior distributions from multiple analyses: Bayesian model averaging, a variant of Bayesian model averaging using leave-one-out predictive densities, and stacking of predictive distributions. We demonstrate these methods by fitting the Hubble Space Telescope (HST) + Spitzer transmission spectrum of the hot Jupiter HD 209458b using models with different cloud and haze prescriptions. All of our ensemble methods lead to uncertainties on retrieved parameters that are larger but more realistic and consistent with physical and chemical expectations. Since they have not typically accounted for model uncertainty, uncertainties of retrieved parameters from HST spectra have likely been underreported. We recommend stacking as the most robust model combination method. Our methods can be used to combine results from independent retrieval codes and from different models within one code. They are also widely applicable to other exoplanet analysis processes, such as combining results from different data reductions.
We present a new code for astrophysical magnetohydrodynamics specifically designed and optimized for high performance and scaling on modern and future supercomputers. We describe a novel hybrid OpenMP/MPI programming model that emerged from a collaboration between Cray, Inc. and the University of Minnesota. This design utilizes MPI-RMA optimized for thread scaling, which allows the code to run extremely efficiently at very high thread counts ideal for the latest generation of multi-core and many-core architectures. Such performance characteristics are needed in the era of “exascale” computing. We describe and demonstrate our high-performance design in detail with the intent that it may be used as a model for other, future astrophysical codes intended for applications demanding exceptional performance.