Search NASA⌕ Search

Engineering topics

Strickland, Michael

Publications and source records attributed to Strickland, Michael.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Bottomonium suppression from the three-loop QCD potential

We compute the suppression of bottomonium in the quark-gluon plasma using the three-loop QCD static potential. The potential describes the spin-averaged bottomonium spectrum below threshold with a less than 1% error. Within potential nonrelativistic quantum chromodynamics and an open quantum systems framework, we compute the evolution of the bottomonium density matrix. The values of the quarkonium transport coefficients are obtained from lattice QCD measurements of the bottomonium in-medium width and thermal mass shift; we additionally include for the first time a vacuum contribution to the dispersive coefficient γ . Using the three-loop potential and the values of the heavy quarkonium transport coefficients, we find that the resulting bottomonium nuclear modification factor is consistent with experimental observations, while at the same time reproducing the lattice measurements of the in-medium width. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Theoretical and experimental constraints for the equation of state of dense and hot matter

Abstract This review aims at providing an extensive discussion of modern constraints relevant for dense and hot strongly interacting matter. It includes theoretical first-principle results from lattice and perturbative QCD, as well as chiral effective field theory results. From the experimental side, it includes heavy-ion collision and low-energy nuclear physics results, as well as observations from neutron stars and their mergers. The validity of different constraints, concerning specific conditions and ranges of applicability, is also provided.

Kumar, Rajesh (ORCID:0000000327463956)↗
DOE OSTI2024

Non-equilibrium Dynamics of the Quark Gluon Plasma (Final Report)

During this funding period I continued my work on far-from-equilibrium formulations of hydrodynamics and bottomonium suppression in the quark-gluon plasma. I published 23 peer-reviewed articles and 9 proceedings contributions. The publications produced during this period have thus far received 705 citations. I delivered 22 talks, seminars, and colloquia about my work, most of which were invited. In addition to the work on hydrodynamics and heavy quarkonia, I also worked on high-order perturbative calculations of N=4 supersymmetric Yang-Mills thermodynamics and contributed to the efforts of the MUSES collaboration whose focus is developing a modular equation-of-state framework for use in modeling neutron stars and neutron star mergers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Bottomonium suppression in 5.02 and 8.16 TeV 𝑝-Pb collisions

Here, we compute the suppression of ϒ⁡(1⁢𝑆), ϒ⁡(2⁢𝑆), and ϒ⁡(3⁢𝑆) states in 𝑝-Pb collisions relative to 𝑝⁢𝑝 collisions, including nuclear parton distribution function (nPDF) effects, coherent energy loss, momentum broadening, and final-state interactions in the quark-gluon plasma. We employ the EPPS21 nPDFs and calculate the uncertainty resulting from variation over the associated error sets. To compute coherent energy loss and momentum broadening, we follow the approach of Arleo, Peigne, and collaborators. The 3+1⁢D viscous hydrodynamical background evolution of the quark-gluon plasma is generated by anisotropic hydrodynamics. The in-medium suppression of bottomonium in the quark-gluon plasma is computed using a next-to-leading-order open quantum system framework formulated within potential nonrelativistic quantum chromodynamics. We find that inclusion of all these effects provides a reasonable description of experimental data from the ALICE, ATLAS, CMS, and LHCb Collaborations for the suppression of ϒ⁡(1⁢𝑆), ϒ⁡(2⁢𝑆), and ϒ⁡(3⁢𝑆) as a function of both transverse momentum and rapidity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Predictions for the sPHENIX physics program

sPHENIX is a next-generation detector experiment at the Relativistic Heavy Ion Collider, designed for a broad set of jet and heavy-flavor probes of the Quark-Gluon Plasma created in heavy ion collisions. In anticipation of the commissioning and first data-taking of the detector in 2023, a RIKEN-BNL Research Center (RBRC) workshop was organized to collect theoretical input and identify compelling aspects of the physics program. Here this paper compiles theoretical predictions from the workshop participants for jet quenching, heavy flavor and quarkonia, cold QCD, and bulk physics measurements at sPHENIX.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Transverse momentum dependent feed-down fractions for bottomonium production

We extract transverse-momentum-dependent feed-down fractions for bottomonium production using a data-driven approach. We use data published by the ATLAS, CMS, and LHCb Collaborations for √𝑠 =7 TeV proton-proton collisions. Based on this collected data, we produce fits to the differential cross sections for the production of both 𝑆- and 𝑃-wave bottomonium states. Combining these fits with branching ratios for excited state decays from the Particle Data Group, we compute the feed-down fractions for both the ϒ⁡(1⁢𝑆) and ϒ⁡(2⁢𝑆) as a function of transverse momentum. Our results indicate a strong dependence on transverse momentum, which is consistent with prior extractions of the feed-down fractions. When evaluated at the average momentum of the states, we find that approximately 75% of ϒ⁡(1⁢𝑆) and ϒ⁡(2⁢𝑆) states are produced directly. Our results for the transverse-momentum-dependent feed-down fractions are provided in tabulated form so that they can be used by other research groups.

Astronomy & Astrophysics↗
DOE OSTI2023

Kaonic Hanbury-Brown–Twiss radii at 200 GeV and 5.02 TeV

Here, we use 3 + 1 dimensional quasiparticle anisotropic hydrodynamics (aHydroQP) to make predictions for kaon Hanbury-Brown–Twiss (HBT) radii in 200 GeV and 5.02 TeV heavy-ion collisions. Using previously determined aHydroQP parameters, we compute kaonic HBT radii and their ratios as a function of the mean transverse momentum of the pair k T . We first consider Au-Au collisions at 200 GeV, finding good agreement between aHydroQP predictions and experimental data up to k T ≈ 0.8 GeV. We then present predictions for kaonic HBT radii and their ratios in 5.02 TeV Pb-Pb collisions. Our aHydroQP predictions do not exhibit a clear k T scaling of the pion and kaon source radii, however, an approximate transverse mass m T scaling is observed, particularly at 200 GeV.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Nuclear Physics Exascale Requirements Review: An Office of Science Review sponsored jointly by Advanced Scientific Computing Research and Nuclear Physics, June 15 - 17, 2016, Gaithersburg, Maryland

Imagine being able to predict — with unprecedented accuracy and precision — the structure of the proton and neutron, and the forces between them, directly from the dynamics of quarks and gluons, and then using this information in calculations of the structure and reactions of atomic nuclei and of the properties of dense neutron stars (NSs). Also imagine discovering new and exotic states of matter, and new laws of nature, by being able to collect more experimental data than we dream possible today, analyzing it in real time to feed back into an experiment, and curating the data with full tracking capabilities and with fully distributed data mining capabilities. Making this vision a reality would improve basic scientific understanding, enabling us to precisely calculate, for example, the spectrum of gravity waves emitted during NS coalescence, and would have important societal applications in nuclear energy research, stockpile stewardship, and other areas. This review presents the components and characteristics of the exascale computing ecosystems necessary to realize this vision.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Nonextensive hydrodynamics of boost-invariant plasmas

We use quasiparticle anisotropic hydrodynamics to study the non-conformal and non-extensive dynamics of a system undergoing boost-invariant Bjorken expansion. To introduce nonextensivity, we use an underlying Tsallis distribution with a time-dependent nonextensivity parameter q. By taking moments of the quasiparticle Boltzmann equation in the relaxation-time approximation, we obtain dynamical equations which allow us to determine the time evolution of all microscopic parameters including q. We compare numerical solutions for bulk observables obtained using the nonextensive evolution with results obtained using quasiparticle anisotropic hydrodynamics with a Boltzmann distribution function (q → 1). We show that the evolution of the temperature, pressure ratio, and scaled energy density, are quite insensitive to which distribution function is assumed. However, we find significant differences in the early-time evolution of the bulk pressure which are observed for even small deviations from the Boltzmann distribution function. Finally, we discuss the existence of non-conformal hydrodynamic attractors for the longitudinal and transverse pressures, the bulk and shear viscous corrections, and the nonextensivity parameter q.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

The complex heavy-quark potential in an anisotropic quark-gluon plasma — Statics and dynamics

We generalize a complex heavy-quark potential model from an isotropic QCD plasma to an anisotropic one by replacing the Debye mass m> D with an anisotropic screening mass depending on the quark pair alignment with respect to the direction of anisotropy. Such an angle-dependent mass is determined by matching the perturbative contributions in the potential model to the exact result obtained in the Hard-Thermal-Loop resummed perturbation theory. An advantage of the resulting potential model is that its angular dependence can be effectively described by using a set of angle-averaged screening masses as proposed in our previous work. Consequently, one could solve a one-dimensional Schrödinger equation with a potential model built by changing the anisotropic screening masses into the corresponding angle-averaged ones, and reproduce the full three-dimensional results for the binding energies and decay widths of low-lying quarkonium bound states to very high accuracy. Finally, turning to dynamics, we demonstrate that the one-dimensional effective potential can accurately describe the time evolution of the vacuum overlaps obtained using the full three-dimensional anisotropic potential. This includes the splitting of different p-wave polarizations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

QTRAJ 1.0: A Lindblad equation solver for heavy-quarkonium dynamics

Here we introduce an open-source package called QTraj that solves the Lindblad equation for heavy-quarkonium dynamics using the quantum trajectories algorithm. The package allows users to simulate the suppression of heavy-quarkonium states using externally-supplied input from 3+1D hydrodynamics simulations. The code uses a split-step pseudo-spectral method for updating the wave-function between jumps, which is implemented using the open-source multi-threaded FFTW3 package. This allows one to have manifestly unitary evolution when using real-valued potentials. In this paper, we provide detailed documentation of QTraj 1.0, installation instructions, and present various tests and benchmarks of the code.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Bulk observables at 5.02 TeV using quasiparticle anisotropic hydrodynamics

Abstract We compare predictions of 3+1D quasiparticle anisotropic hydrodynamics (aHydroQP) for a large set of bulk observables with experimental data collected in 5.02 TeV Pb–Pb collisions. We make predictions for identified hadron spectra, identified hadron average transverse momentum, charged particle multiplicity as a function of pseudorapidity, the kaon-to-pion ( $$K/\pi $$ K / π ) and proton-to-pion ( $$p/\pi $$ p / π ) ratios, identified particle and charged particle elliptic flow, and HBT radii. We compare to data collected by the ALICE collaboration in 5.02 TeV Pb–Pb collisions. We find that, based on available data, these bulk observables are well described by aHydroQP with an assumed initial central temperature of $$T_0=630$$ T 0 = 630 MeV at $$\tau _0 = 0.25$$ τ 0 = 0.25 fm/c and a constant specific shear viscosity of $$\eta /s=0.159$$ η / s = 0.159 , which corresponds to a peak specific bulk viscosity of $$\zeta /s = 0.048$$ ζ / s = 0.048 . In particular, we find that the momentum dependence of the kaon-to-pion ( $$K/\pi $$ K / π ) and proton-to-pion ( $$p/\pi $$ p / π ) ratios reported recently by the ALICE collaboration are extremely well described by aHydroQP in the 0–5% centrality class.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

The imaginary part of the heavy-quark potential from real-time Yang-Mills dynamics

We extract the imaginary part of the heavy-quark potential using classical-statistical simulations of real-time Yang-Mills dynamics in classical thermal equilibrium. The r-dependence of the imaginary part of the potential is extracted by measuring the temporal decay of Wilson loops of spatial length r. We compare our results to continuum expressions obtained using hard thermal loop theory and to semi-analytic lattice perturbation theory calculations using the hard classical loop formalism. We find that, when plotted as a function of m D r, where m D is the hard classical loop Debye mass, the imaginary part of the heavy-quark potential shows little sensitivity to the lattice spacing at small m D r ≲ 1 and agrees well with the semi-analytic hard classical loop result. For large quark-antiquark separations, we quantify the magnitude of the non-perturbative long-range corrections to the imaginary part of the heavy-quark potential. We present our results for a wide range of temperatures, lattice spacings, and lattice volumes. This work sets the stage for extracting the imaginary part of the heavy-quark potential in an expanding non-equilibrium Yang Mills plasma.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

\( \mathcal{N} \) = 4 supersymmetric Yang-Mills thermodynamics to order λ 2

We calculate the resummed perturbative free energy of \( \mathcal{N} \) = 4 supersymmetric Yang-Mills in four spacetime dimensions (SYM 4 , 4 ) through second order in the ’t Hooft coupling λ at finite temperature and zero chemical potential. Our final result is ultraviolet finite and all infrared divergences generated at three-loop level are canceled by summing over SYM 4 , 4 ring diagrams. Non-analytic terms at \( \mathcal{O} \) ( λ 3 / 2 ) and \( \mathcal{O} \) ( λ 2 log λ ) are generated by dressing the A 0 and scalar propagators. The gauge-field Debye mass m D and the scalar thermal mass M D are determined from their corresponding finite-temperature self-energies. Based on this, we obtain the three-loop thermodynamic functions of SYM 4 , 4 to \( \mathcal{O} \) ( λ 2 ). We compare our final result with prior results obtained in the weak- and strong-coupling limits and construct a generalized Padé approximant that interpolates between the weak-coupling result and the large- N c strong-coupling result. Our results suggest that the \( \mathcal{O} \) ( λ 2 ) weak-coupling result for the scaled entropy density is a quantitatively reliable approximation to the scaled entropy density for 0 ≤ λ ≲ 2.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Bottomonium suppression in an open quantum system using the quantum trajectories method

We solve the Lindblad equation describing the Brownian motion of a Coulombic heavy quark-antiquark pair in a strongly coupled quark-gluon plasma using the highly efficient Monte Carlo wave-function method. The Lindblad equation has been derived in the framework of pNRQCD and fully accounts for the quantum and non-Abelian nature of the system. The hydrodynamics of the plasma is realistically implemented through a 3+1D dissipative hydrodynamics code. We compute the bottomonium nuclear modification factor and compare with the most recent LHC data. The computation does not rely on any free parameter, as it depends on two transport coefficients that have been evaluated independently in lattice QCD. Our final results, which include late-time feed down of excited states, agree well with the available data from LHC 5.02 TeV PbPb collisions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Bottomonium suppression and elliptic flow using Heavy Quarkonium Quantum Dynamics

We introduce a framework called Heavy Quarkonium Quantum Dynamics (HQQD) which can be used to compute the dynamical suppression of heavy quarkonia propagating in the quark-gluon plasma using real-time in-medium quantum evolution. Using HQQD we compute large sets of real-time solutions to the Schrödinger equation using a realistic in-medium complex-valued potential. We sample 2 million quarkonia wave packet trajectories and evolve them through the QGP using HQQD to obtain their survival probabilities. The computation is performed using three different HQQD model parameter sets in order to estimate our systematic uncertainty. After taking into account final state feed down we compare our results to existing experimental data for the suppression and elliptic flow of bottomonium states and find that HQQD predictions are good agreement with available data for R AA as a function of N part and p T collected at \( \sqrt{s_{\mathrm{NN}}} \) = 5 . 02 TeV. In the case of v 2 for the various states, we find that the path-length dependence of ϒ(1 s ) suppression results in quite small v 2 for ϒ(1 s ). Our prediction for the integrated elliptic flow for ϒ(1 s ) in the 10 – 90% centrality class, which now includes an estimate of the systematic error, is v 2 [ϒ(1 s )] = 0 . 003 ± 0 . 0007 ± \( {}_{0.0013}^{0.0006} \) . We also find that, due to their increased suppression, excited bottomonium states have a larger elliptic flow. Based on this observation we make predictions for v 2 [ϒ(2 s )] and v 2 [ϒ(3 s )] as a function of centrality and transverse momentum.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Pion interferometry at 200 GeV using anisotropic hydrodynamics

In this paper, we continue our phenomenological studies of heavy-ion collisions using 3 + 1 D anisotropic hydrodynamics (aHydro). In previous works, we compared quasiparticle aHydro (aHydroQP) with ALICE 2.76 TeV Pb-Pb and BNL Relativistic Heavy Ion Collider 200 GeV Au-Au collision results. At both energies, the agreement was quite good between aHydroQP and the experimental data for many observables. In this work, we present comparisons of the Hanbury Brown–Twiss (HBT) radii and their ratios determined using π + π + pairs produced in 200 GeV Au-Au collisions. We first present comparisons with STAR results for the HBT radii and their ratios. We then present comparisons with PHENIX results for the HBT radii and their ratios. In both cases, we find reasonable agreement between aHydroQP predictions and available experimental results for the ratios of HBT radii. At the level of the radii themselves, in some cases quantitative differences on the order of 10–20% remain, which deserve further research.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗