Search NASA⌕ Search

SEARCH · Search NASA

Results for “CALCULATORS”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Full dimensional Franck-Condon factors for the acetylene A{sup ~} {sup 1}A{sub u}—X{sup ~1}Σ{sup +}{sub g} transition. I. Method for calculating polyatomic linear—bent vibrational intensity factors and evaluation of calculated intensities for the gerade vibrational modes in acetylene

Franck-Condon vibrational overlap integrals for the A{sup ~1}A{sub u}—X{sup ~1}Σ{sup +}{sub g} transition in acetylene have been calculated in full dimension in the harmonic normal mode basis. The calculation uses the method of generating functions first developed for polyatomic Franck-Condon factors by Sharp and Rosenstock [J. Chem. Phys. 41(11), 3453–3463 (1964)], and previously applied to acetylene by Watson [J. Mol. Spectrosc. 207(2), 276–284 (2001)] in a reduced-dimension calculation. Because the transition involves a large change in the equilibrium geometry of the electronic states, two different types of corrections to the coordinate transformation are considered to first order: corrections for axis-switching between the Cartesian molecular frames and corrections for the curvilinear nature of the normal modes at large amplitude. The angular factor in the wave function for the out-of-plane component of the trans bending mode, ν{sub 4}{sup ″}, is treated as a rotation, which results in an Eckart constraint on the polar coordinates of the bending modes. To simplify the calculation, the other degenerate bending mode, ν{sub 5}{sup ″}, is integrated in the Cartesian basis and later transformed to the constrained polar coordinate basis, restoring the conventional v and l quantum numbers. An updated A{sup ~}-state harmonic force field obtained recently in the R. W. Field research group is evaluated. The results for transitions involving the gerade vibrational modes are in qualitative agreement with experiment. Calculated results for transitions involving ungerade modes are presented in Paper II of this series [G. B. Park, J. H. Baraban, and R. W. Field, “Full dimensional Franck–Condon factors for the acetylene A{sup ~1}A{sub u}—X{sup ~1}Σ{sup +}{sub g} transition. II. Vibrational overlap factors for levels involving excitation in ungerade modes,” J. Chem. Phys. 141, 134305 (2014)].

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

Do Practical Standard Coupled Cluster Calculations Agree Better than Kohn–Sham Calculations with Currently Available Functionals When Compared to the Best Available Experimental Data for Dissociation Energies of Bonds to 3d Transition Metals?

Coupled-cluster (CC) methods have been extensively used as the high-level approach in quantum electronic structure theory to predict various properties of molecules when experimental results are unavailable. It is often assumed that CC methods, if they include at least up to connected-triple-excitation quasiperturbative corrections to a full treatment of single and double excitations (in particular, CCSD(T)), and a very large basis set, are more accurate than Kohn–Sham (KS) density functional theory (DFT). In the present work, we tested and compared the performance of standard CC and KS methods on bond energy calculations of 20 3d transition metal-containing diatomic molecules against the most reliable experimental data available, as collected in a database called 3dMLBE20. It is found that, although the CCSD(T) and higher levels CC methods have mean unsigned deviations from experiment that are smaller than most exchange-correlation functionals for metal–ligand bond energies of transition metals, the improvement is less than one standard deviation of the mean unsigned deviation. Furthermore, on average, almost half of the 42 exchange-correlation functionals that we tested are closer to experiment than CCSD(T) with the same extended basis set for the same molecule. The results show that, when both relativistic and core–valence correlation effects are considered, even the very high-level (expensive) CC method with single, double, triple, and perturbative quadruple cluster operators, namely, CCSDT(2)Q, averaged over 20 bond energies, gives a mean unsigned deviation (MUD(20) = 4.7 kcal/mol when one correlates only valence, 3p, and 3s electrons of transition metals and only valence electrons of ligands, or 4.6 kcal/mol when one correlates all core electrons except for 1s shells of transition metals, S, and Cl); and that is similar to some good xc functionals (e.g., B97-1 (MUD(20) = 4.5 kcal/mol) and PW6B95 (MUD(20) = 4.9 kcal/mol)) when the same basis set is used. We found that, for both coupled cluster calculations and KS calculations, the T1 diagnostics correlate the errors better than either the M diagnostics or the B1 DFT-based diagnostics. The potential use of practical standard CC methods as a benchmark theory is further confounded by the finding that CC and DFT methods usually have different signs of the error. We conclude that the available experimental data do not provide a justification for using conventional single-reference CC theory calculations to validate or test xc functionals for systems involving 3d transition metals.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2015

Improved methods for Feynman path integral calculations and their application to calculate converged vibrational–rotational partition functions, free energies, enthalpies, entropies, and heat capacities for methane

We present an improved version of our “path-by-path” enhanced same path extrapolation scheme for Feynman path integral (FPI) calculations that permits rapid convergence with discretization errors ranging from O(P{sup −6}) to O(P{sup −12}), where P is the number of path discretization points. We also present two extensions of our importance sampling and stratified sampling schemes for calculating vibrational–rotational partition functions by the FPI method. The first is the use of importance functions for dihedral angles between sets of generalized Jacobi coordinate vectors. The second is an extension of our stratification scheme to allow some strata to be defined based only on coordinate information while other strata are defined based on both the geometry and the energy of the centroid of the Feynman path. These enhanced methods are applied to calculate converged partition functions by FPI methods, and these results are compared to ones obtained earlier by vibrational configuration interaction (VCI) calculations, both calculations being for the Jordan–Gilbert potential energy surface. The earlier VCI calculations are found to agree well (within ∼1.5%) with the new benchmarks. The FPI partition functions presented here are estimated to be converged to within a 2σ statistical uncertainty of between 0.04% and 0.07% for the given potential energy surface for temperatures in the range 300–3000 K and are the most accurately converged partition functions for a given potential energy surface for any molecule with five or more atoms. We also tabulate free energies, enthalpies, entropies, and heat capacities.

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

Methods to Calculate Electronic Excited-State Dynamics for Molecules on Large Metal Clusters with Many States: Ensuring Fast Overlap Calculations and a Robust Choice of Phase

Here, we present an efficient set of methods for propagating excited-state dynamics involving a large number of configuration interaction singles (CIS) or Tamm-Dancoff approximation (TDA) single-reference excited states. Specifically, (i) following Head-Gordon et al., we implement an exact evaluation of the overlap of singly-excited CIS/TDA electronic states at different nuclear geometries using a biorthogonal basis and (ii) we employ a unified protocol for choosing the correct phase for each adiabat at each geometry. For many-electron systems, the combination of these techniques significantly reduces the computational cost of integrating the electronic Schrodinger equation and imposes minimal overhead on top of the underlying electronic structure calculation. As a demonstration, we calculate the electronic excited-state dynamics for a hydrogen molecule scattering off a silver metal cluster, focusing on high-lying excited states, where many electrons can be excited collectively and crossings are plentiful. Interestingly, we find that the high-lying, plasmon-like collective excitation spectrum changes with nuclear dynamics, highlighting the need to simulate non-adiabatic nuclear dynamics and plasmonic excitations simultaneously. In the future, the combination of methods presented here should help theorists build a mechanistic understanding of plasmon-assisted charge transfer and excitation energy relaxation processes near a nanoparticle or metal surface.

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

Inverse load calculation procedure for offshore wind turbines and application to a 5‐MW wind turbine support structure

Abstract A significant number of wind turbines installed today have reached their designed service life of 20 years, and the number will rise continuously. Most of these turbines promise a more economical performance if they operate for more than 20 years. To assess a continued operation, we have to analyze the load‐bearing capacity of the support structure with respect to site‐specific conditions. Such an analysis requires the comparison of the loads used for the design of the support structure with the actual loads experienced. This publication presents the application of a so‐called inverse load calculation to a 5‐MW wind turbine support structure. The inverse load calculation determines external loads derived from a mechanical description of the support structure and from measured structural responses. Using numerical simulations with the software fast , we investigated the influence of wind‐turbine‐specific effects such as the wind turbine control or the dynamic interaction between the loads and the support structure to the presented inverse load calculation procedure. fast is used to study the inverse calculation of simultaneously acting wind and wave loads, which has not been carried out until now. Furthermore, the application of the inverse load calculation procedure to a real 5‐MW wind turbine support structure is demonstrated. In terms of this practical application, setting up the mechanical system for the support structure using measurement data is discussed. The paper presents results for defined load cases and assesses the accuracy of the inversely derived dynamic loads for both the simulations and the practical application. Copyright © 2017 John Wiley & Sons, Ltd.

Pahn, T.↗
DOE OSTI2017

A novel approach to evaluate soil heat flux calculation: An analytical review of nine methods

Abstract There are no direct methods to evaluate calculated soil heat flux (SHF) at the surface ( G 0 ). Instead, validation and cross evaluation of methods for calculating G 0 usually rely on the conventional calorimetric method or the degree of the surface energy balance closure. However, there is uncertainty in the calorimetric method itself, and factors apart from G 0 also contribute to nonclosure of the surface energy balance. Here we used a novel approach to evaluate nine different methods for calculating SHF, including the calorimetric method and methods based on analytical solutions of the heat diffusion equation. The SHF ( G z ) measured by a self‐calibrating SHF plate at a depth of z = 5 cm below the surface (hereafter G m _5 cm ) was deployed as a reference. Each SHF calculation method was assessed by comparing the calculated G z at the same depth (hereafter G c _5 cm ) with G m _5 cm . The calorimetric method and simple measurement method performed best in determining G c _5 cm but still underestimated G m _5 cm by 19% during the daytime. Possible causes for this underestimation include errors and uncertainties in SHF measurements and soil thermal properties, as well as the phase lag between G c _5 cm and G m _5 cm . Our results indicate that the calorimetric method achieves the most accurate SHF estimates if self‐calibrating SHF plates are deployed at two depths (e.g., 5 cm and 10 cm), soil temperature and water content measurements are made in a few depths between the two plates, and soil thermal properties are accurately quantified.

Gao, Zhongming↗
DOE OSTI2016

Challenges in large scale quantum mechanical calculations: Challenges in large scale quantum mechanical calculations

During the past decades, quantum mechanical methods have undergone an amazing transition from pioneering investigations of experts into a wide range of practical applications, made by a vast community of researchers. First principles calculations of systems containing up to a few hundred atoms have become a standard in many branches of science. The sizes of the systems which can be simulated have increased even further during recent years, and quantum-mechanical calculations of systems up to many thousands of atoms are nowadays possible. This opens up new appealing possibilities, in particular for interdisciplinary work, bridging together communities of different needs and sensibilities. In this review we will present the current status of this topic, and will also give an outlook on the vast multitude of applications, challenges and opportunities stimulated by electronic structure calculations, making this field an important working tool and bringing together researchers of many different domains.

Ratcliff, Laura E.↗
DOE OSTI2016

Calculation of Debye-Scherrer diffraction patterns from highly stressed polycrystalline materials

Calculations of Debye-Scherrer diffraction patterns from polycrystalline materials have typically been done in the limit of small deviatoric stresses. Although these methods are well suited for experiments conducted near hydrostatic conditions, more robust models are required to diagnose the large strain anisotropies present in dynamic compression experiments. A method to predict Debye-Scherrer diffraction patterns for arbitrary strains has been presented in the Voigt (iso-strain) limit [Higginbotham, J. Appl. Phys. 115, 174906 (2014)]. Here, we present a method to calculate Debye-Scherrer diffraction patterns from highly stressed polycrystalline samples in the Reuss (iso-stress) limit. This analysis uses elastic constants to calculate lattice strains for all initial crystallite orientations, enabling elastic anisotropy and sample texture effects to be modeled directly. The effects of probing geometry, deviatoric stresses, and sample texture are demonstrated and compared to Voigt limit predictions. An example of shock-compressed polycrystalline diamond is presented to illustrate how this model can be applied and demonstrates the importance of including material strength when interpreting diffraction in dynamic compression experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Solid-state dimer method for calculating solid-solid phase transitions

The dimer method is a minimum mode following algorithm for finding saddle points on a potential energy surface of atomic systems. Here, the dimer method is extended to include the cell degrees of freedom for periodic solid-state systems. Using this method, reaction pathways of solid-solid phase transitions can be determined without having to specify the final state structure or reaction mechanism. Example calculations include concerted phase transitions between CdSe polymorphs and a nucleation and growth mechanism for the A15 to BCC transition in Mo.

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

Bandgap calculations and trends of organometal halide perovskites

Energy production from the Sun requires a stable efficient light absorber. Promising candidates in this respect are organometal perovskites (ABX{sub 3}), which have been intensely investigated during the last years. Here, we have performed electronic structure calculations of 240 perovskites composed of Cs, CH{sub 3}NH{sub 3}, and HC(NH{sub 2}){sub 2} as A-cation, Sn and Pb as B-ion, and a combination of Cl, Br, and I as anions. The calculated gaps span over a region from 0.5 to 5.0 eV. In addition, the trends over bandgaps have been investigated: the bandgap increases with an increase of the electronegativities of the constituent species, while it reduces with an increase of the lattice constants of the system.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Monte Carlo calculations for reference dosimetry of electron beams with the PTW Roos and NE2571 ion chambers

Purpose: To investigate recommendations for reference dosimetry of electron beams and gradient effects for the NE2571 chamber and to provide beam quality conversion factors using Monte Carlo simulations of the PTW Roos and NE2571 ion chambers. Methods: The EGSnrc code system is used to calculate the absorbed dose-to-water and the dose to the gas in fully modeled ion chambers as a function of depth in water. Electron beams are modeled using realistic accelerator simulations as well as beams modeled as collimated point sources from realistic electron beam spectra or monoenergetic electrons. Beam quality conversion factors are calculated with ratios of the doses to water and to the air in the ion chamber in electron beams and a cobalt-60 reference field. The overall ion chamber correction factor is studied using calculations of water-to-air stopping power ratios. Results: The use of an effective point of measurement shift of 1.55 mm from the front face of the PTW Roos chamber, which places the point of measurement inside the chamber cavity, minimizes the difference betweenR{sub 50}, the beam quality specifier, calculated from chamber simulations compared to that obtained using depth-dose calculations in water. A similar shift minimizes the variation of the overall ion chamber correction factor with depth to the practical range and reduces the root-mean-square deviation of a fit to calculated beam quality conversion factors at the reference depth as a function of R{sub 50}. Similarly, an upstream shift of 0.34 r{sub cav} allows a more accurate determination of R{sub 50} from NE2571 chamber calculations and reduces the variation of the overall ion chamber correction factor with depth. The determination of the gradient correction using a shift of 0.22 r{sub cav} optimizes the root-mean-square deviation of a fit to calculated beam quality conversion factors if all beams investigated are considered. However, if only clinical beams are considered, a good fit to results for beam quality conversion factors is obtained without explicitly correcting for gradient effects. The inadequacy of R{sub 50} to uniquely specify beam quality for the accurate selection of k{sub Q} factors is discussed. Systematic uncertainties in beam quality conversion factors are analyzed for the NE2571 chamber and amount to between 0.4% and 1.2% depending on assumptions used. Conclusions: The calculated beam quality conversion factors for the PTW Roos chamber obtained here are in good agreement with literature data. These results characterize the use of an NE2571 ion chamber for reference dosimetry of electron beams even in low-energy beams.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2017

Adaptive local basis set for Kohn–Sham density functional theory in a discontinuous Galerkin framework II: Force, vibration, and molecular dynamics calculations

Recently, we have proposed the adaptive local basis set for electronic structure calculations based on Kohn–Sham density functional theory in a pseudopotential framework. The adaptive local basis set is efficient and systematically improvable for total energy calculations. In this paper, we present the calculation of atomic forces, which can be used for a range of applications such as geometry optimization and molecular dynamics simulation. We demonstrate that, under mild assumptions, the computation of atomic forces can scale nearly linearly with the number of atoms in the system using the adaptive local basis set. We quantify the accuracy of the Hellmann–Feynman forces for a range of physical systems, benchmarked against converged planewave calculations, and find that the adaptive local basis set is efficient for both force and energy calculations, requiring at most a few tens of basis functions per atom to attain accuracies required in practice. Since the adaptive local basis set has implicit dependence on atomic positions, Pulay forces are in general nonzero. However, we find that the Pulay force is numerically small and systematically decreasing with increasing basis completeness, so that the Hellmann–Feynman force is sufficient for basis sizes of a few tens of basis functions per atom. We verify the accuracy of the computed forces in static calculations of quasi-1D and 3D disordered Si systems, vibration calculation of a quasi-1D Si system, and molecular dynamics calculations of H{sub 2} and liquid Al–Si alloy systems, where we show systematic convergence to benchmark planewave results and results from the literature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Evaluation of a new commercial Monte Carlo dose calculation algorithm for electron beams

Purpose: In this report the authors present the validation of a Monte Carlo dose calculation algorithm (XiO EMC from Elekta Software) for electron beams. Methods: Calculated and measured dose distributions were compared for homogeneous water phantoms and for a 3D heterogeneous phantom meant to approximate the geometry of a trachea and spine. Comparisons of measurements and calculated data were performed using 2D and 3D gamma index dose comparison metrics. Results: Measured outputs agree with calculated values within estimated uncertainties for standard and extended SSDs for open applicators, and for cutouts, with the exception of the 17 MeV electron beam at extended SSD for cutout sizes smaller than 5 × 5 cm{sup 2}. Good agreement was obtained between calculated and experimental depth dose curves and dose profiles (minimum number of measurements that pass a 2%/2 mm agreement 2D gamma index criteria for any applicator or energy was 97%). Dose calculations in a heterogeneous phantom agree with radiochromic film measurements (>98% of pixels pass a 3 dimensional 3%/2 mm γ-criteria) provided that the steep dose gradient in the depth direction is considered. Conclusions: Clinically acceptable agreement (at the 2%/2 mm level) between the measurements and calculated data for measurements in water are obtained for this dose calculation algorithm. Radiochromic film is a useful tool to evaluate the accuracy of electron MC treatment planning systems in heterogeneous media.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2014

New Light-Harvesting Materials Using Accurate and Efficient Bandgap Calculations

Electronic bandgap calculations are presented for 2400 experimentally known materials from the Materials Project database and the bandgaps, obtained with different types of functionals within density functional theory and (partial) self-consistent GW approximation, are compared for 20 randomly chosen compounds forming an unconventional set of ternary and quaternary materials. It is shown that the computationally cheap GLLB-SC potential gives results in good agreement (around 15%) with the more advanced and demanding eigenvalue-self-consistent GW. This allows for a high-throughput screening of materials for different applications where the bandgaps are used as descriptors for the efficiency of a photoelectrochemical device. Here, new light harvesting materials are proposed to be used in a one-photon photoelectrochemical device for water splitting by combining the estimation of the bandgaps with the stability analysis using Pourbaix diagrams and with the evaluation of the position of the band edges. Using this methodology, 25 candidate materials are obtained and 5 of them appear to have a realistic possibility of being used as photocatalyst in a one-photon water splitting device.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Steady and unsteady calculations on thermal striping phenomena in triple-parallel jet

The phenomenon of thermal striping is encountered in liquid metal cooled fast reactors (LMFR), in which temperature fluctuation due to convective mixing between hot and cold fluids can lead to a possibility of crack initiation and propagation in the structure due to high cycle thermal fatigue. Using sodium experiments of parallel triple jets configuration performed by Japan Atomic Energy Agency (JAEA) as benchmark, numerical simulations were carried out to evaluate the temperature fluctuation characteristics in fluid and the transfer characteristics of temperature fluctuation from fluid to structure, which is important to assess the potential thermal fatigue damage. In this study, both steady (RANS) and unsteady (URANS, LES) methods were applied to predict the temperature fluctuations of thermal striping. The parametric studies on the effects of mesh density and boundary conditions on the accuracy of the overall solutions were also conducted. The velocity, temperature and temperature fluctuation intensity distribution were compared with the experimental data. As expected, steady calculation has limited success in predicting the thermal–hydraulic characteristics of the thermal striping, highlighting the limitations of the RANS approach in unsteady heat transfer simulations. The unsteady results exhibited reasonably good agreement with experimental results for temperature fluctuation intensity, as well as the average temperature and velocity components at the measurement locations.

Yu, Y. Q.↗