Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Particle correlations & fluctuations”

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 307 records · Page 17

An Analytical Model for the Prediction of a Micro-Dosimeter Response Function

A rapid analytical procedure for the prediction of a micro-dosimeter response function in low Earth orbit (LEO), correlated with the Space Transportation System (STS, shuttle) Tissue Equivalent Proportional Counter (TEPC) measurements is presented. The analytical model takes into consideration the energy loss straggling and chord length distribution of the detector, and is capable of predicting energy deposition fluctuations in a cylindrical micro-volume of arbitrary aspect ratio (height/diameter) by incoming ions through both direct and indirect (ray) events. At any designated (ray traced) target point within the vehicle, the model accepts the differential flux spectrum of Galactic Cosmic Rays (GCR) and/or trapped protons at LEO as input. On a desktop PC, the response function of TEPC for each ion in the GCR/trapped field is computed at the average rate of 30 seconds/ion. The ionizing radiation environment at LEO is represented by O'Neill fs GCR model (2004), covering charged particles in the 1 less than or equal to Z less than or equal to 28. O'Neill's free space GCR model is coupled with the Langley Research Center (LaRC) angular dependent geomagnetic cutoff model to compute the transmission coefficient in LEO. The trapped proton environment is represented by a LaRC developed time dependent procedure which couples the AP8MIN/AP8MAX, Deep River Neutron Monitor (DRNM) and F10.7 solar radio frequency measurements. The albedo neutron environment is represented by the extrapolation of the Atmospheric Ionizing Radiation (AIR) measurements. The charged particle transport calculations correlated with STS 51 and 114 flights are accomplished by using the most recent version (2005) of the LaRC deterministic High charge (Z) and Energy TRaNsport (HZETRN) code. We present the correlations between the TEPC model predictions (response function) and TEPC measured differential/integral spectra in the lineal energy (y) domain for both GCR and trapped protons, with the conclusion that the model correctly accounts for the increase in flux at low y values where energetic ions are the primary contributor. We further discuss that, even with the incorporation of angular dependency in the cutoffs, comparison of the GCR differential/integral flux between STS 51 and 114 TEPC measured data and current calculations indicates that there still exists an underestimation by the simulations at low to mid range y values. This underestimation is partly related the exclusion of the secondary pion particle production from the current version of HZETRN.

Badavi, Francis F.↗

Viscoelastic relaxation and topological fluctuations in glass-forming liquids

In this study, a method for characterizing the topological fluctuations in liquids is proposed. This approach exploits the concept of the weighted gyration tensor of a collection of particles and permits the definition of a local configurational unit (LCU). The first principal axis of the gyration tensor serves as the director of the LCU, which can be tracked and analyzed by molecular dynamics simulations. Analysis of moderately supercooled Kob–Andersen mixtures suggests that orientational relaxation of the LCU closely follows viscoelastic relaxation and exhibits a two-stage behavior. The slow relaxing component of the LCU corresponds to the structural, Maxwellian mechanical relaxation. Additionally, it is found that the mean curvature of the LCUs is approximately zero at the Maxwell relaxation time with the Gaussian curvature being negative. This observation implies that structural relaxation occurs when the configurationally stable and destabilized regions interpenetrate each other in a bicontinuous manner. Finally, the mean and Gaussian curvatures of the LCUs can serve as reduced variables for the shear stress correlation, providing a compelling proof of the close connection between viscoelastic relaxation and topological fluctuations in glass-forming liquids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Relaxation and Aging of Nanosphere Assemblies at a Water–Oil Interface

The relaxation and aging of an assembly of spherical nanoparticles (NPs) at a water-oil interface are characterized in situ by grazing incidence X-ray photon correlation spectroscopy. The dynamics of the interfacial assembly is measured while the interface saturates with NPs. Weak attractions between NPs lead to gel-like structures in the assembly, where the in-plane ordering is inhibited by the broad size distribution of the NPs. Structural rearrangements on the length scale of the NP-NP center-to-center distances proceed by intermittent fluctuations instead of continuous cooperative motions. Here, the coexistence of rapid and slow NP populations is confirmed, as commonly observed in soft glass-forming materials. Dynamics are increasingly slowed as the NPs initially segregate to the locally clustered interface. The structural relaxation of the NPs in these localized clusters is 5 orders of magnitude slower than that of free particles in the bulk. When the interface is nearly saturated, the time for relaxation increases suddenly due to the onset of local jamming, and the dynamics slow exponentially afterward until the system reaches collective jamming by cooperative rearrangements. This investigation provides insights into structural relaxations near the glass transition and the evolution of the structure and dynamics of the assemblies as they transition from an isotropic liquid to a dense disordered film.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Topologically Protected Flatness in Chiral Moiré Heterostructures

The observation of delicate correlated phases in twisted heterostructures of graphene and transition metal dichalcogenides suggests that moiré flat bands are intrinsically resilient against certain types of disorder. Here, we investigate the robustness of moiré flat bands in the chiral limit of the Bistritzer-MacDonald model—applicable to both platforms in certain limits—and demonstrate drastic differences between the first magic angle and higher magic angles in response to chiral symmetric disorder that arise, for instance, from lattice relaxation. We understand these differences using a hidden constant of motion that permits the decomposition of the non-Abelian gauge field induced by interlayer tunnelings into two decoupled Abelian ones. At all magic angles, the resulting effective magnetic field splits into an anomalous contribution and a fluctuating part. The anomalous field maps the moiré flat bands onto a zeroth Dirac Landau level, whose flatness withstands any chiral symmetric perturbation such as nonuniform magnetic fields due to a topological index theorem—thereby underscoring a topological mechanism for band flatness. Only the first magic angle can fully harness this topological protection due to its weak fluctuating magnetic field. In higher magic angles, the amplitude of fluctuations largely exceeds the anomalous contribution, which we find results in a physically meaningless chiral operator and an extremely large sensitivity to microscopic details and an exponential collapse of the single-particle gap. Through numerical simulations, we further study various types of disorder and identify the scattering processes that are enhanced or suppressed in the chiral limit. Interestingly, we find that the topological suppression of disorder broadening persists away from the chiral limit and is further accentuated by isolating a single sublattice polarized flat band in energy. Our analysis suggests the Berry curvature hot spot at the top of the K and K ′ valence band in the transition metal dichalcogenide monolayers is essential for the stability of its moiré flat bands and their correlated states. Published by the American Physical Society 2025

Crépel, Valentin (ORCID:0000000302403412)↗

Kinetic theory of turbulent compressible flows and comparison with classical theory

A complete set of turbulent correlations is given in terms of expansion coefficients in a double series of Hermite polynomials of a two-particle correlation in phase space. Only two of these coefficients, corresponding to Reynolds stress and turbulent heat flux, are shown to appear in gasdynamic equations of turbulence. This is accomplished by incorporating the phase-space correlation in the collision integral of Boltzmann's equation and by deriving generalized Navier-Stokes and Fourier transport relations. In contrast to the classical formalism, the simple expression remains invariant in form whether the flow is compressible or not. Nevertheless, a bilinear transformation of fluctuating quantities shows that the two formalisms are identical with regard to terms of double correlations. The kinetic theory justifies nonexistence of the higher-order correlations characteristic of compressible-turbulence equations in the classical regime.

Tsuge, S.↗

Fluid Physics in a Fluctuating Acceleration Environment

We summarize several aspects of an ongoing investigation of the effects that stochastic residual accelerations (g-jitter) onboard spacecraft can have on experiments conducted in a microgravity environment. The residual acceleration field is modeled as a narrow band noise, characterized by three independent parameters: intensity (g(exp 2)), dominant angular frequency Omega, and characteristic correlation time tau. Realistic values for these parameters are obtained from an analysis of acceleration data corresponding to the SL-J mission, as recorded by the SAMS instruments. We then use the model to address the random motion of a solid particle suspended in an incompressible fluid subjected to such random accelerations. As an extension, the effect of jitter on coarsening of a solid-liquid mixture is briefly discussed, and corrections to diffusion controlled coarsening evaluated. We conclude that jitter will not be significant in the experiment 'Coarsening of solid-liquid mixtures' to be conducted in microgravity. Finally, modifications to the location of onset of instability in systems driven by a random force are discussed by extending the standard reduction to the center manifold to the stochastic case. Results pertaining to time-modulated oscillatory convection are briefly discussed.

Thomson, J. Ross↗

Polyhexatic and polycrystalline states of skyrmion lattices

We report Monte Carlo studies of lattices of up to 10 5 skyrmions treated as particles with negative core energy and repulsive interaction obtained from a microscopic spin model. Temperature dependence of translational and orientational correlations has been investigated for different experimental protocols and initial conditions. Cooling the skyrmion liquid from a fully disordered high-temperature state results in the formation of a skyrmion polycrystal. A perfect skyrmion lattice prepared at T = 0, on raising temperature undergoes a first-order melting transition into a polyhexatic state that consists of large orientationally ordered domains of fluctuating shape. On the further increasing temperature, these domains decrease in size, leading to a fully disordered liquid of skyrmions.

36 MATERIALS SCIENCE↗

Production of neutron-rich heavy nuclei in deep-inelastic 208 Pb + 208 Pb collisions within the stochastic mean-field theory

In deep-inelastic collisions of heavy nuclei, reaction products with a wide range of mass and charge are produced. Such collisions been considered as a possible way to produce superheavy nuclei, as an alternative to fusion reactions. To provide reliable theoretical predictions, it is desired to develop microscopic approaches that correctly and accurately describe nucleon transfer processes in dissipative collisions of heavy nuclei. The purpose of the present work is (1) to investigate the mechanism of nucleon transfers in dissipative collisions of two heavy nuclei, and (2) to explore possible pathways to produce neutron-rich heavy nuclei, through detailed theoretical analyses of fluctuations and correlations in nucleon transfers in 208 Pb + 208 Pb reactions. Three-dimensional time-dependent Hartree-Fock (TDHF) calculations are performed for the collisions of 208 Pb + 208 Pb at 𝐸 c.m. = 832, 936, and 1040 MeV, using the Skyrme SLy4d energy density functional. To calculate fluctuations and correlations in nucleon transfers, we employ the stochastic mean-field (SMF) theory, and the results are compared with another theoretical framework currently available, the time-dependent random phase approximation (TDRPA). Primary and secondary production cross sections are calculated with the SMF theory combined with a statistical model, GEMINI ++ . Using information of nucleon flow across a neck of colliding nuclei in TDHF calculations, we solve quantal diffusion equations for fluctuations and correlations in nucleon transfers based on the SMF theory. From the SMF calculations, we obtain the time evolution of diffusion coefficients as well as fluctuations and correlations in nucleon transfers for a range of initial orbital angular momenta. We compare the results of the SMF calculations with those of TDRPA, showing that TDRPA tends to predict substantially larger fluctuations and correlations in strongly damped collisions of heavy nuclei, which exhibit complex initial angular momentum dependence, while the SMF results provide almost constant (stable) values. Using the obtained fluctuations and correlations, we calculate primary and secondary production cross sections for the 208 Pb + 208 Pb collisions. From the results, we find that both lighter and heavier reaction products as compared to 208 Pb are produced for a wide region in the 𝑁−𝑍 plane as primary products, thanks to the quantal diffusion mechanism in the dissipative collisions. However, we show that cross sections for production of heavy nuclei with 𝑍 ≳ 90 or 𝑁 ≳ 135 are washed out due to secondary particle evaporation and/or fission processes. On the other hand, we find that there remain sizable cross sections for production of neutron-rich nuclei along 𝑁 = 126 with 𝑍< 82, even after secondary disintegration processes. We demonstrate that the secondary production cross sections depend weakly on incident energies, but lower (higher) energy is slightly preferred for production of nuclei with smaller (larger) atomic numbers as compared to 𝑍 = 82. Here, based on the microscopic SMF calculations, it has been shown that deep-inelastic collisions of heavy nuclei, such as 208 Pb + 208 Pb examined in this study, can be a promising means to produce neutron-rich heavy nuclei along 𝑁=126. Discrepancies between the SMF and TDRPA approaches are left unsolved for future investigations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Global gyro-kinetic ion temperature gradient and trapped electron mode turbulence modelling in X -point geometry in negative and positive triangularity

Comparative modelling of ion temperature gradient (ITG)/trapped electron mode (TEM) turbulence in negative (NT) and positive (PT) triangularity plasma shapes was done using the nonlinear global full-f gyrokinetic particle code JOREK-GK in the realistic X-point tokamak geometry including the Scrape Off Layer (SOL) for TCV and DIII-D parameters. A comparison of JOREK-GK code with the gyrokinetic codes GS2 and GENE was done using NT/PT triangularity TCV L-modes parameters showing good agreement between codes in linear growth rates and clear beneficial effect of NT as compared to PT. Global non-linear modelling of the ITG/TEM saturated turbulence for realistic DIII-D NT pulses was done and compared with numerically constructed PT equilibrium with the same plasma profiles. Existence of longer correlation length of density fluctuations in PT compared to NT was demonstrated. Stronger and more sheared zonal flows are generated via Reynolds stress in NT compared to PT. These factors are stabilizing for TEM/ITG turbulence in NT and lead to smaller heat fluxes and heat conductivities in NT compared to PT configuration. Weak dependence of plasma confinement on collisionality and plasma rotation was found in modelling of DIII-D NT shots similar to the experiment. The confinement scaling with normalized ion gyro-radius ρ* was estimated both for NT and PT. Bohm-like scaling was obtained in both configurations, however with better confinement for NT compared to PT which could be favourable factor for reactor size machines.

ITG↗

Concentration and Velocity Gradients in Fluidized Beds

In this work we focus on the height dependence of particle concentration, average velocity components, fluctuations in these velocities and, with the flow turned off, the sedimentation velocity. The latter quantities are measured using Particle Imaging Velocimetry (PIV). The PIV technique uses a 1-megapixel camera to capture two time-displaced images of particles in the bed. The depth of field of the imaging system is approximately 0.5 cm. The camera images a region with characteristic length of 2.6 cm for the small particles and 4.7 cm. for the large particles. The local direction of particle flow is determined by calculating the correlation function for sub-regions of 32 x 32 pixels. The velocity vector map is created from this correlation function using the time between images (we use 15 to 30 ms). The software is sensitive variations of 1/64th of a pixel. We produce velocity maps at various heights, each consisting of 3844 velocities. We break this map into three vertical zones for increased height information. The concentration profile is measured using an expanded (1 cm diameter) linearly polarized HeNe Laser incident on the fluidized bed. A COHU camera (gamma=1, AGC off) with a lens and a polarizer images the transmitted linearly polarized light to minimize the effects of multiply scattered light. The intensity profile (640 X 480 pixels) is well described by a Gaussian fit and the height of the Gaussian is used to characterize the concentration. This value is compared to the heights found for known concentrations. The sedimentation velocity is estimated using by imaging a region near the bottom of the bed and using PIV to measure the velocity as a function of time. With a nearly uniform concentration profile, the time can be converted to height information. The stable fluidized beds are made from large pseudo-monodisperse particles (silica spheres with radii (250-300) microns and (425-500) microns) dispersed in a glycerin/water mix. The Peclet number is sufficiently large that Brownian motion of the particles can be ignored and the Reynolds number sufficiently small that particle inertia is negligible. A packed particle bed is used to randomize and disperse the flowing fluid introduced by a peristaltic pump. The bed itself is a rectangular glass cell 8 cm wide (x), 0.8 cm deep and a height of 30.5 cm (z). The depth of field of the camera is approximately 0.5 cm so depth information is averaged. Over flow fluid is returned to the reservoir making a closed loop system. In these experiments the particles form a sediment approximately 5.7 cm high with the pump off and expand to 22 cm with the pump on. For the smaller particles the pump velocity is .5 millimeters per second and 1.1 millimeters per second for the large particles. At this concentration the bed has a very well defined top where particle concentration rapidly drops to zero.

McClymer, James P.↗

Fundamental processes in the expansion, energization, and coupling of single- and multi-Ion plasmas in space: Laboratory simulation experiments

We have conducted a laboratory investigation into the physics of plasma expansions and their associated energization processes. We studied single- and multi-ion plasma processes in self-expansions, and included light and heavy ions and heavy/light mixtures to encompass the phenomenological regimes of the solar and polar winds and the AMPTE and CRRES chemical release programs. The laboratory experiments provided spatially-distributed time-dependent measurements of total plasma density, temperature, and density fluctuation power spectra with the data confirming the long-theorized electron energization process in an expanding cloud - a result that was impossible to determine in spaceborne experiments (as e.g., in the CRRES program). These results provided the missing link in previous laboratory and spaceborne programs. confirming important elements in our understanding of such solar-terrestrial processes as manifested in expanding plasmas in the solar wind (e.g., CMES) and in ionospheric outflow in plasmaspheric fluctuate refilling after a storm. The energization signatures were seen in an entire series of runs that varied the ion species (Ar', Xe', Kr' and Ne'), and correlative studies included spectral analyses of electrostatic waves collocated with the energized electron distributions. In all cases wave energies were most intense during the times in which the suprathermal populations were present, with wave intensity increasing with the intensity of the suprathermal electron population. This is consistent with theoretical expectations wherein the energization process is directly attributable to wave particle interactions. No resonance conditions were observed, in an overall framework in which the general wave characteristics were broadband with power decreasing with increasing frequency.

Szuszczewicz, E. P.↗

Flow field characteristics of a 127-pin rod bundle with hexagonal spacer grids

A determination of nominal flow phenomena in liquid metal fast reactor (LMFR) fuel assemblies is critical toward generation-IV reactor development. Axially positioned spacer grids are used to maintain the geometry of hexagonal rod bundles and simultaneously introduce perturbations in the flow. Three-dimensional (3D) printed asymmetric honeycomb spacer grids were installed in a prototypical 127-pin LMFR fuel assembly model to study complex fluid dynamics interactions induced by the spacer grid and rods. To characterize flow dynamics in this intricate geometry, time-resolved particle image velocimetry (TR-PIV) using the matched-index-of-refraction method was employed to obtain non-intrusive velocity measurements for three axial planes (one near-wall and two interior planes) at a Reynolds number of 6000. The statistical TR-PIV results compared sub-channel-dependent normalized time-averaged velocity, velocity fluctuations, Reynolds stress, vorticity, and turbulence kinetic energy distributions. TR-PIV line profiles characterized downstream spacer grid flow dynamics. Two-point spatial and spatial–temporal cross-correlation fields revealed local coherent structures and quantified convection velocities of traveling vortices. Spatial–temporal decomposition using dynamic mode decomposition (DMD) applied to the near-wall vorticity fields extracted turbulent structures and flow instabilities in the wake region of the spacer grid, along with their decay and frequency rates. Reduced-order velocity fields from DMD reconstructions identified the most energy-containing coherent structures persistent in the near-wall region. This research provides experimental data sets and analyses of flow behavior in rod bundles with hexagonal spacer grids. In conclusion, the results are critical toward LMFR design and geometry optimization, crucial for the validation of computational fluid dynamics and reduced-order flow models.

42 ENGINEERING↗

Discrete ion stochastic continuum overdamped solvent algorithm for modeling electrolytes

In this paper we develop a methodology for the mesoscale simulation of strong electrolytes. The methodology is an extension of the fluctuating immersed-boundary approach that treats a solute as discrete Lagrangian particles that interact with Eulerian hydrodynamic and electrostatic fields. In both algorithms the immersed-boundary method of Peskin is used for particle-field coupling. Hydrodynamic interactions are taken to be overdamped, with thermal noise incorporated using the fluctuating Stokes equation, including a "dry diffusion" Brownian motion to account for scales not resolved by the coarse-grained model of the solvent. Long-range electrostatic interactions are computed by solving the Poisson equation, with short-range corrections included using an immersed-boundary variant of the classical particle-particle particle-mesh technique. Also included is a short-range repulsive force based on the Weeks-Chandler-Andersen potential. This methodology is validated by comparison to Debye-Hückel theory for ion-ion pair correlation functions, and Debye-Hückel-Onsager theory for conductivity, including the Wien effect for strong electric fields. In each case, good agreement is observed, provided that hydrodynamic interactions at the typical ion-ion separation are resolved by the fluid grid.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Pressure-based process monitoring of direct-ink write material extrusion additive manufacturing

As additive manufacturing (AM) has become a reliable method for creating complex and unique hardware rapidly, the quality assurance of printed parts remains a priority. In situ process monitoring offers an approach for performing quality control while simultaneously minimizing post-production inspection. For extrusion printing processes, direct linkages between extrusion pressure fluctuations and print defects can be established by integrating pressure sensors onto the print head. In this work, the sensitivity of process monitoring is tested using engineered spherical defects. Pressure and force sensors located near an ink reservoir and just before the nozzle are shown to assist in identification of air bubbles, changes in height between the print head and build surface, clogs, and particle aggregates with a detection threshold of 60–70% of the nozzle diameter. Visual evidence of printed bead distortion is quantified using optical image analysis and correlated to pressure measurements. Importantly, this methodology provides an ability to monitor the quality of AM parts produced by extrusion printing methods and can be accomplished using commonly available pressure-sensing equipment.

36 MATERIALS SCIENCE↗

Equilibrium expectations for non-Gaussian fluctuations near a QCD critical point

With the highly anticipated results from the Beam Energy Scan II program at RHIC being recently revealed, an understanding of particle-number fluctuations and their significance as a potential signature of a possible QCD critical point is crucial. Early works that embarked on this endeavor sought to estimate the fluctuations due to the presence of a critical point assuming they stay in equilibrium. From these results came the proposal to focus efforts on higher, non-Gaussian, moments of the event-by-event distributions, in particular of the number of protons. These non-Gaussian moments are especially sensitive to critical fluctuations, as their magnitudes are proportional to high powers of the critical correlation length. As the equation of state provides key input for hydrodynamical simulations of heavy-ion collisions, we estimate equilibrium fluctuations from the BEST equation of state (EoS) that includes critical features from the 3D Ising Model. In particular, the proton factorial cumulants and their dependence on non-universal mapping parameters is investigated within the BEST EoS. Furthermore, the correlation length, as a central quantity for the assessment of fluctuations in the vicinity of a critical point, is also calculated in a consistent manner with the scaling equation of state. An understanding of the equilibrium estimates of proton factorial cumulants will be useful for further comparison to estimates of out-of-equilibrium fluctuations in order to determine the magnitude of the observable fluctuations to be expected in heavyion collision experiments, in which the time spent near a critical point is short.

Karthein, Jamie M. [Massachusetts Institute of Tec↗

The velocity correlation function in cosmic-ray diffusion theory

It is shown that Earl's (1973) eigenvalue sum for the cosmic-ray spatial diffusion coefficient parallel to the mean magnetic field is precisely equivalent to the time integral of the particle-velocity correlation function parallel to the mean field. A derivation due to Kubo (1957) is applied to cosmic-ray pitch-angle scattering, and it is proven that all nine components of the cosmic-ray diffusion tensor can be expressed as integrals over the velocity correlation function. A pitch-angle correlation function is derived, and the effect of long-wavelength turbulence on the velocity correlation function and spatial diffusion coefficients is examined. Application of the velocity-correlation method to a realistic case involving both pitch-angle scattering and appreciable fluctuation in the direction of the local field indicates that long-wavelength turbulence in the local field reduces the parallel diffusion coefficient and places an upper limit on the ratio of the perpendicular to parallel diffusion coefficients.

Forman, M. A.↗

Beyond AdS 2 /dCFT 1 : insertions in two Wilson loops

We consider two-point correlators of local operator insertions in a system of two Wilson-Maldacena loops in $\mathcal{N}$ = 4 supersymmetric Yang-Mills theory on both sides of the AdS/CFT correspondence. On the holographic side the correlator of two Wilson-Maldacena loops is given by a classical string world-sheet which in one phase connects two asymptotically AdS 2 regions and in the other phase is given by two disconnected AdS 2 caps; this configuration breaks supersymmetry as well as conformal invariance. We present a complete systematic account of the string world-sheet fluctuations, including the fermionic sector, and study the behavior of the holographic two-point correlators. On the field theory side we compute certain two-point correlators of local operator insertions by resumming sets of ladder diagrams. Our results demonstrate the efficacy of previously developed methods in tackling this non-conformal, non-susy regime.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Accessing subnuclear fluctuations and saturation with multiplicity dependent J/ψ production in p+p and p+Pb collisions

We study the production of J/ψ vector mesons as a function of charged hadron multiplicity in p+p and p+Pb collisions at LHC energies. We employ the color glass condensate framework, using running coupling Balitsky-Kovchegov evolved dipole amplitudes, to compute gluon and cc¯-pair production. We use fragmentation functions to obtain charged hadrons, and explore two different hadronization schemes for the J/ψ: non-relativistic quantum chromodynamics and the improved color evaporation model. In our framework, event-by-event multiplicity fluctuations of both hadrons and J/ψ are driven by geometric and saturation scale normalization fluctuations. Studying the correlation between J/ψ and hadron multiplicity, we show that the characteristic difference between forward and backward rapidity in p+Pb collisions is a result of different degrees of saturation probed at different rapidities. We demonstrate that experimental data on heavy-flavor production as a function of event activity provide stringent constraints on the fluctuating proton structure

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗