Mars lacks a substantial magnetic field; as a result, the solar wind ablates the Martian atmosphere, making the surface uninhabitable. Therefore, any terraforming attempt will require an artificial Martian magnetic shield. The fundamental challenge of building an artificial magnetosphere is to condense planetary-scale currents and magnetic fields down to the smallest mass possible. Superconducting electromagnets offer a way to do this. However, the underlying physics of superconductors and electromagnets limits this concentration. Based upon these fundamental limitations, we show that the amount of superconducting material is proportional to $B_c^{-2}a^{-3}$, where $B_c$ is the critical magnetic field for the superconductor and $a$ is the loop radius of a solenoid. Since $B_c$ is set by fundamental physics, the only truly adjustable parameter for the design is the loop radius; a larger loop radius minimizes the amount of superconducting material required. This non-intuitive result means that the "intuitive" strategy of building a compact electromagnet and placing it between Mars and the Sun at the first Lagrange point is unfeasible. Considering reasonable limits on $B_c$, the smallest possible loop radius is $\sim$10 km, and the magnetic shield would have a mass of $\sim 10^{19}$ g. Most high-temperature superconductors are constructed of rare elements; given solar system abundances, building a superconductor with $\sim 10^{19}$ g would require mining a solar system body with several times $10^{25}$ g; this is approximately 10% of Mars. We find that the most feasible design is to encircle Mars with a superconducting wire with a loop radius of $\sim$ 3400 km. The resulting wire diameter can be as small as $\sim$5 cm. With this design, the magnetic shield would have a mass of $\sim 10^{12}$ g and would require mining $\sim 10^{18}$ g, or only 0.1\% of Olympus Mons.

We explore a new simulation scheme for partial differential equations (PDE's) called Information Field Dynamics (IFD). Information field dynamics attempts to improve on existing simulation schemes by incorporating Bayesian field inference, which seeks to preserve the maximum amount of information about the field being simulated. The field inference is truly Bayesian and thus depends on a notion of prior belief. Here, we analytically prove that a restricted subset of simulation schemes in IFD are consistent, and thus deliver valid predictions in the limit of high resolutions. This has not previously been done for any IFD schemes. This restricted subset is roughly analogous to traditional fixed-grid numerical PDE solvers, given the additional restriction of translational symmetry. Furthermore, given an arbitrary IFD scheme modelling a PDE, it is a-priori not obvious to what order the scheme is accurate in space and time. For this subset of models, we also derive an easy rule-of-thumb for determining the order of accuracy of the simulation. As with all analytic consistency analysis, an analysis for nontrivial systems is intractable, thus these results are intended as a general indicator of the validity of the approach, and it is hoped that the results will generalize.

In this study we explore a new simulation scheme for partial differential equations known as Information Field Dynamics (IFD). Information field dynamics attempts to improve on existing simulation schemes by incorporating Bayesian field inference into the simulation scheme. The field inference is truly Bayesian and thus depends on a notion of prior belief. A number of results are presented, both theoretical and practical. Many small fixes and results on the general theory are presented, before exploring two general classes of simulation schemes that are possible in the IFD framework. For both, we present a set of theoretical results alongside the development of a prototype scheme. The first class of models corresponds roughly to traditional fixed-grid numerical PDE solvers. The prior Bayesian assumption in these models is that the fields are smooth, and their correlation structure does not vary between locations. For these reasons we call them translation-invariant schemes. We show the requirements for stability of these schemes, but most importantly we prove that these schemes indeed converge to the true behaviour of the field in the limit of high resolutions. Convergence had never been shown for any previous IFD scheme. We also find the error scaling of these codes and show that they implement something very analogous to a high-order finite-difference derivative approximation, which are the most elementary and well-studied numerical schemes. This is an important result, which proves the validity of the IFD approach. The second class of schemes, called the SPH-like schemes are similar to existing Smooth Particle Hydrodynamics codes, in which the simulation grid moves with the flow of the field being modelled.