Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic timing analysis”

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 1,045 records · Page 58

Modeling the dynamics of phospholipids in the fluid phase of liposomes

Here, we present the derivation of a new model to describe neutron spin echo spectroscopy and quasi-elastic neutron scattering data on liposomes. We compare the new model with existing approaches and benchmark it with experimental data. The analysis indicates the importance of including all major contributions in the modeling of the intermediate scattering function. Simultaneous analysis of the experimental data on lipids with full contrast and tail contrast matched samples reveals highly confined lipid tail motion. A comparison of their dynamics demonstrates the statistical independence of tail-motion and height–height correlation of the membrane. A more detailed analysis indicates that the lipid tail relaxation is confined to a potential with cylindrical symmetry, in addition to the undulation and diffusive motion of the liposome. Despite substantial differences in the chemistry of the fatty acid tails, the observation indicates a universal behavior. The analysis of partially deuterated systems confirms the strong contribution of the lipid tail to the intermediate scattering function. Within the time range from 5 to 100 ns, the intermediate scattering function can be described by the height–height correlation function. The existence of the fast-localized tail motion and the contribution of slow translational diffusion of liposomes determine the intermediate scattering function for t < 5 ns and t > 100 ns, respectively. Taking into account the limited time window lowers the bending moduli by a factor of 1.3 (DOPC) to 2 (DMPC) compared to the full range.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spatially varying colloidal phase behavior on multi-dimensional energy landscapes

A method is reported to determine equilibrium concentration profiles and local phase behavior for colloids on multi-dimensional energy landscapes. A general expression is derived based on local particle concentration and osmotic pressure differences that are balanced by forces on colloids due to energy landscape gradients. This analysis is applied to colloidal particles in high frequency AC electric fields within octupolar electrodes, where the energy landscape can be shaped in two dimensions. These results are also directly applicable to any particles having induced dipoles in spatially non-uniform electromagnetic fields. Predictions based on modelling colloids with an effective hard disk equation of state indicate inhomogeneous solid and fluid states coexisting on different shaped energy landscapes including multiple minima. Model predictions show excellent agreement with time-averaged Brownian Dynamic simulations at equilibrium. Findings demonstrate a general approach to understand colloidal phase behavior on energy landscapes due to external fields, which could enable control of colloidal microstructure on morphing energy landscapes and the inverse design of fields to assemble hierarchically structured colloidal materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electrostatic upper-hybrid mode instability driven by a ring electron distribution

Quasi electrostatic fluctuations in the upper-hybrid frequency range are commonly detected in the planetary magnetospheric environment. The origin of such phenomena may relate to the instability driven by a loss-cone feature associated with the electrons populating the dipole-like magnetic field. The present paper carries out a one-dimensional electrostatic particle-in-cell simulation accompanied by a reduced quasilinear kinetic theoretical analysis to investigate the dynamics of the upper-hybrid mode instability driven by an initial ring electron distribution function, which is a form of loss-cone distribution. A favorable comparison is found between the two approaches, which shows that the reduced quasilinear theory, which is grounded in the concept of a model of the particle distribution function that is assumed to maintain a fixed mathematical form except that the macroscopic parameters that define the distribution are allowed to evolve in time, can be an effective tool in the study of plasma instabilities, especially if it is guided by and validated against the more rigorous simulation result.

Physics↗

The effect of local chemical ordering on Frank-Read source activation in a refractory multi-principal element alloy

In this work, we investigate the operation of Frank-Read (FR) sources in a refractory multi-principal element alloy (MPEA). Simulations of discrete dislocation motion in MPEAs is enabled by the development of a phase field dislocation dynamics model that treats the atomic-scale fluctuations in lattice energies across the glide plane, present due to local ordering in the chemical composition within the nominally random MPEA atomic structure. We consider, through simulation, a range of length scales over which ordering occurs, varying from short-range lengths, a few times dislocation core width, to long-range lengths, an order of magnitude longer than the core. Characteristic of this body-centered cubic MPEA, the simulations also include screw/edge character dependence in glide resistance, as informed by atomic scale simulation. The critical stresses to activate the source for the same source size are found to be statistically distributed, as a direct consequence of the underlying variation in lattice energy. Furthermore, analysis of critical state for activating edge and screw FR sources in the MPEA reveals that FR source operation occurs via a two-step mechanism, involving athermal kink-pair formation, unlike the conventional FR source operation in a material with no composition fluctuations. This mechanism lowers the average critical stress required to activate the FR source and causes the statistical dispersion in critical stress to depend on the range of composition ordering. More importantly, it leads to a more severe dependence of source strength on FR source length than predicted by line tension alone.

36 MATERIALS SCIENCE↗

Heat load measurements for the PIP-II pHB650 cryomodule

Phase-3 testing of the pHB650 cryomodule at the PIP-II Injector Test Facility was conducted to evaluate the effectiveness of heat load mitigations performed after earlier phases of testing and to continue pinpointing any sources of unexpectedly high heat loads.. The programme measured HTTS, LTTS, and 2 K isothermal/non-isothermal loads under "standard", "linac", and "simulated dynamic" operating modes, recording data both inside the cryomodule and across the bayonet can circuits. Thermal-acoustic oscillations were eliminated by replacing the original G10 cooldown-valve stem with a stainless-steel stem fitted with wipers. A newly developed Python script automated acquisition of ACNET data, performed real-time heat-load calculations, and generated plots and tables that were posted to the electronic logbook within minutes, vastly reducing manual effort and accelerating feedback between SRF and cryogenics teams. Analysis showed that JT heat-exchanger effectiveness and temperature stratification in the two-phase and relief piping strongly influence the observed loads and helped isolate sources of excess heat. The campaign demonstrates that rigorous pre-test planning, real-time diagnostics, and automated reporting can improve both accuracy and efficiency, providing a template for future PIP-II cryomodule tests and for implementing targeted heat-load mitigations.

Porwisiak, D. [Fermilab; Wroclaw Tech. U.] (ORCID:↗

JIRIAF: JLAB Integrated Research Infrastructure Acros Facilities

The JIRIAF project aims to combine geographically diverse computing facilities into an integrated science infrastructure. This project starts by dynamically evaluating temporarily unallocated or idled compute resources from multiple providers. These resources are integrated to handle additional workloads without affecting local running jobs. This paper describes our approach to launch best-effort batch tasks which exploit these underutilized resources. Our system measures the real-time behavior of jobs running on a machine and learns to distinguish typical performance from outliers. Unsupervised ML techniques are used to analyze hardware-level performance measures, followed by a real-time cross-correlation analysis to determine which applications cause performance degradation. We then ameliorate bad behavior by throttling these processes. We demonstrate that problematic performance interference can be detected and acted on, which makes it possible to continue to share resources between applications and simultaneously maintain high utilization levels in a computing cluster. We relocate the CLAS12 data processing workflow to a remote data center for a case study, preventing file migration and temporal data persistency.

Lawrence, David↗

Glovebox CFD Analysis: Low-Flow, Laminar CFD Analysis of Glovebox Ventilation and Humidity Transport in Ansys Fluent 2025 R1

There are gloveboxes that frequently exceed the 30% relative humidity level required for operations. These exceedances result in alternative operations needing to be utilized. This report documents a computational fluid dynamics (CFD) analysis performed to characterize airflow and humidity transport within a representative glovebox under nominal, low-flow conditions. The analysis focuses on establishing the flow regime, identifying transport-limited regions, and evaluating the implications of laminar flow on humidity removal. A single-phase, laminar flow model with a species transport model was used to represent air and water vapor mixing. Results demonstrate that, at the nominal flow rates, airflow is laminar with weak mixing, which will result in long residence times, tens of minutes long or more, in low-velocity regions. These conclusions suggests that humidity reduction is transport-limited. This analysis provides qualitative insight into airflow structure and humidity transport trends. In addition to characterization of glovebox airflow and humidity transport, results are intended to inform experimental validation and future transient modeling efforts.

42 ENGINEERING↗

Atomic Force Microscopy to Characterize Antimicrobial Peptide-Induced Defects in Model Supported Lipid Bilayers

Antimicrobial peptides (AMPs) interact with bacterial cell membranes through a variety of mechanisms, causing changes extending from nanopore formation to microscale membrane lysis, eventually leading to cell death. Several AMPs also disrupt mammalian cell membranes, despite their significantly different lipid composition and such collateral hemolytic damage hinders the potential therapeutic applicability of the AMP as an anti-microbial. Elucidating the mechanisms underlying the AMP–membrane interactions is challenging due to the variations in the chemical and structural features of the AMPs, the complex compositional variations of cell membranes and the inadequacy of any single experimental technique to comprehensively probe them. (1) Background: Atomic Force Microscopy (AFM) imaging can be used in combination with other techniques to help understand how AMPs alter the orientation and structural organization of the molecules within cell membranes exposed to AMPs. The structure, size, net charge, hydrophobicity and amphipathicity of the AMPs affect how they interact with cell membranes of differing lipid compositions. (2) Methods: Our study examined two different types of AMPs, a 20-amino acid, neutral, α-helical (amphipathic) peptide, alamethicin, and a 13-amino acid, non-α-helical cationic peptide, indolicidin (which intramolecularly folds, creating a hydrophobic core), for their interactions with supported lipid bilayers (SLBs). Robust SLB model membranes on quartz supports, incorporating predominantly anionic lipids representative of bacterial cells, are currently not available and remain to be developed. Therefore, the SLBs of zwitterionic egg phosphatidylcholine (PC), which represents the composition of a mammalian cell membrane, was utilized as the model membrane. This also allows for a comparison with the results obtained from the Quartz Crystal Microbalance with Dissipation (QCM-D) experiments conducted for these peptides interacting with the same zwitterionic SLBs. Further, in the case of alamethicin, because of its neutrality, the lipid charge may be less relevant for understanding its membrane interactions. (3) Results: Using AFM imaging and roughness analysis, we found that alamethicin produced large, unstable defects in the membrane at 5 µM concentrations, and completely removed the bilayer at 10 µM. Indolicidin produced smaller holes in the bilayer at 5 and 10 µM, although they were able to fill in over time. The root-mean-square (RMS) roughness values for the images showed that the surface roughness caused by visible defects peaked after peptide injection and gradually decreased over time. (4) Conclusions: AFM is useful for helping to uncover the dynamic interactions between different AMPs and cell membranes, which can facilitate the selection and design of more efficient AMPs for use in therapeutics and antimicrobial applications.

Swana, Kathleen W.↗

Equation-Free Coarse Control of Distributed Parameter Systems via Local Neural Operators

The control of high-dimensional distributed parameter systems (DPS) remains a challenge when explicit coarse-grained equations are unavailable. Classical equation-free (EF) approaches rely on fine-scale simulators treated as black-box timesteppers. However, repeated simulations for steady-state computation, linearization, and control design are often computationally prohibitive, or the microscopic timestepper may not even be available, leaving us with data as the only resource. We propose a data-driven alternative that uses local neural operators, trained on spatiotemporal microscopic/mesoscopic data, to obtain efficient short-time solution operators. These surrogates are employed within Krylov subspace methods to compute coarse steady and unsteady-states, while also providing Jacobian information in a matrix-free manner. Krylov-Arnoldi iterations then approximate the dominant eigenspectrum, yielding reduced models that capture the open-loop slow dynamics without explicit Jacobian assembly. Both discrete-time Linear Quadratic Regulator (dLQR) and pole-placement (PP) controllers are based on this reduced system and lifted back to the full nonlinear dynamics, thereby closing the feedback loop.

93B52, 93C20, 47N70, 65J15, 65M32, 68T07, 68T20, 6↗

Local conformations and heterogeneities in structures and dynamics of isotactic polypropylene adsorbed onto carbon fiber

Carbon fiber (CF) reinforced polymers (CFRPs) have experienced widespread use in various industries. One of the most important parameters that controls the macroscopic property of CFRPs is the interface between a polymer matrix and CF. There is growing evidence to suggest the formation of a bound polymer layer (BPL), i.e., polymer chains that physically adsorb on a filler surface. However, this interface is always in contact with the thicker part of a polymer matrix, rendering its understanding a difficult task. Therein, we use CF-reinforced isotactic polypropylene (iPP) as a rational CFRP. To characterize the BPL on the CF surface, we extracted it from the CFRP using solvent-rinsing with p-xylene. The physical and thermal properties of the BPL were characterized by differential scanning calorimetry and thermogravimetric analysis, while its microscopic structures and dynamics were probed by small-angle neutron scattering and quasi-elastic neutron scattering (QENS) techniques. Subsequently, we employed atomistic molecular dynamics (MD) simulations to complement the QENS results above the bulk melting temperature and reveal details that were experimentally inaccessible. In this work, we observed that the degree of crystallinity of the BPL was quite lower than the bulk, while the melting temperature of the BPL remained the same as the bulks. Within the given length and time scales probed by QENS and MD, we also observed that most of the bound chains were mobile, with the formation of a high-density region (less than 1 nm in thickness) near the CF surface. The segmental dynamics of the bound chains probed by both QENS and MD were also much faster than those of the free chains, possibly due to the presence of a free surface region at the topmost surface of the BPL. Furthermore, the MD results demonstrated that the backbone chains and side groups lie nearly flat on the CF surface, which is the driving force for the flattening process of the iPP BPL to overcome the conformational entropy loss in the total free energy.

36 MATERIALS SCIENCE↗

From Bulk to Interface: Solvent Exchange Dynamics and Their Role in Ion Transport and the Interfacial Model of Rechargeable Magnesium Batteries

Multivalent battery chemistries have been explored in response to the increasing demand for high-energy rechargeable batteries utilizing sustainable resources. Solvation structures of working cations have been recognized as a key component in the design of electrolytes; however, most structure–property correlations of metal ions in organic electrolytes usually build upon favorable static solvation structures, often overlooking solvent exchange dynamics. We here report the ion solvation structures and solvent exchange rates of magnesium electrolytes in various solvents by using multimodal nuclear magnetic resonance (NMR) analysis and molecular dynamics/density functional theory (MD/DFT) calculations. These magnesium solvation structures and solvent exchange dynamics are correlated to the combined effects of several physicochemical properties of the solvents. Moreover, Mg 2+ transport and interfacial charge transfer efficiency are found to be closely correlated to the solvent exchange rate in the binary electrolytes where the solvent exchange is tunable by the fraction of diluent solvents. Our primary findings are (1) most battery-related solvents undergo ultraslow solvent exchange coordinating to Mg 2+ (with time scales ranging from 0.5 μs to 5 ms), (2) the cation transport mechanism is a mixture of vehicular and structural diffusion even at the ultraslow exchange limit (with faster solvent exchange leading to faster cation transport), and (3) an interfacial model wherein organic-rich regions facilitate desolvation and inorganic regions promote Mg 2+ transport is consistent with our NMR, electrochemistry, and cryogenic X-ray photoelectron spectroscopy (cryo-XPS) results. In conclusion, this observed ultraslow solvent exchange and its importance for ion transport and interfacial properties necessitate the judicious selection of solvents and informed design of electrolyte blends for multivalent electrolytes.

25 ENERGY STORAGE↗

Rain-Induced Stratification of the Equatorial Indian Ocean and Its Potential Feedback to the Atmosphere

Surface freshening through precipitation can act to stably stratify the upper ocean, forming a rain layer (RL). RLs inhibit subsurface vertical mixing, isolating deeper ocean layers from the atmosphere. This process has been studied using observations and idealized simulations. The present ocean modeling study builds upon this body of work by incorporating spatially resolved and realistic atmospheric forcing. Fine-scale observations of the upper ocean collected during the Dynamics of the Madden-Julian Oscillation field campaign are used to verify the General Ocean Turbulence Model (GOTM). Spatiotemporal characteristics of equatorial Indian Ocean RLs are then investigated by forcing a 2D array of GOTM columns with realistic and well-resolved output from an existing regional atmospheric simulation. RL influence on the ocean-atmosphere system is evaluated through analysis of RL-induced modification to surface fluxes and sea surface temperature (SST). This analysis demonstrates that RLs cool the ocean surface on time scales longer than the associated precipitation event. A second simulation with identical atmospheric forcing to that in the first, but with rainfall set to zero, is performed to investigate the role of rain temperature and salinity stratification in maintaining cold SST anomalies within RLs. Approximately one third, or 0.1°C, of the SST reduction within RLs can be attributed to rain effects, while the remainder is attributed to changes in atmospheric temperature and humidity. The prolonged RL-induced SST anomalies enhance SST gradients that have been shown to favor the initiation of atmospheric convection. These findings encourage continued research of RL feedbacks to the atmosphere.

54 ENVIRONMENTAL SCIENCES↗

Turbulent flow characteristics in an 84-pin rod bundle for typical and damaged spacer grids

Hexagonal rod bundles arranged in a tightly packed triangular lattice are extensively used for heat transfer and energy generation applications. Staggered spacer grids are used to maintain the structural integrity of gas-cooled fast reactor (GFR) fuel assemblies, while inducing localized turbulence in flow. Damage to these spacer grids results in a disruption of flow fields within these hexagonal fuel bundles. Experimental flow visualizations are critical to identify the differences in local flow properties that the structural damage may cause. This experimental research investigates the flow-field characteristics at a near-wall and center plane in a prototypical 84-pin GFR fuel assembly. Newly installed typical spacers and spacers subject to naturally occurring damage due to material degradation over prolonged experimentation were investigated. Velocity fields were acquired by utilizing the matched-index-of-refraction method to obtain time-resolved particle image velocimetry measurements for a Reynolds number of 12 000. Reynolds decomposition statistical results divulged differences in the time-averaged velocity, velocity fluctuations, flow anisotropy, and Reynolds stress distributions. Galilean decomposition demarcated the influence of spacer grid damage on the velocity fields. To extract turbulent structures and elucidate mechanisms of flow instabilities, proper orthogonal decomposition analysis was employed. Reduced order flow reconstructions enabled the application of vortex identification algorithms to determine the spatial and statistical characteristics of vortices generated. This research work provides unique experimental data on the spacer grid condition-dependent flow. The results offer a deeper understanding of fluid dynamics behavior to support GFR rod bundle design efforts and computational fluid dynamics model validation.

36 MATERIALS SCIENCE↗

Effect of Hydrogen on Steady-State and Transient Combustion Instability Characteristics

This paper examines the effects of steady-state and transient hydrogen enrichment on thermoacoustic instability in a model gas turbine combustor. Combustion instability, a feedback loop between flame heat release rate oscillations and combustor acoustics, is characterized in a swirl-stabilized flame operated at a range of hydrogen–natural gas fuel blends and heat rates. Measurements of combustor chamber pressure fluctuations and CH* chemiluminescence imaging are used to characterize instability at a range of operating conditions. Steady-state tests show that both mixture heat rate and hydrogen content affect system stability. At a given heat rate, higher levels of hydrogen result in unstable combustion. As heat rate increases, instability occurs at lower concentrations of hydrogen in the fuel. Transient operation was tested in two directions—instability onset and decay—and two hydrogen-addition times—a short time of 1 ms and a longer time of 4 s. Results show that instability onset processes, through the transient addition of hydrogen, are highly repeatable regardless of the timescale of hydrogen addition. Certain instability decay processes are less repeatable, resulting in cases that do not fully transition from unstable to stable combustion despite similar changes in hydrogen fuel flow rate. Flame behavior before, during, and after the transient is characterized using high-speed CH* chemiluminescence imaging. Analysis of the high-speed images shows changes in flame stabilization and dynamics during the onset and decay processes. Finally, the results of this study can have implications for systems that experience variations in fuel composition, particularly in light of growing interest in hydrogen as a renewable fuel.

42 ENGINEERING↗

A New Era of Discovery: The 2023 Long-Range Plan for Nuclear Science (V.1.2)

Nuclear science is the investigation of how protons and neutrons are formed from elementary particles and how the forces between those particles produce both nuclei and the vast variety of nuclear phenomena that occur in the universe. It has evolved into a broad field that addresses profound scientific questions: Where does the mass of visible matter come from? How do stars ignite, live, and die? How do nuclei illuminate the search for new laws of nature? This science points the way to using nuclei to build new technologies that benefit society. The 2015 Nobel Prize in physics was shared by nuclear physicists Art McDonald and Takaaki Kajita for the discovery of neutrino oscillations, which confirmed that neutrinos have mass. Our progress on big questions like this one since 2015 has been remarkable owing to new experimental tools, theoretical breakthroughs, powerful computational techniques, and the talented people who make these innovations possible. Focusing on these new tools, the Facility for Rare Isotope Beams (FRIB) at Michigan State University is already producing exciting results on decays of never-before-produced isotopes a year after it was completed on time and on budget. The energy upgrade of the Continuous Electron Beam Accelerator Facility (CEBAF) at the Thomas Jefferson National Accelerator Facility (Jefferson Lab) was also completed on schedule and on budget—new data from this facility are revealing the spectrum, structure, and dynamics of protons, neutrons, nuclei, and mesons. On the theory front, we can now calculate the distribution of quarks inside the proton from first principles. The implementation of artificial intelligence (AI) and machine learning (ML) techniques has led to improved data analysis and increased efficiency in running experiments and theoretical calculations. The impact of nuclear science goes beyond expanding the frontiers of knowledge about matter in the universe. We simultaneously develop a STEM work force that advances the security, technology, health, and wealth of our nation. Some connections are obvious. Expert scientists trained to work with radioactive nuclei are in demand in nuclear security arenas and are highly sought after by various government agencies and private industries. Graduate students and postdoctoral fellows (postdocs) obtain extensive computational, modeling, and data science skills that are similarly in high demand. Less obvious but equally important is the connection between these trained scientists and success in other professions, including medicine, energy, and entrepreneurial pursuits. The workforce that enables discovery in nuclear science also makes breakthroughs in technologies with tremendous impact on the nation’s economic advancement.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

SUHMO: an adaptive mesh refinement SUbglacial Hydrology MOdel v1.0

Abstract. Water flowing under ice sheets and glaciers can have a strong influence on ice dynamics, particularly through pressure changes, suggesting that a comprehensive ice sheet model should include the effect of basal hydrology. Modeling subglacial hydrology remains a challenge, however, mainly due to the range of spatial and temporal scales involved – from subglacial channels to vast subglacial lakes. Additionally, networks of subglacial drainage channels dynamically evolve over time. To address some of these challenges, we have developed an adaptive mesh refinement (AMR) model based on the Chombo software framework. We extend the model proposed by Sommers et al. (2018) with a small but significant change to accommodate the transition from unresolved to resolved flow features. We handle the strong nonlinearities present in the equations by resorting to an efficient nonlinear full approximation scheme multigrid (FAS-MG) algorithm. We outline the details of the algorithm and present convergence analysis results demonstrating its good performance. Additionally, we present results validating our approach, using test cases from the Subglacial Hydrology Model Intercomparison Project (SHMIP) (de Fleurian et al., 2018). We finish by presenting a more complex, 100 km-by-100 km synthetic test case with peaks and valleys that we use to investigate the effective pressure distribution as the number of AMR levels increases. These preliminary results suggest that a minimum spatial resolution is needed to properly capture channel features, but additional work is required to precisely quantify this and its impact on accurately modeling the coupled ice sheet–hydrology system. The efficiency of our approach, relying on localized refinement, is also demonstrated. Future work will include coupling the SUbglacial Hydrology MOdel (SUHMO) with the BISICLES AMR ice sheet model (Cornford et al., 2013), both built on the same numerical framework.

58 GEOSCIENCES↗

Wave packet dynamics and control in excited states of molecular nitrogen

Wave packet interferometry with vacuum ultraviolet light has been used to probe a complex region of the electronic spectrum of molecular nitrogen, N 2 . Wave packets of Rydberg and valence states were excited by using double pulses of vacuum ultraviolet (VUV), free-electron-laser (FEL) light. These wave packets were composed of contributions from multiple electronic states with a moderate principal quantum number (n ~ 4-9) and a range of vibrational and rotational quantum numbers. The phase relationship of the two FEL pulses varied in time, but as demonstrated previously, a shot-by-shot analysis allows the spectra to be sorted according to the phase between the two pulses. The wave packets were probed by angle-resolved photoionization using an infrared pulse with a variable delay after the pair of excitation pulses. The photoelectron branching fractions and angular distributions display oscillations that depend on both the time delays and the relative phases of the VUV pulses. The combination of frequency, time delay, and phase selection provides significant control over the ionization process and ultimately improves the ability to analyze and assign complex molecular spectra.

74 ATOMIC AND MOLECULAR PHYSICS↗

Effect of causality constraints on Bayesian analyses of heavy-ion collisions

There have long been questions about the limits to the validity of relativistic fluid dynamics and whether it is being used outside its regime of validity in modern simulations of relativistic heavy-ion collisions. An important new tool for answering this question is a causality analysis in the nonlinear regime—if the solutions of the evolution equations do not respect relativistic causality, then they are not a faithful representation of the underlying relativistic theory (in this case, quantum chromodynamics). Using this nonlinear criterion, it has recently been shown that hydrodynamics is indeed being used outside its regime of validity in simulations, at least sometimes. Here we explore the phenomenological implications, particularly the quantitative effects of demanding limits on acausality in modern Bayesian parameter estimation. We find that, while typically only a small fraction of the system's energy is initially in an acausal regime, placing strict limits on the allowed energy fraction significantly changes the preferred properties of the initial condition, which in turn alters the extracted medium properties such as bulk viscosity, where large values are no longer favored. Furthermore, these findings highlight the importance of developing better theoretical descriptions of the early-time, out-of-equilibrium dynamics of relativistic heavy-ion collisions.

Bayesian methods↗