Schwarzschild light bending

A horizontal plane wave of photons enters from the left and meets a non-rotating black hole at the origin. Geometric units G = c = M = 1. Each photon obeys a null geodesic with conserved energy E and angular momentum L; its fate is fixed by the impact parameter b = L / E. Photons with |b| < 3√3 ≈ 5.196 cross the photon sphere at r = 3 and are swallowed (red); photons with |b| > 3√3 are deflected (blue), with photons just above critical looping the photon sphere multiple times before escaping. Drag the sliders to vary the photon count and the impact-parameter range.

Figure 1. Plane wave of null geodesics in the Schwarzschild equatorial plane (M = 1). Method: per-photon velocity-Verlet on the radial Hamiltonian H = p_r²/2 + L²(1−2/r)/(2r²) from shared/js/engine/symplectic.js; angular coordinate φ advances as L/r².

N 41

b_max 9.00

Show b_crit = 3√3 reference lines and critical photon

WHAT TO TRY

Vary each control and watch the rail readouts respond.
Compare the diagnostic plot against the live scene.