Search

Thermal electromagnetic radiation is emitted by the quark-gluon plasma (QGP) throughout its space-time evolution, with production rates that depend characteristically on the temperature. We study this temperature using thermal photons and dileptons using the Trajectum heavy ion code, which is constrained by Bayesian analysis. In addition we present the elliptic flow of both the thermal photons and thermal dileptons including systematic uncertainties corresponding to the model parameter uncertainty. We give a comprehensive overview of the resulting effective temperatures $T_{\rm eff}$, obtained from thermal photon transverse momentum and thermal dilepton invariant mass distributions, as well as the dependence of $T_{\rm eff}$ on various selection criteria of these probes. We conclude that the $T_{\rm eff}$ obtained from thermal photons is mostly insensitive to the temperature of the QGP with a value of $T_{\rm eff} \sim$ 250-300 MeV depending on their transverse momentum, almost independent of collision centrality. Thermal dileptons are much better probes of the QGP temperature as they do not suffer from a blue shift as their invariant mass is used, allowing for a more precise constraint of the QGP temperature during different stages of the evolution of the system. By applying selection criteria on the dilepton transverse momentum and the invariant mass we are able to extract fluid temperatures on average times ranging from late emission ($\langle \tau \rangle = 5.6\,$fm$/c$) to very early emissions ($\langle \tau \rangle < 1.0\,$fm$/c$). Furthermore, we show how these selection criteria can be used to map the elliptic flow of the system all throughout its evolution., Comment: 13 pages, 13 figures. The Trajectum code can be found at https://sites.google.com/view/govertnijs/trajectum and plotting routines can be found at http://wilkevanderschee.nl/trajectum

Advances in CAD and CAM have enabled engineers and design teams to digitally design parts with unprecedented ease. Software solutions now come with a range of modules for optimizing designs for performance requirements, generating instructions for manufacturing, and digitally tracking the entire process from design to procurement in the form of product life-cycle management tools. However, existing solutions force design teams and corporations to take a primarily serial approach where manufacturing and procurement decisions are largely contingent on design, rather than being an integral part of the design process. In this work, we propose a new approach to part making where design, manufacturing, and supply chain requirements and resources can be jointly considered and optimized. We present the Generative Manufacturing compiler that accepts as input the following: 1) An engineering part requirements specification that includes quantities such as loads, domain envelope, mass, and compliance, 2) A business part requirements specification that includes production volume, cost, and lead time, 3) Contextual knowledge about the current manufacturing state such as availability of relevant manufacturing equipment, materials, and workforce, both locally and through the supply chain. Based on these factors, the compiler generates and evaluates manufacturing process alternatives and the optimal derivative designs that are implied by each process, and enables a user guided iterative exploration of the design space. As part of our initial implementation of this compiler, we demonstrate the effectiveness of our approach on examples of a cantilever beam problem and a rocket engine mount problem and showcase its utility in creating and selecting optimal solutions according to the requirements and resources.

Accessing the terahertz (THz) spectral domain through surface-enhanced Raman spectroscopy (SERS) is challenging and opens up the study of low-frequency molecular and electronic excitations. Compared to direct THz probing of heterogenous ensembles, the extreme plasmonic confinement of visible light to deep sub-wavelength scales allows the study of hundreds or even single molecules. We show that self-assembled molecular monolayers of a set of simple aromatic thiols confined inside single-particle plasmonic nanocavities can be distinguished by their low-wavenumber spectral peaks below 200 cm-1, after removal of a bosonic inelastic contribution and an exponential background from the spectrum Developing environment-dependent density-functional-theory simulations of the metal-molecule configuration enables the assignment and classification of their THz vibrations as well as the identification of intermolecular coupling effects and of the influence of the gold surface configuration Furthermore, we show dramatically narrower THz SERS spectra from individual molecules at picocavities, which indicates the possibility to study intrinsic vibrational properties beyond inhomogeneous broadening further supporting the key role of local environment.

The System for Measuring Overlap with Gas (SMOG2) at the LHCb detector enables the study of fixed-target ion-ion collisions at relativistic energies ($\sqrt{s_{\rm NN}}\sim100$ GeV in the centre-of-mass). With input from \textit{ab initio} calculations of the structure of $^{16}$O and $^{20}$Ne, we compute 3+1D hydrodynamic predictions for the anisotropic flow of Pb+Ne and Pb+O collisions, to be tested with upcoming LHCb data. This will allow the detailed study of quark-gluon plasma (QGP) formation as well as experimental tests of the predicted nuclear shapes. Elliptic flow ($v_2$) in Pb+Ne collisions is greatly enhanced compared to the Pb+O baseline due to the shape of $^{20}$Ne, which is deformed in a bowling-pin geometry. Owing to the large $^{208}$Pb radius, this effect is seen in a broad centrality range, a unique feature of this collision configuration. Larger elliptic flow further enhances the quadrangular flow ($v_4$) of Pb+Ne collisions via non-linear coupling, and impacts the sign of the kurtosis of the elliptic flow vector distribution ($c_2\{4\}$). Exploiting the shape of $^{20}$Ne proves thus an ideal method to investigate the formation of QGP in fixed-target experiments at LHCb, and demonstrates the power of SMOG2 as a tool to image nuclear ground states.

Whether or not femto-scale droplets of quark-gluon plasma (QGP) are formed in so-called small systems at high-energy colliders is a pressing question in the phenomenology of the strong interaction. For proton-proton or proton-nucleus collisions the answer is inconclusive due to the large theoretical uncertainties plaguing the description of these processes. While upcoming data on collisions of $^{16}$O nuclei may mitigate these uncertainties in the near future, here we demonstrate the unique possibilities offered by complementing $^{16}$O$^{16}$O data with collisions of $^{20}$Ne ions. We couple both NLEFT and PGCM ab initio descriptions of the structure of $^{20}$Ne and $^{16}$O to hydrodynamic simulations of $^{16}$O$^{16}$O and $^{20}$Ne$^{20}$Ne collisions at high energy. We isolate the imprints of the bowling-pin shape of $^{20}$Ne on the collective flow of hadrons, which can be used to perform quantitative tests of the hydrodynamic QGP paradigm. In particular, we predict that the elliptic flow of $^{20}$Ne$^{20}$Ne collisions is enhanced by as much as 1.170(8)$_{\rm stat.}$(30)$_{\rm syst.}$ for NLEFT and 1.139(6)$_{\rm stat.}$(39)$_{\rm syst.}$ for PGCM relative to $^{16}$O$^{16}$O collisions for the 1% most central events. At the same time, theoretical uncertainties largely cancel when studying relative variations of observables between two systems. This demonstrates a method based on experiments with two light-ion species for precision characterizations of the collective dynamics and its emergence in a small system., Comment: 17 pages, 14 figures. The Trajectum code can be found at https://sites.google.com/view/govertnijs/trajectum and plotting routines can be found at http://wilkevanderschee.nl/trajectum

We report on progress in the nonperturbative understanding of quarkonium dynamics inside a thermal plasma. The time evolution of small-size quarkonium is governed by two-point correlation functions of chromoelectric fields dressed with an adjoint Wilson line, known in this context as generalized gluon distributions (GGDs). The GGDs have been calculated in both weakly and strongly coupled plasmas by using perturbative and holographic methods. Strikingly, the results of our calculations for a strongly coupled plasma indicate that the quarkonium dissociation and recombination rates vanish in the transport descriptions that assume quarkonium undergoes Markovian dynamics. However, this does not imply that the dynamics is trivial. As a starting point to explore the phenomenological consequences of the result at strong coupling, we show a calculation of the $\Upsilon(1S)$ formation probability in time-dependent perturbation theory. This is a first step towards the development of a transport formalism that includes non-Markovian effects, which, depending on how close the as of yet undetermined nonperturbative QCD result of the GGDs is to the strongly coupled $\mathcal{N}=4$ SYM result, could very well dominate over the Markovian ones in quark-gluon plasma produced at RHIC and the LHC., Comment: 4 pages, 1 figure. Contribution to the proceedings of the 30th International Conference on Ultra-relativistic Nucleus-Nucleus Collisions (Quark Matter 2023)

Ultracentral heavy ion collisions due to their exceptionally large multiplicity probe an interesting regime of quark-gluon plasma where the size is (mostly) fixed and fluctuations in the initial condition dominate. Spurred by the recent measurement of the CMS collaboration we investigate the driving factors of the increase of transverse momentum, including a complete analysis of the influence of the QCD equation of state. Particularly interesting is the influence of the centrality selection as well as the initial energy deposition., Comment: 4 pages, 5 figures and an appendix. The Trajectum code can be found at https://sites.google.com/view/govertnijs/trajectum and plotting routines can be found at http://wilkevanderschee.nl/trajectum

Modified Newtonian dynamics (MOND) is a promising alternative to dark matter. To further test the theory, there is a need for fluid- and particle-dynamics simulations. The force in MOND is not a direct particle-particle interaction, but derives from a potential for which a nonlinear partial differential equation (PDE) needs to be solved. Normally, this makes the problem of simulating dynamical evolution computationally expensive. We intend to develop a fast particle-mesh (PM) code for MOND (the AQUAL formalism). We transformed the nonlinear equation for MOND into a system of linear PDEs plus one algebraic equation. An iterative scheme with the fast Fourier transform (FFT) produces successively better numerical approximations. The algorithm was tested for dynamical systems in MOND where analytical solutions are known: the two-body problem, a body with a circular ring, and a spherical distribution of particles in thermal equilibrium in the self-consistent potential. The PM code can accurately calculate the forces at subpixel scale and reproduces the analytical solutions. Four iterations are required for the potential, but when the spatial steps are small compared to the kernel width, one iteration is suffices. The use of a smoothing kernel for the accelerations is inevitable in order to eliminate the self-gravity of the point particles. Our PDE solver is $15$ to $42$ times as slow as a standard Poisson solver. However, the smoothing and particle propagation takes up most of the time above one particle per $10^3$ pixels. The FFTs, the smoothing, and the propagation part in the code can all be parallelized., Comment: 13 pages, 13 figures

We study the generalized gluon distribution that governs the dynamics of quarkonium inside a non-Abelian thermal plasma characterizing its dissociation and recombination rates. This gluon distribution can be written in terms of a correlation function of two chromoelectric fields connected by an adjoint Wilson line. We formulate and calculate this object in $\mathcal{N}=4$ supersymmetric Yang-Mills theory at strong coupling using the AdS/CFT correspondence, allowing for a nonzero center-of-mass velocity $v$ of the heavy quark pair relative to the medium. The effect of a moving medium on the dynamics of the heavy quark pair is described by the simple substitution $T \to \sqrt{\gamma} \, T$ in agreement with previous calculations of other observables at strong coupling, where $T$ is the temperature of the plasma in its rest frame, and $\gamma = (1 - v^2)^{-1/2}$ is the Lorentz boost factor. Such a velocity dependence can be important when the quarkonium momentum is larger than its mass. Contrary to general expectations for open quantum systems weakly coupled with large thermal environments, the contributions to the transition rates that are usually thought of as the leading ones in Markovian descriptions vanish in this strongly coupled plasma. This calls for new theoretical developments to assess the effects of strongly coupled non-Abelian plasmas on in-medium quarkonium dynamics. Finally, we compare our results with those from weakly coupled QCD, and find that the QCD result moves toward the $\mathcal{N}=4$ strongly coupled result as the coupling constant is increased within the regime of applicability of perturbation theory. This behavior makes it even more pressing to develop a non-Markovian description of quarkonium in-medium dynamics., Comment: 12 pages, 3 figures; v2: published version

Suppression of open heavy flavors and quarkonia in heavy-ion collisions is among the most informative probes of the quark-gluon plasma. Interpreting the full wealth of data obtained from the collision events requires a precise understanding of the evolution of heavy quarks and quarkonia as they propagate through the nearly thermal and strongly coupled plasma. Systematic theoretical studies of quarkonium time evolution in the QGP in the regime where the temperature of the QGP is much smaller than the inverse of quarkonium size have only been carried out in the past few years. Such calculations require the evaluation of a gauge-invariant correlator of chromoelectric fields dressed with Wilson lines, which is similar to, but different from, the correlation used to define the well-known arXiv:hep-ph/0605199 heavy quark diffusion coefficient. The origin of this difference has been explained in arXiv:2107.03945, arXiv:2205.04477, arXiv:2306.13127. Here we show the results of the calculation of the analogous correlator in strongly coupled $\mathcal{N}=4$ SYM using the AdS/CFT correspondence at a finite temperature arXiv:2304.03298. While it resembles the open heavy quark case, it has some crucial differences that highlight the relevance of quantum color correlations. We will also discuss the results for the quarkonium transport coefficients obtained from this correlator, thereby establishing the first analytic results at strong coupling in this context. We find that they vanish in the strong coupling limit for an $\mathcal{N}=4$ SYM plasma with a large number of colors., Comment: 6 pages, 2 figures. Contribution to the proceedings of the 11th International Conference on Hard and Electromagnetic Probes of High Energy Nuclear Collisions (Hard Probes 2023)

We present an estimate of the yield of hard probes expected for collisions of the isobars $^{96}_{44}$Ru and $^{96}_{40}$Zr at collision energies reachable at RHIC and the LHC\@. These yields are proportional to the number of binary nucleon-nucleon interactions, which is characteristically different due to the presence of the large neutron skin in $^{96}_{40}$Zr. This provides an independent opportunity to measure the difference between the neutron skin of $^{96}_{44}$Ru and $^{96}_{40}$Zr, which can provide an important constraint on the Equation of State of cold neutron-rich matter., Comment: 5 pages, 6 figures and an appendix. Comments are welcome

Emergent bulk properties of matter governed by the strong nuclear force give rise to physical phenomena across vastly different scales, ranging from the shape of atomic nuclei to the masses and radii of neutron stars. They can be accessed on Earth by measuring the spatial extent of the outer skin made of neutrons that characterises the surface of heavy nuclei. The isotope $^{208}$Pb, owing to its simple structure and neutron excess, has been in this context the target of many dedicated efforts. Here, we determine the neutron skin from measurements of particle distributions and their collective flow in $^{208}$Pb+$^{208}$Pb collisions at ultrarelativistic energy performed at the Large Hadron Collider, which are sensitive to the overall size of the colliding $^{208}$Pb ions. By means of state-of-the-art global analysis tools within the hydrodynamic model of heavy-ion collisions, we infer a neutron skin $\Delta r_{np}=0.217\pm0.058$ fm, consistent with nuclear theory predictions, and competitive in accuracy with a recent determination from parity-violating asymmetries in polarised electron scattering. We establish thus a new experimental method to systematically measure neutron distributions in the ground state of atomic nuclei., Comment: 8 pages, 6 figures. The Trajectum code can be found at https://sites.google.com/view/govertnijs/trajectum. Plotting routines can be found at http://wilkevanderschee.nl/trajectum. v2 expanded appendices, matches published version

We present an extended Bayesian analysis using Trajectum where the initial condition can now include binary scaling. For the far-from-equilibrium evolution before hydrodynamics we introduce an interpolation between free streaming and a holographically inspired evolution that exhibits fast hydrodynamization. We find strong evidence that binary scaling is incompatible with experimental data and find evidence that the holographic far-from-equilibrium evolution is preferred. We end with a discussion on several changes and improvements in the Bayesian framework., Comment: 18 pages, 8 figures. The Trajectum code can be found at https://sites.google.com/view/govertnijs/trajectum. Plotting routines can be found at http://wilkevanderschee.nl/trajectum

Previous studies have shown that a gauge-invariant correlation function of two chromoelectric fields connected by a straight timelike adjoint Wilson line encodes crucial information about quark-gluon plasma (QGP) that determines the dynamics of small-sized quarkonium in the medium. Motivated by the successes of holographic calculations to describe strongly coupled QGP, we calculate the analog gauge-invariant correlation function in strongly coupled $\mathcal{N}=4$ supersymmetric Yang-Mills theory at finite temperature by using the AdS/CFT correspondence. Our results indicate that the transition processes between bound and unbound quarkonium states are suppressed in strongly coupled plasmas, and moreover, the leading contributions to these transition processes vanish in both the quantum Brownian motion and quantum optical limits of open quantum system approaches to quarkonia., Comment: 82 pages, 11 figures

Anti-Stokes photoluminescence (PL) is light emission at a higher photon energy than the excitation, with applications in optical cooling, bioimaging, lasing, and quantum optics. Here, we show how plasmonic nano-cavities activate anti-Stokes PL in WSe$_2$ monolayers through resonant excitation of a dark exciton. The tightly confined plasmonic fields excite the out-of-plane transition dipole of the dark exciton, leading to light emission from the bright exciton at higher energy. Through statistical measurements on hundreds of plasmonic cavities, we show that coupling to the dark exciton is key to achieving a near hundred-fold enhancement of the upconverted PL intensity. This is further corroborated by experiments in which the laser excitation wavelength is tuned across the dark exciton. Finally, we show that an asymmetric nanoparticle shape and precise geometry are key for consistent activation of the dark exciton and efficient PL upconversion. Our work introduces a new excitation channel for anti-Stokes PL in WSe$_2$ and paves the way for large-area substrates providing optical cooling, anti-Stokes lasing, and radiative engineering of excitons.

This White Paper presents the community inputs and scientific conclusions from the Hot and Cold QCD Town Meeting that took place September 23-25, 2022 at MIT, as part of the Nuclear Science Advisory Committee (NSAC) 2023 Long Range Planning process. A total of 424 physicists registered for the meeting. The meeting highlighted progress in Quantum Chromodynamics (QCD) nuclear physics since the 2015 LRP (LRP15) and identified key questions and plausible paths to obtaining answers to those questions, defining priorities for our research over the coming decade. In defining the priority of outstanding physics opportunities for the future, both prospects for the short (~ 5 years) and longer term (5-10 years and beyond) are identified together with the facilities, personnel and other resources needed to maximize the discovery potential and maintain United States leadership in QCD physics worldwide. This White Paper is organized as follows: In the Executive Summary, we detail the Recommendations and Initiatives that were presented and discussed at the Town Meeting, and their supporting rationales. Section 2 highlights major progress and accomplishments of the past seven years. It is followed, in Section 3, by an overview of the physics opportunities for the immediate future, and in relation with the next QCD frontier: the EIC. Section 4 provides an overview of the physics motivations and goals associated with the EIC. Section 5 is devoted to the workforce development and support of diversity, equity and inclusion. This is followed by a dedicated section on computing in Section 6. Section 7 describes the national need for nuclear data science and the relevance to QCD research., Comment: QCD Town Meeting White Paper, as submitted to 2023 NSAC LRP committee on Feb. 28, 2023

Algorithms that involve both forecasting and optimization are at the core of solutions to many difficult real-world problems, such as in supply chains (inventory optimization), traffic, and in the transition towards carbon-free energy generation in battery/load/production scheduling in sustainable energy systems. Typically, in these scenarios we want to solve an optimization problem that depends on unknown future values, which therefore need to be forecast. As both forecasting and optimization are difficult problems in their own right, relatively few research has been done in this area. This paper presents the findings of the ``IEEE-CIS Technical Challenge on Predict+Optimize for Renewable Energy Scheduling," held in 2021. We present a comparison and evaluation of the seven highest-ranked solutions in the competition, to provide researchers with a benchmark problem and to establish the state of the art for this benchmark, with the aim to foster and facilitate research in this area. The competition used data from the Monash Microgrid, as well as weather data and energy market data. It then focused on two main challenges: forecasting renewable energy production and demand, and obtaining an optimal schedule for the activities (lectures) and on-site batteries that lead to the lowest cost of energy. The most accurate forecasts were obtained by gradient-boosted tree and random forest models, and optimization was mostly performed using mixed integer linear and quadratic programming. The winning method predicted different scenarios and optimized over all scenarios jointly using a sample average approximation method.

Driven by the global decarbonization effort, the rapid integration of renewable energy into the conventional electricity grid presents new challenges and opportunities for the battery energy storage system (BESS) participating in the energy market. Energy arbitrage can be a significant source of revenue for the BESS due to the increasing price volatility in the spot market caused by the mismatch between renewable generation and electricity demand. In addition, the Frequency Control Ancillary Services (FCAS) markets established to stabilize the grid can offer higher returns for the BESS due to their capability to respond within milliseconds. Therefore, it is crucial for the BESS to carefully decide how much capacity to assign to each market to maximize the total profit under uncertain market conditions. This paper formulates the bidding problem of the BESS as a Markov Decision Process, which enables the BESS to participate in both the spot market and the FCAS market to maximize profit. Then, Proximal Policy Optimization, a model-free deep reinforcement learning algorithm, is employed to learn the optimal bidding strategy from the dynamic environment of the energy market under a continuous bidding scale. The proposed model is trained and validated using real-world historical data of the Australian National Electricity Market. The results demonstrate that our developed joint bidding strategy in both markets is significantly profitable compared to individual markets.

The rapid adoption of residential solar photovoltaics (PV) has resulted in regular overvoltage events, due to correlated reverse power flows. Currently, PV inverters prevent damage to electronics by curtailing energy production in response to overvoltage. However, this disproportionately affects households at the far end of the feeder, leading to an unfair allocation of the potential value of energy produced. Globally optimizing for fair curtailment requires accurate feeder parameters, which are often unknown. This paper investigates reinforcement learning, which gradually optimizes a fair PV curtailment strategy by interacting with the system. We evaluate six fairness metrics on how well they can be learned compared to an optimal solution oracle. We show that all definitions permit efficient learning, suggesting that reinforcement learning is a promising approach to achieving both safe and fair PV coordination.

Plasmonic nano-antennas can focus light to nanometre length-scales providing intense field enhancements. For the tightest optical confinements (0.5-5 nm) achieved in plasmonic gaps, the gap spacing, refractive index, and facet width play a dominant role in determining the optical properties making tuning through antenna shape challenging. We show here that controlling the surrounding refractive index instead allows both efficient frequency tuning and enhanced in/output-coupling through retardation matching as this allows dark modes to become optically active, improving widespread functionalities.

Molecular vibrations and their dynamics are of outstanding importance for electronic and thermal transport in nanoscale devices as well as for molecular catalysis. The vibrational dynamics of <100 molecules are studied through three-colour time-resolved coherent anti-Stokes Raman spectroscopy (trCARS) using plasmonic nanoantennas. This isolates molecular signals from four-wave mixing (FWM), while using exceptionally low nanowatt powers to avoid molecular damage via single-photon lock-in detection. FWM is found to be strongly suppressed in nm-wide plasmonic gaps compared to plasmonic nanoparticles. The ultrafast vibrational decay rates of biphenyl-4-thiol molecules are accelerated ten-fold by a transient rise in local non-equilibrium temperature excited by enhanced, pulsed optical fields within these plasmonic nanocavities. Separating the contributions of vibrational population decay and dephasing carefully explores the vibrational decay channels of these tightly confined molecules. Such extreme plasmonic enhancement within nanogaps opens up prospects for measuring single-molecule vibrationally-coupled dynamics and diverse molecular optomechanics phenomena.