Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FORMULATION”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Coherent states formulation of polymer field theory

We introduce a stable and efficient complex Langevin (CL) scheme to enable the first direct numerical simulations of the coherent-states (CS) formulation of polymer field theory. In contrast with Edwards’ well-known auxiliary-field (AF) framework, the CS formulation does not contain an embedded nonlinear, non-local, implicit functional of the auxiliary fields, and the action of the field theory has a fully explicit, semi-local, and finite-order polynomial character. In the context of a polymer solution model, we demonstrate that the new CS-CL dynamical scheme for sampling fluctuations in the space of coherent states yields results in good agreement with now-standard AF-CL simulations. The formalism is potentially applicable to a broad range of polymer architectures and may facilitate systematic generation of trial actions for use in coarse-graining and numerical renormalization-group studies.

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

Differentiable but exact formulation of density-functional theory

The universal density functional F of density-functional theory is a complicated and ill-behaved function of the density—in particular, F is not differentiable, making many formal manipulations more complicated. While F has been well characterized in terms of convex analysis as forming a conjugate pair (E, F) with the ground-state energy E via the Hohenberg–Kohn and Lieb variation principles, F is nondifferentiable and subdifferentiable only on a small (but dense) subset of its domain. In this article, we apply a tool from convex analysis, Moreau–Yosida regularization, to construct, for any ε > 0, pairs of conjugate functionals ({sup ε}E, {sup ε}F) that converge to (E, F) pointwise everywhere as ε → 0{sup +}, and such that {sup ε}F is (Fréchet) differentiable. For technical reasons, we limit our attention to molecular electronic systems in a finite but large box. It is noteworthy that no information is lost in the Moreau–Yosida regularization: the physical ground-state energy E(v) is exactly recoverable from the regularized ground-state energy {sup ε}E(v) in a simple way. All concepts and results pertaining to the original (E, F) pair have direct counterparts in results for ({sup ε}E, {sup ε}F). The Moreau–Yosida regularization therefore allows for an exact, differentiable formulation of density-functional theory. In particular, taking advantage of the differentiability of {sup ε}F, a rigorous formulation of Kohn–Sham theory is presented that does not suffer from the noninteracting representability problem in standard Kohn–Sham theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

The Lagrangian formulation of strong-field quantum electrodynamics in a plasma

The Lagrangian formulation of the scalar and spinor quantum electrodynamics in the presence of strong laser fields in a plasma medium is considered. We include the plasma influence in the free Lagrangian analogously to the “Furry picture” and obtain coupled equations of motion for the plasma particles and for the laser propagation. We demonstrate that the strong-field wave (i.e., the laser) satisfies a massive dispersion relation and obtain self-consistently the effective mass of the laser photons. The Lagrangian formulation derived in this paper is the basis for the cross sections calculation of quantum processes taking place in the presence of a plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Inclusion of line tension effect in classical nucleation theory for heterogeneous nucleation: A rigorous thermodynamic formulation and some unique conclusions

A rigorous thermodynamic formulation of the geometric model for heterogeneous nucleation including line tension effect is missing till date due to the associated mathematical hurdles. In this work, we develop a novel thermodynamic formulation based on Classical Nucleation Theory (CNT), which is supposed to illustrate a systematic and a more plausible analysis for the heterogeneous nucleation on a planar surface including the line tension effect. The appreciable range of the critical microscopic contact angle (θ{sub c}), obtained from the generalized Young’s equation and the stability analysis, is θ{sub ∞} < θ{sub c} < θ′ for positive line tension and is θ{sub M} < θ{sub c} < θ{sub ∞} for negative line tension. θ{sub ∞} is the macroscopic contact angle, θ′ is the contact angle for which the Helmholtz free energy has the minimum value for the positive line tension, and θ{sub M} is the local minima of the nondimensional line tension effect for the negative line tension. The shape factor f, which is basically the dimensionless critical free energy barrier, becomes higher for lower values of θ{sub ∞} and higher values of θ{sub c} for positive line tension. The combined effect due to the presence of the triple line and the interfacial areas (f{sup L} + f{sup S}) in shape factor is always within (0, 3.2), resulting f in the range of (0, 1.7) for positive line tension. A formerly presumed appreciable range for θ{sub c}(0 < θ{sub c} < θ{sub ∞}) is found not to be true when the effect of negative line tension is considered for CNT. Estimation based on the property values of some real fluids confirms the relevance of the present analysis.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

An exact formulation of the time-ordered exponential using path-sums

We present the path-sum formulation for the time-ordered exponential of a time-dependent matrix. The path-sum formulation gives the time-ordered exponential as a branched continued fraction of finite depth and breadth. The terms of the path-sum have an elementary interpretation as self-avoiding walks and self-avoiding polygons on a graph. Our result is based on a representation of the time-ordered exponential as the inverse of an operator, the mapping of this inverse to sums of walks on a graphs, and the algebraic structure of sets of walks. We give examples demonstrating our approach. We establish a super-exponential decay bound for the magnitude of the entries of the time-ordered exponential of sparse matrices. We give explicit results for matrices with commonly encountered sparse structures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Beyond the electric-dipole approximation: A formulation and implementation of molecular response theory for the description of absorption of electromagnetic field radiation

We present a formulation of molecular response theory for the description of a quantum mechanical molecular system in the presence of a weak, monochromatic, linearly polarized electromagnetic field without introducing truncated multipolar expansions. The presentation focuses on a description of linear absorption by adopting the energy-loss approach in combination with the complex polarization propagator formulation of response theory. Going beyond the electric-dipole approximation is essential whenever studying electric-dipole-forbidden transitions, and in general, non-dipolar effects become increasingly important when addressing spectroscopies involving higher-energy photons. These two aspects are examined by our study of the near K-edge X-ray absorption fine structure of the alkaline earth metals (Mg, Ca, Sr, Ba, and Ra) as well as the trans-polyenes. In following the series of alkaline earth metals, the sizes of non-dipolar effects are probed with respect to increasing photon energies and a detailed assessment of results is made in terms of studying the pertinent transition electron densities and in particular their spatial extension in comparison with the photon wavelength. Along the series of trans-polyenes, the sizes of non-dipolar effects are probed for X-ray spectroscopies on organic molecules with respect to the spatial extension of the chromophore.

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

A time-dependent formulation of multi-reference perturbation theory

Here, we discuss the time-dependent formulation of perturbation theory in the context of the interacting zeroth-order Hamiltonians that appear in multi-reference situations. As an example, we present a time-dependent formulation and implementation of second-order n-electron valence perturbation theory. The resulting time-dependent n-electron valence second-order perturbation theory (t-NEVPT2) method yields the fully uncontracted n-electron valence perturbation wavefunction and energy, but has a lower computational scaling than the usual contracted variants, and also avoids the construction of high-order density matrices and the diagonalization of metrics.We present results of t-NEVPT2 for the water, nitrogen, carbon, and chromium molecules and outline directions for the future.

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

A momentum-space formulation without partial wave decomposition for scattering of two spin-half particles

A new momentum-space formulation for scattering of two spin-half particles, both either identical or unidentical, is formulated. As basis states the free linear-momentum states are not expanded into the angular-momentum states, the system’s spin states are described by the product of the spin states of the two particles, and the system’s isospin states by the total isospin states of the two particles. We evaluate the Lippmann-Schwinger equations for the T-matrix elements in these basis states. The azimuthal behavior of the potential and of the T-matrix elements leads to a set of coupled integral equations for the T-matrix elements in two variables only, which are the magnitude of the relative momentum and the scattering angle. Some symmetry relations for the potential and the T-matrix elements reduce the number of the integral equations to be solved. A set of six spin operators to express any interaction of two spin-half particles is introduced. We show the spin-averaged differential cross section as being calculated in terms of the solution of the set of the integral equations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Characterization of detonation soot produced during steady and overdriven conditions for three high explosive formulations

The detonation of high explosives (HE) produces a dense fluid of molecular gases and solid carbon. The solid detonation carbon contains various carbon allotropes such as detonation nanodiamonds, onion-like carbon, graphite and amorphous carbon, with the formation of the different forms dependent upon pressure, temperature and the environmental conditions of the detonation. We have collected solid carbon residues from controlled detonations of three HE formulations (Composition B-3, PBX 9501, and PBX 9502). Soot was collected from experiments designed to produce both steady and overdriven conditions, and from detonations in both an ambient (air) atmosphere and in an inert Ar atmosphere. Differences in solid carbon residues were quantified using X-ray photoelectron spectroscopy and carbon isotope measurements. Environmental conditions, HE formulation, and peak pressures influenced the amount of and isotopic composition of the carbon in the soot. Detonations in an Ar atmosphere produced greater amounts of carbon soot with lower δ 13 C values than those in ambient air. Therefore, solid carbon residues continued to evolve after detonation due to excess oxygen in the ambient air detonations. As well, higher peak pressures in overdriven conditions produced less carbon soot with, in general, higher δ 13 C values. Consequently, while overdriven conditions only produced peak pressures for a limited duration, it was enough to influence the composition of the solid carbon residues.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A self-consistent field formulation of excited state mean field theory

In this work, we show that, as in Hartree-Fock theory, the orbitals for excited state mean field theory can be optimized via a self-consistent one-electron equation in which electron-electron repulsion is accounted for through mean field operators. In addition to showing that this excited state ansatz is sufficiently close to a mean field product state to admit a one-electron formulation, this approach brings the orbital optimization speed to within roughly a factor of two of ground state mean field theory. The approach parallels Hartree Fock theory in multiple ways, including the presence of a commutator condition, a one-electron mean-field working equation, and acceleration via direct inversion in the iterative subspace. When combined with a configuration interaction singles Davidson solver for the excitation coefficients, the self-consistent field formulation dramatically reduces the cost of the theory compared to previous approaches based on quasi-Newton descent.

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

Reduced scaling formulation of CASPT2 analytical gradients using the supporting subspace method

We present a reduced scaling and exact reformulation of state specific complete active space second-order perturbation (CASPT2) analytical gradients in terms of the MP2 and Fock derivatives using the supporting subspace method. This work follows naturally from the supporting subspace formulation of the CASPT2 energy in terms of the MP2 energy using dressed orbitals and Fock builds. For a given active space configuration, the terms corresponding to the MP2-gradient can be evaluated with O(N5) operations, while the rest of the calculations can be computed with O(N3) operations using Fock builds, Fock gradients, and linear algebra. When tensor-hyper-contraction is applied simultaneously, the computational cost can be further reduced to O(N4) for a fixed active space size. The new formulation enables efficient implementation of CASPT2 analytical gradients by leveraging the existing graphical processing unit (GPU)-based MP2 and Fock routines. We present benchmark results that demonstrate the accuracy and performance of the new method. Example applications of the new method in ab initio molecular dynamics simulation and constrained geometry optimization are given.

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

INFLUENCE OF SOLAR WIND HEATING FORMULATIONS ON THE PROPERTIES OF SHOCKS IN THE CORONA

One of the challenges in constructing global magnetohydrodynamic (MHD) models of the inner heliosphere for, e.g., space weather forecasting purposes, is to correctly capture the acceleration and expansion of the solar wind. In current models, various ad hoc heating prescriptions are introduced in order to obtain a realistic steady-state solar wind solution. In this work, we demonstrate, by performing MHD simulations of erupting coronal mass ejections (CMEs) on identical solar wind solutions employing different heating formulations, that the dynamics and properties of the CME-driven shocks are significantly altered depending on the applied heating prescription. Furthermore, we show how two popular heating formulations can be altered so as to yield shock properties consistent with theory and available coronal shock observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

FORMULATION OF NON-STEADY-STATE DUST FORMATION PROCESS IN ASTROPHYSICAL ENVIRONMENTS

The non-steady-state formation of small clusters and the growth of grains accompanied by chemical reactions are formulated under the consideration that the collision of key gas species (key molecule) controls the kinetics of dust formation process. The formula allows us to evaluate the size distribution and condensation efficiency of dust formed in astrophysical environments. We apply the formulation to the formation of C and MgSiO{sub 3} grains in the ejecta of supernovae, as an example, to investigate how the non-steady effect influences the formation process, condensation efficiency f{sub con,{sub ∞}}, and average radius a{sub ave,{sub ∞}} of newly formed grains in comparison with the results calculated with the steady-state nucleation rate. We show that the steady-state nucleation rate is a good approximation if the collision timescale of key molecule τ{sub coll} is much smaller than the timescale τ{sub sat} with which the supersaturation ratio increases; otherwise the effect of the non-steady state becomes remarkable, leading to a lower f{sub con,{sub ∞}} and a larger a{sub ave,{sub ∞}}. Examining the results of calculations, we reveal that the steady-state nucleation rate is applicable if the cooling gas satisfies Λ ≡ τ{sub sat}/τ{sub coll} ∼> 30 during the formation of dust, and find that f{sub con,{sub ∞}} and a{sub ave,{sub ∞}} are uniquely determined by Λ{sub on} at the onset time t{sub on} of dust formation. The approximation formulae for f{sub con,{sub ∞}} and a{sub ave,{sub ∞}} as a function of Λ{sub on} could be useful in estimating the mass and typical size of newly formed grains from observed or model-predicted physical properties not only in supernova ejecta but also in mass-loss winds from evolved stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Generalized massive gravity in arbitrary dimensions and its Hamiltonian formulation

We extend the four-dimensional de Rham-Gabadadze-Tolley (dRGT) massive gravity model to a general scalar massive-tensor theory in arbitrary dimensions, coupling a dRGT massive graviton to multiple scalars and allowing for generic kinetic and mass matrix mixing between the massive graviton and the scalars, and derive its Hamiltonian formulation and associated constraint system. When passing to the Hamiltonian formulation, two different sectors arise: a general sector and a special sector. Although obtained via different ways, there are two second class constraints in either of the two sectors, eliminating the BD ghost. However, for the special sector, there are still ghost instabilities except for the case of two dimensions. In particular, for the special sector with one scalar, there is a ''second BD ghost''.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Role of vertex corrections in the matrix formulation of the random phase approximation for the multiorbital Hubbard model

In the framework of a multiorbital Hubbard model description of superconductivity, a matrix formulation of the superconducting pairing interaction that has been widely used is designed to treat spin, charge, and orbital fluctuations within a random phase approximation (RPA). In terms of Feynman diagrams, this takes into account particle-hole ladder and bubble contributions as expected. It turns out, however, that this matrix formulation also generates additional terms which have the diagrammatic structure of vertex corrections. Furthermore we examine these terms and discuss the relationship between the matrix-RPA superconducting pairing interaction and the Feynman diagrams that it sums.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Cluster size convergence of the density matrix embedding theory and its dynamical cluster formulation: A study with an auxiliary-field quantum Monte Carlo solver

In this work, we investigate the cluster size convergence of the energy and observables using two forms of density matrix embedding theory (DMET): the original cluster form (CDMET) and a new formulation motivated by the dynamical cluster approximation (DCA-DMET). Both methods are applied to the half-filled one- and two-dimensional Hubbard models using a sign-problem free auxiliary-field quantum Monte Carlo impurity solver, which allows for the treatment of large impurity clusters of up to 100 sites. While CDMET is more accurate at smaller impurity cluster sizes, DCA-DMET exhibits faster asymptotic convergence towards the thermodynamic limit. We use our two formulations to produce new accurate estimates for the energy and local moment of the two-dimensional Hubbard model for U/t=2,4,6. These results compare favorably with the best data available in the literature, and help resolve earlier uncertainties in the moment for U/t=2.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Gauge invariance of radiative jet functions in the position-space formulation of SCET

In subleading powers of soft-collinear effective theory (SCET), the Lagrangian contains couplings between soft quarks and hard-collinear quarks. Matrix elements of the hard-collinear parts of these couplings are radiative jet functions. In the position-space formulation of SCET, the Lagrangians are constructed from operators that appear to be gauge invariant. Nevertheless, we find violations of gauge invariance arise in the hard-collinear sector because gauge transformations can shift the momentum of a hard-collinear quark field from the hard-collinear sector to the soft sector, where the hard-collinear fields, by definition, have no support. The violations of gauge invariance are manifested in perturbation theory in the hard-collinear sector through the absence of certain Feynman diagrams that would be present in full QCD. A consequence of the absence of these diagrams is that the radiative jet functions that follow directly from the position-space Lagrangians are not gauge invariant, and we demonstrate this through explicit calculations in lower-order perturbation theory. We obtain gauge-invariant Lagrangians by adding to existing position-space Lagrangians terms that are proportional to the soft-quark equation of motion. These gauge-invariant Lagrangians are valid for nonzero, as well as zero, quark masses. We also remark briefly on the gauge invariance of certain Lagrangians that have been constructed in the label-momentum formulation of SCET. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2010

Black hole entropy from the SU(2)-invariant formulation of type I isolated horizons

A detailed analysis of the spherically symmetric isolated horizon system is performed in terms of the connection formulation of general relativity. The system is shown to admit a manifestly SU(2) invariant formulation where the (effective) horizon degrees of freedom are described by an SU(2) Chern-Simons theory. This leads to a more transparent description of the quantum theory in the context of loop quantum gravity and modifications of the form of the horizon entropy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗