Recent discoveries of superconducting phases in the samples of meteorites suggest the possibility of a natural occurrence of superconducting state in space. Superconductors are known to exhibit interesting behaviours when subjected to external magnetic fields, such as levitation. Similar force may act on a superconducting bit in space. The goal of this paper is to quantify this force and assess its effects. Several scenarios in which a superconducting bit can be produced and interact with a magnetic field in space are suggested. The force acting on a superconductor in different conditions is calculated with numerical simulations. The dependence on a magnetic flux density, its gradient, and the geometry and the properties of the superconductor are found. The empirical formulas are derived and used to calculate a magnetic force. The resultant force is extremely weak in all analysed scenarios. It is found that its strength decreases rapidly with the distance from the source of the magnetic flux. Its effect on trajectory of the superconductor is almost negligible. Some possibilities of increasing its strength and the effects are considered.
CERNpublisher-imprint-nameSpringervolume-issue-count12issue-article-count90issue-toc-levels0issue-pricelist-year2021issue-copyright-holderThe Author(s), under exclusive licence to Società Italiana di Fisica and Springer-Verlag GmbH Germany, part of Springer Natureissue-copyright-year2021article-contains-esmNoarticle-numbering-styleContentOnlyarticle-registration-date-year2021article-registration-date-month6article-registration-date-day26article-toc-levels0toc-levels0volume-typeRegularjournal-productNonStandardArchiveJournalnumbering-styleContentOnlyarticle-grants-typeOpenChoicemetadata-grantOpenAccessabstract-grantOpenAccessbodypdf-grantOpenAccessbodyhtml-grantOpenAccessbibliography-grantOpenAccessesm-grantOpenAccessonline-firstfalsepdf-file-referenceBodyRef/PDF/13360_2021_Article_1703.pdfpdf-typeTypesettarget-typeOnlinePDFissue-online-date-year2021issue-online-date-month8issue-online-date-day1issue-typeRegulararticle-typeOriginalPaperjournal-subject-primaryPhysicsjournal-subject-secondaryApplied and Technical Physicsjournal-subject-secondaryTheoretical, Mathematical and Computational Physicsjournal-subject-secondaryComplex Systemsjournal-subject-secondaryCondensed Matter Physicsjournal-subject-secondaryAtomic, Molecular, Optical and Plasma Physicsjournal-subject-collectionPhysics and Astronomyopen-accesstrue
The original version of this article was revised due to a wrong affiliation.
A correction to this article is available online at https://doi.org/10.1140/epjp/s13360-021-02016-4.
Copyright comment
corrected publication 2021
Introduction
Naturally occurring superconducting phases were recently found by Wampler et al. [1] in the samples from two meteorites—iron meteorite Mundrabilla and ureilite GRA 95205 using magnetic field modulated microwave spectroscopy [2]. Some other hints were observed earlier by Guenon et al. [3]. This discovery prompts to analyse the possible effect of superconductivity on the behaviour of celestial bodies and quantify the strength of a magnetic force acting on possible superconducting bodies in space.
The presence of natural superconductors on Earth is rare, and the search for them yielded mixed results [4]. Superconductivity in natural terrestrial minerals was observed in covellite by di Benedetto et al. [5] cooled down below its Tc\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$T_\mathrm{c}$$\end{document} of 1.63 K, and in miassite, with Tc\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$T_\mathrm{c}$$\end{document} of 5.4 K [6]. Some elements naturally present in native form are superconductors too. Notably, Tc\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$T_\mathrm{c}$$\end{document} of tin is 3.72 K and of lead is 7.19 K [7]. It was the alloy of these elements (along with indium), which was detected by Wampler et al. [1]. Another candidate for a space superconductor are some types of diamonds. They are quite abundant, both in stardust with sizes of approximately 2 nm [8] and as larger grains, up to 100 μ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mu $$\end{document}m [9]. Boron-doped diamonds were found to be superconducting at 4 K by Ekimov et al. [10]. Typical temperatures of grains in space were estimated as 5–15 K by Greenberg [11]. This is in line with values of 8–20 K expected in cold dark clouds by Bergin and Tafalla [12]. These temperatures are higher than Tc\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$T_\mathrm{c}$$\end{document} of most of the mentioned superconducting materials, rendering the chances of them exhibiting superconductivity in natural conditions slim. Tc\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$T_\mathrm{c}$$\end{document} of more complex compounds can be significantly higher. Doped C60\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$C_{60}$$\end{document} was observed to be superconducting at 40K [13], and the presence of similar structures in space was described by Maier and Campbell [14].
The most prominent example of a compound superconductor with relatively simple structure and high Tc\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$T_\mathrm{c}$$\end{document} of 39 K is magnesium diboride, which properties were discovered by Nagamatsu et al. [15]. Nanoparticles of MgB2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$_2$$\end{document} can be formed from elemental powders at temperatures above 750 K [16, 17]. The effect of impurities on the quality of a superconductor is non-negligible and was shown by Glowacki et al. [18] to decrease Tc\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$T_\mathrm{c}$$\end{document} even to 22 K. The possibility of formation of MgB2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$_2$$\end{document} during a meteorite impact and subsequent ejection to space cannot be ruled out [19, 20], especially given the presence of both magnesium and boron on the planetary surfaces [21, 22]. Therefore, the existence of a superconducting state in space, while still hypothetical, is possible.
A magnetic force is exerted on a superconductor when it interacts with the external magnetic fluxes. In terms of the experimental and numerical analysis, the most relevant data come from the investigations of magnetic levitation [23]. Significant differences were observed between field cooling (FC) and zero-field cooling (ZFC) conditions [24]. It was also found that the increase of temperature, leading to the decrease of critical current Jc\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$J_\mathrm{c}$$\end{document}, negatively affects the strength of levitation force, as shown by Bernstein et al. [25]. The results of numerical modelling of levitation systems were published by Anischenko et al. [26] and Qin et al. [27]. While providing a useful point of reference, the research on levitation is mostly performed using small-scale systems with permanent magnets serving as the sources of the magnetic flux. The results obtained with such setups are susceptible to the geometry of the magnet and the superconductor, as shown by Yang et al. [28].
Compared to the magnetic fields used during the mentioned experiments, the magnetic fields in space are of much larger scale and generally more uniform in terms of local gradients [29]. Stars can be sources of considerable magnetic fluxes. Typical magnetic flux density on the surface of the Sun ranges from 1.5×10-3\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$1.5\times 10^{-3}$$\end{document} T to 2.5×10-3\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$2.5\times 10^{-3}$$\end{document} T, in sunspots it reaches 1.5×10-1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$1.5\times 10^{-1}$$\end{document} T [30]. Even stronger fluxes (with densities estimated as 1.5×10-1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$1.5\times 10^{-1}$$\end{document} T) are produced by smaller, rapidly rotating M dwarfs [31]. The strongest known sources of magnetic flux in the Universe are magnetars, with magnetic flux density of more than 1010\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$10^{10}$$\end{document} T [32]. In the context of superconductivity, stars and star-related objects are problematic due to their high surface temperature (in the case of magnetars it is more than 106\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$10^6$$\end{document} K, as estimated by Gao et al. [33]), and high energy output, leading to high equilibrium temperatures in their vicinity, dropping below Tc\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$T_\mathrm{c}$$\end{document} only at large distances.
Planetary magnetic fields range greatly in size and strength. In most cases, the source of their magnetic flux is convection in liquid or metallic cores, an effect known as a planetary dynamo [34]. Gas giants can have large magnetospheres and generate relatively strong magnetic fluxes [35]. Equatorial magnetic flux density of Jupiter is reported as 4.17×10-4\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$4.17\times 10^{-4}$$\end{document} T, with maximum values on the surface reaching 2×10-3\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$2\times 10^{-3}$$\end{document} T [36]. Rocky planets, if a liquid core is present and sufficiently active, also can have magnetospheres. Magnetic flux density on the surface of the Earth ranges from approximately 2.3×10-5\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$2.3\times 10^{-5}$$\end{document} T to 6.6×10-5\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$6.6\times 10^{-5}$$\end{document} T and is a subject to constant change [37]. Thus, it is stronger than Saturn’s. Equilibrium temperatures close to the Earth, however, are significantly higher than Tc\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$T_\mathrm{c}$$\end{document} of any known superconductor.
An interesting possibility of interaction between a planetary dynamo and a superconducting bit arises when considering unbound and wide-orbit planets. Stevenson [38] gave an example of a planet ejected from a planetary system and maintaining bodies of water on its surface. Due to the thick atmosphere, such a hypothetical planet would have an effective temperature of approximately 30 K, giving equilibrium temperature around it well below Tc\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$T_\mathrm{c}$$\end{document} of MgB2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$_2$$\end{document}. Meanwhile, high temperature on the surface would be kept thanks to dense and thick atmosphere, nuclear reactions in its interior and residual heat from the time of its formation. The latter implies the presence of the convective currents, which could lead to the generation of a magnetic flux. Another possibility of heating the interior are tidal forces possibly acting on moons of unbound planets [39]. Brown dwarfs also can act as unbound sources of strong magnetic fluxes (with densities up 4.1×10-1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$4.1\times 10^{-1}$$\end{document} T) with relatively low surface temperatures estimated as slightly above 800 K [40].
The goal of this paper is to analyse the behaviour of a small superconducting bit in a magnetic flux that could be generated by a larger celestial body. In order to achieve this goal, a number of numerical simulations are performed to develop simple empirical formulas for a magnetic force acting on a superconductor interacting with a magnetic flux with a given density and its gradient. The formula should also take into consideration the properties of a superconducting material, such as critical current density Jc\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$J_\mathrm{c}$$\end{document} and Tc\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$T_\mathrm{c}$$\end{document}, its geometry. For continuity, the results from numerical calculations used to find the formulas for magnetic force are presented along its formulation in methods. Then, the magnetic force acting on a superconducting bit is calculated in different scenarios. The analysed scenarios include orbiting around a Sun-like star, a magnetar, a brown dwarf, and a rogue planet exhibiting an Earth-like magnetosphere.
Methods
To calculate a magnetic force acting on a superconductor in different conditions, a numerical model is developed in COMSOL Multiphysics. Magnetic field formulation package is employed, which uses H-formulation for calculation of the magnetic fields and current densities, as described further. Several assumptions are made to simplify the numerical model and ease the extraction of dependencies on several factors. The calculations are performed in 2D with circular and rectangular geometries presented schematically in Fig. 1. In the actual model, the environment region is much larger than a superconducting region. It is modelled as a vacuum with relative magnetic permeability μr\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mu _\mathrm{r}$$\end{document} equal to 1. Because of that, the magnetic induction B on the edges of the modelled region can be found as μ0H\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mu _0 H$$\end{document}, where μ0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mu _0$$\end{document} is magnetic permeability of vacuum and H is magnetic field.
General geometry of the numerical model for finding the formulas for magnetic force acting on a superconductor, and the applied boundary conditions
Prescribed value of B is applied as boundary condition on the edges. It is assumed that it has only one nonzero component By\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$B_y$$\end{document}. Analyses are performed for different combinations of background magnetic induction Bb\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$B_b$$\end{document} and its gradient along axis x, ∇By(x)\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nabla B_{y(x)}$$\end{document} The example value of By\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$B_{y}$$\end{document} versus time is presented in Fig. 2a. Initially By\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$B_y$$\end{document} is 0 on both edges. During the first 0.25 s, it is increased to the value Bb\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$B_b$$\end{document} used at a given modelling run, resulting in uniform applied By\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$B_y$$\end{document} across all boundaries. It is maintained as such for 0.75 s to allow for relaxation. Then, ∇By(x)\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nabla B_{y(x)}$$\end{document} is gradually increased. After each increase, the applied By\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$B_y$$\end{document} is constant for the period of 1 s and the total value of magnetic force is read at the end of this time. ∇By(x)\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nabla B_{y(x)}$$\end{document} is changed to have positive and negative values as shown in Fig. 2b to account for the possible hysteretic effects.μ0∂H∂t+∇×EJ=0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} \mu _0\frac{\partial \mathbf{H} }{\partial t}+\nabla \times \mathbf{E} \left( \mathbf{J} \right) =0 \end{aligned}$$\end{document}
Boundary conditions in the numerical model versus time. a Values of By\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$B_{y}$$\end{document} in T prescribed on the edges of a 5-m-wide model region versus time, with the background magnetic induction Bb\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$B_b$$\end{document} of 0.1 T and maximum gradient along axis xBy(x)\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$B_{y(x)}$$\end{document} of 0.01 T. Bl\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$B_\mathrm{l}$$\end{document} is the value prescribed on the left edge, Br\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$B_\mathrm{r}$$\end{document} is the value prescribed on the right edge. b∇By(x)\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nabla B_{y(x)}$$\end{document} in T/m at the same conditions
H-formulation is based on the work by Pecher et al. [41], which in turn is based on Maxwell’s formulas. H-formulation is widely used in modelling of superconductors [42], and it has been successfully applied to model the levitation behaviour of a superconductor by Sass et al. [43]. While applying the formulation, it is assumed that the excitation frequencies are low, allowing to use the quasi-static electromagnetic model. With these assumptions, Faraday’s equation takes form of formula 1. The magnetic field is denoted as H\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{H} $$\end{document}, the electric field is E\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{E} $$\end{document}, t is time, μ0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mu _0$$\end{document} is magnetic permeability of vacuum and J\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{J} $$\end{document} is electric current density.J=∇×H\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} \mathbf{J} =\nabla \times \mathbf{H} \end{aligned}$$\end{document}To satisfy Ampere’s law, J\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{J} $$\end{document} is expressed with formula 2, also being a quasi-static approximation.E=Jσ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} \mathbf{E }=\frac{\mathbf{J }}{{\sigma }} \end{aligned}$$\end{document}Two different formulas are used to find E\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{E} $$\end{document} in different regions of a model, as shown in Fig. 1. In a non-superconducting environment region, representing vacuum surrounding a superconductor, E\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{E} $$\end{document} is found with formula 3. σ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\sigma $$\end{document} here is electrical conductivity. In H-formulation, slightly unphysically, a very small value (less than 1000 S/m) of σ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\sigma $$\end{document} has to be attributed to the environment region to maintain a numerical stability.E=E0|J|-JcJcnJ|J|when|J|≥Jc0when|J|<Jc\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} \mathbf{E }=\left\{ \begin{array}{l@{\quad }ll} E_0\left( \frac{|\mathbf{J }|-J_{c}}{J_{c}}\right) ^{n}\frac{\mathbf{J }}{|J|} &{} \text {when} &{} |\mathbf{J }|\ge J_\mathrm{c}\\ 0 &{} \text {when} &{} |\mathbf{J }|<J_\mathrm{c} \end{array}\right. \end{aligned}$$\end{document}Current–voltage characteristics of a superconductor is highly nonlinear. To account for that in a superconducting region, E\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{E} $$\end{document} is calculated depending on a current density J\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{J} $$\end{document} with the power law (formula 4). In this formula, Jc\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$J_\mathrm{c}$$\end{document} is a critical current density. During the calculations of the magnetic force, it is assumed as constant for a given case, to ease the extraction of the effect of its value on the force. In further calculations, the dependence on the magnetic induction and temperature is introduced. Similarly, the exponent of power law n is assumed to have a constant value of 40. Threshold electric field E0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$E_0$$\end{document} is assumed as 100 μVm-1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathrm {\mu V\, m^{-1}}$$\end{document} [44]. To account for the finite size of the superconducting bit, a constraint setting the total value of the electric current in the superconducting region to 0 is imposed. The numerical mesh used in the model is very dense in the region close to the surface of superconductor, where the shielding electric currents are expected to appear.fL=J×B\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} \mathbf{f} _\mathrm{L}=\mathbf{J} \times \mathbf{B} \end{aligned}$$\end{document}Lorentz force density fL\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{f} _\mathrm{L}$$\end{document} is calculated using formula 5. Since a 2D model is used, only a total force per unit length Fm/L\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{F} _m/L$$\end{document} can be found. This value is the integral of fL\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{f} _\mathrm{L}$$\end{document} over the surface of superconducting region. In order to find the total value of the magnetic force Fm\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$F_m$$\end{document}, Fm/L\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{F} _m/L$$\end{document} has to be multiplied by L, taking into consideration the shape changes along z axis, perpendicular to the xy plane shown in Fig. 1.
The origin of net Fm\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{F} _m$$\end{document} acting on a superconductor is the inhomogeneity of the density of electric current, which arises from the inhomogeneity of Bb\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{B} _b$$\end{document}. Thus, the net Fm\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{F} _m$$\end{document} appears only in a non-uniform background magnetic flux. The strength of Fm\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{F} _m$$\end{document} is affected by several factors. The ones considered in this paper are strength and degree of non-uniformity of Bb\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{B} _b$$\end{document}, material properties of the superconductor represented by Jc\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$J_\mathrm{c}$$\end{document}, and the geometry of the superconducting bit. To quantify their effect, the numerical computations are performed in the wide range of possible values and the resulting Fm\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{F} _m$$\end{document} is calculated. The base values in the computations are: background magnetic flux density of 1 mT, diameter of superconductor (for circular geometry) or side length (for rectangular geometry) of 0.02 m and critical current density Jc\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$J_\mathrm{c}$$\end{document} of 2×108\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$2\times 10^8$$\end{document} A/m2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$^2$$\end{document}.
Dependence of Fm\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{F} _m$$\end{document} on Bb\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{B} _b$$\end{document} and its gradient is presented in Fig. 3. For a given Bb\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{B} _b$$\end{document}, the strength of Fm\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{F} _m$$\end{document} rises linearly with its gradient, as shown in Fig. 3a. No significant dependence of strength and direction of Fm\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{F} _m$$\end{document} on an angle between Bb\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{B} _b$$\end{document} and ∇Bb\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nabla \mathbf{B} _b$$\end{document} is observed. It is important in the considerations of the effect of geometry—different formulas apply for the dimensions in different orientations, as shown further. The net Fm\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{F} _m$$\end{document} always acts towards decreasing Bb\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{B} _b$$\end{document}. In further analysis, the values of Fm\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{F} _m$$\end{document} obtained for a given parameter are collectively considered in terms of a linear coefficient a to help in the representation of dependencies. a of FmL∇B\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\frac{\mathbf{F }_m}{L}\left( \nabla B\right) $$\end{document} increases linearly with Bb\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{B} _b$$\end{document} in a wide range of its strength.
Dependence of |Fm|/L\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$|\mathbf{F} _m|/L$$\end{document} on Bb\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{B} _b$$\end{document} and ∇Bb\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nabla \mathbf{B} _b$$\end{document}. a|Fm|/L\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$|\mathbf{F} _m|/L$$\end{document} as the function of ∇B\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nabla \mathbf{B} $$\end{document} for different Bb\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{B} _b$$\end{document}, b dependence of directional coefficient a—the slope of |Fm|/L\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$|\mathbf{F} _m|/L$$\end{document} on Bb\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{B} _b$$\end{document}. The dependence appears to become nonlinear for very low and high values of Bb\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{B} _b$$\end{document}
Some deviations from the linear dependence of a on Bb\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$B_b$$\end{document} can be seen in Fig. 3b for very low and high values of Bb\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$B_b$$\end{document}. The source of the deviation in the lower range of Bb\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{B} _b$$\end{document} is the reversal of direction of By\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$B_y$$\end{document} occurring in the superconductor region even a when a very small gradient is applied. Additionally, in such conditions the shielding currents flow in the layer much thinner than the resolution of a numerical mesh. Thus, it can be safely assumed that the validity of linear fit extends into the low Bb\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$B_b$$\end{document} range. In the case of high Bb\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$B_b$$\end{document}, a magnetic flux can penetrate the superconductor beyond the surface layer, leading to more complex patterns of the shielding currents. The values of Bb\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$B_b$$\end{document}, in which the deviation is significant, are higher than the ones realistically encountered in space by a superconductor. Thus, the linear fit is sufficient for the purpose of the analysis carried out in this paper. Hysteretic effects are negligible in every analysed case.|a|=b1·tanhb2·Jc+b3b4\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} |a| = b_1\cdot \left( \tanh \left( b_2\cdot J_\mathrm{c} + b_3\right) \right) ^{b_4} \end{aligned}$$\end{document}Values of |Fm|/L\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$|\mathbf{F} _m|/L$$\end{document} versus critical current density Jc\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$J_\mathrm{c}$$\end{document} are presented in Fig. 4. The effect of changing Jc\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$J_\mathrm{c}$$\end{document} even by several orders of magnitude appears to be rather weak in Fig. 4a, but the analysis of the directional coefficients reveals a clear dependence shown in Fig. 4b. A good fit is obtained with formula 6.
The values of fitting parameters for b1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$b_1$$\end{document} to b4\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$b_4$$\end{document}, respectively, are: 0.248, 7.95×10-9\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$7.95\times 10^{-9}$$\end{document}, -0.075 and 1.48×10-2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$1.48\times 10^{-2}$$\end{document} for Jc\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$J_\mathrm{c}$$\end{document} given in A/m2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$^2$$\end{document}. The observed dependence comes from lower local J\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{J} $$\end{document} and the deeper penetration of a superconductor by the magnetic flux, resulting in the lower local values of fL\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$f_\mathrm{L}$$\end{document} acting on superconductor and, ultimately, a smaller net Fm\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{F} _m$$\end{document}. It can be expected that with Jc\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$J_\mathrm{c}$$\end{document} approaching 0, Fm\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{F} _m$$\end{document} will disappear whatsoever, as there will be no shielding currents.
Dependence of |Fm|/L\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$|\mathbf{F} _m|/L$$\end{document} on critical current density Jc\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$J_\mathrm{c}$$\end{document}. a|Fm|/L\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$|\mathbf{F} _m|/L$$\end{document} as the function of ∇Bb\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nabla \mathbf{B} _b$$\end{document} for different Jc\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$J_\mathrm{c}$$\end{document}, b directional coefficient a of |Fm|/L\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$|\mathbf{F} _m|/L$$\end{document} versus ∇Bb\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nabla \mathbf{B} _b$$\end{document} with a hyperbolic tangent fit
The effect of the superconductor geometry on Fm\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{F} _m$$\end{document} is analysed for two geometries—circular and rectangular. Dependence of Fm\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{F} _m$$\end{document} on radius r of a circular superconductor is presented in Fig. 5. Figure 5a shows the results of numerical calculations for different values of r. The value of directional coefficient increases with the square of r, as can be seen in Fig. 5b. In a magnetic field with a constant gradient, the total difference of magnetic flux density across the superconductor rises linearly with increasing r. Similarly, the circumference of the ring in which the shielding currents are flowing also increases linearly, as it is proportional to radius of a circle. The combined effect of these factors gives rise to the observed quadratic dependence.
Dependence of |Fm|/L\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$|\mathbf{F} _m|/L$$\end{document} on r. a|Fm|/L\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$|\mathbf{F} _m|/L$$\end{document} as the function of ∇Bb\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nabla \mathbf{B} _b$$\end{document} for different r, b directional coefficient a of |Fm|/L\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$|\mathbf{F} _m|/L$$\end{document} versus ∇Bb\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nabla \mathbf{B} _b$$\end{document} with a quadratic fit
In the case of the rectangular geometry, the increase of the dimension of the superconductor parallel to ∇Bb\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nabla \mathbf{B} _b$$\end{document} has a very similar effect, as shown in Fig. 6. In the presented analysis, this dimension is the width of the sample and is represented as ratio kw=w/h\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$k_w = w/h$$\end{document}. As in previous case, the increasing difference of Bb\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{B} _b$$\end{document} across the superconductor, and of the surface where the non-symmetric currents appear, combine to form a quadratic dependence (Fig. 6b). The effect of changing the height of the sample (the dimension perpendicular to ∇Bb\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nabla \mathbf{B} _b$$\end{document}) can be described by a logarithmic function defined with formula 7. The results of fitting and calculations for changing height are shown in Fig. 7. The fitted parameters c1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$c_1$$\end{document}, c2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$c_2$$\end{document} and c3\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$c_3$$\end{document} are 0.106, 21.57 and 11.78, respectively.|a|=c1·lnc2kw+c3\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned}&\displaystyle |a| = c_1\cdot \ln \left( \frac{c_2}{k_w} + c_3\right) \end{aligned}$$\end{document}FSL=Cc·|B|·∇B·r2·tanhb2·Jc+b3b4\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned}&\displaystyle \frac{\mathbf{F }_S}{L} = C_\mathrm{c}\cdot |\mathbf{B }|\cdot \nabla \mathbf{B }\cdot r^2\cdot \left( \tanh \left( b_2\cdot J_\mathrm{c} + b_3\right) \right) ^{b_4} \end{aligned}$$\end{document}FSL=Cr·|B|·∇B·w2·tanhb2·Jc+b3b4·lnc2·hw+c3\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned}&\displaystyle \frac{\mathbf{F }_S}{L} = C_{r}\cdot |\mathbf{B }|\cdot \nabla \mathbf{B }\cdot w^2 \cdot \left( \tanh \left( b_2\cdot J_\mathrm{c} + b_3\right) \right) ^{b_4}\cdot \ln \left( \frac{c_2\cdot h}{w}+c_3\right) \end{aligned}$$\end{document}Based on the above considerations, formula 8 is used to find |Fm|/L\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$|\mathbf{F} _m|/L$$\end{document} acting on a circular cross-section of a superconductor. Analogically, formula 9 is used for rectangular cross-sections. The value of constant Cc\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$C_\mathrm{c}$$\end{document} is approximately -2.49×106\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${-2.49\times 10^6}$$\end{document} and Cr\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$C_\mathrm{r}$$\end{document} is -2.67×105\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${-2.67\times 10^5}$$\end{document}, when all parameters are given in SI units.FS=23Cc·|B|·∇B·r3·tanhb2·Jc+b3b4\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned}&\displaystyle \mathbf{F} _S= \frac{2}{3}C_\mathrm{c}\cdot |\mathbf{B} |\cdot \nabla \mathbf{B} \cdot r^3\cdot \left( \tanh \left( b_2\cdot J_\mathrm{c} + b_3\right) \right) ^{b_4} \end{aligned}$$\end{document}FS=Cr·|B|·∇B·w3·tanhb2·Jc+b3b4·lnc2·hw+c3·lnc2·Lw+c3\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned}&\displaystyle \mathbf{F} _S = C_{r}\cdot |\mathbf{B} |\cdot \nabla \mathbf{B} \cdot w^3 \cdot \left( \tanh \left( b_2\cdot J_\mathrm{c} + b_3\right) \right) ^{b_4}\cdot \ln \left( \frac{c_2\cdot h}{w}+c_3\right) \cdot \ln \left( \frac{c_2\cdot L}{w}+c_3\right) \nonumber \\ \end{aligned}$$\end{document}As mentioned, to obtain the formulas for force acting on 3D objects, formulas 8 and 9 have to be integrated along z axis. In the case of a rectangular prism, the logarithmic dependence has to be considered as well. The resultant formulas are 10 for a sphere and 11 for a rectangular prism.
Dependence |Fm|/L\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$|\mathbf{F} _m|/L$$\end{document} on ratio kw\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$k_w$$\end{document} between width and height w/h, with the base height of 0.2 m. a|Fm|/L\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$|\mathbf{F} _m|/L$$\end{document} as the function of ∇Bb\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nabla \mathbf{B} _b$$\end{document} for different kw\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$k_w$$\end{document}, b directional coefficient a of |Fm|/L\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$|\mathbf{F} _m|/L$$\end{document} versus ∇Bb\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nabla \mathbf{B} _b$$\end{document} with a quadratic fit
Dependence of |Fm|/L\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$|\mathbf{F} _m|/L$$\end{document} acting on a rectangular superconductor on 1/kw\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$1/k_w$$\end{document}, with the base width of 0.2 m. a|Fm|/L\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$|\mathbf{F} _m|/L$$\end{document} as the function of ∇Bb\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nabla \mathbf{B} _b$$\end{document} for different 1/kw\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$1/k_w$$\end{document}, b directional coefficient a of Fx/L\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$F_x/L$$\end{document} versus ∇Bb\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nabla \mathbf{B} _b$$\end{document} with a logarithmic fit
The presented formulas are valid for |Bb|\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$|\mathbf{B} _b|$$\end{document} lower than approximately 1 T. They do not consider the effects of field cooling of a superconductor, and the fitted values of parameters are only approximations. Nevertheless, the formulas are sufficiently accurate to provide the realistic values of forces acting on a superconductor in space and to allow the prediction of its behaviour.
In the view of the reasoning presented in Introduction, the most realistic candidate for a superconductor naturally present in space is MgB2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$_2$$\end{document} and its material properties will be used in further calculations. Since it is characterised by low anisotropy [45], Jc\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$J_\mathrm{c}$$\end{document} is assumed to depend only on magnetic flux density B and temperature T. The values are interpolated based on experimental results gathered by Buzea and Yamashita [46] and shown in Fig. 8. The analysed superconductor is assumed to be a spherical piece of pure MgB2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$_2$$\end{document} with the radius of 0.02 m.
Values of critical current density Jc\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$J_\mathrm{c}$$\end{document} of MgB2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$_2$$\end{document}, depending on temperature and magnetic flux density, used in trajectory model (based on data from Buzea and Yamashita [46])
Teq=TbodyR2s1-AB1/4\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} T_{eq} = T_{body}\sqrt{\frac{R}{2s}}\left( 1 - A_B\right) ^{1/4} \end{aligned}$$\end{document}A temperature of a superconducting bit is assumed to be equal to equilibrium temperature Teq\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$T_{eq}$$\end{document} at a given location, found with formula 12. In this formula, Tbody\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$T_{body}$$\end{document} is an average surface temperature of a celestial body interacting with a superconductor, which is treated as a black body. R is the celestial body’s radius, s is the distance between the interacting objects and AB\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$A_B$$\end{document} is Bond albedo of a superconductor, assumed as 0.25. The values of Tbody\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$T_{body}$$\end{document} are listed in Table 1.
Calculation of B\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{B} $$\end{document} is based on a simple numerical model, assuming that a celestial body is a spherical region with an electric current flowing around its axis in angular direction. In other words, the body is treated as a spherical electromagnet. Current density is assumed to be uniform across the volume of the body and its value is selected so the maximum |B|\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$|\mathbf{B} |$$\end{document} on the surface matches the assumed one, according to Table 1. Further calculations are performed only on the surface perpendicular to the axis of the current, where only a single component Bz\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$B_z$$\end{document} of B\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{B} $$\end{document} exists. The obtained values of magnetic flux density are then fitted to find the exact value at a given distance from a body. Gradient is calculated numerically based on the interpolated data.
Table 1 lists several feasible scenarios in which a superconducting bit can interact with bodies generating magnetic flux. Scenario A assumes that a superconductor is interacting with a Sun-like star with the magnetic flux on its surface chosen as the average of a typical range. In scenario B, an interaction with a magnetar is considered, including extremely high magnetic flux densities and temperatures. In scenario C, a relatively low-temperature brown dwarf is selected as a magnetic flux source. Finally, in scenario D a rogue planet with surface magnetic flux density similar to Earth is considered, the major difference being significantly lower effective surface temperature than Earth’s.
FG=GmB·mSd2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} \mathbf{F} _G = G\frac{m_B\cdot m_S}{d^2} \end{aligned}$$\end{document}In each scenario, the values of Fm\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{F} _m$$\end{document} obtained with formula 10 are compared with gravitational force FG\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{F} _G$$\end{document}. Calculation of FG\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{F} _G$$\end{document} is performed using a standard formula (13). G is gravitational constant, d is distance between centre of mass of an orbited body and a superconductor. The ratio Fm/FG\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\left| \mathbf{F} _m\right| /\left| \mathbf{F} _G\right| $$\end{document} is calculated at each respective point. At all considered scenarios, the two forces are acting in opposite directions—FG\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{F} _G$$\end{document} pulls a superconducting bit towards a larger body, while Fm\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{F} _m$$\end{document} repels it.
Results and discussion
The calculated values of Fm\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{F} _m$$\end{document} acting on a superconductor are very small at each scenario. The strongest Fm\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{F} _m$$\end{document} is found with the assumptions from scenario D (rogue planet). The values obtained with this scenario are presented in the following figures. Equilibrium temperature and critical current density are shown in Fig. 11. In all scenarios, the strength of B\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{B} $$\end{document} has a negligible effect on critical current density, and its value is almost completely determined by temperature.
a Equilibrium temperature around a rogue planet, b resultant critical current density
Critical temperature of bulk MgB2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$_2$$\end{document} at the conditions of almost zero magnetic field is 39 K. However, according to the data from Buzea and Yamashita [46], the temperature at which Jc\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$J_\mathrm{c}$$\end{document} of MgB2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$_2$$\end{document} reaches the values higher than 5×106\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$5\times 10^6$$\end{document} A/m2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$^2$$\end{document} is 31.05 K. The considerations shown in Fig. 4 indicate that when Jc\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$J_\mathrm{c}$$\end{document} is below this value, there is almost no magnetic force present. Assuming a uniform temperature distribution of a rogue planet surface, a temperature of 31.05 K is achieved at the distance of approximately 7340 km from its centre, approximately 340 km from its surface (Fig. 9a). Critical current density increases with the distance and decreasing equilibrium temperature, as shown in Fig. 9b.
a Magnetic induction, b gradient of magnetic induction
B\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{B} $$\end{document} and ∇B\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nabla \mathbf{B} $$\end{document} quickly decrease with the distance from the planet centre—magnetic induction decreases as 1/d3\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$1/d^3$$\end{document} and gradient as 1/d4\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$1/d^4$$\end{document}. B\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{B} $$\end{document} is very weak and uniform—maximum strength in the considered region is 5×10-5\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$5\times 10^{-5}$$\end{document} T (Fig. 10a) and maximum ∇B\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nabla \mathbf{B} $$\end{document} is 1.05×10-11\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$1.05\times 10^{-11}$$\end{document} T/m (Fig. 10b). It leads to extremely small values of Fm\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{F} _m$$\end{document}.
a Strength of force due to magnetic interaction between superconductor and a magnetic field of a rogue planet, b ratio between the strength of magnetic force and gravity
As shown in Fig. 11a, despite the relatively low critical current density, Fm\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{F} _m$$\end{document} is strongest in the region closest to the planet, almost immediately after the equilibrium temperature is below the critical temperature. The maximum found value of Fm\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{F} _m$$\end{document} is only 4×10-15\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$4\times 10^{-15}$$\end{document} N for the considered sphere of MgB2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$_2$$\end{document} with r of 0.02 m. With the used assumptions, both the gravitational force and the magnetic interaction force depend on r as r3\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$r^3$$\end{document}. Thus, the ratio between the two will remain the same regardless of the size of a superconducting bit. Assuming that the superconductor is initially on a circular orbit in the region where Fm\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{F} _m$$\end{document} is maximum, the operation of the force will increase the circular orbit radius by 50 nm, making it almost negligible. In the considered geometry, an acceleration due to Fm\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{F} _m$$\end{document} is dependent only on the position of a superconductor and not its mass, giving the maximum value of 4.62×10-14\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$4.62\times 10^{-14}$$\end{document} m/s2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$^2$$\end{document}
The results of calculations for all scenarios are gathered in Table 2. The obtained values of the magnetic interaction force are even lower for other bodies. With the force decreasing as 1/d7\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$1/d^7$$\end{document}, the larger distances of superconducting transition present in other scenarios mean that even the extremely strong sources of magnetic flux such as magnetars fail to exert a meaningful force on a superconductor. On the other hand, such dependence opens some possibilities to increase the effect of superconductivity on a superconductor trajectory.
One of the assumptions taken during the calculations of Fm\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{F} _m$$\end{document} is that the temperature of a superconductor is equal to an equilibrium temperature at a given location. It is not strictly true, as each body possess a certain thermal inertia. Therefore, it is possible for a sufficiently large and fast-moving superconductor passing by an object to maintain a superconducting state in the region with higher magnetic flux density and thus to experience a stronger force. In the considered example of a rogue planet, the closest possible approach to d equal to an assumed atmosphere radius leads to a twofold increase of Fm\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{F} _m$$\end{document} when compared with the case not considering heat capacity of a superconductor.
The assumed uniform, spherical form of a superconductor is rather unlikely. It can be expected that a superconducting body, similarly to the one described at Wampler et al. [1], would contain multiple smaller, disjointed domains with complex shapes. The elongation of the superconductor interface in the direction of a magnetic field gradient can lead to a larger effective Fm\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{F} _m$$\end{document}. Additionally, the presence of multiple domains can lead to more complex interactions, such as the appearance of a time-varying angular momentum, and cause the superconductor to rotate in resonance with the changes in magnetic flux. Such effects can be further strengthened by non-uniformities of the external magnetic flux density.
The effects of superconductivity can become non-negligible at larger time scales when considering a superconductor orbiting a larger body. Since the force is usually directed away from the source of magnetic flux, it can be a mechanism counteracting drag and radiation pressure. Additionally, a superconducting bit may significantly change a magnetic environment in its vicinity. It could be especially important in affecting the behaviour of space dust and the formation of larger bodies, if they contained ferromagnetic particles.
Generally, though, the possible presence of natural superconductivity at small bodies would have a very small effect on their trajectories. The generation of the natural magnetic fluxes is often connected with the movement of conducting masses and an associated emission of heat. It leads to the increase of equilibrium temperature above the critical temperature in the regions where the magnetic fluxes and their gradients are high enough to generate an appreciable force. This does not exclude the possibility of using artificial cooling and artificial superconductors to allow for a larger force generation. In the case of the terrestrial magnetic field, the effect would still be too weak to be practical. On the other hand, Jupiter’s strong magnetosphere and lower temperatures around the planet may make the use of superconducting magnetic sails possible.
Conclusion
The natural existence of a superconducting state in space is possible, and it could have some effect on the behaviour of a small body. However, naturally occurring superconductivity fails to produce a strong magnetic force in viable scenarios. In most cases, a magnetic flux generation is connected with a heat emission. It causes an equilibrium temperature in the vicinity of the magnetic flux source to be above a critical temperature of the superconducting materials, which may be produced and ejected to space by natural processes.
A magnetic flux density and its gradient decrease rapidly with the increasing distance from a body. The derived formulas show the dependence of magnetic force on the distance between the bodies as 1/d7\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$1/d^7$$\end{document}. That makes the force at the regions with sufficiently low equilibrium temperature extremely small, to the point where it can be neglected. Any significant effect of superconductivity on the trajectory of a body would require long time scales or a cold source of strong magnetic flux.
The obtained results do not rule out the possibility of using an artificial cooling and high-temperature superconductors to obtain higher forces at specific conditions. It is planned to explore this possibility in the future works on the subject. Superconducting magnetic sails can find use as the method to maintain an orbit by counteracting drag forces. Superconductivity can also play some role in the formation of larger bodies and behaviour of space dust by generating small regions with relatively strong gradients of magnetic field.
Funding
Open Access funding provided by CERN.
ReferencesJ. Wampler, M. Thiemens, S. Cheng, Y. Zhu, I.K. Schuller, Superconductivity found in meteorites, in Proceedings of the National Academy of Sciences, 2020. ISSN 0027-8424. https://doi.org/10.1073/pnas.1918056117. https://www.pnas.org/content/early/2020/03/17/1918056117RamírezJGBasaranACde la VentaJPereiroJSchullerIKMagnetic field modulated microwave spectroscopy across phase transitions and the search for new superconductors20147790939022014RPPh...77i3902R10.1088/0034-4885/77/9/093902GuenonSRamirezJGBasaranACWamplerJThiemensMSchullerIKSearch for new superconductors: an electro-magnetic phase transition in an iron meteorite inclusion at 117 K2017297–30424010.1007/s10948-016-3708-7R. Wang. Search for natural occurring superconductivity and novel phenomena: magnetic transitions in natural transition metal compounds. PhD thesis, University of Maryland (2017)F. di Benedetto, M. Borgheresi, A. Caneschi, G. Chastanet, C. Cipriani, D. Gatteschi, G. Pratesi, M. Romanelli, R. Sessoli. First evidence of natural superconductivity: covellite. European Journal of Mineralogy, 18(3), 283–287 (2006). ISSN 0935-1221. https://doi.org/10.1127/0935-1221/2006/0018-0283R. Settai, K. Katayama, H. Muranaka, T. Takeuchi, A. Thamizhavel, I. Sheikin, Y. Onuki, Superconducting properties in Rh17S15 under magnetic field and pressure. Journal of Physics and Chemistry of Solids, 71(4): 700–703, (2010). ISSN 0022-3697. https://doi.org/10.1016/j.jpcs.2009.12.070. http://www.sciencedirect.com/science/article/pii/S0022369709004132. 15th International Symposium on Intercalation Compounds - ISIC 15 May 10-15, 2009, BeijingEisensteinJSuperconducting elements1954262772911954RvMP...26..277E10.1103/RevModPhys.26.277TielensAGGMThe molecular universe201385102110812013RvMP...85.1021T10.1103/RevModPhys.85.1021NabieiFBadroJDennenwaldtTOveisiECantoniMHébertCEl GoresyABarratJ-AGilletPA large planetary body inferred from diamond inclusions in a ureilite meteorite20189113272018NatCo...9.1327N10.1038/s41467-018-03808-6ISSN 2041-1723E.A. Ekimov, V.A. Sidorov, E.D. Bauer, N.N. Mel’nik, N.J. Curro, J.D. Thompson, S.M. Stishov, Superconductivity in diamond. Nature 428(6982), 542–545 (2004). https://doi.org/10.1038/nature02449 (ISSN 1476-4687)GreenbergJMInterstellar grain temperatures. effects of grain materials and radiation fields1971122401971A&A....12..240GBerginEATafallaMCold dark clouds: The initial conditions for star formation20074513393962007ARA&A..45..339B10.1146/annurev.astro.45.071206.100404TanigakiKPrassidesKConducting and superconducting properties of alkali-metal C60 fullerides199551515152710.1039/JM9950501515MaierJPCampbellEKFullerenes in space201756184920492910.1002/anie.201612117NagamatsuJNakagawaNMuranakaTZenitaniYAkimitsuJSuperconductivity at 39 k in magnesium diboride2001410682463642001Natur.410...63N10.1038/35065039ISSN 1476-4687C. Cui, D. Liu, Y. Shen, J. Sun, F. Meng, R. Wang, S. Liu, A. Greer, S. Chen, B. Glowacki, Nanoparticles of the superconductor MgB2: structural characterization and in situ study of synthesis kinetics. Acta Materialia, 52(20), 5757–5760, (2004). ISSN 1359-6454. https://doi.org/10.1016/j.actamat.2004.08.033. URL http://www.sciencedirect.com/science/article/pii/S1359645404005038MuralidharMKentaNKoblischkaMRMurakamiMHigh critical current densities in bulk MgB2 fabricated using amorphous boron201521210214121452015PSSAR.212.2141M10.1002/pssa.201532108B. Glowacki, M. Majoros, M. Eisterer, S. Toenies, H. Weber, M. Fukutomi, K. Komori, K. Togano, MgB2 superconductors for applications. Physica C: Superconductivity, 387 (1): 153 – 161, (2003). ISSN 0921-4534. https://doi.org/10.1016/S0921-4534(03)00662-2. http://www.sciencedirect.com/science/article/pii/S0921453403006622. Proceedings of the 3rd Polish-US Workshop on Superconductivity and Magnetism of Advanced MaterialsCoxMACavosieAJBlandPAMiljkovićKWingateMTMicrostructural dynamics of central uplifts: Reidite offset by zircon twins at the Woodleigh impact structure2018461198398610.1130/G45127.1ISSN 0091-7613RubinAEShock, post-shock annealing, and post-annealing shock in ureilites20064111251332006M&PS...41..125R10.1111/j.1945-5100.2006.tb00197.xJ.D. Stephenson, L.J. Hallis, K. Nagashima, S.J. Freeland, Boron enrichment in martian clay. PLOS ONE, 8(6): 1–4 (2013). https://doi.org/10.1371/journal.pone.0064624VanimanDTBishDLChiperaSJFialipsCIWilliam CareyJFeldmanWCMagnesium sulphate salts and the history of water on Mars200443170096636652004Natur.431..663V10.1038/nature02973ISSN 1476-4687BernsteinPNoudemJSuperconducting magnetic levitation: principle, materials, physics and models20203330330012020SuScT..33c3001B10.1088/1361-6668/ab63bdMorandiAFabbriMRibaniPLDennisADurrellJShiYCardwellDThe measurement and modeling of the levitation force between single-grain YBCO bulk superconductors and permanent magnets201828511010.1109/TASC.2018.2822721BernsteinPColsonLDupontLNoudemJInvestigation of the levitation force of field-cooled YBCO and MgB2 disks as functions of temperature20173060650072017SuScT..30f5007B10.1088/1361-6668/aa69ecAnischenkoIPokrovskiiSRudnevIOsipovMModeling of magnetization and levitation force of HTS tapes in magnetic fields of complex configurations2019321010500110.1088/1361-6668/ab2bbeQinMJLiGLiuHKDouSXBrandtEHCalculation of the hysteretic force between a superconductor and a magnet2002660245162002PhRvB..66b4516Q10.1103/PhysRevB.66.024516YangWMChaoXXBianXBLiuPFengYZhangPXZhouLThe effect of magnet size on the levitation force and attractive force of single-domain YBCO bulk superconductors20031677897922003SuScT..16..789Y10.1088/0953-2048/16/7/308DurrerRNeronovACosmological magnetic fields: their generation, evolution and observation2013211622013A&ARv..21...62D10.1007/s00159-013-0062-7ISSN 1432-0754OssendrijverMThe solar dynamo20031142873672003A&ARv..11..287O10.1007/s00159-003-0019-3ISSN 1432-0754ShulyakDReinersAEngelnAMaloLYadavRMorinJKochukhovOStrong dipole magnetic fields in fast rotating fully convective stars20171801842017NatAs...1..184S10.1038/s41550-017-0184ISSN 2397-3366DuncanRCThompsonCFormation of very strongly magnetized neutron stars: implications for gamma-ray bursts1992392L91992ApJ...392L...9D10.1086/186413GaoZWangNPengQ-HThe surface and inner temperatures of magnetars201381210.1017/S1743921312024222S. Labrosse, Thermal and magnetic evolution of the Earth’s core. Physics of the Earth and Planetary Interiors, 140(1): 127 – 143, (2003). ISSN 0031-9201. https://doi.org/10.1016/j.pepi.2003.07.006. http://www.sciencedirect.com/science/article/pii/S0031920103001687. Geophysical and Geochemical Evolution of the Deep EarthKruppNGiant magnetospheres in our solar system: Jupiter and Saturn compared2014221752014A&ARv..22...75K10.1007/s00159-014-0075-xISSN 1432-0754J.E.P. Connerney, S. Kotsiaros, R.J. Oliversen, J.R. Espley, J.L. Joergensen, P.S. Joergensen, J.M.G. Merayo, M. Herceg, J. Bloxham, K.M. Moore, S.J. Bolton, S.M. Levin, A new model of Jupiter’s magnetic field from Juno’s first nine orbits. Geophys. Res. Lett. 45(6), 2590–2596 (2018). https://doi.org/10.1002/2018GL077312A. Chulliat, S. Macmillan, P. Alken, C. Beggan, M. Nair, B. Hamilton, A. Woods, V. Ridley, S. Maus, A. Thomson, The US/UK world magnetic model for 2015-2020: Technical report. National Geophysical Data Center, NOAA, (2015). https://doi.org/10.7289/V5TB14V7StevensonDJLife-sustaining planets in interstellar space?1999400673932321999Natur.400...32S10.1038/21811ISSN 1476-4687DebesJHSigurdssonSThe survival rate of ejected terrestrial planets with moons20076682L167L1702007ApJ...668L.167D10.1086/523103KaoMMHallinanGPinedaJSStevensonDBurgasserAThe strongest magnetic fields on the coolest brown dwarfs20182372252018ApJS..237...25K10.3847/1538-4365/aac2d5R. Pecher, M. McCulloch, S. Chapman, L. Prigozhin, 3d-modelling of bulk type-ii superconductors using unconstrained h-formulation, in Proceeding of the EUCAS 2003 (2003)ShenBGrilliFCoombsTOverview of h-formulation: A versatile tool for modeling electromagnetics in high-temperature superconductor applications2020810040310041410.1109/ACCESS.2020.2996177SassFSoteloGGde Andrade JuniorRSiroisFH-formulation for simulating levitation forces acting on HTS bulks and stacks of 2g coated conductors201528121250122015SuScT..28l5012S10.1088/0953-2048/28/12/125012KvitkovicJPatelSZhangMZhangZPeetzJMarneyAPamidiSEnhanced magnetic field sensing using planar high-temperature superconductor shields20182841510.1109/TASC.2018.2813538ISSN 1051-8223VinodKKumarRGASyamaprasadUProspects for MgB2 superconductors for magnet application2006201R1R1310.1088/0953-2048/20/1/r01BuzeaCYamashitaTReview of the superconducting properties of MgB220011411R115R1462001SuScT..14R.115B10.1088/0953-2048/14/11/201