Arxiv:2003.07985V1

Home , Electrodynamic tether

arXiv:2003.07985v1 [astro-ph.EP] 17 Mar 2020 Keywords: sadtoa xlrtr aevrn rvdn oemis more providing — maneuvering to exploratory due additional consideration as under is tethers electrodynamic using e acltossol emd oflyicueNpuinr Neptunian include fully to made be should calculations new Dynamics Composition il aiu fcec fapoiaey1,weesa10 a whereas con this 12, it; approximately facing of S/C efficiency the maximum with yield dipole, the of plane rati meridian mass (spacecraft-to-tether efficiency capture make to oecmlx h ioei oae eo h qaoilpla equatorial the below located is dipole the complex: more xpae statistics Exoplanet Giants: Gas with associated G those Ice from exploring different in decade are interest next that of the points multiple for are missions There flagship [2]. as NASA by considered orbite standard a t sion. than potential capabilities the mission have greater be and create [1] — exploration properties further thermodynamic as their of well pow as and insertion, propulsion orbital free for provide specific can Encel involving Saturn’s tethers or missions electrodynamic Europa minor Jupiter’s exploring For as such visits Giants. Gas knowledg significant the provide to about helped have Jupiter— at Juno Introduction 1. comp When outside. field the describes s which in — currents to approximation due are fields magnetic Planetary densities. field magnetic and density Gas for detailed and broad been have missions planetary Past Abstract rpitsbitdt Elsevier to submitted Preprint [email protected] oet-rgwr erae ucl ihdsac,tu re thus distance, with quickly decreases work Lorentz-drag ∗ mi addresses: Email author Corresponding h w c ins—rnsadNpue aebeen have Neptune— and —Uranus Giants Ice two The and Galileo as well as Saturn, at —Cassini missions Broad isi h citcpaeisl.Npue—niethe —unlike Neptune itself. plane ecliptic nearly the Uranus in of axis lies rotation The So- early dynamics. intriguing, System the lar in bodies big other with collision mut fwtr moi,admethane. and ammonia, substantial water, contain of amounts Neptune and Uranus and hydrogen while of helium— made primarily are Saturn and —Jupiter Giants. Gas than abundant more hw iglrfaue htaepsil in of signs possible are that features singular Shows ue ins c on,Eetoyai ehr,Magnetic tethers, Electrodynamic moons, Icy Giants, Outer akdydfeetbtenIeadGsGiants Gas and Ice between different Markedly [email protected] J Peláez)(J. aasget htIeGat r much are Giants Ice that suggests Data B huhttescnda ihildfie est rfie,w profiles, density ill-defined with deal can tethers though , b nvria oićiad ard z .inrs3 Madr 3, C.Cisneros Pz. Madrid, Politécnica de Universidad ehrCpueo pccata Neptune at spacecraft of Capture Tether J .Sanmart´ın), R. (J. a .R Sanmart´ın R. J. elAaei eIgneáo Spain Ingenier´ıa of de Academia Real )wl bv oofe ae h / ih pial reach optimally might S/C The case. no-offset a above well o) m ml oueisd h lnt antcmmn vector magnetic-moment planet, the inside volume small ome urn pccatprassa atr ls oteplanet the to close capture at periapsis spacecraft quiring h blt fttest rvd repouso n oe f power and propulsion free provide to tethers of ability the ein ycrns sesdb etn oaigltl dur little rotating Neptune by eased is synchronism venient tto n osdrdtie ioeadqarpl correc quadrupole and dipole detailed consider and otation incpblt hnasadr rie iso.Tte ope Tether mission. orbiter standard a than capability sion ins oprdt yymsin o c ins Presently Giants. Ice for missions flyby to compared Giants, ◦ e shgl fstfo h lntcne,ada ag tilt large at and center, planet the from offset highly is ne, mis- r cause adus, eiinerrwudrdc fcec yaot6.Efficienc 6%. about by efficiency reduce would error meridian iants rdwt aunadJptr h etna antcstruct magnetic Neptunian the Jupiter, and Saturn with ared er o e a,b oc e ntlength unit per force elblwamxmm hr-ici value short-circuit maximum: a below well ntec-oaigmgeie min lsaamtoa fiel motional a plasma E ambient magnetized co-rotating the in hntp ehradmsigcleto 5.Ti requires length This tether [5]. ing collection missing cross and electrons tether attracted thin-tape energetic and/or heating tether periefiin scmae ihJptr[] nefc,t effect, In [3]. could Jupiter tethers field with case, magnetic Saturn compared a as in inefficient As appear Neptune. to mission minor a ayaogtette cuigi oata in Langmuir giant high a its Jupiter because as for expected than act requirements less plasma-bias capture-efficiency to relative Design it —causing and probe. tether current the making along [4], vary end conduc spherical anodic large the a and insulation tether eliminating iny(/-otte asratio, mass (S/C-to-tether ciency antcfil structure field Magnetic ekrfils h / eaievelocity relative S/C The fields. weaker ehrcnutvt) u oomcefcsadproportiona and effects field ohmic to due conductivity)— tether u;cpueefiinyi prxmtl . o ohJupit both [6]. for Saturn 3.5 and approximately is capture-efficiency short-circu mum; particular the near reaching current-density .Peláez J. , m ≡ h olwn icse hte ehr ih eue for used be might tethers whether discusses following The irpiebr-ehrcnetehne urn o by flow current enhanced concept bare-tether disruptive A capture B hti ihyipratfrtte neato sde- as — analysis. interaction present tether the in for tailed important highly is that rtn hc si ergaeorbit. retrograde in moon, is large which one Triton, just has moons— many with Giants other v o ekrfils ehr a vi hs sus with issues, those avoid may tethers fields, weaker For . ′ ∧ b, B ∗ nteSCrfrnefae and frame, reference S/C the in t h ndcsgetsl-dutn oaccommodate to self-adjusting segment anodic the ith B d200 Spain 28040, id ssmlrysal n pccatcpueeffi- spacecraft-capture and small, similarly is L t I t eplnt-vrgdcretdensity current length-averaged keep –to ∧ B ntte current tether on etn hw tiigdifference striking a shows Neptune M B SC au ih euti strong in result might value n loigtelreoffset large the allowing and / eipi hncosn the crossing when periapsis n yial ioelaw dipole a typically and , robtlisrina well as insertion orbital or m aindpnso plasma on depends ration n atr.Calculations capture. ing tions. t erae as decreases ) v iso oNeptune to mission a , ′ ihisrtto axis. rotation its with I ≡ rvnby driven B r ssignificantly is ure σ v eut suggest results y t xrsLorentz exerts − E ac 9 2020 19, March m v pl ( σ tmaxi- it E induces t n the ing B m limit- being o at tor 2 . to l are for he er d Section 2 introduces Neptune’s environmentissues concern- LOWER ELECTRONIC ing incomplete data on the ionospheric plasma electron density DENSITY Ne and magnetic field structure. Section 3 considers parameters e− e− of Neptune as a viable planetary tether environment, which is standard in an operation scheme. In Section 4, a simple ver- Electrons sion of the complex dipole model OTD2 —regularly used in collected the literature on Neptunian magnetics— is used to determine along tether capture-efficiency. Results exceeding the Jupiter/Saturn values segment of are discussed in Section 5 and Conclusions are summarized in positive bias Section 6. PSfrag replacements Negligible ion 2. Neptune environment issues impact rate

Tether operation depends on electron density Ne —appearing Direction of in tether current— and magnetic field in both motional field and current flow Lorentz force. The required and available information for both electron density and magnetic field have drastically different

Plasma characteristics. The short 1989 Voyager 2 ﬂyby, which came contactor very close to Neptune, made getting models from data a difﬁ- cult task. e− e−

2.1. Plasma density issues Voyager 2 provided three types of data that yielded conflict- Figure 1: Interaction tether-ambient plasma ing descriptions of plasma density at Neptune, two involving instruments on board, a third one using radio signals from the A bare tether accommodates density variations by allowing the spacecraft to Earth. First, there are in-situ data from measure- anodic segment to self-adjust by varying its fraction of tether ments by the PLS (plasma) instrument, using 4 modulated-grid length (Figure 1). This could make a broad range of Ne values Faraday cup detectors [7]. Next, a Plasma Wave instrument lead to length-averaged current near the short-circuit maximum receiving wave (particularly whistler wave) data from a Plane- wt ht σt Em (a design reference-value independent of actual × tary Radio Astronomy on-board antenna provided electric-field electron densities), for a convenient range of Lt values, not af- wave data into a wideband waveform receiver [8]. Finally, ra- fecting tether mass if tether-tape width wt is adjusted at a fixed dio signals from Voyager 2 to Earth, following S/C occultations thickness ht . by Neptune, allowed tracking stations in Australia and Japan to In the OML (orbital-motion-limited)electron-collection reg- determine ionospheric Ne profiles versus altitude [9]. imen of bare tethers, the Lorentz force on the tether is propor- 3 The Faraday cups showed maximum density Ne 2 cm− tional to its length-averaged current, ∼ during the brief closest approach at 79 North, r 1.18 RN, or ◦ ≈ about 4,500 km altitude [10] — while both radio-antenna and Iav = wt ht σt Em iav(Lt /L⋆) (1) × tracking-stations showed densities orders of magnitude above. 1/3 2/3 4 3 L⋆ ∝ Em (σt ht /Ne) (2) Tracking data ranged from values of order 10 cm− around 2 3 1,250 km, down to 10 cm− in the 2,500/5,000 km altitude where L⋆ is a characteristic length gauging the bare-tether electron- range [11]. The waveform receiver showed intermediate densi- collection impedance against ohmic resistance and increasing 3 3 ties —around 10 cm− — at in-between altitudes [12]. with decreasing Ne [3]. The dimensionless average-current iav < Modeling of whistler propagation did present some issues 1 (see Figure 2), which vanishes with length ratio Lt /L⋆, ap- [12]. More importantly, regarding conflicting data, the ques- proaches 1 at large Lt /L⋆ following the law: tion was raised whether Voyager 2 had entered Neptune’s iono- iav = 1 L⋆/Lt for Lt /L⋆ > 4. (3) sphere — or whether a cold-density plasma component was de- − tectable by the PLS instrument, which might be interacting with Equation (3) shows effective OML-current accommodation to precipitating particles and not co-rotating flow [10]. Further, drops in electron density. Consider L /L decreasing from 20 radio-tracking data above 2,500 km were problematic because t ⋆ to 4, with tether length L and parameters in equation (2) other of variations in the ionosphere and interplanetary medium that t than N to be constant. The value of N itself would then be the radio signal traverses [11]. e e lower by a 0.09 factor, whereas current i in (3) would decrease Certainly, data from the Voyager 2 flyby did not yield a av from 0.95 to 0.75. full model of the ambient magnetized plasma around Neptune. Tether operation, however, would first involve plasma density 2.2. Magnetic field issues in the lower, no-conflict ionosphere, and secondly deal with ill- Regarding magnetic fields, B, the fields of Gas Giants are defined density profiles by appealing appropriately sized tapes. present due to currents from chargesrepeatedly moving in some 2 1.0 3. Neptune particular S/C-capture regime

0.9 For the hyperbolic orbit of a S/C in Hohmann transfer between heliocentric circular orbits at Earth and Neptune, the ar- 0.8 riving velocity in the Neptune frame is about v∞ = 3.96 km/s,

t 2

h 0.7 t the orbital speciﬁc energy εh being v∞/2. Using that value in

m the general relation between ε and eccentricity e at given peri- E 0.6

t σ apsis rp,

I 0.5 ε = µN (e 1)/(2rp) (5) − yields a hyperbolic eccentricity very close to unity, 0.4 2 eh 1 = v RN/µN = 0.058 (6) 0.3 − ∞ 3 2 0.2 where µN 6835107 km /s is the gravitational constant, and we conveniently≈ set the periapsis very close to the planet (as 0.1 in Jupiter and Saturn analyses [5, 6]) rp RN 24,765 km. Lorentz-drag capture results in barely elliptical≈ ≈ orbits (eccen- 0.0 PSfrag replacements 0 1 2 3 4 5 6 7 8 9 10 tricity ec just below 1). For conditions of interest, calculations L/L ∗ will be carried out using a parabolic orbit throughout, with no sensible change except moving from open to closed. The re- Figure 2: Length-averaged current vs. tether length quired capture-dynamics is, in a sense, weak. When the relative velocity v′ opposes v the Lorentz force small volume deep inside the planet. Outside the planet, at is actually thrust. This is the case for prograde circular orbits “large” distances from the system of steady currents but close beyond a radius as where plasma co-rotating with planet spin to the planet, the field is generally described through the uni- and orbital velocities are equal. versal magnetic-moment vector concept and its dipole-law ap- 2 1/3 ΩaS = µ/aS aS/R ∝ ρpl/Ω (7) proximation. A magnetic moment m = mum of magnitude m ⇒ 3 (gauss meter ) and unit vector um located at a “point” rm gives Among Giants, Neptunep and Saturn have the highest and lowest × a definite magnetic field at a “faraway” point r [13] densities, respectively, while both Ice Giants have spin slower than Gas Giants. This results in values a /R 1.89and a /R 3 s S s N B(r) 3(uρ m)uρ m /ρ (4) ≈ ≈ ≈ · − 3.41, with correspondingvalues —comparedto Jupiter and Earth— in between the two, and larger than either, respectively. For with ρ ρ u r r vector position in planetary frame with = ρ m the parabolic orbits of interest, the radius where drag vanishes displaced origin.≡ For− Saturn, m is at its center (r 0, ρ r) m = along with the relative tangential velocity v is [3] and near parallel to its rotation axis — with Jupiter≈ also having t′ relatively similar conditions. At points in a circular equatorial 3/2 2 1/2 rM /R = √2(aS/R) ∝ ρpl/Ω (8) orbit, the field would then read B(r)= m/r3, and the Lorentz force on a tether, orbiting vertical, is parallel− to the velocity. varying from 3.7 for Saturn to 8.9 for Neptune. For Neptune, a dipole law with rm well off-center and um This makes for important capture differences at Neptune, complex orientation is valid for somewhat large r/R [14] dis- with most of the drag arc from rp to rM rp being irrelevant, N ≫ tances that are of no interest for tether applications due to the because Lorentz drag becomes negligiblewell before r reaching Lorentz force, which is quadratic in the field, decreases with the rM, quite opposite the case with Saturn, which made its retro- inverse 6th power of distance to the planet. In the OTD2 offset- grade capture convenient[6]. This is manifest here in the veloc- 1/2 tilted detailed model [15], with magnetic moment of magnitude ity at the parabolic periapsis, vp = (2µN/RN) 23.5 km/s 3 ≈ m = 0.13G R , the dipole is located 0.19 RN below the equa- being much larger than co-rotation velocity 2πRN /16.1 h × N ≈ torial plane and 0.52 RN radially away from the rotation axis 2.7 km/s, and allows using v v throughout the relevant drag ′ ≈ —in the plane containing axis and dipole— where um is tilted arc. 47◦ with respect to the axis and 22◦ off the meridian plane. Furthermore, because of the large dipole offset —equal mo- For r < 2.5 RN , quadrupole terms in a multipole expan- ment m in capture is greatly more efficient than for a no-offset sion of field B decrease faster than the dipole —and/or addi- case— the S/C might optimally reach periapsis when crossing tional multipoles from differently localized current sources— the meridian plane that contains the dipole center, resulting in and may generally need be considered at differencewith Jupiter the S/C facing the dipole when at periapsis (Figure 3). This and Saturn [16, 17, 18]. Alternative spherical harmonic analysis convenient synchronism is somewhat eased by Neptune having has reached a reasonable quadrupole description, though rele- slow spin and high density —with the high density making S/C vance is dependent on longitude and latitude ranges. For sim- orbital motion fast— as following from the Barker equation for plicity in the present analysis we shall keep the above dipole parabolic orbits, giving time t from periapsis-pass versus radius, term with a reduced approximation. (3vp/2rp)t = √r˜ 1 (2 + r˜), (r˜ r/rp) (9) − ≡ 3 v rP

z y ρ −→m r θ f R 2rP 2rP

r ρ θ PSfrag replacements f R x

PSfrag replacements Figure 4: Full synchronism case

At full synchronism, the symmetry relation between ρ and r at positive and negative values of true anomaly θ leads to equal Lorentz-drag work in symmetric arcs at periapsis. We may thus Figure 3: Magnetic dipole and S/C trajectory rP write the Lorentz-drag work as

r u W˙ Ldr The characteristic time of the S/C motion, for rp R, is thus of WL = W˙ L dt = 2 (14) ≈ ∆t × R dr/dt order Z Z N 3 rp/vp R /2GMpl ∝ 1/√ρpl. (10) the upper limit r , though well below r , is here considered ≈ u M Because of the slow planetq spin and the short range of Lorentz —to include the radial arc— and contributes to drag. We then contribution to its drag work, we will here estimate ρ distances have while ignoring the rotation of Neptune and its rm position. ru 2 mt 3 2 v 2 rdr WL = σt 0.13GRN cos47◦ iav r 6 × − ρ h i R ρ t Z N 2µN(r RN) 4. Capture-efficiency calculation (15)− where we used the radial speed-rate in parabolic arcs,p dr/dt = In the present estimate, we ignore both distance 0.19 R 1/2 N [2µN(r RN )] /r, finally yielding and the 22◦ angle. We consider Neptune magnetics as present- − m ing a 0.55 RN offset and a 47◦ tilt, both having a substantial t 2 WL = σt (0.13G cos47◦) iav r vpRN I(r˜u) (16) effect on capture efficiency. Also, as in [3, 5, 6], we shall con- − ρt h i × sider a S/C approachingNeptune in an equatorial orbit. Writing r˜u 2 d˜r I(r˜u) × (17) B = Bax +Beq, where subscripts ax and eq stand for field com- ≡ 1 ρ˜ 6√r˜ 1 ponents along the Neptunian rotational axis and in the equato- Z − rial plane, respectively, we may rewrite a standard power law ρ˜ andr ˜ keeping a simple relation along the S/C orbital motion W˙ L = v (Lt Iav B) as (see Figure 4). · ∧ 2 2 W˙ L = v (Lt Iav Bax)= Lt Iav (v Bax) (11) ρ˜ = (r˜ + f ) 4 f , (ρ˜ ρ/RN, r˜ r/RN, f = 0.55). · ∧ − · ∧ − ≡ ≡ (18) All three vectors v, Iav and Beq lie in the equatorial plane, thus Using Equation (5-6) in drag-work per unit mass having no joint power effect (Figures 3,4). Having ignored the 2 0.19 RN distance parallel to the Neptunian axis makes the sec- WL µN v∞ eh ec | | = εh εc = (eh ec)= − , (19) ond term in Equation (4) contribute to Bax. MSC − 2rp − 2 eh 1 − 4.1. Full synchronism case and Equation (16) then yields capture efficiency Using Equation (1) with m-subscript omitted, we write MSC MSC eh ec WL < − = | 2 | = Lt Iav = Lt wt ht σt iav EequI (12) mt mt × eh 1 mt v∞/2 × − σt vpRN 2 where unit vectors along tether (u ) and E are both radial at = 2 (0.13G 0.68) iav I(r˜u) (20) I eq ρ v2 × h ir × periapsis and nearly so in a short effective-drag arc in Figure 4. t ∞ Also, using v′ v to write EequI Eeq v Bax in Equation At the lower limit in the integral,r ˜ = 1, Equation (18) gives ≈ ≈ ≈ ∧ l (11), Lorentz power becomes ρ˜l = 1 f = 0.45. Now we consider upper values, ρ˜u = n(1 − − f ) ρñ, with n > 1 yielding forr ˜u, at Equation (18), ˙ mt 2 2 mt 2 3 2 ≡ WL(r,ρ)= σt iav v Bax = σt iavv m cos47◦/ρ − ρt × − ρt × 2 (13) r˜u = rñ 4 f + ρ˜ f , n > 1 (21) ≡ n − 4 q n r˜u ρ˜u I(r˜u) would decrease by a factor of 0.215, and L/L⋆, iav and capture 1.5 1.080 0.675 68.16 efficiency would increase to values 9.30, 0.9, and 11.3, respec- 1.7 1.120 0.765 70.10 tively. 2 1.185 0.90 71.088 3 1.456 1.35 71.752 4.2. Synchronism mismatch 4 1.782 1.80 71.818 If the S/C reaches periapsis beforeor aftercrossing the meridian plane containingthe dipole, Equation (14) is no longervalid, with capture drag-work being the sum over drag work for pos- Table 1: I(r˜ ) for different values of n u itive (negative) true-anomaly θ ranges, which are differently drag-effective. Full work is the same, however, whether the The integral I(r˜u), —which is not singular as moving to a vari- dipole lies at θ positive or negative; for definiteness, set the abler ˜ 1+z2 would manifest— converges rapidly, as shown in dipole at a positive value α. Equation (17) takes now the form ≡ Table 1 and Figure 5, allowing integration to stop at r 1.12rp r˜ r˜ in the parabolic orbit, with a value around 70.1. ≈ u d˜r u d˜r I(r˜u) 6 + 6 = ≡ 1 ˜ √ 1 ˜ √ Z ρ r˜ 1 Z ρ+ r˜ 1 80 − − − r˜u d˜r ρ˜ 6 70 = 6 1 + −6 (25) 1 ρ˜ √r˜ 1 ρ˜+ Z − − 60 where 50 2 2 √ I(r˜u) 40 ρ˜ = (r˜ + f ) 4 f + 2 f (2 r˜)(1 cosα) 2 r˜ 1sinα ± − − − ± − 30 (ρ˜ ρ/RN, r˜ r/hRN) (26)i ≡ ≡ 20 + (-) signs corresponding to distances from the dipole to points 10 in the parabolic orbit of common r value and negative (posi- PSfrag replacements tive) true anomaly. The parenthesis inside the integral starts at 0 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2 a value of 2 at the lower limitr ˜l = 1, and decreases fast asr ˜ r˜ u increases. At that lower limitr ˜l = 1 Equation (26) gives a common Figure 5: I(r˜u). Points correspond to n = 2,3,4 value, For an aluminum tether, Equation (20) then provides maxi- 4 f sin2(α/2) ρ˜l = (1 f ) 1 + (27) mum capture-efficiency MSC/mt = 5.35 iav r. Similar to the − s (1 f )2 Saturn research [6], efficiency can be increased× h i by a gravity − Let us nowconsider upper limits for the integralin (25),r ˜ r˜ , assist from Jupiter through a tail flyby, increasing the heliocen- u = n corresponding in (26) to ρ˜ nρ˜l n > 1. Results for the work tric S/C velocity at the encounter with Neptune —which would − ≡ be greater than the velocity from a direct Hohmann transfer. integrals are given in Table 2, showing that capture-efficiency This reduces the hyperbolic velocity relative to the planet by decrease would be about 1.8 % for α = 5◦. Decrease is 6% for a 1.53 factor, from 3.96 to 2.59 km/s. Capture efficiency then α = 10◦.

increases in Equation (20) by a 2.34 factor to 12.52 iav r. MSC × h i n r˜u ρ˜u I(r˜u) m decrease Regarding the average dimensionless current iav(L/L⋆), the t 1.5 1.100 0.682 68.11 1.8 % expression for L⋆ in Equation (1) reads, for aluminum [19], 1.7 1.144 0.774 69.42 1.7 % 1/3 2/3 (meEm) σt ht 2 1.214 0.909 70.19 1.6% L⋆ 3π (22) ≡ 27/3e Ne 1/3 5 3 2/3 Em ht 10 /cm Table 2: I(r˜u) for different values of n. Small mismatch α = 5◦ L⋆ 2.78km (23) ≈ · 100V/km · 100 µm · N e

Because the drag-arc is short we set Em(r) Emp throughout ≈ the r-integration. Using 5. Discussion of Results

0.68 0.13 4 2 4 2 The result of our analysis is the determination of capture Bp = × 3 (10− Vs/m )= 0.97 10− Vs/m (24) 0.45 × × efficiency, defined as a ratio between mass MSC of captured S/C and mass m of tether required for the capture (includedin M ). we find that E = v B = 22.7 100 V/km. Setting h = t SC mp p p t We found that the large offset of the Neptune dipole resulted 10 µm, N = 103 cm 3, yields L ×= 36.6 km. A tether length e − ⋆ in it definitively exceeding the Saturn capture efficiency, thus L = 73.2 km leads to iav(2)= 0.52 and to an efficiency 6.51 4 3 compensating for the negative effect of the large tilt and the in the case of a Jupiter flyby. For a density of 10 cm− , L⋆ 5 3 3 low value of magnetic moment 0.13 GRN against 0.21 GRS. 10. REFERENCES As in the Saturn case [6], a gravity assist from Jupiter helped an efficient capture. References High efficiency in the presence of a large dipole-offsetwould [1] J. R. Sanmartin, E. C. Lorenzini, Exploration of outer planets using teth- require the S/C to arrive at periapsis —at a point very close to ers for power and propulsion, Journal of Propulsion and Power 21 (3) the planet— lying in the meridian plane of the Neptune dipole. (2005) 573–576. doi:10.2514/1.10772. That synchronismwas foundto be reasonablynecessary, whether [2] Ice giants – pre-decadal survey mission study report (2017), JPL D - 100520 the S/C reached that plane following, or ahead of, arrival at pe- http://www.lpi.usra.edu/icegiants/mission study/Full-Report.pdf, last riapsis, a 10◦ mismatch resulted in a 6 % decrease in efficiency. acces July 2019, (2001). Our calculations used several simplifications. One was ig- [3] J. R. Sanmartin, M. Charro, E. C. Lorenzini, H. Garrett, C. Bombardelli, noring the planet rotation while the S/C moved over the quite C. Bramanti, Electrodynamic tether at Jupiter – I: Capture operation and constraints, IEEE Transactions on Plasma Science 36 (5) (2008) 2450– limited parabolic arc contributing to drag. This arose from dis- 2458. doi:10.1109/TPS.2008.2002580. tinctive properties of Neptune among the 4 Giant Planets. First, [4] J. R. Sanmart´ın, M. Mart´ınez-Sánchez, E. Ahedo, Bare wire anodes for along with Uranus, it has a slow spin as compared with Gas Gi- electrodynamic tether, Journal of Propulsion and Power 9 (3) (1993) 353– ants. Secondly, it has the highest density; for a parabolic orbit 360. doi:10.2514/3.23629. [5] J. R. Sanmart´ın, M. Charro, H. B. Garrett, G. Sánchez-Arriaga, with periapsis close to a planet, the characteristic time of S/C A. Sánchez-Torres, Analysis of tether-mission concept for multiple flybys motion is shorter for higher densities. of moon Europa, Journal of Propulsion and Power 33 (2) (2017) 338— Another simplification involved the dipole model used. The 342. doi:10.2514/1.B36205. [6] J. Sanmart´ın, J. Peláez, I. Carrera-Calvo, Comparative Saturn-versus- dipole lies 0.19RN below the equatorial plane; we ignored this Jupiter tether operation, Journal of Geophysical Research: Space Physics fact. Only the component of field B along the rotational axis 123 (7) (2018) 6026–6030. doi:10.1029/2018ja025574. contributes to drag in an equatorial orbit. Our approximation [7] J. W. Belcher, H. S. Bridge, F. Bagenal, B. Coppi, O. Divers, A. Eviatar, excluded the first term in Equation (4) from contributing to G. S. Gordon, A. J. Lazarus, R. L. McNutt, K. W. Ogilvie, J. D. Richard- son, G. L. Siscoe, E. C. Sittler, J. T. Steinberg, J. D. Sullivan, A. Szabo, drag, because the dipole and S/C laid in the equatorial plane. L. Villanueva, V. M. Vasyliunas, M. Zhang, Plasma observations near With the dipole below the equatorial plane, the first term in neptune: Initial resultas from voyager 2, Science 246 (1989) 1478–1483. Equation (4) would contribute to drag. doi:10.1126/science.246.4936.1478. Furthermore, quadrupole terms, which are more important [8] D. A. Gurnett, W. Kurth, I. H. Cairns, L. Granroth, Whistlers in neptune’s magnetosphere: Evidence of atmospheric lightning, Journal of than in Jupiter and Saturn cases because of the proximity of Geophysical Research: Space Physics 95 (A12) (1990) 20967–20976. the magnetic moment to the planet surface, were fully ignored. doi:10.1029/ja095ia12p20967. Spherical harmonic analysis first provided roughly resolved har- [9] G. F. Lindal, The atmosphere of neptune-an analysis of radio occultation monic coefficients describing the quadrupole [15], then reached data acquired with voyager 2, The Astronomical Journal 103 (1992) 967– 982. doi:10.1086/116119. a consistent description [16], [17]. [10] J. Richardson, J. Belcher, A. Szabo, R. McNutt Jr, The plasma environment of Neptune, Neptune and Triton (1995) 279–340 Edited by D. P. Cruikshank and M. S. Matthew, University of Arizona Press, Tucson. 6. Conclusions [11] J. Bishop, S. K. Atreya, P. N. Romani, G. S. Orton, B. R. Sandel, R. V. The high capture efficiency results suggest new calculations Yelle, The middle and upper atmosphere of Neptune, Neptune and Triton (1995) 427–487. to includeNeptunerotation and to use the OTD2 detailed dipole- [12] D. Gurnett, W. Kurth, Plasma waves and related phenomena in the mag- model, considering further quadrupole corrections; consider- netosphere of Neptune, Neptune and Triton (1995) 389–423. ation of planet rotation would require revising the sensitivity [13] L. D. Landau, E. M. Lifshitz, The classical theory of fields 16, secs. 43 of capture efficiency to synchronism mismatch. The upgraded and 44, Pergamon Press, Oxford (1962). [14] N. F. Ness, M. H. Acuna, L. F. Burlaga, J. E. Connerney, R. P. Lepping, analysis would allow the start of planning a visit to moon Triton F. M. Neubauer, Magnetic fields at Neptune, Science 246 (4936) (1989) mission. 1473–1478. doi:10.1126/science.246.4936.1473. [15] J. Connerney, M. H. Acuna, N. F. Ness, The magnetic field of Neptune, Journal of Geophysical Research: Space Physics 96 (S01) (1991) 19023– 7. Competing Interests Statement 19042. doi:10.1029/91ja01165. [16] N. F. Ness, M. H. Acuña, J. E. Connerney, Neptune’s magnetic field and There are no competing interests amongst the authors. field-geometric properties, Neptune and Triton (1995) 141–168. [17] J. Connerney, Magnetic fields of the outer planets, Journal of Geophysical Research: Planets 98 (E10) (1993) 18659–18679. 8. Funding doi:10.1029/93je00980. This work was supportedby MINECO/AEI and FEDER/EU [18] C. Russell, M. Dougherty, Magnetic fields of the outer planets, Space science reviews 152 (1-4) (2010) 251–269. under Project ESP2017-87271-P. doi:10.1007/s11214-009-9621-7. [19] E. Ahedo, J. Sanmart´ın, Analysis of bare-tether systems for deorbiting low-earth-orbit satellites, Journal of Spacecraft and Rockets 39 (2) (2002) 9. Acknowledgements 198–205. doi:10.2514/2.3820. The authors thank the MINECO/AEI of Spain for their fi- nancial support. Part of J. Pelez’s contribution was made at the University of Colorado at Boulder on a stay fundedby MICINN (ref. PRX19/00470) and the Fulbright Commission, to which he thanks the support received. 6