SEARCH · Search NASA
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.
A Monte Carlo method for angle scattering of electrons in air that accommodates the small-angle multiple scattering and larger-angle single scattering limits is introduced. The algorithm is designed for use in a particle-in-cell simulation of electron transport and electromagnetic wave effects in air. The method is illustrated in example calculations.
Explore the source record for details and available documents.
CO electrooxidation on the Au(111) surface exhibits an onset potential that depends strongly on the pH of the electrolyte. In alkaline solution, the onset potential for CO electrooxidation is found, experimentally, to be 0.5 V lower than in acidic media on the reversible hydrogen electrode scale. This phenomenon is explained here with density functional theory which is used to calculate adsorbate binding energies including the electric double layer of the Au(111)/aqueous interface. Our model consists of a charged Au(111) slab and implicit solvation for the electrolyte. The double reference method is used to determine the potential-dependent CO electrooxidation reaction energetics. A microkinetic model, based upon the calculated reaction energetics, confi rms the lower onset potential in alkaline media. Our results show that there are three factors contributing to the lower onset potential in base: (1) stronger CO adsorption, (2) attraction between adsorbed CO and OH - , and (3) the high concentration of OH - in base.
Cu/SSZ-13 selective catalytic reduction catalysts display activity loss after mild hydrothermal aging without support degradation or loss of isolated Cu ion content. By applying electron paramagnetic resonance spectroscopy to compare hydrated, dehydrated and NH3-saturated catalysts before and after hydrothermal aging, as well as density functional theory calculations, it is found that hydrothermal aging induces Cu relocation, leading to stronger Cu-support interactions. This is particularly evidenced by ZCuOH conversion to Z2Cu, which leads to the generation of a unique structure with two Cu(II) ions positioned in a double 6-membered ring prism. The Cu(II) ions that do not undergo such dramatic changes also appear to stay closer to framework windows as evidenced by the overall increased anisotropy during EPR measurements. In situ EPR measurements under SCR conditions demonstrate that, in comparison to the fresh catalyst, higher percentages of Cu species stay as Cu(II) in the aged catalyst. This suggests that the reduction half-cycle of the SCR redox, i.e., Cu(II) ? Cu(I), becomes slower for the hydrothermally aged catalyst.
Explore the source record for details and available documents.
Abstract Classical turning surfaces of Kohn–Sham potentials separate classically allowed regions (CARs) from classically forbidden regions (CFRs). They are useful for understanding many chemical properties of molecules but need not exist in solids, where the density never decays to zero. At equilibrium geometries, we find that CFRs are absent in perfect metals, rare in covalent semiconductors at equilibrium, but common in ionic and molecular crystals. In all materials, CFRs appear or grow as the internuclear distances are uniformly expanded. They can also appear at a monovacancy in a metal. Calculations with several approximate density functionals and codes confirm these behaviors. A classical picture of conduction suggests that CARs should be connected in metals, and disconnected in wide-gap insulators, and is confirmed in the limits of extreme compression and expansion. Surprisingly, many semiconductors have no CFR at equilibrium, a key finding for density functional construction. Nonetheless, a strong correlation with insulating behavior can still be inferred. Moreover, equilibrium bond lengths for all cases can be estimated from the bond type and the sum of the classical turning radii of the free atoms or ions.
Abstract The repeated slab approach has become a de facto standard to accurately describe surface properties of materials by density functional theory calculations with periodic boundary conditions. For materials exhibiting spontaneous polarization, we show that the conventional scheme of passivation with pseudo hydrogen is unable to realize a charge-neutral surface. The presence of a net surface charge induces via Gauss’s law a macroscopic electric field through the slab and results in poor size convergence with respect to the thickness of the slab. We propose a modified passivation method that accounts for the effect of spontaneous polarization, describes the correct bulk limits and boosts convergence with respect to slab thickness. The robustness, reliability, and superior convergence of energetics and electronic structure achieved by the proposed method are demonstrated using the example of polar ZnO surfaces.
Abstract Point defects in hexagonal boron nitride (hBN) have attracted growing attention as bright single-photon emitters. However, understanding of their atomic structure and radiative properties remains incomplete. Here we study the excited states and radiative lifetimes of over 20 native defects and carbon or oxygen impurities in hBN using ab initio density functional theory and GW plus Bethe-Salpeter equation calculations, generating a large data set of their emission energy, polarization and lifetime. We find a wide variability across quantum emitters, with exciton energies ranging from 0.3 to 4 eV and radiative lifetimes from ns to ms for different defect structures. Through a Bayesian statistical analysis, we identify various high-likelihood charge-neutral defect emitters, among which the native V N N B defect is predicted to possess emission energy and radiative lifetime in agreement with experiments. Our work advances the microscopic understanding of hBN single-photon emitters and introduces a computational framework to characterize and identify quantum emitters in 2D materials.
We investigated the structure-dependent vibrational properties of different Mg(BH 4 ) 2 polymorphs (α, β, γ, and δ phases) with a combination of neutron vibrational spectroscopy (NVS) measurements and density functional theory (DFT) calculations, with emphasis placed on the effects of the local structure and orientation of the BH 4 - anions. DFT simulations closely match the neutron vibrational spectra. The main bands in the low-energy region (20–80 meV) are associated with the BH4 - librational modes. The features in the intermediate energy region (80–120 meV) are attributed to overtones and combination bands arising from the lower-energy modes. The features in the high-energy region (120–200 meV) correspond to the BH 4 - symmetric and asymmetric bending vibrations, of which four peaks located at 140, 142, 160, and 172 meV are especially intense. There are noticeable intensity distribution variations in the vibrational bands for different polymorphs. We can explain these differences using the spatial distribution of BH 4 - anions within various structures. An example of the possible identification of products after the hydrogenation of MgB 2 , using NVS measurements, is presented. Our results provide fundamental insights of benefit to researchers currently studying these promising hydrogen-storage materials.
The generation of accurate kinetic parameters such as mean generation time Λ and effective delayed neutron fraction β eff via Monte Carlo codes is established. Employing these in downstream deterministic codes warrants another step to ensure no additional error is introduced by the low-order transport operator when computing forward and adjoint fluxes for bilinear weighting of these parameters. Another complexity stems from applying superhomogenization (SPH) equivalence in non-fundamental mode approximations, where reference and low-order calculations rely on a 3D full core model. In these cases, SPH factors can optionally be computed for only part of the geometry while preserving reaction rates and K-effective, but the impact of such approximations on kinetics parameters has not been thoroughly studied. This paper aims at studying the preservation of bilinearly-weighted quantities in the Serpent–Griffin calculation procedure. Diffusion and transport evaluations of IPEN/MB-01, Godiva, and Flattop were carried out with the Griffin reactor physics code, testing available modeling options using Serpent-generated multigroup cross sections and equivalence data. Verifying Griffin against Serpent indicates sensitivities to multigroup energy grid selection and regional application of SPH equivalence, introducing significant errors; these were demonstrated to be reduced through the use of a transport method together with a finer energy grid.
Following our previous study of prototypical insertion reactions of energetically asymmetric type with the RPMD (Ring-Polymer Molecular Dynamics) method [Y. Li, Y. Suleimanov, and H. Guo, J. Phys. Chem. Lett. 5, 700 (2014)], we extend it to two other prototypical insertion reactions with much less exothermicity (near thermoneutral), namely, X + H{sub 2} → HX + H where X = C({sup 1}D), S({sup 1}D), in order to assess the accuracy of this method for calculating thermal rate coefficients for this class of reactions. For both chemical reactions, RPMD displays remarkable accuracy and agreement with the previous quantum dynamic results that make it encouraging for the future application of the RPMD to other barrier-less, complex-forming reactions involving polyatomic reactants with any exothermicity.
We report a rigorous full dimensional quantum dynamics algorithm, the multi-layer Lanczos method, for computing vibrational energies and dipole transition intensities of polyatomic molecules without any dynamics approximation. The multi-layer Lanczos method is developed by using a few advanced techniques including the guided spectral transform Lanczos method, multi-layer Lanczos iteration approach, recursive residue generation method, and dipole-wavefunction contraction. The quantum molecular Hamiltonian at the total angular momentum J = 0 is represented in a set of orthogonal polyspherical coordinates so that the large amplitude motions of vibrations are naturally described. In particular, the algorithm is general and problem-independent. An application is illustrated by calculating the infrared vibrational dipole transition spectrum of CH{sub 4} based on the ab initio T8 potential energy surface of Schwenke and Partridge [Spectrochimica Acta, Part A 57, 887 (2001)] and the low-order truncated ab initio dipole moment surfaces of Yurchenko et al. [J. Mol. Spectrosc. 291, 69 (2013)]. A comparison with experiments is made. The algorithm is also applicable for Raman polarizability active spectra.
We investigate atom-diatom reactive collisions, as a preliminary step, in order to assess the possibility of forming Rb{sub 2} molecules in their lowest triplet electronic state by cold collisions of rubidium atoms on the surface of helium nanodroplets. A simple model related to the well-known Rosen treatment of linear triatomic molecules [N. Rosen, J. Chem. Phys. 1, 319 (1933)] in relative coordinates is used, allowing to estimate reactive probabilities for different values of the total angular momentum. The best available full dimensional potential energy surface [Guillon et al., J. Chem. Phys. 136, 174307 (2012)] is employed through the calculations. Noticeable values of the probabilities in the ultracold regime, which numerically fulfill the Wigner threshold law, support the feasibility of the process. The rubidium dimer is mainly produced at high vibrational states, and the reactivity is more efficient for a bosonic helium partner than when the fermion species is considered.
The microscopic dielectric response is a key quantity for electronic materials such as organic semiconductors. Calculations of this response for molecular crystals are currently either expensive or rely on extreme simplifications such as multipole expansions which lack microscopic detail. We present an alternate approach using a microscopic analogue of the Clausius-Mossotti equation, which constructs the dielectric response of a crystal from an eigenvalue decomposition of the dielectric response of individual molecules. This method can potentially be used to examine the effects of defects, disorder, and surfaces on the dielectric properties of molecular solids.
The quantum Monte Carlo (QMC) technique is used to generate accurate energy benchmarks for methane-water clusters containing a single methane monomer and up to 20 water monomers. The benchmarks for each type of cluster are computed for a set of geometries drawn from molecular dynamics simulations. The accuracy of QMC is expected to be comparable with that of coupled-cluster calculations, and this is confirmed by comparisons for the CH{sub 4}-H{sub 2}O dimer. The benchmarks are used to assess the accuracy of the second-order Møller-Plesset (MP2) approximation close to the complete basis-set limit. A recently developed embedded many-body technique is shown to give an efficient procedure for computing basis-set converged MP2 energies for the large clusters. It is found that MP2 values for the methane binding energies and the cohesive energies of the water clusters without methane are in close agreement with the QMC benchmarks, but the agreement is aided by partial cancelation between 2-body and beyond-2-body errors of MP2. The embedding approach allows MP2 to be applied without loss of accuracy to the methane hydrate crystal, and it is shown that the resulting methane binding energy and the cohesive energy of the water lattice agree almost exactly with recently reported QMC values.
We present a formalism and an implementation for calculating spin-orbit couplings (SOCs) within the EOM-CCSD (equation-of-motion coupled-cluster with single and double substitutions) approach. The following variants of EOM-CCSD are considered: EOM-CCSD for excitation energies (EOM-EE-CCSD), EOM-CCSD with spin-flip (EOM-SF-CCSD), EOM-CCSD for ionization potentials (EOM-IP-CCSD) and electron attachment (EOM-EA-CCSD). We employ a perturbative approach in which the SOCs are computed as matrix elements of the respective part of the Breit-Pauli Hamiltonian using zeroth-order non-relativistic wave functions. We follow the expectation-value approach rather than the response-theory formulation for property calculations. Both the full two-electron treatment and the mean-field approximation (a partial account of the two-electron contributions) have been implemented and benchmarked using several small molecules containing elements up to the fourth row of the periodic table. The benchmark results show the excellent performance of the perturbative treatment and the mean-field approximation. When used with an appropriate basis set, the errors with respect to experiment are below 5% for the considered examples. The findings regarding basis-set requirements are in agreement with previous studies. The impact of different correlation treatment in zeroth-order wave functions is analyzed. Overall, the EOM-IP-CCSD, EOM-EA-CCSD, EOM-EE-CCSD, and EOM-SF-CCSD wave functions yield SOCs that agree well with each other (and with the experimental values when available). Using an EOM-CCSD approach that provides a more balanced description of the target states yields more accurate results.