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.
Increasing the applicability of density functional theory. IV. Consequences of ionization-potential improved exchange-correlation potentials
This paper's objective is to create a “consistent” mean-field based Kohn-Sham (KS) density functional theory (DFT) meaning the functional should not only provide good total energy properties, but also the corresponding KS eigenvalues should be accurate approximations to the vertical ionization potentials (VIPs) of the molecule, as the latter condition attests to the viability of the exchange-correlation potential (V{sub XC}). None of the prominently used DFT approaches show these properties: the optimized effective potential V{sub XC} based ab initio dft does. A local, range-separated hybrid potential cam-QTP-00 is introduced as the basis for a “consistent” KS DFT approach. The computed VIPs as the negative of KS eigenvalue have a mean absolute error of 0.8 eV for an extensive set of molecule's electron ionizations, including the core. Barrier heights, equilibrium geometries, and magnetic properties obtained from the potential are in good agreement with experiment. A similar accuracy with less computational efforts can be achieved by using a non-variational global hybrid variant of the QTP-00 approach.
Unsupervised Learning for Equitable DER Control: Preprint
In the context of managing distributed energy resources (DERs) within distribution networks (DNs), this work focuses on the task of developing local controllers. We propose an unsupervised learning framework to train functions that can closely approximate optimal power flow (OPF) solutions. The primary aim is to establish specific conditions under which these learned functions can collectively guide the network towards desired configurations asymptotically, leveraging an incremental control approach. The flexibility of the proposed methodology allows to integrate fairness-driven components into the cost function associated with the OPF problem. This addition seeks to mitigate power curtailment disparities among DERs, thereby promoting equitable power injections across the network. To demonstrate the effectiveness of the proposed approach, power flow simulations are conducted using the IEEE 37-bus feeder. The findings not only showcase the guaranteed system stability but also underscore its improved overall performance.
Electric double layer for a size-asymmetric electrolyte around a spherical colloid
We have studied the structure of a size-asymmetric electrolyte on charged colloids by a density functional perturbation theory. The hard-sphere contribution has been approximated as the direct pair correlation function with the coupling parameter, whereas the electronic contribution has been approximated as the mean-spherical approximation in the bulk phase. The calculated results for the ionic density distributions and mean electrostatic potentials are in very good agreement with the computer simulations over a wide range of colloid sizes and electrolyte concentrations. The present theory provides better structural results than the hypernetted-chain equation based on the mean spherical approximation. We have confirmed that the overcharging appears when the counterions are larger than the coions. The overcharging disappears everywhere when the electrostatic repulsion becomes strong enough, while the charge reversal is observed when the coions are larger than the counterions, and the reversal effect appears for a size-asymmetric electrolyte at high surface charge densities. The charge reversal occurs even for the point of zero charge, mainly due to the depletion force between two ions. The present theory is able to provide interesting insights about the charge reversal and overcharging phenomena occurring at the interface.
Communication: Two-component ring-coupled-cluster computation of the correlation energy in the random-phase approximation
Within the framework of density-functional theory, the correlation energy is computed in the random-phase approximation (RPA) using spinors obtained from a two-component relativistic Kohn–Sham calculation accounting for spin–orbit interactions. Ring-coupled-cluster equations are solved to obtain the two-component RPA correlation energy. Results are presented for the hydrides of the halogens Br, I, and At as well as of the coinage metals Cu, Ag, and Au, based on two-component relativistic exact-decoupling Kohn–Sham calculations.
Systematic theoretical investigation of the zero-field splitting in Gd(III) complexes: Wave function and density functional approaches
The zero-field splitting (ZFS) of the electronic ground state in paramagnetic ions is a sensitive probe of the variations in the electronic and molecular structure with an impact on fields ranging from fundamental physical chemistry to medical applications. A detailed analysis of the ZFS in a series of symmetric Gd(III) complexes is presented in order to establish the applicability and accuracy of computational methods using multiconfigurational complete-active-space self-consistent field wave functions and of density functional theory calculations. The various computational schemes are then applied to larger complexes Gd(III)DOTA(H{sub 2}O){sup −}, Gd(III)DTPA(H{sub 2}O){sup 2−}, and Gd(III)(H{sub 2}O){sub 8}{sup 3+} in order to analyze how the theoretical results compare to experimentally derived parameters. In contrast to approximations based on density functional theory, the multiconfigurational methods produce results for the ZFS of Gd(III) complexes on the correct order of magnitude.
A Short Note on the Scaling Function Constant Problem in the Two-Dimensional Ising Model
We provide a simple derivation of the constant factor in the short-distance asymptotics of the tau-function associated with the 2-point function of the two-dimensional Ising model. This factor was first computed by Tracy (Commun Math Phys 142:297–311, 1991) via an exponential series expansion of the correlation function. Further simplifications in the analysis are due to Tracy and Widom (Commun Math Phys 190:697–721, 1998) using Fredholm determinant representations of the correlation function and Wiener–Hopf approximation results for the underlying resolvent operator. Our method relies on an action integral representation of the tau-function and asymptotic results for the underlying Painlevé-III transcendent from McCoy et al. (J Math Phys 18:1058–1092, 1977).
Modeling and characterizing anisotropic inclusion orientation in heterogeneous material via directional cluster functions and stochastic microstructure reconstruction
We present a framework to model and characterize the microstructure of heterogeneous materials with anisotropic inclusions of secondary phases based on the directional correlation functions of the inclusions. Specifically, we have devised an efficient method to incorporate both directional two-point correlation functions S{sub 2} and directional two-point cluster functions C{sub 2} that contain non-trivial topological connectedness information into the simulated annealing microstructure reconstruction procedure. Our framework is applied to model an anisotropic aluminum alloy and the accuracy of the reconstructed structural models is assessed by quantitative comparison with the actual microstructure obtained via x-ray tomography. We show that incorporation of directional clustering information via C{sub 2} significantly improves the accuracy of the reconstruction. In addition, a set of analytical “basis” correlation functions are introduced to approximate the actual S{sub 2} and C{sub 2} of the material. With the proper choice of basis functions, the anisotropic microstructure can be represented by a handful of parameters including the effective linear sizes of the iron-rich and silicon-rich inclusions along three orthogonal directions. This provides a general and efficient means for heterogeneous material modeling that enables one to significantly reduce the data set required to characterize the anisotropic microstructure.
Structural and orientation effects on electronic energy transfer between silicon quantum dots with dopants and with silver adsorbates
Starting from the atomic structure of silicon quantum dots (QDs), and utilizing ab initio electronic structure calculations within the Förster resonance energy transfer (FRET) treatment, a model has been developed to characterize electronic excitation energy transfer between QDs. Electronic energy transfer rates, K{sub EET}, between selected identical pairs of crystalline silicon quantum dots systems, either bare, doped with Al or P, or adsorbed with Ag and Ag{sub 3}, have been calculated and analyzed to extend previous work on light absorption by QDs. The effects of their size and relative orientation on energy transfer rates for each system have also been considered. Using time-dependent density functional theory and the hybrid functional HSE06, the FRET treatment was employed to model electronic energy transfer rates within the dipole-dipole interaction approximation. Calculations with adsorbed Ag show that: (a) addition of Ag increases rates up to 100 times, (b) addition of Ag{sub 3} increases rates up to 1000 times, (c) collinear alignment of permanent dipoles increases transfer rates by an order of magnitude compared to parallel orientation, and (d) smaller QD-size increases transfer due to greater electronic orbitals overlap. Calculations with dopants show that: (a) p-type and n-type dopants enhance energy transfer up to two orders of magnitude, (b) surface-doping with P and center-doping with Al show the greatest rates, and (c) K{sub EET} is largest for collinear permanent dipoles when the dopant is on the outer surface and for parallel permanent dipoles when the dopant is inside the QD.
Improved energies and local energies with weighted variational Monte Carlo
Neural network parametrizations have increasingly been used to represent the ground and excited states in variational Monte Carlo (VMC) with promising results. However, traditional VMC methods only optimize the wave function in regions of peak probability. The wave function is uncontrolled in the tails of the probability distribution, which can limit the accuracy of the trained wave function. To improve the approximation accuracy in the probability tails, this paper interprets VMC as a gradient flow in the space of wave functions, followed by a projection step. From this perspective, arbitrary probability distributions can be used in the projection step, allowing the user to prioritize accuracy in different regions of state space. Motivated by this theoretical perspective, the paper tests a weighted VMC method on the antiferromagnetic Heisenberg model for a periodic spin chain. Compared to traditional VMC, weighted VMC reduces the error in the ground state energy by a factor of 2, and it reduces the errors in the local energies away from the mode by large factors of 10 2 –10 4 .
Semi-classical approximation description of static properties of nuclei
In this report, we present recent development in semi-classical description of static properties of nuclei. By employing the Wigner transform, we derive simple semiclassical approximations for evaluating properties of nuclear system. We present results of calculations of static properties, using the energy density functional associated with Skyrme effective nucleon–nucleon interaction. In particular, we consider properties of the Wigner distribution function, line of beta stability, Coulomb and symmetry energies, deformation energies, equation of state, nuclear radii, and single-particle level density.
Variational formulation of Budyko-Sellers climate models
A class of simple climate models including those of the Budyko-Sellers type are formulated from a variational principle. A functional is constructed for the zonally averaged mean annual temperature field such that extrema of the functional occur when the climate satisfies the usual energy-balance equation. Local minima of the functional correspond to stable solutions while saddle points correspond to unstable solutions. The technique can be used to construct approximate solutions from trial functions and to carry out finite-amplitude stability analyses. A spectral example is given in explicit detail.
Ensemble density variational methods with self- and ghost-interaction-corrected functionals
Ensemble density functional theory (DFT) offers a way of predicting excited-states energies of atomic and molecular systems without referring to a density response function. Despite a significant theoretical work, practical applications of the proposed approximations have been scarce and they do not allow for a fair judgement of the potential usefulness of ensemble DFT with available functionals. In the paper, we investigate two forms of ensemble density functionals formulated within ensemble DFT framework: the Gross, Oliveira, and Kohn (GOK) functional proposed by Gross et al. [Phys. Rev. A 37, 2809 (1988)] alongside the orbital-dependent eDFT form of the functional introduced by Nagy [J. Phys. B 34, 2363 (2001)] (the acronym eDFT proposed in analogy to eHF – ensemble Hartree-Fock method). Local and semi-local ground-state density functionals are employed in both approaches. Approximate ensemble density functionals contain not only spurious self-interaction but also the so-called ghost-interaction which has no counterpart in the ground-state DFT. We propose how to correct the GOK functional for both kinds of interactions in approximations that go beyond the exact-exchange functional. Numerical applications lead to a conclusion that functionals free of the ghost-interaction by construction, i.e., eDFT, yield much more reliable results than approximate self- and ghost-interaction-corrected GOK functional. Additionally, local density functional corrected for self-interaction employed in the eDFT framework yields excitations energies of the accuracy comparable to that of the uncorrected semi-local eDFT functional.
Structural and opto-electronic properties of 2D AlSb monolayer
We have investigated dielectric function related optical properties such as refractive index, absorption coefficient of two-dimensional hexagonal system of aluminum antimony (AlSb). We have also find structural and electronic properties of AlSb which show direct/indirect band gap with planar structure, employing the density functional theory using the generalized gradient approximation (GGA) given by Perdew-Burke-Ernzerhof (PBE) functional for exchange-correlation potential. The refractive index n(ω) increases with frequency in the near infrared region but in visible region n(ω) increasing after decrease.
High temperature ablative foam
An ablative foam composition is formed of approximately 150 to 250 parts by weight polymeric isocyanate having an isocyanate functionality of 2.6 to 3.2; approximately 15 to 30 parts by weight reactive flame retardant having a hydroxyl number range from 200-260; approximately 10 to 40 parts by weight non-reactive flame retardant; approximately 10 to 40 parts by weight nonhydrolyzable silicone copolymer having a hydroxyl number range from 75-205; and approximately 3 to 16 parts by weight amine initiated polyether resin having an isocyanate functionality greater than or equal to 3.0 and a hydroxyl number range from 400-800.
Constructing tensor network influence functionals for general quantum dynamics
Here, we describe an iterative formalism to compute influence functionals that describe the general quantum dynamics of a subsystem beyond the assumption of linear coupling to a quadratic bath. We use a space-time tensor network representation of the influence functional and investigate its approximability in terms of its bond dimension and time-like entanglement in the tensor network description. We study two numerical models, the spin-boson model and a model of interacting hard-core bosons in a 1D harmonic trap. We find that the influence functional and the intermediates involved in its construction can be efficiently approximated by low bond dimension tensor networks in certain dynamical regimes, which allows the quantum dynamics to be accurately computed for longer times than with direct time evolution methods. However, as one iteratively integrates out the bath, the correlations in the influence functional can first increase before decreasing, indicating that the final compressibility of the influence functional is achieved via non-trivial cancellation.
Investigation of structural, electronic, elastic and optical properties of Cd{sub 1-x-y}Zn{sub x}Hg{sub y}Te alloys
Structural, optical and electronic properties and elastic constants of Cd1{sub -x-y}Zn{sub x} Hg{sub y}Te alloys have been studied by employing the commercial code Castep based on density functional theory. The generalized gradient approximation and local density approximation were utilized as exchange correlation. Using elastic constants for compounds, bulk modulus, band gap, Fermi energy and Kramers–Kronig relations, dielectric constants and the refractive index have been found through calculations. Apart from these, X-ray measurements revealed elastic constants and Vegard’s law. It is seen that results obtained from theory and experiments are all in agreement.
Simplified computational approach for dual-probe heat-pulse method
Two equations are currently available for estimating soil volumetric heat capacity (pc) with the dual-probe heat-pulse (DPHP) method. One is simple but gives only approximate results because it assumes that the DPHP sensor releases an impulse of heat instantaneously. The other explicitly accounts for the finite duration of heating and gives exact results. Unfortunately, the equation that gives exact results involves the exponential integral function, which is not available in most computer spreadsheet software packages or data logger function libraries. In this note we introduce an approximation of the exact equation that contains only simple algebraic functions. The approximation consists of the first five terms of a Taylor series, which are written as a telescoped polynomial for computational purposes. For most applications of the DPHP method, the polynomial approximation gives estimates of pc that are at least an order of magnitude more accurate than estimates obtained from the simple equation based on instantaneous heating.