Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “derived lattice”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Thermoelectric Performance of the 2D Bi 2 Si 2 Te 6 Semiconductor

Bi 2 Si 2 Te 6 , a 2D compound, is a direct band gap semiconductor with an optical band gap of 0.25 eV, and is a promising thermoelectric material. Single-phase Bi 2 Si 2 Te 6 is prepared by a scalable ball-milling and annealing process and the highly densified polycrystalline samples are prepared by spark plasma sintering. Bi 2 Si 2 Te 6 shows a p-type semiconductor transport behavior and exhibits an intrinsically low lattice thermal conductivity of ~0.48 Wm -1 K -1 (cross-plane) at 573 K. The first-principles density functional theory calculations indicate that such low lattice thermal conductivity is derived from the interactions between acoustic phonons and low-lying optical phonons, local vibrations of Bi, the low Debye temperature and strong anharmonicity result from the unique 2D crystal structure and metavalent bonding of Bi 2 Si 2 Te 6 . The Bi 2 Si 2 Te 6 exhibits an optimal figure of merit ZT of ~0.51 at 623 K, which can be further enhanced by the substitution of Bi with Pb. Pb doping leads to a large increase in power factor S 2 σ, from ~4.0 μWcm -1 K- 2 of Bi 2 Si 2 Te 6 to ~8.0 μWcm -1 K -2 of Bi 1.98 Pb 0.02 Si 2 Te 6 at 775 K, owing to the increase in carrier concentration. Moreover, Pb doping induces a further reduction in the lattice thermal conductivity to ~0.38 Wm -1 K -1 (cross-plane) at 623 K in Bi 1.98 Pb 0.02 Si 2 Te 6 ,. The simultaneous optimization of the power factor and lattice thermal conductivity achieves a peak ZT of ~0.90 at 723 K and a high average ZT of ~0.66 at 400–773 K in Bi 1.98 Pb 0.02 Si 2 Te 6 .

36 MATERIALS SCIENCE↗
DOE OSTI2021

Generation of the TSL for Zirconium Hydrides from Ab Initio Methods

Zirconium hydride (ZrHx) is a moderator material used in TRIGA and other reactors that may exist in multiple phases with varying stoichiometry, which include the δ phase and the ϵ phase. Current ENDF/B-VIII.0 ZrHx thermal scattering law (TSL) evaluations do not distinguish between phases. These sub-libraries were generated with the LEAPR module of NJOY using historic phonon spectra derived from a central force model and assume incoherent elastic scattering for both bound hydrogen and zirconium, which neglects the effects of crystal structures important for scattering from zirconium bound in ZrHx. In this work, the TSLs for hydrogen and zirconium bound in δ-ZrHx and ϵ-ZrH2 were generated from phonon spectra derived from modern ab initio lattice dynamics methods and ab initio molecular dynamics. Subsequently, TSLs for hydrogen and zirconium in ZrHx and ZrH2 were generated using the Full Law Analysis Scattering System Hub (FLASSH) code. The built-in generalized coherent elastic routine was used to generate the previously neglected elastic contribution from zirconium for this material. The present TSLs provide both a re-evaluation of the current ZrH sub-libraries and expansion of the set of TSLs available for the examination of neutrons in systems with zirconium hydride, permitting explicit treatment of δ and ϵ phases.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Many-Exciton Quantum Dynamics in a Ruddlesden–Popper Tin Iodide

We present a study on the many-body exciton interactions in a Ruddlesden–Popper tin halide, namely, (PEA) 2 SnI 4 (PEA = phenylethylammonium), using coherent two-dimensional electronic spectroscopy. The optical dephasing times of the third-order polarization observed in these systems are determined by exciton many-body interactions and lattice fluctuations. We investigate the excitation-induced dephasing (EID) and observe a significant reduction of the dephasing time with increasing excitation density as compared to its lead counterpart (PEA) 2 PbI 4 , which we have previously reported in a separate publication [J. Chem. Phys. 2020, 153, 164706]. Surprisingly, we find that the EID interaction parameter is four orders of magnitude higher in (PEA) 2 SnI 4 than that in (PEA) 2 PbI 4 . This increase in the EID rate may be due to exciton localization arising from a more statically disordered lattice in the tin derivative. This is supported by the observation of multiple closely spaced exciton states and the broadening of the linewidth with increasing population time (spectral diffusion), which suggests a static disordered structure relative to the highly dynamic lead-halide. Additionally, we find that the exciton nonlinear coherent lineshape shows evidence of a biexcitonic state with low binding energy (<10 meV) not observed in the lead system. We model the lineshapes based on a stochastic scattering theory that accounts for the interaction with a nonstationary population of dark background excitations. Our study provides evidence of differences in the exciton quantum dynamics between tin- and lead-based Ruddlesden–Popper metal halides (RPMHs) and links them to the exciton–exciton interaction strength and the static disorder aspect of the crystalline structure.

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

Heterostructural Alloy Phase Diagram for (Cd 1-x Zn x ) 3 As 2

Alloying the topological semimetal Cd 3 As 2 with Zn 3 As 2 provides a potential route for controlling the electronic properties. We predict the alloy phase diagram from first-principles calculations, considering that both end members have a crystal structure derived from the antifluorite lattice, but with different arrangements of the unoccupied cation sites. To overcome the limitations of the regular solution approximation and to include short-range order effects, we perform Monte Carlo simulations, parameterize the temperature dependence of the mixing enthalpy ΔH m , and perform thermodynamic integration of the free energy. The resulting phase diagram exhibits features that are unique to heterostructural alloy systems and provides computational predictions of solubility limits and composition ranges that are stable against spinodal decomposition.

36 MATERIALS SCIENCE↗
NASA NTRS1979

A Digital Program for Calculating the Interaction Between Flexible Structures, Unsteady Aerodynamics and Active Controls

A computer program, ISAC, is described which calculates the stability and response of a flexible airplane equipped with active controls. The equations of motion relative to a fixed inertial coordinate system are formulated in terms of the airplane's rigid body motion and its unrestrained normal vibration modes. Unsteady aerodynamic forces are derived from a doublet lattice lifting surface theory. The theoretical basis for the program is briefly explained together with a description of input data and output results.

Peele, E. L.↗
NASA NTRS1996

Subsonic Wing Optimization for Handling Qualities Using ACSYNT

The capability to accurately and rapidly predict aircraft stability derivatives using one comprehensive analysis tool has been created. The PREDAVOR tool has the following capabilities: rapid estimation of stability derivatives using a vortex lattice method, calculation of a longitudinal handling qualities metric, and inherent methodology to optimize a given aircraft configuration for longitudinal handling qualities, including an intuitive graphical interface. The PREDAVOR tool may be applied to both subsonic and supersonic designs, as well as conventional and unconventional, symmetric and asymmetric configurations. The workstation-based tool uses as its model a three-dimensional model of the configuration generated using a computer aided design (CAD) package. The PREDAVOR tool was applied to a Lear Jet Model 23 and the North American XB-70 Valkyrie.

Soban, Danielle Suzanne↗
DOE OSTI2022

Factorization connecting continuum & lattice TMDs

Transverse-momentum-dependent parton distribution functions (TMDs) can be studied from first principles by a perturbative matching onto lattice-calculable quantities: so-called lattice TMDs, which are a class of equal-time correlators that includes quasi-TMDs and TMDs in the Lorentz-invariant approach. We introduce a general correlator that includes as special cases these two Lattice TMDs and continuum TMDs, like the Collins scheme. Then, to facilitate the derivation of a factorization relation between lattice and continuum TMDs, we construct a new scheme, the Large Rapidity (LR) scheme, intermediate between the Collins and quasi-TMDs. The LR and Collins schemes differ only by an order of limits, and can be matched onto one another by a multiplicative kernel. We show that this same matching also holds between quasi and Collins TMDs, which enables us to prove a factorization relation between these quantities to all orders in α s . Our results imply that there is no mixing between various quark flavors or gluons when matching Collins and quasi TMDs, making the lattice calculation of individual flavors and gluon TMDs easier than anticipated. We cross-check these results explicitly at one loop and discuss implications for other physical-to-lattice scheme factorizations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2003

Deriving Laws from Ordering Relations

It took much effort in the early days of non-Euclidean geometry to break away from the mindset that all spaces are flat and that two distinct parallel lines do not cross. Up to that point, all that was known was Euclidean geometry, and it was difficult to imagine anything else. We have suffered a similar handicap brought on by the enormous relevance of Boolean algebra to the problems of our age-logic and set theory. Previously, I demonstrated that the algebra of questions is not Boolean, but rather is described by the free distributive algebra. To get to this stage took much effort, as many obstacles-most self-placed-had to be overcome. As Boolean algebras were all I had ever known, it was almost impossible for me to imagine working with an algebra where elements do not have complements. With this realization, it became very clear that the sum and product rules of probability theory at the most basic level had absolutely nothing to do with the Boolean algebra of logical statements. Instead, a measure of degree of inclusion can be invented for many different partially ordered sets, and the sum and product rules fall out of the associativity and distributivity of the algebra. To reinforce this very important idea, this paper will go over how these constructions are made, while focusing on the underlying assumptions. I will derive the sum and product rules for a distributive lattice in general and demonstrate how this leads to probability theory on the Boolean lattice and is related to the calculus of quantum mechanical amplitudes on the partially ordered set of experimental setups. I will also discuss the rules that can be derived from modular lattices and their relevance to the cross-ratio of projective geometry.

Knuth, Kevin H.↗
NASA NTRS1991

Constraints on the Parent Bodies of Collected Interplanetary Dust Particles

Samples of Interplanetary Dust Particles (IDPs) have now been collected from the stratosphere, from the Earth's ocean beds, and from the ice caps of Greenland and Antarctica. The most likely candidates for the sources of these particles are comets and asteroids. Comparison of the infrared spectra, elemental compositions, and mineralogy of the collected dust with atmospheric entry models and data obtained from cometary probes and telescopic observations has provided important constraints on the possible sources of the various types of collected dust. These constraints lead to the following conclusions. First, most of the deep sea, Greenland, and Antarctic spherules larger than 100 microns are derived from asteroids. Second, the stratospheric IDPs dominated by hydrated layer-lattice silicate minerals are also most likely derived from asteroids. Finally, the stratospheric IDPs dominated by the anhydrous minerals olivine and pyroxene are most likely from comets. The consequences of these parent body assignments are discussed.

Sandford, S. A.↗
DOE OSTI2020

Electronic structures of a diagonally striped lattice: Multiple ( N - 1 ) -fold degenerate flat bands

We investigate the electronic structure of an interesting two-dimensional (2D) diagonally striped lattice (DSL). It consists of arrays of N different types of “atoms” or “molecules” with a fixed sequence in both horizontal and vertical directions. Using a tight-binding model, we show that the DSL (N>2) has the symmetry group Amm2 and Pmm2, and exhibits an oscillatory metallic and insulating phases for the odd- and even- N number, respectively. Some conventional 2D lattices, such as the Lieb lattice, can be related to the derivative DSL via vacancy formation by removing one type of atom in every other (N-1) row of the original DSL. Interestingly, there are multiple (N-1)-fold degenerate flat bands in one group of the derivative DSLs, which affords a unique platform for studying many-body physics. In addition to atomic and molecular lattices, we suggest other artificial DSLs, such as photonic and phononic, to be also constructed with possibly interesting properties.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Lieb-Schultz-Mattis, Luttinger, and 't Hooft - anomaly matching in lattice systems

We analyze lattice Hamiltonian systems whose global symmetries have ’t Hooft anomalies. As is common in the study of anomalies, they are probed by coupling the system to classical background gauge fields. For flat fields (vanishing field strength), the nonzero spatial components of the gauge fields can be thought of as twisted boundary conditions, or equivalently, as topological defects. The symmetries of the twisted Hilbert space and their representations capture the anomalies. We demonstrate this approach with a number of examples. In some of them, the anomalous symmetries are internal symmetries of the lattice system, but they do not act on-site. (We clarify the notion of “on-site action.”) In other cases, the anomalous symmetries involve lattice translations. Using this approach we frame many known and new results in a unified fashion. In this work, we limit ourselves to 1+1d systems with a spatial lattice. In particular, we present a lattice system that flows to the c=1 c = 1 compact boson system with any radius (no BKT transition) with the full internal symmetry of the continuum theory, with its anomalies and its T-duality. As another application, we analyze various spin chain models and phrase their Lieb-Shultz-Mattis theorem as an ’t Hooft anomaly matching condition. We also show in what sense filling constraints like Luttinger theorem can and cannot be viewed as reflecting an anomaly. As a by-product, our understanding allows us to use information from the continuum theory to derive some exact results in lattice model of interest, such as the lattice momenta of the low-energy states.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2000

Interference Lattice-based Loop Nest Tilings for Stencil Computations

A common method for improving performance of stencil operations on structured multi-dimensional discretization grids is loop tiling. Tile shapes and sizes are usually determined heuristically, based on the size of the primary data cache. We provide a lower bound on the numbers of cache misses that must be incurred by any tiling, and a close achievable bound using a particular tiling based on the grid interference lattice. The latter tiling is used to derive highly efficient loop orderings. The total number of cache misses of a code is the sum of (necessary) cold misses and misses caused by elements being dropped from the cache between successive loads (replacement misses). Maximizing temporal locality is equivalent to minimizing replacement misses. Temporal locality of loop nests implementing stencil operations is optimized by tilings that avoid data conflicts. We divide the loop nest iteration space into conflict-free tiles, derived from the cache miss equation. The tiling involves the definition of the grid interference lattice an equivalence class of grid points whose images in main memory map to the same location in the cache-and the construction of a special basis for the lattice. Conflicts only occur on the boundaries of the tiles, unless the tiles are too thin. We show that the surface area of the tiles is bounded for grids of any dimensionality, and for caches of any associativity, provided the eccentricity of the fundamental parallelepiped (the tile spanned by the basis) of the lattice is bounded. Eccentricity is determined by two factors, aspect ratio and skewness. The aspect ratio of the parallelepiped can be bounded by appropriate array padding. The skewness can be bounded by the choice of a proper basis. Combining these two strategies ensures that pathologically thin tiles are avoided. They do not, however, minimize replacement misses per se. The reason is that tile visitation order influences the number of data conflicts on the tile boundaries. If two adjacent tiles are visited successively, there will be no replacement misses on the shared boundary. The iteration space may be covered with pencils larger than the size of the cache while avoiding data conflicts if the pencils are traversed by a scanning-face method. Replacement misses are incurred only on the boundaries of the pencils, and the number of misses is minimized by maximizing the volume of the scanning face, not the volume of the tile. We present an algorithm for constructing the most efficient scanning face for a given grid and stencil operator. In two dimensions it is based on a continued fraction algorithm. In three dimensions it follows Voronoi's successive minima algorithm. We show experimental results of using the scanning face, and compare with canonical loop orderings.

VanderWijngaart, Rob F.↗
DOE OSTI2020

On the pressure and temperature dependence of adsorption densities and other thermodynamic properties in gas shales

Black shale source rocks have become a key natural gas resource in the US and China. Unlike in conventional hydrocarbon reservoirs, a significant fraction of the gas in shales is adsorbed onto the inordinate specific surface areas (SSA) of micropores inside organic matter with densities exceeding those of bulk gas by an order of magnitude. Yet, the distribution of adsorbed versus bulk gas remains poorly understood. Experiments measure excess adsorption isotherms, which are corrected to absolute adsorption amounts by a factor that involves the density, &#x3C1; abs , of the adsorption layer(s). Constant values of &#x3C1; abs are widely used in the literature but, we argue, are inconsistent with the pressure- and temperature-dependence of adsorption behavior. In this work, we propose a more reliable approach that assumes a constant volume of each adsorption layer that can be determined from the SSA of the substrate. We measure nine high-pressure methane excess adsorption isotherms up to 15 MPa for three shale samples from China and Germany at temperatures of 65, 75, and 95 °C, and low-pressure carbon dioxide and nitrogen isotherms to estimate the SSA. We also investigate another 24 isotherms at 35, 50, and 60° C or 8 Chinese shale samples and 3 isotherms at 60, 90, and 120° C for an Argentinian shale sample from the literature. A new variation of the multi-layer Ono-Kondo lattice model is introduced to derive the temperature and pressure dependent adsorption layer densities, which are subsequently used to obtain the absolute adsorption isotherms of interest. Finally, this practical methodology has the potential to significantly improve our estimates of shale gas-in-place resources, and adsorption modeling more broadly.

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

Finite-volume quantization condition from the N/D representation

We propose a new model-independent method for determining hadronic resonances from lattice QCD. The formalism is derived from the general principles of unitarity and analyticity, as encoded in the N/D representation of a partial-wave two-body amplitude. The associated quantization condition relates the finite-volume spectrum to the infinite-volume numerator $\mathscr{N}$, used to reconstruct the scattering amplitude from dispersive relations. Unlike the original Lüscher condition, this new formalism is valid for energies coinciding with the left-hand cuts from arbitrary one- and multiparticle exchanges.

Dawid, Sebastian M. [Univ. of Washington, Seattle,↗
DOE OSTI2020

Tunable Magnetic Exchange between Rare-Earth Metal 5d and Iron 3d States: A Case Study of the Multiple Magnetic Transitions in Gd 6 FeBi 2 and the Solid Solutions Dy 6– x Gd x FeBi 2 (1 ≤ x ≤ 5) with Curie Temperatures in the Range 120–350 K

In this paper, the room-temperature magnet Gd 6 FeBi 2 was comprehensively characterized by means of temperature-dependent single-crystal X-ray diffraction, magnetization measurements, and experimental electron density and electronic structure computations. This work explores the electron-spin effects of Fe on this structure and suggests that Gd6FeBi2 shows structural features and exchange interactions, which are clearly distinguishable from other analogues, including the solid solutions Dy 6–x Gd x FeBi 2 (0 ≤ x ≤ 5). The unique traits of Gd 6 FeBi 2 and its derivatives encompass abnormal variations of lattice parameters with the temperature, as well as the presence of multiple magnetic transitions of complex nature. Based on the comprehensive analyses, it can be suggested that the unique magnetic response in this material is the result from tunable Fe spin states, which are coupled with strong Gd 5d–Fe 3d interactions, which differ from other members of this large, structural family. In the ground state, according to density functional theory calculations, the Fe atom carries two net electrons, whose spins are coupled in an antiparallel fashion to the spins of the Gd electrons. Near the Curie temperature, the Fe net moment is reduced intermittently, coupled with a multiple of magnetic transitions. The magnetic correlations are manifested in unexpected variations of the lattice parameters as a function of temperature, suggesting spin–lattice interactions. Furthermore, this study emphasizes the important role of Gd 5d electrons in tuning Fe 3d-based magnetic contribution, which enables better understanding on the spin coupling in other related compounds and sheds light on the development of new magnetic and spintronic compounds.

36 MATERIALS SCIENCE↗
DOE OSTI2025

MS25: Materials Science-Focused Benchmark Data Set for Machine Learning Interatomic Potentials

Here, we present MS25, a benchmark data set for evaluating machine learning interatomic potentials (MLIPs) across diverse materials-relevant systems including MgO surfaces, liquid water, zeolites, a catalytic Pt surface reaction, high-entropy alloys (HEAs), and disordered Zr-oxides. Five MLIP architectures (MACE, NequIP, Allegro, MTP, and Torch-ANI) are trained and tested, focusing not only on traditional metrics (energies, forces, and stresses) but also explicitly validating derived physical observables such as lattice constants, volumes, and reaction barriers. We find that most models reach comparable accuracy on standard error metrics across the simple systems, although equivariant MLIPs offer 1.5–2× improvements over nonequivariant MLIPs in energy and force error for structurally complex or compositionally disordered environments such as HEAs and Zr–O systems. Our analysis highlights that low errors in energy and force predictions do not guarantee reliable observables, emphasizing the necessity of explicit validation. We demonstrate limitations in cross-framework transferability, as models trained on one zeolite framework (CHA) fail to reliably generalize to predictions of structurally distinct frameworks (e.g., MFI). Size-extensive tests show some dependence on system size for MgO, resulting from forced periodicity. The HEA and Zr–O data sets are identified as challenging tests for future benchmarks and MLIP model architecture developments as they show significant differentiation in error between MLIP architectures and are still relatively difficult at 1000 training images. Moving forward, we recommend that benchmarking efforts shift their focus from marginal accuracy improvements in energy and force errors toward identifying and understanding model failure modes, rigorously assessing transferability, and evaluating how their errors affect observable predictions. For researchers looking to choose an MLIP architecture, we suggest selecting equivariant MLIP architectures if the complexity of the system is a challenge. For simple materials problems, auxiliary features such as integration with molecular dynamics engines, trade-offs between computational data set generation cost vs MLIP inference speed, and framework integration may play a more important decision factor than small differences in error metrics that are unlikely to matter for production-level research.

chemical structure↗
DOE OSTI2023

A j eff = 1/2 Kitaev material on the triangular lattice: the case of NaRuO 2

Motivated by recent reports of a quantum-disordered ground state in the triangular lattice compound NaRuO 2 , we derive a j eff = 1/2 magnetic model for this system by means of first-principles calculations. The pseudospin Hamiltonian is dominated by bond-dependent off-diagonal Γ interactions, complemented by a ferromagnetic Heisenberg exchange and a notably antiferromagnetic Kitaev term. In addition to bilinear interactions, we find a sizable four-spin ring exchange contribution with a strongly anisotropic character, which has been so far overlooked when modeling Kitaev materials. The analysis of the magnetic model, based on the minimization of the classical energy and exact diagonalization of the quantum Hamiltonian, points toward the existence of a rather robust easy-plane ferromagnetic order, which cannot be easily destabilized by physically relevant perturbations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗