Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “theory - 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 937 records · Page 52

Crystal growth, density functional theory, and scintillation properties of Tl 3 LnCl 6 :Ce 3+ and TlLn 2 Cl 7 :Ce 3+ (Ln = Y, Gd)

In this paper we report on the crystal growth, density functional theory (DFT) calculations and scintillation properties of Tl 3 LnCl 6 :Ce and TlLn 2 Cl 7 :Ce(Ln = Y, Gd). Crystals were grown by the Vertical Bridgman technique up to 16 mm in diameter and 25 mm long. Crystals of Tl 3 YCl 6 :Ce and Tl 3 GdCl 6 :Ce belong to the family of the Nesohalides which have the monoclinic crystal structure. Crystals of TlY 2 Cl 7 :Ce and TlGd 2 Cl 7 :Ce belong to the family of the Phyllohalides which have either the monoclinic or the orthorhombic crystal structure. We report the light yields of these Nesohalides and Phyllohalides are typically on the order of 35,000–40,000 ph/MeV.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

A layered 2D computational framework: Theory and applications to nuclear fuel behavior

Nuclear fuel performance computer codes have been developed over the last 50 years to analyze fuel behavior under various operating conditions. Traditionally, these codes used quasi-two-dimensional (also commonly known as 1.5D) representations of the fuel rod, which model the rod using a set of one-dimensional axisymmetric models that represent the behavior at specific axial positions on the rod. Modern fuel performance codes have the ability to investigate full three-dimensional (3D) effects and couple to other physics-based codes for true multiphysics simulations. However, with increasing complexity comes increasing computational costs. Many phenomena of interest involve azimuthally-varying behavior that cannot be represented using the aforementioned quasi-two-dimensional approach, but do not require the use of a full 3D model. To efficiently address these problems, there is a need for a computational framework that provides a compromise between the quasi-two-dimensional and full 3D models. In this paper, we present a new quasi-three-dimensional approach that represents the fuel as a set of 2D planar models that represent the behavior of the fuel cross-section at various axial positions. Presented here are the theory behind the methodology, test cases to illustrate proper implementation, and practical applications of its use in the BISON fuel performance code for a variety of cases in nuclear fuel analysis, such as fuel fracture, axial fuel relocation, and cladding distension and oxidation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Studying QGP with flow: A theory overview

I review recent developments in the phenomenological study of the quark-gluon plasma (QGP) transport properties based on a personal selection of results that were presented at Quark Matter 2019. Additionally, the constraints on the temperature dependence of QGP shear and bulk viscosity are summarized. I discuss new theory advancements towards more realistic 3D dynamical simulations of heavy-ion collisions at finite baryon density. The challenges and opportunities of applying hydrodynamics to small collision systems are highlighted.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

On the canonical formulation of gauge field theories and Poincaré transformations

We address the Hamiltonian formulation of classical gauge field theories while putting forward results some of which are not entirely new, though they do not appear to be well known. We refer in particular to the fact that neither the canonical energy momentum vector ( P μ ) nor the gauge invariant energy momentum vector ( P inv μ ) do generate space-time translations of the gauge field by means of the Poisson brackets: In a general gauge, one has to consider the so-called kinematical energy momentum vector and, in a specific gauge (like the radiation gauge in electrodynamics), one has to consider the Dirac brackets rather than the Poisson brackets. Similar arguments apply to rotations and to Lorentz boosts and are of direct relevance to the “nucleon spin crisis” since the spin of the proton involves a contribution which is due to the angular momentum vector of gluons and thereby requires a proper treatment of the latter. We conclude with some comments on the relationships between the different approaches to quantization (canonical quantization based on the classical Hamiltonian formulation, Gupta-Bleuler, path integrals, BRST, covariant canonical approaches).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Matching the Standard Model to heavy-quark effective theory and nonrelativistic QCD

We find the leading electroweak corrections to the Lagrangians of heavy-quark effective theory and nonrelativistic QCD. These corrections appear in the Wilson coefficients of the two- and four-quark operators and are considered here at one-loop order through $\mathscr{O}(1/m^3)$ and $\mathscr{O}(1/m^2)$, respectively. The two-quark operators through this order include new parity violating terms, which we derive analogously to the parity preserving QCD result at one-loop order.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Effect of single atom Platinum (Pt) doping and facet dependent on the electronic structure and light absorption of Lanthanum Titanium Oxide (La 2 Ti 2 O 7 ): A Density Functional Theory study

Charge generation and separation are two key processes for semiconductor photocatalysis. Here, we use Pt as a single atom catalyst to systematically examine the facet-dependent electronic band structure and light absorption of the layered perovskite-type wide-gap semiconductor, lanthanum titanium oxide (La 2 Ti 2 O 7 , LTO) by means of density functional theory simulations. In this study, it is found that single Pt atom doping of different LTO surfaces (here, (100), (101) and (001)) can not only create states in the bandgap that would promote the formation of recombination centers, but also shift the optical absorption edge to the visible region. Interestingly, the Pt doping forms a heterojunction, with the valence band maximum consisting of Pt states and the conduction band minimum consisting of LTO states, respectively. To study the facet dependent surface activity, acetic acid (CH 3 COOH) was used as a model molecule to investigate the adsorption and charge transfer on the (101), (100) and (001) Pt-LTO surface facets. The results show that the (101) facet could enable stronger adsorption of CH 3 COOH by promoting more electron transfer during the interfacial interaction. Our theoretical findings aim to promote the design and optimization of the single atom catalysts (SACs) for photocatalytic applications and other broad catalysis systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Status and Challenges of Density Functional Theory

Here, we discuss some of the challenges facing density functional theory (DFT) and recent progress in DFT for both ground and excited electronic states. We discuss key aspects of the results we have been able to obtain with the strategy of designing density functionals to have various ingredients and functional forms that are then optimized to accurately predict various types of properties and systems with as much universality as possible. Finally, we make specific recommendations of approximate density functionals that are well suited for particular kinds of applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Specific resonance mode enhancement and suppression using non-uniform polarization of piezoelectric wafers: Theory and experiments

In this work, we present experimental demonstration of specific resonance mode enhancement and suppression in circular and rectangular piezoelectric wafers with engineered non-uniform polarization profiles. The polarization profiles are designed based on the electromechanical impedance response of non-uniformly polarized wafers as obtained from the theory. The circular wafers are fabricated with non-uniform polarization profiles that involve a central polarized region surrounded by an unpolarized region. The radius of the polarization zone is designed based on the condition for specific mode enhancement obtained from the electromechanical impedance response of non-uniformly polarized wafers. We actually show how the condition can not only be used to enhance but also to suppress electromechanical resonances. Two kinds of wafers are designed and fabricated to specifically suppress second and third radial modes respectively. Similarly, rectangular wafers are designed with two different kinds of non-uniform polarization profiles — the first of which enhances the second in-plane extensional mode in the impedance spectrum and the second polarization profile suppresses all the electromechanical resonances pertaining to the in-plane extensional modes and selectively excites only the in-plane bending modes. The proposed approach of using non-uniformly polarized wafers finds application in designing multi-frequency sensors/transducers, frequency-tuned receivers, acoustic beamforming, and other non-traditional applications such as information storage.

36 MATERIALS SCIENCE↗

Thermoosmotic Transport in Nanochannels Grafted with pH-responsive Polyelectrolyte Brushes Modelled Using Augmented Strong Stretching Theory

In this paper, we develop a theory to establish that the thermoosmotic (TOS) effects, induced by the application of an axial temperature gradient, leads to a massive enhancement in liquid transport in nanochannels grafted with charged polyelectrolyte (PE) brushes. Here, we quantify the TOS transport by quantifying the induced electric field and the induced TOS flow field. The different components of the electric field, namely the ionic component, the thermal component, and the osmotic component, as well as the contributions of different ions to these components, are quantified. Furthermore, we express the TOS velocity as a combination of chemiosmotic (COS), thermal, and electroosmotic (EOS) components. The COS and the thermal components augment each other and the overall strength and direction of the TOS flow are dictated by the direction and the relative strength of the EOS component. Most importantly, we compare the cases of brush-grafted nanochannels with those of the brush-free nanochannels of identical surface charge densities: the TOS transport is massively augmented in the brush-grafted nanochannels attributed to the combination of the localization of the EDL (and hence any body force that depends on the EDL charge density) away from the nanochannel wall (i.e., the location of the maximum drag force) and the presence of a possible molecular slip (experienced by the liquid) along the brush surface.

42 ENGINEERING↗

Atomic-Level Structure of Mesoporous Hexagonal Boron Nitride Determined by High-Resolution Solid-State Multinuclear Magnetic Resonance Spectroscopy and Density Functional Theory Calculations

Mesoporous hexagonal boron nitride (p-BN) has received significant attention over the last decade as a promising candidate for water cleaning/pollutant removal and hydrogen storage applications. In this work, high-resolution solid-state NMR spectroscopy and plane-wave density-functional theory (DFT) calculations are used to obtain an atomic-level description of p-BN. 1 H– 15 N or 1 H– 14 N heteronuclear (HETCOR) correlation experiments recorded with either conventional NMR at room temperature or dynamic nuclear polarization surface-enhanced spectroscopy (DNP-SENS) at ca. 100 K reveal NB 2 H, NBH 2 , NBH 3 + species residing on the edges of BN sheets. Ultra-high field 35.2 T 11 B NMR spectroscopy was used to resolve 11 B NMR signals from BN 3 , BN 2 O x (OH) 1–x (x = 0–1), BNO x (OH) 2–x (x = 0–2), BO x (OH) 3–x (x = 0–3), and BO x (OH) 4–x – (x = 0–4). Importantly, 2D 11 B dipolar double-quantum–single-quantum homonuclear correlation spectra reveal that many pore/defect sites are composed of boron oxide/hydroxide clusters connected to the BN framework through BN 2 O units. 1D and 2D 11 B{ 15 N} HETCOR NMR experiments, in addition to plane-wave DFT calculations of nine different structural models, further confirm the assignment of all NMR signals. The detailed structure determination of the pore and edge/defect sites within p-BN should further enable the rational design and development of next-generation p-BN-based materials. In addition, the techniques outlined here should be applicable to determine structure within other porous and/or boron-based materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Rechargeable Alkali-Ion Battery Materials: Theory and Computation

Since its development in the 1970s, the rechargeable alkali-ion battery has proven to be a truly transformative technology, providing portable energy storage for devices ranging from small portable electronics to sizable electric vehicles. Here, we present a review of modern theoretical and computational approaches to the study and design of rechargeable alkali-ion battery materials. Additionally, starting from fundamental thermodynamics and kinetics phenomenological equations, we rigorously derive the theoretical relationships for key battery properties, such as voltage, capacity, alkali diffusivity, and other electrochemically relevant computable quantities. We then present an overview of computational techniques for the study of rechargeable alkali-ion battery materials, followed by a critical review of the literature applying these techniques to yield crucial insights into battery operation and performance. Finally, we provide perspectives on outstanding challenges and opportunities in the theory and computation of rechargeable alkali-ion battery materials.

25 ENERGY STORAGE↗

Forces from Stochastic Density Functional Theory under Nonorthogonal Atom-Centered Basis Sets

We develop a formalism for calculating forces on the nuclei within the linear-scaling stochastic density functional theory (sDFT) in a nonorthogonal atom-centered basis set representation (Fabian et al. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2019, 9, e1412, 10.1002/wcms.1412) and apply it to the Tryptophan Zipper 2 (Trp-zip2) peptide solvated in water. We use an embedded-fragment approach to reduce the statistical errors (fluctuation and systematic bias), where the entire peptide is the main fragment and the remaining 425 water molecules are grouped into small fragments. We analyze the magnitude of the statistical errors in the forces and find that the systematic bias is of the order of 0.065 eV/Å (~1.2 × 10 -3 E h / a 0 ) when 120 stochastic orbitals are used, independently of system size. This magnitude of bias is sufficiently small to ensure that the bond lengths estimated by stochastic DFT (within a Langevin molecular dynamics simulation) will deviate by less than 1% from those predicted by a deterministic calculation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Comprehensive Basis-Set Testing of Extended Symmetry-Adapted Perturbation Theory and Assessment of Mixed-Basis Combinations to Reduce Cost

Hybrid or “extended” symmetry-adapted perturbation theory (XSAPT) replaces traditional SAPT’s treatment of dispersion with better performing alternatives while at the same time extending two-body (dimer) SAPT to a many-body treatment of polarization using a self-consistent charge embedding procedure. The present work presents a systematic study of how XSAPT interaction energies and energy components converge with respect to the choice of Gaussian basis set. Errors can be reduced in a systematic way using correlation-consistent basis sets, with aug-cc-pVTZ results converged within <0.1 kcal/mol. Similar (if slightly less systematic) behavior is obtained using Karlsruhe basis sets at much lower cost, and we introduce new versions with limited augmentation that are even more efficient. Pople-style basis sets, which are more efficient still, often afford good results if a large number of polarization functions are included. The dispersion models used in XSAPT afford much faster basis-set convergence as compared to the perturbative description of dispersion in conventional SAPT, meaning that “compromise” basis sets (such as jun-cc-pVDZ) are no longer required and benchmark-quality results can be obtained using triple-ζ basis sets. The use of diffuse functions proves to be essential, especially for the description of hydrogen bonds. As a result, the “δ(Hartree–Fock)” correction for high-order induction can be performed in double-ζ basis sets without significant loss of accuracy, leading to a mixed-basis approach that offers 4× speedup over the existing (cubic scaling) XSAPT approach.

74 ATOMIC AND MOLECULAR PHYSICS↗

Implementation and Validation of Constrained Density Functional Theory Forces in the CP2K Package

Constrained density functional theory (CDFT) is a powerful tool for the prediction of electron transfer parameters in condensed phase simulations at a reasonable computational cost. In this work we present an extension to CDFT in the popular mixed Gaussian/plane wave electronic structure package CP2K, implementing the additional force terms arising from a constraint based on Hirshfeld charge partitioning. This improves upon the existing Becke partitioning scheme, which is prone to give unphysical atomic charges. We verify this implementation for a variety of systems: electron transfer in (H 2 O) 2 + in vacuum, electron tunnelling between oxygen vacancy centers in solid MgO, and electron self-exchange in aqueous Ru 2+ -Ru 3+ . We find good agreement with previous plane-wave CDFT results for the same systems, but at a significantly lower computational cost, and we discuss the general reliability of condensed phase CDFT calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Computational Protocol to Evaluate Electron–Phonon Interactions Within Density Matrix Perturbation Theory

We present a computational protocol, based on density matrix perturbation theory, to obtain non-adiabatic, frequency-dependent electron–phonon self-energies for molecules and solids. Our approach enables the evaluation of electron–phonon interaction using hybrid functionals, for spin-polarized systems, and the computational overhead to include dynamical and non-adiabatic terms in the evaluation of electron–phonon self-energies is negligible. We discuss results for molecules, as well as pristine and defective solids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An Algorithm for Atom-Centered Lossy Compression of the Atomic Orbital Basis in Density Functional Theory Calculations

Large atomic-orbital (AO) basis sets of at least triple and preferably quadruple-ζ (QZ) size are required to adequately converge Kohn–Sham density functional theory (DFT) calculations toward the complete basis set limit. However, incrementing the cardinal number by one nearly doubles the AO basis dimension, and the computational cost scales as the cube of the AO dimension, so this is very computationally demanding. Here, in this work, we develop and test a threshold-based natural atomic orbital (NAO) scheme in which ϵ-NAOs are obtained as eigenfunctions of atomic blocks of the density matrix in a one-center orthogonalized representation. This enables compression of the AO basis that is optimal for a given threshold, 10 –ϵ , by discarding NAOs with occupation numbers below that threshold. Extensive pilot test calculations using the Hartree–Fock functional and taking the converged density matrix as input suggest that a threshold of 10 –5 can yield a compression factor (ratio of AO to compressed ϵ-NAO dimension) between 2.5 and 4.5 for the QZ pc-3 basis. The errors in relative energies are typically less than 0.1 kcal/mol when the compressed basis is used instead of the uncompressed basis. Between 10 and 100 times smaller errors (i.e., usually less than 0.01 kcal/mol) can be obtained with a threshold 10 –7 , while the compression factor is typically between 2 and 2.5.

basis sets↗

One-Body Properties and Their Perturbative Accuracy with Aufbau Suppressed Coupled Cluster Theory

In this work we derived and implemented the calculation of the one-body reduced density matrix for Aufbau suppressed coupled cluster theory, from which excited state natural orbitals and one-body properties, like atomic populations and dipole moments, are obtained. We utilized the natural orbitals to refine the ASCC solution for simple valence and Rydberg systems, exploring the process of repeatedly solving the ASCC equations in successive natural orbital bases to achieve independence from the starting molecular orbitals. For dipole moments in small molecules where high-level comparison data is available, we find that the accuracy of ASCC essentially matches that of linear response and equation-of-motion coupled cluster as long as care is taken to preserve the response's perturbative completeness.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Photodissociation and Theory to Investigate Uranium Oxide Cluster Cations

Uranium oxide cluster cations of the form U n O m + are produced by laser vaporization of a depleted uranium rod in a pulsed supersonic expansion. Ions are mass-analyzed and mass-selected with a time-of-flight spectrometer and studied with UV laser multiphoton dissociation. Cations of the stoichiometry UO 2 (UO 3 ) n + were observed as photofragments from all photodissociated cluster cations. (UO 3 ) n + clusters were also observed to result from dissociation of larger (UO 3 ) n + clusters, with UO 3 neutral as a common leaving group. Electronic structure calculations were used to investigate the stability of the prominent uranium oxide cluster cations using density functional theory (DFT) with the hybrid B3LYP exchange-correlation functional and at the CCSD(T) level and cc-pVnZ-PP basis sets (n = D,T), including diffuse orbitals as computational expense and availability permitted. Clustering energies, relative energies and dissociation energies of the cations are reported. The lowest energy neutral (UO 3 ) n clusters up to n = 3 are rings, n = 4 and 5 are chains with very low energy rings, and n = 6 is 3D. Finally, the lowest energy structures for UO 2 (UO 3 ) n + are composed of uranyl-like UO 2 + units bound by bridging oxygens to other UO 2 2+ units for n = 2 and 3, and for n = 4 a more complex 3D structure is predicted.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗