Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Formal Methods”

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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

SU-F-T-555: Accurate Stereotactic Cone TMRs Converted from PDDs Scanned with Ray Trace

Purpose: To investigate whether the accuracy of TMRs for stereotactic cones converted from PDDs scanned with Ray Trace can be improved, when compared against the TMRs converted from the traditional PDDs. Methods: Ray Trace measurement in Sun Nuclear 3D Scanner is for accurate scan of small field PDDs. The system detects the center of field at two depths, for example, at 3 and 20 cm in our study, and then performs scan along the line passing the two centers. With both Ray Trace and the traditional method, PDDs for conical cones of 4, 5, 7.5, 10, 12.5, 15, and 17.5 mm diameter (jaws set to 5×5 cm) were obtained for 6X FFF and 10X FFF energies on a Varian Edge linac, using Edge detectors. The formalism of converting PDD to TMR given in Khan’s book (4th Edition, p.161) was applied. Sp values at dmax were obtained by measuring cone Scp and Sc. Continuous direct measurement of TMR by filling/draining water to/from the tank and spot measurement by moving the tank and detector were also performed with the same equipment, using 100 cm SDD. Results: For 6XFFF energy and all the cones, TMRs converted from Ray Trace were very close to the continuous and spot measurement, while TMRs converted from traditional PDDs had larger deviation. Along the central axis beyond dmax, 1.7% of TMR data points calculated from Ray Trace had more 3% deviation from measurement, with maximal deviation of 5.2%. Whereas, 34% of TMR points calculated from traditional PDDs had more than 3% deviation, with maximum of 5.7%. In this initial study, Ray Trace scans for 10XFFF beam were noisy, further measurement is warranted. Conclusion: The Ray Trace could improve the accuracy of PDDs measurement and the calculated TMRs for stereotactic cones, which was within 3% of the measured TMRs.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

Condensation of Pseudomagnons in a Two-Dimensional Anisotropic S = 1 Pseudospin System

Abstract—A 2D anisotropic system of S = 1 centers of the charge triplet type in systems with variable valence or “semi-hard-core” boson systems with a limitation for the occupation of lattice sites n = 0, 1, 2 is studied in the framework of the pseudospin formalism. Assuming that the ground state is a quantum paramagnet, the pseudo-spin wave spectrum and also the conditions of the condensation of pseudomagnons with a phase transition to a superconducting state have been found using the Schwinger boson method.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Twisted 3D N=4 supersymmetric YM on deformed A{sub 3}{sup *} lattice

We study a class of twisted 3D N=4 supersymmetric Yang-Mills (SYM) theory on particular 3-dimensional lattice L{sub 3D} formally denoted as L{sub 3D}{sup su{sub 3}×u{sub 1}} and given by non-trivial fibration L{sub 1D}{sup u{sub 1}}×L{sub 2D}{sup su{sub 3}} with base L{sub 2D}{sup su{sub 3}}=A{sub 2}{sup *}, the weight lattice of SU(3). We first, develop the twisted 3D N=4 SYM in continuum by using superspace method where the scalar supercharge Q is manifestly exhibited. Then, we show how to engineer the 3D lattice L{sub 3D}{sup su{sub 3}×u{sub 1}} that host this theory. After that we build the lattice action S{sub latt} invariant under the following three points: (i) U(N) gauge invariance, (ii) BRST symmetry, (iii) the S{sub 3} point group symmetry of L{sub 3D}{sup su{sub 3}×u{sub 1}}. Other features such as reduction to twisted 2D supersymmetry with 8 supercharges living on L{sub 2D}≡L{sub 2D}{sup su{sub 2}×u{sub 1}}, the extension to twisted maximal 5D SYM with 16 supercharges on lattice L{sub 5D}≡L{sub 5D}{sup su{sub 4}×u{sub 1}} as well as the relation with known results are also given.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Quantum Monte Carlo formalism for dynamical pions and nucleons

In most simulations of nonrelativistic nuclear systems, the wave functions found solving the many-body Schrodinger equations describe the quantum-mechanical amplitudes of the nucleonic degrees of freedom. In those simulations the pionic contributions are encoded in nuclear potentials and electroweak currents, and they determine the low-momentum behavior. In this work we present an alternative quantum Monte Carlo formalism in which both relativistic pions and nonrelativistic nucleons are explicitly included in the quantum-mechanical states of the system. We report the renormalization of the nucleon mass as a function of the momentum cutoff, a Euclidean time density correlation function that deals with the short-time nucleon diffusion, and the pion cloud density and momentum distributions. In the two-nucleon sector we show that the interaction of two static nucleons at large distances reduces to the one-pion exchange potential, and we fit the low-energy constants of the contact interactions to reproduce the binding energy of the deuteron and two neutrons in finite volumes. We show that the method can be readily applied to light-nuclei.

Madeira, Lucas↗
DOE OSTI2014

Calculation of positron binding energies using the generalized any particle propagator theory

We recently extended the electron propagator theory to any type of quantum species based in the framework of the Any-Particle Molecular Orbital (APMO) approach [J. Romero, E. Posada, R. Flores-Moreno, and A. Reyes, J. Chem. Phys. 137, 074105 (2012)]. The generalized any particle molecular orbital propagator theory (APMO/PT) was implemented in its quasiparticle second order version in the LOWDIN code and was applied to calculate nuclear quantum effects in electron binding energies and proton binding energies in molecular systems [M. Díaz-Tinoco, J. Romero, J. V. Ortiz, A. Reyes, and R. Flores-Moreno, J. Chem. Phys. 138, 194108 (2013)]. In this work, we present the derivation of third order quasiparticle APMO/PT methods and we apply them to calculate positron binding energies (PBEs) of atoms and molecules. We calculated the PBEs of anions and some diatomic molecules using the second order, third order, and renormalized third order quasiparticle APMO/PT approaches and compared our results with those previously calculated employing configuration interaction (CI), explicitly correlated and quantum Montecarlo methodologies. We found that renormalized APMO/PT methods can achieve accuracies of ∼0.35 eV for anionic systems, compared to Full-CI results, and provide a quantitative description of positron binding to anionic and highly polar species. Third order APMO/PT approaches display considerable potential to study positron binding to large molecules because of the fifth power scaling with respect to the number of basis sets. In this regard, we present additional PBE calculations of some small polar organic molecules, amino acids and DNA nucleobases. We complement our numerical assessment with formal and numerical analyses of the treatment of electron-positron correlation within the quasiparticle propagator approach.

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

SU-E-T-237: Monte Carlo Dosimetric Characterization of the Mobetron Mobile Linac

Purpose: The aim of this work is to characterize dosimetrically a clinical intraoperative electron beam accelerator, Mobetron (IntraOp Medical, Inc.) in clinical use in our Hospital. Once this first step is completed our purpose is to evaluate shielding requirements for such a device by preparing adequate phase space files. Methods: It is known that electron beam simulation parameters required for state-of-the-art Monte Carlo codes to obtain a good match with measured data, like the mean energy or the FWHM, may not be code-independent due to the different set of process simulated and formalisms involved. Then, to cross-check our results against any issue in the simulation we have compared experimental data (PDD and profiles for electrons in the range 4 to 12 MeV) with simulations performed independently using both Penelope2011 and Geant4 codes. To do so, the geometry and materials of the head of the accelerator have been fully characterized following information provided by the manufacturer. Results: Both simulations agree with experimental data within experimental uncertainties (±1 mm displacement), although small variations (less than 10%) in the mean energy and FWHM are required to match measured values depending on the code used. Conclusion: Independent Monte Carlo simulations were used to obtain an excellent match to measured electron dose distributions. This opens the road to use such data for evaluating shielding requirements which is the main objective of this project.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2022

Overcoming significant challenges in extracting off-stoichiometric thermodynamics using the compound energy formalism through complementary use of experimental and first principles data: A case study of Ba 1-x Sr x FeO 3-δ

The compound energy formalism (CEF) is a powerful framework to describe the thermodynamics of metal oxides as a function of off-stoichiometry, temperature, and composition. The thermodynamic properties are crucial materials design attributes in metal oxide-based oxygen-exchange chemical processes. Despite the richness of information an accurate CEF model provides, a method to determine a unique and accurate fit for oxygen exchange materials remains elusive. This contribution details a method for fitting the CEF model that overcomes the current fitting challenges through three innovations: 1) the combination of density functional theory calculations with experimental data decorrelates excess terms and delineates the enthalpic/entropic contributions to the Gibbs free energy; 2) a systematic determination of the important CEF model terms, removing thermodynamic predetermining human intervention; 3) a self-consistent solution of the starting oxygen offstoichiometry (δ 0 ) of thermogravimetric measurements. Thus, our method enables the reliable extraction of off-stoichiometric metal oxide thermodynamic properties and facilitates rapid materials compositional screening, and reliable process design of systems dependent on off-stoichiometric redox-active metal oxides. We apply this method to a Ba x Sr 1-x FeO 3-δ test case. We find by systematically examining the performance of the CEF model fit with and without each innovation that all three innovations are necessary for an accurate fit. We determined that reduction enthalpy is higher and more sensitive to off-stoichiometry when the Sr fraction is large (139.5 and 185.3 kJ/mol O 2 for SrFeO 3 at δ = 0 and δ = 0.5, respectively vs. nearly constant 83 kJ/mol O 2 for BaFeO 3 ). However, the reduction entropy is mostly insensitive to Sr fraction, but highly dependent on δ suggesting larger contributions in the non-configurational entropy parameters.

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

Quantum effects on dislocation motion from ring-polymer molecular dynamics

Quantum motion of atoms known as zero-point vibration was recently proposed to explain a long-standing discrepancy between theoretically computed and experimentally measured low-temperature plastic strength of iron and possibly other metals with high atomic masses. This finding challenges the traditional notion that quantum motion of atoms is relatively unimportant in solids comprised of heavy atoms. Here we report quantum dynamic simulations of quantum effects on dislocation motion within the exact formalism of Ring-Polymer Molecular Dynamics (RPMD). To extend the reach of quantum atomistic simulations to length and time scales relevant for extended defects in materials, we implemented RPMD in the open-source code LAMMPS thus making the RPMD method widely available to the community. We use our RPMD/LAMMPS approach for direct calculations of dislocation mobility and its effects on the yield strength of α-iron. Our simulation results establish that quantum effects are noticeable at temperatures below 50 K but account for only a modest (≈13% at T = 0 K) overall reduction in the Peierls barrier, at variance with the factor of two reduction predicted earlier based on the more approximate framework of harmonic transition state theory. Our results confirm that zero-point vibrations provide ample additional agitation for atomic motion that increases with decreasing temperature, however its enhancing effect on dislocation mobility is largely offset by an increase in the effective atom size, an effect known as quantum dispersion that has not been accounted for in the previous calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Hadron production at LHC in dipole momentum space

The dipole color approach is the framework that considers the quark-antiquark pair scattering off the target. The rapidity evolution of color dipoles is given by the nonlinear Balitsky-Kovchegov (BK) equation, for which analytical solutions are not yet known. A good way to explore the asymptotic BK solutions is through the traveling wave method of QCD, that uses a correspondence between the BK evolution equation in momentum space and reaction-diffusion physics. Using the traveling wave based AGBS model for the dipole amplitude in momentum space, and within the k{sub t}-factorization formalism, we describe the LHC data on single inclusive hadron yield for p-p collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

A flexible experimental setup for femtosecond time-resolved broad-band ellipsometry and magneto-optics

A versatile experimental setup for femtosecond time-resolved ellipsometry and magneto-optical Kerr effect measurements in the visible light range is described. The apparatus is based on the pump-probe technique and combines a broad-band probing beam with an intense near-infrared pump. According to Fresnel scattering matrix formalism, the analysis of the reflected beam at different polarization states of the incident probe light allows one to determine the diagonal and the off-diagonal elements of the dielectric tensor in the investigated sample. Moreover, the pump-probe method permits to study the dynamics of the dielectric response after a short and intense optical excitation. The performance of the experimental apparatus is tested on CrO{sub 2} single crystals as a benchmark.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

Freely cooling granular gases with short-ranged attractive potentials

We treat the case of an undriven gas of inelastic hard-spheres with short-ranged attractive potentials via an extension of the pseudo-Liouville operator formalism. New evolution equations for the granular temperature and coordination number are obtained. The granular temperature exhibits deviation from both Haff’s law and the case of long-ranged potentials. We verify this departure using soft-sphere discrete element method simulations. Excellent agreement is found for the duration of the simulation even beyond where exclusively binary collisions are expected. Simulations show the emergence of strong spatial-velocity correlations on the length scale of the last peak in the pair-correlation function but do not show strong correlations beyond this length scale. We argue that molecular chaos may remain an adequate approximation if the system is modelled as a Smoluchowski type equation with aggregation and break-up processes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Reduced scaling Hilbert space variational Monte Carlo

We show that for both single-Slater-Jastrow and Jastrow geminal power wave functions the formal cost scaling of Hilbert space variational Monte Carlo can be reduced from fifth to fourth order in the system size, thus bringing it in line with the long-standing scaling of its real space counterpart. While traditional quantum chemistry methods can reduce costs related to the two-electron integral tensor through various tensor decomposition methods, we show that such approaches are ineffective in the presence of Hilbert space Jastrow factors. Instead, we develop a simple semi-stochastic approach that can take similar advantage of the near-sparsity of this four-index tensor. Through demonstrations on alkanes of increasing length, we show that accuracy and overall statistical uncertainty are not meaningfully affected and that a total cost crossover is reached as early as 50 electrons when using a minimal basis. Further study will be needed to assess where the crossover occurs in more compact molecular geometries and larger basis sets and to explore how in that context the crossover can be accelerated.

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

High pressure structural and lattice dynamics study of α-In2Se3

Layered α-In2Se3 has been studied using a combined in situ synchrotron angle-dispersive powder x-ray diffraction and Raman spectroscopy study in a diamond anvil cell up to 60+ GPa, at room temperature. Helium, which remains fairly hydrostatic up to the highest pressure in this study, was used as the pressure-transmitting medium. The results from both experimental methods reveal a pressure-induced structural phase transition from α-In2Se3 to a monoclinic β'-In2Se3 structure at ≈1 GPa, in agreement with previous studies. Based on our detailed measurements using both experimental techniques and the F-f formalism, the β'-In2Se3 structure remains stable up to 45 GPa, without a clear indication of a phase transition toward the previously reported β-In2Se3 phase. Above this pressure, In2Se3 adopts a disordered solid-solution-like orthorhombic structure, phase IV. The results are discussed in comparison with the relevant previous studies of α-In2Se3 under pressure.

Feng, Shiyu↗
DOE OSTI2015

From the ternary Eu(Au/In) 2 and EuAu 4 (Au/In) 2 with remarkable Au/In distributions to a new structure type: The gold-rich Eu 5 Au 16 (Au/In) 6 structure

The ternary Eu(Au/In) 2 (EuAu 0.46 In 1.54 (2) ) (I), EuAu 4 (Au/In) 2 (EuAu 4+x In 2–x with x = 0.75(2) (II), 0.93(2), and 1.03(2)), and Eu 5 Au 16 (Au/In) 6 (Eu 5 Au 17.29 In 4.71(3) ) (III) have been synthesized, and their structures were characterized by single-crystal X-ray diffraction. I and II crystallize with the CeCu 2 -type (Pearson Symbol oI 12; Imma ; Z = 4; a = 4.9018(4) Å; b = 7.8237(5) Å; c = 8.4457(5) Å) and the YbAl 4 Mo 2 -type (tI14; I4/ mmm ; Z = 2; a = 7.1612(7) Å; c = 5.5268(7) Å) and exhibit significant Au/In disorder. I is composed of an Au/In-mixed diamond-related host lattice encapsulating Eu atoms, while the structure of II features ribbons of distorted, squared Au 8 prisms enclosing Eu, Au, and In atoms. Combination of these structural motifs leads to a new structure type as observed for Eu 5 Au 16 (Au/In) 6 (Eu 5 Au 17.29 In 4.71(3) ) (oS108; Cmcm ; Z = 4; a = 7.2283(4) Å; b = 9.0499(6) Å; c = 34.619(2) Å), which formally represents a one-dimensional intergrowth of the series EuAu 2 –“EuAu 4 In 2 ”. The site preferences of the disordered Au/In positions in II were investigated for different hypothetical “EuAu 4 (Au/In) 2 ” models using the projector-augmented wave method and indicate that these structures attempt to optimize the frequencies of the heteroatomic Au–In contacts. Furthermore, a chemical bonding analysis on two “EuAu 5 In” and “EuAu 4 In 2 ” models employed the TB-LMTO-ASA method and reveals that the subtle interplay between the local atomic environments and the bond energies determines the structural and site preferences for these systems.

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

A linear response framework for quantum simulation of bosonic and fermionic correlation functions

Response functions are a fundamental aspect of physics; they represent the link between experimental observations and the underlying quantum many-body state. However, this link is often under-appreciated, as the Lehmann formalism for obtaining response functions in linear response has no direct link to experiment. Within the context of quantum computing, and via a linear response framework, we restore this link by making the experiment an inextricable part of the quantum simulation. This method can be frequency- and momentum-selective, avoids limitations on operators that can be directly measured, and can be more efficient than competing methods. As prototypical examples of response functions, we demonstrate that both bosonic and fermionic Green’s functions can be obtained, and apply these ideas to the study of a charge-density-wave material on the ibm_auckland superconducting quantum computer. The linear response method provides a robust framework for using quantum computers to study systems in physics and chemistry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Designing field-controllable graphene-dot-graphene single molecule switches: A quantum-theoretical proof-of-concept under realistic operating conditions

A theoretical proof of the concept that a particularly designed graphene-based moletronics device, constituted by two semi-infinite graphene subunits, acting as source and drain electrodes, and a central benzenoid ring rotator (a “quantum dot”), could act as a field-controllable molecular switch is outlined and analyzed with the density functional theory approach. Besides the ideal (0 K) case, we also consider the operation of such a device under realistic operating (i.e., finite-temperature) conditions. An in-depth insight into the physics behind device controllability by an external field was gained by thorough analyses of the torsional potential of the dot under various conditions (absence or presence of an external gating field with varying strength), computing the torsional correlation time and transition probabilities within the Bloembergen-Purcell-Pound formalism. Both classical and quantum mechanical tunneling contributions to the intramolecular rotation were considered in the model. The main idea that we put forward in the present study is that intramolecular rotors can be controlled by the gating field even in cases when these groups do not possess a permanent dipole moment (as in cases considered previously by us [I. Petreska et al., J. Chem. Phys. 134, 014708-1–014708-12 (2011)] and also by other groups [P. E. Kornilovitch et al., Phys. Rev. B 66, 245413-1–245413-7 (2002)]). Consequently, one can control the molecular switching properties by an external electrostatic field utilizing even nonpolar intramolecular rotors (i.e., in a more general case than those considered so far). Molecular admittance of the currently considered graphene-based molecular switch under various conditions is analyzed employing non-equilibrium Green’s function formalism, as well as by analysis of frontier molecular orbitals’ behavior.

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

Analysis of Three Body Decays in Quasi-Real Photoproduction

This thesis presents preliminary measurements of the moments of angular distribution for the K?K+ mesonic final state using a quasi-real, linearly polarised photon. Two of the main points of focus for this analysis were the development of the formalism for photoproduced vector-pseudoscalar decay and the tools by which the analysis was carried out. The formalism described how the moments of angular distribution could be extracted from the angular decay variables, and how they related to resonance decays in terms of partial waves. In extracting the moments, it was found that adapted Markov Chain Monte Carlo (MCMC) methods proved to be more effective in terms of extracting results from the data compared to the sole use of gradient descent based fitting algorithms. These measurements are complementary to other decay channels and production mechanisms, with the global aim of probing the existence of exotic mesons. Data from this analysis was taken by the CLAS collaboration at the Thomas Jefferson National Accelerator Facility (JLab) by scattering from a liquid hydrogen target using a 10.6 GeV electron. This work was done within the MesonEx program at CLAS12 that aims to map the spectrum of mesons, and in doing so, gain a greater understanding of Quantum ChromoDynamics (QCD) and the strong force.

Wishart, Robert↗
DOE OSTI2015

Forces and stress in second order Møller-Plesset perturbation theory for condensed phase systems within the resolution-of-identity Gaussian and plane waves approach

The forces acting on the atoms as well as the stress tensor are crucial ingredients for calculating the structural and dynamical properties of systems in the condensed phase. Here, these derivatives of the total energy are evaluated for the second-order Møller-Plesset perturbation energy (MP2) in the framework of the resolution of identity Gaussian and plane waves method, in a way that is fully consistent with how the total energy is computed. This consistency is non-trivial, given the different ways employed to compute Coulomb, exchange, and canonical four center integrals, and allows, for example, for energy conserving dynamics in various ensembles. Based on this formalism, a massively parallel algorithm has been developed for finite and extended system. The designed parallel algorithm displays, with respect to the system size, cubic, quartic, and quintic requirements, respectively, for the memory, communication, and computation. All these requirements are reduced with an increasing number of processes, and the measured performance shows excellent parallel scalability and efficiency up to thousands of nodes. Additionally, the computationally more demanding quintic scaling steps can be accelerated by employing graphics processing units (GPU’s) showing, for large systems, a gain of almost a factor two compared to the standard central processing unit-only case. In this way, the evaluation of the derivatives of the RI-MP2 energy can be performed within a few minutes for systems containing hundreds of atoms and thousands of basis functions. With good time to solution, the implementation thus opens the possibility to perform molecular dynamics (MD) simulations in various ensembles (microcanonical ensemble and isobaric-isothermal ensemble) at the MP2 level of theory. Geometry optimization, full cell relaxation, and energy conserving MD simulations have been performed for a variety of molecular crystals including NH{sub 3}, CO{sub 2}, formic acid, and benzene.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗