Search NASA⌕ Search

SEARCH · Search NASA

Results for “diffuse functions”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Benchmarking Correlation-Consistent Basis Sets for Frequency-Dependent Polarizabilities with Multiresolution Analysis

This paper presents the first converged frequency-dependent HF polarizability results for general molecules, on a set of 89 closed-shell atoms and molecules. The solver employs multiresolution analysis (MRA) in a multiwavelet basis to compute both ground and response states to a guaranteed precision, which are validated against independent numerical grid calculations on atoms and linear molecules. The MRA ground-state energies and response properties are used to evaluate results in correlation-consistent basis sets up to 5Z augmented with either single or double diffuse functions and core-polarization functions. Systematic trends are revealed through consideration of chemical composition as well as the use of machine learning to cluster convergence trends, the latter suggesting the possibility of learning and correcting basis-set error.

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

Learning effective stochastic differential equations from microscopic simulations: Linking stochastic numerics to deep learning

We identify effective stochastic differential equations (SDEs) for coarse observables of fine-grained particle- or agent-based simulations; these SDEs then provide useful coarse surrogate models of the fine scale dynamics. We approximate the drift and diffusivity functions in these effective SDEs through neural networks, which can be thought of as effective stochastic ResNets. The loss function is inspired by, and embodies, the structure of established stochastic numerical integrators (here, Euler–Maruyama and Milstein); our approximations can thus benefit from backward error analysis of these underlying numerical schemes. They also lend themselves naturally to “physics-informed” gray-box identification when approximate coarse models, such as mean field equations, are available. Existing numerical integration schemes for Langevin-type equations and for stochastic partial differential equations can also be used for training; we demonstrate this on a stochastically forced oscillator and the stochastic wave equation. Our approach does not require long trajectories, works on scattered snapshot data, and is designed to naturally handle different time steps per snapshot. We consider both the case where the coarse collective observables are known in advance, as well as the case where they must be found in a data-driven manner.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Band structures and Z 2 invariants of two-dimensional transition metal dichalcogenide monolayers from fully relativistic Dirac-Kohn-Sham theory using Gaussian-type orbitals

Two-dimensional (2D) materials exhibit a wide range of remarkable phenomena, many of which owe their existence to the relativistic spin-orbit coupling (SOC) effects. To understand and predict properties of materials containing heavy elements, such as the transition-metal dichalcogenides (TMDs), relativistic effects must be taken into account in first-principles calculations. We present an all-electron method based on the four-component Dirac Hamiltonian and Gaussian-type orbitals (GTOs) that overcomes complications associated with linear dependencies and ill-conditioned matrices that arise when diffuse functions are included in the basis. Until now, there has been no systematic study of the convergence of GTO basis sets for periodic solids either at the nonrelativistic or the relativistic level. Here, in this study, we provide such a study of relativistic band structures of the 2D TMDs in the hexagonal (2H), tetragonal (1T), and distorted tetragonal (1T') structures, along with a discussion of their SOC-driven properties (Rashba splitting and Z 2 topological invariants). We demonstrate the viability of our approach even when large basis sets with multiple basis functions involving various valence orbitals (denoted triple- and quadruple-ζ) are used in the relativistic regime. Our method does not require the use of pseudopotentials and provides access to all electronic states within the same framework. Our study paves the way for direct studies of material properties, such as the parameters in spin Hamiltonians, that depend heavily on the electron density near atomic nuclei where relativistic and SOC effects are the strongest.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1971

Ablation and radiation coupled viscous hypersonic shock layers, volume 1

The results for a stagnation-line analysis of the radiative heating of a phenolic-nylon ablator are presented. The analysis includes flow field coupling with the ablator surface, equilibrium chemistry, a step-function diffusion model and a coupled line and continuum radiation calculation. This report serves as the documentation, i e. users manual and operating instructions for the computer programs listed in the report.

Engel, C. D.↗
NASA NTRS1985

A potential-energy surface study of the 2A1 and low-lying dissociative states of the methoxy radical

Accurate, ab initio quantum chemical techniques are applied in the present study of low lying bound and dissociative states of the methoxy radical at C3nu conformations, using a double zeta quality basis set that is augmented with polarization and diffuse functions. Excitation energy estimates are obtained for vertical excitation, vertical deexcitation, and system origin. The rate of methoxy photolysis is estimated to be too small to warrant its inclusion in atmospheric models.

Jackels, C. F.↗
NASA NTRS1985

N2(+) bound quartet and sextet state potential energy curves

The N2(+) potential energies have been determined from a complete active space self-consistent field calculation with active 2s and 2p electrons. A (6s 4p 3d 1f) Gaussian basis set was used together with additional higher angular momentum and diffuse functions. The calculated potential energy curves for the states 4Sigma(mu)(+), 4Pi(g), and 6Sigma(g)(+), for which there are no spectroscopic observations, are presented. The corresponding spectroscopic constants have been determined from a polynomial curve fit to the computed energies near the well minima and are shown. The 6Sigma(g)(+) state is found to be significantly bound, with a minimum at 1.72 A.

Partridge, H.↗
NASA NTRS1987

Potential surfaces for O atom-polymer reactions

Ab initio quantum chemistry methods are used to study the energetics of interactions of O atoms with organic compounds. Polyethylene (CH2)n has been chosen as the model system to study the interactions of O(3P) and O(1D) atoms with polymers. In particular, H abstraction is investigated and polyethylene is represented by a C3 (propane) oligomeric model. The gradient method, as implemented in the GRADSCF package of programs, is used to determine the geometries and energies of products and reactants. The saddle point, barrier geometry is determined by minimizing the squares of the gradients of the potential with respect to the internal coordinates. To correctly describe the change in bonding during the reaction at least a two configuration MCSCF (multiconfiguration self consistent field) or GVB (generalized valence bond) wave function has to be used. Basis sets include standard Pople and Dunning sets, however, increased with polarization functions and diffuse p functions on both the C and O atoms. The latter is important due to the O(-) character of the wave function at the saddle point and products. Normal modes and vibrational energy levels are given for the reactants, saddle points and products. Finally, quantitative energetics are obtained by implementing a small CAS (complete active space) approach followed by limited configuration interaction (CI) calculations. Comparisons are made with available experimental data.

Laskowski, B. C.↗
NASA NTRS1987

Effective core potential parameters for first- and second-row atoms

An improved effective core potential (ECP) technique is described and used to give ECP parameters for the atoms of the first two rows of the periodic table. A given basis set is parametrized which allows for a direct comparison with all-electron calculations. Extensive test calculations on first- and second-row molecules using the ECP have been performed, giving excellent agreement with the all-electron results at the SCF, CASSCF, and CI levels of treatment. Correlating and diffuse functions may be added without modifying the ECP parameters. The present ECP descriptions result in CPU time reductions of the order of 50 percent in addition to reduced disk storage.

Pettersson, Lars G. M.↗
NASA NTRS1991

Theoretical study of the low-lying bound states of O2

It is demonstrated that a complete-active-space self-consistent-field (CASSCF) (2p)/MRCI + Q (multireference configuration interaction with a Davidson correction) description in a (13s8p6d 4f2g)/((5s4p3d 2f1g) atomic natural orbits (ANO) basis set supplemented with diffuse functions provides a quantitative description of the six lowest states of O2. The calculated potentials are within 0.05 eV (1.2 kilocal/mol) of accurate experimental results. The importance of substantially expanding the primitive basis set has been investigated, and it is demonstrated that such expansions yield insignificant improvement in the spectroscopic constants. Potential energy curves have also been reported for the weakly bound states of O2. The 5Pi(g) state is estimated to have a D(e) of 0.16 +/- 0.03 eV. The upper bound of D(e) is found to be sufficiently large that the importance of this state as a precursor for the formation of O2 (b 1Sigma(t)(+)) and O(1S) should be reconsidered.

Partridge, Harry↗
NASA NTRS1992

Comment on broadening of water microwave lines by collisions with helium atoms

Earlier theoretical descriptions of interaction forces and collision dynamics in the H2O-He system (Palma et al., 1988; Green et al., 1991) are improved by calculating the interaction at additional orientations and by including diffuse functions in the molecular orbital expansion basis set. The improvements make substantial changes to the potential energy surface, but make only small changes in predicted microwave pressure broadening coefficients. The qualitative temperature dependence and relative widths of different lines appear to be converged with respect to successive improvements to the theoretical description of this system. Remaining discrepancies with available data must reflect problems with the experimental values.

Maluendes, Sergio↗
NASA NTRS2000

Postural control model interpretation of stabilogram diffusion analysis

Collins and De Luca [Collins JJ. De Luca CJ (1993) Exp Brain Res 95: 308-318] introduced a new method known as stabilogram diffusion analysis that provides a quantitative statistical measure of the apparently random variations of center-of-pressure (COP) trajectories recorded during quiet upright stance in humans. This analysis generates a stabilogram diffusion function (SDF) that summarizes the mean square COP displacement as a function of the time interval between COP comparisons. SDFs have a characteristic two-part form that suggests the presence of two different control regimes: a short-term open-loop control behavior and a longer-term closed-loop behavior. This paper demonstrates that a very simple closed-loop control model of upright stance can generate realistic SDFs. The model consists of an inverted pendulum body with torque applied at the ankle joint. This torque includes a random disturbance torque and a control torque. The control torque is a function of the deviation (error signal) between the desired upright body position and the actual body position, and is generated in proportion to the error signal, the derivative of the error signal, and the integral of the error signal [i.e. a proportional, integral and derivative (PID) neural controller]. The control torque is applied with a time delay representing conduction, processing, and muscle activation delays. Variations in the PID parameters and the time delay generate variations in SDFs that mimic real experimental SDFs. This model analysis allows one to interpret experimentally observed changes in SDFs in terms of variations in neural controller and time delay parameters rather than in terms of open-loop versus closed-loop behavior.

NASA Program Biomedical Research and Countermeasur↗
NASA NTRS2014

A Radiation Chemistry Code Based on the Greens Functions of the Diffusion Equation

Ionizing radiation produces several radiolytic species such as.OH, e-aq, and H. when interacting with biological matter. Following their creation, radiolytic species diffuse and chemically react with biological molecules such as DNA. Despite years of research, many questions on the DNA damage by ionizing radiation remains, notably on the indirect effect, i.e. the damage resulting from the reactions of the radiolytic species with DNA. To simulate DNA damage by ionizing radiation, we are developing a step-by-step radiation chemistry code that is based on the Green's functions of the diffusion equation (GFDE), which is able to follow the trajectories of all particles and their reactions with time. In the recent years, simulations based on the GFDE have been used extensively in biochemistry, notably to simulate biochemical networks in time and space and are often used as the "gold standard" to validate diffusion-reaction theories. The exact GFDE for partially diffusion-controlled reactions is difficult to use because of its complex form. Therefore, the radial Green's function, which is much simpler, is often used. Hence, much effort has been devoted to the sampling of the radial Green's functions, for which we have developed a sampling algorithm This algorithm only yields the inter-particle distance vector length after a time step; the sampling of the deviation angle of the inter-particle vector is not taken into consideration. In this work, we show that the radial distribution is predicted by the exact radial Green's function. We also use a technique developed by Clifford et al. to generate the inter-particle vector deviation angles, knowing the inter-particle vector length before and after a time step. The results are compared with those predicted by the exact GFDE and by the analytical angular functions for free diffusion. This first step in the creation of the radiation chemistry code should help the understanding of the contribution of the indirect effect in the formation of DNA damage and double-strand breaks.

Plante, Ianik↗
DOE OSTI2021

Revisiting diffusion-based moisture functions: why do they fail?

Microbial activities in soil are controlled by the balance between physical and biological mechanisms that are currently implied by soil moisture-respiration relationships, but not explicitly represented. Diffusion based moisture functions focus on physical controls, while substrate kinetics based models primarily address biological drivers. How to harmonize these two types of models for robust predictions of soil respiration and other processes has not been addressed. Here we derived a generalizable analytical solution of diffusion-limited nonlinear microbial uptake kinetics, in which microbial metabolic rates were governed by both physical and biological drivers via substrate diffusion and microbial affinity. We proposed a soil trait concept to integrate both physical and biological perspectives into one unified function, and tested it using both laboratory and field based observations from a variety of soils. This approach provides a tractable and scalable approach that incorporates soil heterogeneity and functional complexity in macroscopic model descriptions.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2011

Step-by-Step Simulation of Radiation Chemistry Using Green Functions for Diffusion-Influenced Reactions

Radiolytic species are formed approximately 1 ps after the passage of ionizing radiation through matter. After their formation, they diffuse and chemically react with other radiolytic species and neighboring biological molecules, leading to various oxidative damage. Therefore, the simulation of radiation chemistry is of considerable importance to understand how radiolytic species damage biological molecules [1]. The step-by-step simulation of chemical reactions is difficult, because the radiolytic species are distributed non-homogeneously in the medium. Consequently, computational approaches based on Green functions for diffusion-influenced reactions should be used [2]. Recently, Green functions for more complex type of reactions have been published [3-4]. We have developed exact random variate generators of these Green functions [5], which will allow us to use them in radiation chemistry codes. Moreover, simulating chemistry using the Green functions is which is computationally very demanding, because the probabilities of reactions between each pair of particles should be evaluated at each timestep [2]. This kind of problem is well adapted for General Purpose Graphic Processing Units (GPGPU), which can handle a large number of similar calculations simultaneously. These new developments will allow us to include more complex reactions in chemistry codes, and to improve the calculation time. This code should be of importance to link radiation track structure simulations and DNA damage models.

Plante, Ianik↗
NASA NTRS1983

Supplemental basis functions for the second transition row elements

It is noted that for molecular calculations, the basis sets presented by Huzinaga (1977) need to be augmented by (1) 5p functions to describe the 5s - 5p near degeneracy; (2) a diffuse 4d function to provide for a balanced description of the 5s2 4dn, 5s1 4d(n+1), and 4d(n+2) states of the atom; and (3) a set of 4f functions to correlate the 4d functions. Here, the diffuse 4d function is similar in function to the diffuse 3d function for the first transition row elements recommended by Hay (1977). A table is included giving the optimized values for the diffuse 4d, the 5p, and 4f (STO exponent) functions. The diffuse 4d function and the 5p functions are optimized at the SCF level on the basis of the 5s1 4d(n+1) state (except for Pd, which is optimized for the 4d10 state) and the 5s1 5p1 4dn state, respectively. The table also gives the energies and the atomic symmetries for each of the SCF calculations.

Walch, S. P.↗
NASA NTRS1972

Relation of the diffuse reflectance remission function to the fundamental optical parameters.

The Kubelka-Munk equations describing the diffuse reflectance of a powdered sample were compared to equations obtained using a uniformly-sized rough-surfaced spherical particle model. The comparison resulted in equations relating the remission function and the Kubelka-Munk constants to the index of refraction, the absorption coefficient, and the average particle diameter of a powdered sample. Published experimental results were used to test the equation relating to the remission function to the fundamental optical parameters.

Simmons, E. L.↗
NASA NTRS2011

Step-by-Step Simulation of Radiation of Radiation Chemistry Using Green Functions for Diffusion-Influenced Reactions

The irradiation of biological systems leads to the formation of radiolytic species such as H(raised dot), (raised dot)OH, H2, H2O2, e(sup -)(sub aq), etc.[1]. These species react with neighboring molecules, which result in damage in biological molecules such as DNA. Radiation chemistry is there for every important to understand the radiobiological consequences of radiation[2]. In this work, we discuss an approach based on the exact Green Functions for diffusion-influenced reactions which may be used to simulate radiation chemistry and eventually extended to study more complex systems, including DNA.

Plante, Ianik↗
NASA NTRS2015

Role of Multiple Atmospheric Reflections in Formation of Electron Distribution Function in the Diffuse Aurora Region

The precipitation of high-energy magnetospheric electrons (E greater than 500-600 electronvolts) in the diffuse aurora contributes significant energy flux into Earth's ionosphere. In the diffuse aurora, precipitating electrons initially injected from the plasmasheet via wave-particle interaction processes degrade in the atmosphere toward lower energies and produce secondary electrons via impact ionization of the neutral atmosphere. These initially precipitating electrons of magnetospheric origin can be additionally reflected back into the magnetosphere by the two magnetically conjugated atmospheres, leading to a series of multiple reflections that can greatly influence the initially precipitating flux at the upper ionospheric boundary (700-800 kilometers) and the resultant population of secondary electrons and electrons cascading toward lower energies. We present the solution of the Boltzmann.Landau kinetic equation that uniformly describes the entire electron distribution function in the diffuse aurora, including the affiliated production of secondary electrons (E is less than or equal to 600 electronvolts) and their energy interplay in the magnetosphere and two conjugated ionospheres. This solution takes into account the role of multiple atmospheric reflections of the precipitated electrons that were initially moved into the loss cone via wave.particle interaction processes in Earth's plasmasheet.

Khazanov, George V.↗