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

Exploring Applied Cryptosystems to Formally Verify Security in Cyber-Physical Systems

This project aims to evaluate RSA as a method for public-key encryption for cyber-physical systems (CPS). As technology advances, cyber attacks are increasing, and with them, the need for cybersecurity advances; the average cost for cybercrime in the world was estimated at $6 trillion in 2021. A public-key cryptosystem that has been around since 1977, RSA has recently garnered some critiques for its fragility, computational cost, and lazy implementation. In this project I will review the mathematical derivation of RSA, analyze the practical implications of such mathematical framework for the security of RSA, and propose a formal methods based approach to verify encryption schemes for CPS.

97 MATHEMATICS AND COMPUTING↗

Gamow shell model description of proton scattering on 18 Ne

Background: The structure of weakly bound/unbound nuclei close to particle drip lines is different from that around the valley of beta stability. A comprehensive description of these systems goes beyond the standard shell model (SM) and demands an open quantum system description of the nuclear many-body system. Purpose: For that purpose, we are using the Gamow shell model (GSM) in this study, which provides a fully microscopic description of bound and unbound nuclear states, nuclear decays, and reactions. Here, we formulate the GSM in coupled-channel (GSM-CC) representation to describe low-energy elastic and inelastic scattering of protons on $\mathrm{^{18}Ne}$. Method: The GSM-CC formalism is applied to a translationally invariant Hamiltonian with an effective finite-range two-body interaction. We discuss in detail the GSM-CC formalism in coordinate space and give the description of the novel equivalent potential method for solving the GSM-CC system of integrodifferential equations. This method is then applied for the description of ($p, p'$) reaction cross-sections. The reactions channels are built by GSM wave functions for the ground state $\mathrm{0^+}$ and the first excited $\mathrm{2^+}$ of $\mathrm{^{18}Ne}$ and a proton wave function expanded in different partial waves. The completeness of this basis is verified by comparing GSM and GSM-CC energies of low-energy resonant states in $\mathrm{^{19}Na}$. The differences between the two calculations provide a measure of the missing configurations in the GSM-CC calculation of low-energy states of $\mathrm{^{19}Na}$ due to the restriction on the number of excited states of $\mathrm{^{18}Ne}$. Results: We present the first application of the GSM-CC formalism for the calculation of excited states of $\mathrm{^{18}Ne}$ and $\mathrm{^{19}Na}$, the excitation function, and the elastic/inelastic differential cross-sections in the $\mathrm{^{18}Ne}$ ($p,p'$) reaction at different energies. This is the first unified description of the spectra and reaction cross-sections in the GSM formalism. The method is shown to be both feasible and accurate. The approximate equivalence of GSM and GSM-CC in describing spectra of $\mathrm{^{19}Na}$ has been demonstrated numerically. Conclusions: The GSM in the coupled-channel representation opens a possibility for the unified description of low-energy nuclear structure and reactions using the same Hamiltonian. While both GSM and GSM-CC can describe energies, widths, and wave functions of the many-body states, the GSM-CC can in addition yield reaction cross-sections. The combined application of GSM and GSM-CC to describe energies of resonant states allows to test the exactitude of calculated cross-sections for a given many-body Hamiltonian.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Selecting Essential Information for Biosurveillance—A Multi-Criteria Decision Analysis

The National Strategy for Biosurveillancedefines biosurveillance as ‘‘the process of gathering, integrating, interpreting, and communicating essential information related to all-hazards threats or disease activity affecting human, animal, or plant health to achieve early detection and warning, contribute to overall situational awareness of the health aspects of an incident, and to enable better decision-making at all levels.’’ However, the strategy does not specify how ‘‘essential information’’ is to be identified and integrated into the current biosurveillance enterprise, or what the metrics qualify information as being ‘‘essential’’. Thequestion of data stream identification and selection requires a structured methodology that can systematically evaluate the tradeoffs between the many criteria that need to be taken in account. Multi-Attribute Utility Theory, a type of multi-criteria decision analysis, can provide a well-defined, structured approach that can offer solutions to this problem. While the use of Multi-Attribute Utility Theoryas a practical method to apply formal scientific decision theoretical approaches to complex, multi-criteria problems has been demonstrated in a variety of fields, this method has never been applied to decision support in biosurveillance.We have developed a formalized decision support analytic framework that can facilitate identification of ‘‘essential information’’ for use in biosurveillance systems or processes and we offer this framework to the global BSV community as a tool for optimizing the BSV enterprise. To demonstrate utility, we applied the framework to the problem of evaluating data streams for use in an integrated global infectious disease surveillance system.

47 OTHER INSTRUMENTATION↗

Air-kerma strength determination of a new directional {sup 103}Pd source

Purpose: A new directional {sup 103}Pd planar source array called a CivaSheet™ has been developed by CivaTech Oncology, Inc., for potential use in low-dose-rate (LDR) brachytherapy treatments. The array consists of multiple individual polymer capsules called CivaDots, containing {sup 103}Pd and a gold shield that attenuates the radiation on one side, thus defining a hot and cold side. This novel source requires new methods to establish a source strength metric. The presence of gold material in such close proximity to the active {sup 103}Pd region causes the source spectrum to be significantly different than the energy spectra of seeds normally used in LDR brachytherapy treatments. In this investigation, the authors perform air-kerma strength (S{sub K}) measurements, develop new correction factors for these measurements based on an experimentally verified energy spectrum, and test the robustness of transferring S{sub K} to a well-type ionization chamber. Methods: S{sub K} measurements were performed with the variable-aperture free-air chamber (VAFAC) at the University of Wisconsin Medical Radiation Research Center. Subsequent measurements were then performed in a well-type ionization chamber. To realize the quantity S{sub K} from a directional source with gold material present, new methods and correction factors were considered. Updated correction factors were calculated using the MCNP 6 Monte Carlo code in order to determine S{sub K} with the presence of gold fluorescent energy lines. In addition to S{sub K} measurements, a low-energy high-purity germanium (HPGe) detector was used to experimentally verify the calculated spectrum, a sodium iodide (NaI) scintillating counter was used to verify the azimuthal and polar anisotropy, and a well-type ionization chamber was used to test the feasibility of disseminating S{sub K} values for a directional source within a cylindrically symmetric measurement volume. Results: The UW VAFAC was successfully used to measure the S{sub K} of four CivaDots with reproducibilities within 0.3%. Monte Carlo methods were used to calculate the UW VAFAC correction factors and the calculated spectrum emitted from a CivaDot was experimentally verified with HPGe detector measurements. The well-type ionization chamber showed minimal variation in response (<1.5%) as a function of source positioning angle, indicating that an American Association of Physicists in Medicine (AAPM) Accredited Dosimetry Calibration Laboratory calibrated well chamber would be a suitable device to transfer an S{sub K}-based calibration to a clinical user. S{sub K} per well-chamber ionization current ratios were consistent among the four dots measured. Additionally, the measurements and predictions of anisotropy show uniform emission within the solid angle of the VAFAC, which demonstrates the robustness of the S{sub K} measurement approach. Conclusions: This characterization of a new {sup 103}Pd directional brachytherapy source helps to establish calibration methods that could ultimately be used in the well-established AAPM Task Group 43 formalism. Monte Carlo methods accurately predict the changes in the energy spectrum caused by the fluorescent x-rays produced in the gold shield.

07 ISOTOPE AND RADIATION SOURCES↗

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↗

Theory and Applications of Macroscale Models in Porous Media

Systems dominated by heterogeneity over a multiplicity of scales, like porous media, still challenge our modeling efforts. The presence of disparate length- and time-scales that control dynamical processes in porous media hinders not only models predictive capabilities, but also their computational efficiency. Macrosopic models, i.e., averaged representations of pore-scale processes, are computationally efficient alternatives to microscale models in the study of transport phenomena in porous media at the system, field or device scale (i.e., at a scale much larger than a characteristic pore size). We present an overview of common upscaling methods used to formally derive macroscale equations from pore-scale (mass, momentum and energy) conservation laws. This review includes the volume averaging method, mixture theory, thermodynamically constrained averaging, homogenization, and renormalization group techniques. We apply these methods to a number of specific problems ranging from food processing to human bronchial system, and from diffusion to multiphase flow, to demonstrate the methods generality and flexibility in handling different applications. The primary intent of such an overview is not to provide a thorough review of all currently available upscaling techniques, nor a complete mathematical treatment of the ones presented, but rather a primer on some of the tools available for upscaling, the basic principles they are based upon, and their specific advantages and drawbacks, so to guide the reader in the choice of the most appropriate method for particular applications and of the most relevant technical literature.

42 ENGINEERING↗

Description of ground and excited electronic states by ensemble density functional method with extended active space

An extended variant of the spin-restricted ensemble-referenced Kohn-Sham (REKS) method, the REKS(4,4) method, designed to describe the ground electronic states of strongly multireference systems is modified to enable calculation of excited states within the time-independent variational formalism. The new method, the state-interaction state-averaged REKS(4,4), i.e., SI-SA-REKS(4,4), is capable of describing several excited states of a molecule involving double bond cleavage, polyradical character, or multiple chromophoric units.We demonstrate that the newmethod correctly describes the ground and the lowest singlet excited states of a molecule (ethylene) undergoing double bond cleavage. The applicability of the new method for excitonic states is illustrated with π stacked ethylene and tetracene dimers. We conclude that the new method can describe a wide range of multireference phenomena.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The optical potential on the lattice

The extraction of hadron-hadron scattering parameters from lattice data by using the Luscher approach becomes increasingly complicated in the presence of inelastic channels. We propose a method for the direct extraction of the complex hadron-hadron optical potential on the lattice, which does not require the use of the multi-channel Luscher formalism. Furthermore, this method is applicable without modifications if some inelastic channels contain three or more particles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Hyper Open-Shell Excited Spin States of Transition-Metal Compounds: FeF 2 , FeF 2 ···Ethane, and FeF 2 ···Ethylene

Spin-state energetics are important for understanding properties that involve more than one spin state, for example, catalysis occurring on two or more potential energy surfaces corresponding to different electronic spins. Very often, multiple-spin processes involve transition-metal compounds, and therefore, it is important to understand the electronic structure and energetics of such compounds in different spin states. In this work, we benchmark relative spin-state energies of FeF 2 with respect to the quintet ground spin state using both singleconfigurational and multi-configurational methods, and we examine how they are affected by the binding of ethane and ethylene to the iron center. We also benchmark the binding energies of the complexes. The single-configurational methods used in this work are the Hartree–Fock method, 32 exchange–correlation functionals, and the CCSD(T) coupled-cluster method in both restricted and unrestricted formalisms. The multi-configurational methods that have been used are CASSCF, CASPT2, CASPT3, MRCI, MRCI+Q, and MR-ACPF. The spin-state splitting energies depend on the functional chosen and of the 32 exchange-correlation functionals investigated here we find that for the septet and spin-projected triplet states of FeF 2 , the M06 functional is the best when compared to our best estimates from multi-reference calculations. If all the nine excitation energies are considered, where there are three excited spin states (singlet, triplet, and septet) for each of the three systems (FeF 2 , FeF 2 ···ethane, and FeF 2 ···ethylene), the three best performing functionals are HLE16, SOGGA11-X, and M06-2X. We find that the binding of ethane perturbs the relative spin-state energy of FeF 2 by only a small amount, but the stronger binding of ethylene has a larger effect. For the spin-state splitting energies of FeF 2 using single-reference CCSD(T), we found the predicted results depend very strongly on precisely how the calculations are done, in particular, on the spin-restricted or spin-unrestricted character of the SCF reference state, which can differ even by around 50 kcal/mol for the SCF reference state and the subsequent CCSD(T) calculations. Furthermore, on analyzing the wave functions of both the spinrestricted or spin-unrestricted formalisms, we find that the lowest-energy singlet and triplet states of the complexes, just like FeF 2 in isolation, can have more unpaired electrons than they are usually assumed to have, i.e. they can be hyper open-shell electronic configurations, and that this can significantly lower the energy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Degree of rate control and De Donder relations – An interpretation based on transition state theory

We describe sensitivities of elementary steps as salient parameters in determining the rate determining character of elementary steps within a reaction network and develop a formalism wherein the overall composite reaction is described by an apparent rate-determining step that is a sensitivity-weighted average of the elementary steps that comprise the reaction. Reaction parameters—apparent reaction orders, apparent enthalpy and entropy of activation—of the composite reaction network are determined within the framework we develop by application of transition state theory to the apparent rate-determining step. From this formalism we develop methods for determination of surface coverages by measuring only reaction orders, novel methods for discrimination between proposed mechanisms, a proof that the kinetic degrees of rate control sum to unity, and a proof for the relationship between fractional coverages and thermodynamic degrees of rate control. Two complementary formalisms for identifying rate-limiting transition states are broadly employed in chemical kinetics—De Donder relations based on measuring the thermodynamic driving forces of elementary steps during reaction and degrees of rate control based on knowledge of the kinetics of elementary reaction steps. Here, the formalism developed herein unifies these two strategies to elucidate in the most general case a relationship between the thermodynamic driving forces of elementary steps captured by the reversibilities and the kinetic and thermodynamic degrees of rate control.

42 ENGINEERING↗

On the Question of the Total Energy in the Fermi–Löwdin Orbital Self-Interaction Correction Method

The Fermi-Löwdin orbital self-interaction correction (FLOSIC) formalism is a novel method for implementing the Perdew-Zunger self-interaction correction (PZ-SIC) in density functional theory calculations. In this paper we consider how the use of Fermi orbitals affects total energies and other calculated properties compared to a standard approach to PZ-SIC that utilizes the localization equation conditions. Here, we directly compare the results of the two methods using identical basis sets and numerical techniques in calculations for isolated atoms up to Kr and for a large test set of molecules. We find differences in total energies that increase with increasing atomic number and show that these differences can be traced to a less negative SIC correction for the 1s orbital in FLOSIC. Importantly, energies for highest occupied orbitals and molecular atomization energies are nearly identical in the two methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Fast nonadiabatic dynamics of many-body quantum systems

Modeling many-body quantum systems with strong interactions is one of the core challenges of modern physics. A range of methods has been developed to approach this task, each with its own idiosyncrasies, approximations, and realm of applicability. However, there remain many problems that are intractable for existing methods. In particular, many approaches face a huge computational barrier when modeling large numbers of coupled electrons and ions at finite temperature. Here, we address this shortfall with a new approach to modeling many-body quantum systems. On the basis of the Bohmian trajectory formalism, our new method treats the full particle dynamics with a considerable increase in computational speed. As a result, we are able to perform large-scale simulations of coupled electron-ion systems without using the adiabatic Born-Oppenheimer approximation.

Science & Technology - Other Topics↗

On the feasibility of using the spiral beam formalism for analysis of cardiograms

A contour analysis method based on a spiral beam formalism is proposed for cardiogram classification. A cardiogram contour proximity metric is introduced in spiral beam intensity space. Normal and infarction cardiograms are classified using the proposed method. The method is shown to provide adequate results in most of the cases examined. (paper)

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Zero modes in the light-front coupled-cluster method

The light-front coupled-cluster (LFCC) method is a technique for solving Hamiltonian eigenvalue problems in light-front-quantized field theories. Its primary purpose is to provide a systematic sequence of solvable approximations to the original eigenvalue problem without the truncation of Fock space. Here we discuss the incorporation of zero modes, modes of zero longitudinal momentum, into the formalism of the method. Without zero modes, the light-front vacuum is trivial, and the vacuum expectation value of the field is always zero. The LFCC method with zero modes provides for vacuum structure, in the form of a generalized coherent state of zero modes, as is illustrated here in two-dimensional model field theories. -- Highlights: •Extends the light-front coupled-cluster method to include zero modes. •Illustrates with an analysis of vacuum structure for phi-3, phi-4, and Wick–Cutkosky model field theories. •Demonstrates the applicability of the LFCC method to theories with spontaneous symmetry breaking.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Photon topology

The topology of photons in vacuum is interesting because there are no photons with k = 0, creating a hole in momentum space. We show that while the set of all photons forms a trivial vector bundle $γ$ over this momentum space, the R and L photons form topologically nontrivial subbundles $γ±$ with first Chern numbers ∓2. In contrast, $γ$ has no linearly polarized subbundles, and there is no Chern number associated with linear polarizations. It is a known difficulty that the standard version of Wigner’s little group method produces singular representations of the Poincaré group for massless particles. By considering representations of the Poincaré group on vector bundles we obtain a version of Wigner’s little group method for massless particles which avoids these singularities. Here we show that any massless bundle representation of the Poincaré group can be canonically decomposed into irreducible bundle representations labeled by helicity, which in turn can be associated to smooth irreducible Hilbert space representations. This proves that the R and L photons are globally well defined as particles and that the photon wave function can be uniquely split into R and L components. This formalism offers a method of quantizing the electromagnetic field without invoking discontinuous polarization vectors as in the traditional scheme. We also demonstrate that the spin-Chern number of photons is not a purely topological quantity. Lastly, there has been an extended debate on whether photon angular momentum can be split into spin and orbital parts. Our work explains the precise issues that prevent this splitting. Photons do not admit a spin operator; instead, the angular momentum associated with photons’ internal degree of freedom is described by a helicity-induced subalgebra corresponding to the translational symmetry of $γ$.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Ultrafast disruptive probing: Simultaneously keeping track of tens of reaction pathways

Ultrafast science depends on different implementations of the well-known pump–probe method. Here, we provide a formal description of ultrafast disruptive probing, a method in which the probe pulse disrupts a transient species that may be a metastable ion or a transient state of matter. Disruptive probing has the advantage of allowing for simultaneous tracking of the yield of tens of different processes. Our presentation includes a numerical model and experimental data on multiple products resulting from the strong-field ionization of two different molecules, partially deuterated methanol and norbornene. The correlated enhancement and depletion signals between all the different fragmentation channels offer comprehensive information on photochemical reaction pathways. In combination with ion imaging and/or coincidence momentum imaging or as complementary to atom-specific probing or ultrafast diffraction methods, disruptive probing is a particularly powerful tool for the study of strong-field laser–matter interactions.

47 OTHER INSTRUMENTATION↗

A deterministic solution of the first order linear Boltzmann transport equation in the presence of external magnetic fields

Purpose: Accurate radiotherapy dose calculation algorithms are essential to any successful radiotherapy program, considering the high level of dose conformity and modulation in many of today’s treatment plans. As technology continues to progress, such as is the case with novel MRI-guided radiotherapy systems, the necessity for dose calculation algorithms to accurately predict delivered dose in increasingly challenging scenarios is vital. To this end, a novel deterministic solution has been developed to the first order linear Boltzmann transport equation which accurately calculates x-ray based radiotherapy doses in the presence of magnetic fields. Methods: The deterministic formalism discussed here with the inclusion of magnetic fields is outlined mathematically using a discrete ordinates angular discretization in an attempt to leverage existing deterministic codes. It is compared against the EGSnrc Monte Carlo code, utilizing the emf-macros addition which calculates the effects of electromagnetic fields. This comparison is performed in an inhomogeneous phantom that was designed to present a challenging calculation for deterministic calculations in 0, 0.6, and 3 T magnetic fields oriented parallel and perpendicular to the radiation beam. The accuracy of the formalism discussed here against Monte Carlo was evaluated with a gamma comparison using a standard 2%/2 mm and a more stringent 1%/1 mm criterion for a standard reference 10 × 10 cm{sup 2} field as well as a smaller 2 × 2 cm{sup 2} field. Results: Greater than 99.8% (94.8%) of all points analyzed passed a 2%/2 mm (1%/1 mm) gamma criterion for all magnetic field strengths and orientations investigated. All dosimetric changes resulting from the inclusion of magnetic fields were accurately calculated using the deterministic formalism. However, despite the algorithm’s high degree of accuracy, it is noticed that this formalism was not unconditionally stable using a discrete ordinate angular discretization. Conclusions: The feasibility of including magnetic field effects in a deterministic solution to the first order linear Boltzmann transport equation is shown. The results show a high degree of accuracy when compared against Monte Carlo calculations in all magnetic field strengths and orientations tested.

60 APPLIED LIFE SCIENCES↗

SU-F-T-16: Experimental Determination of Ionization Chamber Correction Factors for In-Phantom Measurements of Reference Air Kerma Rate and Absorbed Water Dose Rate of Brachytherapy 192Ir Source

Purpose: Following the method of in-phantom measurements of reference air kerma rate (Ka) at 100cm and absorbed water dose rate (Dw1) at 1cm of high-dose-rate 192Ir brachytherapy source using 60Co absorbed-dose-to-water calibrated (ND,w,60Co) ionization chamber (IC), we experimentally determined the in-phantom correction factors (kglob) of the PTW30013 (PTW, Freiburg, Germany) IC by comparing the Monte Carlo (MC)-calculated kglob of the other PTW30016 IC. Methods: The Dw1 formalism of in-phantom measurement is: M*ND,w,60Co*(kglob)Dw1, where M is the collected charges, and (kglob)Dw1 the in-phantom Dw1 correction factor. Similarly, Ka is determined by M*ND,w,60Co*(kglob)ka, where (kglob)ka the in-phantom Ka correction factor. Two thimble ICs PTW30013 and another PTW30016 having a ND,w,60Co from the German primary standard laboratory (PTB) were simultaneously exposed to the microselectron 192Ir v2 source at 8cm in a PMMA phantom. A reference well chamber (PTW33004) with a PTB transfer Ka calibration Nka was used for comparing the in-phantom measurements to derive the experimental (kglob)ka factors. We determined the experimental (kglob)Dw1 of the PTW30013 by comparing the PTW30016 measurements with MC-calculated (kglob)Dw1. Results: Ka results of the PTW30016 based on ND,w,60Co and MC-calculated (kglob)ka differ from the well chamber results based on Nka by 1.6% and from the manufacturer by 1.0%. Experimental (kglob)ka factors for the PTW30016 and two other PTW30013 are 0.00683, 0.00681 and 0.00679, and vary <0.5% with 1mm source positioning uncertainty. Experimental (kglob)Dw1 of the PTW30013 ICs are 75.3 and 75.6, and differ by 1.6% from the conversion by dose rate constant from the AAPM report 229. Conclusion: The 1.7% difference between MC and experimental (kglob)ka for the PTW30016 IC is within the PTB 2.5% expanded uncertainty in Ka calibration standard. Using a single IC with ND,w,60Co to calibrate the brachytherapy source and dose output in external radiotherapy is feasible. MC validation of the PTW30013(kglob)Dw1 is warranted.

60 APPLIED LIFE SCIENCES↗