Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SLATER METHOD”

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 91 records · Page 5

Excited electronic states from a variational approach based on symmetry-projected Hartree–Fock configurations

Recent work from our research group has demonstrated that symmetry-projected Hartree–Fock (HF) methods provide a compact representation of molecular ground state wavefunctions based on a superposition of non-orthogonal Slater determinants. The symmetry-projected ansatz can account for static correlations in a computationally efficient way. Here we present a variational extension of this methodology applicable to excited states of the same symmetry as the ground state. Benchmark calculations on the C{sub 2} dimer with a modest basis set, which allows comparison with full configuration interaction results, indicate that this extension provides a high quality description of the low-lying spectrum for the entire dissociation profile. We apply the same methodology to obtain the full low-lying vertical excitation spectrum of formaldehyde, in good agreement with available theoretical and experimental data, as well as to a challenging model C{sub 2v} insertion pathway for BeH{sub 2}. The variational excited state methodology developed in this work has two remarkable traits: it is fully black-box and will be applicable to fairly large systems thanks to its mean-field computational cost.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spectroradiometric considerations for advanced land observing systems

Research aimed at improving the inflight absolute radiometric calibration of advanced land observing systems was initiated. Emphasis was on the satellite sensor calibration program at White Sands. Topics addressed include: absolute radiometric calibration of advanced remote sensing; atmospheric effects on reflected radiation; inflight radiometric calibration; field radiometric methods for reflectance and atmospheric measurement; and calibration of field relectance radiometers.

Slater, P. N.↗

Variations in in-flight absolute radiometric calibration

Variations in the in-flight absolute radiometric calibration of the Coastal Zone Color Scanner and the Thematic Mapper (TM) are reviewed. At short wavelengths, the sensors show a gradual reduction in response, while in the mid-IR the TM shows oscillatory variations. One set of measurements made at White Sands, New Mexico shows anomalous results in TM bands 2 and 4. The results of a reflectance-based and a radiance-based calibration method at White Sands are described. An analysis of the radiance-based method shows the value of such measurements from helicopter altitudes for calibration.

Slater, Philip N.↗

Factorization in large-scale many-body calculations

One approach for solving interacting many-fermion systems is the configuration-interaction method, also sometimes called the interacting shell model, where one finds eigenvalues of the Hamiltonian in a many-body basis of Slater determinants (antisymmetrized products of single-particle wavefunctions). The resulting Hamiltonian matrix is typically very sparse, but for large systems the nonzero matrix elements can nonetheless require terabytes or more of storage. An alternate algorithm, applicable to a broad class of systems with symmetry, in our case rotational invariance, is to exactly factorize both the basis and the interaction using additive/multiplicative quantum numbers; such an algorithm recreates the many-body matrix elements on the fly and can reduce the storage requirements by an order of magnitude or more. Here, we discuss factorization in general and introduce a novel, generalized factorization method, essentially a ‘double-factorization’ which speeds up basis generation and set-up of required arrays. Although we emphasize techniques, we also place factorization in the context of a specific (unpublished) configuration-interaction code, BIGSTICK, which runs both on serial and parallel machines, and discuss the savings in memory due to factorization.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Extending accurate spectral induced polarization measurements into the kHz range: modelling and removal of errors from interactions between the parasitic capacitive coupling and the sample holder

SUMMARY Spectral induced polarization (SIP), describing the measurement of the frequency domain electrical impedance magnitude and phase of porous materials, has been widely used to characterize subsurface hydrological/biogeochemical properties and processes. SIP data collected at frequencies higher than 100 Hz are expected to describe the polarization of small particles providing insights into the physicochemical properties of clays, nanoparticles and microorganisms. However, the phase measurements at these high frequencies are often contaminated by errors due to the parasitic capacitive coupling of the SIP instrument, especially for lower conductivity samples. We developed a model showing the measured phase is the sum of the true sample phase and an error term ωCinZx, where ω is the angular frequency; Cin is the instrument input capacitance and Zx is a measurable impedance function related to the sample holder properties and the reference resistor. Based on this model, a new phase correction method is proposed that results in highly accurate SIP data up to 20 kHz as well as the determination of Cin. We tested the correction method using electric circuits, NaCl fluids and three unconsolidated samples (sand, sand-clay and sand-pyrite mixtures). The corrected phase for the circuit and NaCl fluid experiments showed excellent agreement with the theoretical phase response across the studied frequency range (errors <1 mrad). For unconsolidated samples, removal of errors results in phase spectra more consistent with expected polarization mechanisms, as based on phase peaks recorded for small pyrite and clay particles at high frequencies. These phase peaks could not be identified in the uncorrected data. Our approach can substantially enhance the value of the SIP method for the characterization of fine-grained sediments and rocks.

58 GEOSCIENCES↗

SeaWiFS calibration and algorithm validation

A substantial amount of field research has been completed during this period (1 October 1994 - 31 October 1995). The applicability of different in-flight calibration methods was studied. In addition, various field instruments and new equipment were tested to appraise their usefulness during such calibration activities. The purposes and (expected) results of these field experiments are detailed in separate sections of this report. One section contains a synopsis of calibrations conducted over both land and water targets. A second section contains the results of an experiment conducted to determine the spatial uniformity of surface reflectance at Lake Tahoe, the chosen water calibration site. Advances were also made in the area of radiative transfer software development. An existing exact radiative transfer code, Successive Orders (SO), was modified to provide outputs of the radiance distribution at arbitrary altitudes above a selected calibration target. These modifications were needed to utilize the SO software in radiance-based sensor calibrations. Serious thought was given to the general applicability of various calibration techniques to calibrations conducted over water surfaces. In particular, the practicality of the radiance-based and reflectance-based methods of in-flight radiometric calibration was evaluated.

Parada, Robert J., Jr.↗

Tight-Binding-Based Method for Simulations of Complex Materials with Accompanying Database

The NRL-Tight-Binding (NRL-TB) method has been applied successfully to many monatomic, binary, and a few ternary systems. The method consists of a fit to Density Functional Theory (DFT) results and has transferability beyond the fitted DFT database including performing tight-binding molecular dynamics simulations. We are working on a new website with a detailed guide of how to use the method in performing the fit to DFT results (Fitting Code) and obtaining non-orthogonal TB parameters expressed in a polynomial form as environment-dependent Slater-Koster parameters (on-site terms), and neighbor distance-dependent Hamiltonian and Overlap terms. These parameters are then used in a Static Code to verify or make predictions for large systems often beyond the computational capabilities of standard DFT codes. Our starting point was a website developed at NRL some years ago which NRL no longer supports. We have now been augmenting that website with new examples of expanding the original capabilities of the method such as incorporating in the codes f-electron states. Our website (esd.cos.gmu.edu/tb) contains two separate branches: One named NRL-TB with the full documentation of the method including source codes and examples. The other branch, named ESD-dB, is a library of electronic structure calculations for all the monatomic elements (a subset of the elements also found in a recent book by the author of this proposal. See reference 22) and an approximately 700 hundred compounds in specific crystal structures. We discuss below the progress made so far in improving and expanding the capabilities of our TB method and the plans for further extensions of our approach.

36 MATERIALS SCIENCE↗

Ab initio study of Fe{sub 2}MnZ (Al, Si, Ge) Heusler alloy using GGA approximation

Density functional theory based on FP-LAPW method used to investigate the electronic structure of Fe{sub 2}MnZ, shows that the total spin magnetic moment shows a trend consistent with the Slater–Pauling curve. The Fe and Mn magnetic moment depend on choice of Z element although the magnetic moment of Z element is negative and less than 0.1 μ{sub B}. Spin polarization calculations evidence 100% spin polarization for Fe{sub 2}MnSi. Fe{sub 2}MnAl and Fe{sub 2}MnGe show metallic behavior with 93%, 98% spin polarization.

36 MATERIALS SCIENCE↗

Nonadiabatic Molecular Dynamics Simulations Based on Time-Dependent Density Functional Tight-Binding Method

Nonadiabatic excited-state molecular dynamics underpin many photophysical and photochemical phenomena, such as exciton dynamics, charge separation and transport. In this work, we present an efficient nonadiabatic molecular dynamic (NAMD) simulation method based on time-dependent density functional tight-binding (TDDFTB) theory. Specifically, the adiabatic electronic structure, an essential NAMD input, is described at the TDDFTB level. Here, the nonadiabatic effects originating from the coupled motions of electrons and nuclei are treated by the trajectory surface hopping algorithm. To improve the computational efficiency, nonadiabatic couplings between excited states within the TDDFTB method are derived and implemented using an analytical approach. Further, the time-dependent nonadiabatic coupling scalars are calculated based on the overlap between molecular orbitals rather than the Slater determinants to speed up the simulations. In addition, the electronic decoherence scheme and a state reassigned unavoided crossings algorithm, which has been implemented in the NEXMD software, are used to improve the accuracy of the simulated dynamics and handle trivial unavoided crossings. Finally, the photoinduced nonadiabatic dynamics of a benzene molecule are simulated to demonstrate our implementation. The results for excited state NAMD simulations of benzene molecule based on TDDFTB method compare well that obtained with numerically expensive time-dependent density functional theory. The proposed methodology provides an attractive theoretical simulation tool for predicting the photophysical and photochemical properties of complex materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Calculation of photodetachment cross sections and photoelectron angular distributions of negative ions using density functional theory

Recently, the development of photoelectron velocity map imaging makes it much easier to obtain the photoelectron angular distributions (PADs) experimentally. However, explanations of PADs are only qualitative in most cases, and very limited works have been reported on how to calculate PAD of anions. In the present work, we report a method using the density-functional-theory Kohn-Sham orbitals to calculate the photodetachment cross sections and the anisotropy parameter β. The spherical average over all random molecular orientation is calculated analytically. A program which can handle both the Gaussian type orbital and the Slater type orbital has been coded. The testing calculations on Li{sup −}, C{sup −}, O{sup −}, F{sup −}, CH{sup −}, OH{sup −}, NH{sub 2}{sup −}, O{sub 2}{sup −}, and S{sub 2}{sup −} show that our method is an efficient way to calculate the photodetachment cross section and anisotropy parameter β for anions, thus promising for large systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Three methods for the absolute calibration of the NOAA AVHRR sensors in-flight

Three techniques for the absolute calibration of two reflective channels of the NOAA AVHRR sensors are described and compared. The technique 1 is based on field measurements and refers to another calibrated satellite sensor; technique 2 utilizes the reflectance-based calibration method; and technique 3 estimates the calibration by reference to another satellite sensor. NOAA-9 and NOAA-10 AVHRR images are analyzed. It is observed that the responsivity of the images degrade compared to prelaunch calibration. The effect of degradation in responsivity of the AVHRR sensors on vegetation indices are studied. It is noted that the uncertainties in the calibration methods is about 7-10 percent making them useful for many applications.

Teillet, P. M.↗

SU-E-T-239: Design and Evaluation of a Nozzle Shielding System in Passively Scattered Proton Therapy

Purpose: To design, implement and evaluate a shielding system that will reduce out-of-field dose experienced by the patient and associated electronic systems in passively scattered proton therapy treatment. Methods: A multi-stage neutron shielding system was retrofitted to the Gantry 1 treatment nozzle at Loma Linda University Medical Center. The system uses multiple borated polyethylene plates staged after the primary beam modifying devices to attenuate and absorb neutrons produced by such devices. This arrangement locates increasing levels of shielding between the sources of secondary particles in the nozzle and the patient. Additionally, the design of this shielding structure allows it to be easily retrofitted to an existing proton nozzle system without impacting design or treatment beam characteristics. The effectiveness of the shielding was evaluated both through experimental measurements and Geant4 Monte Carlo simulations. Results: Measurements were completed with Landauer Luxel+ dosimeters that use optically stimulated luminescence and CR-39 to detect fast neutrons, thermal neutrons, protons, photons and beta particles. Measurements of a 250 MeV proton beam indicated that the shielding system reduced out-of-field dose to the patient by almost half with dose equivalent values at 50 and 40 cm from the field edge decreasing from 0.965 and 1.262 mSv/Gy to 0.596 and 0.777 mSv/Gy respectively. The installation of the multi-stage shielding system also reduced dose equivalent experienced by electronic systems installed in the treatment room by up to 80%. Geant4 simulations were also used to evaluate the neutron fluence at various positions in the treatment room as well as provide information on microdosimetry spectra within the patient and treatment room. Conclusion: The shielding system described above proved to be an effective an inexpensive method of reducing out-of-field doses to the patient and electronic systems and can be easily retrofitted to existing passive scattering nozzles.

60 APPLIED LIFE SCIENCES↗

Kβ X-ray Emission Spectra Analysis Using Bayesian Optimization

The Kβ X-ray emission spectrum of 3 d transition metals is rich with electronic and structural information due to strong exchange interactions with the valence shell of the metal, and has become crucial for understanding their spin and oxidation states. The spectrum is commonly treated using crystal-field multiplet theory, a semi-empirical theory that uses tunable parameters to control the strength of the effects present in X-ray emission spectroscopy (XES). However, determining the experimental values of these parameters remains a challenge. We present a methodology that applies Bayesian optimization to crystal-field multiplet theory to determine parameter values. The algorithm is tested on the X-ray emission spectra of a collection of Mn, Co, and Ni oxides. We are able to find optimal values for the four most impactful parameters: Slater−Condon reduction factors F dd , F pd , and G pd , and crystal field splitting 10 Dq . The algorithm produces significantly improved accuracy compared to current analysis methods, and probes interparameter dependencies by modeling the error landscape. This advancement enhances XES analysis by offering an approach of obtaining quantitative electronic structural information on 3 d transition metal valence shells, facilitating applications across various scientific fields.

Bayesian optimization↗

Charge-transfer excited states: Seeking a balanced and efficient wave function ansatz in variational Monte Carlo

We present a simple and efficient wave function ansatz for the treatment of excited charge-transfer states in real-space quantum Monte Carlo methods. Using the recently-introduced variation-after-response method, this ansatz allows a crucial orbital optimization step to be performed beyond a configuration interaction singles expansion, while only requiring calculation of two Slater determinant objects. As a result, we demonstrate this ansatz for the illustrative example of the stretched LiF molecule, for a range of excited states of formaldehyde, and finally for the more challenging ethylene-tetrafluoroethylene molecule.

74 ATOMIC AND MOLECULAR PHYSICS↗

QMCPACK : an open source ab initio quantum Monte Carlo package for the electronic structure of atoms, molecules and solids

QMCPACK is an open source quantum Monte Carlo package for ab-initio electronic structure calculations. It supports calculations of metallic and insulating solids, molecules, atoms, and some model Hamiltonians. Implemented real space quantum Monte Carlo algorithms include variational, diffusion, and reptation Monte Carlo. QMCPACK uses Slater-Jastrow type trial wave functions in conjunction with a sophisticated optimizer capable of optimizing tens of thousands of parameters. The orbital space auxiliary field quantum Monte Carlo method is also implemented, enabling cross validation between different highly accurate methods. The code is specifically optimized for calculations with large numbers of electrons on the latest high performance computing architectures, including multicore central processing unit (CPU) and graphical processing unit (GPU) systems. We detail the program’s capabilities, outline its structure, and give examples of its use in current research calculations. The package is available at http://www.qmcpack.org.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Robust Half-Metallicity in Quaternary Heusler Compounds KSrNZ (Z = O, S, and Se)

The first principle calculations based on density functional theory using full potential linearized augmented plane wave method was performed on the quaternary Heusler alloys KSrNZ (Z = O, S, and Se). Studies on structural properties showed that quaternary Heusler alloys KSrNZ (Z = O, S, and Se) were stable in Y{sub I}-type structure. Furthermore, values of formation energies confirmed that these alloys may be synthesized experimentally. Verification of electronic properties resulted that KSrNO, KSrNS, and KSrNSe compounds were half-metals at equilibrium lattice constants with majority band gaps of 4.61, 3.82, and 3.51 eV, respectively. KSrNZ (Z = O, S, and Se) compounds maintained half-metallic behavior in relatively wide ranges of lattice constant under uniform strain distortion. This indicates that these compounds are promising materials for future spintronic applications. Furthermore, the origin of half-metallic characteristic was verified using electronic density of states. Calculations of magnetic properties showed that the total magnetic moments of KSrNZ (Z = O, S, and Se) were obtained equal to 2 μ{sub B} satisfying the famous Slater-Pauling rule (M{sub tot}= 12-Z{sub tot}).

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Communication: On the stability of ice 0, ice i, and I h

Using ab initio methods, we examine the stability of ice 0, a recently proposed tetragonal form of ice implicated in the homogeneous freezing of water [J. Russo, F. Romano, and H. Tanaka, Nat. Mater. 13, 670 (2014)]. Vibrational frequencies are computed across the complete Brillouin Zone using Density Functional Theory (DFT), to confirm mechanical stability and quantify the free energy of ice 0 relative to ice I h . The robustness of this result is tested via dispersion corrected semi-local and hybrid DFT, and Quantum Monte-Carlo calculation of lattice energies. Results indicate that popular molecular models only slightly overestimate the stability of ice zero. In addition, we study all possible realisations of proton disorder within the ice zero unit cell, and identify the ground state as ferroelectric. Comparisons are made to other low density metastable forms of ice, suggesting that the ice i structure [C. J. Fennel and J. D. Gezelter, J. Chem. Theory Comput. 1, 662 (2005)] may be equally relevant to ice formation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗