Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Classical thermodynamics”

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 433 records · Page 24

Topological kinetic crossover in a nanomagnet array

Ergodic kinetics, which are critical to equilibrium thermodynamics, can be constrained by a system’s topology. Here, we studied a model nanomagnetic array in which such constraints visibly affect the behavior of the magnetic moments. In this system, magnetic excitations connect into thermally active one-dimensional strings whose motion can be imaged in real time. At high temperatures, our data showed the merging, breaking, and reconnecting of strings, resulting in the system transitioning between topologically distinct configurations. Below a crossover temperature, the string motion is dominated by simple changes in length and shape. In this low-temperature regime, the system is energetically stable because of its inability to explore all possible topological configurations. This kinetic crossover suggests a generalizable conception of topologically broken ergodicity and limited equilibration.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Continuous symmetry breaking in a two-dimensional Rydberg array

Spontaneous symmetry breaking underlies much of our classification of phases of matter and their associated transitions. The nature of the underlying symmetry being broken determines many of the qualitative properties of the phase; this is illustrated by the case of discrete versus continuous symmetry breaking. Indeed, in contrast to the discrete case, the breaking of a continuous symmetry leads to the emergence of gapless Goldstone modes controlling, for instance, the thermodynamic stability of the ordered phase. Here, we realize a two-dimensional dipolar XY model that shows a continuous spin-rotational symmetry using a programmable Rydberg quantum simulator. We demonstrate the adiabatic preparation of correlated low-temperature states of both the XY ferromagnet and the XY antiferromagnet. In the ferromagnetic case, we characterize the presence of a long-range XY order, a feature prohibited in the absence of long-range dipolar interaction. In conclusion, our exploration of the many-body physics of XY interactions complements recent works using the Rydberg-blockade mechanism to realize Ising-type interactions showing discrete spin rotation symmetry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Lifetime of Almost Strong Edge-Mode Operators in One-Dimensional, Interacting, Symmetry Protected Topological Phases

Almost strong edge-mode operators arising at the boundaries of certain interacting one-dimensional symmetry protected topological phases with Z 2 symmetry have infinite temperature lifetimes that are nonperturbatively long in the integrability breaking terms, making them promising as bits for quantum information processing. We extract the lifetime of these edge-mode operators for small system sizes as well as in the thermodynamic limit. For the latter, a Lanczos scheme is employed to map the operator dynamics to a one-dimensional tight-binding model of a single particle in Krylov space. We suggest this model to be that of a spatially inhomogeneous Su-Schrieffer-Heeger model with a hopping amplitude that increases away from the boundary, and a dimerization that decreases away from the boundary. We associate this dimerized or staggered structure with the existence of the almost strong mode. Thus, the short time dynamics of the almost strong mode is that of the edge mode of the Su-Schrieffer-Heeger model, while the long time dynamics involves decay due to tunneling out of that mode, followed by chaotic operator spreading. Furthermore, we show that competing scattering processes can lead to interference effects that can significantly enhance the lifetime.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Ethylene oxidation on unpromoted silver catalysts: Reaction pathway and selectivity analysis using DFT calculations

Chemical conversions in catalytic partial oxidation processes of light hydrocarbons are responsible for the production of numerous industrial chemicals, plastics, and intermediates. These processes are relatively expensive to perform, and are typically operated at high thermodynamic inefficiency, so the development of novel, highly efficient catalysts would prove to be very cost effective. Herein, our study focused on surface catalytic mechanisms of the ethylene oxide (EO) formation process. Periodic plane-wave Density Functional Theory (DFT) methods were used to analyze related reaction mechanisms on the Ag(111) surface facet with low coverage. Energetic changes of related species and pathways were calculated. Key surface species are identified to suggest the factors for the observed selectivity during EO formation. Our results are consistent with previous kinetic modeling efforts in the literature which did not employ DFT analysis. Lastly, our study demonstrates how fundamental theoretical investigations and multi-scale modeling techniques are currently impacting the advancement of rational catalyst design and microkinetic modeling techniques in the light hydrocarbon processing industry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effects of cluster expansion on the locations of phase transition boundary as a first step to quantify uncertainty in first principles statistical mechanics framework

Predicting phase diagrams from first principle calculations eliminates the need of tedious experimental trials and errors. Fully automating first principle phase diagram calculations without any sort of human intervention has been a long daunting task and troubling scientific communities for decades. This grand problem remains not fully resolved, largely due to the vastly high-dimensional parameter space associated with density functional theory, cluster expansion, lattice Monte Carlo, and the substantial uncertainty propagating through a set of complex simulations. As a first step to tackle this grand problem, we reported a first demonstration of how sensitive phase boundary locations can be to various cluster expansion fittings and input DFT training data. To the best knowledge of the authors, this study reported the first ever attempt to quantify uncertainty in first principles statistical mechanics framework. In addition, a semi-automated phase transition detection algorithm has been devised in this paper to deal with the associated statistical errors and uncertainties from Monte Carlo method and its predecessor DFT calculations and cluster expansions. This algorithm has been applied in a classical cluster expansion Mg-Cd binary alloy system to detect phase transitions at various locations of phase diagram using different cluster expansions to demonstrate its predictive power and quantify uncertainties. The results suggested that using Chebyshev basis function shifted the transition locations toward the dilute solid solution phase and expanded the phase stability of concentrated ordered phases to wider composition range. The addition of perturbed defect configurations into training data set lowered the order-disorder transition temperature to be away from true transition temperature, suggesting the transition nature is indeed configurational disorder dominated rather than defect assisted. Finally, we found that the weighting has negligible effect on the transition locations, except for the case of using Chebyshev basis function to fit all non-weighted configurations that can be susceptible to Monte Carlo sampling hysteresis.

36 MATERIALS SCIENCE↗

Simulating Effective QED on Quantum Computers

In recent years simulations of chemistry and condensed materials has emerged as one of the preeminent applications of quantum computing, offering an exponential speedup for the solution of the electronic structure for certain strongly correlated electronic systems. To date, most treatments have ignored the question of whether relativistic effects, which are described most generally by quantum electrodynamics (QED), can also be simulated on a quantum computer in polynomial time. Here we show that effective QED, which is equivalent to QED to second order in perturbation theory, can be simulated in polynomial time under reasonable assumptions while properly treating all four components of the wavefunction of the fermionic field. In particular, we provide a detailed analysis of such simulations in position and momentum basis using Trotter-Suzuki formulas. We find that the number of T -gates needed to perform such simulations on a 3 D lattice of n s sites scales at worst as O ( n s 3 / ϵ ) 1 + o ( 1 ) in the thermodynamic limit for position basis simulations and O ( n s 4 + 2 / 3 / ϵ ) 1 + o ( 1 ) in momentum basis. We also find that qubitization scales slightly better with a worst case scaling of O ~ ( n s 2 + 2 / 3 / ϵ ) for lattice eQED and complications in the prepare circuit leads to a slightly worse scaling in momentum basis of O ~ ( n s 5 + 2 / 3 / ϵ ) . We further provide concrete gate counts for simulating a relativistic version of the uniform electron gas that show challenging problems can be simulated using fewer than 10 13 non-Clifford operations and also provide a detailed discussion of how to prepare multi-reference configuration interaction states in effective QED which can provide a reasonable initial guess for the ground state. Finally, we estimate the planewave cutoffs needed to accurately simulate heavy elements such as gold.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Scaling Turbulent Combustion Fields in Explosions

We considered the topic of explosions from spherical high-explosive (HE) charges. We studied how the turbulent combustion fields scale. On the basis of theories of dimensional analysis by Bridgman and similarity theories of Sedov and Barenblatt, we found that all fields scaled with the explosion length scale r0. This included the blast wave, the mean and root mean squared (RMS) profiles of thermodynamic variables, combustion variables, velocities, vorticity, and turbulent Reynolds stresses. This was a consequence of the formulation of the problem and our numerical method, which both satisfied the similarity conditions of Sedov. We performed numerical simulations of 1 g charges and 1 kg charges; the solutions were identical (within roundoff error) when plotted in scaled variables. We also explored scaling laws related to three-phase pyrotechnic explosions. We show that although the scaling formally broke down, the fireball still essentially scaled with the explosion length scale r0. However, the discrete Lagrange particles (DLP) (phase 2) and the heterogeneous continuum model (HCM) of the DLP wakes (phase 3) did not scale with r0, and mean and RMS profiles could differ by a factor of 10 in some regions. This was because the DLP particles and wakes introduced an additional scale that broke the similarity conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Transfer learning as a method to reproduce high-fidelity non-local thermodynamic equilibrium opacities in simulations

Simulations of high-energy density physics often need non-local thermodynamic equilibrium opacity data. These data, however, are expensive to produce at relatively low fidelity. It is even more so at high fidelity such that the opacity calculations can contribute 95 % of the total computation time. This proportion can even reach large proportions. Neural networks can be used to replace the standard calculations of low-fidelity data, and the neural networks can be trained to reproduce artificial, high-fidelity opacity spectra. In this work, it is demonstrated that a novel neural network architecture trained to reproduce high-fidelity krypton spectra through transfer learning can be used in simulations. Further, it is demonstrated that this can be done while achieving a relative per cent error of the peak radiative temperature of the hohlraum of approximately 1 % to 4 % while achieving a 19.4 $\times$ speed up.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Many-Body Perspective of Nuclear Quantum Effects in Aqueous Clusters

Nuclear quantum effects play an important role in the structure and thermodynamics of aqueous systems. Here, by performing a many-body expansion with nuclear-electronic orbital (NEO) theory, we show that proton quantization can give rise to significant energetic contributions for many-body interactions spanning several molecules in single-point energy calculations of water clusters. Although zero-point motion produces a large increase in energy at the one-body level, nuclear quantum effects serve to stabilize higher-order molecular interactions. These results are significant because they demonstrate that nuclear quantum effects play a nontrivial role in many-body interactions of aqueous systems. Our approach also provides a pathway for incorporating nuclear quantum effects into water potential energy surfaces. The NEO approach is advantageous for many-body expansion analyses because it includes nuclear quantum effects directly in the energies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Relativistic spin hydrodynamics with torsion and linear response theory for spin relaxation

Using the second law of local thermodynamics and the first-order Palatini formalism, we formulate relativistic spin hydrodynamics for quantum field theories with Dirac fermions, such as QED and QCD, in a torsionful curved background. We work in a regime where spin density, which is assumed to relax much slower than other non-hydrodynamic modes, is treated as an independent degree of freedom in an extended hydrodynamic description. Spin hydrodynamics in our approach contains only three non-hydrodynamic modes corresponding to a spin vector, whose relaxation time is controlled by a new transport coefficient: the rotational viscosity. We study linear response theory and observe an interesting mode mixing phenomenon between the transverse shear and the spin density modes. We propose several field-theoretical ways to compute the spin relaxation time and the rotational viscosity, via the Green-Kubo formula based on retarded correlation functions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Materials Design Directions for Solar Thermochemical Water Splitting

The sustainable, economical production of molecular hydrogen is a crucial component of a net zero-greenhouse-gas-emissions future. Solar thermochemical water splitting (STWS) offers a renewable route to hydrogen with the potential to help decarbonize several industries, including transportation, manufacturing, mining, metals processing, and electricity generation, as well as provide sustainable hydrogen as a chemical feedstock. STWS uses high temperatures generated from concentrated sunlight or other sustainable means for high-temperature heat to produce hydrogen and oxygen from steam. For example, in its simplest form of a two-step thermochemical cycle, a redox-active metal oxide is heated to ≈1700-2000 K, driving off molecular oxygen while producing oxygen vacancies in the material. The reduced metal oxide then cools (ideally with the extracted heat recuperated for re-use) and, in a separate step, comes into contact with steam, which reacts with oxygen vacancies to produce molecular hydrogen while recovering the original state of the metal oxide. Despite its promising use of the entire solar spectrum to split water thermochemically, the current estimated cost of hydrogen produced via STWS is ≈4-6× the U.S. Department of Energy (DOE) Hydrogen Shot target value of $1/kg. One contributing approach to bridging this cost gap is the design of new materials with improved thermodynamic properties to enable higher efficiencies. The state-of-the-art (SOA) redox-active metal oxide for STWS is ceria (CeO 2 ), due to its close to optimal, although too high, oxygen vacancy formation enthalpy and large configurational and electronic entropy of reduction. However, ceria requires high operating temperatures and its efficiency is insufficient. Therefore, efforts to increase the efficiency of STWS cycles have focused on further optimizing oxygen vacancy formation enthalpies and augmenting the reduction entropy via substitution or doping and materials discovery schemes. Examples of the latter include the perovskites BaCe 0.25 Mn 0.75 O 3 and (Ca,Ce)(Ti,Mn)O 3 . These efforts and others have revealed intuitive chemical principles for the efficient and systematic design of more effective materials, such as the strong correlation between the enthalpies of crystal bond dissociation and solid-state cation reduction with the enthalpy of oxygen vacancy formation, as well as configurational entropy augmentation via the coexistence of two or more redox-active cation sublattices. The purpose of this chapter is to prepare the reader with an up-to-date account of STWS redox-active materials, both the SOA and promising newcomers, as well as to provide chemically intuitive strategies for improving their cycle efficiencies through materials design – in conjunction with ongoing efforts in reactor engineering and gas separations – to reach the cost points for commercial viability. First, we will introduce the thermodynamics of STWS using a two-step, metal-oxide, thermochemical cycle with economics in mind. We also will compare the pros and cons of processes that do or do not involve phase changes. Second, we will describe the qualities that make ceria the SOA STWS redox-active material, as well as its limitations. Third, we will survey some of the most promising candidates to date in the search for materials to supplant ceria, emphasizing the post-ternary, metal-oxide-perovskite alloys. Lastly, we will enumerate and discuss the following materials design directions for STWS redox-active materials: crystal reduction potentials as a proxy for oxygen vacancy formation enthalpies, engineering the electronic and configurational entropy of reduction via f-shells and simultaneous redox, and vetting materials stability via temperature-dependent phase diagrams and melting-point prediction.

08 HYDROGEN↗

Learning crystal field parameters using convolutional neural networks

We present a deep machine learning algorithm to extract crystal field (CF) Stevens parameters from thermodynamic data of rare-earth magnetic materials. The algorithm employs a two-dimensional convolutional neural network (CNN) that is trained on magnetization, magnetic susceptibility and specific heat data that is calculated theoretically within the single-ion approximation and further processed using a standard wavelet transformation. We apply the method to crystal fields of cubic, hexagonal and tetragonal symmetry and for both integer and half-integer total angular momentum values J J of the ground state multiplet. We evaluate its performance on both theoretically generated synthetic and previously published experimental data on CeAgSb _2 2 , PrAgSb _2 2 and PrMg _2 2 Cu _9 9 , and find that it can reliably and accurately extract the CF parameters for all site symmetries and values of J J considered. This demonstrates that CNNs provide an unbiased approach to extracting CF parameters that avoids tedious multi-parameter fitting procedures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Experimental Study of the Fundamental Properties of Warm Dense Mixtures

The aim of the proposed research was to provide a multi-scale study of the properties of warm dense hydrocarbons, by studying the thermodynamic properties through the equation of state and the microscopic properties by x-ray scattering. Understanding the fundamental properties of warm dense mixtures is an intellectual challenge due to the complexity of the system. Unlike liquids or gases, where constituent particle interactions occur through collisions or bonding between valence electrons, in these strongly coupled systems the atoms will be partially ionized and compressed so tightly together that interactions between the inner core electrons can play a role in the systems chemistry. Advances in computational capabilities and development of new theoretical models have been used to predict the properties of mixtures but there is currently no experimental data of the fundamental interaction between the particles in mixtures to test these predictions against. The proposed research was to experimentally investigate the interaction of the different species in hydrocarbons by measuring the compressibility of substances with different carbon and hydrogen ratios and the complexity of the microscopic interactions through elastic and inelastic x-ray scattering. To study the bulk properties of warm dense hydrocarbons we established the Warm Dense Matter Research Laboratory (WDMRL) in the Institute for Shock Physics at Washington State University. The goal was for experiments in the WDMRL to determine shock loading conditions of interest in the hydrocarbon mixtures. We planned on using Hugoniot EOS measurements and a range of carbon and hydrogen concentrations to determine conditions when the EOS of the mixture varied significantly from that of the classical mixing model. Even though shock transit measurements through aluminum foils suggested pressures upto 400GPa, the experiments in the WDMRL were unsuccessful in getting usable shockwave compression data above 100GPa in polystyrene which was below the pressure of interest for hydrocarbon mixtures (>200GPa). To complete the project, we tested a technique using x-ray phase contrast imaging to map the location of tracer layers in a test sample of polycarbonate to record the material motion in dynamically compressed samples. A technique that will be useful for future warm dense matter experiments. These experiments used <200nm gold layers in polycarbonate samples to measure the material velocity and shock speed using x-ray phase contrast imaging. These results were compared to continuum surface measurements performed at ISP and show that the tracer layer technique can measure hydrodynamic properties accurately in dynamically compressed materials.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modeling Equilibrium Solid–Liquid Interfaces under Effective Constant Chemical Potential Using Machine Learning Interatomic Potentials

The chemical potential (μ) of species in solution is essential for understanding various chemical processes at interfaces. Molecular dynamics (MD) simulations, constrained by fixed compositions, cannot maintain constant chemical potential with reference to a targeted concentration or chemical potential under nonequilibrium or dynamic conditions, as solute species can migrate to the interface and deplete (or enrich) the bulk due to solute-interface interactions. In this study, we introduce a simple and computationally efficient approach named iterative quasi-constant chemical potential molecular dynamics (iqCμMD) simulation, which helps simulate targeted molar concentrations of species in solution. iqCμMD overcomes the limitations of conventional MD by adjusting the number of species in the solution to reach a target bulk concentration (chemical potential), which allows simulation of the interface under the bulk conditions comparable to experiment. We demonstrate our approach using machine learning interatomic potential (MLIP)-based MD simulations of the Na 2 SO 4,aq –graphene interface, and to show the transferability of our approach, we also perform classical force field-based MD simulations of NaCl aq –air and NaCl aq –graphite interfaces, which produce comparable results to previous CμMD simulations. Our results also show that the iqCμMD approach efficiently achieves the desired bulk ion concentration within two iterations, and by utilizing MLIPs, we can achieve converged results using relatively small-scale simulations compared to previous CμMD simulations. By combining iqCμMD with MLIP-driven simulations, solid–liquid interfaces can be modeled under an effective constant chemical potential with DFT-level accuracy. Here, we show that iqCμMD offers a robust and simple computational framework for constant chemical potential simulations, as its only requirement is to be able to converge interfacial simulations with a measurable bulk region.

Chemical structure↗

Exciton-driven antiferromagnetic metal in a correlated van der Waals insulator

Abstract Collective excitations of bound electron-hole pairs—known as excitons—are ubiquitous in condensed matter, emerging in systems as diverse as band semiconductors, molecular crystals, and proteins. Recently, their existence in strongly correlated electron materials has attracted increasing interest due to the excitons’ unique coupling to spin and orbital degrees of freedom. The non-equilibrium driving of such dressed quasiparticles offers a promising platform for realizing unconventional many-body phenomena and phases beyond thermodynamic equilibrium. Here, we achieve this in the van der Waals correlated insulator NiPS 3 by photoexciting its newly discovered spin–orbit-entangled excitons that arise from Zhang-Rice states. By monitoring the time evolution of the terahertz conductivity, we observe the coexistence of itinerant carriers produced by exciton dissociation and a long-wavelength antiferromagnetic magnon that coherently precesses in time. These results demonstrate the emergence of a transient metallic state that preserves long-range antiferromagnetism, a phase that cannot be reached by simply tuning the temperature. More broadly, our findings open an avenue toward the exciton-mediated optical manipulation of magnetism.

36 MATERIALS SCIENCE↗

Towards accurate orbital-free simulations: A generalized gradient approximation for the noninteracting free energy density functional

For orbital-free ab initio molecular dynamics, especially on systems in extreme thermodynamic conditions, we provide the first pseudo-potential-adapted generalized gradient approximation (GGA) functional for the non-interacting free energy. This is achieved by systematic finite-temperature extension of our recent LKT ground state non-interacting kinetic energy GGA functional (Phys. Rev. B 98, 041111(R) (2018)). We test the performance of the new functional first via static lattice calculations on crystalline aluminum and silicon. Then we compare deuterium equation of state results against both path-integral Monte Carlo and conventional (orbital-dependent) Kohn-Sham results. The new functional, denoted LKTF, outperforms the previous best semi-local free energy functional, VT84F (Phys. Rev. B 88, 161108(R) (2013)), and provides modestly faster simulations. Finally, we discuss subtleties of identification of kinetic and entropic contributions to non-interacting free-energy functionals obtained by extension from ground state orbital-free kinetic energy functionals.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Computational investigation of water glasses using machine-learning potentials

The molecular origins of water’s anomalous properties have long been a subject of scientific inquiry. The liquid–liquid phase transition hypothesis, which posits the existence of distinct low-density and high-density liquid states separated by a first-order phase transition terminating at a critical point, has gained increasing experimental and computational support and offers a thermodynamically consistent framework for many of water’s anomalies. However, experimental challenges in avoiding crystallization near the postulated liquid–liquid critical point have focused attention to water’s canonical glassy states: low-density and high-density amorphous ice. Here, we use two Deep Potential machine-learning models, trained on the Strongly Constrained and Appropriately Normed density functional and the highly accurate Many-Body Polarizable potential, to conduct an investigation of water’s glassy phenomenology based on quantum mechanical calculations. Despite not being explicitly trained on amorphous ices, both models accurately capture the structure and transformation of the water glasses, including their interconversion along different thermodynamic paths. Isobaric quenching of liquid water at various pressures generates a continuum of intermediate amorphous ices and density fluctuations increase near the liquid–liquid critical pressure. The glass transition temperatures of the amorphous ices produced at different pressures exhibit two distinct branches, corresponding to low-density and high-density amorphous ice behaviors, consistent with experiment and the liquid–liquid transition hypothesis. Extrapolating transformation pressures from isothermal compressions to experimental compression rates brings our simulations into excellent agreement with data. Our findings demonstrate that machine-learning potentials trained on equilibrium phases can effectively model nonequilibrium glassy behavior and pave the way for studying long-timescale, out-of-equilibrium processes with quantum mechanical accuracy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Crystal and magnetic structure of polar oxide HoCrWO 6

In this study, polar magnetic oxide HoCrWO 6 is synthesized and its crystal structure, magnetic structure, and thermodynamic properties are investigated. HoCrWO 6 forms the polar crystal structure (space group Pna 2 1 (#33)) due to the cation ordering of W 6+ and Cr 3+ . There is an antiferromagnetic transition at T N = 24.5 K along with the magnetic entropy change (~5 J.Kg. –1 K –1 at 70 kOe). Neutron diffraction measurement indicates that both Cr and Ho sublattices are ordered with the moment of 2.32(5)μ B and 8.7(4)μ B at 2 K, respectively. While Cr forms A-type collinear antiferromagnetic (AFM) structure with magnetic moment along the b axis, Ho sublattice orders in a non-coplanar AFM arrangement. A comparison with isostructural DyFeWO 6 and DyCrWO 6 indicates that the magnetic structure of this family of compounds is controlled by the presence or absence of e g electrons in the transition metal sublattice.

42 ENGINEERING↗