Search NASA⌕ Search

SEARCH · Search NASA

Results for “Collision Rates”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Porting HEP Parameterized Calorimeter Simulation Code to GPUs

The High Energy Physics (HEP) experiments, such as those at the Large Hadron Collider (LHC), traditionally consume large amounts of CPU cycles for detector simulations and data analysis, but rarely use compute accelerators such as GPUs. As the LHC is upgraded to allow for higher luminosity, resulting in much higher data rates, purely relying on CPUs may not provide enough computing power to support the simulation and data analysis needs. As a proof of concept, we investigate the feasibility of porting a HEP parameterized calorimeter simulation code to GPUs. We have chosen to use FastCaloSim, the ATLAS fast parametrized calorimeter simulation. While FastCaloSim is sufficiently fast such that it does not impose a bottleneck in detector simulations overall, significant speed-ups in the processing of large samples can be achieved from GPU parallelization at both the particle (intra-event) and event levels; this is especially beneficial in conditions expected at the high-luminosity LHC, where extremely high per-event particle multiplicities will result from the many simultaneous proton-proton collisions. We report our experience with porting FastCaloSim to NVIDIA GPUs using CUDA. A preliminary Kokkos implementation of FastCaloSim for portability to other parallel architectures is also described.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Chamber Interior Viewing System (CIVS 2.0)

The NIF facility is used to conduct fusion and high-density scientific research. The experiments are conducted inside of a 10-meter diameter vacuum chamber called the “Target Chamber”. The target chamber is covered in ports to permit adding equipment to the exterior. Some of this equipment extends into the target chamber such as the target positioners and diagnostics. Other equipment sits on the surface and looks inward, such as the CIVS system. CIVS stands for Chamber Interior Viewing System. The CIVS system provides simple black and white digital images of the interior of the target chamber in real time. These images are used for positioner alignment, trouble shooting and observing for positioner collision. The views through these cameras need to be clear enough for the operators to safely conduct these tasks. One of the products of the experiments conducted in the target chamber are neutrons. These neutrons expand outward from the center of the target chamber with great energy. If the neutrons hit the sensor chip on a digital camera they will “kill” the sensor pixel where they contact, and the result will be a permanently bright white pixel on the image. Over time these “dead” pixels will build up and make the camera less and less usable providing a “staticky” image. The current solution of replacing the cameras when they become damaged is unsustainable due to both the cameras being out of production and the increased rate of camera failure based upon the increased number and energy of neutrons being generated. Current supplies of the correct camera type are dwindling with no source of replacement and the projected rate of replacement is unsustainable. Without a redesign, the function provided by the current CIVS system will be lost.

42 ENGINEERING↗
DOE OSTI2020

Deep learning surrogate model for kinetic Landau-fluid closure with collision

In this work, the kinetic Landau-fluid (LF) closure with collision and periodic boundary condition is used in the development of the deep learning (DL) surrogate model. A classical neural network, namely, feedforward neural network or sometimes termed multilayer perceptron, is constructed and trained to learn the kinetic LF closure in the static limit and arbitrary mean free path in configuration space. The preliminary relation between best hyperparameters and critical parameters for data generation is found. Compared with the numerical approach (non-Fourier method) of the LF closure, the deep learning surrogate model shows an order of magnitude of improvement in terms of accuracy. Perhaps most importantly, the surrogate model closure has been integrated for the first time with fluid simulations. Our DL-enabled fluid simulations, for the first time, give the correct Landau damping rate for a wide range of wave vectors, while the Hammett–Perkins closure cannot produce the correct damping rate. We correctly connect the collisionless Hammett–Perkins closure and collisional Braginskii closure to reproduce the intrinsic nonlocal feature of the heat flux with DL techniques. We address the most concerning error accumulation problem and find that simulations with the deep learning surrogate model are as good as, if not better than, simulations with the analytic closure in terms of long-term numerical stability in the linear Landau damping test.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Intrinsic nonlocality of spin- and polarization-resolved probabilities in strong-field quantum electrodynamics

Spin and polarization are central to precision tests of fundamental physics and for interpreting radiation from astrophysical sources and ultraintense laser-matter experiments. Here, focusing on the fundamental process of nonlinear Compton scattering, we demonstrate that a key assumption underlying current strong-field quantum electrodynamics models, i.e., that emission can be treated as an instantaneous random event sampled from a local differential rate, is inconsistent once emission angles, electron spin, and/or photon polarization are resolved. Namely, even in strictly constant and uniform fields , the resulting fully differential distribution is sign indefinite, yielding negative inferred probabilities. The physical reason is that the photon emission probability builds up over a finite length of the electron trajectory, the formation region, during which the electron direction changes by roughly the same small angle that defines the radiation cone. Therefore, we put forward a new method where we integrate over this formation region analytically to obtain a physically consistent electron spin and photon polarization model. We show that the implementation of our model is compatible with existing Monte Carlo and particle-in-cell workflows. Simulations of a GeV-class electron-laser collision accessible at current petawatt facilities and of emission in a pulsarlike magnetic field are shown to reveal spin and polarization patterns that differ even qualitatively from state-of-the-art local models. In particular, our new model predicts substantial angle-dependent circular photon polarization where the well-known collinear-emission approach yields none, and a pronounced helicity bias in the recoiling electrons absent from current predictions. These findings have direct implications for upcoming strong-field QED experiments and for interpreting polarized radiation from extreme astrophysical environments.

astrophysical electromagnetic fields↗
DOE OSTI2025

Investigating the Relative Roles of INPs and CCN in a Simulated Thunderstorm Using a New Immersion Freezing Algorithm

Microphysical processes in deep convective clouds are sensitive to the number concentrations of cloud condensation nuclei (CCN) and ice nucleating particles (INPs), but the effects of INPs are less studied. Modeling studies investigating the effects of INPs and/or CCN on deep convection typically retain a volume‐dependent raindrop freezing relation. The resulting neglect of aerosol accumulation in raindrops via drop collisions has likely produced unrealistic storm responses to INPs in past studies. To address this deficiency, a new immersion freezing algorithm was developed and embedded in a bulk microphysics scheme that freezes both cloud drops and raindrops using the same immersion freezing INP (IF‐INP) activity spectrum based on measurements. Multiple idealized simulations of a single case of deep convection observed during the Clouds, Aerosols, and Complex Terrain Interactions (CACTI) field campaign were conducted, with microphysical differences produced by independently altering IF‐INP temperature dependencies and CCN number concentrations from their observed values. Surface precipitation in all simulations resulted almost exclusively from riming graupel that melted upon descending to the surface. Rainfall and cold pools were substantially and systematically weakened with increased CCN due to decreased graupel riming rates but were relatively insensitive to variations in the magnitude and slope of IF‐INP spectra due to compensating depletion of supercooled liquid water. These compensating processes were a consequence of the accumulation of IF‐INPs in raindrops, encouraging caution in studying IF‐INP effects upon thunderstorms using traditional volume‐dependent drop freezing relationships.

CCN↗
DOE OSTI2022

Single molecule tracking of bacterial cell surface cytochromes reveals dynamics that impact long-distance electron transport

Using a series of multiheme cytochromes, the metal-reducing bacterium Shewanella oneidensis MR-1 can perform extracellular electron transfer (EET) to respire redox-active surfaces, including minerals and electrodes outside the cell. While the role of multiheme cytochromes in transporting electrons across the cell wall is well established, these cytochromes were also recently found to facilitate long-distance (micrometer-scale) redox conduction along outer membranes and across multiple cells bridging electrodes. Recent studies proposed that long-distance conduction arises from the interplay of electron hopping and cytochrome diffusion, which allows collisions and electron exchange between cytochromes along membranes. However, the diffusive dynamics of the multiheme cytochromes have never been observed or quantified in vivo, making it difficult to assess their hypothesized contribution to the collision-exchange mechanism. Here, we use quantum dot labeling, total internal reflection fluorescence microscopy, and single-particle tracking to quantify the lateral diffusive dynamics of the outer membrane-associated decaheme cytochromes MtrC and OmcA, two key components of EET in S. oneidensis. We observe confined diffusion behavior for both quantum dot-labeled MtrC and OmcA along cell surfaces (diffusion coefficients D MtrC = 0.0192 ± 0.0018 µm 2 /s, D OmcA = 0.0125 ± 0.0024 µm 2 /s) and the membrane extensions thought to function as bacterial nanowires. We find that these dynamics can trace a path for electron transport via overlap of cytochrome trajectories, consistent with the long-distance conduction mechanism. The measured dynamics inform kinetic Monte Carlo simulations that combine direct electron hopping and redox molecule diffusion, revealing significant electron transport rates along cells and membrane nanowires.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Detonation structure in the presence of mixture stratification using reaction-resolved simulations

Many investigations of detonation-based combustors have identified reactant mixture inhomogeneity as having a leading-order impact on wave dynamics and combustion efficiency. To examine this phenomenon in a simplified context, an array of two- and three-dimensional channel detonation simulations are conducted in the present work. The reactant mixture consists of stratified fuel and air, wherein the randomly distributed equivalence ratio field features a characteristic stratification length scale. Detailed chemical kinetics are implemented in an adaptive mesh refinement solution framework where the region near the shock front is resolved with Ο (100) cells per representative ZND induction length. The results show that in comparison to baseline cases with uniform reactant mixtures, reactant stratification has a marked impact on the detonation structure. Increasing the stratification length scale increases the size and irregularity of the detonation cells, yielding larger variations in wave speed. Triple point collisions in fuel-rich regions lead to local wave speeds above the notional mean CJ speed, but wave passage through inert regions causes the local wave speed and strength to diminish. Further, conditional statistics show that increasing the stratification length scale increases the variance in pressure and temperature in the primary reaction zone, as well as the variance in heat release over a range of mixture conditions. In addition to the reactant mixture, the impact of the boundary condition behind the detonation is also investigated. The results show that an inflow boundary condition acts to over-drive the wave, leading to higher peak pressures, smaller detonation cells, and increased reactant consumption. On the other hand, cases with a wall behind the wave exhibit weaker waves with lower peak pressures and heat release rates, as well as greater variance in conditional quantities. Comparisons between complementary two- and three-dimensional simulations show reasonable qualitative agreement in wave structure, speed, and conditional statistics.

42 ENGINEERING↗
DOE OSTI2026

Numerical-heating effects in atmospheric pressure streamer discharges simulated with a PIC code

Artificial heating in plasma simulations is a well-known phenomenon which occurs when, among other things, the Debye length is poorly resolved by the simulation mesh. Here, in this work, the degree to which numerical-heating occurs during a simulation of a nanosecond atmospheric pressure streamer discharge is examined. The streamer is simulated using a two-dimensional finite-element, particle-in-cell code Empire, which uses direct simulation Monte Carlo for binary particle interactions. Initially, an estimate of the numerical-heating rate applied to Empire is performed using a simple plasma model. Second, a positive atmospheric pressure streamer discharge simulation is performed to study the effects of numerical heating on plasma density, electron temperature, and streamer velocity. The nominal Debye length is approximately 1 μm and the amount of numerical heating introduced in the simulation is varied by using mesh sizes ranging from 2 μm to 20 μm. A measurable numerical heating quantity is proposed that can be used to estimate the appropriate element size and quantify the numerical-heating that can be expected over the simulation time for an atmospheric pressure streamer. In conclusion while Δx/λ D violations can be an issue it is not likely to be an issue with streamer discharges that are temporally short and occur in environments where collision frequencies are high. This result validates the rationale of grid size choices for a large amount of previously published works where Δx/λ D violation was not clearly addressed. Primary finding of this work is that numerical heating is of minor concern for plasma simulations where electron–neutral collisions are numerous such that multiple collisions can occur within a single plasma period.

Nikic, Dejan [University of New Mexico, Albuquerqu↗
DOE OSTI2022

First-collision source treatment for ray effect mitigation in discrete-ordinate radiation transport solutions

Deterministic transport codes play a fundamental role in the modelling and simulation of neutron transport. One of the most common deterministic methods is the method of discrete ordinates, also known as the S method. While offering significant advantages over other deterministic methods or stochastic methods like Monte Carlo, the method of discrete ordinates suffers from non-physical artifacts in its local solution due to its discretization of angle. These artifacts, referred to as ray effects because of their ray-like appearance, tend to be worse in problems with small sources in areas with little scattering. Significant effort has gone into developing methods to mitigate ray effects, such as the first-collision source treatment, which separates the angular flux into the uncollided and collided fluxes and solving them using non-traditional techniques such as ray tracing. One such code capable of doing this is Lawrence Livermore National Laboratory's deterministic transport code ARDRA. Current ray tracing methods typically trace to a set of points inside a zone to compute an overall flux. However, this approach has significant drawbacks, such as a low order of convergence and not being conservative. Therefore, a new method has been developed that traces instead to a set of points on each of a zone's surfaces and computing the currents, before using these to obtain the flux. A comparison between these two ray tracing methods showed significant advantages to the new surface method, including inherent conservation, a higher convergence rate, and an increase in calculable information like leakage. This work performed under the auspices of the U.S. Department of Energy by Lawrence Liver- more National Laboratory under Contract DE-AC52-07NA27344. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

CFD Simulation of Helium Flow Loop Test Section

A helium flow loop is being assembled at Oak Ridge National Laboratory to analyze heat transfer enhancement for systems such as blanket and divertor components. To efficiently identify optimum geometries for heat transfer enhancement in these applications, simulation work is performed to optimize test section designs that are built and tested in the helium flow loop that operates at 4 MPa and a mass flow rate of 100 g/s. Different ribbed geometries that examine rib shape, rib height, rib orientation, rib spacing, and three dimensional orientation are modeled and simulated in STAR-CCM+ to compare their ability to remove heat and mitigate pressure drop. Following the simulations, models are selected and manufactured for the helium flow loop tests. Simulations initially focus on a hydrodynamic study to determine the appropriate mesh and physics models and then add a heat flux to analyze the heat transfer abilities of the models. The simulations are run in steady state and use a Reynolds-averaged Navier-Stokes k-ε turbulence model. The helium is modeled as an ideal gas. The simulation explores models of geometries that enhance the heat transfer and decrease pressure drop with an overall goal of increasing fluid collision with the wall. Enhanced geometries are simulated to select appropriate designs for manufacturing, and preliminary experimental results are used to validate the simulations. Furthermore, the factors that are being analyzed in the comparison between the experimental and the simulated results include matching thermocouple temperatures, pressure drop, roughness, and fluid velocity.

42 ENGINEERING↗
DOE OSTI2023

Measurements of the Higgs boson production cross section and couplings in the W boson pair decay channel in proton-proton collisions at $\sqrt{s}=13\,\text {Te\hspace{-.08em}V} $

Production cross sections of the standard model Higgs boson decaying to a pair of W bosons are measured in proton-proton collisions at a center-of-mass energy of 13$\,\text {Te\hspace{-.08em}V}$. The analysis targets Higgs bosons produced via gluon fusion, vector boson fusion, and in association with a W or Z boson. Candidate events are required to have at least two charged leptons and moderate missing transverse momentum, targeting events with at least one leptonically decaying W boson originating from the Higgs boson. Results are presented in the form of inclusive and differential cross sections in the simplified template cross section framework, as well as couplings of the Higgs boson to vector bosons and fermions. The data set collected by the CMS detector during 2016–2018 is used, corresponding to an integrated luminosity of 138$\,\text {fb}^{-1}$. The signal strength modifier $\mu $, defined as the ratio of the observed production rate in a given decay channel to the standard model expectation, is measured to be $\mu = 0.95^{+0.10}_{-0.09}$. All results are found to be compatible with the standard model within the uncertainties.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Advancing simulations of coupled electron and phonon nonequilibrium dynamics using adaptive and multirate time integration

Electronic structure calculations in the time domain provide a deeper understanding of nonequilibrium dynamics in materials. The real-time Boltzmann equation (rt-BTE), used in conjunction with accurate interactions computed from first principles, has enabled reliable predictions of coupled electron and lattice dynamics. However, the timescales and system sizes accessible with this approach are still limited, with two main challenges being the different timescales of electron and phonon interactions and the cost of computing collision integrals. As a result, only a few examples of these calculations exist, mainly for two-dimensional (2D) materials. Here we leverage adaptive and multirate time integration methods to achieve a major step forward in solving the coupled rt-BTEs for electrons and phonons. Relative to conventional (non-adaptive) time-stepping, our approach achieves a 10x speedup for a target accuracy, or greater accuracy by 3–6 orders of magnitude for the same computational cost, enabling efficient calculations in both 2D and bulk materials. This efficiency is showcased by computing the coupled electron and lattice dynamics in graphene up to ~100 ps, as well as modeling ultrafast lattice dynamics and thermal diffuse scattering maps in bulk materials (silicon and gallium arsenide). In addition to improved efficiency, our adaptive method can resolve the characteristic rates of different physical processes, thus naturally bridging different timescales. This enables simulations of longer timescales and provides a framework for modeling multiscale dynamics of coupled degrees of freedom in matter. Our work opens new opportunities for quantitative studies of nonequilibrium physics in materials, including driven lattice dynamics with phonons coupled to electrons, spin, and other degrees of freedom.

Yao, Jia [California Institute of Technology (CalT↗
DOE OSTI2023

Charmonium production in $p$Ne collisions at $\sqrt{s_{\scriptscriptstyle \text {NN}}} =68.5$ GeV

The measurement of charmonium states produced in proton-neon (pNe) collisions by the LHCb experiment in its fixed-target configuration is presented. The pro duction of J/ψ and ψ(2S) mesons is studied with a beam of 2.5 TeV protons colliding on gaseous neon targets at rest, corresponding to a nucleon-nucleon centre-of-mass energy $\sqrt{s_{NN}}$ = 68.5 GeV. The data sample corresponds to an integrated luminosity of 21.7 ± 1.4 nb -1 . The J/ψ and ψ(2S) hadrons are reconstructed in μ + μ - final states. The J/ψ production cross-section per target nucleon in the centre-of mass rapidity range y* ϵ [-2.29, 0] is found to be 506 ± 8 ± 46 nb/nucleon. The ratio of J/ψ and D 0 cross-sections is evaluated to (1.06 ± 0.02 ± 0.09)%. The ψ(2S) to J/ψ relative production rate is found to be (1.67 ± 0.27 ± 0.10)% in good agreement with other measurements involving beam and target nuclei of similar sizes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Forward neutrinos from charm at the Large Hadron Collider

The currently operating FASER experiment and the planned Forward Physics Facility (FPF) will detect a large number of neutrinos produced in proton-proton collisions at the LHC. In this work, we estimate neutrino fluxes at these detectors from charm meson decays, which will be particularly important for the ν e and ν τ channels. We make prediction using both the next-to-leading order collinear factorization and the k T -factorization approaches to model the production of charm quarks as well as different schemes to model their hadronization into charm hadrons. In particular, we emphasize that a sophisticated modeling of hadronization involving beam remnants is needed for predictions at FASER and FPF due to the sensitivity to the charm hadron production at low transverse momenta and very forward rapidity. As example, we use the string fragmentation approach implemented in 8. While both standard fragmentation functions and 8 are able to describe LHCb data, we find that 8 predicts significantly higher rate of high energy neutrinos, highlighting the importance of using the correct hadronization model when making predictions. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2026

Mapping Rydberg States of H 2 with the Halfium R-Matrix Method

In this article, we use the Halfium R-matrix method to investigate the Rydberg states of the H 2 molecule up to n = 20, filling the gap above the low-lying bound states already calculated with configuration interaction packages. Moreover, we show that the use of Quantum Defect Theory scaling laws, allows for a comprehensive analysis of the regular patterns resulting from the coupling between Rydberg series and doubly excited states. The results should open the door for more efficient quasi-diabatization of the potential energy curves which is required for calculating cross sections and rate coefficients of the (e + H 2 + ) collisional processes, involved in the plasma modeling for fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Measurements of Higgs boson production cross sections and couplings in the diphoton decay channel at $ \sqrt{\mathrm{s}} $ = 13 TeV

Measurements of Higgs boson production cross sections and couplings in events where the Higgs boson decays into a pair of photons are reported. Events are selected from a sample of proton-proton collisions at s√ = 13 TeV collected by the CMS detector at the LHC from 2016 to 2018, corresponding to an integrated luminosity of 137 fb -1 . Analysis categories enriched in Higgs boson events produced via gluon fusion, vector boson fusion, vector boson associated production, and production associated with top quarks are constructed. The total Higgs boson signal strength, relative to the standard model (SM) prediction, is measured to be 1.12±0.09. Other properties of the Higgs boson are measured, including SM signal strength modifiers, production cross sections, and its couplings to other particles. These include the most precise measurements of gluon fusion and vector boson fusion Higgs boson production in several different kinematic regions, the first measurement of Higgs boson production in association with a top quark pair in five regions of the Higgs boson transverse momentum, and an upper limit on the rate of Higgs boson production in association with a single top quark. All results are found to be in agreement with the SM expectations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Study of the f 0 ( 980 ) and f 0 ( 500 ) Scalar Mesons through the Decay D s + → π + π − e + ν e

Using e + e − collision data corresponding to an integrated luminosity of 7.33 fb − 1 recorded by the BESIII detector at center-of-mass energies between 4.128 and 4.226 GeV, we present an analysis of the decay D s + → π + π − e + ν e , where the D s + is produced via the process e + e − → D s * ± D s ∓ . We observe the f 0 ( 980 ) in the π + π − system and the branching fraction of the decay D s + → f 0 ( 980 ) e + ν e with f 0 ( 980 ) → π + π − measured to be ( 1.72 ± 0.1 3 stat ± 0.1 0 syst ) × 10 − 3 , where the uncertainties are statistical and systematic, respectively. The dynamics of the D s + → f 0 ( 980 ) e + ν e decay are studied with the simple pole parametrization of the hadronic form factor and the Flatté formula describing the f 0 ( 980 ) in the differential decay rate, and the product of the form factor f + f 0 ( 0 ) and the c → s Cabibbo-Kobayashi-Maskawa matrix element | V c s | is determined for the first time to be f + f 0 ( 0 ) | V c s | = 0.504 ± 0.01 7 stat ± 0.03 5 syst . Furthermore, the decay D s + → f 0 ( 500 ) e + ν e is searched for the first time but no signal is found. The upper limit on the branching fraction of D s + → f 0 ( 500 ) e + ν e , f 0 ( 500 ) → π + π − decay is set to be 3.3 × 10 − 4 at 90% confidence level. Published by the American Physical Society 2024

Ablikim, M.↗
DOE OSTI2021

Quarkonium semiclassical transport in quark-gluon plasma: factorization and quantum correction

We study quarkonium transport in the quark-gluon plasma by using the potential nonrelativistic QCD (pNRQCD) effective field theory and the framework of open quantum systems. We argue that the coupling between quarkonium and the thermal bath is weak using separation of scales, so the initial density matrix of the total system factorizes and the time evolution of the subsystem is Markovian. We derive the semiclassical Boltzmann equation for quarkonium by applying a Wigner transform to the Lindblad equation and carrying out a semiclassical expansion. We resum relevant interactions to all orders in the coupling constant at leading power of the nonrelativistic and multipole expansions. The derivation is valid for both weakly coupled and strongly coupled quark-gluon plasmas. We find reaction rates in the transport equation factorize into a quarkonium dipole transition function and a chromoelectric gluon distribution function. For the differential reaction rate, the definition of the momentum dependent chromoelectric gluon distribution function involves staple-shaped Wilson lines. For the inclusive reaction rate, the Wilson lines collapse into a straight line along the real time axis and the distribution becomes momentum independent. The relation between the two Wilson lines is analogous to the relation between the Wilson lines appearing in the gluon parton distribution function (PDF) and the gluon transverse momentum dependent parton distribution function (TMDPDF). The centrality dependence of the quarkonium nuclear modification factor measured by experiments probes the momentum independent distribution while the transverse momentum dependence and measurements of the azimuthal angular anisotropy may be able to probe the momentum dependent one. We discuss one way to indirectly constrain the quarkonium in-medium real potential by using the factorization formula and lattice calculations. The leading quantum correction to the semiclassical transport equation of quarkonium is also worked out. The study can be easily generalized to quarkonium transport in cold nuclear matter, which is relevant for quarkonium production in eA collisions in the future Electron-Ion Collider.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗