Researcher profile

Michael Strickland

Michael Strickland contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
22works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

22 published item(s)

preprint2023arXiv

3D Multi-system Bayesian Calibration with Energy Conservation to Study Rapidity-dependent Dynamics of Nuclear Collisions

Considerable information about the early-stage dynamics of heavy-ion collisions is encoded in the rapidity dependence of measurements. To leverage the large amount of experimental data, we perform a systematic analysis using three-dimensional hydrodynamic simulations of multiple collision systems -- large and small, symmetric and asymmetric. Specifically, we perform fully 3D multi-stage hydrodynamic simulations initialized by a parameterized model for rapidity-dependent energy deposition, which we calibrate on the hadron multiplicity and anisotropic flow coefficients. We utilize Bayesian inference to constrain properties of the early- and late- time dynamics of the system, and highlight the impact of enforcing global energy conservation in our 3D model.

preprint2023arXiv

Bottomonium suppression in the quark-gluon plasma -- From effective field theories to non-unitary quantum evolution

In this proceedings contribution I review recent work which computes the suppression of bottomonium production in heavy-ion collisions using open quantum systems methods applied within the potential non-relativistic quantum chromodynamics (pNRQCD) effective field theory. I discuss how the computation of bottomonium suppression can be reduced to solving a Gorini-Kossakowski-Sudarshan-Lindblad (GKSL) quantum master equation for the evolution of the b-bbar reduced density matrix. The open quantum systems approach used allows one to take into account the non-equilibrium dynamics and decoherence of bottomonium in the quark-gluon plasma. Finally, I present comparisons of phenomenological predictions obtained using a recently obtained next-to-leading-order GKSL equation with ALICE, ATLAS, and CMS experimental data for bottomonium suppression and elliptic flow.

preprint2022arXiv

${\cal N}=4$ supersymmetric Yang-Mills thermodynamics from effective field theory

The free energy density of ${\cal N}=4$ supersymmetric Yang-Mills theory in four space-time dimensions is derived through second order in the 't Hooft coupling $λ$ at finite temperature using effective-field theory methods. The contributions to the free energy density at this order come from the hard scale $T$ and the soft scale $\sqrtλ T$. The effects of the scale $T$ are encoded in the coefficients of an effective three-dimensional field theory that is obtained by dimensional reduction at finite temperature. The effects of the electric scale $\sqrtλ T$ are taken into account by perturbative calculations in the effective theory.

preprint2022arXiv

Effective field theory treatment of ${\cal N}=4$ supersymmetric Yang-Mills thermodynamics

At finite temperature the free energy density of ${\cal N}=4$ supersymmetric Yang-Mills can be calculated using resummed perturbation theory through the order $λ^{5/2}$. Effective field theory methods provide a useful alternative approach to streamline these calculations. In this proceedings contribution, I review recent work with my collaborators where we used effective field theory methods to calculate the free energy density of ${\cal N}=4$ supersymmetric Yang-Mills in four spacetime dimensions through second order in the 't Hooft coupling $λ$. At this order the contributions to the free energy density come from the hard scale $T$ and the soft scale $\sqrtλT$. The contribution from the scale $T$ enters through the coefficients in the effective Lagrangian obtained by dimensional reduction and the effects of the scale $gT$ can be calculated using perturbative methods in the effective theory.

preprint2022arXiv

Long and short distance behavior of the imaginary part of the heavy-quark potential

The imaginary part of the effective heavy-quark potential is related to the total in-medium decay width of heavy quark-antiquark bound states. We extract the static limit of this quantity using classical-statistical simulations of real-time Yang-Mills dynamics by measuring the temporal decay of Wilson loops. By performing the simulations on finer and larger lattices, we are able to show that the nonperturbative results follow the same form as the perturbative ones. 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. We also extract approximations for the short-distance behavior of the classical potential.

preprint2022arXiv

Non-equilibrium attractor in high-temperature QCD plasmas

We establish the existence of a far-from-equilibrium attractor in weakly-coupled gauge theory undergoing 0+1d Bjorken expansion which goes beyond the energy-momentum tensor to the detailed form of the one-particle distribution function. We then demonstrate that the dynamics can be re-scaled at intermediate times and represented by universal exponents. Finally, we assess different procedures for reconstructing the full one-particle distribution function from the energy-momentum tensor along the attractor and discuss implications for the freeze-out procedure used in the phenomenological analysis of ultra-relativistic nuclear collisions

preprint2022arXiv

Non-equilibrium evolution of quarkonium in medium in the open quantum system approach

In this proceedings contribution, I review recent work that aims to provide a more comprehensive and systematic understanding of bottomonium dynamics in the quark-gluon plasma using an open quantum system (OQS) approach that is applied in the framework of the potential non-relativistic QCD (pNRQCD) effective field theory and coupled to realistic hydrodynamical backgrounds that have been tuned to soft hadron observables. I review how the computation of bottomonium suppression can be reduced to solving a Gorini-Kossakowski-Sudarshan-Lindblad (GKSL) equation for the evolution of the b-bbar reduced density matrix, which includes both singlet and octet states plus medium-induced transitions between them at next-to-leading order (NLO) in the binding energy over temperature. Finally, I present comparisons of phenomenological predictions of the NLO OQS+pNRQCD approach and experimental data for bottomonium suppression and elliptic flow in LHC 5.02 TeV Pb-Pb collisions obtained using both smooth and fluctuating hydrodynamic initial conditions.

preprint2022arXiv

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

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.

preprint2022arXiv

QuantumFDTD -- A computational framework for the relativistic Schrödinger equation

We extend the publicly available quantumfdtd code. It was originally intended for solving the time-independent three-dimensional Schrödinger equation via the finite-difference time-domain (FDTD) method and for extracting the ground, first, and second excited states. We (a) include the case of the relativistic Schrödinger equation and (b) add two optimized FFT-based kinetic energy terms for the non-relativistic case. All the three new kinetic terms are computed using Fast Fourier Transform (FFT). We release the resulting code as version 3 of quantumfdtd. Finally, the code now supports arbitrary external file-based potentials and the option to project out distinct parity eigenstates from the solutions. Our goal is quark models used for phenomenological descriptions of QCD bound states, described by the three-dimensional Schrödinger equation. However, we target any field where solving either the non-relativistic or the relativistic three-dimensional Schrödinger equation is required.

preprint2022arXiv

Scheme dependence of two-loop HTLpt-resummed $\text{SYM}_{4,4}$ thermodynamics

The resummed thermodynamics of ${\cal N}=4$ supersymmetric Yang-Mills theory in four space-time dimensions ($\text{SYM}_{4,4}$) has been calculated previously to two loop order within hard thermal loop perturbation theory (HTLpt) using the canonical dimensional regularization (DRG) scheme. Herein, we revisit this calculation using the regularization by dimensional reduction (RDR) scheme. Since the RDR scheme manifestly preserves supersymmetry it is the preferred scheme, however, it is important to assess if and by how much the resummed perturbative results depend on the regularization scheme used. Comparing predictions for the scaled entropy obtained using the DRG and RDR schemes we find that for $λ\lesssim 6$ they are numerically very similar. We then compare the results obtained in both schemes with the strict perturbative result, which is accurate up to order $λ^2$, and a generalized Padé approximant constructed from the known large-$N_c$ weak- and strong-coupling expansions. Comparing the strict perturbative expansion of the two-loop HTLpt result with the perturbative expansion to order $λ^2$, we find that both the DRG and RDR HTLpt calculations result in the same scheme-independent predictions for the coefficients at order $λ$, $λ^{3/2}$, and $λ^2 \logλ$, however, at order $λ^2$ there is a residual regularization scheme dependence.

preprint2022arXiv

Some recent advances in the understanding of ${\cal N}=4$ supersymmetric Yang-Mills thermodynamics

The interest in the thermodynamics of supersymmetric Yang-Mills started after Maldacena proposed the duality between string theory on AdS backgrounds and the large-N limit of SYM theories. One of the motivations to study the thermal properties of ${\cal N}=4$ supersymmetric Yang-Mills in four dimensions (SYM$_{4,4}$) is that at high temperatures, the weak-coupling limit of this theory has many similarities with high temperature quantum chromodynamics (QCD). In this proceedings contribution, we review recent calculations of the resummed perturbative free energy of ${\cal N}=4$ supersymmetric Yang-Mills in four spacetime dimensions through second order in the 't Hooft coupling $λ$ at finite temperature and zero chemical potential. 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.

preprint2022arXiv

The effective complex heavy-quark potential in an anisotropic quark-gluon plasma

We introduce a method for reducing anisotropic heavy-quark potentials to isotropic potentials by using an effective screening mass that depends on the quantum numbers $l$ and $m$ of a given state. We demonstrate that, using the resulting 1D effective potential model, one can solve a 1D Schrödinger equation and reproduce the full 3D results for the energies and binding energies of low-lying heavy-quarkonium bound states to relatively high accuracy. This includes the splitting of different p-wave polarizations. The resulting 1D effective model provides a way to include momentum anisotropy effects in open quantum system simulations of heavy-quarkonium dynamics in the quark-gluon plasma.

preprint2022arXiv

The non-equilibrium attractor for kinetic theory in relaxation time approximation

I demonstrate that the concept of a non-equilibrium attractor can be extended beyond the lowest-order moments typically considered in hydrodynamic treatments. Using a previously obtained exact solution to the relaxation-time approximation Boltzmann equation for a transversally homogeneous and boost-invariant system subject to Bjorken flow, I derive an equation obeyed by all moments of the one-particle distribution function. Using numerical solutions, I show that, similar to the pressure anisotropy, all moments of the distribution function exhibit attractor-like behavior wherein all initial conditions converge to a universal solution after a short time with the exception of moments which are sensitive to modes with zero longitudinal momentum and high transverse momentum. In addition, I compute the exact solution for the distribution function itself on very fine lattices in momentum space and demonstrate that (a) an attractor for the full distribution function exists and (b) solutions with generic initial conditions relax to this solution, first at low momentum and later at high momentum.

preprint2021arXiv

Bottomonium suppression and elliptic flow in heavy-ion collisions

In this proceedings contribution I review recent progress concerning the suppression of bottomonium production in the quark-gluon plasma. Making use of open quantum system methods applied to potential non-relativistic quantum chromodynamics one can show that the dynamics of heavy-quarkonium bound states satisfying the scale hierarchy 1/a_0 >> pi T ~ m_D >> E obey a Lindblad equation whose solution provides the quantum evolution of the heavy-quarkonium reduced density matrix. To solve the resulting Lindblad equation we use a quantum trajectories algorithm which allows one to include all possible angular momentum states of the quark-antiquark probe in a scalable manner. We solve the Lindblad equation using a tuned 3+1D dissipative hydrodynamics code for the background temperature evolution. We then consider a large number of Monte-Carlo sampled bottomonium trajectories embedded in this background. This allows us to extract the centrality- and p_T-dependence of the nuclear suppression factor R_AA[Upsilon] and elliptic flow v_2[Upsilon]. We find good agreement between our model predictions and available sqrt(s_NN) = 5.02 TeV Pb-Pb collision experimental data from the ALICE, ATLAS, and CMS collaborations.

preprint2021arXiv

Bottomonium suppression and flow in heavy-ion collisions

The strong suppression of bottomonia production in ultra-relativistic heavy-ion collisions is a smoking gun for the creation of a deconfined quark-gluon plasma (QGP). In this proceedings contribution, I review recent work that aims to provide a more comprehensive and systematic understanding of bottomonium dynamics in the QGP through the use of pNRQCD and an open quantum systems approach. This approach allows one to evolve the heavy-quarkonium reduced density matrix, taking into account non-unitary effective Hamiltonian evolution of the wave-function and quantum jumps between different angular momentum and color states. In the case of a strong coupled QGP in which E << T,m_D << 1/a_0, the corresponding evolution equation is Markovian and can therefore be mapped to a Lindblad evolution equation. To solve the resulting Lindblad equation, we make use of a stochastic unraveling called the quantum trajectories algorithm and couple the non-abelian quantum evolution to a realistic 3+1D viscous hydrodynamical background. Using a large number of Monte-Carlo sampled bottomonium trajectories, we make predictions for bottomonium R_AA and elliptic flow as a function of centrality and transverse momentum and compare to data collected by the ALICE, ATLAS, and CMS collaborations.

preprint2021arXiv

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.

preprint2020arXiv

An improved anisotropic hydrodynamics ansatz

We introduce an improved form for the anisotropic hydrodynamics distribution function which explicitly takes into account the free-streaming and equilibrating contributions separately. We demonstrate that with this improvement one can better reproduce exact results available in the literature for the evolution of moments of the distribution function, in particular, for moments which contain no powers of the longitudinal momentum in their definition (m=0 moments). Using the resulting dynamical equations, we extract the non-equilibrium attractor associated with our improved aHydro ansatz and demonstrate that the improvement also allows one to better reproduce the exact dynamical attractor obtained using kinetic theory in the relaxation time approximation, particularly at early rescaled times and for m=0 moments.

preprint2020arXiv

Non-equilibrium attractor in high-temperature QCD plasmas

We establish the existence of a far-from-equilibrium attractor in weakly-coupled gauge theory undergoing one-dimensional Bjorken expansion. We demonstrate that the resulting far-from-equilibrium evolution is insensitive to certain features of the initial condition, including both the initial momentum-space anisotropy and initial occupancy. We find that this insensitivity extends beyond the energy-momentum tensor to the detailed form of the one-particle distribution function. Based on our results, we assess different procedures for reconstructing the full one-particle distribution function from the energy-momentum tensor along the attractor and discuss implications for the freeze-out procedure used in the phenomenological analysis of ultra-relativistic nuclear collisions.

preprint2020arXiv

Pseudothermalization of the quark-gluon plasma

In this proceedings contribution I review recent work in kinetic theory which demonstrates that, for system undergoing Bjorken expansion, there exists an attractor in all moments of the one-particle distribution function. I discuss how this attractor emerges in both exact solutions obtained in relaxation time approximation (RTA) and the effective kinetic theory approach to high-temperature quantum chromodynamics (QCD). The QCD effective kinetic theory collisional kernel used includes both elastic (2 <-> 2) and LPM-resummed inelastic (2 <-> 1) contributions. The results obtained indicate that a pseudothermal attractor exists in both RTA and QCD kinetic theory and that their respective attractors can be extended to early times when the system is far from equilibrium. Finally, I discuss how knowledge of the QCD effective kinetic theory attractor can be used to assess different hydrodynamic freeze-out prescriptions used in heavy-ion phenomenology. The results obtained show that improved freeze-out prescriptions such as anisotropic hydrodynamics perform better in conditions corresponding to those generated in high-multiplicity pA and pp collisions, e.g. short lifetime and high inverse Reynolds number.

preprint2020arXiv

Two-loop HTL-resummed thermodynamics for N=4 supersymmetric Yang-Mills theory

We compute the two-loop hard-thermal-loop (HTL) resummed thermodynamic potential for N=4 supersymmetric Yang-Mills (SYM). Our final result is manifestly gauge-invariant and was renormalized using only simple vacuum energy, gluon mass, scalar mass, and quark mass counter terms. The HTL mass parameters m_D, M_D, and m_q are then determined self-consistently using a variational prescription which results in a set of coupled gap equations. Based on this, we obtain the two-loop HTL-resummed thermodynamic functions of N=4 SYM. We compare our final result with known results obtained in the weak- and strong-coupling limits. We also compare to previously obtained approximately self-consistent HTL resummations and Padé approximants. We find that the two-loop HTL resummed results for the scaled entropy density is a quantitatively reliable approximation to the scaled entropy density for 0 <= lambda <~ 2 and is in agreement with previous approximately self-consistent HTL resummation results for lambda <~ 6.

preprint2019arXiv

Exact solution for the non-equilibrium attractor in number-conserving relaxation time approximation

We extend previous studies of the conformal 0+1d kinetic non-equilibrium attractor in relaxation time approximation by enforcing number conservation through the introduction of a dynamical fugacity (chemical potential). We derive two coupled integral equations for the effective temperature and fugacity which are then solved numerically to obtain the exact solution. The resulting solutions exhibit convergence to a unique non-equilibrium attractor when the scaled moments of the distribution function are plotted as a function of the rescaled time w = tau/tau_eq. This occurs even though the system is out of chemical equilibrium at late times. In addition, compared to the case where number conservation was not imposed, we find that the moments converge to their respective attractors more quickly, particularly for moments with m=0. Finally, we compare the resulting attractor moments with predictions from different hydrodynamic frameworks.

preprint2019arXiv

Photon production and elliptic flow from momentum-anisotropic quark-gluon plasma

The emission of real photons from a momentum-anisotropic quark-gluon plasma (QGP) is affected by both the collective flow of the radiating medium and the modification of local rest frame emission rate due to the anisotropic momentum distribution of partonic degrees of freedom. In this paper, we first calculate the photon production rate from an ellipsoidally momentum-anisotropic QGP including hard contributions from Compton scattering and quark pair annihilation and soft contribution calculated using the hard thermal loop (HTL) approximation. We introduce a parametrization of the nonequilibrium rate in order to facilitate its further application in yield and flow calculations. We convolve the anisotropic photon rate with the space-time evolution of QGP provided by 3+1d anisotropic hydrodynamics (aHydro) to obtain the yield and the elliptic flow coefficient $v_2$ of photons from QGP generated at Pb-Pb collisions at LHC at 2.76 TeV and Au-Au collisions at RHIC at 200 GeV. We investigate the effects of various parameters on the results. In particular we analyze the sensitivity of results to initial momentum anisotropy.