Search NASA⌕ Search

SEARCH · Search NASA

Results for “Wave 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 1,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Theoretical calculations of rotationally inelastic collisions of He with NaK(A {sup 1}Σ{sup +}): Transfer of population, orientation, and alignment

We have performed extensive calculations to investigate thermal energy, rotationally inelastic collisions of NaK (A{sup 1}Σ{sup +}) with He. We determined a potential energy surface using a multi-reference configuration interaction wave function as implemented by the GAMESS electronic structure code, and we have performed coupled channel scattering calculations using the Arthurs and Dalgarno formalism. We also calculate the Grawert coefficients B{sub λ}(j, j′) for each j → j′ transition. These coefficients are used to determine the probability that orientation and alignment are preserved in collisions taking place in a cell environment. The calculations include all rotational levels with j or j′ between 0 and 50, and total (translational and rotational) energies in the range 0.0002–0.0025 a.u. (∼44–550 cm{sup −1}). The calculated cross sections for transitions with even values of Δj tend to be larger than those for transitions with odd Δj, in agreement with the recent experiments of Wolfe et al. (J. Chem. Phys. 134, 174301 (2011)). The calculations of the energy dependence of the cross sections and the calculations of the fraction of orientation and alignment preserved in collisions also exhibit distinctly different behaviors for odd and even values of Δj. The calculations also indicate that the average fraction of orientation or alignment preserved in a transition becomes larger as j increases. We interpret this behavior using the semiclassical model of Derouard, which also leads to a simple way of visualizing the distribution of the angles between the initial and final angular momentum vectors j and j′. Finally, we compare the exact quantum results for j → j′ transitions with results based on the simpler, energy sudden approximation. That approximation is shown to be quite accurate.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Quantum theory of atoms in molecules/charge-charge flux-dipole flux models for fundamental vibrational intensity changes on H-bond formation of water and hydrogen fluoride

The Quantum Theory of Atoms In Molecules/Charge-Charge Flux-Dipole Flux (QTAIM/CCFDF) model has been used to investigate the electronic structure variations associated with intensity changes on dimerization for the vibrations of the water and hydrogen fluoride dimers as well as in the water-hydrogen fluoride complex. QCISD/cc-pVTZ wave functions applied in the QTAIM/CCFDF model accurately provide the fundamental band intensities of water and its dimer predicting symmetric and antisymmetric stretching intensity increases for the donor unit of 159 and 47 km mol{sup −1} on H-bond formation compared with the experimental values of 141 and 53 km mol{sup −1}. The symmetric stretching of the proton donor water in the dimer has intensity contributions parallel and perpendicular to its C{sub 2v} axis. The largest calculated increase of 107 km mol{sup −1} is perpendicular to this axis and owes to equilibrium atomic charge displacements on vibration. Charge flux decreases occurring parallel and perpendicular to this axis result in 42 and 40 km mol{sup −1} total intensity increases for the symmetric and antisymmetric stretches, respectively. These decreases in charge flux result in intensity enhancements because of the interaction contributions to the intensities between charge flux and the other quantities. Even though dipole flux contributions are much smaller than the charge and charge flux ones in both monomer and dimer water they are important for calculating the total intensity values for their stretching vibrations since the charge-charge flux interaction term cancels the charge and charge flux contributions. The QTAIM/CCFDF hydrogen-bonded stretching intensity strengthening of 321 km mol{sup −1} on HF dimerization and 592 km mol{sup −1} on HF:H{sub 2}O complexation can essentially be explained by charge, charge flux and their interaction cross term. Atomic contributions to the intensities are also calculated. The bridge hydrogen atomic contributions alone explain 145, 237, and 574 km mol{sup −1} of the H-bond stretching intensity enhancements for the water and HF dimers and their heterodimer compared with total increments of 149, 321, and 592 km mol{sup −1}, respectively.

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

General Description of Fission Observables: GEF Model Code

The GEF (“GEneral description of Fission observables”) model code is documented. It describes the observables for spontaneous fission, neutron-induced fission and, more generally, for fission of a compound nucleus from any other entrance channel, with given excitation energy and angular momentum. The GEF model is applicable for a wide range of isotopes from Z = 80 to Z = 112 and beyond, up to excitation energies of about 100 MeV. The results of the GEF model are compared with fission barriers, fission probabilities, fission-fragment mass- and nuclide distributions, isomeric ratios, total kinetic energies, and prompt-neutron and prompt-gamma yields and energy spectra from neutron-induced and spontaneous fission. Derived properties of delayed neutrons and decay heat are also considered. The GEF model is based on a general approach to nuclear fission that explains a great part of the complex appearance of fission observables on the basis of fundamental laws of physics and general properties of microscopic systems and mathematical objects. The topographic theorem is used to estimate the fission-barrier heights from theoretical macroscopic saddle-point and ground-state masses and experimental ground-state masses. Motivated by the theoretically predicted early localisation of nucleonic wave functions in a necked-in shape, the properties of the relevant fragment shells are extracted. These are used to determine the depths and the widths of the fission valleys corresponding to the different fission channels and to describe the fission-fragment distributions and deformations at scission by a statistical approach. A modified composite nuclear-level-density formula is proposed. It respects some features in the superfluid regime that are in accordance with new experimental findings and with theoretical expectations. These are a constant-temperature behaviour that is consistent with a considerably increased heat capacity and an increased pairing condensation energy that is consistent with the collective enhancement of the level density. The exchange of excitation energy and nucleons between the nascent fragments on the way from saddle to scission is estimated according to statistical mechanics. As a result, excitation energy and unpaired nucleons are predominantly transferred to the heavy fragment in the superfluid regime. This description reproduces some rather peculiar observed features of the prompt-neutron multiplicities and of the even-odd effect in fission-fragment Z distributions. For completeness, some conventional descriptions are used for calculating pre-equilibrium emission, fission probabilities and statistical emission of neutrons and gamma radiation from the excited fragments. Preference is given to simple models that can also be applied to exotic nuclei compared to more sophisticated models that need precise empirical input of nuclear properties, e.g. spectroscopic information. The approach reveals a high degree of regularity and provides a considerable insight into the physics of the fission process. Fission observables can be calculated with a precision that complies with the needs for applications in nuclear technology without specific adjustments to measured data of individual systems. The GEF executable runs out of the box with no need for entering any empirical data. This unique feature is of valuable importance, because the number of systems and energies of potential significance for fundamental and applied science will never be possible to be measured. The relevance of the approach for examining the consistency of experimental results and for evaluating nuclear data is demonstrated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2014

Intervalence charge transfer luminescence: Interplay between anomalous and 5d − 4f emissions in Yb-doped fluorite-type crystals

In this paper, we report the existence of intervalence charge transfer (IVCT) luminescence in Yb-doped fluorite-type crystals associated with Yb{sup 2+}–Yb{sup 3+} mixed valence pairs. By means of embedded cluster, wave function theory ab initio calculations, we show that the widely studied, very broad band, anomalous emission of Yb{sup 2+}-doped CaF{sub 2} and SrF{sub 2}, usually associated with impurity-trapped excitons, is, rather, an IVCT luminescence associated with Yb{sup 2+}–Yb{sup 3+} mixed valence pairs. The IVCT luminescence is very efficiently excited by a two-photon upconversion mechanism where each photon provokes the same strong 4f{sup 14}–1A{sub 1g}→ 4f{sup 13}({sup 2}F{sub 7/2})5de{sub g}–1T{sub 1u} absorption in the Yb{sup 2+} part of the pair: the first one, from the pair ground state; the second one, from an excited state of the pair whose Yb{sup 3+} moiety is in the higher 4f{sup 13}({sup 2}F{sub 5/2}) multiplet. The Yb{sup 2+}–Yb{sup 3+} → Yb{sup 3+}–Yb{sup 2+} IVCT emission consists of an Yb{sup 2+} 5de{sub g} → Yb{sup 3+} 4f{sub 7/2} charge transfer accompanied by a 4f{sub 7/2} → 4f{sub 5/2} deexcitation within the Yb{sup 2+} 4f{sup 13} subshell: [{sup 2}F{sub 5/2}5de{sub g},{sup 2}F{sub 7/2}] → [{sup 2}F{sub 7/2},4f{sup 14}]. The IVCT vertical transition leaves the oxidized and reduced moieties of the pair after electron transfer very far from their equilibrium structures; this explains the unexpectedly large band width of the emission band and its low peak energy, because the large reorganization energies are subtracted from the normal emission. The IVCT energy diagrams resulting from the quantum mechanical calculations explain the different luminescent properties of Yb-doped CaF{sub 2}, SrF{sub 2}, BaF{sub 2}, and SrCl{sub 2}: the presence of IVCT luminescence in Yb-doped CaF{sub 2} and SrF{sub 2}; its coexistence with regular 5d-4f emission in SrF{sub 2}; its absence in BaF{sub 2} and SrCl{sub 2}; the quenching of all emissions in BaF{sub 2}; and the presence of additional 5d–4f emissions in SrCl{sub 2} which are absent in SrF{sub 2}. They also allow to interpret and reproduce recent experiments on transient photoluminescence enhancement in Yb{sup 2+}-doped CaF{sub 2} and SrF{sub 2}, the appearance of Yb{sup 2+} 4f–5d absorption bands in the excitation spectra of the IR Yb{sup 3+} emission in partly reduced CaF{sub 2}:Yb{sup 3+} samples, and to identify the broadband observed in the excitation spectrum of the so far called anomalous emission of SrF{sub 2}:Yb{sup 2+} as an IVCT absorption, which corresponds to an Yb{sup 2+} 4f{sub 5/2} → Yb{sup 3+} 4f{sub 7/2} electron transfer.

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

Science Requirements and Conceptual Design for a Polarized Medium Energy Electron-Ion Collider at Jlab

Researchers have envisioned an electron-ion collider with ion species up to heavy ions, high polarization of electrons and light ions, and a well-matched center-of-mass energy range as an ideal gluon microscope to explore new frontiers of nuclear science. In its most recent Long Range Plan, the Nuclear Science Advisory Committee (NSAC) of the US Department of Energy and the National Science Foundation endorsed such a collider in the form of a 'half-recommendation.' As a response to this science need, Jefferson Lab and its user community have been engaged in feasibility studies of a medium energy polarized electron-ion collider (MEIC), cost-effectively utilizing Jefferson Lab's already existing Continuous Electron Beam Accelerator Facility (CEBAF). In close collaboration, this community of nuclear physicists and accelerator scientists has rigorously explored the science case and design concept for this envisioned grand instrument of science. An electron-ion collider embodies the vision of reaching the next frontier in Quantum Chromodynamics - understanding the behavior of hadrons as complex bound states of quarks and gluons. Whereas the 12 GeV Upgrade of CEBAF will map the valence-quark components of the nucleon and nuclear wave functions in detail, an electron-ion collider will determine the largely unknown role sea quarks play and for the first time study the glue that binds all atomic nuclei. The MEIC will allow nuclear scientists to map the spin and spatial structure of quarks and gluons in nucleons, to discover the collective effects of gluons in nuclei, and to understand the emergence of hadrons from quarks and gluons. The proposed electron-ion collider at Jefferson Lab will collide a highly polarized electron beam originating from the CEBAF recirculating superconducting radiofrequency (SRF) linear accelerator (linac) with highly polarized light-ion beams or unpolarized light- to heavy-ion beams from a new ion accelerator and storage complex. Since the very beginning, the design studies at Jefferson Lab have focused on achieving high collider performance, particularly ultrahigh luminosities up to 10{sup 34} cm{sup -2}s{sup -1} per detector with large acceptance, while maintaining high polarization for both the electron and light-ion beams. These are the two key performance requirements of a future electron-ion collider facility as articulated by the NSAC Long Range Plan. In MEIC, a new ion complex is designed specifically to deliver ion beams that match the high bunch repetition and highly polarized electron beam from CEBAF. During the last two years, both development of the science case and optimization of the machine design point toward a medium-energy electron-ion collider as the topmost goal for Jefferson Lab. The MEIC, with relatively compact collider rings, can deliver a luminosity above 10{sup 34} cm{sup -2}s{sup -1} at a center-of-mass energy up to 65 GeV. It offers an electron energy up to 11 GeV, a proton energy up to 100 GeV, and corresponding energies per nucleon for heavy ions with the same magnetic rigidity. This design choice balances the scope of the science program, collider capabilities, accelerator technology innovation, and total project cost. An energy upgrade could be implemented in the future by adding two large collider rings housed in another large tunnel to push the center-of-mass energy up to or exceeding 140 GeV. After careful consideration of an alternative electron energy recovery linac on ion storage ring approach, a ring-ring collider scenario at high bunch repetition frequency was found to offer fully competitive performance while eliminating the uncertainties of challenging R&D on ampere-class polarized electron sources and many-pass energy-recovery linacs (ERLs). The essential new elements of an MEIC facility at Jefferson Lab are an electron storage ring and an entirely new, modern ion acceleration and storage complex. For the high-current electron collider ring, the upgraded 12 GeV CEBAF SRF linac will serve as a full-energy injector, and, if needed, provide top-off refilling. The CEBAF fixed-target nuclear physics program can be simultaneously operated since the filling time of the electron ring is very short. The ion complex for MEIC consists of sources for polarized light ions and unpolarized light to heavy ions, an SRF ion linac with proton energy up to 280 MeV, a 3 GeV prebooster synchrotron, a large booster synchrotron for proton energy up to 20 GeV, and a medium-energy collider ring with energy up to 100 GeV. The ion complex can accelerate other species of ions with corresponding energies at each accelerating stage. There are three collision points planned for MEIC. Two of them are for collisions with medium-energy ions; the third is for low energy ion beams stored in a dedicated low-energy compact storage ring, as a possible follow-on project.

43 PARTICLE ACCELERATORS↗
DOE OSTI2016

Rigorous Ab Initio Quantum Embedding for Quantum Chemistry Using Green’s Function Theory: Screened Interaction, Nonlocal Self-Energy Relaxation, Orbital Basis, and Chemical Accuracy

We present a detailed discussion of the self-energy embedding theory (SEET), which is a quantum embedding scheme allowing us to describe a chosen subsystem very accurately while keeping the description of the environment at a lower level. We apply SEET to molecular examples where our chosen subsystem is made out of a set of strongly correlated orbitals while the weakly correlated orbitals constitute an environment. Consequently, a highly accurate method is used to calculate the self-energy for the system, while a lower-level method is employed to find the self-energy for the environment. Such a self-energy separation is very general, and to make the SEET procedure applicable to multiple systems, a detailed and practical procedure for the evaluation of the system and environment self-energy is necessary. We list all of the intricacies for one of the possible procedures while focusing our discussion on many practical implementation aspects such as the choice of best orbital basis, impurity solver, and many steps necessary to reach chemical accuracy. Finally, on a set of carefully chosen molecular examples, we demonstrate that SEET, which is a controlled, systematically improvable Green’s function method, can be as accurate as established wave function quantum chemistry methods.

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

Is the Inversion of Phosphorus Trihalides (PF 3 , PCl 3 , PBr 3 , and PI 3 ) a Diradical Process?

This work explores possible reaction paths for the inversion of a series of trigonal pyramidal phosphorus trihalides, PF 3 , PCl 3 , PBr 3 , and PI 3 , and it especially addresses the question of whether and when the bonding of the lowest-energy species along the inversion paths should be described as a hyper open-shell diradical. The various paths for inversion are calculated using a single-reference method within the framework of Kohn-Sham density functional theory and also with multi-reference wave function methods. Our calculated results using both kinds of methods show that, for all the halogens studied (F, Cl, Br, and I), the lowestenergy singlet path for the inversion occurs by the formation of a C 2v transition structure rather than a D 3 h transition structure. This geometrical preference agrees with what has been inferred previously based on closed-shell singlet calculations. But in the present study, we examined not only closed-shell singlet transition states but also open-shell singlet states and triplet states for calculating stationary points and inversion paths, and for some of the phosphorus trihalides, we found that paths involving open-shell configurations are lower in energy than those restricted to closed-shell configurations. We analyzed the changes along the paths in terms of hybridization and orientation of the frontier orbitals and in terms of locally avoided crossings, and the extent of diradical character was quantified by calculating the effective number of unpaired electrons. Even for the singlet inversion path that goes via a D 3 h structure, the barrier for PF 3 , PCl 3 , and PBr 3 is higher for a closed-shell singlet spin state than for the open-shell singlet configuration. Furthermore, the energy of the triplet D 3 h structure is below even the open-shell D 3 h singlet for PCl 3 , PBr 3 , and PI 3 . This necessitates rethinking the role of open-shell states in nominally closedshell processes.

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

Atomistic Defect Makes a Phase Plate for the Generation and High-Angular Splitting of Electron Vortex Beams

Topological defects in solid-state materials by breaking the translational symmetry offer emerging properties that are not present in their parental phases. For example, edge dislocations—the 2π phase-winding topological defects—in antiferromagnetic NiO crystals can exhibit ferromagnetic behaviors. We study how these defects could give rise to topological orders when they interact with a high-energy electron beam. To probe this interaction, we formed a coherent electron nanobeam in a scanning transmission electron microscope and recorded the far-field transmitted patterns as the beam steps through the edge dislocation core in [001] NiO. Surprisingly, we found the amplitude patterns of the <020> Bragg disks evolve in a similar manner to the evolution of an annular solar eclipse. Using the ptychographic technique, we recovered the missing phase information in the diffraction plane and revealed the topological phase vortices in the diffracted beams. Through atomic topological defects, the wave function of electrons can be converted from plane wave to electron vortex. Technologically, this approach provides a feasible route for the fabrication of phase plates that can generate electron vortex beams with an angular separation that is three orders of magnitude larger than what traditional nanofabrication technology can offer. This advance will enable the collection of magnetic circular dichroism spectra with high spatial resolution and high efficiency, boosting the understanding of the relationship between symmetry breaking and magnetic property of individual topological defect at the atomic scale.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Quantum Monte Carlo calculation of the binding energy of the beryllium dimer

The accurate calculation of the binding energy of the beryllium dimer is a challenging theoretical problem. In this study, the binding energy of Be{sub 2} is calculated using the diffusion Monte Carlo (DMC) method, using single Slater determinant and multiconfigurational trial functions. DMC calculations using single-determinant trial wave functions of orbitals obtained from density functional theory calculations overestimate the binding energy, while DMC calculations using Hartree-Fock or CAS(4,8), complete active space trial functions significantly underestimate the binding energy. In order to obtain an accurate value of the binding energy of Be{sub 2} from DMC calculations, it is necessary to employ trial functions that include excitations outside the valence space. Our best estimate DMC result for the binding energy of Be{sub 2}, obtained by using configuration interaction trial functions and extrapolating in the threshold for the configurations retained in the trial function, is 908 cm{sup −1}, only slightly below the 935 cm{sup −1} value derived from experiment.

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

Roadmap on methods and software for electronic structure based simulations in chemistry and materials

This Roadmap article provides a succinct, comprehensive overview of the state of electronic structure methods and software for molecular and materials simulations. Seventeen distinct sections collect insights by 51 leading scientists in the field. Each contribution addresses the status of a particular area, as well as current challenges and anticipated future advances, with a particular eye towards software related aspects and providing key references for further reading. Foundational sections cover density functional theory and its implementation in real-world simulation frameworks, Green's function based many-body perturbation theory, wave-function based and stochastic electronic structure approaches, relativistic effects and semiempirical electronic structure theory approaches. Subsequent sections cover nuclear quantum effects, real-time propagation of the electronic structure, challenges for computational spectroscopy simulations, and exploration of complex potential energy surfaces. The final sections summarize practical aspects, including computational workflows for complex simulation tasks, the impact of current and future high-performance computing architectures, software engineering practices, education and training to maintain and broaden the community, as well as the status of and needs for electronic structure based modeling from the vantage point of industry environments. Overall, the field of electronic structure software and method development continues to unlock immense opportunities for future scientific discovery, based on the growing ability of computations to reveal complex phenomena, processes and properties that are determined by the make-up of matter at the atomic scale, with high precision.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Quantum Monte Carlo calculations of light nuclei with local chiral two- and three-nucleon interactions

Local chiral effective field theory interactions have recently been developed and used in the context of quantum Monte Carlo few- and many-body methods for nuclear physics. In this paper, we go over detailed features of local chiral nucleon-nucleon interactions and examine their effect on properties of the deuteron, paying special attention to the perturbativeness of the expansion. We then turn to three-nucleon interactions, focusing on operator ambiguities and their interplay with regulator effects. We then discuss the nuclear Green's function Monte Carlo method, going over both wave-function correlations and approximations for the two- and three-body propagators. Finally, following this, we present a range of results on light nuclei: Binding energies and distribution functions are contrasted and compared, starting from several different microscopic interactions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Detection of magnetic circular dichroism in amorphous materials utilizing a single-crystalline overlayer

Physicists are fascinated with topological defects in solid-state materials, because by breaking the translational symmetry they offer emerging properties that are not present in their parental phases. For example, edge dislocations—the 2π phase-winding topological defects—in antiferromagnetic NiO crystals can exhibit ferromagnetic behaviors. Herein, we study how these defects could give rise to exotic topological orders when they interact with a high energy electron beam. To probe this interaction, we formed a coherent electron nanobeam in a scanning transmission electron microscope and recorded the far-field transmitted patterns as the beam steps through the edge dislocation core in [001] NiO. Surprisingly, we found the amplitude patterns of the <020> Bragg disks evolve in a similar manner to the evolution of an annular solar eclipse. Using the ptychographic technique, we recovered the missing phase information in the diffraction plane and revealed the topological phase vortices in the diffracted beams. Through atomic topological defects, the wave function of electrons can be converted from plane wave to electron vortex. This approach provides a new perspective for boosting the collection efficiency of magnetic circular dichroism spectra with high spatial resolution and understanding the relationship between symmetry breaking and exotic property of individual topological defect at atomic level.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Effects of spin-orbit coupling and very large supercells on the description of acceptors in CdTe

Predicting accurate shallow donor and acceptor levels in semiconductors has been quite challenging using periodic boundary conditions as implemented in current density functional theory codes. The reason is that the wave functions associated with the shallow centers are quite extended and are not fully contained in the supercell (typically with a few hundred atoms) so impurities in the periodically repeated image cells interact with each other. Errors of ~0.1 eV are expected, and these are of the same order of magnitude as the ionization energies themselves. In the case of acceptors in CdTe, this problem is exacerbated by the strong spin-orbit coupling that split the Te-related states at the top of the valence band, making the calculations at least 8 times more expensive. CdTe is an important solar-cell material with record high efficiency of 22%. One of the main limiting factors to increasing the efficiency towards the theoretical limit of ~30% is the often reported very low hole concentration. Pushing to the limit of computational capability by using very large supercells with spin-orbit coupling we report the results of hybrid functional calculations of group-V acceptors in CdTe. We show that extrapolation to the dilute limit leads to an interpretation of the experimental data that is qualitatively different from previous DFT and hybrid functional calculation reports. We find that the group-V impurities indeed behave as shallow acceptors and that the corresponding compensating AX-centers are unstable and do not limit p-type doping. We address the differences between our results and previous theoretical predictions and show that our calculated ionization energies predict hole concentrations that are in excellent agreement with recent temperature-dependent Hall measurements on high-quality single-crystal samples.

14 SOLAR ENERGY↗
DOE OSTI2025

Overset-Grid Method with Smooth Orbital Partitioning for Molecular Scattering Calculations

To solve molecular photoionization and electron scattering problems, we use an overset-grid representation of electronic continuum functions, which has an extended central spherical grid that overlaps small spherical grids (subgrids) centered on each atom of a polyatomic molecule. Here, in this work, we present an improved algorithm that smoothly partitions the total wave function between the central grid and the atomic subgrids. The smooth partitioning allows one to use approximately one-fourth the number of partial waves on the central grid compared to our previous implementation with switching functions. The resulting numerical method for treating electron scattering and photoionization of polyatomic molecules combines the accuracy and flexibility of pure numerical grid representations with the rapid convergence of hybrid combinations of atom-centered basis-set expansions and grid methods. The overset-grid representation is implemented using the complex Kohn variational principle for scattering and photoionization amplitudes. The faster convergence with respect to the number of central grid partial waves is demonstrated and accuracy is verified by comparisons with the previous implementation and with far more computationally demanding single-center numerical expansions in electron-molecule scattering and photoionization calculations on the neon dimer (Ne 2 ) system, carbon tetrafluoride (CF 4 ) molecule, and the pyridine (C 5 H 5 N) molecule in the static-exchange approximation.

Molecules↗
DOE OSTI2018

Configurational forces in electronic structure calculations using Kohn-Sham density functional theory

Here we derive the expressions for configurational forces in Kohn-Sham density functional theory, which correspond to the generalized variational force computed as the derivative of the Kohn-Sham energy functional with respect to the position of a material point x. These configurational forces that result from the inner variations of the Kohn-Sham energy functional provide a unified framework to compute atomic forces as well as stress tensor for geometry optimization. Importantly, owing to the variational nature of the formulation, these configurational forces inherently account for the Pulay corrections. The formulation presented in this work treats both pseudopotential and all-electron calculations in a single framework, and employs a local variational real-space formulation of Kohn-Sham density functional theory (DFT) expressed in terms of the nonorthogonal wave functions that is amenable to reduced-order scaling techniques. We demonstrate the accuracy and performance of the proposed configurational force approach on benchmark all-electron and pseudopotential calculations conducted using higher-order finite-element discretization. To this end, we examine the rates of convergence of the finite-element discretization in the computed forces and stresses for various materials systems, and, further, verify the accuracy from finite differencing the energy. Wherever applicable, we also compare the forces and stresses with those obtained from Kohn-Sham DFT calculations employing plane-wave basis (pseudopotential calculations) and Gaussian basis (all-electron calculations). Finally, we verify the accuracy of the forces on large materials systems involving a metallic aluminum nanocluster containing 666 atoms and an alkane chain containing 902 atoms, where the Kohn-Sham electronic ground state is computed using a reduced-order scaling subspace projection technique

36 MATERIALS SCIENCE↗
NASA NTRS1981

Functional methods for waves in random media

Some basic ideas in functional methods for waves in random media are illustrated through a simple random differential equation. These methods are then generalized to solve certain random parabolic equations via an exponential representation given by the Feynman-Kac formula. It is shown that these functional methods are applicable to a number of problems in random wave propagation. They include the forward-scattering approximation in Gaussian white-noise media; the solution of the optical beam propagation problem by a phase-integral method; the high-frequency scattering by bounded random media, and a derivation of approximate moment equations from the functional integral representation.

Chow, P. L.↗
NASA NTRS1981

Functional methods for waves in random media

Some basic ideas in functional methods for waves in random media are illustrated through a simple random differential equation. These methods are then generalized to solve certain random parabolic equations via an exponential representation given by the Feynman-Kac formula. It is shown that these functional methods are applicable to a number of problems in random wave propagation. They include the forward-scattering approximation in Gaussian white-noise media; the solution of the optical beam propagation problem by a phase-integral method; the high-frequency scattering by bounded random media; and a derivation of approximate moment equations from the functional integral representation.

Chow, P. L.↗
DOE OSTI2016

Dielectric function and plasmons in graphene: A self-consistent-field calculation within a Markovian master equation formalism

We introduce a method for calculating the dielectric function of nanostructures with an arbitrary band dispersion and Bloch wave functions. The linear response of a dissipative electronic system to an external electromagnetic field is calculated by a self-consistent-field approach within a Markovian master equation formalism (SCF-MMEF) coupled with full-wave electromagnetic equations. The SCF-MMEF accurately accounts for several concurrent scattering mechanisms. The method captures interband electron-hole-pair generation, as well as the interband and intraband electron scattering with phonons and impurities. We employ the SCF-MMEF to calculate the dielectric function, complex conductivity, and loss function for supported graphene. From the loss-function maximum, we obtain plasmon dispersion and propagation length for different substrate types [nonpolar diamondlike carbon (DLC) and polar SiO 2 and hBN], impurity densities, carrier densities, and temperatures. Plasmons on the two polar substrates are suppressed below the highest surface phonon energy, while the spectrum is broad on the nonpolar DLC. Plasmon propagation lengths are comparable on polar and nonpolar substrates and are on the order of tens of nanometers, considerably shorter than previously reported. As a result, they improve with fewer impurities, at lower temperatures, and at higher carrier densities.

77 NANOSCIENCE AND NANOTECHNOLOGY↗