Search

Abstract: Following the reconstruction of the TEXT tokamak at Huazhong University of Science and Technology in China, renamed as J-TEXT, a plethora of experimental and theoretical investigations has been conducted to elucidate the intricacies of turbulent transport within the tokamak configuration. These endeavors encompass not only the J-TEXT device’s experimental advancements but also delve into critical issues pertinent to the optimization of future fusion devices and reactors. The research includes topics on the suppression of turbulence, flow drive and damping, density limit, non-local transport, intrinsic toroidal flow, turbulence and flow with magnetic islands, turbulent transport in the stochastic layer, and turbulence and zonal flow with energetic particles or helium ash. Several important achievements have been made in the last few years, which will be further elaborated upon in this comprehensive review.

Since the first plasma realized in 2020, a series of key systems on HL-3 (known as HL-2M before) tokamak have been equipped/upgraded, including in-vessel components (the first wall, lower divertor, and toroidal cryogenic/water-cooling/baking/glow discharge systems, etc.), auxiliary heating system of 11 MW, and 28 diagnostic systems (to measure the plasma density electron temperature, radiation, magnetic field, etc.). Magnet field systems were commissioned firstly for divertor plasma discharges. During the 2nd experimental campaign of HL-3 tokamak, several great progresses have been achieved. Firstly, the successful operation with plasma current larger than 1 MA was achieved under a divertor configuration. Secondly, the advanced divertor concept with two distinct snowflake configurations was realized. It is found that the distribution of ion saturation current and heat flux on bottom plate becomes wide due to magnetic surface expansion, demonstrating the advantage of such configuration in the heat flux mitigation. In addition, using the combination of NBI, ECRH and LHCD, the standard sawtoothing high confinement mode of megampere plasma was firstly accessed on the HL-3. The successful commissioning of HL-3 is beneficial for the initial operation of ITER.

Since the last IAEA-FEC in 2021, significant progress on the development of long pulse steady state scenario and its related key physics and technologies have been achieved, including the reproducible 403 s long-pulse steady-state H-mode plasma with pure radio frequency (RF) power heating. A thousand-second time scale (∼1056 s) fully non-inductive plasma with high injected energy up to 1.73 GJ has also been achieved. The EAST operational regime of high βP has been significantly extended (H98y2 > 1.3, βP ∼ 4.0, βN ∼ 2.4 and ne/nGW ∼ 1.0) using RF and neutral beam injection (NBI). The full edge localized mode suppression using the n = 4 resonant magnetic perturbations has been achieved in ITER-like standard type-I ELMy H-mode plasmas with q95 ≈ 3.1 on EAST, extrapolating favorably to the ITER baseline scenario. The sustained large ELM control and stable partial detachment have been achieved with Ne seeding. The underlying physics of plasma-beta effect for error field penetration, where toroidal effect dominates, is disclosed by comparing the results in cylindrical theory and MARS-Q simulation in EAST. Breakdown and plasma initiation at low toroidal electric fields (

We study the effect of strain on the magnetic properties and magnetization configurations in nanogranular FexGe1-x films (x = 0.53 ±0.05) with and without B20 FeGe nanocrystals surrounded by an amorphous structure. Relaxed films on amorphous silicon nitride membranes reveal a disordered skyrmion phase while films near and on top of a rigid substrate favor ferromagnetism and an anisotropic hybridization of Fe d levels and spin-polarized Ge sp band states. The weakly coupled topological states emerge at room temperature and become more abundant at cryogenic temperatures without showing indications of pinning at defects or confinement to individual grains. These results demonstrate the possibility to control magnetic exchange and topological magnetism by strain and inform magnetoelasticity-mediated voltage control of topological phases in amorphous quantum materials.

In this study, we analyse ‘magneto-Stokes’ flow, a fundamental magnetohydrodynamic (MHD) flow that shares the cylindrical-annular geometry of the Taylor–Couette cell but uses applied electromagnetic forces to circulate a free-surface layer of electrolyte at low Reynolds numbers. The first complete, analytical solution for time-dependent magneto-Stokes flow is presented and validated with coupled laboratory and numerical experiments. Three regimes are distinguished (shallow-layer, transitional and deep-layer flow regimes), and their influence on the efficiency of microscale mixing is clarified. The solution in the shallow-layer limit belongs to a newly identified class of MHD potential flows, and thus induces mixing without the aid of axial vorticity. We show that these shallow-layer magneto-Stokes flows can still augment mixing in distinct Taylor dispersion and advection-dominated mixing regimes. The existence of enhanced mixing across all three distinguished flow regimes is predicted by asymptotic scaling laws and supported by three-dimensional numerical simulations. Mixing enhancement is initiated with the least electromagnetic forcing in channels with order-unity depth-to-gap-width ratios. If the strength of the electromagnetic forcing is not a constraint, then shallow-layer flows can still yield the shortest mixing times in the advection-dominated limit. Our robust description of momentum evolution and mixing of passive tracers makes the annular magneto-Stokes system fit for use as an MHD reference flow.

This paper presents a machine learning-based airfoil shape optimization framework designed to increase flutter resistance and reduce drag. Using the National Advisory Committee for Aeronautics airfoil as the base design and a Hicks–Henne bump function, we employ multi-objective Bayesian optimization and harmonic balance-based flutter prediction. The optimization process yields a Pareto front revealing trade-off relationships between the flutter speed index and drag coefficient. The optimized airfoils, resembling those of evolved marine animals, outperform the base design in terms of flutter resistance and drag. These results demonstrate the framework's potential to enhance aircraft performance and safety by addressing aeroelastic factors.

Spot fires pose a major risk and add to the already complex physics, which makes fire spread so hard to predict, especially in the wildland urban interface. Firebrands can not only cross fuel breaks and thwart other suppression efforts but also directly damage infrastructure and block evacuation routes. Transport models and computational fluid dynamics tools often make simplifications when predicting spot fire risk, but there is a relative lack of experimental data to validate such parameterizations. To this end, we present a field experiment performed at the University of California Berkeley Blodgett Research Forest in California where we recorded the flame and firebrands emanating from a nighttime hand-drawn pile fire using high-frequency imaging. We used image-processing to characterize the fire intensity and turbulence as well as particle tracking velocimetry to measure ejected firebrand kinematics as they are lofted by the plume. We further collected embers that settled around the fire at varying distances and measured their size, shape, density, and settling distributions. We also examine existing physics-based time-averaged models of firebrand lofting and note discrepancies between such models, often used due to their speed and simplicity, and our experimental observations. Finally, we discuss some implications our observations could have on future modeling efforts by considering the time-dependent fire dynamics, intermittency in the plume turbulence, and in the firebrand generation rate. To the best of our knowledge, these are the first in situ observations of firebrand generation and lofting from representative fuels, addressing a major source of data gap and uncertainty in the wildland fire literature.

Laboratory experiments and particle-resolved simulations are employed to investigate the settling dynamics of a pair of rigidly connected spherical particles of unequal size. They yield a detailed picture of the transient evolution and the terminal values of the aggregate's orientation angle and its settling and drift velocities as functions of the aspect ratio and the Galileo number $Ga$ , which denotes the ratio of buoyancy and viscous forces acting on the aggregate. At low to moderate values of $Ga$ , the aggregate's orientation and velocity converge to their terminal values monotonically, whereas for higher $Ga$ -values the aggregate tends to undergo a more complex motion. If the aggregate assumes an asymmetric terminal orientation, it displays a non-zero terminal drift velocity. For diameter ratios much larger than one and small $Ga$ , the terminal orientation of the aggregate becomes approximately vertical, whereas when $Ga$ is sufficiently large for flow separation to occur, the aggregate orients itself such that the smaller sphere is located at the separation line. Empirical scaling laws are obtained for the terminal settling velocity and orientation angle as functions of the aspect ratio and $Ga$ for diameter ratios from 1 to 4 and particle-to-fluid density ratios from 1.3 to 5. An analysis of the accompanying flow field shows the formation of vortical structures exhibiting complex topologies in the aggregate's wake, and indicates the formation of a horizontal pressure gradient across the larger sphere, which represents the main reason for the emergence of the drift velocity.

A real-time capable core Ion Cyclotron Range of Frequencies (ICRF) heating model on NSTX and WEST is developed. The model is based on two nonlinear regression algorithms, the random forest ensemble of decision trees and the multilayer perceptron neural network. The algorithms are trained on TORIC ICRF spectrum solver simulations of the expected flat-top operation scenarios in NSTX and WEST assuming Maxwellian plasmas. The surrogate models are shown to successfully capture the multi-species core ICRF power absorption predicted by the original model for the high harmonic fast wave and the ion cyclotron minority heating schemes while reducing the computational time by six orders of magnitude. Although these models can be expanded, the achieved regression scoring, computational efficiency and increased model robustness suggest these strategies can be implemented into integrated modeling frameworks for real-time control applications.

Abstract: Long-wavelength density fluctuations ( k ρ i 0.97, indicating the prevalence of density ‘voids’ at inner radii and density ‘blobs’ at outer radii and outside of the separatrix. The turbulence intensity flux ⟨ v ~ r n ~ 2 ⟩ is calculated to characterize turbulence spreading at the plasma edge. During ECH/NBI power ramps and at counter- I p injected torque, ⟨ v ~ r n ~ 2 ⟩ is directed inward inside the separatrix, which is evidence of inward spreading of turbulence intensity from the edge gradient region caused by the inner propagation of density ‘voids’. Significantly weaker ⟨ v ~ r n ~ 2 ⟩ is observed with co- I p torque. A correlation between co- I p torque, turbulence intensity δ n / n at ρ = 0.97, and increased srape-off layer (SOL) heat flux decay length λ q is found in the torque scan scenario, showing that edge turbulence plays a material role in determining the SOL conditions and heat flux width.

Charge conservation and the Pauli exclusion principle result from fundamental symmetries in the standard model of particle physics, and are typically taken as axiomatic. High-precision tests for small violations of these symmetries could point to new physics. Here we consider three models for violation of these processes, which would produce detectable ionization in the high-purity germanium detectors of the Majorana Demonstrator experiment. Using a 37.5 kg yr exposure, we report a lower limit on the electron mean lifetime, improving the previous best limit for the e→νeνe¯νe decay channel by more than an order of magnitude. We also present searches for two types of violation of the Pauli exclusion principle, setting limits on the probability of an electron to be found in a symmetric quantum state.

Landauer’s principle states that erasing a bit of information at fixed temperature T costs at least k B T ln ⁡ 2 units of work. Here we investigate erasure at varying temperature, to which Landauer’s result does not apply. We formulate bit erasure as a stochastic nonequilibrium process involving a compression of configuration space, with physical and logical states associated in a symmetric way. Erasure starts and ends at temperature T, but temperature can otherwise vary with time in an arbitrary way. Defined in this way, erasure is governed by a set of nonequilibrium fluctuation relations that show that varying-temperature erasure can done with less work than k B T ln ⁡ 2 . As a result, erasure and the complementary process of bit randomization can be combined to form a work-producing engine cycle.

In this article we investigate novel signatures of radiation reaction via the angular deflection of an electron beam colliding at 90 degrees with an intense laser pulse. Due to the radiation reaction effect, the electrons can be deflected towards the beam axis for plane wave backgrounds, which is not possible in the absence of radiation reaction effects. The magnitude and size of the deflection angle can be controlled by tailoring the laser pulse shapes. The effect is first derived analytically using the Landau-Lifshitz equation, which allows to determine the important scaling behavior with laser intensity and particle energy. We then move on to full scale 3D Monte Carlo simulations to verify the effect is observable with present day laser technology. We investigate the opportunities for an indirect observation of laser depletion in such side scattering scenarios.

Laser-driven neutron sources are routinely produced by the interaction of laser-accelerated protons with a converter. They present complementary characteristics to those of conventional accelerator-based neutron sources (e.g. short pulse durations, enabling novel applications like radiography). We present here results from an experiment aimed at performing a global characterization of the neutrons produced using the Titan laser at the Jupiter Laser Facility (Livermore, USA), where protons were accelerated from 23 m thick plastic targets and directed onto a LiF converter to produce neutrons. For this purpose, several diagnostics were used to measure these neutron emissions, such as CR-39, activation foils, time-of-flight detectors and direct measurement of residual activity in the LiF converters. The use of these different, independently operating diagnostics enables comparison of the various measurements performed to provide a robust characterization. These measurements led to a neutron yield of neutrons per shot with a modest angular dependence, close to that simulated.

Experimental studies of the dynamics of shear flow and turbulence spreading at the edge of tokamak plasmas are reported. Scans of line-averaged density and plasma current are carried out while approaching the Greenwald density limit on the J-TEXT tokamak. In all scans, when the Greenwald fraction f G = n ¯ / n G = n ¯ / ( I p / π a 2 ) increases, a common feature of enhanced turbulence spreading and edge cooling is found. The result suggests that turbulence spreading is a good indicator of edge cooling, indeed better than turbulent particle transport is. The normalized turbulence spreading power increases significantly when the normalized E × B shearing rate decreases. This indicates that turbulence spreading becomes prominent when the shearing rate is weaker than the turbulence scattering rate. The asymmetry between positive/negative (blobs/holes) spreading events, turbulence spreading power and shear flow are discussed. These results elucidate the important effects of interaction between shear flow and turbulence spreading on plasma edge cooling.

In the context of condensed matter physics, the Sommerfeld parameter describes the enhancement or suppression of free-carrier charge density in the vicinity of a charged center. The Sommerfeld parameter is known for three-dimensional systems and is integral to the description of trap-assisted recombination in solids. Here we derive the Sommerfeld parameter in one and two dimensions and compare with the results in three dimensions. We provide an approximate analytical expression for the Sommerfeld parameter in two dimensions. Our results indicate that the effect of the Sommerfeld parameter is to suppress trap-assisted recombination in decreased dimensionality.

Antiferromagnets are promising platforms for transduction and transmission of quantum information via magnons—the quanta of spin waves—and they offer advantages over ferromagnets in regard to dissipation, speed of response and robustness to external fields. Recently, transduction was shown in a van der Waals antiferromagnet, where strong spin-exciton coupling enables readout of the amplitude and phase of coherent magnons by photons of visible light. This discovery shifts the focus of research to transmission, specifically to exploring the non-local interactions that enable magnon wave packets to propagate. Here we demonstrate that magnon propagation is mediated by long-range dipole–dipole interaction. This coupling is an inevitable consequence of fundamental electrodynamics and, as such, will likely mediate the propagation of spin at long wavelengths in the entire class of van der Waals magnets currently under investigation. Successfully identifying the mechanism of spin propagation provides a set of optimization rules, as well as caveats, that are essential for any future applications of these promising systems.

Direct laser acceleration of electrons during a high-energy, picosecond laser interaction with an underdense plasma has been demonstrated to be substantially enhanced by controlling the laser focusing geometry. Experiments using the OMEGA EP facility measured electrons accelerated to maximum energies exceeding 120 times the ponderomotive energy under certain laser focusing, pulse energy, and plasma density conditions. Two-dimensional particle-in-cell simulations show that the laser focusing conditions alter the laser field evolution, channel fields generation, and electron oscillation, all of which contribute to the final electron energies. The optimal laser focusing condition occurs when the transverse oscillation amplitude of the accelerated electron in the channel fields matches the laser beam width, resulting in efficient energy gain. Through this observation, a simple model was developed to calculate the optimal laser focal spot size in more general conditions and is validated by experimental data.

Anomalous transport of multi-species plasma is considered with the generalized Hasegawa–Wakatani model [A. R. Knyazev and S. I. Krasheninnikov, Phys. Plasmas 31, 012502 (2024)] further extended to incorporate the Finite Larmor Radius (FLR) effects. By introducing the “associated” enstrophy, it is shown that with no FLR effects (where anomalous transport of all ion species is described as a transport of passive scalars in the turbulent fields of the electrostatic potential and electron density fluctuations) the fluctuating densities of ion species converge to the state where they are linearly proportional to electron density and vorticity fluctuations, which confirm previous numerical findings of [A. R. Knyazev and S. I. Krasheninnikov, Phys. Plasmas 31, 012502 (2024)]. However, in contrast to the “cold” ion approximation, with the FLR effects included, both the plasma turbulence and the dynamics of all ion species become interconnected. Therefore, for simplicity, the FLR effects in this work were considered only for a small “trace” impurity fraction. It is found that for light (neon) “trace” impurity, the FLR effects reduce both anomalous flux and density fluctuations. However, for heavy (tungsten) “trace” impurity, the FLR effects exhibit non-monotonic impact on anomalous transport.

This paper investigates the transport of drugs delivered by direct injection into the cerebrospinal fluid (CSF) that fills the intrathecal space surrounding the spinal cord. Because of the small drug diffusivity, the dispersion of neutrally buoyant drugs has been shown in previous work to rely mainly on the mean Lagrangian flow associated with the CSF oscillatory motion. Attention is given here to effects of buoyancy, arising when the drug density differs from the CSF density. For the typical density differences found in applications, the associated Richardson number is shown to be of order unity, so that the Lagrangian drift includes a buoyancy-induced component that depends on the spatial distribution of the drug, resulting in a slowly evolving cycle-averaged flow problem that can be analysed with two-time scale methods. The asymptotic analysis leads to a nonlinear integro-differential equation for the spatiotemporal solute evolution that describes accurately drug dispersion at a fraction of the cost involved in direct numerical simulations of the oscillatory flow. The model equation is used to predict drug dispersion of positively and negatively buoyant drugs in an anatomically correct spinal canal, with separate attention given to drug delivery via bolus injection and constant infusion.

Abstract: The self-organization is one of the most interesting phenomena in the non-equilibrium complex system, generating ordered structures of different sizes and durations. In tokamak plasmas, various self-organized phenomena have been reported, and two of them, coexisting in the near-marginal (interaction dominant) regime, are avalanches and the E × B staircase. Avalanches mean the ballistic flux propagation event through successive interactions as it propagates, and the E × B staircase means a globally ordered pattern of self-organized zonal flow layers. Various models have been suggested to understand their characteristics and relation, but experimental researches have been mostly limited to the demonstration of their existence. Here we report detailed analyses of their dynamics and statistics and explain their relation. Avalanches influence the formation and the width distribution of the E × B staircase, while the E × B staircase confines avalanches within its mesoscopic width until dissipated or penetrated. Our perspective to consider them the self-organization phenomena enhances our fundamental understanding of them as well as links our findings with the self-organization of mesoscopic structures in various complex systems.

Abstract: We report on comprehensive experimental studies of turbulence spreading in edge plasmas. These studies demonstrate the relation of turbulence spreading and entrainment to intermittent convective density fluctuation events or bursts (i.e. blobs and holes). The non-diffusive character of turbulence spreading is thus elucidated. The turbulence spreading velocity (or mean jet velocity) manifests a linear correlation with the skewness of density fluctuations, and increases with the auto-correlation time of density fluctuations. Turbulence spreading by positive density fluctuations is outward, while spreading by negative density fluctuations is inward. The degree of symmetry breaking between outward propagating blobs and inward propagating holes increases with the amplitude of density fluctuations. Thus, blob-hole asymmetry emerges as crucial to turbulence spreading. These results highlight the important role of intermittent convective events in conveying the spreading of turbulence, and constitute a fundamental challenge to existing diffusive models of spreading.

The Exabiome project seeks to improve the understanding of microbiomes through the development of methods for accelerating metagenomic science using exascale computing. This article gives an overview of scientific impact of the three components of the project: metagenome assembly, protein family detection, and comparative analysis of metagenomes. Exabiome developed MetaHipMer, the only metagenome assembler capable of scaling to full exascale systems. MetaHipMer has enabled ground-breaking assemblies on the Frontier supercomputer, with many scientific benefits, such as the discovery of rare species and viral genomes. To investigate protein families, Exabiome developed two exascale tools, PASTIS and HipMCL. Together, these can utilize exascale resources to understand the functional diversity of billions of dark matter proteins and novel protein families. For comparative analysis, Exabiome developed kmerprof, a tool that can be used to compare huge metagenomes for many different scientific purposes, for example, grouping human microbiomes according to body location.

Earthquakes present worldwide risk to economic and human safety. The 2023 earthquakes in Turkiye provided a reminder of the potential for catastrophic consequences with 50,700 deaths and 15.7 million people affected. The ability to predict ground motions and infrastructure damage for earthquakes continues to be a challenging problem for scientists and engineers. Until now, estimates of ground motions have been performed empirically by looking at sparse data from past earthquakes. This approach can provide statistical information on intensity amplitudes but cannot inform site-specific ground motions essential to developing the most effective resilience. Interest has grown in large-scale computational models to simulate earthquakes at regional scale. The U.S. Department of Energy EarthQuake SIMulation (EQSIM) framework was developed for regional-scale earthquake simulations at unprecedented fidelity, taking advantage of emerging GPU-accelerated systems. This article describes the EQSIM workflow and demonstrates regional-scale simulations with the new computational capability available to scientists in their quest to mitigate future disasters.

Models have historically represented fractured porous media with continuum descriptions that characterize the media using bulk parameters. The impact of small-scale features is not captured in these models, although they may be controlling the performance of subsurface applications. Pore-scale models can simulate processes in small-scale features by representing the pore space geometry explicitly but are computationally expensive for large domains. The alternative multiscale approach entails the combination of pore-scale and continuum-scale descriptions in a single framework. We use Chombo-Crunch, a computational capability that discretizes complex geometries with an adaptive embedded boundary method to contrast these two approaches. Chombo-Crunch takes advantage of recent computational performance and memory bandwidth improvements resulting from the emergence of exascale computing resources. These combined improvements enable the efficient simulation of reactive transport in fractured media with a high degree of fidelity and the ability to capture the control small-scale processes exert on the overall medium evolution.

Pressurized water reactors are designed to operate in a single-phase flow. However, during a flow loss or other off-design conditions liquid temperature may exceed saturation temperature and, if a continuous film of gas forms, a boiling casualty may result. If superhydrophobic surfaces are introduced among the fuel cell assembly, vapor bubbles show an affinity to these surfaces and gas may coalesce to escape faster, resulting in a larger margin to reach critical heat flux. In the present study, we consider air and liquid water mixture examining the overall flow dynamics in a case with no bulk liquid flow, reminiscent of a case with coolant pump failure. The parameter ranges examined span Reynolds number based on gas superficial velocity up to 5300 and produced bubble with bond numbers from 1 to 270, and Weber numbers based on estimated terminal velocity from 0.6 to 150. The Morton number was fixed at 2.63×10-11. As the flow with high gas volume fraction becomes optically opaque and is sensitive to intrusive instrumentation, a custom-built photon-counting dual energy threshold X-ray computed tomography system is employed for measurement of the time average void fraction within the bundle non-intrusively. Two dual plane wire mesh sensors upstream and downstream of the rod bundle are employed to obtain comparison void fraction, velocity profiles and bubble size distributions. Additionally, traditional pressure-based gas holdup measurements are employed to calculate time- and volume-averaged void fraction. The data show that the presence of superhydrophobic surface being present in the rod bundle results in a significantly lower gas volume fraction and higher localized void fraction at the superhydrophobic coating, as compared to those in a similar rod bundle without superhydrophobic internals. Graphical abstract: (Figure presented.)

The work presents numerical simulation studies of the role that dynamic plasma recycling on the main wall and divertor target surfaces plays in transient edge plasma transport phenomena, such as edge localized modes (ELMs). The studies are performed by coupling the edge plasma transport code UEDGE [Rognlien et al., J. Nucl. Mater. 196–198, 347 (1992)] and the wall reaction–diffusion transport code FACE [Smirnov et al., Fusion Sci. Technol. 71, 75 (2017)]. The two-dimensional, time-dependent, two-way coupling of the codes, in a realistic tokamak geometry, is accomplished using the Integrated Plasma Simulator framework [Elwasif et al., in 18th Euromicro Conference on Parallel, Distributed and Network-Based Processing (PDP 2010), Pisa, Italy (IEEE, 2010), pp. 419–427] for all modeled material plasma boundaries. The simulations show that dynamic plasma recycling has substantially different characteristics on the main wall and on the divertor plates. It is demonstrated that during an ELM cycle the outer wall can dynamically absorb and release a number of particles comparable to that expelled by the ELM from the core plasma, by far exceeding the dynamic retention capacity of the divertor surfaces. The resulting evolution of the edge and divertor plasma conditions during an ELM cycle is analyzed.

Abstract: In this work, we propose a reduced model with a dynamical critical gradient to study the formation, propagation, and conversion of the transport barriers. In contrast to the commonly adopted static critical gradient, an evolving critical gradient self-consistently softens the profile stiffness, so as to facilitate the generation of transport barriers. This is especially crucial to the internal transport barrier (ITB) formation. Numerically, we show that the inhomogeneity of turbulent and neoclassical diffusivities can induce the global wave front propagation of the transport barrier. When the heating power ramps quickly, the ITB propagates unidirectionally to the edge region and converts into an edge transport barrier. For slow power ramping, the propagating ITB will bifurcate into bidirectional wavefronts and finally convert into a steady double transport barrier state. Our model uncovers the vital role of a dynamical ‘profile-stiffness’ in depicting the global dynamics of the transport barrier.

The introduction of localized electronic states into a metal can alter its physical properties, for example enabling exotic metal physics including heavy fermion and strange metal behaviour. A common source of localized states in such systems are partially filled 4f and 5f shells because of the inherently compact nature of those orbitals. The interaction of electrons in these orbitals with the conduction sea is well described by the Kondo framework. However, there have also been observations of Kondo-like behaviour in 3d transition metal oxides and in 4d- and 5d-containing van der Waals heterostructures. This calls for a broader consideration of the physical requirements for Kondo systems. Here we show transport and thermodynamic hallmarks of heavy fermion and strange metal behaviour that arise in the kagome metal Ni3In, wherein the source of localized states is destructive interference-induced band flattening of partially filled Ni 3d states. With magnetic field and pressure tuning, we also find evidence that the system is proximate to quantum criticality, extending the analogy to f-electron Kondo lattices. These observations highlight the role of hopping frustration in metallic systems as a potential source for strong correlations. Additionally, this suggests a lattice-driven approach to realizing correlated metals with non-trivial band topology.

The concept and theory of potential vorticity in drift wave turbulence are extended to the case of an inhomogeneous magnetic field. A one-field magnetic potential vorticity conserving equation is derived via the use of conservative gyrokinetics. The similarity between the corresponding systems for drift wave turbulence and shallow water theory is discussed in detail. Zonal flow physics in an inhomogeneous magnetic field is discussed. In particular, a Charney–Drazin type nonacceleration theorem is derived from the novel system, which conserves magnetic potential vorticity. Extensions of the turbulent equipartition theory to the transport of magnetic potential vorticity are proposed.

Hydrodynamic equations for a one-component plasma are derived as a unification of the Euler equations with long-range Coulomb interaction. By using a variational principle, these equations self-consistently unify thermodynamics, dispersion laws, nonlinear motion, and conservation laws. In the moderate and strong coupling limits, it is argued that these equations work down to the length scale of the interparticle spacing. The use of a variational principle also ensures that closure is achieved self-consistently. Hydrodynamic equations are evaluated in both the Eulerian frame, where the fluid variables depend on the position in the laboratory, and the Lagrangian frame, where they depend on the position in some reference state, such as the initial position. Each frame has its advantages and our final theory combines elements of both. The properties of longitudinal and transverse dispersion laws are calculated for the hydrodynamic equations. A simple step function approximation for the pair distribution function enables simple calculations that reveal the structure of the equations of motion. The obtained dispersion laws are compared to molecular dynamics simulations and the theory of quasilocalized charge approximation. The action, which gives excellent agreement for both longitudinal and transverse dispersion laws for a wide range of coupling strengths, is elucidated. Agreement with numerical experiments shows that such a hydrodynamic approach can be used to accurately describe a one-component plasma at very small length scales comparable to the average interparticle spacing. The validity of this approach suggests considering nonlinear flows and other systems with long-range interactions in the future.

The family of transition-metal dipnictides has been of theoretical and experimental interest because this family hosts topological states and extremely large magnetoresistance (MR). Recently,TaAs2, a member of this family, has been predicted to support a topological crystalline insulating state. Here, by using high-resolution angle-resolved photoemission spectroscopy (ARPES), we reveal both closed and open pockets in the metallic Fermi surface (FS) and linearly dispersive bands on the (2‾01) surface, along with the presence of extreme MR observed from magneto-transport measurements. A comparison of the ARPES results with first-principles computations shows that the linearly dispersive bands on the measured surface ofTaAs2are trivial bulk bands. The absence of symmetry-protected surface state on the (2‾01) surface indicates its topologically dark nature. The presence of open FS features suggests that the open-orbit fermiology could contribute to the extremely large MR ofTaAs2.

The trade-off between robustness and tunability is a central challenge in the pursuit of quantum simulation and fault-tolerant quantum computation. In particular, quantum architectures are often designed to achieve high coherence at the expense of tunability. Many current qubit designs have fixed energy levels and consequently limited types of controllable interactions. Here by adiabatically transforming fixed-frequency superconducting circuits into modifiable Floquet qubits, we demonstrate an XXZ Heisenberg interaction with fully adjustable anisotropy. This interaction model can act as the primitive for an expressive set of quantum operations, but is also the basis for quantum simulations of spin systems. To illustrate the robustness and versatility of our Floquet protocol, we tailor the Heisenberg Hamiltonian and implement two-qubit iSWAP, CZ and SWAP gates with good estimated fidelities. In addition, we implement a Heisenberg interaction between higher energy levels and employ it to construct a three-qubit CCZ gate, also with a competitive fidelity. Our protocol applies to multiple fixed-frequency high-coherence platforms, providing a collection of interactions for high-performance quantum information processing. It also establishes the potential of the Floquet framework as a tool for exploring quantum electrodynamics and optimal control.

BOUT++ turbulence simulations were performed to investigate the impact of turbulence spreading on the edge localized mode (ELM) size and divertor heat flux width $({\lambda _q})$ broadening in small ELM regimes. This study is motivated by EAST experiments. BOUT++ linear simulations of a pedestal radial electric field (Er) scan show that the dominant toroidal number mode (n) shifts from high-n to low-n, with a narrow mode spectrum, and the maximum linear growth rate increases as the pedestal Er well deepens. The nonlinear simulations show that as the net E × B pedestal flow increases, the pressure fluctuation level and its inward penetration beyond the top of the pedestal both increase. This leads to a transition from small ELMs to large ELMs. Both inward and outward turbulence spreading are sensitive to the scrape-off-layer (SOL) plasma profiles. The inward turbulence spreading increases for the steep SOL profiles, leading to increasing pedestal energy loss in the small ELM regime. The SOL width $({\lambda _q})$ is significantly broadened progressing from the ELM-free to small ELM regime, due to the onset of strong radial turbulent transport. The extent of the SOL width $({\lambda _q})$ broadening depends strongly on outward turbulence spreading. The fluctuation energy intensity flux ${\varGamma _\varepsilon }$ at the separatrix can be enhanced by increasing either pedestal Er flow shear or local SOL pressure gradient. The ${\lambda _q}$ is broadened as the fluctuation energy intensity flux ${\varGamma _\varepsilon }$ at the last close flux surface (LCFS) increases. Local SOL E × B flow shear will restrain outward turbulence spreading and the associated heat flux width broadening. Operating in H-mode with small ELMs has the potential to solve two critical problems: reducing the ELM size and broadening the SOL width.

Interactions between fluid flow and elastic structures are important in many naturally occurring and engineered systems. This review collects and organizes recent theoretical and experimental developments in understanding fluid-structure interactions at low Reynolds numbers. Particular attention is given to the motion of objects moving in close proximity to deformable soft materials and the ensuing interplay between fluid flow and elastic deformation. We discuss how this interplay can be understood in terms of forces and torques, and harnessed in applications such as microrheometry, tribology, and soft robotics. We then discuss the interaction of soft and wet objects close to contact, where intermolecular forces and surface roughness effects become important and are sources of complexity and opportunity.

Abstract: We present a multi-scale model of quasi-mode evolution in a stochastic magnetic field. The similarity between a quasi-mode and a ballooning mode enables us to address the challenges arising from the disparate geometries in the theories of ballooning modes in the presence of resonant magnetic perturbations. We obtain useful insights into our understanding of ballooning mode dynamics in a stochastic background. To maintain quasi-neutrality at all scales, the beat between the quasi-mode and the stochastic magnetic field drives microturbulence, which drives the turbulent background that promotes mixing and damps the quasi-mode. As a result of the broad mode structure of the quasi-mode, the turbulent viscosity and the turbulent diffusivity produced by the microturbulence are larger than those in our related study on resistive interchange modes. The stochastic magnetic field can also enhance the effective plasma inertia and reduce the effective drive, thereby slowing the mode growth. A nontrivial correlation between the microturbulence and the magnetic perturbations is shown to develop. This could account for the reduction in the Jensen-Shannon complexity of pedestal turbulence in the RMP ELM suppression phase observed in recent experiments. Directions for future experimental and theoretical studies are suggested.

A simple two-dimensional fluid–structure interaction problem, involving viscous oscillatory flow in a channel separated by an elastic membrane from a fluid-filled slender cavity, is analysed to shed light on the flow dynamics pertaining to syringomyelia, a neurological disorder characterized by the appearance of a large tubular cavity (syrinx) within the spinal cord. The focus is on configurations in which the velocity induced in the cavity, representing the syrinx, is comparable to that found in the channel, representing the subarachnoid space surrounding the spinal cord, both flows being coupled through a linear elastic equation describing the membrane deformation. An asymptotic analysis for small stroke lengths leads to closed-form expressions for the leading-order oscillatory flow, and also for the stationary flow associated with the first-order corrections, the latter involving a steady distribution of transmembrane pressure. The magnitude of the induced flow is found to depend strongly on the frequency, with the result that for channel flow rates of non-sinusoidal waveform, as those found in the spinal canal, higher harmonics can dominate the sloshing motion in the cavity, in agreement with previous in vivo observations. Under some conditions, the cycle-averaged transmembrane pressure, also showing a marked dependence on the frequency, changes sign on increasing the cavity transverse dimension (i.e. orthogonal to the cord axis), underscoring the importance of cavity size in connection with the underlying hydrodynamics. The analytic results presented here can be instrumental in guiding future numerical investigations, needed to clarify the pathogenesis of syringomyelia cavities.

Elucidating the phase diagram of solid hydrogen is a key objective in condensed matter physics. Several decades ago, it was proposed that at low temperatures and high pressures, solid hydrogen would be a metal with a high superconducting transition temperature. This transition to a metallic state can happen through the closing of the energy gap in the molecular solid or through a transition to an atomic solid. Recent experiments have managed to reach pressures in the range of 400-500 GPa, providing valuable insights. There is strong evidence suggesting that metallization via either of these mechanisms occurs within this pressure range. Computational and experimental studies have identified multiple promising crystal phases, but the limited accuracy of calculations and the limited capabilities of experiments prevent us from determining unequivocally the observed phase or phases. Therefore, it is crucial to investigate the superconducting properties of all the candidate phases. Recently, we reported the superconducting properties of theC2/c-24,Cmca-12,Cmca-4 andI41/amd-2 phases, including anharmonic effects. Here, we report the effects of anisotropy on superconducting properties using Eliashberg theory. Then, we investigate the superconducting properties of deuterium and estimate the size of the isotope effect for each phase. We find that the isotope effect on superconductivity is diminished by anharmonicity in theC2/c-24 andCmca-12 phases and enlarged in theCmca-4 andI41/amd-2 phases. Our anharmonic calculations of theC2/c-24 phase of deuterium agree closely with the most recent experiment by Loubeyreet al(2022Phys. Rev. Lett.29035501), indicating that theC2/c-24 phase remains the leading candidate in this pressure range, and has a strong anharmonic character. These characteristics can serve to distinguish among crystal phases in experiment. Furthermore, expanding our understanding of superconductivity in pure hydrogen holds significance in the study of high-Tchydrides.

Exascale supercomputers are first-of-their-kind instruments with potentially paradigm-shifting capabilities. However, enabling complex applications at scale and high performance can be difficult, requiring hard work and deep technical understanding. Specific programming models are often needed to achieve the desired performance, and, early on, systems are undergoing debugging and rapid development. This makes it challenging to run on new systems quickly following their deployment. The Application Integration (AppInt) area of the Exascale Computing Project was designed to be the integration point between applications, supporting software, system environments, high-performance computing facilities, and vendors. We describe how the AppInt team addressed these challenges while also promoting the use of portable and sustainable programming models and helping harden systems prior to general availability.

The Exascale Computing Project software deployment effort developed and advanced DevOps capabilities. One goal was to enable robust continuous integration (CI) workflows that span the protected high-performance computing environments found within many of the U.S. Department of Energy's (DOE's) national laboratories. This article highlights several challenges encountered with enabling automation, such as charging models for CI jobs and meeting individualized security requirements that revolve around strongly associating running code with a human identity. It also describes how the Jacamar CI tool evolved to meet later requirements and became a key aspect of the solutions currently offered. Derived from this experience, we offer a conceptual framework for understanding current and future CI challenges at DOE facilities and offer suggestions for long-term solutions.

Computational and data-enabled science and engineering are revolutionizing advances throughout science and society, at all scales of computing. For example, teams in the U.S. Department of Energy's Exascale Computing Project have been tackling new frontiers in modeling, simulation, and analysis by exploiting unprecedented exascale computing capabilities-building an advanced software ecosystem that supports next-generation applications and addresses disruptive changes in computer architectures. However, concerns are growing about the productivity of the developers of scientific software. Members of the Interoperable Design of Extreme-scale Application Software project serve as catalysts to address these challenges through fostering software communities, incubating and curating methodologies and resources, and disseminating knowledge to advance developer productivity and software sustainability. This article discusses how these synergistic activities are advancing scientific discovery-mitigating technical risks by building a firmer foundation for reproducible, sustainable science at all scales of computing, from laptops to clusters to exascale and beyond.

We consider a constructive definition of the multivariate Pareto that factorizes the random vector into a radial component and an independent angular component. The former follows a univariate Pareto distribution, and the latter is defined on the surface of the positive orthant of the infinity norm unit hypercube. We propose a method for inferring the distribution of the angular component by identifying its support as the limit of the positive orthant of the unit p-norm spheres and introduce a projected gamma family of distributions defined through the normalization of a vector of independent random gammas to the space. This serves to construct a flexible family of distributions obtained as a Dirichlet process mixture of projected gammas. For model assessment, we discuss scoring methods appropriate to distributions on the unit hypercube. In particular, working with the energy score criterion, we develop a kernel metric that produces a proper scoring rule and presents a simulation study to compare different modeling choices using the proposed metric. Using our approach, we describe the dependence structure of extreme values in the integrated vapor transport (IVT), data describing the flow of atmospheric moisture along the coast of California. We find clear but heterogeneous geographical dependence.

The U.S. Exascale Computing Project (ECP) has succeeded in preparing applications to run efficiently on the first reported exascale supercomputers in the world. To achieve this, it modernized the whole leadership software stack, from libraries to simulation codes. In this article, we contrast selected leadership software before and after the ECP. We discuss how sustainable research software development for leadership computing can embrace the conversation with the hardware vendors, leadership computing facilities, software community, and domain scientists who are the application developers and integrators of software products. We elaborate on how software needs to take portability as a central design principle and to benefit from interdependent teams; we also demonstrate how moving to programming languages with high momentum, like modern C++, can help improve the sustainability, interoperability, and performance of research software. Finally, we showcase how cross-institutional efforts can enable algorithm advances that are beyond incremental performance optimization.