This image shows the electric current density in the orbital plane of a binary neutron star merger about 2 milliseconds after the stars have come into contact. During this phase of the evolution, matter in the neutron stars is violently stirred by the Kelvin-Helmholtz instability generating ultra-strong magnetic fields of up to 10^17 Gauss (roughly 100 quadrillion times the Earth's magnetic field). This simulation was run on Aurora using the team's INCITE allocation. Image: Dr. Eduardo Mario Gutiérrez, Penn State.