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 19 records

Influence of the complete basis set approximation, tight weighted-core, and diffuse functions on the DLPNO-CCSD(T1) atomization energies of neutral H,C,O-compounds

The impact of complete basis set extrapolation schemes (CBS), diffuse functions, and tight weighted core functions on enthalpies of formation predicted via the DLPNO-CCSD(T1) reduced Feller Peterson-Dixon approach has been examined for neutral H,C,O-compounds. All tested three-point (TZ/QZ/5Z) extrapolation schemes result in mean unsigned deviation (MUD) below 2 kJ mol -1 relative to the experiment. Here, the two-point QZ/5Z and TZ/QZ CBS 1/$l^3_{max}$ extrapolation schemes are inferior to their inverse power counterpart (1/($l_{max}$ + 1/2) 4 ) by 1.3 and 4.3 kJ mol -1 . The CBS extrapolated frozen core atomization energies are insensitive (within 1 kJ mol -1 ) to augmentation of the basis set with tight weighted core functions. The core-valence correlation effects converge already at triple-ζ, although double- ζ/triple- ζ CBS extrapolation performs better and is recommended. The effect of diffuse function augmentation converges slowly, and cannot be reproduced with double- ζ or triple- ζ calculations as these are plagued with basis set superposition and incompleteness errors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The importance of diffuse f functions for transition metals

The importance of a diffuse f function for determining the dissociation energy (De) of Cu2 and the 3F-3D energy separation in Ni atom is investigated. It is found that the diffuse f contributes at most 0.05 eV to the De of Cu2 when added to a basis containing tight f functions and a flexibly contracted d basis. The diffuse f function is found to decrease the 3F-3D separation in Ni, but by substantially less than the tight f functions.

Langhoff, S. R.↗

A photometric function for diffuse reflection by particulate materials

A photometric function is proposed to describe the diffuse reflection of radiation by particulate materials. Both multiple scattering and the dominant effects of particle shadowing are included and the function is verified by comparisons with the photometries of laboratory surfaces. Brightness measurements of planetary and other diffusely scattering surfaces can be used to calculate the brightness for geometries other than those used in the measurements and for which the Minnaert function does not apply. The measurements also can be directly related to such surface characteristics as particle size, single-particle albedo, and compactness.

Meador, W. E.↗

Diffusion-convection function of cosmic rays

The fundamental properties and some numerical results of the solution of the diffusion equation of an impulsive cosmic-ray point source in an uniform, unbounded and spherically symmetrical moving medium is presented. The diffusion-convection(D-C) function is an elementary composite function of the solution of the D-C equation for the particles injected impulsively from a diffusive point source into the medium. It is the analytic solution derived by the dimensional method for the propagation equation of solar cosmic rays in the heliosphere, i.e. the interplanetary space. Because of the introduction of convection effect of solar wind, a nonhomogeneous term appears in the propagation equation, it is difficult to express its solution in terms of the ordinary special functions. The research made so far has led to a solution containing only the first order approximation of the convection effect.

Zhang, G.↗

Polarizabilities of neutral atoms and atomic ions with a noble gas electron configuration

Atomic polarizabilities play an important role in the development of force fields for molecular simulations, as well as for the development of qualitative concepts of atomic and molecular behavior. Coupled cluster theory at the coupled cluster singles doubles triples level with very large correlation-consistent basis sets with extended diffuse functions has been used to predict the polarizabilities of the atomic neutrals, mono-cations and mono-anions with a noble gas configuration. Additional corrections for scalar relativistic and spin–orbit effects were also included for the electron configurations of Kr, Xe, and Rn. The results are in excellent agreement with experiment or with other high level calculations where available. The current results for most of these species represent the best available values for the polarizabilities. The results show that the polarizability of H - is very difficult to calculate without extremely diffuse functions. The polarizability of H - is the largest value, 34.05 Å 3 , calculated for all species in the current study. The polarizabilities of the remaining halogen anions are also the best available values. The polarizabilities of the halogen anions (excluding F - ) and H - have a linear correlation with the electron affinity of the neutral atom. Spin– orbit effects, even for closed shell species, cannot be ignored for quantitative accuracy, and the inclusion of spin–orbit effects for Fr + , Rn, and At - increases the polarizability by 4%, 6%, and 15%, respectively.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Dissociation Energies of He2, HeH, and ArH; A Bond Function Study

The bond energies and bond lengths are determined for He2, HeH, and ArH at the CCSD(T) level using both atom-centered basis sets and those that include bond functions. The addition of bond functions dramatically improves the rate of convergence of the results with respect to the size of the atom-centered basis set; with bond functions, triple zeta atom-centered basis set, outperform quintuple zeta basis sets without bond functions. The addition of bond functions also reduces the number of diffuse functions that must be added to the atom-centered sets. Employing bond functions appear to offer a very cost effective method of computing the interaction between weakly bound systems, especially for He.

Bauschlicher, Charles W., Jr.↗

Effect of the Reconnection Electric Field on Electron Distribution Functions in the Diffusion Region of Magnetotail Reconnection

Electron distribution functions in the electron diffusion region during symmetric magnetic reconnection are investigated by means of theory and fully kinetic simulations. Crescent‐like striations are formed in distribution functions in the velocity plane perpendicular to the magnetic field. Using an analytical current sheet, we solve the equation of motion for electrons, and derive the shape of a crescent distribution, as a function of the distance from the neutral line, field gradients, and the reconnection electric field. Each crescent is tilted in the velocity plane because of the acceleration by the reconnection electric field, and multiple stripes appear due to multiple meandering bounces. Applying the theory to distribution functions observed in Earth's magnetotail, we deduce the amplitude of the reconnection electric field.

Bessho, N.↗

The binding of atomic hydrogen on graphene from density functional theory and diffusion Monte Carlo calculations

In this research, density functional theory (DFT) and diffusion Monte Carlo (DMC) methods are used to calculate the binding energy of a H atom chemisorbed on the graphene surface. The DMC value of the binding energy is about 16% smaller in magnitude than the Perdew–Burke–Ernzerhof (PBE) result. The inclusion of exact exchange through the use of the Heyd–Scuseria–Ernzerhof functional brings the DFT value of the binding energy closer in line with the DMC result. It is also found that there are significant differences in the charge distributions determined using PBE and DMC approaches.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The application of far IR lasers to the investigation of the optical properties of materials

The effect of diffusion on the saturation and absorption properties of molecular gases is considered. The analysis starts with a simple derivation of the multi-level rate equations. Incorporating the diffusion term, calculations are carried out both for the homogeneous and inhomogeneous broadening cases. Both graphical and analytical solutions are obtained for the absorption coefficient and the saturation intensity. It is shown that the effect of diffusion is characterized by a diffusion function, which varies inversely with the radius of the beam at low pressures and small beam sizes. The analysis is done in terms of two parameters k1a and k2a, which are simple functions of the intensity of the beam I, the effective relaxation rate w, and the effective diffusion constant D the a is designated as the radius of the beam.

Stafsudd, O. M.↗

Electron Distribution Functions in the Diffusion Region of Asymmetric Magnetic Reconnection

We study electron distribution functions in a diffusion region of antiparallel asymmetric reconnection by means of particle-in-cell simulations and analytical theory. At the electron stagnation point, the electron distribution comprises a crescent-shaped population and a core component. The crescent-shaped distribution is due to electrons coming from the magnetosheath toward the stagnation point and accelerated mainly by electric field normal to the current sheet. Only a part of magnetosheath electrons can reach the stagnation point and form the crescent-shaped distribution that has a boundary of a parabolic curve. The penetration length of magnetosheath electrons into the magnetosphere is derived. We expect that satellite observations can detect crescent-shaped electron distributions during magnetopause reconnection.

Electron↗

Cation–Ligand Interactions Dictate Salt Partitioning and Diffusivity in Ligand-Functionalized Polymer Membranes

Membranes are an attractive alternative to current thermal separations due to their scalability and energy efficiency in desalinating water. Unfortunately, many of the conventional membrane materials available today are unable to differentiate between ionic solutes, especially alkali cations, compromising their use in ion–ion separations. Inspired by the ion-specific interactions exhibited by biological ion channels, recent research efforts have focused on synthesizing and characterizing new polymeric materials that incorporate ligands into polymer networks to bias solubility and/or diffusivity of one cationic species over another. Despite these efforts, little is known about the influence of incorporating ligands into polymer membranes on solubility and diffusivity of the complexing species. In this study, we first build a qualitative model of salt partitioning, diffusivity, and permeability in generic cation-complexing ligand-functionalized polymer membranes. Next, to validate our model and hypotheses, we perform atomistic molecular dynamics simulations of a 12-crown-4-functionalized membrane in the presence of alkali halide salts at low concentration. Generally, cation complexation enhances cation solubility but decreases diffusivity. Interestingly, the reduction in diffusivity is predicted to be larger than the enhancement in solubility for materials which operate by the mechanisms proposed in our physical picture, ultimately resulting in a reduction in the permeability of the selectively complexing ion.

36 MATERIALS SCIENCE↗

Computations of ion diffusion coefficients from the Boltzmann-Fokker-Planck equation

The Boltzmann-Fokker-Planck equation is solved with the Chapman-Enskog method of analysis for the velocity distribution functions of helium, carbon, nitrogen, and oxygen. The analysis is a perturbation scheme based on the assumption of a collision-dominated gas, and the calculations are carried out to first order. The elements considered are treated as trace constituents in an electron-proton gas. From the resulting distribution functions, diffusion coefficients are computed which are found to be 20-30% less than those obtained by Chapman and Burgers. In addition, it is shown that the return current of cold electrons needed to maintain quasi-neutrality in a plasma with a temperature gradient contributes a term in the thermal diffusion coefficient omitted erroneously in previous works. This added term resolves the longstanding controversy over the discrepancy between the coefficients of Chapman and Burgers, which are seen to be completely equivalent in the light of this analysis. The viscosity coefficient for an electron-proton gas is also computed and found to be 7% less than that obtained by Braginskii.

Roussel-Dupre, R.↗

Reliable Spectroscopic Constants for CCH-, NH2- and Their Isotopomers from an Accurate Potential Energy Function

Accurate quartic force fields have been determined for the CCH- and NH2- molecular anions using the singles and doubles coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations, CCSD(T). Very large one-particle basis sets have been used including diffuse functions and up through g-type functions. Correlation of the nitrogen and carbon core electrons has been included, as well as other "small" effects, such as the diagonal Born-Oppenheimer correction, and basis set extrapolation, and corrections for higher-order correlation effects and scalar relativistic effects. Fundamental vibrational frequencies have been computed using standard second-order perturbation theory as well as variational methods. Comparison with the available experimental data is presented and discussed. The implications of our research for the astronomical observation of molecular anions will be discussed.

Lee, Timothy J.↗

FunDiff: diffusion models over function spaces for physics-informed generative modeling

Recent advances in generative modeling-particularly diffusion models and flow matching-have been widely used for synthesizing discrete data such as images and videos. However, adapting these models to physical applications remains challenging, as the quantities of interest are continuous functions governed by complex physical laws. To address this, we introduce FunDiff, an efficient and robust framework for generative modeling in function spaces. FunDiff combines a latent diffusion process with a function autoencoder architecture to handle input functions with varying discretizations, generates continuous functions that can be evaluated at arbitrary locations, and seamlessly incorporate physical priors. These priors are enforced through architectural constraints or physics-informed loss functions, ensuring that generated samples satisfy fundamental physical laws. We theoretically establish minimax optimality guarantees for density estimation in function spaces, demonstrating that diffusion-based estimators achieve optimal convergence rates under suitable regularity conditions. We further demonstrate the practical effectiveness of FunDiff across diverse applications in fluid dynamics and solid mechanics. Empirical results indicate that our method can generate physically consistent samples with high fidelity to the target distribution, and exhibit robustness to noisy and low-resolution data.

Wang, Sifan [Yale University, New Haven, CT (Unite↗

Comprehensive Basis-Set Testing of Extended Symmetry-Adapted Perturbation Theory and Assessment of Mixed-Basis Combinations to Reduce Cost

Hybrid or “extended” symmetry-adapted perturbation theory (XSAPT) replaces traditional SAPT’s treatment of dispersion with better performing alternatives while at the same time extending two-body (dimer) SAPT to a many-body treatment of polarization using a self-consistent charge embedding procedure. The present work presents a systematic study of how XSAPT interaction energies and energy components converge with respect to the choice of Gaussian basis set. Errors can be reduced in a systematic way using correlation-consistent basis sets, with aug-cc-pVTZ results converged within <0.1 kcal/mol. Similar (if slightly less systematic) behavior is obtained using Karlsruhe basis sets at much lower cost, and we introduce new versions with limited augmentation that are even more efficient. Pople-style basis sets, which are more efficient still, often afford good results if a large number of polarization functions are included. The dispersion models used in XSAPT afford much faster basis-set convergence as compared to the perturbative description of dispersion in conventional SAPT, meaning that “compromise” basis sets (such as jun-cc-pVDZ) are no longer required and benchmark-quality results can be obtained using triple-ζ basis sets. The use of diffuse functions proves to be essential, especially for the description of hydrogen bonds. As a result, the “δ(Hartree–Fock)” correction for high-order induction can be performed in double-ζ basis sets without significant loss of accuracy, leading to a mixed-basis approach that offers 4× speedup over the existing (cubic scaling) XSAPT approach.

74 ATOMIC AND MOLECULAR PHYSICS↗

The effect of hydrogen-bonding on flavin’s infrared absorption spectrum

Cluster and continuum solvation computational models are employed to model the effect of hydrogen bonding interactions on the vibrational modes of lumiflavin. Calculated spectra were compared to experimental Fourier-transform infrared (FTIR) spectra in the diagnostic 1450–1800 cm –1 range, where intense V C=C , V C=N , V C 2 =O , and V C 4 =O stretching modes of flavin’s isoalloxazine ring are found. Local mode analysis is used to describe the strength of hydrogen-bonding in cluster models. The computations indicate that V C=C and V C=N mode frequencies are relatively insensitive to intermolecular interactions while the V C 2 =O and V C 4 =O modes are sensitive to direct (and also indirect for V C 2 =O ) hydrogen-bonding interactions. Although flavin is neutral, basis sets without the diffuse functions provide incorrect relative frequencies and intensities. The 6-31+G* basis set is found to be adequate for this system, and there is limited benefit to considering larger basis sets. Calculated vibrational mode frequencies agree with experimentally determined frequencies in solution when cluster models with multiple water molecules are used. Accurate simulation of relative FTIR band intensities, on the other hand, requires a continuum (or possibly quantum mechanical/molecular mechanical) model that accounts for long-range electrostatic effects. Lastly, an experimental peak at ca. 1624 cm –1 that is typically assigned to the V C 2 =O vibrational stretching mode has a complicated shape that suggests multiple underlying contributions. Our calculations show that this band has contributions from both the C 6 –C 7 and C 2 = O stretching vibrations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Effective Fragment Molecular Orbital Method: Achieving High Scalability and Accuracy for Large Systems

The effective fragment molecular orbital (EFMO) method has been developed to predict the total energy of a very large molecular system accurately (with respect to the underlying quantum mechanical method) and efficiently by taking advantage of the locality of strong chemical interactions and employing a two-level hierarchical parallelism. The accuracy of the EFMO method is partly attributed to the accurate and robust intermolecular interaction prediction between distant fragments, in particular, the many-body polarization and dispersion effects, which require the generation of static and dynamic polarizability tensors by solving the coupled perturbed Hartree–Fock (CPHF) and time-dependent HF (TDHF) equations, respectively. Solving the CPHF and TDHF equations is the main EFMO computational bottleneck due to the inefficient (serial) and I/O-intensive implementation of the CPHF and TDHF solvers. In this work, the efficiency and scalability of the EFMO method are significantly improved with a new CPU memory-based implementation for solving the CPHF and TDHF equations that are parallelized by either message passing interface (MPI) or hybrid MPI/OpenMP. Here, the accuracy of the EFMO method is demonstrated for both covalently bonded systems and noncovalently bound molecular clusters by systematically examining the effects of basis sets and a key distance-related cutoff parameter, R cut . R cut determines whether a fragment pair (dimer) is treated by the chosen ab initio method or calculated using the effective fragment potential (EFP) method (separated dimers). Decreasing the value of Rcut increases the number of separated (EFP) dimers, thereby decreasing the computational effort. It is demonstrated that excellent accuracy (<1 kcal/mol error per fragment) can be achieved when using a sufficiently large basis set with diffuse functions coupled with a small R cut value. With the new parallel implementation, the total EFMO wall time is substantially reduced, especially with a high number of MPI ranks. Given a sufficient workload, nearly ideal strong scaling is achieved for the CPHF and TDHF parts of the calculation. For the first time, EFMO calculations with the inclusion of long-range polarization and dispersion interactions on a hydrated mesoporous silica nanoparticle with explicit water solvent molecules (more than 15k atoms) are achieved on a massively parallel supercomputer using nearly 1000 physical nodes. In addition, EFMO calculations on the carbinolamine formation step of an amine-catalyzed aldol reaction at the nanoscale with explicit solvent effects are presented.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗