IRON: Ion Routing Optimization Toolchain

• IRON Toolchain is a modular numerical workflow that generates time-dependent electrode voltages for efficient, low-excitation ion shuttling in segmented rf traps.
• It integrates electrostatic field solving, unconstrained quadratic optimization, digital waveform postprocessing, and dynamic simulation to ensure precision and adherence to hardware constraints.
• Performance benchmarks and example protocols demonstrate its capability to achieve sub-nanometer residuals and frequency stability, supporting scalable trapped-ion quantum processor designs.

The IRON (Ion Routing Optimization and eNvelope-shaping) toolchain is a numerical workflow for designing and validating high-fidelity, low-excitation transport protocols for trapped-ion qubits in segmented, radiofrequency (rf) trap architectures. The toolchain systematically generates time-dependent electrode voltages for complex traps, supporting arbitrary geometries including linear arrays, junctions, and multi-zone devices. By combining a boundary-element electrostatic solver, unconstrained quadratic optimization, digital waveform postprocessing, and full dynamical simulation, IRON enables efficient voltage protocol generation with rigorous control over hardware constraints and motional excitation, establishing a scalable numerical foundation for modern and next-generation trapped-ion quantum processors (Conta et al., 13 Jan 2026).

1. Modular Workflow of the IRON Toolchain

IRON consists of four key modules:

1. Electrostatic Field Solver: Utilizes a 3D CAD surface mesh of trap electrodes as input, applying Dirichlet boundary conditions (1 V on a single electrode, 0 V elsewhere) and solving Laplace’s equation,

$\Delta \phi_n(\mathbf{r}) = 0$

typically via a boundary-element method. The result is a unit-voltage potential $\phi_n(\mathbf{r})$ for each electrode, cached for fast field interrogation at arbitrary points.

1. Waveform Optimizer: Accepts the precomputed potentials, a prescribed shuttling trajectory $\{\mathbf{r}_{w,t}\}$, desired secular frequencies $\omega_{\mathrm{ref},u}$, orientations, and hardware constraints ($V_{\max}$, bandwidth). At each support point, $\phi_n(\mathbf{r})$ is expanded in a local multipole basis (regular solid harmonics up to order $L$) for efficient field and curvature calculation. The optimizer assembles a quadratic cost functional incorporating penalties for well position, confinement fidelity, voltage magnitude, slew rate, and fixed-voltage anchoring, then solves for voltages $\{V_{n,t}\}$ via linear algebra or unconstrained minimization.
2. Waveform Postprocessing: Processes the discrete voltages to match hardware (e.g., AWG) constraints. Steps include spline interpolation, time remapping to nonlinear profiles, grid resampling, and compensation for transfer-function filtering by inverting the finite impulse response (FIR) filter with Tikhonov-regularized least-squares.
3. Dynamical Simulator: Propagates the ion’s motion under the time-dependent, postprocessed voltages. Forces are computed at each timestep using

$$F(\mathbf{r},t) = Q \left[ \mathbf{E}_{\mathrm{rf}}(\mathbf{r}) + \sum_n V_{\mathrm{pre},n}(t) \, \mathbf{e}_n(\mathbf{r}) \right]$$

and motions are integrated (e.g., via velocity-Verlet), extracting residual motional excitation and secular frequency drifts.

This modular architecture enables rapid prototyping and reusability across different trap designs and shuttling protocols.

2. Mathematical Modeling of Trap Potentials

The electrostatic environment is governed by Laplace’s equation with Dirichlet voltage conditions imposed per electrode:

$$\Phi(\mathbf{r}, t) = \sum_{n=1}^N V_n(t) \, \phi_n(\mathbf{r})$$

where each $\phi_n(\mathbf{r})$0 is the unit response from electrode $\phi_n(\mathbf{r})$1. For computational efficiency, each $\phi_n(\mathbf{r})$2 is expanded about each path support point $\phi_n(\mathbf{r})$3 in a local solid-harmonic basis:

$\phi_n(\mathbf{r})$4

with $\phi_n(\mathbf{r})$5. Coefficients $\phi_n(\mathbf{r})$6 are extracted by sampling $\phi_n(\mathbf{r})$7 on a design sphere of radius $\phi_n(\mathbf{r})$8 using spherical-design or Fibonacci grids and Gram-matrix orthogonalization. Analytical differentiation of these expansions yields local fields and Hessians required for optimization objectives.

3. Optimization Problem and Cost Functional

Voltage protocols are encoded as $\phi_n(\mathbf{r})$9. The total quadratic cost functional is composed of five penalty terms:

• Well-position residual ($\{\mathbf{r}_{w,t}\}$0): Ensures proximity of the moving potential minimum to the target trajectory.
• Hessian/confinement residual ($\{\mathbf{r}_{w,t}\}$1): Matches the potential curvature to target secular frequencies and orientations.
• Voltage-magnitude penalty ($\{\mathbf{r}_{w,t}\}$2): Restricts voltage amplitudes.
• Voltage-slew/bandwidth penalty ($\{\mathbf{r}_{w,t}\}$3): Limits stepwise voltage variation.
• Fixed-voltage anchoring ($\{\mathbf{r}_{w,t}\}$4): Optionally constrains certain electrodes.

Explicitly,

$\{\mathbf{r}_{w,t}\}$5

$\{\mathbf{r}_{w,t}\}$6

Because the total cost is quadratic, the optimal voltages arise from solving $\{\mathbf{r}_{w,t}\}$7 with band-diagonal $\{\mathbf{r}_{w,t}\}$8 or by unconstrained minimization (e.g., BFGS).

4. Waveform Generation and Hardware Interface

Optimized discrete voltages are spline-interpolated to produce continuous profiles $\{\mathbf{r}_{w,t}\}$9 for $\omega_{\mathrm{ref},u}$0, possibly remapped nonlinearly (e.g., sinusoidal), grid-resampled to hardware time steps, and corrected for known AWG bandwidth limitations. The pre-ramp voltage $\omega_{\mathrm{ref},u}$1 is obtained by solving

$\omega_{\mathrm{ref},u}$2

where $\omega_{\mathrm{ref},u}$3 is the FIR filter, $\omega_{\mathrm{ref},u}$4 is the finite-difference operator, and $\omega_{\mathrm{ref},u}$5 is a regularization parameter. This inversion ensures that, after convolution through hardware and bandwidth limitations, the actual delivered waveform reproduces the intended voltage trajectory to within $\omega_{\mathrm{ref},u}$6 relative accuracy.

5. Dynamical Integration and Validation

Ion motion under real voltage protocols is simulated by integrating

$\omega_{\mathrm{ref},u}$7

using the velocity-Verlet method with fine time steps. Instantaneous motional modes are monitored by extracting local Hessians,

$\omega_{\mathrm{ref},u}$8

where $\omega_{\mathrm{ref},u}$9 are the principal curvatures.

Validation benchmarks against measured secular frequencies in standard geometries (e.g., Ruster-et-al. linear traps) show agreement within $V_{\max}$0 across design spheres $V_{\max}$1 from $V_{\max}$2 to $V_{\max}$3, well inside spectroscopic accuracy limits. Quartic multipole expansions with $V_{\max}$4 sphere samples agree to within $V_{\max}$5 of full ($V_{\max}$6) expansions. Mesh substructure effects on multipole coefficients remain $V_{\max}$7 for minimum face areas $V_{\max}$8.

6. Performance and Example Protocols

IRON demonstrates rapid prototyping performance. Representative timings (40 electrodes) are: 230 s for unit-potential solves (NullSpace), 11 s for multipole expansions (40 electrodes $V_{\max}$9 300 path points $\phi_n(\mathbf{r})$0 $\phi_n(\mathbf{r})$1), and 0.9 s for solution of a $\phi_n(\mathbf{r})$2 linear system (sparse CG).

Application examples confirm protocol fidelity:

• Linear five-segment shuttling: $\phi_n(\mathbf{r})$3 μm per segment, $\phi_n(\mathbf{r})$4 steps, with axial frequency $\phi_n(\mathbf{r})$5 MHz, radial $\phi_n(\mathbf{r})$6 MHz. Final residuals below 10 nm (axial), $\phi_n(\mathbf{r})$7 nm (radial), $\phi_n(\mathbf{r})$8, and principal-axis rotations $\phi_n(\mathbf{r})$9 mrad.
• X-junction corner shuttling: 90° bend, 50 μm radius, 36 dc electrodes, $L$0 steps, voltages $L$1 V. Axial $L$2 drops by $L$3, radial confinement exceeding axial, and residuals below 50 nm (axial) and 10 μm (vertical) in the worst region.

These protocols illustrate IRON’s capability for robust, low-excitation shuttling in arbitrary segmented traps.

7. Significance and Extensibility

IRON supplies an extensible and numerically efficient foundation for transport protocol design in segmented ion-trap quantum processors. Its general architecture accommodates flexible trap geometries, complex constraint sets, and hardware-aware waveform generation. By enabling rapid, reliable numerical prototyping, IRON supports the scalable deployment of fast, low-excitation shuttling necessary for large-scale quantum information processing (Conta et al., 13 Jan 2026). The explicit algorithmic steps and parameterizations provided enable full re-implementation for research and experimental validation of quantum transport protocols.