Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “equilibrium potential”

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 181 records · Page 10

Global minimization via classical tunneling assisted by collective force field formation

Simple elements interacting in networks can give rise to intricate emergent behaviors. Examples such as synchronization and phase transitions often apply in many contexts, as many different systems may reduce to the same effective model. Here, we demonstrate such a behavior in a model inspired by memristors. When weakly driven, the system is described by movement in an effective potential, but when strongly driven, instabilities cause escapes from local minima, which can be interpreted as an unstable tunneling mechanism. We dub this collective and nonperturbative effect a “Lyapunov force,” which steers the system toward the global minimum of the potential function, even if the full system has a constellation of equilibrium points growing exponentially with the system size. This mechanism is appealing for its physical relevance in nanoscale physics and for its possible applications in optimization, Monte Carlo schemes, and machine learning.

97 MATHEMATICS AND COMPUTING↗

Random matrix theory of the isospectral twirling

We present a systematic construction of probes into the dynamics of isospectral ensembles of Hamiltonians by the notion of Isospectral twirling, expanding the scopes and methods of ref. [1]. The relevant ensembles of Hamiltonians are those defined by salient spectral probability distributions. The Gaussian Unitary Ensembles (GUE) describes a class of quantum chaotic Hamiltonians, while spectra corresponding to the Poisson and Gaussian Diagonal Ensemble (GDE) describe non chaotic, integrable dynamics. We compute the Isospectral twirling of several classes of important quantities in the analysis of quantum many-body systems: Frame potentials, Loschmidt Echos, OTOCs, Entanglement, Tripartite mutual information, coherence, distance to equilibrium states, work in quantum batteries and extension to CP-maps. Moreover, we perform averages in these ensembles by random matrix theory and show how these quantities clearly separate chaotic quantum dynamics from non chaotic ones.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Out-of-equilibrium dynamics of a grid-like Fe( II ) spin crossover dimer triggered by a two-photon excitation

The application of two-photon excitation (TPE) in the study of light-responsive materials holds immense potential due to its deeper penetration and reduced photodamage. Despite these benefits, TPE has been underutilised in the investigation of the photoinduced spin crossover (SCO) phenomenon. Here, we employ TPE to delve into the out-of-equilibrium dynamics of a SCO Fe II dimer of the form [Fe II (HL) 2 ] 2 (BF 4 ) 4 ·2MeCN (HL = 3,5-bis{6-(2,2'-bipyridyl)}pyrazole). Optical transient absorption (OTA) spectroscopy in solution proves that the same dynamics take place under both one-photon excitation (OPE) and TPE. The results show the emergence of the photoinduced high spin state in less than 2 ps and with a lifetime of 147 ns. Time-resolved photocrystallography (TRXRD) reveals a single molecular reorganisation within the first 500 ps following TPE. Additionally, variable temperature single crystal X-ray diffraction (VTSCXRD) and magnetic susceptibility measurements confirm that the thermal transition is silenced by the solvent. While the results of the OTA and TRXRD utilising TPE are intriguing, the high pump fluencies required to excite enough metal centres to the high spin state may impair its practical application. Nonetheless, this study sheds light on the potential of TPE for the investigation of the out-of-equilibrium dynamics of SCO complexes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Coherent phonon flatband generated in GaAs/AlAs superlattices via layer-selective optical pumping

Flatbands, characterized by their dispersionless energy levels in electronic, magnetic, and phononic systems, hold substantial potential for advancements in electronics and quantum information processing. Most flatbands exist in thermal equilibrium and cannot be easily created or annihilated externally, limiting their flexibility as switchable knobs for use in microelectronics and quantum applications. In our work, we demonstrate the generation of a coherent phonon flatband in a GaAs/AlAs superlattice using 800 nm femtosecond laser pulses. This coherent phonon flatband does not correspond to a phonon eigenmode at equilibrium and exhibits strong coupling with two branches of coherently excited longitudinal phonon modes. With molecular dynamics simulations, we show more generally that the coherent phonon flatband can be induced by coherently and spatially modulated optical excitations of superlattice structures. Our results highlight a pathway for coherent phonon flatband creation in the time domain that can be generalized to various superlattice systems, potentially inspiring the realization of coherent flatband generation of other quasiparticles.

Condensed-matter physics↗

Computational Fermi Level Engineering and Doping-Type Conversion of Mg:Ga 2 O 3 via Three-Step Synthesis Process

Gallium oxide (Ga 2 O 3 ) is being actively explored for electronics that can operate at high power, temperature, and frequency as well as for deep-ultraviolet optoelectronics and other applications due to its ultra-wide bandgap (UWBG) and low projected fabrication cost of large-size and high-quality crystals. Efficient n-type doping of monoclinic beta-phase of Ga 2 O 3 has been achieved, but p-type doping faces fundamental obstacles due to compensation, deep acceptor levels, and the polaron transport mechanism of free holes. However, aside from the challenges of achieving p-type conductivity, plenty of opportunity exists to engineer the position of the Fermi level for improved design of Ga 2 O 3 -based devices. We use first-principles defect theory and defect equilibrium calculations to simulate a three-step growth-annealing-quench synthesis protocol for hydrogen-assisted Mg doping in ß-Ga 2 O 3 . The simulations take into account the gas phase equilibrium between H 2 , O 2 , and H 2 O, which determines the H chemical potential. We predict Ga 2 O 3 doping-type conversion to a net p-type regime after growth under reducing conditions in the presence of H2 followed by O-rich annealing, which is a similar process to Mg acceptor activation by H removal in GaN. For equilibrium annealing with re-equilibration of compensating O vacancies, there is an optimal temperature that maximizes the Ga 2 O 3 net acceptor density for a given Mg doping level; the acceptor density is further increased in the non-equilibrium annealing scenario without re-equilibration. After quenching to operating temperature, the Ga 2 O 3 Fermi level drops below mid-gap down to about 1.5 eV above the valence band maximum, creating a significant number of uncompensated neutral Mg Ga 0 acceptors. The resulting free hole concentration in Ga 2 O 3 is very low even at elevated operating temperature (~10 8 cm -3 at 400°C) due to the deep energy level of these Mg acceptors, and hole conductivity is further impeded by the polaron hopping mechanism. However, the Fermi-level reduction and suppression of free electron density in this doping-type converted (NA>ND) Ga 2 O 3 material are important for improved designs of Ga 2 O 3 electronic devices. These results illustrate the power of computational predictions not only for new materials but also for their synthesis science.

36 MATERIALS SCIENCE↗

Global biomass supply modeling for long-run management of the climate system

Bioenergy is projected to have a prominent, valuable, and maybe essential, role in climate management. However, there is significant variation in projected bioenergy deployment results, as well as concerns about the potential environmental and social implications of supplying biomass. Bioenergy deployment projections are market equilibrium solutions from integrated modeling, yet little is known about the underlying modeling of the supply of biomass as a feedstock for energy use in these modeling frameworks. We undertake a novel diagnostic analysis with ten global models to elucidate, compare, and assess how biomass is supplied within the models used to inform long-run climate management. With experiments that isolate and reveal biomass supply modeling behavior and characteristics (costs, emissions, land use, market effects), we learn about biomass supply tendencies and differences. The insights provide a new level of modeling transparency and understanding of estimated global biomass supplies that informs evaluation of the potential for bioenergy in managing the climate and interpretation of integrated modeling. For each model, we characterize the potential distributions of global biomass supply across regions and feedstock types for increasing levels of quantity supplied, as well as some of the potential societal externalities of supplying biomass. We also evaluate the biomass supply implications of managing these externalities. Finally, we interpret biomass market results from integrated modeling in terms of our new understanding of biomass supply. Overall, we find little consensus between models on where biomass could be cost-effectively produced and the implications. We also reveal model specific biomass supply narratives, with results providing new insights into integrated modeling bioenergy outcomes and differences. The analysis finds that many integrated models are considering and managing emissions and land use externalities of supplying biomass and estimating that environmental and societal trade-offs in the form of land emissions, land conversion, and higher agricultural prices are cost-effective, and to some degree a reality of using biomass, to address climate change.

09 BIOMASS FUELS↗

Downstream Separation of Formic Acid with Anion-Exchange Resin from Electrocatalytic Carbon Dioxide (CO 2 ) Conversion: Adsorption, Kinetics, and Equilibrium Modeling

The electrocatalytic carbon dioxide reduction reaction (CO 2 RR) that produces liquid formic acid (FA) is one potential route to minimize the CO 2 concentration in the atmosphere by utilizing CO 2 as a feedstock. However, one of the main disadvantages of this route is the high energy demand required for purification through conventional distillation. This method of purification has been shown to result in a large carbon footprint in life-cycle analysis (LCA). Therefore, finding an alternative purification technique that can integrate with an electrochemical CO 2 RR process with the potential to lower the carbon footprint would be of great interest. Under the experimental conditions of the CO 2 RR, FA remains in its basic form as the formate anion in the catholyte. Anion exchange resins with their basic properties can separate the formate anion from the catholyte mixture as an alternative to energy intensive distillation. In this work, three anion exchange resins with increasing basicity order: Amberlite IRA-96 < Amberlite IRA-910 < Ambersep 900 are tested to separate FA in the presence of aqueous solution of potassium bicarbonate catholyte with varying levels of resin and FA. Kinetics and equilibrium studies data for the FA adsorption are interpreted using several kinetics and isotherm models. The kinetics data fit well with a pseudo-first-order model at high initial FA concentrations and a pseudo-second-order model at low initial FA concentrations. The experimental data can be best explained with the Freundlich isotherm model. The fitted results show that strongly basic Ambersep 900 is more effective at separating FA, with a maximum adsorption capacity of 336.7 mg/g compared to Amberlite IRA-96 (275.2 mg/g) and Amberlite IRA-910 (209.2 mg/g) in neutral water. However, weakly basic Amberlite IRA-96 (110.8 mg/g) can separate FA more efficiently from potassium bicarbonate catholyte than Amberlite IRA-910 (99.9 mg/g) and Ambersep 900 (40.0 mg/g). As a result, competitive adsorption occurs in the presence of potassium bicarbonate catholyte, but the formic acid adsorption capacity is still sufficient for adequate separation.

20 FOSSIL-FUELED POWER PLANTS↗

Assessment of Performance of Density Functionals for Predicting Potential Energy Curves in Hydrogen Storage Applications

We report the availability of accurate computational tools for modeling and simulation is vital to accelerate the discovery of materials capable of storing hydrogen (H 2 ) under given parameters of pressure swing and temperature. Previously, we compiled the H2Bind275 data set consisting of equilibrium geometries and assessed the performance of 55 density functionals over this data set (Veccham, S. P.; Head-Gordon, M. J. Chem. Theory Comput. 2020, 16, 4963-4982). As it is crucial for computational tools to accurately model the entire potential energy curve (PEC), in addition to the equilibrium geometry, we extended this data set with 389 new data points to include two compressed and three elongated geometries along 78 PECs for H 2 binding, forming the H2Bind78 × 7 data set. By assessing the performance of 55 density functionals on this significantly larger and more comprehensive H2Bind78 × 7 data set, we identified the best performing density functionals for H 2 binding applications: PBE0-DH, ωB97X-V, ωB97M-V, and DSD-PBEPBE-D3(BJ). The addition of Hartree-Fock exchange improves the performance of density functionals, albeit not uniformly throughout the PEC. We recommend the usage of ωB97X-V and ωB97M-V density functionals as they offer good performance for both geometries and energies. In addition, we also identified B97M-V and B97M-rV as the best semilocal density functionals for predicting H 2 binding energy at its equilibrium geometry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mapping the kinetic evolution of metastable grain boundaries under non-equilibrium processing

The kinetic evolution of a multiplicity of metastable grain boundaries (GBs) under fast driving conditions are studied by atomistic modeling. Assisted with an enhanced statistical analysis, the energetic evolution of GBs over a broad metastability-temperature space is mapped out, wherein two distinct regimes—an ageing regime and a rejuvenating regime—are retrieved with high fidelity. By comparing the results under various conditions (e.g. random perturbations, isothermal annealing, and fast heating/cooling), it is shown that such ageing/rejuvenating mechanism map is universal, irrespective of the actual stimuli used to elicit the metastable GBs. Here, the ageing/rejuvenating phenomena are demonstrated to stem from the energy imbalance of uphill climbing and downhill dropping during sequential transitions in the system's potential energy landscape. Without the necessity of introducing free parameters, such model can reconcile experimentally measured hardness variation of nanocrystalline metals subjected to femto-second laser irradiation, and it therefore provides a novel perspective on achieving a plurality of interfacial states and facilitating previously inaccessible property regimes.

36 MATERIALS SCIENCE↗

Part I: Predicting performance of Purolite A532E resins for remediation of comingled contaminants in groundwater

Ion exchange (IX) resins are used in pump-and-treat (P&T) facilities to remove soluble groundwater contaminants. However, natural anions present at concentrations orders of magnitude higher than contaminants can compete for IX sites and impact resin lifecycles. Here, the Hanford Site’s 200 West Area P&T facility (Washington State, USA) was selected as a case study because it currently uses two IX resins: Purolite® A532E (A532E) to remove pertechnetate (TcO 4 - ) and DOWEX 21K (DOWEX) to remove uranium from groundwater. Nitrate (NO 3 - ), sulfate (SO 4 2- ), chloride (Cl - ), and carbonate (CO 3 2- ) anions have been identified to potentially compete for A532E and DOWEX IX sites. Hanford-relevant anion groundwater concentrations were used to design a series of laboratory-scale batch experiments to evaluate the impact of competing anions on resin performance and potential kinetic effects. These data are then modeled to obtain Cl--normalized equilibrium exchange coefficients (K) needed to predict IX resin performance. The work is presented in two parts, with IX performance evaluated for A532E in Part I and DOWEX in Part II. Part I results demonstrate that TcO 4 - uptake is not impacted by NO 3 - , SO 4 2- , Cl - , CO 3 2- (as HCO 3 - ) and U(VI) carbonate anions, with K TcO4-/Cl- > 4,000, likely due to the high selectivity of A532E trihexylammonium sites for the large, weakly hydrated TcO 4 - anion. Other anion K values were K NO3-/Cl- = 20, K SO4--/Cl- = 0.2, K HCO3-/Cl- = 0.09, K U/Cl- = 370–1000. These K values provide conservative parameters for predicting A532E performance, and demonstrate that, under these test conditions, A532E will remove TcO 4 - from current and future influent streams to meet groundwater treatment objectives.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

STEPs-SOL, a Peptoid Force Field Parameterization to Include Solvent Effects

As peptoids (N-substituted glycines) continue to gain popularity as a class of biomimetic polymers, the importance and demand for accurate force fields in molecular simulations also grow. Building on the vacuum-optimized Systematic and Extensible Force Field for Peptoids (STEPs) force field, here we present STEPs-SOL, a novel peptoid force field parametrization that effectively incorporates solvent effects to enhance the accuracy of peptoid simulations. The development of STEPs-SOL is based on the need for precise electrostatic modeling achieved through solvent-specific partial charge optimization. Here, our systematic approach significantly improves agreement with experimental measurements, reducing the mean absolute error in cis/trans ratio predictions (ΔG c/t ) by an average of 38% across multiple peptoid residues and solvent environments. This improved parametrization addresses computational challenges associated with nonbonded energies while maintaining a workflow that relies on high-level quantum mechanical data rather than depending solely on limited experimental equilibrium properties. By evaluating the effects of conformational bias in restrained electrostatic potential (RESP) charge generation and examining their impact on peptoid conformations in various solvents, we enhance our understanding of peptoid structural dynamics while providing a more accurate modeling framework.

force field↗

Nonequilibrium Solvent Effects during Photodissociation in Liquids: Dynamical Energy Surfaces, Caging and Chemical Identity

In the gas phase, potential energy surfaces can be used to provide insight into the details of photochemical reaction dynamics. In solution, however, it is unclear what potential energy surfaces, if any, can be used to describe even simple chemical reactions such as the photodissociation of a diatomic solute. In this paper, we use mixed quantum/classical (MQC) molecular dynamics (MD) to study the photodissociation of $Na$ $^+_2$ in both liquid Ar and liquid tetrahydrofuran (THF). We examine both the gas-phase potential surfaces and potentials of mean force (PMF), which assume that the solvent remains at equilibrium with the solute throughout the photodissociation process and show that neither resemble a nonequilibrium dynamical energy surface that is generated by taking the time integral of work. For the photodissociation of $Na$ $^+_2$ in liquid Ar, the dynamical energy surface shows clear signatures of solvent caging, and the degree of caging is directly related to the mass of the solvent atoms. For $Na$ $^+_2$ in liquid THF, local specific interactions between the solute and solvent lead to changes in chemical identity that create a kinetic trap that effectively prevents the molecule from dissociating. Here, the results show that nonequilibrium effects play an important role even in simple solution-phase reactions, requiring the use of dynamical energy surface to understand such chemical events.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Temperature and density dependent pair potential for deuterium under shock

Large-scale classical molecular dynamics (CMD) simulations naturally include the microscopic physics necessary for atomistic modeling of shock release at the ablator-fuel interface in an inertial confinement fusion (ICF) capsule. Here, the multi-megabar shocks utilized in ICF experiments can drive the deuterium fuel from ambient to electron volt temperatures (T) and multi-fold compression. Modeling interatomic interactions over such an extreme range of conditions is challenging for empirical bond order potentials. We generate a pair potential for deuterium with explicit temperature and mass density dependence from ab initio density functional theory molecular dynamics using the iterative Boltzmann inversion method. This potential accurately reproduces the radial distribution functions and pressures from DFT in CMD equilibrium simulations across a wide range of thermodynamic conditions, yet fails to return the expected Hugoniot relations when used in direct CMD shock simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Phase equilibrium of liquid water and hexagonal ice from enhanced sampling molecular dynamics simulations

We study the phase equilibrium between liquid water and ice Ih modeled by the TIP4P/Ice interatomic potential using enhanced sampling molecular dynamics simulations. Our approach is based on the calculation of ice Ih-liquid free energy differences from simulations that visit reversibly both phases. The reversible interconversion is achieved by introducing a static bias potential as a function of an order parameter. The order parameter was tailored to crystallize the hexagonal diamond structure of oxygen in ice Ih. Furthermore, we analyze the effect of the system size on the ice Ih-liquid free energy differences, and we obtain a melting temperature of 270 K in the thermodynamic limit. This result is in agreement with estimates from thermodynamic integration (272 K) and coexistence simulations (270 K). Since the order parameter does not include information about the coordinates of the protons, the spontaneously formed solid configurations contain proton disorder as expected for ice Ih.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Structure of the electrical double layer at the ice–water interface

The surface of ice in contact with water contains sites that undergo deprotonation and protonation and can act as adsorption sites for aqueous ions. Therefore, an electrical double layer should form at this interface and existing models for describing the electrical double layer at metal oxide–water interfaces should be able to be modified to describe the surface charge, surface potential, and ionic occupancy at the ice–water interface. I used a surface complexation model along with literature measurements of the zeta potential of ice in brines of various strength and pH to constrain equilibrium constants. I then made predictions of ion site occupancy, surface charge density, and partitioning of counterions between the Stern and diffuse layers. The equilibrium constant for cation adsorption is more than 5 orders of magnitude larger than the other constants, indicating that this reaction dominates even at low salinity. Deprotonated OH sites are predicted to be slightly more abundant than dangling O sites, consistent with previous work. Surface charge densities are on the order of ±0.001 C/m2 and are always negative at the moderate pH values of interest to atmospheric and geophysical applications (6–9). In this pH range, over 99% of the counterions are contained in the Stern layer. This suggests that diffuse layer polarization will not occur because the ionic concentrations in the diffuse layer are nearly identical to those in the bulk electrolyte and that electrical conduction and polarization in the Stern layer will be negligible due to reduced ion mobility.

Daigle, Hugh (ORCID:0000000260628321)↗

Overcoming time and complexity limitations in molecular dynamics investigations of equilibrium melting

Abstract A hybrid Monte-Carlo molecular-dynamics method for determining solidus and liquidus compositions in multicomponent systems is presented that overcomes both the time limitations in conventional molecular dynamics that prevent the evolution of distinct solid and liquid compositions via diffusion and the complexity challenge that prevents use of thermodynamic assessment in systems of many components. This hybrid method is validated in the Cu–Ni system against an independent assessment of solidus and liquidus compositions based on the regular solution model. Strategies for efficient mapping of different phase diagrams, based on the thermodynamic parameter T 0 , the temperature at which two phases of the composition X 0 have equal free energies, are presented and then demonstrated for the copper-nickel fully miscible system and the gold–silicon eutectic system. A calculation of the solidus and liquidus sampled during the equilibrium melting of equiatomic CrMnFeCoNi is performed, indicating that this method has potential to be extended to the study of many component alloys.

Au-Si↗

Heat transport in liquid water from first-principles and deep neural network simulations

In this work, we compute the thermal conductivity of water within linear response theory from equilibrium molecular dynamics simulations, by adopting two different approaches. In one, the potential energy surface (PES) is derived on the fly from the electronic ground state of density functional theory (DFT) and the corresponding analytical expression is used for the energy flux. In the other, the PES is represented by a deep neural network (DNN) trained on DFT data, whereby the PES has an explicit local decomposition and the energy flux takes a particularly simple expression. By virtue of a gauge invariance principle, established by Marcolongo, Umari, and Baroni, the two approaches should be equivalent if the PES were reproduced accurately by the DNN model. We test this hypothesis by calculating the thermal conductivity, at the GGA (PBE) level of theory, using the direct formulation and its DNN proxy, finding that both approaches yield the same conductivity, in excess of the experimental value by approximately 60%. Besides being numerically much more efficient than its direct DFT counterpart, the DNN scheme has the advantage of being easily applicable to more sophisticated DFT approximations, such as meta-GGA and hybrid functionals, for which it would be hard to derive analytically the expression of the energy flux. We find in this way that a DNN model, trained on meta-GGA (SCAN) data, reduces the deviation from experiment of the predicted thermal conductivity by about 50%, leaving the question open as to whether the residual error is due to deficiencies of the functional, to a neglect of nuclear quantum effects in the atomic dynamics, or, likely, to a combination of the two.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

General Method for Determining Redox Potentials without Electrolyte

A novel method to determine redox potentials without electrolyte is presented here. The method is based on a new ability to determine the dissociation constant, $K°_d$, for ion pairs formed between any radical anion and any inert electrolyte counterion. These dissociation constants can be used to determine relative shifts of redox potential as a function of electrolyte concentration, connecting referenced potentials determined with electrochemistry (with 0.1 M electrolyte) to electrolyte-free values. Pulse radiolysis created radical anions enabling determination of equilibrium constants for electron transfer between anions of donor and acceptor molecules as a function of electrolyte concentration in THF. The measurements determined “composite equilibrium constants”, $K_{eqC}$, which contain information about the dissociation constant for the electrolyte cations, X + , with the radical anions of both the donor, $K°_d$(D –• ,X + ) and the acceptor, $K°_d$(A –• ,X + ). Dissociation constants were obtained for a selection of radical anions with tetrabutylammonium (TBA + ). The electrolyte was found to shift the reduction potentials of small molecules 1-methylpyrene and trans-stilbene by close to +130 mV whereas oligo-fluorenes and polyfluorenes experienced shifts of only (+25 ± 6) mV due to charge delocalization weakening the ion pair. These shifts for reduction of aromatic hydrocarbon molecules are smaller than shifts of +232 and +451 mV seen previously for benzophenone radical anion with TBA + and Na + respectively where the charge on the radical anion is localized largely on one C=O bond, thus forming a more tightly bound ion pair.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗