IMF Reconstruction using Machine Learning:

Proof of Concept

Ramiz A. Qudsi

qudsira@bu.edu

(https://qudsiramiz.github.io)

https://slides.com/qudsi/magnetore_bu

Multiscale nature of IMF

Coleman (ApJ, 1968)

Bale et al. (Natur, 2019)

Competing Theories of Turbulence

Critical effects of magnetic structure on plasma dynamics:

  • Propagation of energetic particles
  • 3D context of kinetic processes like magnetic reconnection

How do we get the full 3D structure of magnetic field?

Put a magnetometer at every point in space

Launch a constellation of spacecraft with magnetometers on them

Compute values at rest of the points by interpolation

Conventional interpolation are not best suited

Gaussian Processes in Machine Learning

Hub

Spacecrafts

$$d_i$$

22

Some relevant scales:

d_i \mathrm{(1AU)} \sim 100 ~\mathrm{km} \\ V_{sw} \sim 500 ~\mathrm{km/sec} \\ X_{sim} \mathrm{(box size)} \sim 40 d_i\\ \hspace{3.1cm}\sim 4\times 10^3 \mathrm{km}\\ d_{spc} \sim [1, 11]~d_i\\ \hspace{1.2cm}\sim \mathrm{ [10^2, 10^3]~ km} \\ f[min,max] \sim V_{sw}/(2\times d_{spc, max, min}) \\ \hspace{0.8cm}\sim [0.25, 2.5]~\mathrm{Hz}\\ \nu (\mathrm{sampling~rate}) \sim 40 \mathrm{Hz}\\ \mathrm{(just~because~of~box~size~and~stuff)}\\ \mathrm{time~to~move~across~the~simulation~box} \sim 10 ~\mathrm{sec}

Gaussian Process Regression Interpolation

Gaussian Process Regression:

$$m(\mathbf{x}) = \mathbb{E}[f(\mathbf{x})]$$

Mean function

$$k(\mathbf{x}, \mathbf{x'}) = \mathbb{E}[(f(\mathbf{x}) - m(\mathbf{x}))(f(\mathbf{x'}) - m(\mathbf{x'}))]$$

Covariance function

$$f(\mathbf{x}) \sim \mathcal{GP}\left(m(\mathbf{x}), k(\mathbf{x}, \mathbf{x'})\right)$$

Gaussian Processes

Constant

Linear

RBF

Matern

Kernels

Results

One component at a time

Maruca 2021, Frontiers in A & SS

All components at a time

Qudsi et. al., in prep

Single component input

Vector input

Different configurations

simulation data

8 spc

24 spc

8 spc

24 spc

8 spc

16 spc

24 spc

34 spc

24 spc*

22 spc*

Moving Forward

What's next?

Explore the impact of relative spacecraft separation and arrangement

Incorporate the zero divergence condition on magnetic field

Explore alternative algorithms/methods

Deep Gaussian Processes

 

Neural Network as GP

Find ways to compare the reconstructed images for different styles

Thank You!

https://slides.com/qudsi/magnetore_bu

Magnetore

By Ramiz Qudsi

Magnetore

Discussion of work done towards MagneToRE: the Magnetic Topology Reconstruction Explorer

  • 176