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

Spin liquid state in a rare-earth hyperkagome lattice

We report quantum fluctuations enhanced by frustration and subtle interplay between competing degrees of freedom offer an ideal ground to realize novel states with fractional quantum numbers in quantum materials that defy standard theoretical paradigms. Quantum spin liquid (QSL) is a highly entangled state wherein frustration-induced strong quantum fluctuations preclude symmetry-breaking phase transitions down to zero temperature without any order parameter. Experimental realizations of QSL in quantum materials with spin dimensionality greater than one is very rare. Here, we present our thermodynamic, nuclear magnetic resonance, muon spin relaxation, and inelastic neutron scattering studies of a rare-earth hyperkagome compound Li 3 Yb 3 Te 2 O 12 in which Yb 3+ ions constitute a three-dimensional spin lattice without any detectable disorder. Our comprehensive experiments evince neither signature of magnetic ordering nor spin freezing down to 38 mK that suggest the realization of dynamic liquid-like ground state in this antiferromagnet. The ground state of this material is interpreted by a low energy J eff = 1/2 degrees of freedom with short-range spin correlations. The present results demonstrate a viable basis to explore spin-orbit driven enigmatic correlated quantum states in a class of rare-earth-based three-dimensional frustrated magnets that may open avenues in theoretical and experimental search for spin liquids.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Probing entanglement in a 2D hard-core Bose–Hubbard lattice

Entanglement and its propagation are central to understanding many physical properties of quantum systems. Notably, within closed quantum many-body systems, entanglement is believed to yield emergent thermodynamic behaviour. However, a universal understanding remains challenging owing to the non-integrability and computational intractability of most large-scale quantum systems. Quantum hardware platforms provide a means to study the formation and scaling of entanglement in interacting many-body systems. Here we use a controllable 4 × 4 array of superconducting qubits to emulate a 2D hard-core Bose–Hubbard (HCBH) lattice. We generate superposition states by simultaneously driving all lattice sites and extract correlation lengths and entanglement entropy across its many-body energy spectrum. We observe volume-law entanglement scaling for states at the centre of the spectrum and a crossover to the onset of area-law scaling near its edges.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Data-driven studies of magnetic two-dimensional materials

We use a data-driven approach to study the magnetic and thermodynamic properties of van der Waals (vdW) layered materials. We investigate monolayers of the form A 2 B 2 X 6 , based on the known material Cr 2 Ge 2 Te 6 , using density functional theory (DFT) calculations and machine learning methods to determine their magnetic properties, such as magnetic order and magnetic moment. We also examine formation energies and use them as a proxy for chemical stability. We show that machine learning tools, combined with DFT calculations, can provide a computationally efficient means to predict properties of such two-dimensional (2D) magnetic materials. Our data analytics approach provides insights into the microscopic origins of magnetic ordering in these systems. For instance, we find that the X site strongly affects the magnetic coupling between neighboring A sites, which drives the magnetic ordering. Our approach opens new ways for rapid discovery of chemically stable vdW materials that exhibit magnetic behavior.

42 ENGINEERING↗
DOE OSTI2023

Quantifying uncertainty in analysis of shockless dynamic compression experiments on platinum. I. Inverse Lagrangian analysis

Absolute measurements of solid-material compressibility by magnetically driven shockless dynamic compression experiments to multi-megabar pressures have the potential to greatly improve the accuracy and precision of pressure calibration standards for use in diamond anvil cell experiments. Here, to this end, we apply characteristics-based inverse Lagrangian analysis (ILA) to 11 sets of ramp-compression data on pure platinum (Pt) metal and then reduce the resulting weighted-mean stress–strain curve to the principal isentrope and room-temperature isotherm using simple models for yield stress and Grüneisen parameter. We introduce several improvements to methods for ILA and quasi-isentrope reduction, the latter including calculation of corrections in wave speed instead of stress and pressure to render results largely independent of initial yield stress while enforcing thermodynamic consistency near zero pressure. More importantly, we quantify in detail the propagation of experimental uncertainty through ILA and model uncertainty through quasi-isentrope reduction, considering all potential sources of error except the electrode and window material models used in ILA. Compared to previous approaches, we find larger uncertainty in longitudinal stress. Monte Carlo analysis demonstrates that uncertainty in the yield-stress model constitutes by far the largest contribution to uncertainty in quasi-isentrope reduction corrections. We present a new room-temperature isotherm for Pt up to 444 GPa, with 1-sigma uncertainty at that pressure of just under ±1.2%; the latter is about a factor of three smaller than uncertainty previously reported for multi-megabar ramp-compression experiments on Pt. The result is well represented by a Vinet-form compression curve with (isothermal) bulk modulus K 0 = 270.3 ± 3.8 GPa, pressure derivative K$^{'}_{0}$= 5.66 ± 0.10, and correlation coefficient $R_{K_{0},K^{'}_{0}}$= –0.843.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Modeling dissociation of ionic liquids with electrolyte NRTL model

Dissociation of ionic liquids (ILs) is critical to electrochemical process performance by affecting ionic conductivity. Prior studies of thermodynamic modeling either assumed no dissociation or complete dissociation that were not consistent with experimental observations. This work presented an eNRTL model for three aqueous imidazolium-based IL systems: [Emim][EtSO 4 ], [Emim][TFA], and [Emim][TfO]. The model accurately correlates the IL dissociation and VLE data in the entire concentration range from pure IL to infinite dilution aqueous solution. The system non-ideality is dominated by the dissociation chemistry. Dissociation behavior is driven by strong short-range interaction of the molecular IL-water pair relative to the dissociated IL-water pair, consistent with prior studies indicating that molecular IL and water can form an energetically favorable complex. The interaction energy parameters calculated in this work qualitatively agrees with the hydrogen bonding ability of the anions previously reported. This research not only provides a model that well describes the IL dissociation but also demonstrates the physical significance of eNRTL model.

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

A Computational Framework to Accelerate the Discovery of Perovskites for Solar Thermochemical Hydrogen Production: Identification of Gd Perovskite Oxide Redox Mediators

A high-throughput computational framework to identify novel multinary perovskite redox mediators is presented, and this framework is applied to discover the Gd-containing perovskite oxide compositions Gd 2 BB'O 6 , GdA'B 2 O 6 , and GdA'BB'O 6 that split water. The computational scheme uses a sequence of empirical approaches to evaluate the stabilities, electronic properties, and oxygen vacancy thermodynamics of these materials, including contributions to the enthalpies and entropies of reduction, ΔH TR and ΔS TR . This scheme uses the machine-learned descriptor τ to identify compositions that are likely stable as perovskites, the bond valence method to estimate the magnitude and phase of BO 6 octahedral tilting and provide accurate initial estimates of perovskite geometries, and density functional theory including magnetic- and defect-sampling to predict STCH-relevant properties. Eighty-three promising STCH candidate perovskite oxides down-selected from 4392 Gd-containing compositions are reported, three of which are referred to experimental collaborators for characterization and exhibit STCH activity. Our results demonstrate that the high-throughput computational scheme described herein—which is used to evaluate Gd-containing compositions but can be applied to any multinary perovskite oxide compositional space(s) of interest—accelerates the discovery of novel STCH active redox mediators with reasonable computational expense.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Tuning the quantumness of simple Bose systems: A universal phase diagram

We present a comprehensive theoretical study of the phase diagram of a system of many Bose particles interacting with a two-body central potential of the so-called Lennard-Jones form. First-principles path-integral computations are carried out, providing essentially exact numerical results on the thermodynamic properties. The theoretical model used here provides a realistic and remarkably general framework for describing simple Bose systems ranging from crystals to normal fluids to superfluids and gases. The interplay between particle interactions on the one hand and quantum indistinguishability and delocalization on the other hand is characterized by a single quantumness parameter, which can be tuned to engineer and explore different regimes. Taking advantage of the rare combination of the versatility of the many-body Hamiltonian and the possibility for exact computations, we systematically investigate the phases of the systems as a function of pressure (P) and temperature (T), as well as the quantumness parameter. Here, we show how the topology of the phase diagram evolves from the known case of 4 He, as the system is made more (and less) quantum, and compare our predictions with available results from mean-field theory. Possible realization and observation of the phases and physical regimes predicted here are discussed in various experimental systems, including hypothetical muonic matter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Machine-learning accelerated geometry optimization in molecular simulation

Geometry optimization is an important part of both computational materials and surface science because it is the path to finding ground state atomic structures and reaction pathways. These properties are used in the estimation of thermodynamic and kinetic properties of molecular and crystal structures. This process is slow at the quantum level of theory because it involves an iterative calculation of forces using quantum chemical codes such as density functional theory (DFT), which are computationally expensive and which limit the speed of the optimization algorithms. It would be highly advantageous to accelerate this process because then one could do either the same amount of work in less time or more work in the same time. Here, we provide a neural network (NN) ensemble based active learning method to accelerate the local geometry optimization for multiple configurations simultaneously. We illustrate the acceleration on several case studies including bare metal surfaces, surfaces with adsorbates, and nudged elastic band for two reactions. In all cases, the accelerated method requires fewer DFT calculations than the standard method. In addition, we provide an Atomic Simulation Environment (ASE)-optimizer Python package to make the usage of the NN ensemble active learning for geometry optimization easier.

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

Pt Particles on a Dynamic TiO 2 Support in Near-Ambient Conditions−Disentangling Size, Pressure, and Support Effects

Platinum particles on reducible oxides are known to form complex and highly dynamic catalyst systems at elevated pressures and temperatures, often adopting active structures that differ from those found at room temperature and under ultrahigh vacuum (UHV). Here, we study the oxidation and structural evolution of subnanometer Pt clusters and nanoparticles supported on rutile TiO 2 (110) across an oxygen pressure range from UHV to 0.1 mbar, using near-ambient pressure X-ray photoelectron spectroscopy (NAP-XPS), scanning tunneling microscopy (STM) under UHV and NAP conditions, and low-energy ion scattering (LEIS). Our results reveal distinct differences in oxidation behavior and thermal stability between Pt nanoparticles and clusters, which are further modulated by the support stoichiometry and oxygen pressure. Small Pt clusters become oxidized even at room temperature but are susceptible to accelerated sintering in 0.1 mbar O 2 at elevated temperatures. In contrast, well-crystallized Pt nanoparticles on near-stoichiometric TiO 2 show weaker oxidation. On a reduced, defective TiO 2 support, Pt instead quickly becomes deeply buried by new titania layers, which are formed during support reoxidation. This process appears to result primarily from interactions of the support with the gas phase, unlike the classical, selflimited encapsulation that is induced by the strong metal−support interaction (SMSI). Finally, we address the full complexity of real catalysts in a direct side-by-side comparison of the single-crystalline model system with a Pt-loaded TiO 2 powder catalyst (P25). We conclude that the stoichiometry of the model supports must be carefully chosen and controlled to accurately reproduce the expected state of powder supports during redox reactions.

metal nanoparticles↗
DOE OSTI2021

Enhancing the formation of ionic defects to study the ice Ih/XI transition with molecular dynamics simulations

Ice Ih, the common form of ice in the biosphere, contains proton disorder. Its proton-ordered counterpart, ice XI, is thermodynamically stable below 72 K. However, the formation of ice XI is kinetically hindered, and experimentally it is obtained by doping with KOH. Doping creates ionic defects that promote the migration of protons and the associated change in proton configuration. In this article, we mimic the effect of doping with a bias potential that enhances the formation of ionic defects in molecular dynamics simulations. The recombination of the ions thus formed proceeds through fast migration of the hydroxide along hydrogen bond loops, providing a physical and expedite way to change the proton configuration. Here, a key ingredient of this approach is a machine learning potential trained with density functional theory data and capable of modelling molecular dissociation. We exemplify the usefulness of this idea by studying the order-disorder transition using an appropriate order parameter that distinguishes the proton environments in ice Ih and XI. We calculate the changes in free energy, enthalpy, and entropy associated with the transition. Our estimated entropy agrees with experiment within the error bars of the calculation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Thermodynamic behavior of correlated electron-hole fluids in van der Waals heterostructures

Abstract Coupled two-dimensional electron-hole bilayers provide a unique platform to study strongly correlated Bose-Fermi mixtures in condensed matter. Electrons and holes in spatially separated layers can bind to form interlayer excitons, composite Bosons expected to support high-temperature exciton condensates. The interlayer excitons can also interact strongly with excess charge carriers when electron and hole densities are unequal. Here, we use optical spectroscopy to quantitatively probe the local thermodynamic properties of strongly correlated electron-hole fluids in MoSe 2 /hBN/WSe 2 heterostructures. We observe a discontinuity in the electron and hole chemical potentials at matched electron and hole densities, a definitive signature of an excitonic insulator ground state. The excitonic insulator is stable up to a Mott density of ~0.8 × 10 12 cm −2 and has a thermal ionization temperature of ~70 K. The density dependence of the electron, hole, and exciton chemical potentials reveals strong correlation effects across the phase diagram. Compared with a non-interacting uniform charge distribution, the correlation effects lead to significant attractive exciton-exciton and exciton-charge interactions in the electron-hole fluid. Our work highlights the unique quantum behavior that can emerge in strongly correlated electron-hole systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Pressure-Induced Phase Transitions in Bilayer La 3 Ni 2 O 7

La 3 Ni 2 O 7 differs from expectations in part because at ambient pressure, it exists in two polymorphs: an unconventional crystal structure with alternating layers of single- and triple-layered nickel−oxygen octahedra (i.e., La 2 NiO 4 plus La 4 Ni 3 O 10 , with the so-called 1313 crystal structure) and a classical double-layer Ruddlesden−Popper phase (i.e., La 3 Ni 2 O 7 , with the so-called 2222 crystal structure, which dominates at high pressure). In this study, we report the pressure-dependent structural and electrical resistive properties of single crystals of the classical double-layered La 3 Ni 2 O 7 , grown at slightly elevated pressure using the floating-zone method. Structural characterization under pressures up to 15.4 GPa reveals a gradual transition from orthorhombic to tetragonal symmetry that is completed at 12−14 GPa. Additionally, we present pressure and field-dependent electrical resistance measurements up to 27.4 GPa, from which we construct a phase diagram. Although a resistive transition at 80 K is observed at high pressures, the resistive hallmarks of superconductivity are not.

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

Progress report on analytic and numerical studies of x-ray-induced impulse

In a variety of high-energy-density (HED) systems, x-rays of a given energy are used to generate shockwaves, bulk motion, and impulse in materials. As radiation energy is deposited within the material, the material is heated. The number of photons of a certain wavelength absorbed is determined by the spectral intensity of the radiation and the material’s opacity evaluated at that wavelength. This heating results in a pressure increase dictated by the material’s equation of state. Depending on the intensity of the absorbed radiation, it may also cause the material to change phase into a liquid, gas, or plasma. The increased pressure drives the heated surface layer to blow off, imparting impulse and sending a compressive wave into the bulk of the material. Additionally, the compression wave interacts with the solid boundary of the material, resulting in a tensile wave that may cause the material to spall. The impulse generated by the deposition of x-ray energy within the sample can be modeled using purely analytical methods, e.g. the Bethe, Bade, Averell, and Yost (BBAY) model. However, the blow-off process is rather complicated, and proper modeling efforts must account for material ejected by spallation, vaporization, jetting, and plasma ablation. For this reason, analytical models have an unclosed term describing the final energy of the blown-off material Ef(z). Prior modeling efforts have arbitrarily fixed this at some value or modeled it using a limiting set of thermodynamic assumptions. The work we are currently performing uses validated simulations using sophisticated photon transport, equation of state, and strength models/data to provide a fit for Ef(z) that is useful for predictive calculation of impulse. We will apply our methodology and show results for different materials and x-ray sources. This work is particularly useful for the design of experiments studying x-ray impulse generation. The present report is outlined as follows. Section 1.2 describes a series of HED experiments investigating x-ray-generated impulse in materials. Section 1.3 describes the computational approach we employ in this study, and presents validation results against the aforementioned experiments. Section 1.4 introduces analytical models for impulse generation, as well as our method for utilizing impulse from simulations to close the models. We discuss the concept of impulse-spectrum sensitivity, it’s application to uncertainty quantification, and derive a very useful analytical expression for it in 1.5. We then summarize recent progress in this project and discuss future work in section 1.6.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Multiple and nonlocal cation redox in Ca–Ce–Ti–Mn oxide perovskites for solar thermochemical applications

Modeling-driven design of redox-active off-stoichiometric oxides for solar thermochemical H 2 production (STCH) seldom has resulted in empirical demonstration of competitive materials. Here, we report the theoretical prediction and experimental evidence that the perovskite Ca 2/3 Ce 1/3 Ti 1/3 Mn 2/3 O 3 is synthesizable with high phase purity, stable, and has desirable redox thermodynamics for STCH, with a predicted average neutral oxygen vacancy (VO) formation energy, E v = 3.30 eV. Flow reactor experiments suggest potentially comparable or greater H 2 production capacity than recent promising Sr–La–Mn–Al and Ba–Ce–Mn metal oxide perovskites. Utilizing quantum-based modeling of a solid solution on both A and B sub-lattices, we predict the impact of nearest-neighbor composition on E v and determine that A-site Ce 4+ reduction dominates the redox-activity of Ca 2/3 Ce 1/3 Ti 1/3 Mn 2/3 O 3 . X-ray absorption spectroscopy measurements provide evidence that supports these predictions and reversible Ce 4+ -to-Ce 3+ reduction. Our models predict that Ce 4+ reduces even when it is not nearest-neighbor to the V O , suggesting that refinement of Ce stoichiometry has the possibility of further enhancing performance.

08 HYDROGEN↗
DOE OSTI2021

Direct numerical simulation of supercritical oxy-methane mixing layers with CO 2 substituted counterparts

Direct numerical simulations (DNS) of temporally developing, three dimensional, CH 4 /CO 2 , CH 4 /O 2 and CO 2 /O 2 mixing layers, are conducted at a supercritical pressure of 300 atm. To effectively model the supercritical regime, the employed formulation includes the compressible form of the governing equations, the cubic Peng-Robinson equation of state and a generalized formulation for heat and mass flux vectors derived from non-equilibrium thermodynamics and fluctuation theory. A linear inviscid stability analysis is also performed for each case, to determine its most unstable wavelength. Flow visualizations reveal the presence of high density gradient magnitude regions for all three mixing layers, with conditional averages indicating increased presence of heavier fluid species within these regions. No significant departures are observed from perfect gas behavior, with compressibility factors very close to unity for all three mixing cases. Applicability of presumed probability density function methods (PDF) is examined for the three supercritical mixing layers. An a priori analysis is also conducted to investigate various simplifying assumptions employed in modeling various subgrid scale (SGS) flux models. Two additional terms are identified in the large eddy simulations (LES) equations, the gradient of SGS contribution of pressure in the momentum equation and the gradient of SGS contribution of heat flux in energy equation, whose magnitudes are similar and comparable with their respective resolved terms. The performance of the scale similarity model to represent these additional terms is investigated. Lastly, the performance of Smagorinsky, gradient and scale similarity models are also investigated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Thermal transport exceeding bulk heat conduction due to nonthermal micro/nanoscale phonon populations

While classical size effects usually lead to a reduced effective thermal conductivity, we report here that nonthermal phonon populations produced by a micro/nanoscale heat source can lead to enhanced heat conduction, exceeding the prediction from Fourier's law. We study nondiffusive thermal transport by phonons at small distances within the framework of the Boltzmann transport equation (BTE) and demonstrate that the transport is significantly affected by the distribution of phonons emitted by the source. We discuss analytical solutions of the steady-state BTE for a source with a sinusoidal spatial profile, as well as for a three-dimensional Gaussian “hot spot,” and provide numerical results for single crystal silicon at room temperature. If a micro/nanoscale heat source produces a thermal phonon distribution, it gets hotter than that predicted by the heat diffusion equation; however, if the source predominantly produces low-frequency acoustic phonons with long mean free paths, it may get significantly cooler than that predicted by the heat equation, yielding an enhanced heat transport beyond bulk heat conduction.

42 ENGINEERING↗
DOE OSTI2023

Staggered scheme for the compressible fluctuating hydrodynamics of multispecies fluid mixtures

Here, we present a numerical formulation for the solution of nonisothermal, compressible Navier-Stokes equations with thermal fluctuations to describe mesoscale transport phenomena in multispecies fluid mixtures. The novelty of our numerical method is the use of staggered grid momenta along with a finite volume discretization of the thermodynamic variables to solve the resulting stochastic partial differential equations. The key advantages of the numerical scheme are that it significantly simplifies the discretization of diffusive and stochastic momentum fluxes into a more compact form, and it provides an unambiguous prescription of boundary conditions involving pressure. The staggered grid scheme more accurately reproduces the equilibrium static structure factor of hydrodynamic fluctuations in gas mixtures compared to a collocated scheme described previously by Balakrishnan et al. [Phys. Rev. E 89, 013017 (2014)1539-375510.1103/PhysRevE.89.013017]. The numerical method is tested for ideal noble gases mixtures under various nonequilibrium conditions, such as applied thermal and concentration gradients, to assess the role of cross-diffusion effects, such as Soret and Dufour, on the long-ranged correlations of hydrodynamic fluctuations, which are also more accurately reproduced compared to the collocated scheme. We numerically study giant nonequilibrium fluctuations driven by concentration gradients and fluctuation-driven Rayleigh-Taylor instability in gas mixtures. Wherever applicable, excellent agreement is observed with theory and measurements from the direct simulation Monte Carlo method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Phonon second harmonic generation in NaBr studied by inelastic neutron scattering and computer simulation

The phenomenon of second harmonic generation (SHG) was found for phonons in anharmonic NaBr by inelastic neutron scattering. The temperature dependence of this phonon SHG was measured from 300 K to 650 K. At 300 K the second harmonic (SH) is seen as a high-energy branch around 33 meV, nearly independent of $\overrightarrow{Q}$. The temperature effective potential (TDEP) method and classical molecular dynamics (MD) simulation with machine learning interatomic potential were able to reproduce the SH, and showed that SHG occurs with the flat transverse optical (TO) phonon branch. A classical model of a nonlinear medium explains the intensity and lifetime of the SH, compared to those of the TO modes. Also successful was a quantum model based on the Heisenberg-Langevin equation for interacting phonons coupled to a thermal bath, which also predicts a spectral distribution of the SH. In conclusion, the measured temperature dependence of the intensity of the second harmonic showed that it follows the Planck distribution of a one-phonon quasiparticle, and not two TO phonons.

36 MATERIALS SCIENCE↗