AETHERIA

Project Title: AETHERIA Real-Time GPU Gravitational Dynamics, Observational Photometry & Sensor-Fused Kinematic Workstation 1. Executive Summary Investigating major galactic mergers and non-equilibrium secular evolution typically requires batched execution on High-Performance Computing (HPC) clusters or reliance on static, pre-rendered kinematic visualizers. AETHERIA is an interactive astrophysics simulation and observational diagnostics workstation developed for the YPAE Space Simulation Hackathon.

Engineered from scratch using Python, Taichi Lang (GPU-JIT compiled kernels), and OpenCV, AETHERIA resolves gravitational equations of motion for a 40,000-particle composite system at 60 FPS on consumer-grade hardware. It directly links dynamic collision mechanics with an empirical observational pipeline—incorporating logarithmic surface brightness profiling (μ(R)), multi-phase interstellar medium (ISM) extinction models, a live core distance scale, and a dual-zone optical hardware sensing interface architected for future expansion into an IMU-instrumented cyber-physical sensor glove.

Theoretical & Practical Motivation 1. The HPC Latency Bottleneck in Pedagogy: Educational and analytical workflows frequently decouple orbital dynamics from observational verification. Simulating tidal morphology across multi-gigayear timescales on supercomputing clusters prevents real-time parameter sweeps during lectures and laboratory sessions. 2. Absence of Synthetic Observational Diagnostics: Many interactive particle visualizers treat stars as massless, uniform point emitters, ignoring physical observables such as the Sérsic profile index (n), effective half-light radius (Reff ), differential tidal shear, and optical depth (τ λ). 3. Immersive Physical Interfacing: Traditional desktop peripherals (mouse and keyboard) impose rigid planar constraints. Modern planetariums, science centers, and academic presentations require intuitive, multi-degree-of-freedom spatial interfaces to manipulate celestial mechanics directly.

Deep Astrophysical Implementation A. Gravitational Kinematics & Regularized Potential Solvers The simulation solves stellar orbital trajectories within the joint gravitational potential of two Supermassive Black Holes (SMBHs) undergoing mutual orbital decay. To prevent unphysical numerical singularities during close perigalacticon passages, point-source acceleration uses a regularized Plummer-softening kernel: a_i = - SUM_{j=0}^{1} [ (G * M_j) / ( ||r_i - r_j||^2 + epsilon^2 )^(3/2) ] * (r_i - r_j) Softening Parameter: epsilon = 0.045 kpc, maintaining collisionless stellar dynamics while mitigating numerical scattering artifacts. Symplectic Integration: Particle trajectories advance through a customized GPU Euler-Verlet integration scheme, tracking the conservation of energy and angular momentum across inbound parabolic plunges, tidal disruption, and post-merger relaxation. Initial Kinematics: Disk particles are seeded in centrifugal equilibrium governed by circular orbital velocity: v_circ(r) = sqrt( (G * M_SMBH) / (r + r_soft) ). Tidal Tensor & Gravitational Shear: Instantaneous tidal strain on orbiting disks is monitored via the tidal field gradient: T_ij ≈ d^2(Phi) / (dx_i * dx_j) ~ (G * M_SMBH) / R^3 This calculates when the tidal boundary is breached and resonant gravitational torques pull stars into extended tidal tails. Collision Stage Progression: [Inbound Plunge] --> [Perigalacticon Flyby] --> [Tidal Tail Ejection] --> [Apogalacticon Turnaround] --> [Coalescence]

B. Stellar Demographics & Radiative Transfer Approximation Rather than populating galaxies with uniform particles, AETHERIA instantiates 40,000 distinct bodies categorized into three astrophysical populations:

Nuclear Bulge Population (22%): Dynamically hot, pressure-supported old stars (M stars ​ ∼1.5M_stars ) exhibiting steep velocity dispersion profiles and high central density, radiating in warm yellow-orange wavelengths (T eff ∼4,500 K). Starburst Disk Arms (50%): Rotationally supported young stellar clusters and OB associations (Teff >10,000 K) defined by logarithmic spiral density functions: r(theta) = r_0 * exp( theta / tan(psi) ) subject to differential rotation and tidal torque.Interstellar Dust Lanes (28%): Cold molecular gas and dust clouds that do not emit light, but attenuate background flux according to an empirical optical depth model: I_observed = I_intrinsic * exp( -tau_lambda ) Simulating extinction lanes and interstellar reddening, cold dust particles cast absorption signatures across the bright stellar background.C. Observational Photometry & Morphology MetricsLogarithmic Surface Brightness Profile mu(R): Evaluates projected luminous flux across 32 logarithmic radial annuli (0.015 kpc<= R <= 1.5 kpc): mu(R) = -2.5 * log10( Sigma(R) ) + C [mag / arcsec^2] This captures the morphological transformation from an exponential disk profile (I(R) = I_0 \exp(-R/h)$) into a collisionally relaxed, pressure-supported spheroid governed by the Sérsic Profile / de Vaucouleurs' law (n = 4): I(R) = I_e * exp( -b_n * [ (R / R_e)^(1/n) - 1 ] ) Dynamic Effective Radius ($R_{\text{eff}}$): Computes the half-light radius enclosing 50% of total cumulative stellar radiance in real time, charting spatial expansion driven by orbital angular momentum redistribution during tidal ejection.Direct SMBH Separation Scale Ruler: Computes instantaneous Euclidean separation in kiloparsecs ($\text{kpc}$) and relative orbital velocity ($V_{\text{rel}}$ in $\text{km/s}$), rendering a calibrated metric ruler with transverse ticks directly between both SMBH attractors in the 3D viewport. Non-Invasive Optical Tracking Engine (Current Implementation) AETHERIA converts standard monocular camera hardware into a real-time spatial telemetry tracker without requiring cloud inference or heavyweight external neural network runtimes:

Dual-Color-Space Skin Isolation: Merges normalized HSV and YCrCb chrominance/luminance thresholding to deliver stable segmentation under varying ambient illuminations: Mask = (H_min <= H <= H_max) AND (Cr_min <= Cr <= Cr_max) AND (Cb_min <= Cb <= Cb_max) Adaptive Spatial Cone Masking: Dynamically suppresses facial centroids within the central upper bounding cone (x∈[0.25W,0.75W],y∈[0,0.5H]), preventing background interference. Topological Convexity & Solidity Analysis: Computes contour hulls and defect depths: Solidity = Area_contour / Area_hull Distinguishes an **Open Palm** ($\text{Solidity} < 0.80, \text{defects} \ge 2 \rightarrow \text{Simulation Run}$) from a **Closed Fist** ($\text{Solidity} \ge 0.80 \rightarrow \text{Simulation Hold}$). A spatial boundary threshold across the upper 25% of the left frame triggers an automated epoch reset.Quadratic Deadzone Rate Vectoring: Evaluates right-hand centroid displacement (Delta x, Delta y) from the calibrated neutral anchor, filtering sensor jitter through a 15% deadzone and mapping output via a polynomial transfer curve: f(Delta) = sign(Delta) * |Delta|^1.3

Impact & Future Roadmap A. Pedagogical & Open-Science Impact 1. Accessible Supercomputing: Replaces expensive cluster queues with an open-source, 60 FPS GPU-native simulation framework that runs on commodity laptops, democratizing computational astrophysics. 2. Turnkey Museum & Planetarium Installation: Provides an interactive kiosk configuration for public outreach centers, allowing users to physically drive galactic interactions and see how parameters shape final morphologies.

Cyber-Physical Sensor Glove Integration (Next Milestone) While the camera pipeline offers immediate, zero-cost accessibility, our next hardware phase bridges optical tracking with a custom embedded sensor glove. Sub-Millimeter Micro-Gesture Parameter Modulation:

Flexion-Driven Sweep Control: Individual finger contractions map directly to scalar parameters (e.g., index-finger flexion scales impact velocity vimp , thumb-middle pinch adjusts mass ratio M2 /M 1). 6-DoF Core Manipulation: Allows researchers to "grab" individual galactic nuclei, altering collision inclination angles (i) and impact vectors via wrist pitch and yaw. Closed-Loop Gravitational Haptic Feedback: Mounts Linear Resonant Actuators (LRAs) on glove fingertips to deliver tactile feedback proportional to calculated tidal strain forces ($T_{ij}$) and gravitational acceleration during perigalacticon passages.

Computational & Simulation RoadmapBarnes-Hut Hierarchical Octree: Transitioning from the current dominant-attractor model to a fully self-gravitating Barnes-Hut octree (O(N \log N)) using Taichi sparse memory allocators (ti.root.pointer), scaling active particles to $N > 200,000. Smoothed Particle Hydrodynamics (SPH): Implementing density kernels (W(r, h)) and Navier-Stokes pressure solvers to simulate interstellar gas shock-heating, starburst initiation, and ram-pressure stripping. Astronomical FITS Pipeline: Adding export utilities to package volumetric data into standardized FITS files with synthetic WCS headers, enabling direct comparison with observations from the Hubble and James Webb Space Telescopes.