Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Theoretical plasma physics”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Measurement of helicon waves with phase contrast imaging on DIII-D – A theoretical feasibility study

A DIII-D high-beta H-mode discharge, with I p = 850 kA, B t = 2.1 T n e = 4 10 19 m –3 , has been designed to validate full wave modeling of helicon waves by optimizing the expected response of the Phase Contrast Imaging diagnostic. Helicon waves have been predicted to have high current drive efficiency off-axis without facing the accessibility issues of lower-hybrid waves. To test these predictions experimentally, DIII-D has recently commissioned a high-power helicon antenna. To confidently predict the behavior of helicon waves in future devices, measurements of their fundamental properties and validation against models will be essential. Phase contrast imaging (PCI) is an absolutely calibrated internal reference interferometer able to measure density fluctuations with radial wavenumbers k R between 1.5 cm –1 and 20 cm –1 . For helicon waves 2 cm –1 < k R < 10 cm –1 is expected, allowing PCI to measure their envelope and wavenumber spectrum. This makes PCI a powerful tool for the validation of state-of-the-art models, like the AORSA full wave code. AORSA is used to compute the density perturbations measured by the PCI with 2D calculations corresponding to 11 different toroidal mode numbers combined to resolve the trajectory of the helicon wave in 3D. This is necessary because the waves travel 120 degrees toroidally from the antenna to the PCI. Here, a cold plasma finite element model (CPFEM) [4] is used to predict propagation through the scrape-off layer. The result of the CPFEM model is connected to AORSA by creating an artificial Gaussian antenna on the last closed flux surface. PCI shows best results for waves with small vertical wavenumbers k z . Modeling the helicon waves for several past DIII-D experiments shows that k z is minimized if the intersection of the helicon and the PCI laser beams occurs in the midplane. For such an optimized scenario the predicted signal level is two orders of magnitude larger than the background density fluctuations arising from broadband turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Center for Tokamak Transients Simulations

The major goals for the Center for Tokamak Transients Simulations (CTTS) were to develop detailed numerical simulations of disruptive transients in tokamak experiments and to use them in theoretical studies of the thermal-quench (TQ) and current-quench (CQ) phases of disruption. The University of Wisconsin-Madison component of CTTS contributed relevant development for the NIMROD code (https://nimrodteam.org) and applied it to disruptions that involve global vertical displacement. Code-development contributions include boundary conditions that represent sheath effects which occur when tokamak plasma is in contact with the surrounding vessel wall during vertical displacement events (VDEs). We also implemented a reduced model of energetic runaway electrons (REs) to simulate the effect of these energetic electrons on the CQ phase of tokamak disruptions. Our computations of VDEs include 2D and 3D verification through benchmarks with two other widely used macroscale plasma simulation codes, M3D-C1 and JOREK. Over nonlinear evolution to the final termination of plasma current, the 2D simulation results on magnetic axis location, toroidal current, and halo current track each other well. With the 3D benchmark computation, all three codes found the same qualitative behavior of the asymmetric instabilities having increasing growth rate after contact with the wall, followed by nonlinear excitation of other toroidal harmonics. Our NIMROD computations for an idealized configuration show agreement with an analytical assessment by V. D. Pustovitov [Nucl. Fusion 55, 113032 (2015)] that the plasma and non-ideal wall remain in force-balance such that net forces only result from magnetic stresses over the outer surface of the wall. Our study of the influence of boundary conditions on VDEs shows that extended-MHD simulations of VDEs are sensitive to electron energy transport to the wall and that boundary conditions on flow are only significant to the extent that they influence the electron energy transport. Resistive-MHD NIMROD simulations of the tokamak discharges in the Madison Symmetric Torus show that magnetic perturbations of poloidal wavenumber m=3 lead to chaotic magnetic topology over the edge region, which helps deconfine energetic electrons. A second area of RE study shows that the coupled resistive MHD/reduced RE model reproduces the resistive-hose particle-beam instability. With respect to educational opportunities, this effort involved and supported four graduate students at the University of Wisconsin-Madison and one postdoctoral associate. The results of this study have been disseminated through journal publications, conference presentations, technical reports, and PhD dissertations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Focusing of laser-accelerated ions with spherically curved targets

A theoretical model for laser-accelerated ion focusing from spherically curved targets using a Gaussian self-similar solution is presented. This model describes the evolution of the focal location and focal radius with final ion energy and the initial radial boundary of the ions, the latter being a function of both target geometry and the ion acceleration radius for a given ion energy. The theory is supported by particle-in-cell simulations of a variety of target shapes with varying radii of curvature and target opening angles, as well as variations in the injected electron beam radius and energy spectrum. The theory and simulations suggest that the focal location varies linearly with the radius of curvature, with the square root of the ratio of ion energy to effective electron temperature, and monotonically increases with the energy-dependent ratio of the initial ion radial boundary to the radius of curvature. Considering ponderomotive acceleration, this suggests that the focal length should scale inversely with intensity, $(I_Lλ_L^2)$ −1/4 , suggesting a 10× increase in intensity will reduce the focal distance by ̃1.8× for a given ion energy.

Physics - Plasma physics↗
DOE OSTI2025

Numerical simulation of compressible fluid-dynamics in the chamber of inertial fusion energy systems

Here, this paper aims to establish new and innovative modeling capabilities for analyzing chambers in Inertial Fusion Energy (IFE) systems. IFE is emerging as a promising method to achieve fusion power production, but several challenges must be overcome to develop an IFE pilot plant or deploy commercial IFE systems. These challenges are both theoretical and technical, encompassing a deeper understanding of the underlying physical phenomena and the development of new technologies and materials. One of the needs is to develop mathematical models to describe IFE systems and numerical tools to simulate them. This paper contributes to this endeavor by presenting a new OpenFOAM solver for IFE systems, focusing on gas dynamics in their chambers. The analysis and development of chamber designs will play a significant role in the transition from single-shot experiments to high-repetition rates, as there is a need to protect the chamber walls from the intense radiation fields produced by fusion reactions. A promising design option, normally referred to as thick wall chamber design, consists in using lithium or molten salt jet arrays within the chamber. A critical phenomenon is the venting of high-pressure gases from the center to the external part of the chamber, passing through the blanket array. This process involves the propagation and attenuation of strong pressure waves, requiring suitable modeling approaches for compressible fluid-dynamics. The solver proposed in this work implements a multi-material hydrodynamics model tailored to accurately describe the non-linear propagation of pressure waves while avoiding numerical oscillation issues typical of high-velocity compressible simulation. This solver is verified against numerical test cases, validated against experimental data, and applied to the analysis of the High-Yield Lithium-Injection Fusion-Energy (HYLIFE-I) concept. The relevance of this paper is threefold. Firstly, it contributes to developing and testing modeling approaches for compressible fluid-dynamics phenomena, with specific focus on the new and unexplored topic of IFE thick-liquid-wall blanket modeling. Secondly, it marks one of the first applications of the OpenFOAM library in the research field of IFE systems. Finally, the investigated problem is of practical interest for IFE developers, as it provides useful indications about relevant phenomena in pressure wave propagation in the chamber of these systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Medium-induced modification of groomed and ungroomed jet mass and angularities in Pb–Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The ALICE Collaboration presents a new suite of jet substructure measurements in Pb–Pb and pp collisions at a center-of-mass energy per nucleon pair $\sqrt{s_{NN}}$ = 5.02 TeV. These measurements provide access to the internal structure of jets via the momentum and angle of their constituents, probing how the quark–gluon plasma modifies jets, an effect known as jet quenching. Jet grooming additionally removes soft wide-angle radiation to enhance perturbative accuracy and reduce experimental uncertainties. We report the groomed and ungroomed jet mass m jet and jet angularities λ$^{κ}_{α}$ using κ = 1 and α > 0. Charged-particle jets are reconstructed at midrapidity using the anti-k T algorithm with resolution parameter R = 0.2. A narrowing of the jet mass and angularity distributions in Pb–Pb collisions with respect to pp is observed and is enhanced for groomed results, confirming modification of the jet core. By using consistent jet definitions and kinematic cuts between the mass and angularities for the first time, previous inconsistencies in the interpretation of quenching measurements are resolved, rectifying a hurdle for understanding how jet quenching arises from first principles and highlighting the importance of a well-controlled baseline. These results are compared with a variety of theoretical models of jet quenching, providing constraints on jet energy-loss mechanisms in the quark–gluon plasma.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Weibel-mediated filamentary structures observed in the ICF context

Here, in light of novel and past experimental results, we demonstrate how Weibel-mediated filamentary structures can develop in the expanding plasma plume of a laser-irradiated foil. The transverse ballistic cooling that occurs during the quasi-spherical plasma expansion naturally drives an electron pressure anisotropy, resulting in the growth of electron current filaments. This effect competes with electron–ion Coulomb collisions, which tend to isotropize the electron distribution function. Based on theoretical and particle-in-cell modeling, we provide estimates of the dominant wavelength and amplitude of the self-generated magnetic fluctuations, which are found to explain experimental data obtained at the OMEGA and Laser Megajoule facilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Observation of $t\bar{t}$ Production in Pb+Pb Collisions at $\sqrt{s_{NN}}$ = 5.02 TeV with the ATLAS Detector

Top-quark pair production is observed in lead–lead (Pb+Pb) collisions at $\sqrt{s_{NN}}$ = 5.02 TeV at the Large Hadron Collider with the ATLAS detector. The data sample was recorded in 2015 and 2018, amounting to an integrated luminosity of 1.9 nb −1 . Events with exactly one electron and one muon and at least two jets are selected. Top-quark pair production is measured with an observed (expected) significance of 5.0 (4.1) standard deviations. The measured top-quark pair production cross section is 𝜎 𝑡⁢$\tilde{i}$ = 3.6⁢$^{+1.0}_{−0.9}$⁢(stat)⁢ $^{+0.8}_{−0.5}$⁢(syst) μ⁢b, with a total relative uncertainty of 31%, and is consistent with theoretical predictions using a range of different nuclear parton distribution functions. The observation of this process consolidates the evidence of the existence of all quark flavors in the preequilibrium stage of the quark-gluon plasma at very high energy densities, similar to the conditions present in the early Universe.

Quark-gluon plasma↗
DOE OSTI2022

Tritium Transport Modeling at the Pore Scale in Ceramic Breeder Materials Using TMAP8

Fusion reactors depend on the blanket material to breed and release tritium at the same rate or faster than it is consumed by the fusion reaction. Cellular ceramic breeders (CCBs) are dense materials that can maintain a high tritium breeding ratio while promoting tritium release because of highly connected pores. Assessing the tritium breeding capabilities of these materials requires a combination of extensive experimental and modeling efforts. In this work, we develop and calibrate a multiphysics model of tritium transport. This novel model accounts for ceramic and pore diffusion, trapping and detrapping, and several surface reactions at the pore surface. We perform a sensitivity analysis and calibrate the model by comparing its predictions against experimental measurements of deuterium absorption. The calibrated model is then used to model tritium absorption in samples with different pore microstructures to investigate the effect of pore interconnectivity on tritium absorption. The model is part of the development of the multiscale, multiphysics framework for tritium transport [i.e., the Tritium Migration Analysis Program (TMAP8)], which is itself built on top of the finite-element multiphysics framework multiphysics object-oriented simulation environment (MOOSE). This study demonstrates some of TMAP8’s capabilities and is the first step toward assessing the tritium breeding capabilities of ceramic breeder material designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Measurement of electrons from beauty-hadron decays in $pp$ and Pb-Pb collisions at $\sqrt{s_{NN}}$ = $5.02$ TeV

The production of electrons from beauty-hadron decays was measured at midrapidity in proton-proton ($pp$) and central Pb-Pb collisions at center-of-mass energy per nucleon-nucleon pair $\sqrt{s_{NN}}$ = 5.02 TeV, using the ALICE detector at the LHC. The cross section measured in p p collisions in the transverse momentum interval $2$ < $p_T$ < $26$ $\mathrm{GeV}$ / $c$ was compared with models based on perturbative quantum chromodynamics calculations. The yield in the 10% most central Pb-Pb collisions, measured in the interval $2$ < $p_T$ < $26$ $\mathrm{GeV}$ / $c$, was used to compute the nuclear modification factor $R_{AA}$, extrapolating the $pp$ reference cross section to $p_T$ larger than 8 $\mathrm{GeV}$ / $c$. The measured $R_{AA}$ shows significant suppression of the yield of electrons from beauty-hadron decays at high $p_T$ and does not show a significant dependence above 8 $\mathrm{GeV}$ / $c$ within uncertainties. The results are described by several theoretical models based on different implementations of the interaction of heavy quarks with a quark-gluon plasma, which predict a smaller energy loss for beauty quarks compared to light and charm quarks.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Model-free Rayleigh weight from x-ray Thomson scattering measurements

X-ray Thomson scattering (XRTS) has emerged as a powerful tool for the diagnostics of matter under extreme conditions. In principle, it gives one access to important system parameters such as the temperature, density, and ionization state, but the interpretation of the measured XRTS intensity usually relies on theoretical models and approximations. In this context, a key property is given by the Rayleigh weight that describes the electronic localization around the ions. Here, we show that it is possible to extract the Rayleigh weight directly from the experimental data without the need for any model calculations or simulations. As a practical application, we consider an experimental measurement of strongly compressed Be at the National Ignition Facility [Döppner et al., Nature 618, 270–275 (2023)]. We demonstrate that experimental results for the Rayleigh weight open up new avenues for the interpretation of XRTS experiments by matching the measurement with ab initio simulations such as density functional theory or path integral Monte Carlo. Interestingly, this new procedure leads to significantly lower density compared to previously used chemical models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Measurement of muon pairs produced via γ γ scattering in nonultraperipheral Pb + Pb collisions at s N N = 5.02 TeV with the ATLAS detector

Results of a measurement of dimuon photoproduction in nonultraperipheral Pb + Pb collisions at √ s NN = 5.02 TeV are presented. The measurement uses ATLAS data from the 2015 and 2018 Pb + Pb data-taking periods at the LHC with an integrated luminosity of 1.94 nb –1 . The γγ → μ + μ – pairs are identified via selections on pair momentum asymmetry and acoplanarity. Differential cross sections for dimuon production are measured in different centrality, average muon momentum, and pair rapidity intervals as functions of acoplanarity and k ⊥ , the transverse momentum kick of one muon relative to the other. Measurements are also made as a function of the rapidity separation of the muons and the angle of the muon pair relative to the second-order event plane to test whether magnetic fields generated in the quark-gluon plasma affect the measured muons. A prior observation of a centrality-dependent broadening of the acoplanarity distribution is confirmed. Furthermore, the improved precision of the measurement reveals a depletion in the number of pairs having small acoplanarity or k ⊥ values in more central collisions. The acoplanarity distributions in a given centrality interval are observed to vary with the mean p T of the muons in the pair, but the k ⊥ distributions do not. Comparisons with recent theoretical predictions are made. The predicted trends associated with effects of magnetic fields on the dimuons are not observed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Soft-hard framework with exact four-momentum conservation for small systems

A new framework, called x-scape, for the combined study of both hard and soft transverse momentum sectors in high-energy proton-proton (𝑝−𝑝) and proton-nucleus (𝑝−𝐴) collisions is set up. A dynamical initial state is set up using the 3d-Glauber model with transverse locations of hotspots within each incoming nucleon. A hard scattering that emanates from two colliding hotspots is carried out using the Pythia generator. Initial state radiation from the incoming hard partons is carried out in a new module called I-matter, which includes the longitudinal location of initial splits. The energy-momentum of both the initial hard partons and their associated beam remnants is removed from the hot spots, depleting the energy-momentum available for the formation of the bulk medium. Outgoing showers are simulated using the matter generator, and results are presented for both cases, allowing for and not allowing for energy loss. First comparisons between this hard-soft model and single inclusive hadron and jet data from 𝑝−𝑝 and minimum bias 𝑝−Pb collisions are presented. Single hadron spectra in 𝑝−𝑝 are used to carry out a limited (in number of parameters) Bayesian calibration of the model. Fair comparisons with data are indicative of the utility of this new framework. Theoretical studies of the correlation between jet 𝑝 𝑇 and event activity at mid and forward rapidity are carried out.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Production of $Υ(nS)$ mesons in Pb + Pb and $pp$ collisions at 5.02 TeV

A measurement of the production of vector bottomonium states, $Υ$(1S), $Υ$(2S), and $Υ$(3S), in Pb + Pb and $pp$ collisions at a center-of-mass energy per nucleon pair of 5.02 TeV is presented. The data correspond to integrated luminosities of 1.38 nb –1 of Pb + Pb data collected in 2018, 0.44 nb –1 of Pb + Pb data collected in 2015, and 0.26 fb –1 of $pp$ data collected in 2017 by the ATLAS detector at the Large Hadron Collider. The measurements are performed in the dimuon decay channel for transverse momentum ${p}_{T}^{μμ}$ < 30 GeV, absolute rapidity |$y^{μμ}$| < 1.5, and Pb + Pb event centrality 0–80%. The production rates of the three bottomonium states in Pb + Pb collisions are compared with those in $pp$ collisions to extract the nuclear modification factors as functions of event centrality, ${p}_{T}^{μμ}$, and $|y^{μμ}|$. In addition, the suppression of the excited states relative to the ground state is studied. The results are compared with theoretical model calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Comparison of inclusive and photon-tagged jet suppression in 5.02 TeV Pb+Pb collisions with ATLAS

Parton energy loss in the quark–gluon plasma (QGP) is studied with a measurement of photon-tagged jet production in 1.7 nb –1 of Pb+Pb data and 260 pb –1 of pp data, both at $\sqrt{^S{\text{NN}}}$ = 5.02 TeV, with the ATLAS detector. The process $pp$ → γ +jet+X and its analogue in Pb+Pb collisions is measured in events containing an isolated photon with transverse momentum ($p_{\text{T}}$) above 50 GeV and reported as a function of jet $p_{\text{T}}$. This selection results in a sample of jets with a steeply falling $p_{\text{T}}$ distribution that are mostly initiated by the showering of quarks. The $pp$ and Pb+Pb measurements are used to report the nuclear modification factor, $R_{\text{AA}}$, and the fractional energy loss, $S_{\text{loss}}$, for photon-tagged jets. In addition, the results are compared with the analogous ones for inclusive jets, which have a significantly smaller quark-initiated fraction. The $R_{\text{AA}}$ and $S_{\text{loss}}$ values are found to be significantly different between those for photon-tagged jets and inclusive jets, demonstrating that energy loss in the QGP is sensitive to the colour-charge of the initiating parton. Finally, the results are also compared with a variety of theoretical models of colour-charge-dependent energy loss.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Energy dependence of polarized 𝛾⁢𝛾 → 𝑒 + ⁢ 𝑒− in peripheral Au+Au collisions at $\sqrt{s_{NN}}$ = 54.4 and 200 GeV with the STAR experiment at RHIC

We report the differential yields at mid-rapidity of the Breit-Wheeler process (𝛾⁢𝛾 → 𝑒 + ⁢𝑒 − ) in peripheral Au+Aucollisions at $\sqrt{s_{NN}}$ = 54.4 and 200 GeV with the STAR experiment at the Relativistic Heavy Ion Collider (RHIC), as a function of energy $\sqrt{s_{NN}}$, 𝑒 + ⁢𝑒 − transverse momentum 𝑝 T , 𝑝$^{2}_{T}$, invariant mass 𝑀 𝑒⁢𝑒 , and azimuthal angle. In the invariant mass range of 0.4 < 𝑀 𝑒⁢𝑒 < 2.6GeV/𝑐 2 at low transverse momentum (𝑝 T < 0.15GeV/𝑐), the yields increase while the pair √⟨𝑝$^{2}_{T}$⟩ decreases with increasing $\sqrt{s_{NN}}$, a feature that is correctly predicted by the QED calculation. Here, the energy dependencies of the measured quantities are sensitive to the nuclear form factor, infrared divergence and photon polarization. The data are compiled and used to extract the charge radius of the Au nucleus.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Single-Shot Reconstruction of Electron Beam Longitudinal Phase Space in a Laser Wakefield Accelerator

We report on a single-shot longitudinal phase-space reconstruction diagnostic for electron beams in a laser wakefield accelerator via the experimental observation of distinct periodic modulations in the angularly resolved spectra. Such modulated angular spectra arise as a result of the direct interaction between the ultrarelativistic electron beam and the laser driver in the presence of the wakefield. A constrained theoretical model for the coupled oscillator, assisted by a genetic algorithm, can recreate the experimental electron spectra and, thus, fully reconstructs the longitudinal phase-space distribution of the electron beam with a temporal resolution of approximately 1.3 fs. In particular, it reveals the slice energy spread of the electron beam, which is important to measure for applications such as x-ray free electron lasers. In our experiment, the root-mean-square energy spread retrieved is bounded at 9.9 MeV, corresponding to a 0.9%–3.0% relative spread, despite the overall GeV energy beam having approximately 100% relative energy spread.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Centrality dependence of Lévy-stable two-pion Bose-Einstein correlations in $\sqrt{𝑠{𝑁⁢𝑁}}$ = 200 GeV Au + Au collisions

The PHENIX experiment measured the centrality dependence of two-pion Bose-Einstein correlation functions in $\sqrt{𝑠{𝑁⁢𝑁}}$ = 200 GeV Au + Au collisions at the Relativistic Heavy Ion Collider at Brookhaven National Laboratory. The data are well represented by Lévy-stable source distributions. The extracted source parameters are the correlation-strength parameter 𝜆, the Lévy index of stability 𝛼, and the Lévy-scale parameter 𝑅 as a function of transverse mass 𝑚 𝑇 and centrality. The 𝜆⁡(𝑚 𝑇 ) parameter is constant at larger values of 𝑚 𝑇 , but decreases as 𝑚 𝑇 decreases. The Lévy-scale parameter 𝑅⁡(𝑚 𝑇 ) decreases with 𝑚 𝑇 and exhibits proportionality to the length scale of the nuclear overlap region. The Lévy exponent 𝛼⁡(𝑚 𝑇 ) is independent of 𝑚 𝑇 within uncertainties in each investigated centrality bin, but shows a clear centrality dependence. At all centralities, the Lévy exponent 𝛼 is significantly different from that of Gaussian (𝛼 = 2) or Cauchy (𝛼 = 1) source distributions. Comparisons to the predictions of Monte-Carlo simulations of resonance-decay chains show that, in all but the most peripheral centrality class (50%–60%), the obtained results are inconsistent with the measurements, unless a significant reduction of the in-medium mass of the 𝜂′ meson is included. Finally, in each centrality class, the best value of the in-medium 𝜂′ mass is compared to the mass of the 𝜂 meson, as well as to several theoretical predictions that consider restoration of U 𝐴⁢ (1) symmetry in hot hadronic matter.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Spatially resolved polarization swings in the supermassive binary black hole candidate OJ 287 with first Event Horizon Telescope observations

We present the first Event Horizon Telescope 1.3 mm observations of the supermassive binary black hole candidate OJ 287. The observations achieved an unprecedented angular resolution of 18 μas and reveal significant structural and polarization variability over just five days, marking the shortest timescale on which such changes have been directly imaged in this source. The inner jet exhibits a twisted ridgeline structure, with features displaying apparent superluminal motions up to about 22 c. The linear polarization maps reveal three main polarized features whose electric-vector position angles (EVPAs) change substantially over the time span of our observations, including a component with a radial polarization consistent with being produced by a recollimation shock. Most notably, we directly resolved two innermost jet components whose EVPAs rotate in opposite directions. The faster component, moving at 2.4 ± 0.9 μas/day (17.4 ± 6.5 c), exhibits counterclockwise EVPA swings of roughly 3.7° per day, while the slower component, with a proper motion of 1.4 ± 0.3 μas/day (10.2 ± 2.2 c), rotates clockwise at approximately 2.5° per day. Previous studies inferred helical magnetic fields in AGN jets from time-resolved or integrated polarization variability but lacked the angular resolution to directly image this effect. Our results provide spatially resolved evidence that a helical magnetic field threads the jet’s collimation and acceleration zone, ruling out models based on the superposition of unresolved components. Our analysis suggests that propagating shocks interact with a Kelvin–Helmholtz plasma instability, illuminating different phases of the helical magnetic field and producing the observed polarization spatial and temporal variability. Moreover, our model naturally accounts for the more rapid polarization rotation observed in the faster moving component. Our model predicts even more rapid swings in polarization, which could be tested with future observations featuring a more densely sampled time coverage.

OJ 287↗