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 127 records · Page 7

Digital implementation of vibration suppression controllers for LSS

A synthesis procedure for positivity controllers is described, and both simulation and experimental results for vibration suppression control experiments at the NASA Marshall ACES and the NASA Langley MiniMAST ground test facilities are presented. Positivity-designed controllers provide general stability robustness to modeling errors, coupled with controller phase compensation which provides stability margins required for actual implementation. A multivariable scaling method of positivity-designed controllers is used for the on-line tuning of nominal digital controllers. A phase compensation method uses model information to predict system measurements and controller states for compensating the nominal controllers for unmodeled time delays in the system.

Bosse, A.↗

SUPIN: A Tool for the Aerodynamic Design and Analysis of Supersonic Inlets

The Supersonic Inlet Design and Analysis Tool (SUPIN) is a computational program that performs geometric modeling and aerodynamic design and analysis of a set of supersonic inlets whose types include the axisymmetric pitot, three-dimensional pitot, axisymmetric, two-dimensional, two-dimensional bifurcated-duct, and streamline-traced inlets. The inlets are modeled by dividing the inlet flowpath into components that start at the freestream and ends at the engine face. The primary components of the inlet include the external supersonic diffuser, throat section, and subsonic diffuser. Each component is characterized by a set of inputs that include geometric and aerodynamic factors. Examples of the geometric factors are angles, lengths, area ratios, and coordinates that describe the geometry of the components. Examples of the aerodynamic factors are Mach numbers, flow ratios, and pressure ratios that set constraints or desired aerodynamic properties of the component. These input factors are specified within a text-based input data file that is read into SUPIN. SUPIN is a Fortran 95 program that performs the inlet design and analysis as a batch or non-interactive process. SUPIN writes output data files describing the inlet geometry and aerodynamic performance. SUPIN uses analytic, empirical, and numerical methods to perform the inlet design and analysis. The inlet geometry is constructed using planar and three-dimensional geometry constructs and is represented by coordinates, angles, areas, profiles, and surfaces. Surfaces are output in the form of a Plot3D surface grid file. SUPIN does have some capability to automatically generate planar and three-dimensional, multi-block, structured grids for computational fluid dynamics (CFD) analysis. The aerodynamic analysis performed by SUPIN computes quasi-one-dimensional properties, such as Mach numbers, pressures, and temperatures, through the inlet flowpath. Planar aerodynamic properties, such as those for shock and Mach waves, are computed for the supersonic compression of the external supersonic diffuser and throat section. The inlet performance is characterized by the inlet flow rates, total pressure recovery, and inlet drag. The primary intent of this document is to describe the usage of SUPIN. The methods used within SUPIN are documented and may serve as a useful reference for understanding the aerodynamic design and analysis of supersonic inlets.

supersonic inlets↗

On extending Kohn-Sham density functionals to systems with fractional number of electrons

Here, we analyze four ways of formulating the Kohn-Sham (KS) density functionals with a fractional number of electrons, through extending the constrained search space from the Kohn-Sham and the generalized Kohn-Sham (GKS) non-interacting $\mathcal{v}$-representable density domain for integer systems to four different sets of densities for fractional systems. In particular, these density sets are (I) ensemble interacting N -representable densities, (II) ensemble non-interacting N -representable densities, (III) non-interacting densities by the Janak construction, and (IV) non-interacting densities whose composing orbitals satisfy the Aufbau occupation principle. By proving the equivalence of the underlying first order reduced density matrices associated with these densities, we show that sets (I), (II), and (III) are equivalent, and all reduce to the Janak construction. Moreover, for functionals with the ensemble v-representable assumption at the minimizer, (III) reduces to (IV) and thus justifies the previous use of the Aufbau protocol within the (G)KS framework in the study of the ground state of fractional electron systems, as defined in the grand canonical ensemble at zero temperature. By further analyzing the Aufbau solution for different density functional approximations (DFAs) in the (G)KS scheme, we rigorously prove that there can be one and only one fractional occupation for the Hartree Fock functional, while there can be multiple fractional occupations for general DFAs in the presence of degeneracy. This has been confirmed by numerical calculations using the local density approximation as a representative of general DFAs. This work thus clarifies important issues on density functional theory calculations for fractional electron systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Microscopically based energy density functionals for nuclei using the density matrix expansion. II. Full optimization and validation

Background: Energy density functional methods provide a generic framework to compute properties of atomic nuclei starting from models of nuclear potentials and the rules of quantum mechanics. Until now, the overwhelming majority of functionals have been constructed either from empirical nuclear potentials such as the Skyrme or Gogny forces, or from systematic gradient-like expansions in the spirit of the density functional theory for atoms. Purpose: In this study, we seek to obtain a usable form of the nuclear energy density functional that is rooted in the modern theory of nuclear forces. We thus consider a functional obtained from the density matrix expansion of local nuclear potentials from chiral effective field theory. We propose a parametrization of this functional carefully calibrated and validated on selected ground-state properties that is suitable for large-scale calculations of nuclear properties. Methods: Our energy functional comprises two main components. The first component is a non-local functional of the density and corresponds to the direct part (Hartree term) of the expectation value of local chiral potentials on a Slater determinant. Contributions to the mean field and the energy of this term are computed by expanding the spatial, finite-range components of the chiral potential onto Gaussian functions. The second component is a local functional of the density and is obtained by applying the density matrix expansion to the exchange part (Fock term) of the expectation value of the local chiral potential. We apply the UNEDF2 optimization protocol to determine the coupling constants of this energy functional. Results: We obtain a set of microscopically constrained functionals for local chiral potentials from leading order up to next-to-next-to-leading order with and without three-body forces and contributions from Δ excitations. These functionals are validated on the calculation of nuclear and neutron matter, nuclear mass tables, single-particle shell structure in closed-shell nuclei, and the fission barrier of 240 Pu . Quantitatively, they perform noticeably better than the more phenomenological Skyrme functionals. Conclusions: The inclusion of higher-order terms in the chiral perturbation expansion seems to produce a systematic improvement in predicting nuclear binding energies while the impact on other observables is not really significant. In conclusion, this result is especially promising since all the fits have been performed at the single-reference level of the energy density functional approach, where important collective correlations such as center-of-mass correction, rotational correction, or zero-point vibrational energies have not been taken into account yet.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Magneto-electronic and Thermodynamic Properties of Quaternary NiFeMnZ (Z = In, Sn) New Spin Gapless Semiconductors

This paper deals with investigation of magneto-electronic and thermodynamic properties of new spin gapless semiconductor quaternary NiFeMnZ (Z = In, Sn) Heusler alloys. In addition, electronic properties such as band structure and density of states were studied in detail. The obtained results are based on first-principles calculations, using the full potential linearized augmented plane wave (FP-LAPW) method within well-known generalized gradient approximation (GGA) coupling with spin-orbit calculation (SOC) in framework of the density functional theory (DFT) implemented in the WIEN2k package. As results, the alloys were found to be new spin gapless semiconductors (SGSs), and both compounds in LiMgPdSn-type crystal structure are stable. The magnetic moment of NiFeMnIn and NiFeMnSn was 4.39 and 4.99 μ{sub B} per formula unit, respectively. This result is in agreement with the Slater-Pauling rule (M{sub tot} = 24 −Z{sub tot}). Furthermore, the thermodynamic properties were also predicted. All these multifunctionalities lead to crucial interest for developing future spintronics device generations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Ferromagnetism in half-metallic quaternary FeVTiAl Heusler compound

The electronic structure and magnetic properties of FeVTiAl quaternary Heusler alloy have been investigated within the density functional theory framework. The material was found completely spin-polarized half-metallic Ferromagnet in the ground state with F-43m structure. The structural stability was further confirmed by calculating different elastic constants in the cubic phase. Present study predicts an energy band gap of 0.72 eV calculated in localized minority spin channel at an equilibrium lattice parameter of 6.0Å. The calculated total spin magnetic moment of 2 µ{sub B}/f.u. is in agreement with the Slater-Pauling rule for full Heusler alloys.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

First-Principle Study of Half-Metallic Ferrimagnet Behavior in Titanium-Based Heusler Alloys Ti{sub 2}FeZ (Z = Al, Ga, and In)

Using density functional theory with the full-potential linearized augmented plane-wave method (FP-LAPW), we have study the structural, electronic, and magnetic properties of Ti{sub 2}FeZ (Z = Al, Ga, and In) alloys with Hg{sub 2}CuTi-type structure. The magnetic stabilities reveal that all our compounds exhibit ferrimagnetic (FiM) behaviors. The electronic structure report the existence of a gap energy equal to 0.56, 0.60, and 0.64 eV for Ti{sub 2}FeAl, Ti{sub 2}FeGa, and Ti{sub 2}FeIn, respectively, in the spin-down state and divulge metallic intersections at the Fermi level for the spin-up state. These results indicate that our compounds have a half-metallic (HM) nature. In addition to this, the total magnetic moments are in agreement with the obtained one by the Slater-Pauling rule (M{sub tot} = Z{sub tot} − 18), which indicates the 100% spin polarization for these compounds.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Half-metallic Ferromagnetism in Novel Rh{sub 2}-based Full-Heusler Alloys Rh{sub 2}FeZ (Z = Ga and In)

In this approach, we have employed the ab initio method of full-potential linearized plane waves plus local orbital based on the framework of density functional theory, within the spin generalized gradient approximation plus U (GGA+U), where U is the Hubbard correlation term for exchange-correlation potential. This method is used to predict the equilibrium structural parameters, the electronic structure, and the magnetic properties of Rh{sub 2}FeGa and Rh{sub 2}FeIn full-Heusler alloys. The structural properties result show that both Rh{sub 2}FeGa and Rh{sub 2}FeIn compounds are stable in ferromagnetic states. The spin-polarized electronic structure demonstrate that both full-Heusler compounds are half-metals with half-metallic band gap (E{sub HM}) of 0.326 eV for Rh{sub 2}FeGa alloy and 0.245 eV for Rh{sub 2}FeIn alloy, where the apparition of an indirect gap from Γ to X in minority spin is shown by these Heusler alloys. Furthermore, calculations of total magnetic moment give the value of 5 µ{sub B}, which agrees with Slater-Pauling rule; also, the atomic magnetic moment of both Rh and Fe transition metals have the same sign which depicts the ferromagnetic behavior.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Atomic emission spectroscopy

The relationship between the Slater-Condon theory and the conditions within the atom as revealed by experimental data was investigated. The first spectrum of Si, Rb, Cl, Br, I, Ne, Ar, and Xe-136 and the second spectrum of As, Cu, and P were determined. Methods for assessing the phase stability of fringe counting interferometers and the design of an autoranging scanning system for digitizing the output of an infrared spectrometer and recording it on magnetic tape are described.

Andrew, K. H.↗

SU-F-T-211: Evaluation of a Dual Focusing Magnet System for the Treatment of Small Proton Targets

Purpose: To investigate magnetic focusing for small volume proton targets using a doublet combination of quadrupole rare earth permanent magnet Halbach cylinder assemblies Methods: Monte Carlo computer simulations were performed using the Geant4 toolkit to compare dose depositions of proton beams transported through two focusing magnets or in their absence. Proton beams with energies of 127 MeV and initial diameters of 5, 8 and 10 mm were delivered through two identical focusing magnets similar to those currently in experimental use at Loma Linda University Medical Center. Analogous experiments used optimized configurations based on the simulation results. Dose was measured by a diode detector and Gafchromic EBT3 film and compared to simulation data. Based on results from the experimental data, an additional set of simulations was performed with an initial beam diameter of 18 mm and a two differing length magnets (40mm & 68mm). Results: Experimental data matched well with Monte Carlo simulations. However, under conditions necessary to produce circular beam spots at target depth, magnetically focused beams using two identical 40 mm length magnets did not meet all of our performance criteria of circular beam spots, improved peak to entrance (P/E) dose ratios and dose delivery efficiencies. The simulations using the longer 68 mm 2nd magnet yielded better results with 34% better P/E dose ratio and 20–50% better dose delivery efficiencies when compared to unfocused 10 mm beams. Conclusion: While magnetic focusing using two magnets with identical focusing power did not yield desired results, ongoing Monte Carlo simulations suggest that increasing the length of the 2nd magnet to 68 mm could improve P/E dose ratios and dose efficiencies. Future work includes additional experimental validation of the longer 2nd magnet setup as well as experiments with triplet magnet systems. This project was sponsored with funding from the Department of Defense (DOD# W81XWH-BAA-10-1).

60 APPLIED LIFE SCIENCES↗

X-ray production cross sections at incident photon energies across the M{sub i} (i=1-5) edges of {sub 90}Th

The X-ray production cross sections for the M{sub k} (k= ξ, δ, α, β, ζ, γ, m{sub 1}, m{sub 2}) groups of X-rays have been evaluated at incident photon energies across the M{sub i} (i =1-5) edges of {sub 90}Th using the relativistic Hartree-Fock-Slater model based photoionisation cross sections and recently reported values of the M-shell X-ray emission rates, fluorescence and Coster Kronig yields. Further, the energies of the prominent (M{sub i}-S{sub j}) (S{sub j}=N{sub j}, O{sub j} and i =1-3, j =1-7) resonant Raman scattered (RRS) peaks at different incident photon energies have also been evaluated using the neutral-atom electron binding energies (E{sub sj}) based on the relaxed orbital relativistic Hartree-Fock-Slater model.

74 ATOMIC AND MOLECULAR PHYSICS↗

Sub-system self-consistency in coupled cluster theory

Here, in this article, we provide numerical evidence indicating that the single-reference coupled-cluster (CC) energies can be calculated alternatively to their copybook definition. We demonstrate that the CC energy can be reconstructed by diagonalizing the effective Hamiltonians describing correlated sub-systems of the many-body system. In the extreme case, we provide numerical evidence that the CC energy can be reproduced through the diagonalization of the effective Hamiltonian describing sub-system composed of a single electron. These properties of the CC formalism can be exploited to design protocols to define effective interactions in sub-systems used as probes to calculate the energy of the entire system and introduce a new type of self-consistency for approximate CC approaches.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The accuracies of effective interactions in downfolding coupled-cluster approaches for small-dimensionality active spaces

Here, this paper evaluates the accuracy of the Hermitian form of the downfolding procedure using the double unitary coupled cluster (DUCC) ansatz on the benchmark systems of linear chains of hydrogen atoms, H6 and H8. The computational infrastructure employs the occupation-number-representation codes to construct the matrix representation of arbitrary second-quantized operators, allowing for the exact representation of exponentials of various operators. The tests demonstrate that external amplitudes from standard single-reference coupled cluster methods that sufficiently describe external (out-of-active-space) correlations reliably parameterize the Hermitian downfolded effective Hamiltonians in the DUCC formalism. The results show that this approach can overcome the problems associated with losing the variational character of corresponding energies in the corresponding SR-CC theories.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Density Functional Study of Half-Metallic Ferromagnetism in the New Heusler Compounds CaYO{sub 2} (Y = Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, and Zn)

The structural, electronic, and magnetic properties of full-Heusler alloys CaYO{sub 2} (Y = Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, and Zn) were investigated using full-potential linearized augmented plane wave method based on density functional theory. It was confirmed that these compounds can be synthesized experimentally due to their negative formation energies. The results predicted that all compounds in AlCu{sub 2}Mn-type structure were more stable than CuHg{sub 2}Ti-type structure. The band structure calculations showed that alloys CaYO{sub 2} (Y = Cr, Fe, Ni, and Cu) in AlCu{sub 2}Mn-type structure and CaYO{sub 2} (Y = Ti, V, Cr, Co, and Ni) in CuHg{sub 2}Ti-type structure have half metallic behavior. The origin of half-metallic characteristic was investigated using electronic density of states for CaFeO{sub 2}. Among these alloys, CaCrO{sub 2} was half-metal in a wide range of lattice constants indicating its resistance against destructive factors of half-metallicity such as temperature and pressure. The total magnetic moments for half-metallic alloys were integer values and obeyed the Slater-Pauling rules. The CaYO{sub 2} alloys with half-metallic character may be good candidates for practical applications in spintronics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Tunable-fidelity wave functions for the ab initio description of scattering and reactions

Here, the no-core shell model (NCSM) is an ab initio method that solves the nuclear many-body problem by expanding the many-particle wave function into a (typically) harmonic oscillator basis and minimizing the energy to obtain the expansion coefficients. Extensions of the NCSM, such as its coupling with microscopic-cluster basis states, further allow for an ab initio treatment of light-ion nuclear reactions of interest for both astrophysics and nuclear technology applications. A downside of the method is the exponential scaling of the basis size with increasing number of nucleons and excitation quanta, which limits its applicability to mass A ≲ 16 nuclei, except for variants where the basis is further down-selected via some truncation scheme. We consider a basis selection method for the NCSM that was first introduced in the context of the large-scale shell model and captures the essential degrees of freedom of the nuclear wave function leading to a favorable complexity scaling for calculations and enabling ab initio reaction calculations in sd-shell nuclei. The particle configurations within the NCSM basis are ordered based on their contribution to the first moment of the Hamiltonian matrix that results from the projection onto the many-body basis. The truncation scheme then consists in retaining only the lowest-first-moment configurations, which typically contain only few many-body basis states (Slater determinants). As the energy threshold above which configurations are disregarded is increased, the size of the basis becomes an almost-continuous variable, allowing for tunable fidelity in the obtained wave functions. The resulting wave functions can then be used directly in ab initio reaction calculations. We present calculations for 7 Li and n + 12 C scattering using nucleon-nucleon interactions derived from chiral effective field theory and softened using the similarity renormalization group method. The obtained energy levels invariably demonstrate exponential convergence with the size of the basis, and we find improved convergence in scattering calculations. To demonstrate the possibilities enabled by the approach, we also present a first calculation for the scattering of neutrons from 24 Mg. The method presented in this work appears promising for future studies of nuclei with mass A > 16, opening multiple future research directions impacting both nuclear astrophysics and nuclear technology applications.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

SU-E-T-72: Commissioning of a Standardized SRS Cone Set: Determination of the Bolus Gap Factors in a Passively Scattered Proton Beam

Purpose: To determine the proton output factors for an SRS cone set using standardized apertures and varied range compensators (bolus blanks); specifically, to determine the best method for modeling the bolus gap factor (BGF) and eliminate the need for patient specific calibrations. Methods: A Standard Imaging A-16 chamber was placed in a Plastic Water phantom to measure the change in dose/MU with different treatment combinations for a proton SRS cone, using standardized apertures and range compensators. Measurements were made with all apertures in the SRS cone set, with four different range compensator thicknesses and five different air gaps between the end of the SRS cone and the surface of the phantom. The chamber was located at iso-center and maintained at a constant depth at the center of modulation for all measurements. Each aperture was placed in the cone to measure the change in MU needed to maintain constant dose at the chamber, as the air gap was increased with different thicknesses of bolus. Results: The dose/MU varied significantly with decreasing aperture size, increasing bolus thickness, or increasing air gap. The measured data was fitted with the lowest order polynomials that accurately described the data, to create a model for determining the change in output for any potential combination of devices used to treat a patient. For a given standardized aperture, the BGF could be described by its constituent factors: the bolus thickness factor (BTF) and the nozzle extension factor (NEF). Conclusion: The methods used to model the dose at the calibration point could be used to accurately predict the change in output for SRS proton beams due to the BGF, eliminating the need for patient specific calibrations. This method for modeling SRS treatments could also be applied to model other treatments using passively scattered proton beams.

07 ISOTOPE AND RADIATION SOURCES↗

Efficient core-excited state orbital perspective on calculating x-ray absorption transitions in determinant framework

X-ray absorption spectroscopy (XAS) is an explicit probe of the unoccupied electronic structure of materials and an invaluable tool for fingerprinting various electronic properties and phenomena. Computational methods capable of simulating and analyzing such spectra are therefore in high demand for complementing experimental results and for extracting valuable insights therefrom. In particular, a recently proposed first-principles approach titled many-body XAS (MBXAS), which approximates the final (initial) state as a Slater determinant constructed from Kohn-Sham (KS) orbitals optimized in the absence (presence) of the relevant core electron has shown promising prospects in evaluating the transition amplitudes. In this article, we show that the MBXAS approach can be rederived using a transition operator expressed entirely in the basis of core-excited state KS orbitals and that this reformulation offers substantial practical and conceptual advantages. In addition to circumventing previous issues of convergence with respect to the number of unoccupied ground-state orbitals, the aforementioned representation reduces the computational expense by rendering the calculation of such orbitals unnecessary altogether. Furthermore, the reformulated approach also provides a direct pathway for comparing the many-body approximation with the so-called single-particle treatment and indicates the relative importance in observed XAS intensity of the relaxation of the valence occupied subspace induced by the core excitation. Finally, using the core-excited state basis, we define auxiliary orbitals for x-ray absorption and demonstrate their utility in explaining the spectral intensity by contrasting them with single-particle approximations to the excited state.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

A Density Functional Theory Investigations of Half-Heusler Compounds RhVZ (Z = P, As, Sb)

We have studied the half-Heusler compounds RhVZ (Z = P, As, Sb) using Density Functional Theory (DFT). The method of the Full Potential Linearly Augmented Plane Wave (FP-LAPW) is employed in the Wien2k package for calculation of structural, elastic, mechanical, electronic, magnetic, and optical properties. Lattice constants are found in the range 5.67 − 5.80 Å. Analysis of elastic properties shows that two of the compounds namely RhVP and RhVAs are ductile while RhVSb is brittle in nature. All three compounds are half-metals as revealed from the band structures and Density of States (DoS) calculations. In all of them, spin-down channels have the small band gaps, while spin-up channels are conducting. These compounds follow the Slater-Pauling 18 (M{sub Tot} = Z{sub Tot} − 18) electron rule with the total magnetic moments in the range 1 − 2 μ{sub B}. Optical properties like dielectric function, refractive index, reflectivity, conductivity, and absorption coefficient are calculated and discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗