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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

PV Performance Modeling Methods and Practices: Results from the 4th PV Performance Modeling Collaborative Workshop.

In 2014, the IEA PVPS Task 13 added the PVPMC as a formal activity to its technical work plan for 2014-2017. The goal of this activity is to expand the reach of the PVPMC to a broader international audience and help to reduce PV performance modeling uncertainties worldwide. One of the main deliverables of this activity is to host one or more PVPMC workshops outside the US to foster more international participation within this collaborative group. This report reviews the results of the first in a series of these joint IEA PVPS Task 13/PVPMC workshops. The 4th PV Performance Modeling Collaborative Workshop was held in Cologne, Germany at the headquarters of TÜV Rheinland on October 22-23, 2015.

14 SOLAR ENERGY↗
DOE OSTI2017

New Methods in Non-Perturbative QCD

In this work, we investigate the properties of quantum chromodynamics (QCD), by using newly developing mathematics and physics formalisms. Almost all of the mass in the visible universe emerges from a quantum chromodynamics (QCD), which has a completely negligible microscopic mass content. An intimately related issue in QCD is the quark confinement problem. Answers to non-perturbative questions in QCD remained largely elusive despite much effort over the years. It is also believed that the usual perturbation theory is inadequate to address these kinds of problems. Perturbation theory gives a divergent asymptotic series (even when the theory is properly renormalized), and there are non-perturbative phenomena which never appear at any order in perturbation theory. Recently, a fascinating bridge between perturbation theory and non-perturbative effects has been found: a formalism called resurgence theory in mathematics tells us that perturbative data and non-perturbative data are intimately related. Translating this to the language of quantum field theory, it turns out that non-perturbative information is present in a coded form in perturbation theory and it can be decoded. We take advantage of this feature, which is particularly useful to understand some unresolved mysteries of QCD from first principles. In particular, we use: a) Circle compactifications which provide a semi-classical window to study confinement and mass gap problems, and calculable prototypes of the deconfinement phase transition; b) Resurgence theory and transseries which provide a unified framework for perturbative and non-perturbative expansion; c) Analytic continuation of path integrals and Lefschetz thimbles which may be useful to address sign problem in QCD at finite density.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

Displacement Fields and Self-Energies of Circular and Polygonal Dislocation Loops in Homogeneous and Layered Anisotropic Solids

There are large classes of materials problems that involve the solutions of stress, displacement, and strain energy of dislocation loops in elastically anisotropic solids, including increasingly detailed investigations of the generation and evolution of irradiation induced defect clusters ranging in sizes from the micro- to meso-scopic length scales. Based on a two-dimensional Fourier transform and Stroh formalism that are ideal for homogeneous and layered anisotropic solids, we have developed robust and computationally efficient methods to calculate the displacement fields for circular and polygonal dislocation loops. Using the homogeneous nature of the Green tensor of order -1, we have shown that the displacement and stress fields of dislocation loops can be obtained by numerical quadrature of a line integral. In addition, it is shown that the sextuple integrals associated with the strain energy of loops can be represented by the product of a pre-factor containing elastic anisotropy effects and a universal term that is singular and equal to that for elastic isotropic case. Furthermore, we have found that the self-energy pre-factor of prismatic loops is identical to the effective modulus of normal contact, and the pre-factor of shear loops differs from the effective indentation modulus in shear by only a few percent. These results provide a convenient method for examining dislocation reaction energetic and efficient procedures for numerical computation of local displacements and stresses of dislocation loops, both of which play integral roles in quantitative defect analyses within combined experimental–theoretical investigations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Multireference Equation-of-Motion Driven Similarity Renormalization Group: Theoretical Foundations and Applications to Ionized States

We present a formulation and implementation of an equation-of-motion (EOM) extension of the multireference driven similarity renormalization group (MR-DSRG) formalism for ionization potentials (IP-EOM-DSRG). The IP-EOM-DSRG formalism results in a Hermitian generalized eigenvalue problem, delivering accurate ionization potentials for strongly correlated systems. The EOM step scales as O(N 5 ) with the basis set size N, allowing for efficient calculation of spectroscopic properties, such as transition energies and intensities. The IP-EOM-DSRG formalism is combined with three truncation schemes of the parent MR-DSRG theory: an iterative nonperturbative method with up to two-body excitations [MR-LDSRG(2)] and second- and third-order perturbative approximations [DSRG-MRPT2/3]. We benchmark these variants by computing (1) the vertical valence ionization potentials of a series of small molecules at both equilibrium and stretched geometries; (2) the spectroscopic constants of several low-lying electronic states of the OH, CN, N 2 + , and CO + radicals; and (3) the binding curves of low-lying electronic states of the CN radical. A comparison with experimental data and theoretical results shows that all three IP-EOM-DSRG methods accurately reproduce the vertical ionization potentials and spectroscopic constants of these systems. Notably, the DSRG-MRPT3 and MR-LDSRG(2) versions outperform several state-of-the-art multireference methods of comparable or higher cost.

Hamiltonians↗
DOE OSTI2015

A hybrid framework of first principles molecular orbital calculations and a three-dimensional integral equation theory for molecular liquids: Multi-center molecular Ornstein–Zernike self-consistent field approach

In this study, we reported the development of a new quantum mechanics/molecular mechanics (QM/MM)-type framework to describe chemical processes in solution by combining standard molecular-orbital calculations with a three-dimensional formalism of integral equation theory for molecular liquids (multi-center molecular Ornstein–Zernike (MC-MOZ) method). The theoretical procedure is very similar to the 3D-reference interaction site model self-consistent field (RISM-SCF) approach. Since the MC-MOZ method is highly parallelized for computation, the present approach has the potential to be one of the most efficient procedures to treat chemical processes in solution. Benchmark tests to check the validity of this approach were performed for two solute (solute water and formaldehyde) systems and a simple S{sub N}2 reaction (Cl{sup −} + CH{sub 3}Cl → ClCH{sub 3} + Cl{sup −}) in aqueous solution. The results for solute molecular properties and solvation structures obtained by the present approach were in reasonable agreement with those obtained by other hybrid frameworks and experiments. In particular, the results of the proposed approach are in excellent agreements with those of 3D-RISM-SCF.

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

A full-potential approach to the relativistic single-site Green's function

One major purpose of studying the single-site scattering problem is to obtain the scattering matrices and differential equation solutions indispensable to multiple scattering theory (MST) calculations. On the other hand, the single-site scattering itself is also appealing because it reveals the physical environment experienced by electrons around the scattering center. In this study, we demonstrate a new formalism to calculate the relativistic full-potential single-site Green's function. We implement this method to calculate the single-site density of states and electron charge densities. Lastly, the code is rigorously tested and with the help of Krein's theorem, the relativistic effects and full potential effects in group V elements and noble metals are thoroughly investigated.

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

Unimodular F ( R ) gravity

We extend the formalism of the Einstein-Hilbert unimodular gravity in the context of modified F ( R ) gravity. After appropriately modifying the Friedmann-Robertson-Walker metric in a way that it becomes compatible to the unimodular condition of having a constant metric determinant, we derive the equations of motion of the unimodular F ( R ) gravity by using the metric formalism of modified gravity with Lagrange multiplier constraint. The resulting equations are studied in frames of reconstruction method, which enables us to realize various cosmological scenarios, which was impossible to realize in the standard Einstein-Hilbert unimodular gravity. Several unimodular F ( R ) inflationary scenarios are presented, and in some cases, concordance with Planck and BICEP2 observational data can be achieved.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

A comparative analysis of US state-level policies and programs to advance energy justice

Abstract The US energy system is undergoing massive changes that have environmental, technological, and societal implications. Decisions and actions taken now will dramatically alter the trajectory of our future energy system. Though new technologies can lessen the dangers of anthropogenic climate change, an effective and sustainable transition requires addressing social justice issues as a priority. To this end, many states have implemented policies and programs that address energy justice. A formal survey of state-level energy justice policies and programs is needed to better understand the methods used to assist communities targeted and the impacts of these efforts. Such an analysis should include an analysis of the metrics used to judge the effectiveness of the enacted policies and programs. In this paper, we analyze the energy justice policies and programs of seven different US states. From this analysis, we identify three important aspects of state energy policies. First, most policies and programs do not articulate their goals, nor cite metrics to judge their successes. This lack of transparency and accountability is a major roadblock to truly just solutions. Second, states focus on a very narrow range of solutions for energy justice issues. The most prevalent being utility bill assistance, which does not address the systemic and enduring challenges faced by many low-income or historically marginalized communities. Lastly, comprehensive energy policy does not acknowledge and address housing, employment, education, and healthcare inequities that exacerbate energy system inequities. The work presented in this paper sheds light on the progress of state-level energy justice policies and programs across the United States, and we argue more work needs to be done to understand best practices in energy justice policymaking.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2025

Dark matter substructure or source model systematics? A case study of cluster lens Abell S1063

Mapping the small-scale structure of the universe through gravitational lensing is a promising tool for probing the particle nature of dark matter. Curved Arc Basis (CAB) has been proposed as a local lensing formalism in galaxy clusters, with the potential to detect low-mass dark matter substructure. In this work, we analyse the cluster lens Abell S1063 in search of dark matter substructure with the CAB formalism, using multiband imaging data from James Webb Space Telescope ( JWST ). We use two different source modelling methods: shapelets and pixel-based source reconstruction based on Delaunay triangulation. We find that source modelling systematics from shapelets result in a disagreement between CAB parameters measured from different filters. Source modelling with Delaunay significantly alleviates this systematic, as seen in the improvement in agreement across filters. We also find that inadequate complexity in source modelling can result in convincing spurious detections of dark matter substructure from strong gravitational lenses, as seen by our $\Delta \text{BIC} > 20$ measurement of a $M \sim 10^{10}$ ${\rm M}_{\odot }$ subhalo with shapelets, a spurious detection that is not reproduced with Delaunay source modelling. We demonstrate that multiband analysis with different JWST filters is key for disentangling source and lens model systematics from dark matter substructure detections.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Real Space Quantum Cluster Formulation for the Typical Medium Theory of Anderson Localization

We develop a real space cluster extension of the typical medium theory (cluster-TMT) to study Anderson localization. By construction, the cluster-TMT approach is formally equivalent to the real space cluster extension of the dynamical mean field theory. Applying the developed method to the 3D Anderson model with a box disorder distribution, we demonstrate that cluster-TMT successfully captures the localization phenomena in all disorder regimes. As a function of the cluster size, our method obtains the correct critical disorder strength for the Anderson localization in 3D, and systematically recovers the re-entrance behavior of the mobility edge. From a general perspective, our developed methodology offers the potential to study Anderson localization at surfaces within quantum embedding theory. This opens the door to studying the interplay between topology and Anderson localization from first principles.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Auxiliary field Quantum Monte Carlo for dilute neutrons on the lattice

We employ constrained path Auxiliary Field Quantum Monte Carlo (AFQMC) in the pursuit of studying physical nuclear systems using a lattice formalism. Since AFQMC has been widely used in the study of condensed-matter systems such as the Hubbard model, we benchmark our method against published results for both one- and two-dimensional Hubbard model calculations. We then turn our attention to cold atomic and nuclear systems. We use an onsite contact interaction that can be tuned in order to reproduce the known scattering length and effective range of a given interaction. Developing this machinery allows us to extend our calculations to study nuclear systems within a lattice formalism. In conclusion, we perform initial calculations for a range of nuclear systems from two- to few-body neutron systems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

An Attractive Way to Correct for Missing Singles Excitations in Unitary Coupled Cluster Doubles Theory

Coupled cluster methods based exclusively on double excitations are comparatively “cheap” and interesting model chemistries, as they are typically able to capture the bulk of the dynamic electron correlation effects. The trade-off in such approximations is that the effect of neglected excitations, particularly single excitations, can be considerable. Using standard and electron-pair-restricted T 2 operators to define two flavors of unitary coupled cluster doubles (UCCD) methods, we investigate the extent to which missing single excitations can be recovered from low-order corrections in many-body perturbation theory (MBPT) within the unitary coupled cluster (UCC) formalism. Here, our analysis includes the derivations of finite-order UCC energy functionals, which are used as a basis to define perturbative estimates of missed single excitations. This leads to the novel UCCD[4S] and UCCD[6S] methods, which consider energy corrections for missing single excitations through fourth- and sixth-order in MBPT, respectively. We also apply the same methodology to the electron-pair-restricted ansatz, but the improvements are only marginal. Our findings show that augmenting UCCD with these post hoc perturbative corrections can lead to UCCSD-quality results.

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

Liquid crystal free energy relaxation by a theoretically informed Monte Carlo method using a finite element quadrature approach

A theoretically informed Monte Carlo method is proposed for Monte Carlo simulation of liquid crystals on the basis of theoretical representations in terms of coarse-grained free energy functionals. The free energy functional is described in the framework of the Landau-de Gennes formalism. A piecewise finite element discretization is used to approximate the alignment field, thereby providing an excellent geometrical representation of curved interfaces and accurate integration of the free energy. The method is suitable for situations where the free energy functional includes highly non-linear terms, including chirality or high-order deformation modes. Finally, the validity of the method is established by comparing the results of Monte Carlo simulations to traditional Ginzburg-Landau minimizations of the free energy using a finite difference scheme, and its usefulness is demonstrated in the context of simulations of chiral liquid crystal droplets with and without nanoparticle inclusions.

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

On the accuracy of the chemically significant eigenvalue method

We study the accuracy and convergence properties of the chemically significant eigenvalues method as proposed by Georgievskii et al. [J. Phys. Chem. A 117, 12146-12154 (2013)] and its close relative, dominant subspace truncation, for reduction of the energy-grained master equation. We formally derive the connection between both reduction techniques and provide hard error bounds for the accuracy of the latter which confirm the empirically excellent accuracy and convergence properties but also unveil practically relevant cases in which both methods are bound to fall short. We propose the use of balanced truncation as an effective alternative in these cases.

Chemistry↗
DOE OSTI2013

A simple method for one-loop renormalization in curved space-time

We present a simple method for deriving the renormalization counterterms from the components of the energy-momentum tensor in curved space-time. This method allows control over the finite parts of the counterterms and provides explicit expressions for each term separately. As an example, the method is used for the self-interacting scalar field in a Friedmann-Robertson-Walker metric in the adiabatic approximation, where we calculate the renormalized equation of motion for the field and the renormalized components of the energy-momentum tensor to fourth adiabatic order while including interactions to one-loop order. Within this formalism the trace anomaly, including contributions from interactions, is shown to have a simple derivation. We compare our results to those obtained by two standard methods, finding agreement with the Schwinger-DeWitt expansion but disagreement with adiabatic subtractions for interacting theories.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Multi-Hadron Observables from Lattice Quantum Chromodynamics

We describe formal work that relates the nite-volume spectrum in a quantum eld theory to scattering and decay amplitudes. This is of particular relevance to numerical calculations performed using Lattice Quantum Chromodynamics (LQCD). Correlators calculated using LQCD can only be determined on the Euclidean time axis. For this reason the standard method of determining scattering amplitudes via the Lehmann-Symanzik-Zimmermann reduction formula cannot be employed. By contrast, the nite-volume spectrum is directly accessible in LQCD calculations. Formalism for relating the spectrum to physical scattering observables is thus highly desirable. In this thesis we develop tools for extracting physical information from LQCD for four types of observables. First we analyze systems with multiple, strongly-coupled two-scalar channels. Here we accommodate both identical and nonidentical scalars, and in the latter case allow for degenerate as well as nondegenerate particle masses. Using relativistic eld theory, and summing to all orders in perturbation theory, we derive a result relating the nite-volume spectrum to the two-to-two scattering amplitudes of the coupled-channel theory. This generalizes the formalism of Martin L uscher for the case of single-channel scattering. Second we consider the weak decay of a single particle into multiple, coupled two-scalar channels. We show how the nite-volume matrix element extracted in LQCD is related to matrix elements of asymptotic two-particle states, and thus to decay amplitudes. This generalizes work by Laurent Lellouch and Martin L uscher. Third we extend the method for extracting matrix elements by considering currents which insert energy, momentum and angular momentum. This allows one to extract transition matrix elements and form factors from LQCD. Finally we look beyond two-particle systems to those with three-particles in asymptotic states. Working again to all orders in relativistic eld theory, we derive a relation between the spectrum and an in nite-volume three-to-three scattering quantity. This nal analysis is the most complicated of the four, because the all-orders summation is more di cult for this system, and also because a number of new technical issues arise in analyzing the contributing diagrams.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

Time-sliced perturbation theory for large scale structure I: general formalism

We present a new analytic approach to describe large scale structure formation in the mildly non-linear regime. The central object of the method is the time-dependent probability distribution function generating correlators of the cosmological observables at a given moment of time. Expanding the distribution function around the Gaussian weight we formulate a perturbative technique to calculate non-linear corrections to cosmological correlators, similar to the diagrammatic expansion in a three-dimensional Euclidean quantum field theory, with time playing the role of an external parameter. For the physically relevant case of cold dark matter in an Einstein-de Sitter universe, the time evolution of the distribution function can be found exactly and is encapsulated by a time-dependent coupling constant controlling the perturbative expansion. We show that all building blocks of the expansion are free from spurious infrared enhanced contributions that plague the standard cosmological perturbation theory. This paves the way towards the systematic resummation of infrared effects in large scale structure formation. We also argue that the approach proposed here provides a natural framework to account for the influence of short-scale dynamics on larger scales along the lines of effective field theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

Analytic derivative couplings in time-dependent density functional theory: Quadratic response theory versus pseudo-wavefunction approach

We revisit the formalism for analytic derivative couplings between excited states in time-dependent density functional theory (TDDFT). We derive and implement these couplings using quadratic response theory, then numerically compare this response-theory formulation to couplings implemented previously based on a pseudo-wavefunction formalism and direct differentiation of the Kohn-Sham determinant. Numerical results, including comparison to full configuration interaction calculations, suggest that the two approaches perform equally well for many molecular systems, provided that the underlying DFT method affords accurate potential energy surfaces. The response contributions are found to be important for certain systems with high symmetry, but can be calculated with only a moderate increase in computational cost beyond what is required for the pseudo-wavefunction approach. In the case of spin-flip TDDFT, we provide a formal proof that the derivative couplings obtained using response theory are identical to those obtained from the pseudo-wavefunction formulation, which validates our previous implementation based on the latter formalism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗