Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “permutation”

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

A survey on application of representation theory to molecular vibration

Representations Theory is used extensively in many of the physical sciences as every physical system has a symmetry group G. Various differential equations determine the vibration of a molecule, and the symmetry group of the molecule acts on the space of solutions of these equations. In this paper we use CH{sub 4} (methane) molecule, which has four hydrogen atoms at the corners of a regular tetrahedron, and a carbon atom at the center of the tetrahedron. The four hydrogen atoms in CH{sub 4} are permuted by the action of the symmetry group and this action fixes the carbon atom. At each of the 5 vertices, we assign three unit vectors, called the standard basis vectors in directions of the three edges which are joined to the vertex. The symmetry group G of the molecules permutes the 15 standard basis vectors, so we may regard Q{sup 15} as a GG By expressing Q{sup 15} as a direct sum of irreducible GG-modules, the problem of finding the normal modes of vibration is reduced to that of computing the eigenvectors of some small matrices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Generalized constructive tree weights

The Loop Vertex Expansion (LVE) is a quantum field theory (QFT) method which explicitly computes the Borel sum of Feynman perturbation series. This LVE relies in a crucial way on symmetric tree weights which define a measure on the set of spanning trees of any connected graph. In this paper we generalize this method by defining new tree weights. They depend on the choice of a partition of a set of vertices of the graph, and when the partition is non-trivial, they are no longer symmetric under permutation of vertices. Nevertheless we prove they have the required positivity property to lead to a convergent LVE; in fact we formulate this positivity property precisely for the first time. Our generalized tree weights are inspired by the Brydges-Battle-Federbush work on cluster expansions and could be particularly suited to the computation of connected functions in QFT. Several concrete examples are explicitly given.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2014

Communication: A benchmark-quality, full-dimensional ab initio potential energy surface for Ar-HOCO

A full-dimensional, global ab initio potential energy surface (PES) for the Ar-HOCO system is presented. The PES consists of a previous intramolecular ab initio PES for HOCO [J. Li, C. Xie, J. Ma, Y. Wang, R. Dawes, D. Xie, J. M. Bowman, and H. Guo, J. Phys. Chem. A 116, 5057 (2012)], plus a new permutationally invariant interaction potential based on fitting 12 432 UCCSD(T)-F12a/aVDZ counterpoise-corrected energies. The latter has a total rms fitting error of about 25 cm{sup −1} for fitted interaction energies up to roughly 12 000 cm{sup −1}. Two additional fits are presented. One is a novel very compact permutational invariant representation, which contains terms only involving the Ar-atom distances. The rms fitting error for this fit is 193 cm{sup −1}. The other fit is the widely used pairwise one. The pairwise fit to the entire data set has an rms fitting error of 427 cm{sup −1}. All of these potentials are used in preliminary classical trajectory calculations of energy transfer with a focus on comparisons with the results using the benchmark potential.

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

Complete spectrum of the infinite-U Hubbard ring using group theory

We present a full analytical solution of the multiconfigurational strongly correlated mixed-valence problem corresponding to the N-Hubbard ring filled with N−1 electrons, and infinite on-site repulsion. While the eigenvalues and the eigenstates of the model are known already, analytical determination of their degeneracy is presented here for the first time. The full solution, including degeneracy count, is achieved for each spin configuration by mapping the Hubbard model into a set of Hückel-annulene problems for rings of variable size. The number and size of these effective Hückel annulenes, both crucial to obtain Hubbard states and their degeneracy, are determined by solving a well-known combinatorial enumeration problem, the necklace problem for N−1 beads and two colors, within each subgroup of the C{sub N−1} permutation group. Symmetry-adapted solution of the necklace enumeration problem is finally achieved by means of the subduction of coset representation technique [S. Fujita, Theor. Chim. Acta 76, 247 (1989)], which provides a general and elegant strategy to solve the one-hole infinite-U Hubbard problem, including degeneracy count, for any ring size. The proposed group theoretical strategy to solve the infinite-U Hubbard problem for N−1 electrons is easily generalized to the case of arbitrary electron count L, by analyzing the permutation group C{sub L} and all its subgroups.

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

A nine-dimensional ab initio global potential energy surface for the H{sub 2}O{sup +} + H{sub 2} → H{sub 3}O{sup +} + H reaction

An accurate full-dimensional global potential energy surface (PES) is developed for the title reaction. While the long-range interactions in the reactant asymptote are represented by an analytical expression, the interaction region of the PES is fit to more than 81 000 of ab initio points at the UCCSD(T)-F12b/AVTZ level using the permutation invariant polynomial neural network approach. Fully symmetric with respect to permutation of all four hydrogen atoms, the PES provides a faithful representation of the ab initio points, with a root mean square error of 1.8 meV or 15 cm{sup −1}. The reaction path for this exoergic reaction features an attractive and barrierless entrance channel, a submerged saddle point, a shallow H{sub 4}O{sup +} well, and a barrierless exit channel. The rate coefficients for the title reaction and kinetic isotope effect have been determined on this PES using quasi-classical trajectories, and they are in good agreement with available experimental data. It is further shown that the H{sub 2}O{sup +} rotational enhancement of reactivity observed experimentally can be traced to the submerged saddle point. Using our recently proposed Sudden Vector Projection model, we demonstrate that a rotational degree of freedom of the H{sub 2}O{sup +} reactant is strongly coupled with the reaction coordinate at this saddle point, thus unraveling the origin of the pronounced mode specificity in this reaction.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Probing the folded state and mechanical unfolding pathways of T4 lysozyme using all-atom and coarse-grained molecular simulation

The Bacteriophage T4 Lysozyme (T4L) is a prototype modular protein comprised of an N-terminal and a C-domain domain, which was extensively studied to understand the folding/unfolding mechanism of modular proteins. To offer detailed structural and dynamic insights to the folded-state stability and the mechanical unfolding behaviors of T4L, we have performed extensive equilibrium and steered molecular dynamics simulations of both the wild-type (WT) and a circular permutation (CP) variant of T4L using all-atom and coarse-grained force fields. Our all-atom and coarse-grained simulations of the folded state have consistently found greater stability of the C-domain than the N-domain in isolation, which is in agreement with past thermostatic studies of T4L. While the all-atom simulation cannot fully explain the mechanical unfolding behaviors of the WT and the CP variant observed in an optical tweezers study, the coarse-grained simulations based on the Go model or a modified elastic network model (mENM) are in qualitative agreement with the experimental finding of greater unfolding cooperativity in the WT than the CP variant. Interestingly, the two coarse-grained models predict different structural mechanisms for the observed change in cooperativity between the WT and the CP variant—while the Go model predicts minor modification of the unfolding pathways by circular permutation (i.e., preserving the general order that the N-domain unfolds before the C-domain), the mENM predicts a dramatic change in unfolding pathways (e.g., different order of N/C-domain unfolding in the WT and the CP variant). Based on our simulations, we have analyzed the limitations of and the key differences between these models and offered testable predictions for future experiments to resolve the structural mechanism for cooperative folding/unfolding of T4L.

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

Magnetohydrodynamic turbulence: Observation and experiment

We provide a tutorial on the paradigms and tools of magnetohydrodynamic (MHD) turbulence. The principal paradigm is that of a turbulent cascade from large scales to small, resulting in power law behavior for the frequency power spectrum for magnetic fluctuations E{sub B}(f). We will describe five useful statistical tools for MHD turbulence in the time domain: the temporal autocorrelation function, the frequency power spectrum, the probability distribution function of temporal increments, the temporal structure function, and the permutation entropy. Each of these tools will be illustrated with an example taken from MHD fluctuations in the solar wind. A single dataset from the Wind satellite will be used to illustrate all five temporal statistical tools.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Asymmetric band offsets in silicon heterojunction solar cells: Impact on device performance

Amorphous/crystalline silicon interfaces feature considerably larger valence than conduction band offsets. In this article, we analyze the impact of such band offset asymmetry on the performance of silicon heterojunction solar cells. To this end, we use silicon suboxides as passivation layers—inserted between substrate and (front or rear) contacts—since such layers enable intentionally exacerbated band-offset asymmetry. Investigating all topologically possible passivation layer permutations and focussing on light and dark current-voltage characteristics, we confirm that to avoid fill factor losses, wider-bandgap silicon oxide films (of at least several nanometer thin) should be avoided in hole-collecting contacts. As a consequence, device implementation of such films as window layers—without degraded carrier collection—demands electron collection at the front and hole collection at the rear. Furthermore, at elevated operating temperatures, once possible carrier transport barriers are overcome by thermionic (field) emission, the device performance is mainly dictated by the passivation of its surfaces. In this context, compared to the standard amorphous silicon layers, the wide-bandgap oxide layers applied here passivate remarkably better at these temperatures, which may represent an additional benefit under practical operation conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Communication: Fitting potential energy surfaces with fundamental invariant neural network

A more flexible neural network (NN) method using the fundamental invariants (FIs) as the input vector is proposed in the construction of potential energy surfaces for molecular systems involving identical atoms. Mathematically, FIs finitely generate the permutation invariant polynomial (PIP) ring. In combination with NN, fundamental invariant neural network (FI-NN) can approximate any function to arbitrary accuracy. Because FI-NN minimizes the size of input permutation invariant polynomials, it can efficiently reduce the evaluation time of potential energy, in particular for polyatomic systems. In this work, we provide the FIs for all possible molecular systems up to five atoms. Potential energy surfaces for OH{sub 3} and CH{sub 4} were constructed with FI-NN, with the accuracy confirmed by full-dimensional quantum dynamic scattering and bound state calculations.

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

Learning molecular energies using localized graph kernels

We report that recent machine learning methods make it possible to model potential energy of atomic configurations with chemical-level accuracy (as calculated from ab initio calculations) and at speeds suitable for molecular dynamics simulation. Best performance is achieved when the known physical constraints are encoded in the machine learning models. For example, the atomic energy is invariant under global translations and rotations; it is also invariant to permutations of same-species atoms. Although simple to state, these symmetries are complicated to encode into machine learning algorithms. In this paper, we present a machine learning approach based on graph theory that naturally incorporates translation, rotation, and permutation symmetries. Specifically, we use a random walk graph kernel to measure the similarity of two adjacency matrices, each of which represents a local atomic environment. This Graph Approximated Energy (GRAPE) approach is flexible and admits many possible extensions. Finally, we benchmark a simple version of GRAPE by predicting atomization energies on a standard dataset of organic molecules.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Neural network based quasi-diabatic Hamiltonians with symmetry adaptation and a correct description of conical intersections

Previously, we demonstrated in a report that artificial neural networks (NNs) can be used to generate quasidiabatic Hamiltonians (H d ) that are capable of representing adiabatic energies, energy gradients, and derivative couplings. In this work, two additional issues are addressed. First, symmetry-adapted functions such as permutation invariant polynomials are introduced to account for complete nuclear permutation inversion symmetry. Second, a partially diagonalized representation is introduced to facilitate a better description of near degeneracy points. The diabatization of 1, 2 1 A states of NH 3 is used as an example. The NN fitting findings are compared to that of a previous fitting with symmetry adapted polynomials.

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

Masked Particle Modeling on Sets: Towards Self-Supervised High Energy Physics Foundation Models

Abstract We propose masked particle modeling (MPM) as a self-supervised method for learning generic, transferable, and reusable representations on unordered sets of inputs for use in high energy physics (HEP) scientific data. This work provides a novel scheme to perform masked modeling based pre-training to learn permutation invariant functions on sets. More generally, this work provides a step towards building large foundation models for HEP that can be generically pre-trained with self-supervised learning and later fine-tuned for a variety of down-stream tasks. In MPM, particles in a set are masked and the training objective is to recover their identity, as defined by a discretized token representation of a pre-trained vector quantized variational autoencoder. We study the efficacy of the method in samples of high energy jets at collider physics experiments, including studies on the impact of discretization, permutation invariance, and ordering. We also study the fine-tuning capability of the model, showing that it can be adapted to tasks such as supervised and weakly supervised jet classification, and that the model can transfer efficiently with small fine-tuning data sets to new classes and new data domains.

Heinrich, Lukas (ORCID:0000000240487584)↗
DOE OSTI2024

Topological Quantum Synchronization of Fractionalized Spins

The gapped symmetric phase of the Affleck-Kennedy-Lieb-Tasaki model exhibits fractionalized spins at the ends of an open chain. Here, we show that breaking SU(2) symmetry and applying a global spin-lowering dissipator achieves synchronization of these fractionalized spins. Additional local dissipators ensure convergence to the ground state manifold. In order to understand which aspects of this synchronization are robust within the entire Haldane-gap phase, we reduce the biquadratic term, which eliminates the need for an external field but destabilizes synchronization. Within the ground state subspace, stability is regained using only the global lowering dissipator. These results demonstrate that fractionalized degrees of freedom can be synchronized in extended systems with a significant degree of robustness arising from topological protection. A direct consequence is that permutation symmetries are not required for the dynamics to be synchronized, representing a clear advantage of topological synchronization compared to synchronization induced by permutation symmetries.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Statistical perspective on embrittling potency for intergranular fracture

Embrittling potency is a thermodynamic metric that assesses the influence of solute segregation to a grain boundary (GB) on intergranular fracture. Historically, authors of studies have reported embrittling potency as a single scalar value, assuming a single segregation site of importance at a GB and a particular cleavage plane. However, the topography of intergranular fracture surfaces is not generally known a priori. Accordingly, we, in this study, present a statistical ensemble approach to compute embrittling potency, where many free surface (FS) permutations are systematically considered to model fracture of a GB. The result is a statistical description of the thermodynamics of GB embrittlement. As a specific example, embrittling potency distributions are presented for Cr segregation to sites at two Ni $\langle 111 \rangle$ symmetric tilt GBs using atomistic simulations. We show that the average embrittling potency for a particular GB site, considering an ensemble of FS permutations, is not equal to the embrittling potency computed using the lowest energy pair of FSs. A mean GB embrittlement is proposed, considering both the likelihood of formation of a particular FS and the probability of solute occupancy at each GB site, to compare the relative embrittling behavior of two distinct GBs.

36 MATERIALS SCIENCE↗
DOE OSTI2023

An over-approximation of entropy for elemental multiconfigurational ground state electronic structures

The ability to predict the chemical and physical properties of a material is directly related to the structure and interactions of its electrons. For materials comprised of f-block elements (the lanthanides and actinides found in the last two rows of the periodic table), the complexity of electronic structure has presented great difficulty in understanding, modeling, and predicting material properties. The complexity of multiconfigurational ground state electronic structures is illustrated herein by the combinatorics of electron permutations within individual and cumulative occupancy configurations. A non-integer orbital occupancy representation of multiconfigurational ground states is described for superposition mixing between multiple near-energy degenerate occupancy configurations and generalized in such a way that established ground states are returned by approximation for elements with less-complex electronic structures. By considering the occupancy configurations as statistical mechanics macrostates, and the permutations of electrons as statistical mechanics microstates within those macrostates, an over-approximation of entropy for multiconfigurational elemental ground state electronic structures has been calculated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

An under-approximation of entropy for elemental multiconfigurational ground state electronic structures

A combinatorial approach has been applied to the allowable permutations of quantum electronic configurations under the constraints of Hund's rule for established ground state configurations toward an under-approximation of electronic structure entropy. Combined with a previously reported over-approximation, the approximations are used in conjunction in an attempt to bracket the upper and lower entropy limits for multiconfigurational ground state electronic structure entropy and compared to known standard molar entropies for the elements. This formality has been used for the application of a classical statistical mechanics methodology to be applied to the discrete sets of quantum mechanical states of Pu in order to calculate orbital occupancies in Pu's multiconfigurational ground state. Without consideration of the relative energies of various possible electronic configurations contributing to the multiconfigurational ground state, the calculations are performed under a general energy degeneracy assumption weighted to the number of permutations for specific configurations. The number of configurations assumed to significantly contribute is gradually constrained in order to approach a low-order approximation of orbital occupancies in Pu that are then compared to experimental and other calculated results from the literature.

Beaux, II, Miles F. (ORCID:000000032192626X)↗
DOE OSTI2015

Estimated dose rates to members of the public from external exposure to patients with {sup 131}I thyroid treatment

Purpose: Estimated dose rates that may result from exposure to patients who had been administered iodine-131 ({sup 131}I) as part of medical therapy were calculated. These effective dose rate estimates were compared with simplified assumptions under United States Nuclear Regulatory Commission Regulatory Guide 8.39, which does not consider body tissue attenuation nor time-dependent redistribution and excretion of the administered {sup 131}I. Methods: Dose rates were estimated for members of the public potentially exposed to external irradiation from patients recently treated with {sup 131}I. Tissue attenuation and iodine biokinetics were considered in the patient in a larger comprehensive effort to improve external dose rate estimates. The external dose rate estimates are based on Monte Carlo simulations using the Phantom with Movable Arms and Legs (PIMAL), previously developed by Oak Ridge National Laboratory and the United States Nuclear Regulatory Commission. PIMAL was employed to model the relative positions of the {sup 131}I patient and members of the public in three exposure scenarios: (1) traveling on a bus in a total of six seated or standing permutations, (2) two nursing home cases where a caregiver is seated at 30 cm from the patient’s bedside and a nursing home resident seated 250 cm away from the patient in an adjacent bed, and (3) two hotel cases where the patient and a guest are in adjacent rooms with beds on opposite sides of the common wall, with the patient and guest both in bed and either seated back-to-back or lying head to head. The biokinetic model predictions of the retention and distribution of {sup 131}I in the patient assumed a single voiding of urinary bladder contents that occurred during the trip at 2, 4, or 8 h after {sup 131}I administration for the public transportation cases, continuous first-order voiding for the nursing home cases, and regular periodic voiding at 4, 8, or 12 h after administration for the hotel room cases. Organ specific activities of {sup 131}I in the thyroid, bladder, and combined remaining tissues were calculated as a function of time after administration. Exposures to members of the public were considered for {sup 131}I patients with normal thyroid uptake (peak thyroid uptake of ∼27% of administered {sup 131}I), differentiated thyroid cancer (DTC, 5% uptake), and hyperthyroidism (80% uptake). Results: The scenario with the patient seated behind the member of the public yielded the highest dose rate estimate of seated public transportation exposure cases. The dose rate to the adjacent room guest was highest for the exposure scenario in which the hotel guest and patient are seated by a factor of ∼4 for the normal and differentiated thyroid cancer uptake cases and by a factor of ∼3 for the hyperthyroid case. Conclusions: It was determined that for all modeled cases, the DTC case yielded the lowest external dose rates, whereas the hyperthyroid case yielded the highest dose rates. In estimating external dose to members of the public from patients with {sup 131}I therapy, consideration must be given to (patient- and case-specific) administered {sup 131}I activities and duration of exposure for a more complete estimate. The method implemented here included a detailed calculation model, which provides a means to determine dose rate estimates for a range of scenarios. The method was demonstrated for variations of three scenarios, showing how dose rates are expected to vary with uptake, voiding pattern, and patient location.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Invariant exchange perturbation theory for multicenter systems: Time-dependent perturbations

A formalism of exchange perturbation theory (EPT) is developed for the case of interactions that explicitly depend on time. Corrections to the wave function obtained in any order of perturbation theory and represented in an invariant form include exchange contributions due to intercenter electron permutations in complex multicenter systems. For collisions of atomic systems with an arbitrary type of interaction, general expressions are obtained for the transfer (T) and scattering (S) matrices in which intercenter electron permutations between overlapping nonorthogonal states belonging to different centers (atoms) are consistently taken into account. The problem of collision of alpha particles with lithium atoms accompanied by the redistribution of electrons between centers is considered. The differential and total charge-exchange cross sections of lithium are calculated.

74 ATOMIC AND MOLECULAR PHYSICS↗