Study of Jupiter's Interior with Quadratic Monte Carlo Simulations

We construct models for Jupiter's interior that match the gravity data obtained by the Juno and Galileo spacecraft. To generate ensembles of models, we introduce a novel quadratic Monte Carlo technique, which is more efficient in confining fitness landscapes than the affine invariant method that relies on linear stretch moves. We compare how long it takes the ensembles of walkers in both methods to travel to the most relevant parameter region. Once there, we compare the autocorrelation time and error bars of the two methods. For a ring potential and the 2d Rosenbrock function, we find that our quadratic Monte Carlo technique is significantly more efficient. Furthermore, we modified the walk moves by adding a scaling factor. We provide the source code and examples so that this method can be applied elsewhere. Here we employ our method to generate five-layer models for Jupiter's interior that include winds and a prominent dilute core, which allows us to match the planet's even and odd gravity harmonics. We compare predictions from the different model ensembles and analyze how much an increase in the temperature at 1 bar and ad hoc change to the equation of state affect the inferred amount of heavy elements in the atmosphere and in the planet overall.