Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “functional approximation”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2012

First-principles calculations of the electronic structure, phase transition and properties of ZrSiO4 polymorphs

First-principles periodic density functional theory (DFT) calculations have been performed to understand the electronic structure, chemical bonding, phase transition, and physical properties of the mineral zircon (in the chemical composition of ZrSiO4) and its high pressure phase reidite. Temperature effect on phase transition and thermal–mechanical properties such as heat capacity and bulk modulus have been studied by combining the equation of states obtained from DFT calculations with the quasi-harmonic Debye model to take into account the entropy contribution to free energy. Local density approximation (LDA) and generalized gradient approximation (GGA) DFT functionals have been systematically compared in predicting the structure and property of this material. It is found that the LDA functional provides a better description of the equilibrium structure and bulk modulus, while GGA predicts a transition pressure closer to experimental values. Both functionals correctly predict the relative stability of the two phases, with GGA giving slightly larger energy differences. The calculated band structures show that both zircon and reidite have indirect bandgaps and the reidite phase has a narrower bandgap than the zircon phase. The electronic density of states and atomic charges analyses show that bonding in the high-pressure reidite phase has a stronger covalent character.

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

Adsorption of Water Monomer and Clusters on Platinum(111) Terrace and Related Steps and Kinks II. Surface Diffusion

Surface diffusion of water monomer, dimer, and trimer on the (111) terrace, (221) and (322) stepped, and (763) and (854) kinked surfaces of platinum was studied by density functional theory using the PW91 approximation to the energy functional. Monomer diffusion on the terrace is facile, with an activation barrier of 0.20 eV, while dimer and trimer diffusions are restricted due to their high activation barriers of 0.43 and 0.48 eV, respectively. During monomer diffusion on the terrace the O–Pt distance increases by 0.54 Å, about 23% of the initial distance of 2.34 Å. The calculated rate of monomer diffusion hops is in good agreement with the onset temperature of diffusion measurements of Daschbach et al., J. Chem. Phys., 120 (2004) 1516. Alternative monomer diffusion pathways, in which the molecule rolls or flips, were also found. These pathways have diffusion barriers of 0.22 eV. During dimer diffusion on the terrace, the donor molecule rises 0.4 Å at the saddle point, while the acceptor rises by only 0.03 Å. Monomer diffusion up to steps and kinks, with activation barriers of 0.11–0.13 eV, facilitate chain formation on top of step edges. The energy landscape of monomer diffusion from terrace to step to kink sites is downhill with a maximum activation barrier of 0.26 eV. A model for water adsorption is presented inwhichmonomer diffusion leads to concurrent formation of terrace clusters and population of steps/kinks, the latter consistent with the STMmeasurements ofMorgenstern et al., Phys. Rev. Lett., 77 (1996) 703.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2022

On the mechanism of intermolecular nitrogen-atom transfer from a lattice-isolated diruthenium nitride intermediate

Catalyst confinement within microporous media provides the opportunity to site isolate reactive intermediates, enforce intermolecular functionalization chemistry by co-localizing reactive intermediates and substrates in molecular-scale interstices, and harness non-covalent host–guest interactions to achieve selectivities that are complementary to those accessible in solution. As part of an ongoing program to develop synthetically useful nitrogen-atom transfer (NAT) catalysts, we have demonstrated intermolecular benzylic amination of toluene at a Ru 2 nitride intermediate confined within the interstices of a Ru2-based metal–organic framework (MOF), Ru 3 (btc) 2 X 3 (btc = 1,3,5-benzenetricarboxylate, i.e., Ru-HKUST-1 for X = Cl). Nitride confinement within the extended MOF lattice enabled intermolecular C–H functionalization of benzylic C–H bonds in preference to nitride dimerization, which was encountered with soluble molecular analogues. Detailed study of the kinetic isotope effects (KIEs, i.e., k H /k D ) of C–H amination, assayed both as intramolecular effects using partially labeled toluene and as intermolecular effects using a mixture of per-labeled and unlabeled toluene, provided evidence for restricted substrate mobility on the time scale of interstitial NAT. Analysis of these KIEs as a function of material mesoporosity provided approximate experimental values for functionalization in the absence of mass transport barriers. Here, we disclose a combined experimental and computational investigation of the mechanism of NAT from a Ru 2 nitride to the C–H bond of toluene. Computed kinetic isotope effects for a H-atom abstraction (HAA)/radical rebound (RR) mechanism are in good agreement with experimental data obtained for C–H amination at the rapid diffusion limit. These results provide the first detailed analysis of the mechanism of intermolecular NAT to a C–H bond, bolster the use of KIEs as a probe of confinement effects on NAT within MOF lattices, and provide mechanistic insights unavailable by experiment because rate-determining mass transport obscured the underlying chemical kinetics.

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

Accelerating self-consistent field iterations in Kohn-Sham density functional theory using a low-rank approximation of the dielectric matrix

We present an efficient preconditioning technique for accelerating the fixed-point iteration in real-space Kohn-Sham density functional theory (DFT) calculations. The preconditioner uses a low-rank approximation of the dielectric matrix (LRDM) based on Gâteaux derivatives of the residual of fixed-point iteration along appropriately chosen direction functions. We develop a computationally efficient method to evaluate these Gâteaux derivatives in conjunction with the Chebyshev filtered subspace iteration procedure, an approach widely used in large-scale Kohn-Sham DFT calculations. Further, we propose a variant of LRDM preconditioner based on adaptive accumulation of low-rank approximations from previous self-consistent field iterations, and also extend the LRDM preconditioner to spin-polarized Kohn-Sham DFT calculations. We demonstrate the robustness and efficiency of the LRDM preconditioner against other widely used preconditioners on a range of benchmark systems with sizes ranging from ~100 to 1100 atoms (~500–20,000 electrons). The benchmark systems include various combinations of metal-insulating-semiconducting heterogeneous material systems, nanoparticles with localized d orbitals near the Fermi energy, nanofilm with metal dopants, and magnetic systems. In all benchmark systems, the LRDM preconditioner converges robustly within 20–30 iterations. In contrast, other widely used preconditioners show slow convergence in many cases, as well as divergence of the fixed-point iteration in some cases. Lastly, we demonstrate the computational efficiency afforded by the LRDM method, with up to 3.4-fold reduction in computational cost for the total ground-state calculation compared to other preconditioners.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Linear-scaling time-dependent density-functional theory beyond the Tamm-Dancoff approximation: Obtaining efficiency and accuracy with in situ optimised local orbitals

We present a solution of the full time-dependent density-functional theory (TDDFT) eigenvalue equation in the linear response formalism exhibiting a linear-scaling computational complexity with system size, without relying on the simplifying Tamm-Dancoff approximation (TDA). The implementation relies on representing the occupied and unoccupied subspaces with two different sets of in situ optimised localised functions, yielding a very compact and efficient representation of the transition density matrix of the excitation with the accuracy associated with a systematic basis set. The TDDFT eigenvalue equation is solved using a preconditioned conjugate gradient algorithm that is very memory-efficient. The algorithm is validated on a small test molecule and a good agreement with results obtained from standard quantum chemistry packages is found, with the preconditioner yielding a significant improvement in convergence rates. The method developed in this work is then used to reproduce experimental results of the absorption spectrum of bacteriochlorophyll in an organic solvent, where it is demonstrated that the TDA fails to reproduce the main features of the low energy spectrum, while the full TDDFT equation yields results in good qualitative agreement with experimental data. Furthermore, the need for explicitly including parts of the solvent into the TDDFT calculations is highlighted, making the treatment of large system sizes necessary that are well within reach of the capabilities of the algorithm introduced here. Finally, the linear-scaling properties of the algorithm are demonstrated by computing the lowest excitation energy of bacteriochlorophyll in solution. The largest systems considered in this work are of the same order of magnitude as a variety of widely studied pigment-protein complexes, opening up the possibility of studying their properties without having to resort to any semiclassical approximations to parts of the protein environment.

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

Hartree approximation to the one loop quantum gravitational correction to the graviton mode function on de Sitter

We use the Hartree approximation to the Einstein equation on de Sitter background to solve for the one loop correction to the graviton mode function. This should give a reasonable approximation to how the ensemble of inflationary gravitons affects a single external graviton. At late times we find that the one loop correction to the plane wave mode function u(η,k) goes like GH{sup 2}ln (a)/a{sup 2}, where a is the inflationary scale factor. One consequence is that the one loop corrections to the ''electric'' components of the linearized Weyl tensor grow compared to the tree order result.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

ORNL_AISD-Ex: Quantum chemical prediction of UV/Vis absorption spectra for over 10 million organic molecules

We performed calculations of electronic excitation energies and associated oscillator strengths based on the time-dependent density-functional tight-binding (TD-DFTB) method [1]. The SMILES (Simplified molecular-input line-entry system) strings of the molecules from the AISD HOMO-LUMO database [2] were converted to a 3D atomistic structure and stored in a PDB file after preliminary geometry optimization using the Merck Molecular Force Field (MMFF94) in RDKit [3,4]. The primary information stored in the PDB file archive consists of Cartesian coordinates for each atom of the molecule in their 3D location in space, along with summary information about the structure, sequence, and experiment. We then performed molecular geometry optimization using the density-functional tight-binding (DFTB) method [5] in the electronic ground state, followed by single-point excited states calculations, as described below. We note that, since RDKit employs a random choice for the generation of molecular conformers, the molecular geometries obtained in this dataset could be different from the ones that were generated when the AISD HOMO-LUMO dataset was generated. The computed excitation energies and associated oscillator strengths can be converted to predict UV/Vis absorption spectra, where excitation energies correspond to absorption peak positions, and oscillator strengths are a good measure of the probability of absorption of visible or UV light in transitions between electronic ground and excited states. The conversion of SMILES strings to 3D Cartesian coordinates of fully DFTB-optimized molecules was successful for 10,502,904 out of 10,502,917 molecules. For these molecules, both geometry optimizations and excited states calculations were successful. The DFTB calculations did not complete for 13 molecules of the original AISD HOMO-LUMO dataset. We still provide information about the geometry of these molecules. The molecules are diverse for chemical compositions (which span 5 non-hydrogen elements: oxygen, carbon, nitrogen, fluorine, sulfur) and molecular size (the smallest molecule contains 5 non-hydrogen atoms, and the largest molecule contains 71 non-hydrogen atoms). The DFTB method [5] is an approximation to density functional theory (DFT), utilizing a minimal basis set in conjunction with a two-center approximation to the electronic Hamiltonian and overlap matrix elements. The DFTB total energy is the sum of an electronic and a repulsive energy contribution, and their calculation requires optimized electronic parameters and diatomic repulsive potential energy functions. All DFTB calculations were performed using the DFTB+ code [6] (version 21.2) and the wrapper for DFTB+ in the Atomic Simulation Environment (ASE) (version 3.22.1) [7], which performed an internal conversion of Cartesian coordinates from PDB to the .gen file format. For the geometry optimizations on the electronic ground state potential energy surface of the molecules, we have chosen the third-order DFTB (DFTB3) method [5c] and employed the matching 3ob set of electronic parameters and repulsive potentials [8]. The empirical γ-damping for hydrogen bond correction, and Grimme's D3 empirical dispersion correction with Becke-Johnson damping (D3(BJ)) [9] dispersion correction was included to improve the description of non-covalent interactions. For excited states single-point energy calculations, we employed the TD-DFTB method in conjunction with the DFTB2 method [5b] and the matching mio [5b,10] and halorg [11] parameter sets. We opted to request the simultaneous calculation of 50 excited states for singlet transition to investigate sufficient number of excited states, based on linear response theory using the Casida equation [Ref: T. A. Niehaus, S. Suhai, F. Della Sala, P Lugli, M. Elstner, G. Seifert, and Th. Frauenheim. Tight-binding approach to time-dependent density-functional response theory. Phys. Rev. B, 63:085108, 2001] and the ARPACK diagonalizer [R. B. Lehoucq, D. C. Sorensen, and C. Yang. Arpack users guide: Solution of large-scale eigenvalue problems by implicitly restarted arnoldi methods, 1997. 46, 51]. The dataset contains 1001 tar.gz files. Tar files are named as “ornl_aisd_ex_1.tar.gz†through “ornl_aisd_ex_1000.tar.gzâ€. Additionally, the 13 failed molecules are in “ornl_aisd_ex_unprocessed.tar.gzâ€. Except for the tar files listed below, each tar file contains 10,500 molecules. Tar files numbered 34, 121, 128, 352, 360, 429, 495, 509, 518, 627, 676, 668, and 862 contain 10,499 molecules each. The last tar file numbered 1000 contains 13,417 molecules. The total size of the uncompressed dataset is over 283 Gigabytes. The code for calculating the electronic excitation energies and statistical analysis of the dataset is provided at the following GitLab repository: https://github.com/ORNL/Analysis-of-Large-Scale-Molecular-Datasets-with-Python Calculating the UV spectrum of a molecule requires performing 3 main operations: 1. Converting the smiles string representation of a molecule into a geometric structure where each atom is assigned XYZ coordinates. The geometric structure is written to the file smiles.pdb. 2. Using smiles.pdb to compute the relaxed geometry of the molecule, which corresponds with the position of the atoms at the position of equilibrium at the ground state. This generates the files band.out, detailed.out, and geo_end.gen. 3. Using geo_end.gen to calculate the UV spectrum of the molecule which is written into the file EXC.DAT. Every molecule in the dataset has its own directory. The files contained in each molecule directory are as follows: 1. geo_end.gen 2. detailed.out 3. band.out 4. EXC.DAT 5. smiles.pdb REFERENCES [1] Niehaus, T. A.; Suhai, S.; Della Salla, F.; Lugli, P.; Elstner, M.; Seifert, G.; Frauenheim, Th. Tight-binding approach to time-dependent density-functional response theory. Phys. Rev. B, 2001, 63, 085108/1-9. [2] Blanchard, A.; Gounley, J.; Metha, K.; Yoo, P.; Irle, S. AISD HOMO-LUMO. DOI: 10.13139/ORNLNCCS/1869409 [3] RDKit: Cheminformatics and Machine Learning Software. 2013, [http://www.rdkit.org] [4] Tosco, P.; Stiefl, N. and Landrum, G. Bringing the MMFF force field to the RDKit: implementation and validation. J Cheminform. 2014, 6, 1–4. [5] a) Porezag, D.; Frauenheim, T.; Kohler, T.; Seifert, G.; Kaschner, Construction of tight-binding-like potentials on the basis of density-functional theory: Application to carbon, R. Phys. Rev. B 1995, 51, 12947-12957; b) Elstner, M.; Porezag, D.; Jungnickel, G.; Elsner, J.; Haugk, M.; Frauenheim, Th.; Suhai, S.; Seifert, G.; Phys. Rev. B 1998, 58, 7260-7268; c) Gaus, M.; Cui, Q.; Elstner, M. DFTB3: Extension of the Self-Consistent-Charge Density-Functional Tight-Binding Method (SCC-DFTB), J. Chem. Theory Comput. 2011, 7, 931-948; d) Cui, Q.; Elstner, M. Density functional tight binding: values of semi-empirical methods in an ab initio era, Phys. Chem. Chem. Phys. 2014, 16, 14368-14377. [6] Hourahine, B. et al. DFTB+, a software package for efficient approximate density functional theory based atomistic simulations, J. Chem. Phys. 2020, 152, 124101/1-19. [7] Larsen, A. H. et al. The atomic simulation environment—a Python library for working with atoms. J. Phys.: Cond. Matter 2017, 29, 273002. [8] Kubillus, M.; Kubar, T.; Gaus, M.; Rezac, J.; Elstner, M. Parameterization of the DFTB3 Method for Br, Ca, Cl, F, I, K, and Na in Organic and Biological Systems, J. Chem. Theory Comput. 2015, 11, 332-342. [9] Brandenburg, J. G.; Grimme, S. Accurate Modeling of Organic Molecular Crystals by Dispersion-Corrected Density Functional Tight Binding (DFTB), J. Phys. Chem. Lett. 2014, 5, 1785−1789. [10] a) Niehaus, T. A.; Elstner, M.; Frauenheim, Th.; Suhai, S. Application of an approximate density-functional method to sulfur containing compounds. J. Mol. Struct.: THEOCHEM 2001, 541, 185-94; b) Elstner, M.; Hobza, P.; Frauenheim, Th.; Suhai, S.; Kaxiras, E. Hydrogen bonding and stacking interactions of nucleic acid base pairs: A density-functional-theory based treatment. J. Chem. Phys. 2001, 114, 5149-55. [11] Kubar, T.; Bodrog, Z.; Gaus, M.; Köhler, C.; Aradi, B.; Frauenheim, Th.; Elstner, M. Parametrization of the SCC-DFTB Method for Halogens. J. Chem. Theory Comput. 2013, 9, 2939-49.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Single particle spectral function in iron pnictide superconductors within two band model

The newly discovered iron based superconductors have shown a great deal of interest in the condensed matter physics community around the world to explore the experimental and theoretical aspects of electronic properties of these materials to underline the technological potential. These materials have a layered tetragonal crystal structure and generally show small anisotropy. Based on ARPES data to understand the band structure of these systems, we present a theoretical tight binding model and numerical computation of the single particle spectral function within two orbital per site for iron pnictide superconductors. The two band tight binding model Hamiltonian containing various orbitals hopping energies, intra- and inter- band electronic correlations and Hund’s coupling energy in Fe 3d orbitals has been used. The expressions of single particle spectral function have been obtained by employing the Green’s function equation of motion approach within BCS-mean-field approximation. The spectral function is numerically computed at different k-points of Brillouin zone in extended s-wave paring symmetry. It is noticed that the behavior of electronic states is different at different k-points of Brillouin zone and highly influenced by onsite Coulomb interactions. Further, It is predicted that the presence of onsite coulomb correlation suppress the spectral weight close to Fermi level in iron pnictide systems. On the basis of numerical computation we have compared our theoretical results with recent angle resolved photoemission spectroscopic ARPES data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Composite-particle decay widths by the generator coordinate method

We study the feasibility of applying the Generator Coordinate Method (GCM) of self-consistent mean-field theory to calculate decay widths of composite particles to composite-particle final states. The main question is how well the GCM can approximate continuum wave functions in the decay channels. The analysis is straightforward under the assumption that the GCM wave functions are separable into internal and Gaussian center-of-mass wave functions. Two methods are examined for calculating decay widths. In one method, the density of final states is computed entirely in the GCM framework. In the other method, it is determined by matching the GCM wave function to an asymptotic scattering wave function. Both methods are applied to a numerical example and are found to agree within their determined uncertainties.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Random Phase Approximation Correlation Energy Using Real-Space Density Functional Perturbation Theory

We present a real-space method for computing the random phase approximation (RPA) correlation energy within Kohn–Sham density functional theory, leveraging the low-rank nature of the frequency-dependent density response operator. In particular, we employ a cubic-scaling formalism based on density functional perturbation theory that circumvents the calculation of the response function matrix, instead relying on the ability to compute its product with a vector through the solution of the associated Sternheimer linear systems. We develop a large-scale parallel implementation of this formalism using the subspace iteration method in conjunction with the spectral quadrature method while employing the Kronecker product-based method for the application of the Coulomb operator and the conjugate orthogonal conjugate gradient method for the solution of the linear systems. We demonstrate convergence with respect to key parameters and verify the method’s accuracy by comparing with plane-wave results. We show that the framework achieves good strong scaling to many thousands of processors, reducing the time to solution for a lithium hydride system with 128 electrons to around 150 s on 4608 processors.

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

The Predictive Power of Exact Constraints and Appropriate Norms in Density Functional Theory

Ground-state Kohn-Sham density functional theory provides, in principle, the exact ground-state energy and electronic spin densities of real interacting electrons in a static external potential. In practice, the exact density functional for the exchange-correlation (xc) energy must be approximated in a computationally efficient way. About 20 mathematical properties of the exact xc functional are known. In this work, we review and discuss these known constraints on the xc energy and hole. By analyzing a sequence of increasingly sophisticated density functional approximations (DFAs), we argue that (a) the satisfaction of more exact constraints and appropriate norms makes a functional more predictive over the immense space of many-electron systems and (b) fitting to bonded systems yields an interpolative DFA that may not extrapolate well to systems unlike those in the fitting set. We discuss both how the class of well-described systems has grown along with constraint satisfaction and the possibilities for future functional development.

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

AN EFFICIENT APPROXIMATION TO THE LIKELIHOOD FOR GRAVITATIONAL WAVE STOCHASTIC BACKGROUND DETECTION USING PULSAR TIMING DATA

Direct detection of gravitational waves by pulsar timing arrays will become feasible over the next few years. In the low frequency regime (10{sup -7} Hz-10{sup -9} Hz), we expect that a superposition of gravitational waves from many sources will manifest itself as an isotropic stochastic gravitational wave background. Currently, a number of techniques exist to detect such a signal; however, many detection methods are computationally challenging. Here we introduce an approximation to the full likelihood function for a pulsar timing array that results in computational savings proportional to the square of the number of pulsars in the array. Through a series of simulations we show that the approximate likelihood function reproduces results obtained from the full likelihood function. We further show, both analytically and through simulations, that, on average, this approximate likelihood function gives unbiased parameter estimates for astrophysically realistic stochastic background amplitudes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

On Coarse Grid Correction Methods in Krylov Subspaces

Two approaches using coarse grid correction in the course of a certain Krylov iterative process are presented. The aim of the correction is to accelerate the iterations. These approaches are based on an approximation of the function sought for by simple basis functions having finite supports. Additional acceleration can be achieved if the iterative process is restarted and the approximate solution is refined. In this case, the resulting process turns out to be a two-level preconditioned method. The influence of different parameters of the iterative process on its convergence is demonstrated by numerical results.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Role of electron correlation on the adenine dimer interaction for non-equilibrium geometries: a benchmark Quantum Monte Carlo study

The accurate description of non-covalent interactions is critical for understanding the structure, dynamics, and eventual function of biomolecules. The adenine dimer serves as a benchmark system for computational methods due to its role in nucleic acid structures and its rich conformational landscape. In this study, we employ benchmark diffusion quantum Monte Carlo (DMC) methods to investigate the relative energies and role of electron correlation on a set of adenine dimer conformations generated via a search of the potential energy landscape using the global optimizer algorithm. Relative DMC energies are compared against a wide range of density functional theory (DFT) approximation results. We find that although most of the DFT functionals perform well for low-energy structures, their accuracy varies significantly for higher-energy conformations, including stacked and T-shaped structures. A large fraction of the variation is due to the treatment of the van der Waals interaction. BLYP, B3LYP, and PBE0 significantly improve with added D4 dispersion, while the recent r2SCAN-D4 and ωB97M-V functionals show the least scatter and closest agreement with the DMC. These findings highlight the delicate nature of these interactions in biomolecular systems and provide guidance for simulations of their structure and dynamics and for the development of machine learned interatomic potentials.

Washburn, Laurel [ORNL] (ORCID:0000000324179335)↗
DOE OSTI2012

Microscopic nuclear mass table with high-performance computing

We overview the methodology behind the large-scale mass table calculations based on the nuclear density functional theory (DFT) with Skyrme energy density functionals (EDFs). The calculations employing massively parallel computers are done in a large configuration space of an axially deformed harmonic oscillator. Nuclear mass tables tabulating global nuclear properties such as binding energies, radii, shape deformations, and pairing gaps are obtained for several Skyrme EDFs augmented by a mixed-type pairing functional, using an approximate particle number projection before variation. Specialty visualization tools have been developed to analyze the results. As an illustrative example of our current capabilities, we show some results for a wide range of even-even nuclei with Z ≤ 120 and N ≤ 300.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

Cohesive energy and structural parameters of binary oxides of groups IIA and IIIB from diffusion quantum Monte Carlo

We have applied the diffusion quantum Monte Carlo (DMC) method to calculate the cohesive energy and the structural parameters of the binary oxides CaO, SrO, BaO, Sc{sub 2}O{sub 3}, Y{sub 2}O{sub 3}, and La{sub 2}O{sub 3}. The aim of our calculations is to systematically quantify the accuracy of the DMC method to study this type of metal oxides. The DMC results were compared with local, semi-local, and hybrid Density Functional Theory (DFT) approximations as well as with experimental measurements. The DMC method yields cohesive energies for these oxides with a mean absolute deviation from experimental measurements of 0.18(2) eV, while with local, semi-local, and hybrid DFT approximations, the deviation is 3.06, 0.94, and 1.23 eV, respectively. For lattice constants, the mean absolute deviations in DMC, local, semi-local, and hybrid DFT approximations are 0.017(1), 0.07, 0.05, and 0.04 Å, respectively. DMC is a highly accurate method, outperforming the DFT approximations in describing the cohesive energies and structural parameters of these binary oxides.

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

Modified Steklov Functions and Numerical Differentiation Formulas

We consider an approximation method based on Steklov functions of the first and second order. We obtain estimates for the norms in the space C of continuous periodic functions and clarify how they connect with numerical differentiation formulas.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Analytical Approximations for Beamstrahlung at Very High Energy Electron-Positron Colliders

Among the many effects that occur in beam-beam electron-positron collisions at TeV energies, emission of hard synchrotron radiation, or beamstrahlung, has special importance. Beamstrahlung determines the energy spectrum of the most energetic electrons, positrons, and photons and supplies the initial condition for the calculation of all other QED processes. In this paper, we show that the description of beamstrahlung simplifies in the limit of large quantum parameter $Υ$, which is realized in 10 TeV collider designs. The beamstrahlung spectra for electrons and photons are given in terms of universal functions. We supply approximations to these functions that will be useful for more general studies of the beam-beam interaction at very high energies.

Accelerator Physics (physics.acc-ph)↗