Search NASA⌕ Search

SEARCH · Search NASA

Results for “Intrinsic computation”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Advection algorithms for quantum neutrino moment transport

Neutrino transport in compact objects is an inherently challenging multidimensional problem. Here, this difficulty is compounded if one includes flavor transformation—an intrinsically quantum phenomenon requiring one to follow the coherence between flavors and thus necessitating the introduction of complex numbers. To reduce the computational burden, simulations of compact objects that include neutrino transport often make use of momentum-angle-integrated moments (the lowest order ones being commonly referred to as the energy density and flux) and these quantities can be generalized to include neutrino flavor, i.e., they become quantum moments. Numerous finite-volume approaches to solving the moment evolution equations for classical neutrino transport have been developed based on solving a Riemann problem at cell interfaces. In this paper we describe our generalization of a Riemann solver for quantum moments, specifically decomposing complex numbers in terms of a (signed) magnitude and phase instead of real and imaginary parts. We then test our new algorithm in numerous cases showing a neutrino fast flavor instability, varying from toy models with analytic solutions to snapshots from neutron star merger simulations. Compared to previous algorithms for neutrino transport with flavor mixing, we find uniformly smaller growth rates of the flavor transformation along with concomitantly larger length-scales, and that the results are a better match with the growth rates seen from multiangle codes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Far-forward neutrinos at the Large Hadron Collider

We present a new calculation of the energy distribution of high-energy neutrinos from the decay of charm and bottom hadrons produced at the Large Hadron Collider (LHC). In the kinematical region of very forward rapidities, heavy-flavor production and decay is a source of tau neutrinos that leads to thousands of charged-current tau neutrino events in a 1 m long, 1 m radius lead neutrino detector at a distance of 480 m from the interaction region. In our computation, next-to-leading order QCD radiative corrections are accounted for in the production cross-sections. Non-perturbative intrinsic-k T effects are approximated by a simple phenomenological model introducing a Gaussian k T -smearing of the parton distribution functions, which might also mimic perturbative effects due to multiple initial-state soft-gluon emissions. The transition from partonic to hadronic states is described by phenomenological fragmentation functions. To study the effect of various input parameters, theoretical predictions for $D^{±}_s$ production are compared with LHCb data on double-differential cross-sections in transverse momentum and rapidity. The uncertainties related to the choice of the input parameter values, ultimately affecting the predictions of the tau neutrino event distributions, are discussed. We consider a 3+1 neutrino mixing scenario to illustrate the potential for a neutrino experiment to constrain the 3+1 parameter space using tau neutrinos and antineutrinos. We find large theoretical uncertainties in the predictions of the neutrino fluxes in the far-forward region. Untangling the effects of tau neutrino oscillations into sterile neutrinos and distinguishing a 3+1 scenario from the standard scenario with three active neutrino flavours, will be challenging due to the large theoretical uncertainties from QCD.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Niobium in the clean limit: An intrinsic type-I superconductor

Niobium is one of the most researched superconductors, both theoretically and experimentally. It is enormously significant in all branches of superconducting applications, from powerful magnets to quantum computing. It is therefore imperative to understand its fundamental properties in great detail. In this work, we use the results of recent microscopic calculations of anisotropic electronic, phonon, and superconducting properties, and apply a thermodynamic criterion for the type of superconductivity, more accurate and straightforward than a conventional Ginzburg-Landau parameter κ-based delineation, to show that pure niobium is a type-I superconductor in the clean limit. However, disorder (impurities, defects, strain, stress) pushes it to become a type-II superconductor.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Response Functions of Correlated Systems within Green's Function Theory

Why use Green's functions as the fundamental variable? Wave-function (Psi) methods are king for high-fidelity and Density-functional (Rho) methods are very efficient (Kohn-Sham). Goldilocks principle: Green's function (G) methods straddle the Rho and Psi methods, intermediate in both accuracy and efficiency. Also, when interest lies in excitations & 2-particle properties: G-methods are natural - intrinsic to the theory.

DMFT↗
DOE OSTI2025

Bigpicc: a graph-based approach to identifying carcinogenic gene combinations from mutation data

Abstract Genome data from cancer patients represents relationships between the presence of a gene mutation and cancer occurrence in a patient. Different types of cancer in human are thought to be caused by combinations of two to nine gene mutations. Identifying these combinations through traditional exhaustive search requires the amount of computation that scales exponentially with the combination size and in most cases is intractable even for cutting-edge supercomputers. We propose a parameter-free heuristic approach that leverages the intrinsic topology of gene-patient mutations to identify carcinogenic combinations. The biological relevance of the identified combinations is measured by using them to predict the presence of tumor in previously unseen samples. The resulting classifiers for 16 cancer types perform on par with exhaustive search results, and score the average of 80.1% sensitivity and 91.6% specificity for the best choice of hit range per cancer type. Our approach is able to find higher-hit carcinogenic combinations targeting which would take years of computations using exhaustive search.

Biochemistry & Molecular Biology↗
DOE OSTI2020

Single Cobalt Sites Dispersed in Hierarchically Porous Nanofiber Networks for Durable and High-Power PGM-Free Cathodes in Fuel Cells

Increasing catalytic activity and durability of atomically dispersed metal–nitrogen–carbon (M–N–C) catalysts for the oxygen reduction reaction (ORR) cathode in proton-exchange-membrane fuel cells remains a grand challenge. In this study, a high-power and durable Co–N–C nanofiber catalyst synthesized through electrospinning cobalt-doped zeolitic imidazolate frameworks into selected polyacrylonitrile and poly(vinylpyrrolidone) polymers is reported. The distinct porous fibrous morphology and hierarchical structures play a vital role in boosting electrode performance by exposing more accessible active sites, providing facile electron conductivity, and facilitating the mass transport of reactant. The enhanced intrinsic activity is attributed to the extra graphitic N dopants surrounding the CoN 4 moieties. The highly graphitized carbon matrix in the catalyst is beneficial for enhancing the carbon corrosion resistance, thereby promoting catalyst stability. The unique nanoscale X-ray computed tomography verifies the well-distributed ionomer coverage throughout the fibrous carbon network in the catalyst. The membrane electrode assembly achieves a power density of 0.40 W cm –2 in a practical H 2 /air cell (1.0 bar) and demonstrates significantly enhanced durability under accelerated stability tests. The combination of the intrinsic activity and stability of single Co sites, along with unique catalyst architecture, provide new insight into designing efficient PGM-free electrodes with improved performance and durability.

25 ENERGY STORAGE↗
DOE OSTI2025

Theoretical and experimental quantification of Suzuki segregation enthalpy and strengthening mechanisms in a binary alloy

Solute segregation to planar defects in metallic alloys has been shown to drastically alter mechanical properties. While various works using first-principles and thermodynamic calculations have studied the fundamental driving forces for solute segregation via the Suzuki criterion, planar defect energy, or a comparison of energies of the HCP-like phase and FCC matrix, a quantitative experimental and computational comparison of equilibrium composition and segregation enthalpies has not yet been reported. In this work, we predict the equilibrium composition and segregation enthalpy to intrinsic stacking faults in a Ni-60Co (at.%) alloy and compare the results to two independent experimental methods. We observed that Co segregates to the innermost two planes of the intrinsic stacking fault, and we found that the experimental segregation enrichment, measured from transmission electron microscopy energy dispersive X-ray spectroscopy, of the faults is 6.8 at.% Co, which is 2.2 at.% less than the predicted value at the same temperature. We also find that the segregation enthalpy measured from the composition profile is −21.1 ± 6.4 meV/atom and separately from differential scanning calorimetry segregation enthalpy is −33.2 meV/atom, whereas the predicted enthalpy is −31 ± 1 meV/atom. Based on these results, we determine that segregation occurs very rapidly, within 8 min at temperatures as low as 36% of the homologous solidus temperature. Furthermore, this analysis provides an overview of the possible dislocation mechanisms responsible for strengthening effects due to solute segregation, and concludes that changes in room temperature hardness from local phase transformation is likely tied to post-segregation room temperature equilibrium partial separation distance.

Ab initio calculation↗
DOE OSTI2022

On-shell correlators and color-kinematics duality in curved symmetric spacetimes

We define a perturbatively calculable quantity — the on-shell correlator — which furnishes a unified description of particle dynamics in curved spacetime. Specializing to the case of flat and anti-de Sitter space, on-shell correlators coincide precisely with on-shell scattering amplitudes and boundary correlators, respectively. Remarkably, we find that symmetric manifolds admit a generalization of on-shell kinematics in which the corresponding momenta are literally the isometry generators of the spacetime acting on the external kinematic data. These isometric momenta are intrinsically non-commutative but exhibit on-shell conditions that are identical to those of flat space, thus providing a common language for computing and representing on-shell correlators which is agnostic about the underlying geometry. Afterwards, we compute tree-level on-shell correlators for biadjoint scalar (BAS) theory and the nonlinear sigma model (NLSM) and learn that color-kinematics duality is manifested at the level of fields under a mapping of the color algebra to the algebra of gauged isometries on the spacetime manifold. Last but not least, we present a field theoretic derivation of the fundamental BCJ relations for on-shell correlators following from the existence of certain conserved currents in BAS theory and the NLSM.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

SAXS of murine amelogenin identifies a persistent dimeric species from pH 5.0 to 8.0

Amelogenin is an intrinsically disordered protein essential to tooth enamel formation in mammals. Here, using advanced small angle X-ray scattering (SAXS) capabilities at synchrotrons and computational models, we revisited measuring the quaternary structure of murine amelogenin as a function of pH and phosphorylation at serine-16. The SAXS data shows that at the pH extremes, amelogenin exists as an extended monomer at pH 3.0 (R g = 38.4 Å) and nanospheres at pH 8.0 (R g = 84.0 Å), consistent with multiple previous observations. At pH 5.0 and above there was no evidence for a significant population of monomeric species. Instead, at pH 5.0 ~ 80% of the population is a heterogenous dimeric species that increases to ~ 100% at pH 5.5. The dimer population was observed at all pH > 5 conditions in dynamic equilibrium with a species in the pentamer range at pH < 6.5 and nanospheres at pH 8.0. At pH 8.0 ~ 40% of the amelogenin still remained in the dimeric state. In general, serine-16 phosphorylation of amelogenin appears to modestly stabilize the population of the dimeric species.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Analysis of the bonding in tetrahedrane and phosphorus-substituted tetrahedranes

The bonding structures of tetrahedrane, phosphatetrahedrane, diphosphatetrahedrane and triphosphatetrahedrane are studied by employing an intrinsic quasi-atomic orbital analysis. Ethane, cyclopropane and tetrahedral P 4 are employed as reference systems. The orbital analysis is paired with the computation of strain energies via isodesmic reactions. The results show that the increase in geometric strain upon transition from ethane to cyclopropane to tetrahedrane weakens the CC bonds, despite leading to shorter C–C interatomic distances. With the increase in strain, the orbitals centered on C and involved in the bonding of the cage structure are observed to have elevated p-character, and the orbital structure of C deviates from sp 3 hybridization. The systematic substitution of CH groups by P atoms in the cage structure of tetrahedrane leads to stronger CC bonds, larger angles in the cage structures of the resulting phosphatetrahedranes, lower p-character in the orbitals involved in the bonding of the cages, and lower strain energies. It is found that P is more amenable to strained molecular arrangements than is C, and that the propensity of a given atom to hybridize s and p functions, or the lack thereof, has implications in the stability of molecules with strained geometries. The combination of the calculations presented here with the existing literature provides insight into the apparent propensity of tetrahedrane and P 4 to ‘break’ their tetrahedral structures. Trends in the bonding interactions, such as bond strengths, s- and p-orbital characters and charge transfer are identified and related to the strain energy observed in each of the analyzed systems.

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

Scalable Deep-Learning-Accelerated Topology Optimization for Additively Manufactured Materials

Topology optimization (TO) is a popular and powerful computational approach for designing novel structures, materials, and devices. Two computational challenges have limited the applicability of TO to a variety of industrial applications. First, a TO problem often involves a large number of design variables to guarantee sufficient expressive power. Second, many TO problems require a large number of expensive physical model simulations, and those simulations cannot be parallelized. To address these issues, we propose a general scalable deep-learning (DL) based TO framework, referred to as SDL-TO, which utilizes parallel schemes in high performance computing (HPC) to accelerate the TO process for designing additively manufactured (AM) materials. Unlike the existing studies of DL for TO, our framework accelerates TO by learning the iterative history data and simultaneously training on the mapping between the given design and its gradient. The surrogate gradient is learned by utilizing parallel computing on multiple CPUs incorporated with a distributed DL training on multiple GPUs. The learned TO gradient enables a fast online update scheme instead of an expensive update based on the physical simulator or solver. Using a local sampling strategy, we achieve to reduce the intrinsic high dimensionality of the design space and improve the training accuracy and the scalability of the SDL-TO framework. The method is demonstrated by benchmark examples and AM materials design for heat conduction. The proposed SDL-TO framework shows competitive performance compared to the baseline methods but significantly reduces the computational cost by a speed up of around 8.6x over the standard TO implementation.

Bi, Sirui↗
DOE OSTI2025

Computational toolkit for predicting thickness of 2D materials using machine learning and autogenerated dataset by large language model

The thickness of 2D materials not only plays a crucial role in determining the performance of nanoelectronic and optoelectronic devices but also introduces complexities in predicting volume-dependent properties, such as energy storage capacity, due to the intrinsic vacuum within these materials. Although a plethora of experimental techniques, including but not limited to optical contrast, Raman spectroscopy, nonlinear optical spectroscopy, near-field optical imaging, and hyperspectral imaging, facilitate the measurement of 2D material thickness, comprehensive data for many materials remain elusive. Over the past decade, the exponential proliferation of 2D materials and their heterostructures has outstripped the capabilities of conventional experimental and computational approaches. In this evolving landscape, machine learning (ML) has emerged as an indispensable tool, offering a scalable approach to augment these traditional methodologies. Addressing the critical gap, we introduce THICK2D—Thickness Hierarchy Inference and Calculation Kit for 2D Materials. This Python-based computational framework harnesses an autogenerated thickness database, developed using large language models, and advanced ML algorithms to facilitate the rapid and scalable estimation of material thickness, relying solely on crystallographic data. To demonstrate the utility and robustness of THICK2D, we successfully used the toolkit to predict the thickness of more than 8000 2D-based materials, sourced from two extensive 2D materials databases. THICK2D is disseminated as an open-source utility, accessible on GitHub at https://github.com/gmp007/THICK2D, and archived on Zenodo at https://10.5281/zenodo.11216648.

Ekuma, Chinedu E. (ORCID:0000000258527556)↗
DOE OSTI2021

Clock transitions guard against spin decoherence in singlet fission

Short coherence times present a primary obstacle in quantum computing and sensing applications. In atomic systems, clock transitions (CTs), formed from avoided crossings in an applied Zeeman field, can substantially increase coherence times. We show how CTs can dampen intrinsic and extrinsic sources of quantum noise in molecules. Conical intersections between two periodic potentials form CTs in electron paramagnetic resonance experiments of the spin-polarized singlet fission photoproduct. We report on a pair of CTs for a two-chromophore molecule in terms of the Zeeman field strength, molecular orientation relative to the field, and molecular geometry.

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

CryoDRGN-AI: neural ab initio reconstruction of challenging cryo-EM and cryo-ET datasets

Proteins and other biomolecules form dynamic macromolecular machines that are tightly orchestrated to move, bind, and perform chemistry. Cryo-electron microscopy (cryo-EM) and cryo-electron tomography (cryo-ET) can access the intrinsic heterogeneity of these complexes and are therefore key tools for understanding their function. However, 3D reconstruction of the collected imaging data presents a challenging computational problem, especially without any starting information, a setting termed ab initio reconstruction. Here, in this study, we introduce cryoDRGN-AI, a method leveraging an expressive neural representation and combining an exhaustive search strategy with gradient-based optimization to process challenging heterogeneous datasets. Using cryoDRGN-AI, we reveal new conformational states in large datasets, reconstruct previously unresolved motions from unfiltered datasets, and demonstrate ab initio reconstruction of biomolecular complexes from in situ data. With this expressive and scalable model for structure determination, we hope to unlock the full potential of cryo-EM and cryo-ET as a high-throughput tool for structural biology and discovery.

Levy, Axel [Stanford Univ., CA (United States); SL↗
DOE OSTI2025

A direct detection method of galaxy intrinsic ellipticity-gravitational shear correlation in non-linear regimes using self-calibration

Intrinsic alignment (IA) of galaxies is a challenging source of contamination in the Cosmic shear (GG) signals. The galaxy intrinsic ellipticity-gravitational shear (IG) correlation is generally the most dominant component of such contamination for cross-correlating redshift bins. One of the most effective techniques to mitigate such contamination is the self-calibration (SC) method which extracts the IG correlation and allows for its removal from the GG signal. In a photometric survey, the SC method first extracts the galaxy number density-galaxy intrinsic ellipticity (gI) correlation from the observed galaxy-galaxy lensing correlation using the redshift dependence of lens-source pairs. The IG correlation is computed through a scaling relation using the gI correlation and other lensing observables. The applicability of the SC method has so far been focused on the linear IA scales and the linear galaxy bias. We extend the SC method beyond the linear regime by modifying its scaling relation which can account for the non-linear galaxy bias model and various IA models. In this study, we provide a framework to detect the IG correlation for the redshift bins for source galaxies for the proposed year 1 survey of the Rubin Legacy Survey of Space and Time (LSST Y1). We tested the method for the tidal alignment and tidal torquing (TATT) model of IA and we found that the scaling relation is accurate within 10% and 20% for cross-correlating and auto-correlating redshift bins, respectively. Hence the suppression of IG contamination in observed GG correlation can be accomplished with a factor of 10 and 5, for cross-correlating and auto-correlating redshift bins, respectively. We tested the method's robustness and found that the suppression of IG contamination by a factor of 5 is still achievable for all combinations of cross-correlating bins even with the inclusion of a moderate amount of uncertainties on IA and bias parameters, respectively. We also make available, a branch of the code FAST-PT to provide gI correlations up to 1-loop order term used by the new SC method.

gravitational lensing↗
DOE OSTI2020

Higher-order phonon scattering: advancing the quantum theory of phonon linewidth, thermal conductivity and thermal radiative properties

Phonon scattering plays a central role in the quantum theory of phonon linewidth, which in turn governs important properties including infrared spectra, Raman spectra, lattice thermal conductivity, thermal radiative properties, and also signifi- cantly affects other important processes such as hot electron relaxation. Since Maradudin and Fein’s classic work in 1962, three-phonon scattering had been considered as the dominant intrinsic phonon scattering mechanism and has seen tremendous advances. However, the role of the higher-order four-phonon scattering had been persistently unclear and so was ignored. The tremendous complexity of the formalism and computational challenges stood in the way, prohibiting the direct and quantitative treatment of four-phonon scattering. In 2016, a rigorous four-phonon scattering formalism was developed, and the prediction was realized using empirical potentials. In 2017, the method was extended using first-principles calculated force constants, and the thermal conductivities of boron arsenides (BAs), Si and diamond were predicted. The predictions for BAs were later confirmed by several independent experiments. Four-phonon scattering has since been investigated in a range of materials and established as an important intrinsic scattering mechanism for thermal transport and radiative properties. Specifically, four-phonon scattering is important when the fourth-order scattering potential or phase space becomes relatively large. The former scenario includes: (i) nearly all materials when the temperature is high; (ii) strongly anharmonic (low thermal conductivity) materials, including most rocksalt compounds, halides, hydrides, chalcogenides and oxides. The latter scenario includes: (iii) materials with large acoustic–optical phonon band gaps, such as XY compounds with a large atomic mass ratio between X and Y; (iv) two- dimensional materials with reflection symmetry, such as single-layer graphene, single-layer boron nitride and carbon nanotubes; and (v) phonons with a large density of states, such as optical phonons, which are important for Raman, infrared and thermal radiative properties. Four-phonon scattering is expected to gain broad interest in various technologically important materials for thermoelectrics, thermal barrier coatings, thermal energy storage, phase change, nuclear power, ultra-high temperature ceramics, infrared spectra, Raman spectra, radiative transport, hot electron relaxation and radiative cooling. Four-phonon scattering has been, and will continue to be, established as an important intrinsic phonon scattering mechanism beyond three-phonon scattering. The prediction of four-phonon scattering will transition from a breakthrough to a new routine in the next decade.

Feng, Tianli↗
DOE OSTI2026

Selective Chemical Looping Combustion of Terminal Alkynes in Mixtures with Alkenes

The selective combustion of terminal alkynes in mixtures with alkenes is demonstrated during anaerobic reduction half-cycles on bulk bismuth oxide (Bi 2 O 3 ) as an approach to remove alkynes, which act as inhibitors in olefin polymerization. Bi 2 O 3 combusts phenylacetylene in styrene, 3-methylphenylacetylene in 3-methylstyrene, propyne in propylene, 1-hexyne in 1-hexene, and 1-octyne in 1-octene, with alkyne combustion selectivities exceeding 96%. Near unity reaction orders for hydrocarbon consumption during reduction half-cycles are consistent with combustion pathways initiated by rate-determining initial C–H activation, which drive selective alkyne combustion through intrinsic differences in the first-order rate constants for alkyne and alkene combustion rather than preferential adsorption of alkynes on Bi 2 O 3 surfaces. Computational assessments of initial C–H activation pathways for alkynes and alkenes on (010) α-Bi 2 O 3 surfaces using density functional theory illustrate that heterolytic transition states which form proton-carbanion pairs on Bi–O sites kinetically favor the activation of alkynes rather than alkenes due to differences in C–H bond acidity, and the barrier for heterolytic C–H activation is dictated in part by the sum of the molecular deprotonation energy and the energy to bind an R – carbanion to a Bi site in its transition-state geometry. Finally, these heterolytic reactivity channels during selective chemical looping combustion present novel routes for purifying olefin gas streams containing alkyne impurities.

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

DS-TIDE: Harnessing Dynamical Systems for Efficient Time-Independent Differential Equation Solving

Time-Independent Differential Equations (TIDEs) are central to modeling equilibrium behavior across a wide range of scientific and engineering domains, from electrostatics to porous media flow. Conventional numerical solvers offer reliable solutions but incur significant computational costs due to fine-grained discretization and iterative procedures. Machine learning-based approaches address this by replacing iterative solving processes with one-time inference; however, their sophisticated models require extensive training resources that often exceed those of traditional solvers. Consequently, designing a TIDE solver that achieves high accuracy, broad applicability, and exceptional computational efficiency remains a fundamental challenge. In this paper, we propose DS-TIDE, a novel hardware solver that is inspired by, and subsequently leverages, the intrinsic connection between Dynamical Systems (DS) and Differential Equations (DEs) to efficiently and accurately solve TIDEs. DS-TIDE employs a CMOS-compatible DS-based processor, whose physical states evolve under carefully designed DE-driven dynamics and naturally converge to equilibrium -- the solution of the target TIDE -- within ~1µs on a ~1-watt DS-TIDE processor. To enhance expressivity, DS-TIDE incorporates Heterogeneous Dynamics with Temporal Layering (HDTL), which solves TIDEs through a three-stage DS evolution -- conditioning, solving, and decoding -- each governed by specialized dynamics. The entire evolution process is analogous to an infinitely deep neural network temporally unrolled, offering the system the capability of representing complex equations. Furthermore, DS-TIDE is equipped with an on-device DS-DE Auto-Alignment mechanism that dynamically adapts intrinsic hardware dynamics within milliseconds, effectively aligning the system’s dynamics to diverse target DEs. Experimental results across TIDEs from a wide range of scientific and engineering domains demonstrate that DS-TIDE achieves ~10^3× speedup, ~10^5× energy savings, and competitive or superior accuracy compared to state-of-the-art numerical and ML-based solvers.

Liu, Chuan↗