Design of a system for controlling a levitating sphere in superfluid 3He at extremely low temperatures

We present a new mechanical probe to study the properties of superfluid 3He at microkelvin temperatures down to 100 μK. The setup consists of a set of coils for levitating a superconducting sphere and controlling its motion in a wide variety of regimes. In particular, the realisation of motion of a levitating body at a uniform velocity presents both an experimental challenge and a promising direction into the study of the edge states in topological superfluid 3He-B. We include the theoretical study of the device stability and simulations to illustrate the capabilities of the control system.


Design of a system for controlling a levitating sphere in superfluid 3 He at extremely low temperatures
Manuel Arrayás 1* , José L. Trueba 1 , Carlos Uriarte 1 & Dmitry E. Zmeev 2 We present a new mechanical probe to study the properties of superfluid 3 He at microkelvin temperatures down to 100 μK. The setup consists of a set of coils for levitating a superconducting sphere and controlling its motion in a wide variety of regimes. In particular, the realisation of motion of a levitating body at a uniform velocity presents both an experimental challenge and a promising direction into the study of the edge states in topological superfluid 3 He-B. We include the theoretical study of the device stability and simulations to illustrate the capabilities of the control system.
Starting from Andronikashvili's experiment with the torsional pendulum 1 , the measurement of drag forces on oscillating mechanical resonators immersed in liquid helium has been one of the main tools in the study of its superfluid properties. On resonance, the drag force is equal and opposite in sign to the driving force, applied to the oscillator, and is known to the experimenter. Multiple properties of the superfluid can be successfully inferred from the behaviour of the drag force. On the other hand, oscillatory motion necessarily imposes significant limits on the experimental capabilities. Certain phenomena only manifest when a mechanical probe moves uniformly in the superfluid and cannot be observed using oscillators. Examples include generation of homogeneous and isotropic turbulence produced by a tightly-fitted grid towed through a channel [2][3][4] and observation of superfluid currents beyond the Landau velocity 5 .
Despite being the simplest form of motion, uniform motion is notoriously difficult to implement in a lowtemperature experiment, mainly due to the limitations on heat release and available space 2,3,6 . Another obstacle to realisation of a successful experiment with a uniformly moving body is that the drag force acting upon it is a lot more difficult to measure than the drag force on an oscillator in resonance 6 .
Recently, the uniform motion technique was used to discover anomalously long lifetimes of quasiparticles bound to the surfaces of superfluid 3 He-B 7 . This time was found to be about 6 ms, much longer than any other time scales in the system 8 . The discovery of the long lifetime of the quasiparticles also explains another recently observed anomalous phenomenon-exceptionally high apparent thermal conductivity in narrow channels of superfluid 3 He: most likely it is the long-lived surface quasiparticles that carry significant energy along the walls of the channel 9 . In the experiments 5,7 , the surface-bound quasiparticles could be promoted into bulk by accelerating a goalpost-shaped superconducting wire, with the crossbar moving parallel to itself. In this case different parts of the posts moved at different velocities. Also, the crossbar was necessarily connected to a large 'bath' of surface-bound excitations on the walls of the experimental container and the states emptied into bulk by acceleration could be refilled from this 'bath' by diffusion. A more favourable topology and well-defined geometry are required for more quantitative experiments. One such experiment is determination of the spectrum of surface excitations, potentially containing Majorana-like fermions 10 , by accelerating the body to various velocities. An example of a probe suitable for these experiments is a levitated sphere.
In this work we present the design of an experimental setup consisting of a levitation and control system for moving a body in a wide range of motion regimes, from oscillatory to uniform, inside superfluid 3 He at temperatures down to 100 μK. We have made the design to work for a superconducting spherical object made of indium. It is important to note that a similar experiment has been already performed in 4 He, with a magnetic sphere being levitated between the plates of a superconducting niobium capacitor 11 . A large voltage was applied in order to charge the sphere, before the niobium was cooled down below its critical temperature. Then the DC voltage was switched off and oscillatory motion was induced by applying an oscillatory voltage amplitude up to a few volts. This arrangement, however, is not suitable for superfluid 3 He experiments, mainly due to the high magnetic field in the experimental region during precooling before nuclear demagnetisation: this field then will be trapped by the superconductor. Moreover one of our main goals is to be able to move the body with a uniform velocity. www.nature.com/scientificreports/ Our levitating system comprises a set of superconducting coils. The main two concentric rings at the bottom of the cell filled with 3 He produce the levitating force, which cancels gravity. The control system is provided by a set of lateral coils, which allow the control of the motion in the horizontal plane. In order to enable the control, the system must sit in the stability region of the magnetic and gravitational energies, which imposes some strict conditions on the final design. Our results prove that we can meet these conditions in a range of currents and masses of the objects which can levitate inside the experimental cell. In this work we provide such analysis of the stability conditions and perform simulations of the motion of the sphere in the superfluid.

Results
The setup. Our experimental technique is based on the setup shown in Fig. 1. Six superconducting coils are placed around the experimental cell, providing the levitation force and motion control of a superconducting sphere inside the cell. The cell, 12 mm in diameter, is filled with superfluid 3 He at temperatures below 200 μK and contains vibrating wire bolometers for detecting energy losses resulting from the motion of the sphere in the superfluid. Two coils are placed at the bottom within a nested arrangement, with their symmetry axes on the vertical direction z, giving the main contribution to the levitating force. The four lateral coils placed on the sides have their axes pointing in the x and y directions, respectively. They have equal radii of 6 mm and play the main role in the motion control of the sphere.
The inner bottom coil has a radius of R 1 = 6 mm while the outer has the radius of R 2 = 15 mm. We can create magnetic fields in opposite directions controlling the direction of the currents I 1 and I 2 so we can obtain the right levitation conditions as we will discuss in the next section.

Stability.
A superconducting body exhibits diamagnetic behaviour. Diamagnetic objects are repelled by magnetic fields and can be levitated. However, in order to achieve stable levitation, the total energy must posses a minimum at the critical point, where gravitational and magnetic forces balance out. The stability conditions, as will be shown below, impose certain constraints on our design. On the other hand, we have limitations for the maximum currents in the coils, dictated by the critical current, and for the maximum field, dictated by the critical field of the superconducting material of the levitating object, which is 28 mT in indium.
The differential amount of work done by the magnetic field on the body is −dM · B , while the induced magnetic moment reads 12 Note that the factor α appears as a geometric factor and V is the volume of the body. In the case of a superconducting spherical body 13 when the radius is bigger than the penetration length a good approximation is α = 1/2 . Due to the perfect diamagnetic behaviour, the magnetic susceptibility is χ = −1 . Integrating the magnetic work as the field increases from zero to its final value, and taking into account the gravitational energy, we get the total energy where m is the mass of the body and g is the gravitational acceleration. The points where levitation occurs are called critical points, as we have ∇U = 0 . For the equilibrium to be stable, the energy must be at a minimum, i.e. the energy will increase in all the directions from the critical points. In mathematical form, the conditions are 14 (2) U(r) = mgz + 1 4 www.nature.com/scientificreports/ We can perform an exact analysis of the stability. Up to small corrections, the field created by a coil can be approximated by the one created by a ring of radius a with n turns carrying a current I as 12 with k = 4ar/[(r + a) 2 + z 2 ] as the argument for the elliptic integrals of the first and second type and B r , B z the radial and vertical components of the field. Using this expression we can compute analytically the magnetic field created by the bottom coils. The stability conditions (3) depend only on the magnetic field configuration as we can see from (2). Thus, in order to achieve a stable levitation, we need to fulfill the condition ∇U = 0 , in the region where ∂ 2 x B 2 > 0, ∂ 2 y B 2 > 0, ∂ 2 z B 2 > 0 , here B is the modulus of the magnetic field. An important consequence of the critical condition of the force being null is that it involves only the density ρ of the levitated body as it reads We stress that a solution for an object made from a type I superconductor and reasonable sizes of coils does not exist, therefore we resort to using a light plastic sphere covered with indium, which has a relatively high critical field. The sphere is shown in the inset to Fig. 1. The coating thickness must be bigger than the penetration length which for indium is 60 nm 15 .
In Fig. 2 we have plotted again a vertical cross section showing the different stability regions when only the bottom coils are switched on. The parameters used are n = 250 turns for both coils, R 1 = 6 mm, I 1 = 2 A for the inner coil, and R 2 = 15 mm, I 2 = 1.7 A (in the opposite direction) for the outer coil. The superconducting body is a hollow sphere made of indium, with internal radius R in = 500 µ m , external radius R out = 550 µ m and a total mass of 1.26 mg. We show the unstable region in blue, while green and red represent the stable region, that is the region where ∂ 2 i B 2 > 0 . The red region is where the magnetic force in the z direction is larger than the gravitational force, while in the green region the gravitational force is larger. Thus, the interface between the green and the red regions is the critical boundary, i.e. the zero force surface. In other words, the levitated body in equilibrium will be positioned at that interface. We observe that close to the centre there is an almost horizontal flat region of stable levitation of about 14 mm which is about 13 times the diameter of the sphere, providing the possibility to control the motion of the sphere in the horizontal direction with the lateral coils.
The controlled motion of the sphere. When the sphere is displaced out of the equilibrium position, it starts to oscillate. As the viscosity is very low, we need to suppress these oscillations by introducing a control system. The control system also will allow us to move the sphere with a constant velocity, which is crucial for the studies of surface excitations in 3 He-B.
The motion around the equilibrium position in the horizontal direction are nonlinear harmonic oscillations whose frequency depends on the amplitude of the motion, but also on the intensity of the magnetic field produced by the bottom coils. In order to prevent the oscillations and to move the sphere at a (quasi-)constant velocity, we can apply a pre-calculated force 3 , by running a controlled current in the lateral coils.  www.nature.com/scientificreports/ As an example, in order to move the sphere from x 1 = 0 mm to x 2 = +2 mm, with a uniform velocity of v 0 = 0.02 m/s we need to use a time-dependent current in the lateral coils of the form: Here the current I L (t) is in amps and t is in seconds. The motion comprises three stages: acceleration, uniform motion and deceleration corresponding to the three time intervals in Eq. 7. The first stage, from t = 0 s to t 1 = 0.01 s, represents the acceleration period of time for reducing the oscillations. The second stage, from t 1 to t 2 , is the time interval in which uniform motion is induced, where t 2 = 0.12 s. The third stage, from t = t 2 to t = 0.13 s, represents the time interval where the decelerating pulse is applied to stop the motion of the sphere. The lateral coils are fed with currents producing the field in the same direction, i.e. considering the coils in the y-axis, looking in the positive direction, both coils would carry clockwise currents to move the sphere in the positive y direction. The result can be seen in Fig. 3.
We can also use a similar control to suppress the vertical oscillations. These oscillations are produced due to the inhomogeneity of the magnetic field in addition to the distortion of the stability region produced by the side coils. The vertical control is implemented by adding a small synchronised time-dependent current to the bottom coils. The current in the inner coil is shown in Eq. 8 and the current applied to the outer coil is shown in Eq. 9. We can see in Fig. 4 that the oscillations are also removed in the vertical direction.
:  Detection of the sphere position. We also need to implement the position detection system to monitor the motion of the sphere and to confirm the validity of our design. For this purpose we will utilise further two coils with axes along x and y. These detection coils will be run at an AC current with a frequency by far exceeding all characteristic mechanical frequencies. The change in the position of the diamagnetic sphere will affect the inductance of the coils which will be registered as a change in their impedance.

Discussion
As a proof of concept, we have developed a levitation and control system to explore properties of superfluid 3 He at microkelvin temperatures. Importantly, this design will allow uniform motion of a levitating body inside the superfluid for the first time. Equally, the same setup can be used to probe the superfluid in a wide range of frequencies of oscillatory motion.
The design consists of a pair of nested coils at the bottom of the cell, which provide the levitation force. The same coils play a role in the control of the sphere's motion. Two pairs of opposing lateral coils in a perpendicular arrangement complete the control system. As shown in this work, the control system allows the uniform motion of the sphere in the horizontal plane. The work presented includes a detailed study of the stability of the set-up and the motion of the sphere in the horizontal plane at uniform velocity. We include the time dependence of the current in the coils needed to control the spurious oscillations of the sphere.
We believe that this levitronic device will open up a new trend in the low temperature research and the exploration of the mechanical properties of quantum fluids as well as topological defects within superfluids.

Methods
Here in this section we include a detailed summary of the experimental setup which is scattered through the text, together with the pseudocode to calculate the stability region. The goal is to facilitate the construction of the device and reproducibility of the calculations needed for the levitation of the sphere.
The data in Table 1 correspond to the coils arrangement shown in the Fig. 1. www.nature.com/scientificreports/ A practical realisation of the hollow sphere would include a plastic ball with a layer of indium deposited on the surface. As long as the mean density of the resulting metallised ball is the same (1.82 g cm −3 ), the parameters of the setup should hold the levitation and stability conditions discussed in this work.
In the implmentation of the stability analysis we used central differences for computing the second order derivatives. The domain is discretized in a N × N mesh, with N = 900 , the step is r = 4.5 × 10 −5 in the radial direction and z = 3.33 × 10 −5 in the vertical direction. Here we included the pseudocode of our implementation: We have developed a Python version of this pseudocode which will be provided freely on request.