AthenaK: Portable GRMHD & NR Framework

• AthenaK is a performance-portable reimplementation of Athena++ that leverages the Kokkos programming model to run efficiently on diverse hardware including CPUs and GPUs.
• It employs advanced Godunov finite-volume methods, block-based adaptive mesh refinement, and staggered constrained transport to ensure robustness in astrophysical fluid dynamics and relativistic simulations.
• The framework demonstrates high scaling efficiency in exascale environments and supports a wide range of physics modules from GRMHD to Z4c numerical relativity for cutting-edge scientific applications.

to=arxiv_search.search ചികിത_json {"^^^^^^^^6query6 performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6^^^^^^^^", "^^^^^^^^6max_results6 ^^^^6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6query6^^^^} สูตรบาคาร่า to=arxiv_search.search 平台开号 to=arxiv_search.search code 天天中彩票中奖_json {"^^^^^^^^6query6 OR abs:\6^^^^^^^^"AthenaK\"", "^^^^^^^^6max_results6 ^^^^6max_results6query6^^^^} to=arxiv_search.search 天天中彩票不能买_json {"^^^^^^^^6query6 "^^^^^^^^6max_results6 ^^^^6max_results6query6^^^^} AthenaK is a performance-portable, open-source reimplementation of the Athena++ block-structured adaptive mesh refinement framework in C++ using the Kokkos programming model. It preserves the Athena-family emphasis on Godunov finite-volume methods, staggered constrained transport, and block-based AMR, while extending the framework to Newtonian, special relativistic, and general relativistic hydrodynamics and magnetohydrodynamics, GR radiation transport, particle modules, and—through companion developments—a Z^^^^6ti:\6^^^^c numerical relativity solver and GRMHD in dynamical spacetimes. The “K” denotes Kokkos-based performance portability across CPUs, NVIDIA and AMD GPUs, and ARM systems, and the code is explicitly oriented toward exascale workloads (&&&^^^^6query6^^^^&&&).

^^^^6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6^^^^. Lineage, scope, and design goals

AthenaK descends directly from Athena++, but it is not a minor port. The framework keeps the core Athena++ ideas—block-based AMR, Godunov finite-volume solvers, flux correction at refinement boundaries, and staggered constrained transport—while rewriting the implementation around a device-first execution model. In this reformulation, blocks are grouped into MeshBlockPack aggregates to reduce GPU kernel-launch overhead, dependent variables are stored in device-resident Kokkos views, and ghost-zone communication, prolongation, restriction, and boundary-condition application are performed inside device kernels before MPI exchange. The host is retained primarily for mesh-tree creation and refinement, task scheduling, and I/O (&&&^^^^6query6^^^^&&&).

The implemented physics coverage is correspondingly broad. The base framework includes Newtonian hydrodynamics and magnetohydrodynamics, special relativistic HD and MHD, fixed-metric GRHD and GRMHD, GR radiation hydrodynamics and magnetohydrodynamics, and particle modules for Lagrangian tracers and charged test particles. Companion developments add a Z^^^^6ti:\6^^^^c solver for the Einstein equations and a Valencia-formulation GRMHD solver in dynamical spacetimes, enabling binary black hole and binary neutron star applications within the same AMR/Kokkos infrastructure (&&&^^^^6query6^^^^&&&).

At the same time, AthenaK is deliberately narrower than some general-purpose mesh frameworks. In the ^^^^6max_results6query6max_results6ti:\6^^^^ framework paper, the code is restricted to Cartesian coordinates; characteristic reconstruction is not yet implemented; and the corner-transport upwind integrator is not included. Those constraints are architectural rather than accidental: the code prioritizes a compact, high-throughput, performance-portable implementation of rectangular-block astrophysical fluid dynamics and numerical relativity (&&&^^^^6query6^^^^&&&).

^^^^6max_results6^^^^. Framework architecture and numerical methods

AthenaK organizes the mesh as a binary tree, or an octree in three dimensions, with two-to-one refinement and a uniform timestep across refinement levels. Conservation across AMR interfaces is enforced by flux correction. Physics modules are independent objects with their own device-resident state and task lists, which avoids a combinatorial explosion of monolithic task graphs when multiple modules are combined in one run (&&&^^^^6query6^^^^&&&).

The numerical core is a reconstruct–solve–advance finite-volume scheme. Reconstruction options include PLM, classic PPM^^^^6ti:\6^^^^, extremum-preserving PPMX, and WENO-Z; Riemann solvers include LLF, HLLE, HLLC, HLLD, and Roe; and time integration uses method-of-lines SSP Runge–Kutta schemes from first to fourth order. In MHD, AthenaK uses the Athena++ upwind constrained-transport algorithm with face-centered magnetic fields and edge-centered electric fields, preserving PRESERVED_PLACEHOLDER_^^^^6query6^^^^ to roundoff. For stiff multi-fluid couplings, the framework includes IMEX second- and third-order schemes, including the positivity-preserving IMEX^^^^6max_results6^^^^ method of Krapp et al. as summarized in the framework description (&&&^^^^6query6^^^^&&&).

A representative conservative form used throughout the code family is

PRESERVED_PLACEHOLDER_^^^^6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6^^^^

with module-specific definitions of PRESERVED_PLACEHOLDER_^^^^6max_results6^^^^, PRESERVED_PLACEHOLDER_^^^^6query6^^^^, and PRESERVED_PLACEHOLDER_^^^^6ti:\6^^^^. In fixed-metric GRMHD, AthenaK follows the conservative HARM-like formulation summarized in the framework paper; in dynamical-spacetime GRMHD it uses the Valencia conservative Eulerian form, evolving densitized variables PRESERVED_PLACEHOLDER_^^^^6 OR abs:\6^^^^ with metric source terms handled consistently in the ^^^^6query6^^^^+^^^^6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6^^^^ split (&&&^^^^6query6^^^^&&&).

Robustness in extreme flows is strongly tied to first-order flux correction. If an RK substep would generate negative density or pressure, or fail primitive recovery, AthenaK recomputes local interface fluxes with first-order donor-cell reconstruction and a more diffusive LLF or HLLE flux. In the dynamical-spacetime GRMHD work this FOFC mechanism, optionally combined with a relaxed discrete maximum principle, is central to reducing reliance on artificial atmospheres and stabilizing neutron-star merger calculations (Fields et al., 2024).

Radiation transport adds another layer of specialization. The finite-solid-angle GR radiation module, originally developed in Athena++ and ported to AthenaK, discretizes the specific intensity in angle, advances it explicitly in space and angle, and couples matter and radiation through a locally implicit source update. The method is designed for stationary spacetimes and highly anisotropic radiation fields, where moment closures such as M^^^^6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6^^^^ are inadequate (White et al., 2023).

^^^^6query6^^^^. Relativistic extensions and numerical relativity

AthenaK’s numerical relativity module evolves the vacuum Einstein equations in the Z^^^^6ti:\6^^^^c formulation. In the implementation described by Zhu et al., Z^^^^6ti:\6^^^^c is combined with the moving-puncture gauge, oct-tree AMR, and cell-centered finite-difference operators designed for GPUs. The module demonstrates sixth-order spatial convergence on linear waves when four ghost cells are used, and fourth-order self-convergence for the gravitational waveform in an equal-mass, non-spinning binary black hole calibration problem. New refinement criteria based on $\chi_{\min}$ and $|\nabla\chi|$ focus resolution near punctures and around smaller holes in high-mass-ratio binaries, while wave-zone derefinement controls protect extraction accuracy at large radius (Zhu et al., 2024).

The GRMHD extension to dynamical spacetimes uses the Valencia conservative Eulerian formulation,

$$\partial_t(\sqrt{\gamma}\,\mathbf{U}) + \partial_i(\sqrt{\gamma}\,\mathbf{F}^i) = \sqrt{\gamma}\,\mathbf{S},$$

with conservative variables $\mathbf{U} = [D, S_i, \tau, B^i]$. In this formulation, spacetime variables are cell-centered rather than vertex-centered, which avoids interpolations between metric and fluid layouts and improves GRMHD performance on block-based AMR meshes. Validation spans magnetized shock tubes, cylindrical blast waves, magnetic-loop advection, oscillating and free-evolving TOV stars, and SANE accretion around a Kerr black hole (Fields et al., 2024).

These developments culminated quickly in production applications. The GW^^^^6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6 OR abs:\6query6AthenaK6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6ti:\6^^^^ target simulation performed with AthenaK uses the moving-puncture approach, Z^^^^6ti:\6^^^^c evolution, sixth-order finite differencing in space, and a low-storage fourth-order Runge–Kutta scheme in time. Initial data are constructed with TwoPunctures and matched to the LazEV-targeted parameters of Lovelace et al. (^^^^6max_results6query6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv66^^^^); horizon quantities are computed with Einstein Toolkit tools; and waveforms are extracted both from finite-radius PRESERVED_PLACEHOLDER_^^^^6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6query6^^^^ and via Cauchy Characteristic Extraction to future null infinity PRESERVED_PLACEHOLDER_^^^^6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6^^^^. For the final remnant, the code reports PRESERVED_PLACEHOLDER_^^^^6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6max_results6^^^^, PRESERVED_PLACEHOLDER_^^^^6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6query6^^^^, and PRESERVED_PLACEHOLDER_^^^^6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6ti:\6^^^^, with agreement against the highest-resolution SpEC result at the level of PRESERVED_PLACEHOLDER_^^^^6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6 OR abs:\6^^^^ in mass, PRESERVED_PLACEHOLDER_^^^^6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv66^^^^ in spin, and about PRESERVED_PLACEHOLDER_^^^^6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6^^^^7 in recoil (&&&^^^^6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6^^^^&&&).

AthenaK’s relativistic reach is not confined to vacuum or compact-binary waveforms. The code has also been used for fixed-spacetime GRMHD simulations of thin accretion disks in Cartesian Kerr–Schild coordinates, where it validated a gravity-dominated geodesic plunge model and found a non-zero PRESERVED_PLACEHOLDER_^^^^6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6^^^^8 angular-momentum drop across the plunging region, together with an order-of-magnitude rise in the local PRESERVED_PLACEHOLDER_^^^^6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6^^^^9 parameter inside the ISCO (&&&^^^^6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6max_results6^^^^&&&).

^^^^6ti:\6^^^^. Performance portability and scaling

AthenaK’s central technical claim is that high-order astrophysical fluid dynamics, GRMHD, and numerical relativity can be implemented once and run efficiently across heterogeneous architectures. The framework paper reports over one billion cell updates per second for three-dimensional hydrodynamics on a single NVIDIA Grace Hopper processor and a typical parallel efficiency of PRESERVED_PLACEHOLDER_^^^^6max_results6query6^^^^ on 6^^^^6 OR abs:\6^^^^,^^^^6 OR abs:\6query66^^^^ AMD GPUs on Frontier (&&&^^^^6query6^^^^&&&).

The reported performance metrics span both framework-level and application-level measurements.

Context	Reported result	Source
Single Grace Hopper processor	over one billion cell updates per second for ^^^^6query6^^^^D hydrodynamics	(&&&^^^^6query6^^^^&&&)
Framework weak scaling on Frontier	typical parallel efficiency of PRESERVED_PLACEHOLDER_^^^^6max_results6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6^^^^ on 6^^^^6 OR abs:\6^^^^,^^^^6 OR abs:\6query66^^^^ AMD GPUs	(&&&^^^^6query6^^^^&&&)
NR module weak scaling on Frontier	PRESERVED_PLACEHOLDER_^^^^6max_results6max_results6^^^^ efficiency on up to 6^^^^6 OR abs:\6^^^^,^^^^6 OR abs:\6query66^^^^ AMD MI^^^^6max_results6 OR abs:\6query6^^^^X GPUs, relative to ^^^^6ti:\6^^^^ GPUs	(Zhu et al., 2024)
NR module strong scaling	PRESERVED_PLACEHOLDER_^^^^6max_results6query6^^^^ on AMD MI^^^^6max_results6 OR abs:\6query6^^^^X and PRESERVED_PLACEHOLDER_^^^^6max_results6ti:\6^^^^ on NVIDIA A^^^^6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6query6query6^^^^ GPUs	(Zhu et al., 2024)
Dynamical-spacetime GRMHD weak scaling	PRESERVED_PLACEHOLDER_^^^^6max_results6 OR abs:\6^^^^ to ^^^^6query6max_results6^^^^,768 GPUs and PRESERVED_PLACEHOLDER_^^^^6max_results66^^^^ to 6^^^^6 OR abs:\6^^^^,^^^^6 OR abs:\6query66^^^^ GPUs	(Fields et al., 2024)
SuperMUC-NG single-node AthenaK study	PRESERVED_PLACEHOLDER_^^^^6max_results6^^^^7 throughput gain and PRESERVED_PLACEHOLDER_^^^^6max_results6^^^^8 energy-efficiency gain for GPU offload at mesh-block edge length PRESERVED_PLACEHOLDER_^^^^6max_results6^^^^9	(&&&^^^^6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6^^^^9&&&)

Several details qualify these headline numbers. First, mesh-block granularity matters. In the SuperMUC-NG node-level study, AthenaK’s peak GPU throughput for the GW^^^^6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6 OR abs:\6query6AthenaK6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6ti:\6^^^^ BBH workload is PRESERVED_PLACEHOLDER_^^^^6query6query6^^^^ zone-cycles/s and the peak energy efficiency is PRESERVED_PLACEHOLDER_^^^^6query6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6^^^^ zone-cycles/J, but those gains saturate only once the mesh-block edge length reaches about ^^^^6max_results6ti:\6^^^^. At edge lengths PRESERVED_PLACEHOLDER_^^^^6query6max_results6^^^^, GPU occupancy and launch overhead erode both throughput and per-Joule efficiency (&&&^^^^6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6^^^^9&&&).

Second, end-to-end production timings are now concrete rather than aspirational. The GW^^^^6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6 OR abs:\6query6AthenaK6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6ti:\6^^^^ production run on Aurora used ^^^^6query6max_results6^^^^ nodes, or ^^^^6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6AthenaK6max_results6^^^^ GPUs, and required about ^^^^6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6query6query6^^^^ wall-clock hours for a ^^^^6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6query6^^^^-orbit BBH evolution with CCE worldtube output and post-processing. In the paper’s interpretation, modern GPU-accelerated systems are roughly PRESERVED_PLACEHOLDER_^^^^6query6query6^^^^ faster than leadership-scale machines available in ^^^^6max_results6query6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv66^^^^ for this class of workflow (&&&^^^^6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6^^^^&&&).

^^^^6 OR abs:\6^^^^. Scientific applications

AthenaK has been used in a notably diverse set of workflows. In binary black hole merger modeling, the GW^^^^6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6 OR abs:\6query6AthenaK6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6ti:\6^^^^ calculation established a full “Cauchy-to-infinity” pipeline: TwoPunctures initial data, Z^^^^6ti:\6^^^^c/moving-puncture evolution with AMR, isolated-horizon remnant analysis, finite-radius PRESERVED_PLACEHOLDER_^^^^6query6ti:\6^^^^ extraction, and CCE waveforms at PRESERVED_PLACEHOLDER_^^^^6query6 OR abs:\6^^^^. The dominant PRESERVED_PLACEHOLDER_^^^^6query66^^^^ mode agrees with SXS:BBH:^^^^6query6query6query6 OR abs:\6^^^^ and RIT:BBH:^^^^6query6query66max_results6^^^^ at the level of maximum dephasing PRESERVED_PLACEHOLDER_^^^^6query6^^^^7 and PRESERVED_PLACEHOLDER_^^^^6query6^^^^8 at merger, respectively, with nearly flat amplitude differences of PRESERVED_PLACEHOLDER_^^^^6query6^^^^9 through the inspiral. The resulting AthenaK waveform was then used in a bilby re-analysis of GW^^^^6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6 OR abs:\6query6AthenaK6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6ti:\6^^^^, yielding PRESERVED_PLACEHOLDER_^^^^6ti:\6query6^^^^, PRESERVED_PLACEHOLDER_^^^^6ti:\6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6^^^^, and PRESERVED_PLACEHOLDER_^^^^6ti:\6max_results6^^^^, broadly consistent with semi-analytic LVK models (&&&^^^^6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6^^^^&&&).

In compact-object GRMHD, the code now spans both black-hole accretion and neutron-star mergers. The plunging-region study around a Schwarzschild black hole used AthenaK for a ^^^^6max_results6query6^^^^,^^^^6query6query6query6^^^^ PRESERVED_PLACEHOLDER_^^^^6ti:\6query6^^^^ thin-disk evolution and found that a gravity-dominated geodesic plunge remains accurate provided non-adiabatic heating is included, plausibly from grid-scale magnetic dissipation in a mid-plane current sheet. The conclusion that constant PRESERVED_PLACEHOLDER_^^^^6ti:\6ti:\6^^^^-disk models are physically inappropriate inside the ISCO follows directly from the measured rise in entropy and magnetic stress across the plunge (&&&^^^^6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6max_results6^^^^&&&). In magnetized binary neutron star mergers, AthenaK has been used with the SFHo tabulated EOS, trapped-lepton advection, WENOZ reconstruction, HLLE, and upwind constrained transport at resolutions down to PRESERVED_PLACEHOLDER_^^^^6ti:\6 OR abs:\6^^^^. Between the two highest resolutions the orbital dephasing over more than seven orbits is only PRESERVED_PLACEHOLDER_^^^^6ti:\66^^^^ radians, magnetic fields are amplified and a polar funnel forms, but baryonic pollution prevents a magnetically dominated outflow by about ^^^^6query6query6^^^^ ms after merger (&&&^^^^6max_results6ti:\6^^^^&&&).

A related neutron-star remnant application used AthenaK to study magnetic flux emergence from a twisted toroidal flux tube in a hot massive post-merger star. In that controlled fixed-metric GRMHD experiment, magnetic buoyant emergence occurs only for extremely large fields significantly exceeding PRESERVED_PLACEHOLDER_^^^^6ti:\6^^^^7, while more typical fields around PRESERVED_PLACEHOLDER_^^^^6ti:\6^^^^8 are dominated by hydrodynamic effects. The reported outcome places strong limitations on remnant-driven magnetically powered outflows in binary neutron star remnants (&&&^^^^6max_results6 OR abs:\6^^^^&&&).

The code is also used well outside numerical relativity. AthenaK has run fully three-dimensional two-fluid MHD turbulence calculations in partially ionized media at PRESERVED_PLACEHOLDER_^^^^6ti:\6^^^^9, showing how velocity, density, and magnetic-field statistics change across neutral–ion coupling regimes and how decoupling occurs over a range of scales rather than at a single ambipolar-diffusion scale (&&&^^^^6max_results66^^^^&&&). It has supported a ^^^^6ti:\6query6^^^^96PRESERVED_PLACEHOLDER_^^^^6 OR abs:\6query6^^^^ multiphase interstellar turbulence simulation on Frontier, with in situ synthetic dust-polarization map generation, reproducing Planck ^^^^6query6 OR abs:\6query6^^^^ GHz inertial-range slopes PRESERVED_PLACEHOLDER_^^^^6 OR abs:\6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6^^^^, PRESERVED_PLACEHOLDER_^^^^6 OR abs:\6max_results6^^^^, and PRESERVED_PLACEHOLDER_^^^^6 OR abs:\6query6^^^^, together with PRESERVED_PLACEHOLDER_^^^^6 OR abs:\6ti:\6^^^^ and PRESERVED_PLACEHOLDER_^^^^6 OR abs:\6 OR abs:\6^^^^ (&&&^^^^6max_results6^^^^7&&&). It has also been used for controlled isothermal turbulence experiments comparing Fourier-space driving and point-source driving, where the canonical turbulence parameter PRESERVED_PLACEHOLDER_^^^^6 OR abs:\66^^^^ was shown to be degenerate with both the spatial locality and temporal correlation of the forcing (&&&^^^^6max_results6^^^^8&&&).

6. Reproducibility, limitations, and current open problems

AthenaK is public infrastructure. The framework is open source under BSD-^^^^6query6^^^^ at https://github.com/IAS-Astrophysics/athenak, uses a CMake-based build, selects physics and algorithms at runtime through input files, writes MPI-IO output in a compact binary format, and provides Python translators to HDF^^^^6 OR abs:\6^^^^. For the GW^^^^6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6 OR abs:\6query6AthenaK6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6ti:\6^^^^ application, the authors additionally released a step-by-step tutorial with all necessary input files, analysis scripts, and SpECTRE CCE converters at https://github.com/dradice/athenak-tutorial-gw^^^^6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6 OR abs:\6query6AthenaK6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6ti:\6^^^^ (&&&^^^^6query6^^^^&&&).

The code’s current limitations are explicit in the literature. At the framework level, AthenaK remains restricted to Cartesian coordinates, does not yet implement characteristic reconstruction or CTU, and has not yet ported some advanced diffusion machinery such as operator-split implicit solvers and super-time-stepping. In the ^^^^6max_results6query6max_results6ti:\6^^^^ framework paper, tabulated microphysical EOS support is described as under test rather than mature base functionality (&&&^^^^6query6^^^^&&&). In the GW^^^^6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6 OR abs:\6query6AthenaK6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6ti:\6^^^^ production study, only a single AthenaK resolution is used, so the error budget is established by cross-validation against SXS and LazEV rather than a formal convergence sequence; gauge differences and minor AMR backscattering noise around PRESERVED_PLACEHOLDER_^^^^6 OR abs:\6^^^^7 in some modes are noted explicitly, and more detailed gauge-invariant waveform comparisons are deferred (&&&^^^^6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6^^^^&&&).

More fundamentally, some AthenaK applications have exposed limitations that are physical and numerical rather than purely infrastructural. In turbulent radiative mixing layers, the apparent resolution independence of total cooling is shown to arise from a cancellation between numerical dissipation and numerical viscosity; the paper argues that this is “remarkable, and perhaps fortuitous,” and recommends resolving the turbulent Field length PRESERVED_PLACEHOLDER_^^^^6 OR abs:\6^^^^8, defined by PRESERVED_PLACEHOLDER_^^^^6 OR abs:\6^^^^9, rather than treating converged $\chi_{\min}$^^^^6query6^^^^ as sufficient evidence of a well-resolved layer (&&&^^^^6query6max_results6^^^^&&&). In magnetized BNS mergers, the failure to achieve $\chi_{\min}$^^^^6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6^^^^ in the polar funnel is traced partly to baryon pollution linked to artificial heating of stellar surfaces, and the authors identify improved Riemann solvers, neutrino cooling and heating, and longer-term evolution as plausible next steps; AthenaK developments under test include HLLD and neutrino transport modules labeled M^^^^6AthenaK performance portable Athena++ AthenaK numerical relativity GRMHD arXiv6^^^^ and full Boltzmann via FP_N (&&&^^^^6max_results6ti:\6^^^^&&&).

Taken together, these papers define AthenaK less as a single application code than as a unified numerical platform. Its distinguishing feature is not merely GPU acceleration, but the combination of block-structured AMR, constrained transport, relativistic fluid solvers, Z^^^^6ti:\6^^^^c numerical relativity, and Kokkos-based performance portability in one code base. The published record already shows that this combination is sufficient for many-orbit BBH waveform production, GPU BNS merger calculations, GR radiation transport, exascale interstellar turbulence, and controlled studies of numerical resolution effects across astrophysical fluid dynamics (&&&^^^^6query6^^^^&&&).