Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase shift keying”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Mechanistic insights to drive catalytic hydrogenation of formamide intermediates to methanol via deaminative hydrogenation

Amine-promoted hydrogenation of CO 2 to methanol typically proceeds via a formamide intermediate when amines are used as additives or if the hydrogenation is performed in carbon capture solvents. The catalysts used for the hydrogenation of the formamide intermediate dictate the selectivity of the products formed: 1) Deoxygenative hydrogenation (C–O bond cleavage) resulting in N-methylation of amine and deactivation of the solvent, 2) Deaminative hydrogenation (C–N bond cleavage) resulting in formation of methanol and regeneration of the solvent. To date, catalytic reductions of CO 2 with amine promoters suffer from poor selectively for methanol which we attribute to the limiting formamide intermediate, though to date, the conditions that favor C–N cleavage have yet to be fully understood. To better understand the reactivity of the formamide intermediates, a range of heterogenous catalysts were used to study the hydrogenation of formamide. Well-known gas phase CO 2 hydrogenation catalysts catalyze the hydrogenation of formamide to N-methyl product via C–O bond cleavage. However, the selectivity can be readily shifted to selective C–N bond cleavage by addition of an additive with sufficient basicity for both homogenous and heterogeneous catalytic systems. The base additive shifts the selectivity by deprotonating a hemiaminal intermediate formed in situ during the formamide hydrogenation. This prevents dehydration process leading to N-methylated product, which is a key capture solvent deactivation pathway that hinders amine use in carbon capture, utilization, and storage (CCUS). The findings from this study provide a roadmap on how to improve the selectivity of known heterogenous catalysts, enabling catalytic reduction of captured CO 2 to methanol.

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

Bonding and Acidity of the Formal Hydride in the Prototypical Au 9 (PPh 3 ) 8 H 2+ Nanocluster

The role of hydrogen atoms as surface ligands on metal nanoclusters is of profound importance but remains difficult to directly study. While hydrogen atoms often appear to be incorporated formally as hydrides, evidence suggests that they donate electrons to the cluster's delocalized superatomic orbitals and may consequently behave as acidic protons that play key roles in synthetic or catalytic mechanisms. Here we directly test this assertion for the prototypical Au 9 (PPh 3 ) 8 H 2+ nanocluster, formed by addition of a hydride to the well-characterized Au 9 (PPh 3 ) 8 3+ . Using gas-phase infrared spectroscopy, we were able to unambiguously isolate Au 9 (PPh 3 ) 8 H 2+ and Au 9 (PPh 3 ) 8 D 2+ , revealing an Au–H stretching mode at 1528 cm –1 that shifts to 1038 cm –1 upon deuteration. This shift is greater than the maximum expected for a typical harmonic potential, suggesting a potential governing cluster-H bonding that has some square-well character consistent with the hydrogen nucleus behaving as a metal atom in the cluster core. Complexing this cluster with very weak bases reveals a redshift of 37 cm –1 in the Au–H vibration, consistent with those typically seen for moderately acidic groups in gas phase molecules and providing an estimate of the acidity of Au 9 (PPh 3 ) 8 H 2+ , at least with regard to its surface reactivity.

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

Liquid Salt Combined-Cycle Pilot Plant Design

The work described in this report is responsive to the Office of Fossil Energy program ‘Energy Storage for Fossil Power Generation.’ This Phase I report has been prepared by Pintail Power LLC, with support from Nexant ECA, Electric Power Research Institute (EPRI) and Southern Company Services as a deliverable for the U.S. Department of Energy for NETL Award DE-FE-00320016. The Liquid Salt Combined Cycle™ (LSCC™) technology provides large-scale energy storage integrated with Fossil Electric Generating Units (FEGUs) to meet critical needs in the energy transition by providing: • the lowest cost large-scale storage for time-shifting of renewable energy, • superior fuel efficiency to reduce GHGs from dispatchable resources, • flexible capacity and ramping to balance variability of wind and solar resources, • essential grid stability services to assure reliability of a low-carbon grid. The LSCC approach: • employs equipment that has already been proven in utility service, • uses safe, non-toxic, non-degrading, perpetual-life storage medium, • leverages and repurposes existing FEGU assets, • expands the value stack of energy storage to reduce market, financing, and commodity risks. Pintail Power has developed the LSCC technology to meet the need for reliable, efficient, and cost-effective integration of Variable Renewable Energy (VRE) into a low-carbon electric grid by coupling proven thermal energy storage with proven gas turbines, steam turbines, and heat transfer equipment. This novel approach is intended to address the key issues facing the grid and operators of renewable and fossil generating units including: • Overgeneration and curtailment of renewables, • Need for fast ramping dispatchable resources, • Improved efficiency and flexibility of fossil units, • Additional peaking capacity to support electrification of transportation and heating, • Provision of reliability services to support high penetration of VRE, especially synchronous inertia and fast frequency response. A Technology Readiness assessment by EPRI confirmed that LSCC technology consists of commercially proven hardware used in industrial and utility applications. Although the novel LSCC approach has not yet been demonstrated as a complete system, interfaces between major components have been conservatively specified. A Phase III pilot is planned to demonstrate equipment integration and operation. The patented innovation is removal of the evaporator section from the exhaust heat recovery system, with the evaporation performed by stored energy in a separate steam generator. This arrangement couples renewable and fossil power generation via long-duration energy storage to deliver cost, performance, and operational synergies, including superior charging and discharging flexibility, reduced fuel consumption and lower CO 2 emissions compared to conventional Combined Cycle Power Plants, and low-cost, large-scale energy storage. The LSCC technology is composed of proven equipment integrated with gas turbine exhaust heat in a novel system. During charging, electric heaters raise the salt temperature as it flows from the Cold Salt Tank to the Hot Salt Tank. During discharging, hot salt produces steam from feedwater that is heated with gas turbine exhaust, which also superheats steam to drive a steam turbine. LSCC technology can be added to any combustion-turbine to integrate renewable energy, provide needed grid services, and increase the value of fossil electric generating units based on the technology’s following attributes: • Long-duration storage enables time-shifting of VRE to avoid curtailment and impairment of renewable assets. • Long storage duration combined with fast-charging capability increases arbitrage opportunities by storing more energy when the price is low and discharging more hours when the price is high. • Long storage duration allows resource adequacy to be supplied across multiple days to increase reliability and reduce risk. • The stored energy reduces fuel heat rate and GHG emissions, and increases merit, so the LSCC dispatches earlier and longer to increase the plant’s capacity factor and asset value. • The stored energy enables pre-heating and startup of the steam cycle, without operating the gas turbine, to enable fast startup and ramping when dispatched for discharge. • The steam turbine can operate without the gas turbine so it can provide valuable synchronous inertia during charging without consuming fuel. • Fast frequency response and regulation services can be provided during charging using solid-state heater and pump controls to vary the charge power input in response to grid signals. • The LSCC system can be configured for resilience including black start, islanded/micro-grid operation, and even self-recharging of storage using either gas turbine power or gas turbine exhaust heat. The commercialization plan is to add LSCC technology to existing simple cycle gas turbine power plants with the 50MW GE LM6000 aero-derivative gas turbine as the reference design basis. A Techno-economic assessment of the reference design evaluated the benefits (Levelized Avoided Cost of Energy) and costs (Levelized Cost of Energy). The plant definition included all major systems and budgetary vendor quotes. Pintail Power and NexantECA developed the overall cost estimate for the LSCC plant up to the total plant cost level, following the DOE-NETL cost estimate guidelines at AACE Class 3 (-20%/+30%). This includes the equipment cost, bulk material, direct and indirect labor costs to arrive at the bare erected cost. Engineering costs are factored from the BEC and added to it to arrive at the EPC cost. Process and project contingencies were then factored from the EPC cost and rolled-up to yield the total plant cost of $\$$184 million for 1746 MWh of discharge electricity. • At $\$$105/kWh, the reference plant costs less than any of the Energy Storage Systems evaluated by PNNL in 2020 for the Energy Storage Grand Challenge. Operations and Maintenance cost estimates were scaled from combined cycle practice, assuming that the LSCC unit was co-located with and sharing some labor expense with other units, to arrive at $\$$2.2 million per year. Plant economics were evaluated using prices from the ERCOT Day-Ahead Market for calendar year 2019 (excluding the market disruptions from the COVID pandemic and the February 2020 deep freeze event). Assuming economic dispatch in the ERCOT Day-Ahead market, the reference plant capacity factor would have discharged for 2777 hours at 91.9 MW, a 31.66% capacity factor, with a marginal cost of $\$$25.59/MWh, and a LACE of $\$$82.41/MWh. Fixed charges were calculated according to EIA guidelines to arrive at an LCOE of $\$$83.48. The benefit-to-cost ratio of 0.99 suggests that the reference plant would have been cost-effective and competitive in the market. EPRI interviewed selected utilities to gauge the need for, applicability of and interest in the LSCC system. Several utilities are currently managing increased load growth along with the inclusion of increasing levels of renewable generation, putting pressure on conventional generation by requiring increased turndown requirements and ultimately lower capacity factors. All of the utilities interviewed have CO 2 reduction targets in the 2030-2050 timeframe that will severely limit the participation of fossil generation and require better utilization of carbon free generation. While there is limited opportunity for storage in the current markets, the utilities interviewed stated that there will be a substantial need for long duration energy storage in the future given the expected trends. Utilizing an energy storage system will generally be preferred over new gas capacity in some cases, with the capabilities of the LSCC system being a potential option for retrofit to existing simple cycle gas turbine units, allowing them to deliver greater participation in the market with lower carbon intensity. A technology gap assessment and technology maturation plan identified a pilot-scale demonstration as the final step before commercialization. Key gaps to be addressed during the Phase II FEED (Front-End Engineering Design) are component selection and design, commissioning procedures, and operational procedures and the control system for LSCC charging and discharging. The project team has been expanded to include Wood Group PLC as EPC. The proposed Phase II work leads to a pilot-scale engineering demonstration (TRL 6) to be conducted at Southern Company’s Plant Rowan, where the prototype system will perform “all the functions that will be required of the operational system.” The proposed pilot will facilitate commercialization (TRL-9) by scale-up to utility-scale systems integrated with peaking GTs or directly to facility scale systems using industrial GTs. The conceptual design for the pilot plant focuses on the novel integration aspects of LSCC technology. A slipstream of gas turbine exhaust will feed a waste heat recovery unit coupled to a molten salt steam generator heated by stored energy. The pilot is intended to demonstrate all key operating modes of the LSCC technology during charging, discharging and standby. The pilot equipment will be approximately one-seventh scale of the LM6000 commercial target and is expected to have commercial off-ramp potential for facility-scale applications.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2024

Quantitative phase imaging by gradient retardance optical microscopy

Abstract Quantitative phase imaging (QPI) has become a vital tool in bioimaging, offering precise measurements of wavefront distortion and, thus, of key cellular metabolism metrics, such as dry mass and density. However, only a few QPI applications have been demonstrated in optically thick specimens, where scattering increases background and reduces contrast. Building upon the concept of structured illumination interferometry, we introduce Gradient Retardance Optical Microscopy (GROM) for QPI of both thin and thick samples. GROM transforms any standard Differential Interference Contrast (DIC) microscope into a QPI platform by incorporating a liquid crystal retarder into the illumination path, enabling independent phase-shifting of the DIC microscope's sheared beams. GROM greatly simplifies related configurations, reduces costs, and eradicates energy losses in parallel imaging modalities, such as fluorescence. We successfully tested GROM on a diverse range of specimens, from microbes and red blood cells to optically thick (~ 300 μm) plant roots without fixation or clearing.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

Superionic conduction in solid polymer electrolytes – decoupling ion transport from segmental relaxation

Solvent-free, solid polymer electrolytes (SPEs) are promising candidates for next-generation, electrochemical energy storage systems due to their potential to enhance safety and performance, enable flexible device architectures, and streamline manufacturing processes. Conventional SPEs suffer from limited ionic conductivity due to the strong coupling between ion transport and (generally slow) polymer segmental relaxation. The realization of superionic conduction in SPEs, in which ions move faster than the structural relaxation of the polymers, requires a shift in design principles to promote this type of decoupled ion motion. In this perspective, we discuss how polymer architecture, ion–ion correlations, and ion–polymer interactions can unlock superionic behavior. We highlight several key design features, such as crystallinity, bulky side groups, high molecular weight, and percolating ionic aggregation, with a focus on creating low-barrier transport pathways in various polymer systems. We also demonstrate opportunities to combine polymer chemistry and data science through high-throughput and automated screening approaches to reveal how phase behavior, ion dynamics, and ionic interactions govern transport, thereby potentially enabling data-driven discovery of superionic polymer electrolyte materials.

Yang, Mengying [Univ. of Delaware, Newark, DE (Uni↗
DOE OSTI2026

Secondary-Side Active Rectifier Synchronization and Control in LCC-LC Compensated Inductive Power Transfer Systems

This paper proposes a communication-less synchronization and output power control strategy for LCC-LC compensated inductive power transfer systems with a secondaryside active rectifier. By replacing the passive rectifier with an active bridge, the proposed architecture eliminates the need for a separate secondary-side DC-DC converter and enables direct regulation of battery power. The key challenge in such a system is achieving robust synchronization of the secondary active bridge with the primary inverter without access to primary-side signals. To address this issue, the secondary-side variables are analyzed and the current through the secondary series inductor is identified as a load-independent synchronization variable with a fixed phase relationship to the primary excitation. A fixed-frequency second order generalized (SOGI)-based phase locked loop (PLL) is used to estimate the switching phase and frequency, and a phase-shift-based output power controller is developed using a phasor-transformer model of the active bridges. A modified control formulation is further introduced to remove phase-dependent loop-gain variation and simplify controller design. Simulation results validate the proposed method under both fixed and varying switching-frequency conditions, demonstrating successful frequency tracking from 82 kHz to 88 kHz and accurate output power regulation over a wide operating range. The proposed method offers a compact and effective solution for high-power IPT systems by enabling secondary-side synchronization and control without wireless communication or an additional DC-DC conversion stage.

Fernandes, Arnold Anthony [ORNL] (ORCID:0009000631↗
DOE OSTI2022

Bridging knowledge gaps in liquid- and vapor-fed CO 2 electrolysis through active electrode area

Increased use of gas diffusion electrodes for CO 2 electroreduction widens the experimental phase space that was previously inaccessible using foil electrodes, raising fundamental questions over the impacts of key variables that translate between liquid- and vapor-fed CO 2 electrolysis systems. This work focuses on studying the interplay of current-potential profiles and electrochemically active surface area (ECSA) by implementing a Cu nanoflower catalyst morphology. The results show decreased overpotentials as much as 460 and 174 mV for foil and gas diffusion electrodes, respectively, while maintaining or improving multi-carbon product current density. Furthermore, these overpotential shifts and product activities normalized by ECSA lead to current-potential relationships akin to those of the Tafel description, which are found through a continuum model to be useful for describing the roughness dependence for both liquid- and vapor-fed systems. This analysis establishes a holistic approach for establishing catalyst design criteria to improve materials development for CO 2 .

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Strong-Weak Duality via Jordan-Wigner Transformation: Using Fermionic Methods for Strongly Correlated su(2) Spin Systems

The Jordan-Wigner transformation establishes a duality between $su(2)$ and fermionic algebras. We present qualitative arguments and numerical evidence that when mapping spins to fermions, the transformation makes strong correlation weaker, as demonstrated by the Hartree-Fock approximation to the transformed Hamiltonian. This result can be rationalized in terms of rank reduction of spin shift terms when transformed to fermions. Conversely, the mapping of fermions to qubits makes strong correlation stronger, complicating its solution when one uses qubit-based correlators. The presence of string operators poses challenges to the implementation of quantum chemistry methods on classical computers, but these can be dealt with using established techniques of low computational cost. Here, our proof of principle results for XXZ and J1-J2 Heisenberg (in 1D and 2D) indicate that the JW transformed fermionic Hamiltonian has reduced complexity in key regions of their phase diagrams, and provides a better starting point for addressing challenging spin problems.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

New Methods for Predicting Non-Born-Oppenheimer Chemistry

Current methods for modeling non-adiabatic molecular dynamics face fundamental limitations when treating geometric phase effects: quantum mechanical phenomena where nuclear wavepackets acquire phase shifts when encircling conical intersections. Existing approaches either neglect these effects entirely or rely on potential energy surfaces arising from the Born-Oppenheimer approximation, which introduce artificial singularities and can overestimate geometric phase contributions. This project developed a new theoretical framework based on exact factorization (XF) methods to overcome these limitations. We derived mathematical formulations for hybrid quantum-classical XF dynamics that selectively treat critical nuclear degrees of freedom quantum mechanically while propagating others classically. This approach addresses the computational intractability that has previously limited exact methods to toy systems. Key innovations include a new approach to systematically identifying nuclear coordinates requiring quantum treatment, as well as novel implementation strategies that interface with existing quantum chemistry codes. The project also developed a proof-of-concept code for treating Jahn-Teller systems and creation of educational materials on non-adiabatic dynamics geared at the graduate level. The theoretical framework developed will enable future systematically improvable calculations of nuclear quantum effects in realistic molecular systems, filling a critical gap in non-adiabatic dynamics methods. This foundation supports future development of predictive tools for designing energy-relevant photochemical processes where quantum coherence effects may be exploited to control reaction outcomes.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2023

Heat pump water heater enhanced with phase change materials thermal energy storage: Modeling study

A promising solution to improve the first hour rating (FHR) of a heat pump water heater (HPWH) involves employing a secondary tank which contains phase change material (PCM) capsules. To better understand the influence of PCM thermal storage on the HPWH operational performance, a dynamic model was developed to simulate and analyze the behavior of a newly developed HPWH technology that incorporates PCM storage into a standard HPWH to optimize key parameters such as the uniform energy factor and FHR. Mathematical models of several key components of the proposed HPWH-PCM integrated thermal energy storage (TES) system, e.g., water heater tank, PCM TES tank, evaporator, compressor, and expansion valve, have been elaborated. Also, a model-based control co-simulation platform was developed to integrate a embed PCM storage HPWH dynamic model with a control model for better supporting control design, analysis, verification, and validation. The model accuracy has been validated through comparing simulation results with lab test results, with a mean average percentage error of <5.5% for most of the selected performance variables. In addition, using the developed co-simulation platform, the demand response control strategy was studied to evaluate the load flexibility of the combined HPWH-PCM storage system by shifting system power usage to outside of the 3.5 h peak load period.

25 ENERGY STORAGE↗
DOE OSTI2024

DG-IMEX method for a two-moment model for radiation transport in the $\mathscr{O}$($v$/$c$) limit

Here, we consider neutral particle systems described by moments of a phase-space density and propose a realizability-preserving numerical method to evolve a spectral two-moment model for particles interacting with a background fluid moving with nonrelativistic velocities. The system of nonlinear moment equations, with special relativistic corrections to $\mathscr{O}$($v$/$c$), expresses a balance between phase-space advection and collisions and includes velocity-dependent terms that account for spatial advection, Doppler shift, and angular aberration. The model is conservative for the correct $\mathscr{O}$($v$/$c$) Eulerian-frame number density and is consistent, to $\mathscr{O}$($v$/$c$), with Eulerian-frame energy and momentum conservation. This model is closely related to the one promoted by Lowrie et al. and similar to models currently used to study transport phenomena in large-scale simulations of astrophysical environments. The proposed numerical method is designed to preserve moment realizability, which guarantees that the moments correspond to a nonnegative phase-space density. The realizability-preserving scheme consists of the following key components: (i) a strong stability-preserving implicit-explicit (IMEX) time-integration method; (ii) a discontinuous Galerkin (DG) phase-space discretization with carefully constructed numerical uxes; (iii) a realizability-preserving implicit collision update; and(iv) a realizability-enforcing limiter. In time integration, nonlinearity of the moment model necessitates solution of nonlinear equations, which we formulate as fixed-point problems and solve with tailored iterative solvers that preserve moment realizability with guaranteed global convergence. We also analyze the simultaneous Eulerian-frame number and energy conservation properties of the semi-discrete DG scheme and propose a "spectral redistribution" scheme that promotes Eulerian-frame energy conservation. Through numerical experiments, we demonstrate the accuracy and robustness of this DG-IMEX method and investigate its Eulerian-frame energy conservation properties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Hofmann Stability Charts Revisited for PIP-II: From Classical Theory to Assumption-Free and ML-Driven Maps

The Hofmann stability chart remains a standard for visualizing parametric resonances in space-charge–dominated linacs, but its use typically relies on non-oscillatory Vlasov dispersion relations with simplifying assumptions (continuous focusing, KV phase space, linear optics, limited transverse–longitudinal coupling). We revisit the chart for the PIP-II linac along three tracks. (1) We reproduce the conventional maps in the (νz/νx, νx/ν0x) plane for relevant εz/εx, providing a validated reference. (2) We remove key assumptions by deriving stability surfaces directly from multi-particle tracking with realistic lattice discreteness, RF defocusing, solenoid/quad optics, and bunched-beam dynamics; local tunes and early-time growth rates are estimated from envelope oscillations and projected to the same coordinates. These assumption-reduced maps recover the canonical stopbands while revealing shifts and broadenings driven by tune modulation, non-KV distributions, and transverse–longitudinal coupling at PIP-II intensities. (3) We train a compact machine-learning surrogate that emulates the growth surface from zero-current optics, tune depression, emittance ratio, bunching factor, and selected lattice descriptors, enabling rapid scans and online working-point selection. We compare the three representations on representative PIP-II sections and discuss implications for commissioning guard bands, resonance avoidance, and routine operations.

Pathak, Abhishek [Fermilab] (ORCID:000000021704208↗
DOE OSTI2023

On the Electride Nature of Na–hP4

Early quantum mechanical models suggested that pressure drives solids towards free-electron metal behavior where the ions are locked into simple close-packed structures. The prediction and subsequent discovery of high-pressure electrides (HPEs), compounds assuming open structures where the valence electrons are localized in interstitial voids, required a paradigm shift. Our quantum chemical calculations on the iconic insulating Na-hP4 HPE show that increasing density causes a 3s→3pd electronic transition due to Pauli repulsion between the 1s2s and 3s states, and orthogonality of the 3pd states to the core. The large lobes of the resulting Na-pd hybrid orbitals point towards the center of an 11-membered penta-capped trigonal prism and overlap constructively, forming multicentered bonds, which are responsible for the emergence of the interstitial charge localization in Na-hP4. As a result, these multicentered bonds facilitate the increased density of this phase, which is key for its stabilization under pressure.

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

Equipartition and the Temperature of Maximum Density of TIP4P/2005 Water

Here, we simulate TIP4P/2005 water in the temperature range of 257 to 318 K with time-steps δ = 0.25, 0.50, 1.00, 2.00, and 4.00 fs. The density–temperature behavior obtained using 0.25 or 0.50 fs is in excellent agreement with each other but differs from those obtained using time steps that have been shown earlier to lead to a breakdown of equipartition. For δt = 0.25 or 0.50 fs, the temperature of maximum density (TMD) is 277.15 K and the density value is in close agreement with experiments. For δt = 1.00 fs, the TMD is 277.15 K, but the density value is shifted higher. For the other time steps considered here, the TMD is shifted to progressively lower values for longer time steps, a trend that holds for different thermostat/barostat combinations. Enhancing the water–water dispersion interaction, as has been recommended for simulating disordered proteins in TIP4P/2005, degrades the description of the liquid–vapor phase envelope. We present a simple physically transparent explanation that highlights the separation of the time scales between translational and rotational motion. We also develop a metric, χ, that we term the equipartition anomaly, to detect equipartition violations in simulations that include molecules that are treated as rigid objects. Calculating χ is shown to be straightforward and sensitive to equipartition violations. A key takeaway from this study is that using sufficiently short time steps (≤0.5 fs) to preserve equipartition is essential for obtaining meaningful liquid water properties and for producing reliable simulation data, as correct ensemble sampling is fundamental to ensure reproducibility across codes and simulation algorithms.

Asthagiri, Dilipkumar N. [Oak Ridge National Labor↗
DOE OSTI2023

Attosecond Time Delay Trends across the Isoelectronic Noble Gas Sequence

The analysis and measurement of Wigner time delays can provide detailed information about the electronic environment within and around atomic and molecular systems, with one the key differences being the lack of a long-range potential after a halogen ion undergoes photoionization. In this work, we use relativistic random-phase approximation to calculate the average Wigner delay from the highest occupied subshells of the atomic pairings (2p, 2s in Fluorine, Neon), (3p, 3s in Chlorine, Argon), (4p, 4s, 3d, in Bromine, Krypton), and (5p, 5s, 4d in Iodine, Xenon). The qualitative behaviors of the Wigner delays between the isoelectronic pairings were found to be similar in nature, with the only large differences occurring at photoelectron energies less than 20 eV and around Cooper minima. Interestingly, the relative shift in Wigner time delays between negatively charged halogens and noble gases decreases as atomic mass increases. All atomic pairings show large differences at low energies, with noble gas atoms showing large positive Wigner delays, while negatively charged halogen ions show negative delays. The implications for photoionization studies in halide-containing molecules is also discussed.

Grafstrom, Brock↗
DOE OSTI2025

Dataset for manuscript "Equipartition and the temperature of maximum density of TIP4P/2005 water"

We simulate TIP4P/2005 water in the temperature range of 257 K to 318 K with time-steps 0.25, 0.50, 1.00, 2.00, and 4.00 fs. The density-temperature behavior obtained using 0.25 or 0.50 fs are in excellent agreement with each other but differ from those obtained using time-steps that have been shown earlier to lead to a breakdown of equipartition. The temperature of maximum density (TMD) is 277.15 K with time-step 0.25 or 0.50 fs, but is shifted to progressively lower values for longer time-steps, a trend that holds for different thermostat/barostat combinations. Enhancing the water-water dispersion interaction, as has been recommended for simulating disordered proteins in TIP4P/2005, degrades the description of the liquid-vapor phase envelope. We present a simple physically transparent reasoning to highlight the separation of the time-scales between translational and rotational motion. We also develop a metric, Chi, that we term the equipartition anomaly, to detect equipartition violations in simulations that include molecules that are treated as rigid objects. Calculating Chi is shown to be straightforward and sensitive to equipartition violations. A key takeaway from this study is that using sufficiently short time-steps (less than or equal to 0.5 fs) to preserve equipartition is essential for obtaining meaningful liquid water properties and for producing reliable simulation data, as correct-ensemble sampling is fundamental to ensure reproducibility across codes and simulation alogrithms. The included dataset provides the raw data used in the preparation of the graphs noted in the manuscript.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Temperature measurement of Quark-Gluon plasma at different stages

In a Quark-Gluon Plasma (QGP), the fundamental building blocks of matter, quarks and gluons, are under extreme conditions of temperature and density. A QGP could exist in the early stages of the Universe, and in various objects and events in the cosmos. The thermodynamic and hydrodynamic properties of the QGP are described by Quantum Chromodynamics (QCD) and can be studied in heavy-ion collisions. Despite being a key thermodynamic parameter, the QGP temperature is still poorly known. Thermal lepton pairs (e + e − and μ + μ − ) are ideal penetrating probes of the true temperature of the emitting source, since their invariant-mass spectra suffer neither from strong final-state interactions nor from blue-shift effects due to rapid expansion. Here we measure the QGP temperature using thermal e+e− production at the Relativistic Heavy Ion Collider (RHIC). The average temperature from the low-mass region (in-medium ρ 0 vector-meson dominant) is (2.01 ± 0.23) × 10 12 K, consistent with the chemical freeze-out temperature from statistical models and the phase transition temperature from Lattice QCD. The average temperature from the intermediate mass region (above the ρ0 mass, QGP dominant) is significantly higher at (3.25 ± 0.60) × 10 12 K. This work provides essential experimental thermodynamic measurements to map out the QCD phase diagram and understand the properties of matter under extreme conditions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Asymmetric Fraunhofer spectra in a topological insulator-based Josephson junction

Josephson junctions with topological insulators as their weak link (S–TI–S junctions) are predicted to host Majorana fermions, which are key to creating qubits for topologically protected quantum computing. But the details of the S–TI–S current-phase relation and its interplay with magnetic fields are not well understood. We fabricate a Bi 2 Se 3 junction with NbTi leads and measure the Fraunhofer patterns of the junction with applied in-plane fields. We observe that asymmetric Fraunhofer patterns appear in the resistance maps of B z vs B x,y , with aperiodic node spacings. These asymmetric patterns appear even at zero parallel field and for temperatures up to 1 K. The anomalous features are compared to asymmetric Fraunhofer patterns expected for finite Cooper pair momentum shifts as well as geometric effects. We show that the geometric effects can dominate, independent of in-plane field magnitude. Furthermore, these results are important for differentiating geometrical phase shifts from those caused by Cooper pair momentum shifting, Majorana mode signatures, or other unconventional superconducting behavior.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗