Keywords: orbital attitude control, electromagnetic interference, machine learning scheduling, Kalman filter, magnetorquers, high-power compute platforms, spacecraft systems co-design, distributed training.
1. INTRODUCTION
The deployment of large-scale machine learning compute infrastructure in Earth orbit has become technically and economically plausible within the past several years, driven by reusable launch vehicle economics, the modular architecture of modern GPU/NPU clusters, and the emergence of commercial orbital platform services. Proposed orbital data center concepts [1,2] anticipate continuous power draws of 10-100 MW and higher, enabled by large solar array deployments or space-based solar power architectures [3].
The systems engineering of orbital compute platforms at this power level introduces a class of interference problem that has no precedent in spacecraft design history. All previous spacecraft — including the International Space Station, which draws approximately 84 kW at peak [4] — operate at power levels where the electromagnetic effects of internal power distribution are negligible relative to attitude control system authority. This negligibility underpins the standard practice of designing power distribution and attitude control as independent subsystems with no required coordination between them.
This independence assumption fails at megawatt-scale compute platform power levels. The failure mechanism is specific to machine learning training workloads, which are characterized by sharp transient current demands — training burst events — rather than the steady or slowly-varying power draws of conventional spacecraft subsystems. A training burst event in a large distributed NPU cluster draws tens of kiloamperes over seconds to tens of seconds. This rapid current change generates a time-varying magnetic dipole moment proportional to the current-area product of the bus geometry. At the power levels and bus geometries relevant to orbital AI infrastructure, this spurious dipole moment is comparable to or greater than the attitude control authority of the magnetorquer system responsible for reaction wheel desaturation.
The consequence is direct: a training burst event can torque the orbital platform, misalign thermal radiator panels, desaturate reaction wheels, and in extreme cases induce structural loading inconsistent with orbital platform design margins. None of these consequences are captured by any existing spacecraft EMI standard [5,6] or attitude control specification, because no existing standard contemplates a spacecraft subsystem capable of generating this magnitude of internal magnetic disturbance.
This paper makes three contributions. First, we formalize the coupling mechanism between ML compute scheduling and orbital attitude control, deriving the interference threshold equation and quantifying the failure margin at representative orbital compute platform parameters. Second, we propose HERALD — a co-designed scheduling and attitude control architecture that enforces the derived interference constraint as a hard scheduling invariant while preserving near-optimal training throughput. Third, we extend the HERALD framework to handle the electromagnetic contributions of power beaming rectenna systems, which represent an additional uncounted interference source in solar-power-beaming orbital platform architectures, and to coordinate plasma phased-array fleet shielding across multi-node deployments.
The paper is organized as follows. Section 2 reviews the relevant background in orbital attitude control, spacecraft EMI standards, and ML training schedulers. Section 3 derives the interference coupling equations and quantifies the failure regime. Section 4 presents the HERALD architecture and state-space formulation. Section 5 addresses the rectenna harmonic interference extension. Section 6 presents the multi-node plasma phased-array coordination protocol. Section 7 discusses limitations and implementation requirements. Section 8 concludes.
2. BACKGROUND
2.1 Orbital Attitude Control with Magnetorquers
Orbital spacecraft attitude control typically uses a combination of reaction wheels for fine attitude control and magnetorquers (current-carrying coils or rods that interact with Earth's geomagnetic field) for reaction wheel desaturation [7,8]. The magnetorquer produces a torque by interacting with the local geomagnetic field B_env:
τ_control = M_control × B_env (1)
where M_control is the magnetic dipole moment commanded by the attitude control system [A·m²] and × denotes the vector cross product. The maximum attitude control torque available from the magnetorquer system is bounded by the maximum achievable dipole moment M_auth — the magnetorquer authority:
||M_control|| ≤ M_auth (2)
For a representative LEO orbital platform at 500 km altitude with Earth's magnetic field strength B_env ≈ 40 μT, an MTQ800-class magnetorquer array achieves M_auth ≈ 200,000 A·m² [9]. The attitude control response time τ_control — the time over which the magnetorquer can effect a meaningful attitude correction — is typically 1-10 seconds for reaction wheel desaturation in LEO [10].
2.2 Spacecraft Internal EMI Standards
MIL-STD-461 [5] and ECSS-E-ST-20-07 [6] specify conducted and radiated emission limits for spacecraft electrical systems. These standards were developed for spacecraft with power levels of tens to hundreds of kilowatts and focus on interference with sensitive scientific instruments and communication systems. They do not address the generation of attitude-relevant magnetic disturbance torques by internal power distribution transients, because no prior spacecraft has operated at power levels where such disturbances are significant relative to attitude control authority.
The applicable MIL-STD-461 limits for conducted emissions on power leads specify maximum current noise spectral density in the frequency range 30 Hz to 10 kHz. Training burst events produce current transients with characteristic frequencies of 0.1-1 Hz — below the lower limit of the MIL-STD-461 conducted emission specification. This gap in the standards reflects the absence of any prior spacecraft with comparable internal power transients at these frequencies.
2.3 Machine Learning Training Schedulers
Distributed ML training on GPU/NPU clusters is managed by schedulers that allocate compute resources to training jobs, manage data pipeline throughput, and coordinate gradient synchronization across nodes [11,12]. The power draw of a training cluster is determined by the compute utilization profile of the scheduled jobs — periods of high utilization (training burst events) are interspersed with periods of lower utilization (data loading, gradient synchronization, checkpointing).
The current draw of a modern GPU during a training burst scales approximately as:
I_burst(t) = P_TDP / V_bus · f_util(t) (3)
where P_TDP is the thermal design power of the GPU, V_bus is the bus voltage, and f_util(t) ∈ [0,1] is the utilization fraction at time t. For a cluster of N_GPU GPUs at bus voltage V_bus, the total cluster current during a burst event is:
I_cluster(t) = N_GPU · P_TDP · f_util(t) / V_bus (4)
The rate of current change during burst initiation — the critical parameter for attitude control interference — depends on the scheduler's burst initiation protocol. Standard schedulers initiate training bursts as fast as the hardware allows, typically achieving full utilization within 100-500 ms. This produces dI/dt values in the range 10^3 - 10^6 A/s for large clusters, depending on bus voltage and cluster size.
2.4 Kalman Filter Attitude Estimation
The extended Kalman filter (EKF) is the standard state estimator for spacecraft attitude control [13]. The EKF maintains a state estimate x_t and error covariance P_t, updated by the prediction-correction cycle:
x_{t|t-1} = f(x_{t-1}, u_t) [prediction] (5)
x_t = x_{t|t-1} + K_t(z_t − h(x_{t|t-1})) [correction] (6)
where f is the state transition function, u_t is the control input, z_t is the measurement vector, h is the measurement function, and K_t is the Kalman gain. For standard spacecraft attitude control, the state vector includes attitude quaternion q, angular velocity ω, and gyroscope bias b_g. HERALD extends this standard formulation to include bus current state and its derivative, creating a coupled estimator that jointly tracks attitude dynamics and compute load dynamics. This extension is the central technical contribution of the HERALD architecture.
3. THE COUPLING MECHANISM: DERIVATION AND QUANTIFICATION
3.1 Spurious Dipole Moment from Bus Current Transients
A current-carrying conductor loop of area A carrying current I produces a magnetic dipole moment:
M = I · A · n̂ (7)
where n̂ is the unit normal to the loop plane. For a spacecraft DC bus, the effective loop area A_eff is determined by the physical routing of the bus conductors and the geometry of the return current path. For a centralized bus architecture with conductors routed along a spacecraft truss of characteristic dimension L, A_eff ≈ L² for a roughly rectangular current loop. For a distributed per-rack bus architecture, A_eff is reduced by the constraint that each rack's current loop is small relative to the total bus geometry.
The spurious dipole moment during a training burst event is:
M_spurious(t) = I_cluster(t) · A_eff (8)
The rate of change of spurious dipole moment during burst initiation is:
dM_spurious/dt = A_eff · dI_cluster/dt (9)
3.2 The Interference Threshold
Attitude control interference becomes significant when the spurious dipole moment competes with the attitude control system's commanded moment. We define the interference threshold as the condition under which the spurious moment exceeds a fraction ε_int of the magnetorquer authority:
M_spurious ≥ ε_int · M_auth (10)
For ε_int = 0.1 (10% interference threshold — the level at which attitude perturbations become measurable in attitude sensor data), the maximum allowable cluster current is:
I_max = ε_int · M_auth / A_eff (11)
The maximum allowable rate of current change during burst initiation is:
dI/dt|_max = ε_int · M_auth / (A_eff · τ_control) (12)
Equation (12) is the fundamental scheduling constraint. Any burst initiation sequence that produces dI/dt > dI/dt|_max during the attitude control response window τ_control will generate attitude perturbations inconsistent with platform pointing requirements.
3.3 Numerical Evaluation at Representative Platform Parameters
Table 1 evaluates equations (11) and (12) at parameters representative of a 40 MW orbital AI compute platform, derived from published data for ISS-heritage bus architecture [4], commercially available magnetorquer systems [9], and Starcloud/Lumen Orbit modular cluster concepts [1,2].
| Parameter | Symbol | Value | Source/Basis |
|---|---|---|---|
| Platform compute power | P_compute | 40 MW | Starcloud-class modular cluster |
| DC bus voltage | V_bus | 400 V | ISS heritage; scalable to 4 kV HVDC |
| Cluster current (training burst) | I_cluster | 10,000–100,000 A | At 400 V; range reflects utilization variation |
| Effective loop area (centralized bus) | A_eff | 20 m² | Compact truss routing, conservative estimate |
| Effective loop area (distributed) | A_eff | < 2 m² per rack | Per-rack feeders reduce loop area |
| Magnetorquer authority | M_auth | 200,000 A·m² | 10–20 MTQ800-class rods with ferromagnetic cores |
| Attitude control response time | τ_control | 5 s | LEO B-field ~40 μT; conservative desaturation |
| Interference threshold (10%) | ε_int · M_auth | 20,000 A·m² | From equation (10) |
| Maximum allowable current | I_max | 1,000–10,000 A | From equation (11); range = centralized/distributed |
| Maximum allowable dI/dt | dI/dt|_max | 200–2,000 A/s | From equation (12); range = centralized/distributed |
Table 1. HERALD constraint parameters at 40 MW orbital compute platform scale.
The interference ratio M_spurious/M_auth reaches 10–500 at centralized bus architecture — a factor of 100–5,000 above the interference threshold. This is a design-critical coupling, not a second-order effect. The interference ratio at aggressive parameter values — 500 MW·m² of spurious dipole moment against 200,000 A·m² of control authority — represents a factor of 2,500 above threshold. Standard attitude control algorithms operating without knowledge of the compute load would experience sustained uncompensated disturbance torques, producing attitude errors potentially exceeding pointing requirements by orders of magnitude for platform-wide training runs.
KEY FINDING: At 40 MW scale with centralized bus architecture, training burst events produce spurious magnetic dipole moments up to 2,500 times the magnetorquer interference threshold. This is not a perturbation to be corrected — it is a dominant torque input that the attitude control system has no visibility into under the current decoupled design paradigm.
4. THE HERALD ARCHITECTURE
4.1 Design Principles
HERALD (Harmonic EM-Resolved Attitude-Load Dispatcher) addresses the coupling problem through joint co-design of the ML training scheduler and the attitude control estimator. Three design principles guide the architecture:
Principle 1 — Predict, don't react. The attitude control system should have advance knowledge of planned training bursts, not discover their electromagnetic effects after the fact. This requires feeding the training job queue forward into the attitude estimator.
Principle 2 — Enforce constraints at the scheduler, not the actuator. The attitude control system should not be required to compensate for burst-induced disturbances — it has limited bandwidth and authority. Instead, the scheduler should be prevented from initiating bursts that would require compensation. The dI/dt constraint is a scheduling constraint, not an attitude control compensation problem.
Principle 3 — Co-design the state vector. The Kalman filter estimator should maintain joint state over attitude dynamics and bus current dynamics. A combined state vector enables optimal estimation of both subsystems with explicit representation of their coupling.
4.2 Extended State Vector
The HERALD state vector extends the standard attitude Kalman filter to include bus current state:
x = [q, ω, b_g, I_bus, dI_bus/dt, M_residual, I_rect, H_rect]ᵀ (13)
where:
- q ∈ SO(3) — attitude quaternion [4 components]
- ω ∈ ℜ³ — angular velocity [rad/s]
- b_g ∈ ℜ³ — gyroscope bias [rad/s]
- I_bus ∈ ℜ — instantaneous DC bus current [A]
- dI_bus/dt ∈ ℜ — bus current rate of change [A/s]
- M_residual ∈ ℜ³ — residual magnetic dipole after magnetorquer compensation [A·m²]
- I_rect ∈ ℜ — power beaming rectenna switching current [A]
- H_rect ∈ ℜᴷ — rectenna switching harmonic content vector [K harmonics]
The inclusion of I_rect and H_rect in the state vector is required because power beaming rectenna systems — a primary power source for orbital compute platforms — produce switching transients at 5-20 kHz whose harmonics extend into the magnetorquer bandwidth. These harmonics are an independent, broadband interference source not captured by the training burst model alone (Section 5).
4.3 State Transition Model
The state transition model couples attitude dynamics and bus current dynamics through the spurious dipole moment term:
q_{t+1} = q_t ⊗ Δq(ω_t, τ_total, Δt) (14)
ω_{t+1} = ω_t + J⁻¹(τ_total − ω_t × Jω_t) · Δt (15)
where J is the platform moment of inertia tensor, ⊗ denotes quaternion multiplication, and τ_total is the total torque:
τ_total = τ_control + τ_spurious + τ_disturbance (16)
τ_spurious = M_spurious × B_env = (I_bus · A_eff + I_rect · A_rect) × B_env (17)
The spurious torque τ_spurious is now an explicit term in the attitude dynamics model, computed from the bus current state and the known bus geometry parameters A_eff and A_rect. This makes the spurious torque a predicted disturbance (compensated by the Kalman filter) rather than an unmodeled noise term.
The bus current dynamics are modeled as:
I_{bus,t+1} = I_{bus,t} + (dI_bus/dt)_t · Δt + w_I (18)
(dI_bus/dt)_{t+1} = f_scheduler(job_queue_t, I_{bus,t}) + w_{dI} (19)
where w_I and w_{dI} are process noise terms and f_scheduler is the HERALD dispatch function (Section 4.5) that predicts the current rate of change from the pending job queue.
4.4 Measurement Model
The HERALD measurement vector includes standard attitude sensors augmented by current measurement:
z_t = [q_star, ω_gyro, B_measured, I_bus_measured, I_rect_measured]ᵀ (20)
where q_star is the star tracker attitude measurement, ω_gyro is the gyroscope angular velocity measurement, B_measured is the magnetometer measurement of the local magnetic field (including contributions from all current loops), and I_bus_measured, I_rect_measured are direct current measurements from bus current sensors.
The magnetometer measurement model includes contributions from both the geomagnetic field and the platform's internal current loops:
B_measured = B_env + B_spurious + v_B (21)
B_spurious = μ₀/(4π) · [3(M_total · r̂)r̂ − M_total] / r³ (22)
where M_total = M_control + M_spurious + M_residual is the total magnetic moment of the platform, r is the distance from the dipole to the magnetometer, and v_B is measurement noise. The inclusion of B_spurious in the measurement model allows the Kalman filter to use magnetometer measurements to refine estimates of M_residual — the residual dipole after commanded magnetorquer compensation.
4.5 The HERALD Dispatch Algorithm
The HERALD dispatch algorithm is a constrained scheduler that enforces the dI/dt constraint derived in equation (12) while optimizing training throughput. The algorithm operates as follows:
HERALD Dispatch Algorithm:
Input: job_queue (ordered list of pending training jobs with resource requirements)
x_t (current Kalman state estimate including I_bus, dI_bus/dt)
dI_max = ε_int · M_auth / (A_eff · τ_control) [constraint, from Eq. 12]
For each candidate job j in job_queue:
- Predict current trajectory if j is initiated at t: I_predicted(t') = I_bus,t + ΔI_j(t'-t) for t' ∈ [t, t+T_ramp_j] where ΔI_j is the current ramp profile for job j
- Compute predicted dI/dt over ramp window: dI_predicted/dt = max |dI_predicted(t')/dt| for t' ∈ [t, t+T_ramp_j]
- Check constraint: if dI_predicted/dt > dI_max: defer j; compute earliest feasible initiation time t_j* t_j* = t + (dI_predicted/dt - dI_max) · T_ramp_j / ΔI_j_total else: initiate j; update I_bus forecast
- Update Kalman state with initiated job's current profile as known input
The dispatch algorithm produces a smooth current envelope that satisfies the dI/dt constraint at every point. Jobs are not cancelled — they are deferred to the earliest feasible initiation time. The throughput cost of this deferral depends on the burst frequency and the tightness of the constraint; Section 7.1 analyzes this cost quantitatively.
The joint optimization objective, incorporating training staleness [14] as an additional scheduling signal:
J = min Σⱼ w₁ · staleness(j) + w₂ · delay(j, t_j*) (23)
subject to dI/dt ≤ dI/dt|_max for all t.
Here staleness(j) is the gradient staleness of job j (inversely proportional to training urgency) and delay(j, t_j*) is the deferral time imposed by the constraint. Jobs with low staleness — whose gradients are current — are prioritized for available current budget. Jobs with high staleness tolerate deferral better, allowing the scheduler to smooth current demand while preserving training quality.
5. RECTENNA HARMONIC INTERFERENCE EXTENSION
5.1 Power Beaming as an Uncounted EM Source
Orbital compute platforms at megawatt scale require power sources beyond what solar arrays alone can provide at reasonable panel area and mass. Space-based solar power beaming — transmitting power from a dedicated solar collection platform to the compute node via microwave or laser — is a candidate primary power architecture [3,15].
The rectenna (rectifying antenna) system at the receiving node converts incident microwave energy to DC power through a diode rectification process. The diode rectification process produces switching transients in the DC conversion stage at the rectenna's fundamental switching frequency f_switch and its harmonics. For a bridge rectifier topology, the fundamental harmonic is at twice the microwave carrier frequency divided by the rectification stage count — typically 5-20 kHz for practical implementations. The harmonic series extends to several MHz.
This harmonic current injection into the DC bus is an EM interference source independent of the training load. It was identified as a design gap in the HERALD architecture because:
- The rectenna switching frequency (5-20 kHz) falls above the 0.1-1 Hz training burst frequency and below the magnetometer sampling rate (typically 1-10 Hz), placing its fundamental frequency in a band not covered by the training load model.
- The harmonic content is broadband and stochastic, unlike the predictable training burst current profile. It cannot be fed forward from the job queue and must be estimated from measurements.
- Rectenna current amplitude is proportional to received power, which varies with pointing accuracy, beam path geometry, and atmospheric conditions. It is not predictable from the training scheduler alone.
5.2 Harmonic Separation Filter
The HERALD state vector includes I_rect and H_rect to enable real-time estimation of the rectenna contribution to M_spurious. The harmonic separation filter decomposes the total measured bus current into training load and rectenna components:
I_bus(t) = I_train(t) + I_rect(t) + I_noise(t) (24)
The training load component I_train(t) is predicted by the HERALD dispatch algorithm and is known a priori. The rectenna component I_rect(t) has a known spectral structure (harmonics at f_switch, 2f_switch, 3f_switch, ...) but unknown amplitude. The separation is achieved by bandpass filtering:
I_rect(t) = Σₖ aₖ · cos(2πk·f_switch·t + φₖ) (25)
where aₖ and φₖ are the amplitude and phase of the k-th harmonic, estimated by the Kalman filter from the magnetometer measurements using the known harmonic structure as a constraint.
The requirement on the LC filter between the rectenna and the DC bus prevents rectenna switching harmonics from propagating to the magnetorquer control bandwidth (DC to 100 Hz):
f_cutoff,LC ≤ 50 Hz (at least 100× below minimum switching frequency of 5 kHz) (26)
The LC filter reduces the rectenna harmonic current injection into the attitude-relevant band by approximately 60 dB (factor of 1,000 in current amplitude) for the fundamental harmonic and more for higher harmonics. The residual rectenna contribution below the filter cutoff is modeled in the Kalman state as a slowly-varying DC term.
6. MULTI-NODE PLASMA PHASED-ARRAY COORDINATION
6.1 Fleet-Scale Electromagnetic Coordination
For orbital AI compute fleets consisting of multiple nodes operating in formation, each node's magnetoplasma thruster ring — used for both station-keeping and active particle shielding during solar energetic particle events — represents an additional electromagnetic coupling between nodes. Independent operation of plasma emission systems across a multi-node fleet creates potential for standing wave interference patterns in the combined magnetic field geometry that can focus charged particles toward the fleet rather than deflecting them. We term this the anti-trap requirement.
The anti-trap condition requires coordinated phase assignment across all fleet nodes. HERALD extends to handle this coordination as a fourth scheduling output, alongside burst throttling, attitude coupling, and gradient staleness management.
6.2 The Anti-Trap Phase Assignment Problem
For a fleet of N nodes with plasma emission systems, the combined magnetic field at position r is:
B_total(r,t) = Σᵢ₌₁ᴺ Bᵢ(r) · cos(ωt + φᵢ) (27)
where Bᵢ(r) is the field contribution from node i at position r, ω is the plasma oscillation frequency, and φᵢ is the phase offset for node i. The anti-trap condition requires:
∇B_total · r̂_outward > 0 for all threat directions (28)
B_total(r_inter-node) > B_total(r_node) for all node positions (29)
Condition (28) ensures the magnetic field gradient points outward — deflecting incoming charged particles away from the fleet. Condition (29) ensures there is no magnetic saddle point between nodes that would channel radiation toward the fleet center.
For a symmetric N-node fleet in a regular geometric arrangement, the optimal phase assignment satisfying both conditions is the uniform phase distribution:
φᵢ* = (2π/N) · i for i = 0, 1, ..., N-1 (30)
This distributes the phase uniformly around the unit circle, producing a combined field geometry with outward-pointing gradient in all directions and no inter-node saddle points.
For asymmetric fleet geometries — irregular spacing, different node power levels, or degraded nodes — HERALD solves the phase assignment as a real-time convex optimization:
{φᵢ*} = argmin Σⱼ∈threat_directions max(0, −∇B_total(rⱼ) · r̂ⱼ) (31)
The objective minimizes the number of threat directions with inward-pointing field gradients. This is a convex problem in the phase variables {φᵢ} when the node positions and field models are fixed, solvable by standard interior-point methods in under 10 ms on modest hardware.
6.3 Storm Mode Integration
During a Carrington-level solar energetic particle event (proton fluence > 10^10 cm^(-2)/min), the plasma phased-array switches to maximum-power collective shielding. HERALD suspends all non-critical compute to maximize available bus current for plasma emission. The coupled optimization objective (equation 23) gains a third term:
J_storm = min w₁·staleness + w₂·delay + w₃·(1 − P_shield) (32)
where P_shield is the shielding effectiveness of the current plasma configuration. During storm mode, w₃ >> w₁, w₂ — shielding takes priority over training throughput. HERALD enforces this priority shift as a constitutional scheduling constraint: no training burst may be initiated during storm mode if it would reduce available plasma bus current below the minimum shielding threshold.
The plasma bus and compute bus are electrically isolated via per-node galvanic optical isolators, ensuring that storm-mode plasma priority does not interact with the dI/dt attitude control constraint. The two constraints operate on independent electrical subsystems.
7. DISCUSSION AND IMPLEMENTATION CONSIDERATIONS
7.1 Training Throughput Cost of the dI/dt Constraint
The dI/dt constraint defers training burst initiation when the predicted current ramp would violate equation (12). The throughput cost depends on the ratio of the unconstrained dI/dt to the constraint threshold dI/dt|_max and the characteristic burst initiation time T_ramp.
For a representative 40 MW platform with centralized bus architecture (dI/dt|_max = 200 A/s from Table 1) and a cluster of 10,000 GPUs at 400 V with T_ramp = 500 ms (hardware-limited ramp time), the unconstrained dI/dt is:
dI/dt_unconstrained = ΔI_cluster / T_ramp = 50,000 A / 0.5 s = 100,000 A/s (33)
The constraint requires dI/dt ≤ 200 A/s, so the burst must be initiated over a ramp time of:
T_ramp,constrained = ΔI_cluster / dI/dt|_max = 50,000 / 200 = 250 s (34)
This 250-second ramp time, compared to the unconstrained 0.5 seconds, represents a significant change in burst initiation protocol. However, the throughput impact depends on burst frequency, not ramp time per burst. For training jobs with characteristic duration of hours to days, a 250-second ramp constitutes less than 1% of total job time and does not meaningfully reduce training throughput. For very short jobs (duration < 10 minutes), the ramp time becomes a meaningful fraction of job time. The HERALD scheduler should prioritize short jobs for initiation during periods when the current budget is already near target level (low dI/dt required), and defer long burst initiations to periods of low residual current.
Switching to distributed bus architecture (A_eff < 2 m² per rack) increases dI/dt|_max by a factor of 10 (equation 12), reducing the constrained ramp time to 25 seconds. This is a substantial improvement and is the primary motivation for recommending distributed bus topology for orbital compute platforms above approximately 10 MW.
7.2 Bus Topology Decision
The interference analysis establishes a clear preference ordering for bus topology:
- Centralized DC bus (single high-current feeder): A_eff is maximized (~20 m²), constraint is tightest (dI/dt|_max ~200 A/s), ramp time penalty is largest. Not recommended above 10 MW.
- Distributed per-rack feeders: A_eff is minimized (<2 m²), constraint is relaxed 10×, ramp time penalty is manageable. Preferred architecture above 10 MW.
- High-voltage DC (HVDC) at 4 kV: Reduces cluster current by 10× at the same power level (I_cluster = P/V_bus), reducing M_spurious by 10×. Compatible with either topology; recommended for platforms above 40 MW.
7.3 Correlated SEP Event Failure Mode
The HERALD Kalman filter assumes that measurement noise terms w_I and v_B are independent. During a Carrington-level solar energetic particle event, this assumption fails: correlated multi-bit upsets across sensor nodes produce correlated measurement errors that violate the independence assumption and cause the Kalman filter to diverge.
The mitigation is a storm-mode protocol that switches the Kalman measurement update from digital sensor readings to analog majority-voting photonic inter-die links during declared storm conditions. Photonic interconnects are immune to charge deposition from ionizing particles — the light signal propagating in a silicon waveguide is not affected by electron-hole pair generation in the surrounding silicon. This provides a radiation-immune measurement channel that maintains filter observability during the worst SEP events.
Storm mode is declared when the particle flux sensor network detects fluence rate exceeding 10^8 cm^(-2)·s^(-1) — a threshold that provides approximately 10-30 minutes of warning before a Carrington-class event reaches peak intensity [16].
7.4 Single-Event Latch-Up on Shared Bus
A heavy-ion strike on a power FET in a shared DC bus node can latch the affected FET into a high-current state, injecting a large current spike into the bus. This is a radiation-induced latch-up (SEL) event [17]. On a shared bus, the latch-up current spike propagates to all nodes, generating a spurious magnetic dipole moment substantially larger than a training burst.
Mitigation requires per-node galvanic isolation with optical triggering: each node's connection to the shared bus passes through a solid-state switch with optical control signal. A SEL detection circuit monitors each node's bus current and opens the isolation switch within microseconds of detecting anomalous current. The optical triggering ensures that a SEL event on one node's electronics cannot propagate a false trip signal to other nodes through a shared control bus.
7.5 Optical Attitude Determination and Thruster-Based Correction as a HERALD Backup Layer
The HERALD architecture described in this paper addresses the compute-attitude coupling problem primarily through prevention: by constraining the rate of current change at burst initiation, the scheduler eliminates the dominant source of spurious magnetic disturbance before it reaches the attitude control system. Prevention, however, is not a complete solution. Residual disturbances from imperfect current prediction, rectenna harmonic injection, and unmodeled coupling pathways will produce small but non-negligible attitude errors that the magnetorquer system alone may be too slow to correct within mission pointing requirements.
A complementary second layer — optical attitude determination combined with fast thruster correction — addresses this residual error. Where HERALD prevents the disturbance through scheduling, optical sensors detect any residual attitude deviation that propagates through the prevention layer, and microthrusters correct it on a timescale shorter than the magnetorquer response window.
Three optical sensor modalities are relevant to this backup layer. Star trackers provide the highest accuracy attitude reference available without ground contact, achieving 1-5 arcsecond pointing knowledge under nominal operating conditions and maintaining function through periods of elevated solar activity with appropriate radiation hardening. Sun sensors provide a coarser but highly radiation-tolerant backup reference, particularly valuable during solar energetic particle events when star tracker performance may degrade. Horizon sensors close the low-frequency drift loop, compensating for the slow secular attitude errors that accumulate between burst events and are below the detection threshold of the HERALD Kalman filter's attitude state vector.
Thruster correction authority requires two complementary technologies operating at different timescales. Cold-gas thrusters provide the fast impulsive correction needed to arrest attitude excursions within the burst initiation window — their response time of tens of milliseconds is an order of magnitude faster than magnetorquer desaturation. Their limitation is delta-V budget: cold-gas systems carry finite propellant and cannot sustain continuous correction over century-scale mission durations. Ion microthrusters address the sustained correction requirement, providing precise low-thrust attitude maintenance with orders-of-magnitude better specific impulse than cold-gas at the cost of slower transient response. The optimal architecture combines both: cold-gas for fast burst-induced correction events, ion thrusters for continuous low-level drift compensation between events.
Together, this optical-plus-thruster backup layer and the HERALD scheduling constraint form a two-layer attitude stability architecture. The first layer eliminates the dominant disturbance source through prediction and constraint. The second layer detects and corrects whatever the first layer misses. The combination provides attitude stability guarantees robust to both predicted electromagnetic disturbances and the residual unmodeled coupling that any real megawatt-scale orbital compute platform will inevitably produce.
8. CONCLUSION
We have identified and formalized a previously uncharacterized coupling between machine learning compute schedulers and orbital attitude control systems that becomes design-critical at megawatt-scale orbital compute platform power levels. The coupling mechanism — training burst events on the platform DC bus generating spurious magnetic dipole moments that compete with magnetorquer attitude control authority — is invisible to standard spacecraft EMI analysis, which does not address disturbance torques generated by internal current distribution at sub-10 Hz frequencies.
The derived interference threshold (equation 12) provides a quantitative scheduling constraint: dI/dt ≤ ε_int · M_auth / (A_eff · τ_control). At representative 40 MW platform parameters with centralized bus architecture, this constraint requires training bursts to ramp current at 200 A/s or less — a factor of 500 more slowly than hardware-limited burst initiation.
The HERALD architecture enforces this constraint through a co-designed Kalman filter that jointly estimates attitude dynamics and bus current dynamics, driving a constrained scheduler that defers burst initiation to the earliest feasible time within the constraint envelope. The HERALD framework extends to handle power beaming rectenna harmonic interference through a harmonic separation filter and LC filter specification, and to coordinate multi-node plasma phased-array fleet shielding through a convex phase assignment optimization. These extensions address electromagnetic interference sources beyond the training load that were not accounted for in any prior orbital platform design.
Three design recommendations emerge from this work for orbital compute platform architects:
- Adopt distributed per-rack bus topology for any platform above 10 MW. Centralized bus architecture makes the attitude coupling unmanageable at higher power levels without unacceptable throughput penalties.
- Co-design the ML training scheduler and the attitude control system from the beginning of the platform design process. Retrofitting scheduler-side current smoothing after the attitude control system is designed will produce suboptimal solutions because the constraint parameters depend on bus geometry decisions made during attitude control system design.
- Include power beaming rectenna switching harmonics in the electromagnetic compatibility specification from the first design review. The rectenna is an independent interference source that the training scheduler cannot predict or control.
The HERALD problem — compute affecting physics affecting mission safety — is likely to recur as orbital compute platforms scale. The framework introduced here provides a template for identifying and formalizing such cross-domain couplings before they become mission failures.