The Quality of the Embedding Potential Is Decisive for Minimal Quantum Region Size in Embedding Calculations: The Case of the Green Fluorescent Protein
Explore the source record for details and available documents.
Nåbo, Lina J.↗
SEARCH · Search NASA
Results for “Calculation”
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.
First principles calculations in support of Pu aging: calculating the effects of lattice imperfections on thermodynamics
Plutonium (Pu) has, in theory, well defined crystal structures: its atoms are arranged in regular spatial patterns. But Pu is radioactive, and as its nuclei decay those regular spatial patterns are interrupted. The interruptions are lattice imperfections, which are known to affect how materials respond to their environment. To adequately model Pu, we need to know which lattice imperfections are present, how they interact with each other, and how they affect the material’s response to its environment. The work presented here aims to use density functional theory (DFT) calculations to begin to answer the latter, in particular, how individual lattice imperfections affect measurable effects including thermal expansion (the change in volume in response to a change in temperature), heat capacity (the amount of thermal energy needed to change a material’s temperature), and elastic moduli (a material’s resistance to applied stresses). Pu poses many computational challenges. The $\textit{f}$ electrons require special attention. Of all the elements, Pu has the largest number of electrons that must be included in the calculations. Thermal effects require calculating the phonons (the lattice vibrations), which for systems containing lattice imperfections demand large, complex computational cells - but computational resources limit the size and complexity. With careful restructuring of how the calculations are performed, all these challenges have been met to enable calculations that provide insight into how lattice imperfections affect Pu’s response to its environment. Reported here are the computational challenges and the advances developed to meet them, along with the first results showing the strong effect that one prototype lattice imperfection (an interstitial Pu atom in a delta-phase Pu lattice) has on Pu’s response to its environment.
36 MATERIALS SCIENCE↗
ICE Calculator 2.0: Final Report for Phase 1 of the National Initiative to Update the Interruption Cost Estimate (ICE) Calculator
In 2021, Berkeley Lab and Resource Innovations, Inc. launched the “ICE 2.0 Initiative” – a national study to refresh the underlying data and enhance the functionality of the ICE Calculator. The Initiative involves Berkeley Lab contracting with sponsoring utilities to administer identical, updated and comprehensive interruption cost surveys to statistically representative samples of each utility’s customers. Berkeley Lab and Resource Innovations then pool the survey results across the utilities and use them to update the analytical engines that drive the ICE Calculator. The ICE 2.0 Initiative is being conducted in phases. Each phase involves the administration of interruption cost surveys to the customers of sponsoring utilities, followed by an update to the ICE Calculator based on analysis of the pooled survey results. This report describes the activities and findings from Phase 1 of the ICE 2.0 Initiative. Phase 1 was sponsored by eight utilities: American Electric Power, Commonwealth Edison, Dominion Energy, Duke Energy, DTE Electric, Exelon, National Grid, and Puget Sound Energy. Phase 1 involved 11 customer interruption cost survey activities representing a total of 24 electricity distribution service territories, 23 of them located in the Eastern and Midwestern regions of the U.S. and one located in the Pacific Northwest. ICE 2.0 vs. 1.0 Comparison This memorandum compares customer power interruption costs estimated using the recently updated Interruption Cost Estimate (ICE) Calculator (“ICE 2.0”) to the original ICE Calculator (“ICE 1.0”). ICE 1.0 was developed in 2009 based on 15 independent power interruption cost surveys conducted by 10 electric utilities between 1989 and 2012. ICE 2.0 was developed in 2025 through a national initiative based on a consistent set of power interruption cost surveys and 11 surveying efforts conducted across 24 electric utility service territories between 2022 and 2024.
32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
Lithium formate ion clusters formation during electrospray ionization: Evidence of magic number clusters by mass spectrometry and ab initio calculations
Small cationic and anionic clusters of lithium formate were generated by electrospray ionization and their fragmentations were studied by tandem mass spectrometry (collision-induced dissociation with N{sub 2}). Singly as well as multiply charged clusters were formed in both positive and negative ion modes with the general formulae, (HCOOLi){sub n}Li{sup +}, (HCOOLi){sub n}Li{sub m}{sup m+}, (HCOOLi){sub n}HCOO{sup −}, and (HCOOLi){sub n}(HCOO){sub m}{sup m−}. Several magic number cluster (MNC) ions were observed in both the positive and negative ion modes although more predominant in the positive ion mode with (HCOOLi){sub 3}Li{sup +} being the most abundant and stable cluster ion. Fragmentations of singly charged positive clusters proceed first by the loss of a dimer unit ((HCOOLi){sub 2}) followed by the loss of monomer units (HCOOLi) although the former remains the dominant dissociation process. In the case of positive cluster ions, all fragmentations lead to the magic cluster (HCOOLi){sub 3}Li{sup +} as the most abundant fragment ion at higher collision energies which then fragments further to dimer and monomer ions at lower abundances. In the negative ion mode, however, singly charged clusters dissociated via sequential loss of monomer units. Multiply charged clusters in both positive and negative ion modes dissociated mainly via Coulomb repulsion. Quantum chemical calculations performed for smaller cluster ions showed that the trimer ion has a closed ring structure similar to the phenalenylium structure with three closed rings connected to the central lithium ion. Further additions of monomer units result in similar symmetric structures for hexamer and nonamer cluster ions. Thermochemical calculations show that trimer cluster ion is relatively more stable than neighboring cluster ions, supporting the experimental observation of a magic number cluster with enhanced stability.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Multireference configuration interaction calculations of the first six ionization potentials of the uranium atom
The first 6 ionization potentials (IPs) of the uranium atom have been calculated using multireference configuration interaction (MRCI+Q) with extrapolations to the complete basis set limit using new all-electron correlation consistent basis sets. The latter was carried out with the third-order Douglas-Kroll-Hess Hamiltonian. Correlation down through the 5s5p5d electrons has been taken into account, as well as contributions to the IPs due to the Lamb shift. Spin-orbit coupling contributions calculated at the 4-component Kramers restricted configuration interaction level, as well as the Gaunt term computed at the Dirac-Hartree-Fock level, were added to the best scalar relativistic results. The final ionization potentials are expected to be accurate to at least 5 kcal/mol (0.2 eV) and thus more reliable than the current experimental values of IP{sub 3} through IP{sub 6}.
74 ATOMIC AND MOLECULAR PHYSICS↗
Self-consistent calculation of the Sommerfeld enhancement
A calculation of the Sommerfeld enhancement is presented and applied to the problem of s-wave non-relativistic dark matter annihilation. The difference from previous computations in the literature is that the effect of the underlying short-range scattering process is consistently included together with the long-range force in the effective QM Schrödinger problem. Our procedure satisfies partial-wave unitarity where previous calculations fail. We provide analytic results for some potentials of phenomenological relevance.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
An efficient approximation of the supercell approach to the calculation of the full phonon spectrum
Not Available
Duong, Thien C. (ORCID:000000018694814X)↗
GPU implementation of the linear scaling three dimensional fragment method for large scale electronic structure calculations
LS3DF, namely linear scaling three-dimensional fragment method, is an efficient linear scaling ab initio total energy electronic structure calculation code based on a divide-and-conquer strategy. In this paper, we present our GPU implementation of the LS3DF code. Our test results show that the GPU code can calculate systems with about ten thousand atoms fully self-consistently in the order of 10 min using thousands of computing nodes. This makes the electronic structure calculations of 10,000-atom nanosystems routine work. This speed is 4.5–6 times faster than the CPU calculations using the same number of nodes on the Titan machine in the Oak Ridge leadership computing facility (OLCF). Such speedup is achieved by (a) carefully re-designing of the computationally heavy kernels; (b) redesign of the communication pattern for heterogeneous supercomputers.
Computer Science↗
Parton distributions and lattice-QCD calculations: Toward 3D structure
The strong force which binds hadrons is described by the theory of quantum chromodynamics (QCD). Determining the character and manifestations of QCD is one of the most important and challenging outstanding issues necessary for a comprehensive understanding of the structure of hadrons. Within the context of the QCD parton picture, the parton distribution functions (PDFs) have been remarkably successful in describing a wide variety of processes. However, these PDFs have generally been confined to the description of collinear partons within the hadron. New experiments and facilities provide the opportunity to additionally explore the transverse structure of hadrons which is described by generalized parton distributions (GPDs) and transverse-momentum-dependent parton distribution functions (TMD PDFs). In our previous report Lin et al. (2018), we compared and contrasted the two main approaches used to determine the collinear PDFs: the first based on perturbative QCD factorization theorems, and the second based on lattice-QCD calculations. In the present report, we provide an update of recent progress on the collinear PDFs, and also expand the scope to encompass the generalized PDFs (GPDs and TMD PDFs). We review the current state of the various calculations, and consider what new data might be available in the near future. We also examine how a shared effort can foster dialog between the PDF and lattice-QCD communities, and yield improvements for these generalized PDFs.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Radiation resistance of oxide dispersion strengthened alloys: Perspectives from in situ observations and rate theory calculations
Here, in situ ion irradiation and rate theory calculations were employed to directly compare the radiation resistance of an oxide dispersion strengthened alloy with that of a conventional ferritic/martensitic alloy. Compared to the rapid buildup of dislocation loops, loop growth, and formation of network dislocations in the conventional ferritic/martensitic alloy, the superior radiation resistance of the oxide dispersion strengthened alloy is manifested by its stable dislocation structure under the same irradiation conditions. Thus, the results are consistent with rate theory calculations, which show that high-density nanoparticles can significantly reduce freely migrating defects and suppress the buildup of clustered defects.
36 MATERIALS SCIENCE↗
A calculator to quantify cover crop effects on soil health and productivity
Not Available
Jian, Jinshi↗
Calculation of Electrochemical Reorganization Energies for Redox Molecules at Self-Assembled Monolayer Modified Electrodes
Electrochemical electron transfer reactions play an important role in energy conversion processes with many technological applications. Electrodes modified by self-assembled monolayers (SAMs) are useful because the double layer effects are reduced. An important quantity for calculating the electron transfer rate constant is the reorganization energy, which is associated with changes in solute geometry and solvent configuration. In this Letter, an approach for calculating the electrochemical solvent reorganization energy for a redox molecule attached to or near a SAM modified electrode is presented. This integral equations formalism polarizable continuum model (IEF-PCM) approach accounts for the detailed electronic structure of the molecule, as well as the contributions from the electrode, SAM, and electronic and inertial solvent responses. The calculated total reorganization energies are in good agreement with experimental data for a series of metal complex in aqueous solution. This approach will be useful for calculating electron transfer rate constants for molecular electrocatalysts. This work was supported as part of the Center for Molecular Electrocatalysis, an Energy Frontier Research Center funded by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences.
Ghosh, Soumya↗
New XDM-corrected potential energy surfaces for Ar–NO(X{sup 2}Π): A comparison with CCSD(T) calculations and experiments
We report new potential energy surfaces for the ground state Ar–NO(X{sup 2}Π) van der Waals system calculated using the unrestricted Hartree-Fock (UHF) method with the addition of the Becke-Roussel correlation functional and exchange-hole dipole moment dispersion correction (XDM). We compare UHFBR-XDM surfaces and those previously reported by Alexander from coupled cluster CCSD(T) calculations [J. Chem. Phys. 111, 7426 (1999)]. The bound states of Ar–NO have been investigated with these new UHFBR-XDM surfaces, including relative energy-level spacing, adiabatic bender states and wave functions, and spectroscopic data. These results have been found to be in good agreement with calculations based on the CCSD(T) PESs. These new PESs are used to investigate the inelastic scattering of NO(X) by Ar. Full close-coupling integral cross sections at collision energies of 442 cm{sup −1}, 1774 cm{sup −1} and differential cross sections at collision energy of 530 cm{sup −1} were determined for transitions out of the lowest NO(X) rotational level (j = ω = 1/2,f). These cross sections are in good agreement with those calculated with CCSD(T) and accordingly in good agreement with the most recent initial and final state resolved experimental data. The UHFBR-XDM scheme yields high-quality potential surfaces with computational cost comparable to the Hartree-Fock method and our results may serve as a benchmark for application of this scheme to collisions between larger molecules.
74 ATOMIC AND MOLECULAR PHYSICS↗
Native defects in Tl{sub 6}SI{sub 4}: Density functional calculations
Tl{sub 6}SI{sub 4} is a promising room-temperature semiconductor radiation detection material. Here, we report density functional calculations of native defects and dielectric properties of Tl{sub 6}SI{sub 4}. Formation energies and defect levels of native point defects and defect complexes are calculated. Donor-acceptor defect complexes are shown to be abundant in Tl{sub 6}SI{sub 4}. High resistivity can be obtained by Fermi level pinning by native donor and acceptor defects. Deep donors that are detrimental to electron transport are identified and methods to mitigate such problem are discussed. Furthermore, we show that mixed ionic-covalent character of Tl{sub 6}SI{sub 4} gives rise to enhanced Born effective charges and large static dielectric constant, which provides effective screening of charged defects and impurities.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Quantum mechanical calculations of state-to-state cross sections and rate constants for the F + DCl → Cl + DF reaction
We present accurate state-to-state quantum wave packet calculations of integral cross sections and rate constants for the title reaction. Calculations are carried out on the best available ground 1{sup 2}A′ global adiabatic potential energy surface of Deskevich et al. [J. Chem. Phys. 124, 224303 (2006)]. Converged state-to-state reaction cross sections have been calculated for collision energies up to 0.5 eV and different initial rotational and vibrational excitations, DCl(v = 0, j = 0 − 1; v = 1, j = 0). Also, initial-state resolved rate constants of the title reaction have been calculated in a temperature range of 100-400 K. It is found that the initial rotational excitation of the DCl molecule does not enhance reactivity, in contract to the reaction with the isotopologue HCl in which initial rotational excitation produces an important enhancement. These differences between the isotopologue reactions are analyzed in detail and attributed to the presence of resonances for HCl(v = 0, j), absent in the case of DCl(v = 0, j). For vibrational excited DCl(v = 1, j), however, the reaction cross section increases noticeably, what is also explained by another resonance.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Fragment approach to constrained density functional theory calculations using Daubechies wavelets
In a recent paper, we presented a linear scaling Kohn-Sham density functional theory (DFT) code based on Daubechies wavelets, where a minimal set of localized support functions are optimized in situ and therefore adapted to the chemical properties of the molecular system. Thanks to the systematically controllable accuracy of the underlying basis set, this approach is able to provide an optimal contracted basis for a given system: accuracies for ground state energies and atomic forces are of the same quality as an uncontracted, cubic scaling approach. This basis set offers, by construction, a natural subset where the density matrix of the system can be projected. In this paper, we demonstrate the flexibility of this minimal basis formalism in providing a basis set that can be reused as-is, i.e., without reoptimization, for charge-constrained DFT calculations within a fragment approach. Support functions, represented in the underlying wavelet grid, of the template fragments are roto-translated with high numerical precision to the required positions and used as projectors for the charge weight function. We demonstrate the interest of this approach to express highly precise and efficient calculations for preparing diabatic states and for the computational setup of systems in complex environments.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Towards a systematic assessment of errors in diffusion Monte Carlo calculations of semiconductors: Case study of zinc selenide and zinc oxide
The fixed node diffusion Monte Carlo (DMC) method has attracted interest in recent years as a way to calculate properties of solid materials with high accuracy. However, the framework for the calculation of properties such as total energies, atomization energies, and excited state energies is not yet fully established. Several outstanding questions remain as to the effect of pseudopotentials, the magnitude of the fixed node error, and the size of supercell finite size effects. Here, we consider in detail the semiconductors ZnSe and ZnO and carry out systematic studies to assess the magnitude of the energy differences arising from controlled and uncontrolled approximations in DMC. The former include time step errors and supercell finite size effects for ground and optically excited states, and the latter include pseudopotentials, the pseudopotential localization approximation, and the fixed node approximation. We find that for these compounds, the errors can be controlled to good precision using modern computational resources and that quantum Monte Carlo calculations using Dirac-Fock pseudopotentials can offer good estimates of both cohesive energy and the gap of these systems. We do however observe differences in calculated optical gaps that arise when different pseudopotentials are used.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗