Search NASA⌕ Search

SEARCH · Search NASA

Results for “hybrid regimes”

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 37 records · Page 2

Electron kinetics in a positive column of AC discharges in a dynamic regime

Abstract We have performed hybrid kinetic-fluid simulations of a positive column in alternating current (AC) argon discharges over a range of driving frequencies f and gas pressure p for the conditions when the spatial nonlocality of the electron energy distribution function (EEDF) is substantial. Our simulations confirmed that the most efficient conditions of plasma maintenance are observed in the dynamic regime when time modulations of mean electron energy (temperature) are substantial. The minimal values of the root mean square electric field and the electron temperature have been observed at f/p values of about 3 kHz Torr −1 in a tube of radius R = 1 cm. The ionization rate and plasma density reached maximal values under these conditions. The numerical solution of a kinetic equation allowed accounting for the kinetic effects associated with spatial and temporal nonlocality of the EEDF. Using the kinetic energy of electrons as an independent variable, we solved an anisotropic tensor diffusion equation in phase space. We clarified the role of different flux components during electron diffusion in phase space over surfaces of constant total energy. We have shown that the kinetic theory uncovers a more exciting and rich physics than the classical ambipolar diffusion (Schottky) model. Non-monotonic radial distributions of excitation rates, metastable densities, and plasma density have been observed in our simulations at pR > 6 Torr cm. The predicted off-axis plasma density peak in the dynamic regime has never been observed in experiments so far. We hope our results stimulate further experimental studies of the AC positive column. The kinetic analysis could help uncover new physics even for such a well-known plasma object as a positive column in noble gases.

Physics↗

Vortex Matter in Confined Superconductors and Mesoscopic Hybrid Heterostructures. Final report

This program studied the fundamental features that underlie the behavior of hybrid ferromagnet/superconductor heterostructures both in the regime of magnetic coupling and proximity effect. In ferromagnet/superconductor hybrid structures, the source of superconducting correlations and ferromagnetism are spatially separated. Therefore, the interplay between ferromagnetism and superconductivity results in physical effects that are not easily observed in bulk superconductors and allows to tune the local superconducting and vortex properties, realizing, thus, a new class of materials. In the case of proximity effect the ferromagnet and the superconductor are in contact and leakage of superconducting pairs is allowed into the ferromagnetic layer. Typically, the decay length for the superconducting pair correlations in the ferromagnet is of the order of one nanometer due to the typically large ferromagnetic exchange field. However, a suitable magnetic structure in the ferromagnetic material can be used to align the spins of the superconducting pairs allowing much longer penetration of superconductivity in the magnetic material. These superconducting pairs with aligned spins can carry spin and provide the basis for superconducting spintronics devices. In the case of magnetic coupling the effect of a spatial inhomogeneous stray field produced by the ferromagnet substantially modifies the nucleation of the order parameter in the superconducting layer and can lead to a significant change in the local superconducting properties. Therefore, magnetically coupled planar ferromagnet-superconductor hybrid structures offer new avenues for manipulation of the superconductivity at the nanoscale and convenient means to control vortex dynamics.

36 MATERIALS SCIENCE↗

A Holistic DSMC Transport Database for Re-Entry and Ablation Modeling

Hybrid simulation frameworks combining Computational Fluid Dynamics (CFD) and Direct Simulation Monte Carlo (DSMC) are frequently employed to efficiently perform high-fidelity solutions of environments containing combined continuum/rarified flow. The use of DSMC, a stochastic, particle-based method, is necessary for high-Knudsen flow where continuum-based assumptions governing CFD break down. However, the DSMC methodology is generally very computationally inefficient to model the continuum regime. In a CFD/DSMC hybrid approach, obtaining an accurate, high-fidelity solution hinges on the consistent treatment of transport properties and the used thermo-chemical models employed within the two solvers. In principle, in regions where CFD and DSMC are both employed, the same gas mixture under the same conditions should have the same properties, regardless of simulation type. Observed differences should be due to non-equilibrium processes, rather than differences in physical models. While the transport models governing CFD and DSMC simulations are starkly different, they can effectively be linked via their use of reduced Chapman-Enskog collision integrals. In CFD, these integrals are typically stored as fitted polynomial expressions and used to directly compute gas transport properties via mixing rules or the full Chapman-Enskog formulation. In DSMC, they can be used to derive the collision parameters needed for the phenomenological collision cross-section models that govern particle interactions, via a Nelder-Mead optimization scheme. The goal of this work is to provide a unified DSMC transport database encompassing the vast majority of known gas species encountered during atmospheric entry, on Earth or any other Solar body. This goal is largely possible due to recently performed ab-initio quantum chemistry calculations. Combined with other high-fidelity literature sources, the planned database will consist of collision integral data for over 200 neutral and ionized species and over 17000 binary collisions. From these collision integrals, Nelder-Mead optimization is used to compute Variable Soft Sphere (VSS) collision model parameters for DSMC, fitted from 300 K to 20000 K. Initial comparisons of transport properties of relevant equilibrium gas mixtures show great agreement between CFD and DSMC-derived results. The completed database will be able to be readily applied to model binary collisions of any gas mixture containing the included species over the specified temperature range, making it a valuable tool for future planetary probe modeling efforts. An example is shown below. Equilibrium mixture transport properties for a 19-species Titan atmospheric model [4] are computed using both fitted VSS parameters and the original CFD collision integral values. Deviations in computed properties between the two approaches is less than 5% for the entire temperature range.

M R Gosma↗

Hybrid Quantum Networks with Discrete Polarizations and Continuous Quadrature Variables

Final Scientific/Technical Report for DOE Project DE-SC0022069. Quantum communication and computation systems have evolved along two largely independent paradigms: discrete-variable (DV) systems that encode information in qubits such as photon polarizations or photon-number states, and continuous-variable (CV) systems that encode information in optical field quadratures. Each approach offers distinct advantages—DVs provide low error rates and compatibility with single-photon platforms, while CVs support deterministic operations and efficient quantum state manipulation. A fundamental challenge for building a scalable Quantum Internet lies in interfacing these two regimes into a unified hybrid architecture that can coherently distribute and process quantum information across heterogeneous quantum nodes. This project aims to develop and demonstrate a hybrid optical quantum network that seamlessly integrates DV and CV systems. Specifically, we investigate a new class of hybrid entanglement between the discrete polarizations of single photons and the continuous quadrature variables of optical cat states, overcoming incompatibilities in existing DV–CV demonstrations. Using this new entanglement resource, the team investigates a multi-node hybrid quantum local area network (Q-LAN) testbed capable of DV–CV entanglement generation, swapping, and distribution across fiber links.

74 ATOMIC AND MOLECULAR PHYSICS↗

Recommended DSMC Collision Model Parameters for Planetary Entry

Hybrid simulation frameworks combining Computational Fluid Dynamics (CFD) and Direct Simulation Monte Carlo (DSMC) are frequently employed to efficiently perform high-fidelity simulations of environments consisting of both continuum and rarified flow. DSMC is a stochastic, particle-based method which solves the fundamental Boltzmann equation and is therefore necessary for high-Knudsen flow where continuum-based assumptions governing CFD break down. However, the DSMC methodology is generally computationally inefficient to model the continuum regime. In a CFD/DSMC hybrid approach, obtaining an accurate, high-fidelity solution hinges on the consistent treatment of transport properties and the thermo-chemical models employed within the two solvers. In principle, in regions where CFD and DSMC are both employed, the same gas mixture under the same conditions should have the same properties, regardless of simulation type. Observed differences should be due to non-equilibrium processes, rather than differences in physical models. The goal of this work is to provide a comprehensive DSMC transport database encompassing the vast majority of known gas species encountered during Earth or other planetary atmospheric entry. This goal is largely possible due to recently performed ab initio quantum chemistry calculations. Combined with other high-fidelity data, the planned database will consist of collision integral data for over 200 neutral and ionized species and over 20000 binary collisions. From these collision integrals, Nelder-Mead optimization is used to compute collision-specific Variable Soft Sphere (VSS) collision model parameters, fitted from 300 K to 20000 K. Initial comparisons of transport properties of relevant equilibrium gas mixtures show great agreement between CFD and DSMC-derived results. The completed database can be readily applied to model binary collisions of any gas mixture containing the included species over the specified temperature range, making it a valuable tool for future planetary probe modeling efforts. An example is shown below in Fig. 1. Equilibrium mixture transport properties for a 35-species mixture composed originally of 10% air and 90% pyrolysis species of a carbon-phenolic ablator material [4] are computed using both fitted VSS Parameters and the original CFD collision integral values. Deviations in computed properties between the two approaches are less than 5% for the entire temperature range.

M. R. Gosma↗

Transition to whistler mediated magnetic reconnection

The transition in the magnetic reconnection rate from the resistive magnetohydrodynamic (MHD) regime where the Alfen wave controls reconnection to a regime in which the ions become unmagnetized and the whistler wave mediates reconnection is explored with 2-D hybrid simulations. In the whistler regime the electrons carry the currents while the ions provide a neutralizing background. A simple physical picture is presented illustrating the role of the whistler mediated reconnection is calculated analytically. The development of an out-of-plane component of the magnetic field is an observable signature of whistler driven reconnection.

Mandt, M. E.↗

Valley-selective optical Stark effect of exciton-polaritons in a monolayer semiconductor

Selective breaking of degenerate energy levels is a well-known tool for coherent manipulation of spin states. Though most simply achieved with magnetic fields, polarization-sensitive optical methods provide high-speed alternatives. Exploiting the optical selection rules of transition metal dichalcogenide monolayers, the optical Stark effect allows for ultrafast manipulation of valley-coherent excitons. Compared to excitons in these materials, microcavity exciton-polaritons offer a promising alternative for valley manipulation, with longer lifetimes, enhanced valley coherence, and operation across wider temperature ranges. Here, we show valley-selective control of polariton energies in WS 2 using the optical Stark effect, extending coherent valley manipulation to the hybrid light-matter regime. Ultrafast pump-probe measurements reveal polariton spectra with strong polarization contrast originating from valley-selective energy shifts. This demonstration of valley degeneracy breaking at picosecond timescales establishes a method for coherent control of valley phenomena in exciton-polaritons.

36 MATERIALS SCIENCE↗

Investigation of plasma instabilities in the polar cusp

During the last six months, considerable progress was made in studying the excitation of electromagnetic waves in the whistler frequency range by an anisotropic or gyrating electron beam. A paper entitled 'Electron Cyclotron Wave Generation by Relativistic Electrons' was published in the Journal of Geophysical Research. Another paper entitled 'Electron Beam Excitation of Upstream Waves in the Whistler Mode Frequency Range' was submitted for publication in Journal of Geophysical Research. This paper is in collaboration with Dr. C. W. Smith at Bartol Research Institute. In this paper, it was shown that an anisotropic electron beam (or gyrating electron beam) is capable of generating both left-hand and right-hand polarized electromagnetic waves in the whistler frequency range. Our earlier paper 'Electromagnetic Components of Auroral Hiss and Lower Hybrid Waves in the Polar Magnetosphere' was accepted for publication in the AGU Chapman Conference on Micro and Meso Scale Phenomena in Space Plasmas. Electromagnetic waves in the lower hybrid and whistler waves regime were identified and a mechanism of how these waves are generated was suggested.

Wong, H. K.↗

Kinetic physics of the mirror instability

Kinetic mechanisms for the growth and saturation of the mirror instability are described using one-dimensional hybrid simulations. Two parameter regimes are considered. In the first regime, a relatively small ion anisotropy excites a slowly growing instability that produces small-amplitude waves; most ions respond to the waves as an adiabatic fluid. In the second regime a large anisotropy excites a rapidly growing instability that generates large-amplitude waves; the response of many ions in this case is nonadiabatic. The difference in ion response is due to the relative importance of two ion populations, resonant and nonresonant. Resonant ions, those ions with low velocities parallel to the background magnetic field, contribute to the growth of the instability as a result of their gyrointeractions with the noncoplanar component of the waves electric field and respond to the mirror waves nonadiabatically. Nonreasonant ions, those with large parallel velocities, respond as an adiabatic fluid. In both regimes, ion anisotropy is reduced by means of the magnetic mirror force; in the second regime, the anisotropy reduction reduces the free energy available for wave growth and leads to the saturation of the mirror instability.

Mckean, M. E.↗

Analysis of Effectiveness of Phoenix Entry Reaction Control System

Interaction between the external flowfield and the reaction control system (RCS) thruster plumes of the Phoenix capsule during entry has been investigated. The analysis covered rarefied, transitional, hypersonic and supersonic flight regimes. Performance of pitch, yaw and roll control authority channels was evaluated, with specific emphasis on the yaw channel due to its low nominal yaw control authority. Because Phoenix had already been constructed and its RCS could not be modified before flight, an assessment of RCS efficacy along the trajectory was needed to determine possible issues and to make necessary software changes. Effectiveness of the system at various regimes was evaluated using a hybrid DSMC-CFD technique, based on DSMC Analysis Code (DAC) code and General Aerodynamic Simulation Program (GASP), the LAURA (Langley Aerothermal Upwind Relaxation Algorithm) code, and the FUN3D (Fully Unstructured 3D) code. Results of the analysis at hypersonic and supersonic conditions suggest a significant aero-RCS interference which reduced the efficacy of the thrusters and could likely produce control reversal. Very little aero-RCS interference was predicted in rarefied and transitional regimes. A recommendation was made to the project to widen controller system deadbands to minimize (if not eliminate) the use of RCS thrusters through hypersonic and supersonic flight regimes, where their performance would be uncertain.

Dyakonov, Artem A.↗

Edge localized mode characteristics and divertor heat flux during stationary and transient phase for CFETR hybrid scenario

The study of edge localized mode (ELM) behavior, stationary heat flux and transient heat flux during the ELM crash phase is performed for a China Fusion Engineering Test Reactor (CFETR) 1 GW hybrid mode operation scenario ($R$ = 7.2 m, ${B_{\text{T}}}$ = 6.5 T, ${I_{\text{p}}}$ = 13.78 MA). Modeling and simulation start with a scenario obtained by multi-code integrated modeling on the one modeling framework for integrated tasks framework. Linear stability and nonlinear simulations of ELM dynamics are carried out using the BOUT++ six-field reduced magnetohydrodynamic module, which show a much smaller ELM energy loss (Δ ELM ~ 0.13%) compared to that of a Type-I ELM. Parametric analysis of the weak linear growth rate and small ELM energy loss characteristics shed light on physics corresponding to a grassy ELM regime for CFETR 1 GW hybrid scenario. The transient heat flux on the divertor target during this small ELM phase is investigated using BOUT++. Here we found that upstream radial transport in the scrape-off-layer (SOL) induced by small ELMs is weak, which keeps it in the drift-dominated region. However, the heat flux width is still broadened to ${\lambda _{\text{q}}}$ = 4.64 mm by the increase of separatrix temperature during ELM nonlinear evolution. The impact of transient peak heat load and ELM energy fluence on tungsten melting and net erosion rate of divertor target is evaluated for the first-time using physics-based transport that connects the pedestal with the SOL. Energy fluence caused by a single ELM pulse is below the tungsten melting limit, while tungsten erosion would exceed the material requirements. We conclude that external mitigation methods, such as divertor detachment and advanced divertor geometry are likely needed for the steady state operation of CFETR.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Pt Particle Size Affects Both the Charge Separation and Water Reduction Efficiencies of CdS–Pt Nanorod Photocatalysts for Light Driven H 2 Generation

Decreasing the metal catalyst size into nanoclusters or even single atom is an emerging direction of developing more efficient and cost-effective photocatalytic systems. Because the catalyst particle size affects both the catalyst activity and light driven charge separation efficiency, their effects on the overall photocatalytic efficiency are still poorly understood. Herein, using a well-defined semiconductor–metal heterostructure with Pt nanoparticle catalysts selectively grown on the apexes of CdS nanorods (NRs), we study the effect of the Pt catalyst size on light driven H 2 generation quantum efficiency (QE H 2 ). With the increase of the Pt catalyst size from 0.7 ± 0.3 to 3.0 ± 0.8 nm, the QEH 2 of CdS–Pt increases from 0.5 ± 0.2% to 38.3 ± 5.1%, by nearly 2 orders of magnitude. Transient absorption spectroscopy measurement reveals that the electron transfer rate from the CdS NR to the Pt tip increases with the Pt diameter following a scaling law of d 5.6 , giving rise to the increase of electron transfer efficiency at larger Pt sizes. Additionally, the observed trend can be understood by a simplified kinetic model that assumes the overall efficiency is the product of the quantum efficiencies of charge separation (including hole transfer, electron transfer, and hole scavenging) and water reduction steps, and for CdS–Pt NRs, the quantum efficiencies of electron transfer and water reduction steps increase with the Pt sizes. Our findings suggest the importance of improving the quantum efficiencies of both charge separation and catalysis in designing efficient semiconductor–metal hybrid photocatalysts, especially in the regime of small metal particle sizes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multi-messenger dynamic imaging of laser-driven shocks in water using a plasma wakefield accelerator

Understanding dense matter hydrodynamics is critical for predicting plasma behavior in environments relevant to laser-driven inertial confinement fusion. Traditional diagnostic sources face limitations in brightness, spatiotemporal resolution, and in their ability to detect relevant electromagnetic fields. In this work, we present a dual-probe, multi-messenger laser wakefield accelerator platform combining ultrafast X-rays and relativistic electron beams at 1 Hz, to interrogate a free-flowing water target in vacuum, heated by an intense 200 ps laser pulse. This scheme enables high-repetition-rate tracking the evolution of the interaction using both particle types. Betatron X-rays reveal a cylindrically symmetric shock compression morphology assisted by low-density vapor, resembling foam-layer-assisted fusion targets. The synchronized electron beam detects time-evolving electromagnetic fields, uncovering charge separation and ion species differentiation during plasma expansion – phenomena not captured by photons or hydrodynamic simulations. We show that combining both probes provides complementary insights spanning kinetic to hydrodynamic regimes, highlighting the need for hybrid physics models to accurately predict fusion-relevant plasma behavior.

imaging techniques↗

Sub-nanosecond time-resolved radiation measurement using x-ray focusing crystal spectrometers

Here, in this paper, we describe a technique using a crystal spectrometer, a silicon-diode detector, and a filtered photoconductive detector to monitor photon energies in the L-shell (0.9–1 keV) and K-shell regimes for nickel and copper hybrid X-pinch x-ray sources. The detectors, system cabling, and an 8 GHz digital oscilloscope in combination enable time resolution better than 200 ps for photoconductive detectors and 700 ps for silicon-diode detectors of the K- and L-shell radiation signals, respectively. We substantially improve the relative timing of signals obtained using the oscilloscope by using an x-ray streak camera with a crystal spectrometer to monitor the L-shell line spectra and, separately, the K-shell line spectra relative to the continuum burst to better than 17 ps time resolution. This combination of instruments enabled and validated a new method by which plasma conditions in nickel and copper X-pinches can be assessed immediately before and after the ~30 ps continuum x-ray burst produced by 370 kA hybrid X-pinches. In general, the method described here can be applied to observe otherwise highly filter-absorbed radiation in the presence of a broad spectrum of higher energy radiation by combining x-ray crystals and detectors.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Analog Quantum Simulation of the Dynamics of Open Quantum Systems with Quantum Dots and Microelectronic Circuits

We introduce a general setup for the analog quantum simulation of the dynamics of open quantum systems based on semiconductor quantum dots electrically connected to a chain of quantum RLC electronic circuits. The dots are chosen to be in the regime of spin-charge hybridization to enhance their sensitivity to the RLC circuits while mitigating the detrimental effects of unwanted noise. In this context, we establish an experimentally realizable map between the hybrid system and a qubit coupled to thermal harmonic environments of arbitrary complexity that enables the analog quantum simulation of open quantum systems. We assess the utility of the simulator by numerically exact emulations that indicate that the experimental setup can faithfully mimic the intended target even in the presence of its natural inherent noise. We further provide a detailed analysis of the physical requirements on the quantum dots and the RLC circuits needed to experimentally realize this proposal that indicates that the simulator can be created with existing technology. The approach can exactly capture the effects of highly structured non-Markovian quantum environments typical of photosynthesis and chemical dynamics, and offer clear potential advantages over conventional and even quantum computation. The proposal opens a general path for effective quantum dynamics simulations based on semiconductor quantum dots.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Impact of Nanoparticle Size and Surface Chemistry on Peptoid Self-Assembly

Self-assembled organic nanomaterials can be generated by bottom-up assembly pathways where the structure is controlled by the organic sequence and altered using pH, temperature, and solvation. In contrast, self-assembled structures based on inorganic nanoparticles typically rely on physical packing and drying effects to achieve uniform superlattices. By combining these two chemistries to access inorganic–organic nanostructures, we aim to understand the key factors that govern the assembly pathway and structural outcomes in hybrid systems. In this work, we outline two assembly regimes between quantum dots (QDs) and reversibly binding peptoids. These regimes can be accessed by changing the solubility and size of the hybrid (peptoid-QD) monomer unit. The hybrid monomers are prepared via ligand exchange and assembled, and the resulting assemblies are studied using ex-situ transmission electron microscopy as a function of assembly time. In aqueous conditions, QDs were found to stabilize certain morphologies of peptoid intermediates and generate a final product consisting of multilayers of small peptoid sheets linked by QDs. The QDs were also seen to facilitate or inhibit assembly in organic solvents based on the relative hydrophobicity of the surface ligands, which ultimately dictated the solubility of the hybrid monomer unit. Increasing the size of the QDs led to large hybrid sheets with regions of highly ordered square-packed QDs. A second, smaller QD species can also be integrated to create binary hybrid lattices. Furthermore, these results create a set of design principles for controlling the structure and structural evolution of hybrid peptoid-QD assemblies and contribute to the predictive synthesis of complex hybrid matter.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Ghent Hybrid model in NuWro: a new neutrino single-pion production model in the GeV regime

Neutrino-induced single-pion production constitutes an essential interaction channel in modern neutrino oscillation experiments, with its products building up a significant fraction of the observable hadronic final states. Frameworks of oscillation analyses strongly rely on Monte Carlo neutrino event generators, which provide theoretical predictions of neutrino interactions on nuclear targets. Thus, it is crucial to integrate state-of-the-art single-pion production models with Monte Carlo simulations to prepare for the upcoming systematics-dominated landscape of neutrino measurements. In this work, we present the implementation of the Ghent Hybrid model for neutrino-induced single-pion production in the NuWro Monte Carlo event generator. The interaction dynamics includes coherently-added contributions from nucleon resonances and a non-resonant background, merged into the pythia branching predictions in the deep-inelastic regime, as instrumented by NuWro. This neutrino-nucleon interaction model is fully incorporated into the nuclear framework of the generator, allowing it to account for the influence of both initial- and final-state nuclear medium effects. We compare the predictions of this integrated implementation with recent pion production data from accelerator-based neutrino experiments. The results of the novel model show improved agreement of the generator predictions with the data and point to the significance of the refined treatment of the description of pion-production processes beyond the ∆ region.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Optimizing the critical temperature and superfluid density of a metal-superconductor bilayer

A promising path to realizing higher superconducting transition temperatures 𝑇 c is the strategic engineering of artificial heterostructures. For example, quantum materials could, in principle, be coupled with other materials to produce a more robust superconducting state. Here, in this work, we add numerical support to the hypothesis that a strongly interacting superconductor weakened by phase fluctuations can boost its 𝑇 c by hybridizing the system with a metal. Using determinant quantum Monte Carlo, we simulate a two-dimensional bilayer composed of an attractive Hubbard model and a metallic layer in two regimes of the interaction strength −|𝑈|. In the strongly interacting regime, we find that increasing the interlayer hybridization 𝑡 ⊥ results in a nonmonotonic enhancement of 𝑇 c , with an optimal value comparable to the maximum 𝑇 c observed in the single-layer attractive Hubbard model, confirming trends inferred from other approaches. In the intermediate coupling regime, when −|𝑈| is close to the value associated with the maximum 𝑇 c of the single-layer model, increasing 𝑡 ⊥ tends to decrease 𝑇 c , implying that the correlated layer was already optimally tuned. Importantly, we demonstrate that the mechanism behind these trends is related to enhancement in the superfluid stiffness, as was initially proposed by Kivelson [Phys. B: Condens. Matter 318, 61 (2002)].

Zhang, Yutan [Univ. of California, Davis, CA (Unit↗