Past planetary missions have been broad and detailed for Gas Giants, compared to flyby missions for Ice Giants. Presently, a mission to Neptune using electrodynamic tethers is under consideration due to the ability of tethers to provide free propulsion and power for orbital insertion as well as additional exploratory maneuvering --- providing more mission capability than a standard orbiter mission. Tether operation depends on plasma density and magnetic field $\mathbf{B}$, though tethers can deal with ill-defined density profiles, with the anodic segment self-adjusting to accommodate densities. Planetary magnetic fields are due to currents in some small volume inside the planet, magnetic-moment vector, and typically a dipole law approximation --- which describes the field outside. When compared with Saturn and Jupiter, the Neptunian magnetic structure is significantly more complex: the dipole is located below the equatorial plane, is highly offset from the planet center, and at large tilt with its rotation axis. Lorentz-drag work decreases quickly with distance, thus requiring spacecraft periapsis at capture close to the planet and allowing the large offset to make capture efficiency (spacecraft-to-tether mass ratio) well above a no-offset case. The S/C might optimally reach periapsis when crossing the meridian plane of the dipole, with the S/C facing it; this convenient synchronism is eased by Neptune rotating little during capture. Calculations yield maximum efficiency of approximately 12, whereas a $10^{\circ}$ meridian error would reduce efficiency by about 6%. Efficiency results suggest new calculations should be made to fully include Neptunian rotation and consider detailed dipole and quadrupole corrections.

The need for the extra dimension in Kustaanheimo-Stiefel (KS) regularization is explained by the topology of the Hopf fibration, which defines the geometry and structure of KS space. A trajectory in Cartesian space is represented by a four-dimensional manifold, called the fundamental manifold. Based on geometric and topological aspects classical concepts of stability are translated to KS language. The separation between manifolds of solutions generalizes the concept of Lyapunov stability. The dimension-raising nature of the fibration transforms fixed points, limit cycles, attractive sets, and Poincar\'e sections to higher-dimensional subspaces. From these concepts chaotic systems are studied. In strongly perturbed problems the numerical error can break the topological structure of KS space: points in a fiber are no longer transformed to the same point in Cartesian space. An observer in three dimensions will see orbits departing from the same initial conditions but diverging in time. This apparent randomness of the integration can only be understood in four dimensions. The concept of topological stability results in a simple method for estimating the time scale in which numerical simulations can be trusted. Ideally all trajectories departing from the same fiber should be KS transformed to a unique trajectory in three-dimensional space, because the fundamental manifold that they constitute is unique. By monitoring how trajectories departing from one fiber separate from the fundamental manifold a critical time, equivalent to the Lyapunov time, is estimated. These concepts are tested on N-body examples: the Pythagorean problem, and an example of field stars interacting with a binary.

The quark mass dependences of light element binding energies and nuclear scattering lengths are derived using chiral perturbation theory in combination with non-perturbative methods. In particular, we present new, improved values for the quark mass dependence of meson resonances that enter the nuclear force. A detailed analysis of the theoretical uncertainties arising in this determination is presented. As an application we derive from a comparison of observed and calculated primordial deuterium and helium abundances a stringent limit on the variation of the light quark mass, $\delta m_q/m_q = 0.02 \pm 0.04$. Inclusion of the neutron lifetime modification under the assumption of a variation of the Higgs vacuum expectation value that translates into changing quark, electron, and weak gauge boson masses, leads to a stronger limit, $|\delta m_q/m_q| < 0.009$.