Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “equilibrium modeling”

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,315 records · Page 73

Designing molecular qubits: computational insights into first-row and group 6 transition metal complexes

In the realm of optically addressable qubits, a previously synthesized and characterized Cr( IV ) pseudo-tetrahedral complex, featuring four strongly donating ligands surrounding the chromium center, has demonstrated potential as a qubit candidate. This study proposes analogs of this complex through a metal substitution strategy, extending the investigation to different complexes based on metal centers selected from first-row and Group 6 transition metals. Computational modeling based on multiconfigurational methods CASPT2 and MC-PDFT was utilized to calculate energy gaps between ground and excited electronic spin states, and zero-field splitting parameters. Simulations were applied to each equilibrium geometry and related deformations based on vibrational modes. All results align with previous experimental findings, but also show that qubits based on V and Ti centers could be more electronically stable than the Cr one, suggesting a lower electronic features dependency from their related geometry. In some cases geometrical deformations provide changes in relative energy gaps between triplet and singlet excited state, that could potentially swap, offering a different initialization process, and some inspiration for ligand design based on such deformations. Additionally, this study identifies an unsynthesized Ti( II ) compound as a promising candidate for molecular qubits. This finding highlights the role of computational multireference methods in the rational design of qubit systems.

Sauza-de la Vega, Arturo [Univ. of Chicago, IL (Un↗

Generalized impurity pinch in partially magnetized multi-ion plasma

In a two-ion-species plasma with disparate ion masses, heavy ions tend to concentrate in the low-temperature region of collisionally magnetized plasma and in the high-temperature region of collisionally unmagnetized plasma, respectively. Moreover, collisional magnetization can be determined as the ratio of the light ion gyrofrequency to the collision frequency of light and heavy ion species, and the behavior of this effect in the intermediate regime of partially magnetized plasma is predominantly dependent on this Hall parameter. Multi-ion cross-field transport has been described before in the collisionally magnetized plasma regime, and generalized pinch relations, which describe densities of ion species in equilibrium in that plasma, are found in the literature. Here, the role of collisional magnetization and Larmor magnetization in multi-ion collisional transport is clarified, and generalized pinch relations are extended to the partially magnetized regime in which the ion Hall parameter may be small, as long as electrons remain collisionally magnetized. Equilibrium ion density profiles have the same dependence on external forces and on each other regardless of collisional magnetization of ions. The expansion of the range of validity of multi-ion collisional transport models makes them applicable to a wider range of laboratory plasma conditions. In particular, ion density profiles evolve sufficiently fast for radial impurity transport to be observable around stagnation on MagLIF, leading to expulsion of heavy ion impurities from the hotspot as long as plasma becomes sufficiently collisionally magnetized during the implosion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Excited quantum anomalous and spin Hall effect: dissociation of flat-bands-enabled excitonic insulator state

Quantum anomalous Hall effect (QAHE) and quantum spin Hall effect (QSHE) are two interesting physical manifestations of 2D materials that have an intrinsic nontrivial band topology. In principle, they are ground-state equilibrium properties characterized by Fermi level lying in a topological gap, below which all the occupied bands are summed to a non-zero topological invariant. In this work, we propose theoretical concepts and models of ‘excited’ QAHE (EQAHE) and EQSHE generated by dissociation of an excitonic insulator (EI) state with complete population inversion (CPI), a unique many-body ground state enabled by two yin-yang flat bands (FBs) of opposite chirality hosted in a diatomic Kagome lattice. The two FBs have a trivial gap in between, i.e. the system is a trivial insulator in the single-particle ground-state, but nontrivial gaps above and below, so that upon photoexcitation the quasi-Fermi levels of both electrons and holes will lie in a nontrivial gap achieved by the CPI-EI state, as demonstrated by exact diagonalization calculations. Then dissociation of singlet and triplet EI state will lead to EQAHE and EQSHE, respectively. Realizations of yin-yang FBs in real materials are also discussed.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Direct prediction of saturated neoclassical tearing modes in slab using an equilibrium approach

We demonstrate for the first time that the nonlinear saturation of neoclassical tearing modes (NTMs) can be found directly using a variational principle based on Taylor relaxation, without needing to simulate the intermediate, resistivity-dependent dynamics. As in previous investigations of classical tearing mode saturation (Loizu et al 2020 Phys. Plasmas 27 070701; Loizu and Bonfiglio 2023 J. Plasma Phys. 89 905890507), we make use of Stepped Pressure Equilibrium Code (SPEC) (Hudson et al 2012 Phys. Plasmas 19 112502), an equilibrium solver based on the variational principle of the multi-region relaxed magnetohydrodynamics (MHDs), featuring stepped pressure profiles and arbitrary magnetic topology. We work in slab geometry and employ a simple bootstrap current model J bs = C$\boldsymbol{\nabla}$p to study the bootstrap-driven tearing modes, scanning over the asymptotic matching parameter Δ' and bootstrap current strength. Saturated island widths produced by SPEC agree well with the predictions of an initial value resistive MHDs code (Huang and Bhattacharjee 2016 Astrophys. J. 818 20) while being orders of magnitude faster to calculate. Additionally, we observe good agreement with a simple analytical modified Rutherford equation, without requiring any fitting coefficients. The match is obtained for both linearly unstable classical tearing modes in the presence of bootstrap current, and NTMs, which are linearly stable but nonlinear-unstable due to the effects of the bootstrap current.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multi-messenger heavy-ion collision physics

We report that this work studies the production of direct photons in relativistic nuclear collisions, along with the production of hadrons. Radiation from the very first instants to the final moments of the evolution is included. The hybrid model used here describes all stages of relativistic heavy-ion collisions. Chronologically, those are an initial state reflecting the collision of nuclei described within the Color Glass Condensate effective theory; a pre-equilibrium phase based on non-equilibrium linear response; relativistic viscous hydrodynamics, and a hadronic afterburner. The effect of the preequilibrium phase on both photonic and hadronic observables is highlighted for the first time. The potential of photon observables - spectrum, differential elliptic and triangular flow - to reveal the chemical equilibration time is studied. Finally, we consider "small collision systems", including proton+nucleus collisions and collisions of light nuclei, as probed by hadronic and electromagnetic observables. We demonstrate how photon production can signal the formation of quark-gluon plasma in such small systems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Leaching Kinetics of Rare Earth Elements from Fire Clay Seam Coal

Recovery of rare earth elements (REEs) from coal samples collected from the Fire Clay coal seam using diluted mineral acid solutions was investigated. The initial processing step was coal recovery using conventional froth flotation which concentrated the REEs in tailing material resulting in an upgrade to values around 700 ppm on a dry whole mass basis. Leaching experiments were performed on the flotation tailings material using a 1.2 M sulfuric acid solution adjusted to a temperature of 75 °C to study the extractability of REEs from coal material. The effect of particle size, leaching time, leaching temperature, and solid concentration on REE leaching recovery were evaluated. The kinetic data obtained from leaching over a range of temperatures suggested that the leaching process follows the shrinking core model with possibly a mixed control mechanism that may be a result of several heterogenous materials leaching simultaneously. Leaching recovery increased rapidly at the beginning of the reaction then slowed as the system reached equilibrium. The apparent activation energy determined from test data obtained over a range of temperatures using 1 M sulfuric acid was 36 kJ/mol for the first 20 min of reaction time and 27 kJ/mol for the leaching period between 20 and 120 min. The leaching of light REEs during the initial stage was determined to be driven by a chemical reaction, followed by the formation of a product layer, which required lower activation energy in the later stage of leaching. In regards to the heavy REEs, the major mechanism for leaching is desorption and the product layer formation does not affect the heavy REEs significantly.

58 GEOSCIENCES↗

On the Utility of Transmission Color Analysis i: Differentiating Super-Earths and Sub-Neptunes

The majority of exoplanets found to date have been discovered via the transit method, and transmission spectroscopy represents the primary method of studying these distant worlds. Currently, in-depth atmospheric characterization of transiting exoplanets entails the use of spectrographs on large telescopes, requiring significant observing time to study each planet. Previous studies have demonstrated trends for solar system worlds using color–color photometry of reflectance spectra, as well as trends within transmission spectra for hot Jupiters. Building on these concepts, we have investigated the use of transmission color photometric analysis for efficient, coarse categorization of exoplanets and for assessing the nature of these worlds, with a focus on resolving the bulk composition degeneracy to aid in discriminating super-Earths and sub-Neptunes. We present our methodology and first results, including spectrum models, model comparison frameworks, and wave band selection criteria. We present our results for different transmission “color” metrics, filter selection methods, and numbers of filters. Assuming noise-free spectra of isothermal atmospheres in chemical equilibrium, with our pipeline, we are able to constrain atmospheric mean molecular weight in order to distinguish between super-Earth and sub-Neptune atmospheres with >90% overall accuracy using specific low-resolution filter combinations, . We also found that increasing the number of filters does not substantially impact this performance. This method could allow for broad characterization of large numbers of planets much more efficiently than current methods permit, enabling population- and system-level studies. Additionally, data collected via this method could inform follow-up observing time by large telescopes for more detailed studies of worlds of interest.

47 OTHER INSTRUMENTATION↗

Near-continuum, hypersonic oxygen flow over a double cone simulated by direct simulation Monte Carlo informed from quantum chemistry

A large-scale, fully resolved direct simulation Monte Carlo (DSMC) computation of a non-equilibrium, reactive flow of pure oxygen over a double cone is presented. Under the simulated near-continuum conditions, the computational demands are shown to be significant because of the wide range of length scales that must be resolved. Therefore, robust grid adaption capabilities and efficient parallelization of the Stochastic PArallel Rarefied-gas Time-accurate Analyzer (SPARTA) code that is utilized in this work are essential. The thermochemical and transport collision models were selected for efficiency and simplicity. First-principles data, obtained from the highly accurate direct molecular simulation method, were used to inform the collision models’ parameters. Importantly, because SPARTA implements molecular collision models using collision-specific energies, the resulting macroscopic relaxation rates were evaluated a posteriori via zero-dimensional heat bath simulations. The comparisons of surface properties, namely heat flux and pressure, show very close agreement with previous computational fluid dynamics (CFD) results. Differences with the measurements were found to be similar to the CFD simulations. The unresolved discrepancy with the measurements could be due to inconsistent free stream conditions with the actual experimental data or missing physical phenomena altogether, for example atomic and molecular oxygen electronically excited states, three-dimensional effects, or more complex gas–surface interactions. As shown in this work, the advantages of obtaining a DSMC particle solution for these flows reside in the method's ability to be directly informed from first principles and to seamlessly describe internal energy non-equilibrium for all modes. With the advent of exascale computing and beyond, particle methods will be an increasingly important tool to verify the validity of physical assumptions in reduced-order models via fully resolved, experimental-scale simulations, down to the level of molecular-level distributions.

Mechanics↗

HydraGNN_Predictive_GFM_2024 - Ensemble of predictive graph foundation models for ground state atomistic materials modeling

We provide the ensemble of fifteen pre-trained graph foundation models (GFMs) for atomistic materials modeling applications. Each one of the fifteen GFMs has been trained on five open-source datasets that (once aggregated) amount to over 154 million atomistic structures, which cover over two-thirds of the natural elements of the periodic table and that comprises a broad set of organic and inorganic compounds. This vast set of atomistic structures comprises ground state configurations that are dynamically stable (i.e., equilibrated structures with atomic forces approximately close to zero values) as well as dynamically unstable structures (i.e., non-equilibrium structures with non-negligible non-zero values of atomic forces). The ensemble of datasets aggregated does NOT include excited states. The datasets have been curated to remove atomistic structures with spectral norm of the force tensor above 100 eV/angstrom. Moreover, a linear term of the energy was computed for each dataset using a linear regression model that uses the chemical concentration of each natural element as regressor. The linear term predicted by the linear regression model has been subtracted from each original energy value to perform a re-alignment of the energy values across different electronic structures approximation theories performed to generate the diverse multi-source, multi-fidelity datasets. The folder "ADIOS_files" contains the set of pre-processed datasets in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used for the development and training of GFMs in this work. The "ADIOS_files" directory contains 6 sub-directories named as follows: - ANI1x-v3.bp - MPTrj-v3.bp - OC2020-20M-v3.bp - OC2020-v3.bp - OC2022-v3.bp - qm7x-v3.bp Each sub-directory contains the pre-processed datasets converted in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used to the development, training, and performance testing of the ensemble go predictive graph foundation models. Each GFM was developed using HydraGNN (https://github.com/ORNL/HydraGNN) as underlying graph neural network (GNN) architecture. The multi-task learning (MTL) capability of HydraGNN was used to simultaneously train the GFMs on labeled values for direct predictions of energy (a total system property of an atomistic structure that measures the chemical stability) and atomic forces (an atomic level property of an atomistic structure that measures the dynamical stability). The hyper parameters of the GFM have been tuned using scalable hyperparameter optimization (HPO) algorithms implemented in the software DeepHyper (https://github.com/deephyper/deephyper). The pre-training of each HPO trial was performed using distributed data parallelism (DDP) to scale the training across 128 compute nodes of the exascale OLCF supercomputer Frontier. Each HPO trial was trained only for 10 epochs and an early stopping was performed to avoid wasting significant computational resources on GNN architectures that were clearly underperforming. For each HPO trial, the 'omnistat' tool developed by (AMD Research - Advanced Micro Device) was used to measure the total energy consumption in kWh. The ensemble of GFMs was obtained by selecting the fifteen best performing HPO trials. Four models have been selected for their clear advantage in accuracy, and these are the GFMs with IDs 229, 156, 147, 260. Additional eleven models have been selected based on judicious balance between accuracy and energy consumption needed for training, and these are the GFMs with IDs 165, 78, 137, 1, 175, 171, 181, 67, 179, 167, 351. Each selected GFM of the ensemble was continued to cumulate a total of at most 30 epochs. In some cases, the total number of epochs actually performed was les than 30 due to two combined factors: (1) the size of the GFM (i.e., the number of model parameters to train) and (2) the total wall-clock time for which the computational resources could be allocated on OLCF-Frontier. The "Ensemble_of_models" directory contains 15 sub-directories named as follows: - gfm_0.229 - gfm_0.156 - gfm_0.147 - gfm_0.260 - gfm_0.165 - gfm_0.78 - gfm_0.137 - gfm_0.1 - gfm_0.175 - gfm_0.171 - gfm_0.181 - gfm_0.67 - gfm_0.179 - gfm_0.167 - gfm_0.351 Each one of these sub-directories refers to one of the fifteen HPO trials that have been selected to continue the pre-training with at most 30 epochs. With each sub-directory associated with a specific HPO trial, the following files can be found: - config.json: file for argument parsing to develop and train an HydraGNN architecture - gfm_0.ID_epoch_N.pk: file with model parameters for HPO ID trial after N epochs of training The ensemble of fifteen GFM architectures was used for (1) ensemble averaging to stabilize the predictions of energy and atomic forces after pre-training for post-processing analysis and (2) ensemble uncertainty quantification (UQ). The code used to develop, pre-train, and load the pre-trained models for post-processing analysis is available on the ORNL-GitHub at the following link: https://github.com/ORNL/HydraGNN/tree/Predictive_GFM_2024

36 MATERIALS SCIENCE↗

Dynamic Heterogeneity and Kovacs’ Memory Effects in Supercooled Water

Understanding the properties of supercooled water is important for developing a comprehensive theory for liquid water and amorphous ices. Because of rapid crystallization for deeply supercooled water, experiments on it are typically carried out under conditions in which the temperature and/or pressure are rapidly changing. As a result, information on the structural relaxation kinetics of supercooled water as it approaches (metastable) equilibrium is useful for interpreting results obtained in this experimentally challenging region of phase space. Here, we used infrared spectroscopy and the fast time resolution obtained by transiently heating nanoscale water films to investigate relaxation kinetics (aging) in supercooled water. When the structural relaxation of the water films was followed using a temperature jump protocol analogous to the classic experiments of Kovacs, similar memory effects were observed. In particular, after suitable aging at one temperature, water’s structure displayed an extremum versus the number of heat pulses upon changing to a second temperature before eventually relaxing to a steady-state structure characteristic of that temperature. A random double well model based on the idea of dynamic heterogeneity in supercooled water accounts for the observations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Optimizing complexity in the kinetic modelling of integrated flue gas purification for pressurized oxy-combustion

A 5-step mechanism and a corresponding kinetic model are presented for integrated flue gas purification, where SOx and NOx are removed during pressurized oxy-coal combustion. Three coupled subsystems (‘blocks”) of gas–liquid transformations are considered: 1) oxidation of NO to NO 2 and dissolution of NO 2 in water; 2) dissolution of SO 2 in the water; and 3) interactions between S- and N- containing substances in the liquid phase. The model has been formulated in terms of only species that can be experimentally measured (model of optimal complexity). Subsystem I is limiting regarding NO removal. Subsystem II is very fast in comparison with the other subsystems and can be treated separately. Analytical expressions for all independent subsystems have been found and compared with computational results for the whole system. The temporal domain of validity of all these subsystems is presented. Finally, the validity of the NO 2 quasi-steady-state assumption is verified, and it is reported that variables such as NO, NO 2 and HNO 3 can be considered decoupled from the whole system since they appear in only Subsystem I, while HSO 3 - and HNO 2 are coupled. The presented kinetic model was verified experimentally by determining the 1:1 ratio of HNO 2 and HNO 3 , in gas–liquid NO-O 2 interaction. Additionally, new experimental facts on kinetic orders and stoichiometry between HNO 2 and HSO 3 - have been used to optimize the model. It is concluded that all NO x - SO x interactions can be essentially grasped by a sequence of only 4 irreversible reactions considering, in addition, the 5th reaction, i.e. the fast ‘SO 2 -water’ equilibrium.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Grain boundary metastability controls irradiation resistance in nanocrystalline metals

Grain boundaries (GBs) in polycrystalline materials are powerful sinks for irradiation defects. While standard theories assume that a GB’s efficiency as a sink is defined solely by its character before irradiation, recent evidence conclusively shows that the irradiation sink efficiency is a highly dynamic property controlled by the intrinsic metastability of GBs under far-from-equilibrium irradiation conditions. In this paper, we reveal that the denuded (i.e., defect-free) zone, typically the signature of a strong sink, can collapse as irradiation damage accumulates. We propose a radiation damage evolution model that captures this behavior based on the emergence of a series of irradiation defect-enabled metastable GB microstate changes that dynamically alter the ability of the GB to absorb further damage. We show that these microstate changes control further defect absorption and give rise to the formation of a defect network that manifests itself as a net Nye-tensor signal detectable via lattice curvature experiments.

36 MATERIALS SCIENCE↗

Self-Inhibition Phenomena in Cu 3 Pt Oxidation by CO 2

Here, this study investigates the oxidation behavior of Cu 3 Pt(100) in CO 2 using a combination of ambient-pressure X-ray photoelectron spectroscopy, mass spectroscopy, and density functional theory modeling. Our in situ measurements reveal the simultaneous oxidation and reduction of Cu 2 O due to the opposing effects of atomic oxygen and CO generated from dissociative CO 2 adsorption, leading to a dynamic equilibrium state of simultaneously occurring redox reactions. Complementary atomistic calculations elucidate the inhibitory effects of subsurface Pt enrichment and the counteracting roles of CO 2 and CO in surface oxidation and reduction. These results provide mechanistic insights into the dissociative pathway of CO 2 molecules and dynamic evolution of surface composition and reactivity of Cu-based alloy catalysts in CO 2 -rich environments, with broader implications for tuning gas–surface reactions by manipulating gas reactants or solid surface composition.

36 MATERIALS SCIENCE↗

Direct observation of a dynamical glass transition in a nanomagnetic artificial Hopfield network

Spin glasses, generally defined as disordered systems with randomized competing interactions, are a widely investigated complex system. Theoretical models describing spin glasses are broadly used in other complex systems, such as those describing brain function, error-correcting codes5 or stock-market dynamics. This wide interest in spin glasses provides strong motivation to generate an artificial spin glass within the framework of artificial spin ice systems. Here we present the experimental realization of an artificial spin glass consisting of dipolar coupled single-domain Ising-type nanomagnets arranged onto an interaction network that replicates the aspects of a Hopfield neural network. Using cryogenic X-ray photoemission electron microscopy (XPEEM), we performed temperature-dependent imaging of thermally driven moment fluctuations within these networks and observed characteristic features of a two-dimensional Ising spin glass. Specifically, the temperature dependence of the spin glass correlation function follows a power-law trend predicted from theoretical models on two-dimensional spin glasses. As a result, we observe clear signatures of the hard-to-observe rugged spin glass free energy in the form of sub-aging, out-of-equilibrium autocorrelations and a transition from stable to unstable dynamics.

36 MATERIALS SCIENCE↗

Dislocation-Grain Boundary Interaction Dataset for FCC Cu

Interactions between dislocations and grain boundaries play a major role in controlling the strength and ductility of structural materials. Experimentally, assessing and probing geometric and stress-based criteria at the local level for dislocation transmission through grain boundaries remains challenging. Therefore, there have been many efforts to systematically generate datasets of dislocation-grain boundary interactions (DGI) via computational models such as molecular dynamics simulations. So far, most DGI datasets have focused only on the subset of nominal minimum-energy grain boundary structures, which limits their applicability, especially to materials processed far from equilibrium. We present a comprehensive database of dislocation-grain boundary interactions for edge, screw, and 60° mixed dislocation with 330 <110> and 257 <112> symmetric tilt grain boundaries (total of 587) in FCC Cu consisting of 73 minimum-energy grain boundary structures and 514 metastable structures. The dataset contains the outcomes for 5234 unique interactions for various dislocation types, grain boundary structures, and applied shear stresses.

36 MATERIALS SCIENCE↗

Multi-scale computational screening and mechanistic insights of cyclic amines as solvents for improved lignocellulosic biomass processing

A computational screening workflow for the efficient deconstruction of cellulose, lignin and hemicellulose fractions of lignocellulosic biomass using cyclic amines as solvents. Lignocellulosic biomass is a promising feedstock for production of affordable fuels and chemicals from renewable resources. Effective solubilization and subsequent deconstruction of its cellulose, hemicellulose, and lignin fractions is essential for the viability of future biorefineries. This study used quantum chemistry-based equilibrium thermodynamics methods to evaluate the potential of 650 cyclic amines to solubilize cellulose, hemicellulose, and lignin. The activity coefficients of solvent - biopolymer interactions were predicted using the COSMO-RS (COnductor-like Screening MOdel for Real Solvents) method and used to identify cyclic amines that can efficiently dissolve and extract selective fractions of biopolymers during biomass pretreatment. Among the 650 cyclic amines, 1-piperazineethanmaine was predicted to be an effective solvent for extracting all three polymers and was experimentally shown to achieve the highest lignin removal (97.1%). Non-covalent interaction, reduced density gradient and quantum chemical calculations were performed to elucidate the dissolution mechanism of lignin, cellulose and hemicellulose and gain further molecular level insights into the interactions between the cyclic amines and biomass polymers that promote efficient solubilization and extraction. These analyses indicated that 1-piperazineethanmaine and 1-methylimidazole make noncovalent van der Waals, electrostatic interactions and hydrogen bonding with lignin, leading to enhanced lignin removal, while the strong intramolecular hydrogen bonding interactions in cellulose and hemicellulose result in weaker solvent-biopolymer interactions. Overall, the computational approach provided an efficient method for identifying cyclic amines tailored for optimal biomass pretreatment and resulted in the identification of a potential new class of solvents for effective biomass pretreatment.

Kumar, Nikhil↗

Nonlinear modeling of ELM mitigation with RMP on HL-2A

Abstract Nonlinear modeling of mitigation of the edge localized mode (ELM) with resonant magnetic perturbation (RMP) is performed for the HL-2A tokamak, utilizing the three-dimensional (3D) magnetohydrodynamic code JOREK. Based on the 3D equilibrium established after application of the n = 1 ( n is the toroidal mode number) RMP at 4.9 kAt coil current with odd parity, ELM mitigation is successfully simulated consistent with the experimental result. Nonlinear simulations show strong mode coupling among toroidal Fourier harmonics, allowing redistribution of the magnetic energy such that the most unstable toroidal mode saturates at a lower level. This magnetic energy cascade offers an explanation of the RMP-induced ELM mitigation achieved in HL-2A. Detailed examination of the simulation results shows persistent resonant field screening even during the ELM mitigation phase. Finite plasma resistivity however does enable partial penetration of the resonant field thus modifying the edge magnetic topology and characteristics of the edge transport. Plasma radial profiles undergo pronounced changes around the pedestal region, when the magnetic energy of the most unstable toroidal mode reaches the maximum value. Systematic scans of the applied RMP coil current with the JOREK simulations find a threshold value of around 4.5 kAt required for achieving the ELM mitigation on HL-2A.

Physics↗

Recombination of 𝐵 𝑐 mesons in ultrarelativistic heavy-ion collisions

High-energy heavy-ion collisions have been suggested as a favorable environment for the production of 𝐵 𝑐 mesons, due to a much larger abundance of charm and bottom quarks compared to elementary reactions. Motivated by recent CMS data for 𝐵$^{+}_{𝑐}$ production in Pb-Pb(5.02 TeV) collisions at the Large Hadron Collider (LHC), we deploy a previously developed transport approach for charmonia and bottomonia to evaluate the kinetics of 𝐵 𝑐 mesons throughout the fireball formed in these reactions. The main inputs to our approach are two transport parameters: the 𝐵 𝑐 's reaction rate and equilibrium limit. Both quantities are determined by previous calculations via a combination of charm and bottom sectors. In-medium binding energies of 𝐵 𝑐 mesons are calculated from a thermodynamic 𝑇 matrix with a lattice-QCD constrained potential, and figure in their inelastic reaction rates. Temperature-dependent equilibrium limits include charm- and bottom-quark fugacities based on their initial production. We compute the centrality dependence of inclusive 𝐵 𝑐 production and transverse-momentum (𝑝 𝑇 ) spectra using two different recombination models: instantaneous coalescence and resonance recombination. The main uncertainty in the resulting nuclear modification factors, 𝑅 AA , is currently associated with the 𝐵 𝑐 cross section in elementary 𝑝⁢𝑝 collisions, caused by the uncertainty in the branching ratio for the 𝐵$^{−}_{𝑐}$ → 𝐽/𝜓⁢𝜇 − ⁢$\overline{𝜈}$ decay. Our results indicate a large enhancement of the 𝑅 AA at low 𝑝 𝑇 , with significant regeneration contributions up to 𝑝 𝑇 ≃ 20GeV. Comparisons to CMS data are carried out but firm conclusions will require a more accurate value of the branching ratio, or alternative channels to measure the 𝐵 𝑐 production in 𝑝⁢ collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗