Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mathematical physics”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A caloritronics-based Mott neuristor

Abstract Machine learning imitates the basic features of biological neural networks at a software level. A strong effort is currently being made to mimic neurons and synapses with hardware components, an approach known as neuromorphic computing. While recent advances in resistive switching have provided a path to emulate synapses at the 10 nm scale, a scalable neuron analogue is yet to be found. Here, we show how heat transfer can be utilized to mimic neuron functionalities in Mott nanodevices. We use the Joule heating created by current spikes to trigger the insulator-to-metal transition in a biased VO 2 nanogap. We show that thermal dynamics allow the implementation of the basic neuron functionalities: activity, leaky integrate-and-fire, volatility and rate coding. This approach could enable neuromorphic hardware to take full advantage of the rapid advances in memristive synapses, allowing for much denser and complex neural networks.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Three-State Majority-vote Model on Scale-Free Networks and the Unitary Relation for Critical Exponents

Abstract We investigate the three-state majority-vote model for opinion dynamics on scale-free and regular networks. In this model, an individual selects an opinion equal to the opinion of the majority of its neighbors with probability 1 − q , and different to it with probability q . The parameter q is called the noise parameter of the model. We build a network of interactions where z neighbors are selected by each added site in the system, a preferential attachment network with degree distribution k − λ , where λ = 3 for a large number of nodes N . In this work, z is called the growth parameter. Using finite-size scaling analysis, we obtain that the critical exponents $$\beta /\bar{\nu }$$ β / ν ¯ and $$\gamma /\bar{\nu }$$ γ / ν ¯ associated with the magnetization and the susceptibility, respectively. Using Monte Carlo simulations, we calculate the critical noise parameter q c as a function of z for the scale-free networks and obtain the phase diagram of the model. We find that the critical exponents add up to unity when using a special volumetric scaling, regardless of the dimension of the network of interactions. We verify this result by obtaining the critical noise and the critical exponents for the two and three-state majority-vote model on cubic lattice networks.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Golden ratio in venation patterns of dragonfly wings

Abstract The vein pattern in insect wings allows this lightweight structure to carry multiple biological functions. Here, an investigation of the angular distribution of the vein struts in dragonfly wings revealed that the golden angle or golden ratio dominates the venation patterns. We find that the golden angle dominates the intervein angles in regions where thin veins and membranes demand strength reinforcement. A golden ratio partition method has thus been developed that explains a set of preferred intervein angles in distorted polygon-shaped venation cells throughout the venation pattern in dragonfly wings. These observations provide new evidence that the wing structure is spatially optimized, by the golden rule in nature, for supporting biomechanical functions of dragonfly wings.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Forces that control self-organization of chemically-propelled Janus tori

Control of the individual and collective behavior of self-propelled synthetic micro-objects has immediate application for nanotechnology, robotics, and precision medicine. Despite significant progress in the synthesis and characterization of self-propelled Janus (two-faced) particles, predictive understanding of their behavior remains challenging, especially if the particles have anisotropic forms. Here, by using molecular simulation, we describe the interactions of chemically-propelled microtori near a wall. The results show that a torus hovers at a certain distance from the wall due to a combination of gravity and hydrodynamic flows generated by the chemical activity. Moreover, electrostatic dipolar interactions between the torus and the wall result in a spontaneous tilt and horizontal translation, in a qualitative agreement with experiment. Simulations of the dynamics of two and four tori near a wall provide evidence for the formation of stable self-propelled bound states. Our results illustrate that self-organization at the microscale occurs due to a combination of multiple factors, including hydrodynamic, chemical, electrostatic and steric interactions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Variational quantum algorithms

Applications such as simulating complicated quantum systems or solving large-scale linear algebra problems are very challenging for classical computers, owing to the extremely high computational cost. Quantum computers promise a solution, although fault-tolerant quantum computers will probably not be available in the near future. Current quantum devices have serious constraints, including limited numbers of qubits and noise processes that limit circuit depth. Variational quantum algorithms (VQAs), which use a classical optimizer to train a parameterized quantum circuit, have emerged as a leading strategy to address these constraints. VQAs have now been proposed for essentially all applications that researchers have envisaged for quantum computers, and they appear to be the best hope for obtaining quantum advantage. Nevertheless, challenges remain, including the trainability, accuracy and efficiency of VQAs. Here we overview the field of VQAs, discuss strategies to overcome their challenges and finally, we highlight the exciting prospects for using them to obtain quantum advantage.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Quantum approximate multi-objective optimization

The goal of multi-objective optimization is to understand optimal trade-offs between competing objective functions by finding the Pareto front, that is, the set of all Pareto-optimal solutions, where no objective can be improved without degrading another one. Multi-objective optimization can be challenging classically, even if the corresponding single-objective optimization problems are efficiently solvable. Thus, multi-objective optimization represents a compelling problem class to analyze with quantum computers. Here we use a low-depth quantum approximate optimization algorithm to approximate the optimal Pareto front of certain multi-objective weighted maximum-cut problems. We demonstrate its performance on an IBM Quantum computer, as well as with matrix product state numerical simulation, and show its potential to outperform classical approaches.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

The SRG/eROSITA All-Sky Survey: Constraints on the structure growth from cluster number counts

Recent advancements in methods used in wide-area surveys have demonstrated the reliability of the number density of galaxy clusters as a viable tool for precision cosmology. Beyond testing the current cosmological paradigm, cluster number counts can also be used to investigate the discrepancies currently affecting cosmological measurements. In particular, cosmological studies based on cosmic shear and other large-scale structure probes routinely find a value for the amplitude of the fluctuations in the universe S​ 8 = σ 8 (Ω m /0.3) 0.5 smaller than the one inferred from the primary cosmic microwave background. In this work, we investigate this tension by measuring structure evolution across cosmic time as probed by the number counts of massive halos with the first SRG/eROSITA All-Sky Survey cluster catalog in the western Galactic hemisphere, complemented with the overlapping Dark Energy Survey Year-3, Kilo-Degree Survey, and Hyper Suprime-Cam data for weak lensing mass calibration, by implementing two different parameterizations and a model-agnostic method. In the first model, we measured the cosmic linear growth index as γ = 1.19 ± 0.21, which is in tension with the standard value of γ = 0.55 but in good statistical agreement with other large-scale structure probes. The second model is a phenomenological scenario in which we rescale the linear matter power spectrum at low redshift to investigate a potential reduction of structure formation, and it provided similar results. Finally, in a third strategy, we considered a standard ΛCDM cosmology, but we separated the cluster catalog into five redshift bins, measuring the cosmological parameters in each and inferring the evolution of the structure formation, finding hints of a reduction. Interestingly, the S ​8 value inferred from the number counts of the cluster eRASS1 when we add a degree of freedom to the matter power spectrum recovers the value inferred by cosmic shear studies. The observed reduction in the growth rate or systematic uncertainties associated with various measurements may account for the discrepancy in the S 8 ​ values suggested between cosmic shear probes and eROSITA cluster number counts and Planck CMB measurements.

cosmological parameters↗
DOE OSTI2025

The occurrence and impact of carbon-oxygen shell mergers in massive stars

In their final stages before undergoing a core-collapse supernova, massive stars may experience mergers between internal shells where carbon (C) and oxygen (O) are consumed as fuels for nuclear burning. This interaction, known as a C-O shell merger, can dramatically alter the internal structure of the star, leading to peculiar nucleosynthesis and potentially influencing the supernova explosion and the propagation of the subsequent supernova shock. Our understanding of the frequency and consequences of C-O shell mergers remains limited. This study aims to identify, for the first time, early diagnostics in the stellar structure that lead to C-O shell mergers in more advanced stages. We also assess their role in shaping the chemical abundances in the most metal poor stars of the Galaxy. We analyzed a set of 209 stellar evolution models available in the literature, with different initial progenitor masses and metallicities. We then compared the nucleosynthetic yields from a subset of these models with the abundances of odd-Z elements in metal-poor stars. We find that the occurrence of C-O shell mergers in stellar models can be predicted with a good approximation based on the outcomes of the central He burning phase, specifically, from the CO core mass (M CO ) and the 12 C central mass fraction (X C12 ): 90% of models with a C-O merger have X C12 <0.277 and M CO <4.90 M ⊙ , with average values of M CO = 4.02 M ⊙ and X C12 = 0.176. The quantities X C12 and M CO are indirectly affected from several stellar properties, including the initial stellar mass and metallicity. Additionally, we confirm that the Sc-rich and K-rich yields from models with C-O mergers would solve the long-standing underproduction of these elements in massive stars. Our results emphasize the crucial role of C-O shell mergers in enriching the interstellar medium, particularly in the production of odd-Z elements. This highlights the necessity of further investigations to refine their influence on presupernova stellar properties and their broader impact on Galactic chemical evolution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Re-evaluation of the TSL for Yttrium Hydride

Yttrium hydride (YH x ) is of interest as a high-temperature moderator material in advanced nuclear reactor systems because of its superior ability to retain hydrogen at elevated temperatures. Thermal neutron scattering laws (TSL) for hydrogen bound in yttrium hydride (H-YH 2 ) and yttrium bound in yttrium hydride (Y-YH 2 ) were previously evaluated by Naval Nuclear Laboratory using the ab initio approach and released in ENDF/B-VIII.0. In that work, density functional theory, incorporating the generalized gradient approximation (GGA) for the exchange-correlation energy, was used to simulate the face-centered cubic structure of YH 2 and calculate the interatomic Hellmann-Feynman forces for a 2×2×2 supercell containing 96 atoms. Lattice dynamics calculations using PHONON were used to determine the phonon density of states. The calculated phonon density of states for H and Y in YH 2 were then used to prepare H-YH 2 and Y-YH 2 TSL evaluations, in the incoherent approximation, using the LEAPR module of NJOY. In addition, elastic scattering was assumed to be incoherent for both H and Y. While the incoherent elastic scattering approximation is appropriate for H-YH 2 , it introduces an undesirable approximation for Y-YH 2 . In this work, we re-evaluate the TSL for Y-YH 2 using FLASSH (Full Law Analysis Scattering System Hub). Y-YH 2 is evaluated using the FLASSH generalized coherent elastic scattering capability in order to capture the Bragg peaks associated with the YH 2 crystal structure which were neglected in the prior NJOY-based evaluation due to limitations in LEAPR. An experimental approach to validate the Y-YH 2 TSL using neutron transmission measurements is discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Validation of Thermal Scattering Laws for Light Water at Elevated Temperatures with Diffusion Experiments

Light water is the most important neutron moderator in many reactor applications. Thermal neutron scattering kernels for H bound in H 2 O are represented by ENDF-format thermal scattering laws (TSLs) evaluated at discrete temperatures T . Nuclear data evaluations are conventionally validated utilizing a combination of critical benchmarks and experimental cross section data. Existing public critical benchmarks may be inadequate for testing water TSLs over the full range of T of interest in reactor applications, and experimental thermal scattering cross section data for water is sparse at elevated T . In this work, MC21 is used to simulate the decay of the equilibrium thermal neutron flux in light-water spheres of several radii. The fundamental-mode time eigenvalue $\mathcal{a}$ is calculated at 22 °C and 227 °C (at saturation pressure) for each sphere using three different H-H 2 O TSLs. Fitting $\mathcal{a}$ to a polynomial function of geometric buckling allows calculation of the thermal neutron diffusion length L . Extrapolation lengths are treated as a function of the transport mean free path and geometry. Experimental L data from 25 publications at 49 temperatures (from 10 °C to 295 °C), computed by several different time-dependent and space-dependent decay methods, is used to develop an empirical fit of L vs. T . The MC21-calculated L for the TSLs tested are within ±1% of the predicted values at 22 °C and 227 °C. This validation approach, which may be repeated at arbitrary T , constitutes an integral benchmark specific to and characterizing the detailed physics of the thermal scattering kernel applied.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Using Prompt Neutron Decay Constant Measurements to Obtain Additional Kinetics Information

Prompt neutron decay constant measurements are typically used to infer the prompt neutron lifetime in a chain-reacting system, but these measurements include much more information about the kinetics of the system. This work focuses specifically on the Rossi-α method of measuring the prompt neutron decay constant. The Rossi-α system of equations consists of five unknown parameters. With five unknowns and three equations, the system of equations solved by the Rossi-α measurement is under-determined. It has been postulated that a well posed problem can be formulated using independent measurement configurations on a single system. The external source strength, the effective delayed neutron fractions, and Diven’s parameter do not typically change for different configurations of a single experiment. If good detector geometry is chosen, the detector efficiency does not change for multiple measurements of a single experiment using Rossi-α measurements. By making these assumptions, there are only two unique variables. This knowledge can be used to solve for most of the unique parameters using two independent, unique Rossi-α measurements.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Resolved Resonance Region Analysis of 206 Pb, 207 Pb, and 208 Pb for Next Generation Lead-Cooled Fast Systems

Recently, great interest has been generated in using lead as a coolant for fast neutron systems and as a result it is important to investigate the ENDF/B-VIII.0 isotopic evaluations that comprise stable le ad. To this end, resonance parameters for 206 Pb, 207 Pb, and 208 Pb were re-evaluated because their resolved reso­nance regions extend beyond 0.5 MeV meaning resonance parameters used in reconstructing cross sections and elastic scattering angular distributions impact fast systems. The impact of resonance parameters is demon­strated by the differences between the evaluations in predicting experimental results from the fully modeled RPI Quasi-Differential Scattering Experiment via MCNPv6.2. In addition, MCNP KCODE calculations of lead-sensitive fast spectra critical benchmarks showed variations of keff on average of 400 pcm, caused solely from differences in elastic scattering angular distributions in 208 Pb. Re-evaluation entailed fitting data with the program SAMMY.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Advanced Modeling and Simulation Methods for Evaluation of Thermal Neutron Scattering Materials

With the rise of interest in thermal neutron scattering data for advanced reactor, criticality safety, and shielding applications, new experimental data are required for evaluation of new materials or for re-evaluation (or validations) of previously evaluated materials. New experimental data are evaluated in a three-step process: (1) computing the phonon characteristics, (2) computing the dynamic structure factor (DSF) from the data, and (3) using the experimental setup to simulate the experimental data. All three steps have challenges, ranging from the need for a sufficiently general material simulation code—a processing code that can compute the corresponding DSF—to having a detailed layout of the instrument/beamline/facility where the data were measured. Whereas phonon characteristics of materials can be computed using various methods (molecular dynamics, density functional theory, etc.), a high-fidelity computation of the DSF and the simulation of the experiment based on the DSF is vital to the accuracy of the evaluation. The latter two steps can be achieved by using the two corresponding code systems developed by instrument scientists at the Spallation Neutron Source (SNS) at Oak Ridge National Laboratory: (1) OCLIMAX, a program that calculates the dynamic structure factor from DFT and MD simulation results, and (2) MCViNE, a Monte Carlo neutron ray-tracing program designed to simulate neutron scattering experiments. Recently, polyethylene and yttrium hydride were measured at the Wide Angular-Range Chopper (ARCS) and SEQUOIA instrument stations of the SNS. These experiments are simulated using the density functional theory code, the Cambridge Serial Total Energy Package (CASTEP), to compute its phonon characteristics (eigenvalues/vectors and PDOS), which is then processed using OCLIMAX to yield the DSF, and finally the data at each instrument station are simulated by the MCViNE for comparison to the measured data for evaluation. For comparison to conventional evaluation methods, the scattering data processed from OCLIMAX are compared against those processed from the LEAPR module of NJOY, and the results from MCViNE simulations are compared against previously used simplified beamline models implemented in the Monte Carlo N-Particle (MCNP) code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Extended Lagrangian Born–Oppenheimer molecular dynamics using a Krylov subspace approximation

It is shown how the electronic equations of motion in extended Lagrangian Born–Oppenheimer molecular dynamics simulations can be integrated using low-rank approximations of the inverse Jacobian kernel. This kernel determines the metric tensor in the harmonic oscillator extension of the Lagrangian that drives the evolution of the electronic degrees of freedom. The proposed kernel approximation is derived from a pseudoinverse of a low-rank estimate of the Jacobian, which is expressed in terms of a generalized set of directional derivatives with directions that are given from a Krylov subspace approximation. The approach allows a tunable and adaptive approximation that can take advantage of efficient preconditioning techniques. The proposed kernel approximation for the integration of the electronic equations of motion makes it possible to apply extended Lagrangian first-principles molecular dynamics simulations to a broader range of problems, including reactive chemical systems with numerically sensitive and unsteady charge solutions. This can be achieved without requiring exact full calculations of the inverse Jacobian kernel in each time step or relying on iterative non-linear self-consistent field optimization of the electronic ground state prior to the force evaluations as in regular direct Born–Oppenheimer molecular dynamics. We note the low-rank approximation of the Jacobian is directly related to Broyden’s class of quasi-Newton algorithms and Jacobian-free Newton–Krylov methods and provides a complementary formulation for the solution of nonlinear systems of equations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

NWChem: Past, present, and future

Specialized computational chemistry packages have permanently reshaped the landscape of chemical sciences by providing tools to support and guide the experimental effort and for prediction of chemical and materials properties. In this regard, a special role has been played by electronic structure packages where complex chemical and materials processes can be modeled using first-principle-driven methodologies. Over the last few decades, the rapid development of computing technologies and a tremendous increase in computational power has offered a unique chance to study complex chemical transformations using sophisticated and predictive many-body techniques to describe correlated behavior of electrons in molecular and condensed phase systems at different levels of theory. In enabling these simulations, a critical role has been played by novel parallel algorithms capable of taking advantage of computational resources to address polynomial scaling of electronic structure methods. NWChem was among the first electronic structure codes that focused on delivering scalable performance for electronic structure simulations. Herein, we briefly review the NWChem suite of computational codes including its history, design principles, parallel tools, current capabilities, outreach and outlook.

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

An Entropy-Maximization Approach to Automated Training Set Generation for Interatomic Potentials

Machine learning-based interatomic potentials are currently garnering a lot of attention as they strive to achieve the accuracy of electronic structure methods at the computational cost of empirical potentials. Given their generic functional forms, the transferability of these potentials is highly dependent on the quality of the training set, the generation of which can be highly labor-intensive. Good training sets should at once contain a very diverse set of configurations while avoiding redundancies that incur cost without providing benefits. We formalize these requirements in a local entropy-maximization framework and propose an automated sampling scheme to sample from this objective function. We show that this approach generates much more diverse training sets than unbiased sampling and is competitive with hand-crafted training sets.

74 ATOMIC AND MOLECULAR PHYSICS↗