Access the full text.
Sign up today, get DeepDyve free for 14 days.
References for this paper are not available at this time. We will be adding them shortly, thank you for your patience.
Breakthrough Starshot is an initiative to prove ultra-fast light-driven nanocrafts, and lay the foundations for a rst launch to Alpha Centauri within the next generation. Along the way, the project could generate important supplementary bene ts to solar system exploration. A number of hard engineering challenges remain to be solved before these missions can become a reality. A system model has been formulated as part of the Starshot systems engineering work. This paper presents the model and describes how it computes cost-optimal point designs. Three point designs are computed: A 0:2 c mission to Alpha Centauri, a 0:01 c solar system precursor mission, and a ground-based test facility based on a vacuum tunnel. All assume that the photon pressure from a 1:06 μm wavelength beam accelerates a circular dielectric sail. The 0:2 c point design assumes $0:01=W lasers, $500=m optics, and $50=kWh energy storage to achieve $8.0B capital cost for the ground-based beam director. In contrast, the energy needed to accelerate each sail costs $6M. Beam director capital cost is minimized by a 4:1 m diameter sail that is accelerated for 9 min. The 0:01 c point design assumes $1=W lasers, $10k/m optics, and $100=kWh energy storage to achieve $517M capital cost for the beam director and $8k energy cost to accelerate each 19 cm diameter sail. The ground-based test facility assumes $100=W lasers, $1M/m optics, $500=kWh energy storage, and $10k/m vacuum tunnel. To reach 20 km s , fast enough to escape the solar system from Earth, takes 0:4 km of vacuum tunnel, 22 kW of lasers, and a 0:6 m diameter telescope, all of which costs $5M. The system model predicts that, ultimately, Starshot can scale to propel probes faster than 0:9 c. Keywords: Breakthrough Starshot, beam-driven sail, beamed energy propulsion, model-based systems engineering 1. Objectives Robert Forward, also of Hughes Research Labs, proposed to use the photon pressure from lasers to accelerate sails Breakthrough Starshot. Send 1 gram of scienti c instru- to speeds that bridge interstellar distances [2]. In his 1984 mentation to the Centauri System to study it. Image its paper [3], Forward suggests coherently combining many planets, look for life and transmit the results to Earth. Do lasers into a phased array to reach high enough power and so by using a beam emitted from the Earth to accelerate large enough primary optic sizes. He correctly identi es a sail carrying the instrumentation to 0:2 c. The capital the key gure of merit as the sail acceleration produced cost of the equipment shall be less than $10B. per unit laser power. This gure of merit favors a sail material that is thin, light, re ective, and has low beam Systems Engineering Work. Ensure that Starshot engi- absorptance. However, Forward considers only metallic neering activities amount to a mission that ful lls the sails, which are absorption-limited to uxes of less than Breakthrough Starshot objectives. 10 kW m . The resulting 3:6 km diameter, 1 t sail accel- System Model. Replace physical experiments with simu- erates to 0:11 c over a period of three years using a 65 GW lations in cases where it saves time and money to do so. laser whose primary optic is 1000 km in diameter. Verify Breakthrough Starshot feasibility and estimate per- In his 1986 paper [4], Forward makes key conceptual formance. Design, optimize, trade-o, and analyze alter- progress by suggesting that sails be made from multilayer natives. Generate and quantify requirements. Model the dielectrics instead of metals. Dielectrics have more than impact of changes. ve orders of magnitude lower absorption than metals be- cause their electrons are bound electrons that cannot ab- 2. Context sorb photons except at discrete wavelengths. In compari- son, metals absorb energy from an incident beam via the The rst operational laser was built by Maiman of mechanism of Joule heating. The metal's softening tem- Hughes Research Labs in 1960 [1]. Only two years later, perature limits the irradiance on the sail, in turn limiting its acceleration. Switching from metals to dielectrics pro- Systems Director, Breakthrough Starshot foundly increases the irradiance and acceleration of the Email address: kevin@parkinresearch.com (Kevin L. G. Parkin) Preprint submitted to arXiv July 27, 2018 arXiv:1805.01306v3 [astro-ph.IM] 26 Jul 2018 sail, which in turn decreases the timescale, beam size, 3. Formulation and cost. Using eight quarter-wave thick layers of dia- The system model is formulated around the propaga- mond, Forward's 6:4 cm diameter, 4:5 mg sail accelerates tion of a beam from a ground-level beamer to a sailcraft to 0:015 c over a period of 2:6 s using a 1 GW laser whose in space above it, as shown in g. 1. The sailcraft begins primary optic is 100 m in diameter. at a given initial displacement above the beamer. This The dielectric sail concept was improved in 1989 when displacement, in combination with the beamer diameter, Landis pointed out that a single quarter-wave thick dielec- is used by a beam propagation model to determine the tric lm accelerates faster than a multilayer lm despite fraction of transmitted power that reaches the sailcraft. A having lower re ectance [5]. In unpublished work, Landis material/optical model calculates how much of the power went on to nd that because of the way that re ectance that is incident on the sailcraft is re ected or absorbed. A scales with thickness, the true optimum sail thickness is relativistic equation of motion then translates this power somewhat less than a quarter wavelength. His result is into an acceleration. The equation of motion is numeri- con rmed later in this paper. cally integrated forward in time until the sailcraft reaches In work leading directly to the formation of Break- its desired cruise velocity. The last photons arriving at the through Starshot in early 2016, a Harvard team led by sailcraft are traced back in space and time to determine Loeb [6] used a simple spreadsheet model to verify the when beam cuto occurs at the beamer. The system pa- feasibility and performance of a beam-driven sail and to rameters of beamer diameter, beamer power, and sailcraft produce point designs. A UCSB team led by Lubin [7] de- diameter are optimized to ensure that the sailcraft actually veloped a comprehensive analytical model that sought to reaches cruise velocity and does so using a minimum-cost describe the leading order system behavior. To use closed- beamer. This cost optimization reduces the dimensionality form equations and avoid numerical trajectory integration, of the model because the beamer diameter, beamer power, both the Harvard and UCSB models used simplifying ap- and sailcraft diameter are no longer inputs. proximations including top-hat beams. To the extent that cost was minimized, it was minimized by manual experi- 3.1. Goubau Beam mentation with the model's input values. In the general case, beam propagation from a ground- The present system model was formulated in March level beamer to a space-based sailcraft involves models 2016 and set out to verify the assumptions and results of representing a phased array, its elements, and transat- the earlier models. Key questions at this stage were: Are mospheric propagation. This system model simpli es the the earlier models correct? What diameter is the sailcraft beamer to an eective primary optic that transmits an and beamer ? How much smaller/cheaper can the beamer idealized beam. The system model makes no representa- be if the sailcraft continues to accelerate until the beam be- tions about the technologies used in the array elements, comes too weak? How cheap must the lasers/microwaves the element sizes, or how the phased array as a whole is be for the beamer to cost less than $10B? implemented. Nor does the system model make detailed estimates of atmospheric attenuation. Instead, all these factors are rolled into a user-provided value for stored to sail including payload beam director transmitted power eciency, and a user-provided value for atmospheric attenuation. The system model does calcu- late the idealized beam transfer eciency, which varies the most of these three eciency factors. 𝑧 , 𝑡 , …………………. SAILCRAFT, BEAM CUTOFF 2 2 TRAJECTORY (1-DOF RK45, RELATIVISTIC) MATERIAL/OPTICAL (STRATIFIED LAYER) 𝐷 , 𝑧 ,𝑡 ,𝑃 , 𝜂 , 𝜌 , 𝑅 , 𝑇 , 𝐴 SAILCRAFT, BEAM INITIALIZATION 𝑠 1 1 𝑏 𝑎 𝑧 BEAM (IDEALIZED GOUBAU) BEAMER 𝐷 $/m , $/Watt, $/kWh COST MODEL Figure 1: System model 2 Sailcraft Unlike Gaussian beams and top-hat beams, Goubau Equations (1) to (4) enable the rapid recalculation of beams [8, 9] describe near-optimal energy transfer between power reaching the sail on each time-step of the trajectory nite optics [10]. For this reason, Goubau beams are integration via used in the context of wireless power transfer [11]. When P = P ; (5) b a b 1 high transfer eciency is needed, the beam pro le resem- bles a Gaussian, and when low transfer eciency can be where is the eciency factor that accounts for atmo- tolerated, the pro le resembles a top-hat. At intermedi- spheric attenuation via absorption and scattering. For ate eciencies, the beam resembles a truncated Gaussian. consistent accounting throughout this paper, P is de ned But which transfer eciency minimizes the beamer cost? to be the laser power that is transmitted by the beamer. P Assumptions about the answer to this question are often is therefore the fraction of transmitted power that is des- wrong. tined to reach the sailcraft. It varies monotonically with Referring to g. 1, the system model represents the , a desirable feature when using numerical optimization. beamer and sailcraft as two areas that are perpendicular to a common axis. Transfer eciency is a function of a single 3.2. Equation of Motion dimensionless parameter that depends on the product The equation of motion relates the power that is inci- of optic diameters D and D , the distance between them s b dent upon the sailcraft to the sailcraft's acceleration. It is z, and the beam wavelength : deduced by requiring momentum to be conserved through z 8z the interaction of the sailcraft and beam photons. The 2 = : (1) D D derivation presented here extends the approach of Kulka- A A s b s b rni [13] to include a dielectric sailcraft having nite trans- There is no closed-form solution for beam power transfer mittance. Also, the sailcraft thermally re-emits absorbed eciency ( ); however, it is closely approximated by [12]: beam energy in the forward and backward directions, and this is included here because it contributes a nonzero drag (a) if a > 1:21748051194181 (a) = (2) force as seen from the beamer frame. (a) otherwise In the beamer (observer) frame, conservation of mo- 1 mentum can be expressed as 8 4 2 (a) = a + a 4a b + 4b 8b + 4 (3) 4b dp dp p s + = 0; (6) 2 6 10 dt dt a a 7a (a) = + ; (4) where the relativistic sailcraft momentum p is given by 2 32 4608 s 2 p = m c : (7) s 0 2 a where a ( ) and b e . c is the speed of light, m is the sailcraft rest mass, and is the Lorentz factor 1= 1 . Newton's second law shows the time rate of change of sailcraft momentum to be equal the apparent force F acting on the sailcraft as seen from the beamer frame, 0:8 dp F = = m c : (8) s 0 dt 0:6 Reflected 0:4 𝑛 ,𝑝 𝑟 𝑟 0:2 Incident Transmitted 𝑛 ,𝑝 𝑛 ,𝑝 𝑠 0 𝑡 0 0 5 10 15 20 Re-emitted backward Re-emitted forward 𝑛 , 𝑝 e r 𝑛 , 𝑝 e 0 Figure 2: Goubau beam power transmission eciency calculated us- ing eq. (2) Figure 3: Control volume De ned here as Goubau de nes it [8, 9]. Others de ne it dier- ently [10, 11]. The time rate of change of photon momentum is calculated using a control volume that co-moves with the sailcraft, ( ) b shown in g. 3. Photons cross the control boundary with Finally, substituting eqs. (6) to (13) into this expression rates n and momenta p. In the gure, the arrows are drawn yields the equation of motion: such that the rates and momenta are positive quantities. P 1 Of the photons which are destined to strike the sailcraft F = (A + 2R) ; (15) c 1 + (not all or even most of them), let us say that n is the rate at which they are transmitted by the beamer. The where P = cp n represents the beam power that is des- b 0 b sailcraft recedes from the beamer at speed . Thus, the tined to strike the sailcraft as measured in the beamer rate at which photons strike the sailcraft, n , is given by s frame. At any point in the trajectory, P is calculated from eq. (5). Equation (15) can be simpli ed by recogniz- n = n (1 ) : (9) s b ing that in the sailcraft rest frame, the sailcraft receives a power P given by [14] Consistent with the other rates, this rate is the number of s photons per unit beamer time, as opposed to sailcraft time. P = P : (16) This is important because there is time dilation between b 1 + the beamer and sailcraft frames. The rate n at which photons are re ected is related to the re ectance R by Thus, n = Rn : (10) s r s F = (A + 2R) : (17) This model asserts that the sail is in thermal equilibrium Combining eq. (15) with eq. (8) yields another form of the because it re-emits the energy that it absorbs and has no equation of motion that is more readily integrated, thermal inertia. Thus, the rate n at which photons are transmitted is given by P 1 = (A + 2R) ; (18) E 1 + n = (1 R A) n ; (11) t s where E = m c is the sailcraft rest energy. 0 0 where A is the absorptance. Absorbed photons heat the If the derivation of force is repeated for only the re- sailcraft. To maintain thermal equilibrium, their energy emitted forward and backward photons shown in g. 3, is re-radiated in the forward and backward directions as then a drag-like force is found: photons spanning a range of wavelengths. For accounting purposes, this model pretends that the incident photons F = ; (19) that would have been absorbed are instead re ected and transmitted in equal measure, in a similar way to the other where the absorbed power P is given by photons. Thus, the photons that are re-emitted are as- 0 0 signed momentum p in the backward direction and p in r 0 P = AP : (20) a s the forward direction. The rate n at which photons are Physically, the thermally-radiated photons blueshift in the re-emitted is related to the absorptance by forward direction and redshift in the backward direction. This unbalances photon momenta between the two direc- n = n : (12) e s tions and manifests as an apparent drag that acts on the Photon momenta are unchanged by the moving control sailcraft from the beamer frame, yet is not felt by the sail- volume, but they are aected by frame conversions within. craft in its rest frame. In the context of dust particles that Photons that are re ected from the sail are Doppler shifted orbit stars, this force is known as Poynting-Robertson drag twice: First in going from the beamer to sailcraft frame, [15]. For a non-absorbing sailcraft, there is no drag, but for then after re ection in going back from the sailcraft to the a perfectly absorbing sailcraft, the equations of Kulkarni beamer frame. The photon is redshifted both times be- [13, 16] overpredict the sailcraft's acceleration by a factor cause the sailcraft is receding from the beamer and vice of (1 + ), which is 20% at 0:2 c. versa. In the beamer frame, an incident photon has mo- This derivation has said that photons strike the sail- mentum p . This is related to the beamer wavelength craft at a rate that is slowed by the classical factor of 0 0 by p = , where h is Planck's constant. Thus, the ratio (1 ), so where do the missing photons go? The situa- of re ected to incident momentum is given by the inverse tion is clari ed by drawing a spacetime diagram, shown in square of the relativistic Doppler factor: g. 4. The missing photons simply catch up with the sail- craft after the beam is turned o. Thus, if a sailcraft accel- p 1 r 0 erates from rest until it reaches the desired cruise velocity = = : (13) p 1 + 0 r at time t and displacement Z , then the beam duration t s b is given by Summing the contributions at the boundary in g. 3, dp t = t : (21) = n p n p + n p + n p n p : (14) b s s 0 t 0 r r e r e 0 dt 4 𝑡 Recalling the equation of motion, eq. (18), and simplifying with A R and 1 yields, Extra time Δ𝑡 = / ; (25) because ! 1 and E / where is sail mass density and is its thickness. Thus, it can now be seen that the 100 nm thick sail is better and accelerates 2.5 times faster than the alternative, all other things being equal. The question of optimum sail thickness and perfor- mance is further complicated because dielectric layer re- ectance varies with thickness. The Starshot system model implements a complex characteristic matrix method as de- scribed by Macleod [17] to calculate the re ectance, trans- mittance and absorptance of an arbitrary but locally- at assembly of thin lms. 0.15 0.25 0.5 0.75 1.0 0.35 0.3 Figure 4: Spacetime diagram 0.25 0.2 The sailcraft equilibrium temperature T is estimated by 0.15 using the Stefan-Boltzmann equation, 0.1 0.05 T = ; (22) 1.8 where is the Stefan-Boltzmann constant and " is the 1.6 2-sided total hemispherical emittance. This means that 1.4 for a temperature-limited sailcraft, P is constant. In the s 1.2 sailcraft rest frame, the equation of motion is simply 0.8 0 s 0.6 F = (A + 2R) : (23) 0.4 0 0.2 Thus, F is constant for a temperature-limited sailcraft. Equating eq. (23) with eq. (17) yields simply 0.1 1 0 δ/λ [-] F = F : (24) Figure 5: Figures of merit predicted by the strati ed layer model Thus, the force is the same in the sailcraft and beamer frames. This is consistent with a pure force that does not In g. 5, the model is used to plot the re ectance of add net heat to the sailcraft. Other forms of the equation a single-layer dielectric lm having a relative refractive of motion produce impure forces that are frame-dependent index of 2 under illumination by free-space plane waves and imply sailcraft heating. oriented at normal incidence to the lm. Thus, the in- The acceleration varies between sailcraft and beamer cident waves have half their free-space wavelength within frames because of the factor of introduced by eq. (8). the lm, and the maximum re ectance occurs when the This variation is seen in the trajectory integration results, lm is one quarter of this thickness, three quarters, ve presented a little later. quarters, and so on. As already shown, system perfor- mance is maximized by choosing the thickness at which 3.3. Strati ed Layer Optical Model R= is maximum. Equation (10) shows this to occur at Which is better, a 1000 nm thick sail with 99.9% re- about =7, with the subsequent maximum having more ectance, or a 100 nm thick sail of the same mass density, than three times worse performance. For illumination at but with 25% re ectance? To de ne `better', a gure of a free-space wavelength of 1:06 μm, this corresponds to an merit is needed. For simple design purposes, the appro- optimum sail thickness of 76 nm. priate gure of merit is not the sail re ectance or the sail thickness, but the sail acceleration per unit beam power. BEAM ON SAILCRAFT Rλ/δ [-] R [-] δ=λ/7 4. Key Trade from the beamer to the sailcraft via increased beamer di- ameter and decreased beam duration. It is always possible to vary beam diameter, power, A simple model for beamer capex C is given by and duration such that a sailcraft reaches its desired cruise velocity. Thus, the key trade lies in reaching the desired C = k + k P + k Q : (26) a l 1;max e 0 cruise velocity at minimum cost. Beamer diameter D , peak transmitted power P , and b 1;max Lasers pulse duration t are dependent variables of the system model, as explained in the next section. Stored energy Q is the integral of the wallplug power drawn over the pulse duration, and this power draw is not constant, as shown later. Factors k , k , and k are independent user-supplied a l e values for cost per unit area, cost per unit power, and cost per unit energy stored. They are technology gures of minimum merit. By choosing a high laser cost factor of $1000=W, for capex example, the cost-optimum solution moves toward longer beam duration, reduced average power, and larger diam- eter, to improve energy transfer eciency at long range. At present, the cost model encompases the beamer only and not the sailcraft or other elements of the system. The beamer is expected to be the largest capital expense. Storage Optics Figure 6: Key trade 5. Solution Procedure Benford [18] presents an analysis for cost-optimized A solution procedure for the system model is shown in beam driven sail missions in which the beamer's primary g. 7. At its top levels, the system model performs nested optic diameter is traded vs. power in order to minimize optimizations to ensure that the sailcraft reaches its cruise capex : A smaller diameter saves money on optics but velocity and that it does so with system elements whose reduces transfer eciency, so power is increased to com- speci cations minimize the beamer capex. pensate. Conversely, lower power saves money on sources but transfer eciency has to be improved to compensate, + ± ± 𝛽 ,𝐷 , 𝑃 , 𝐷 Inputs 𝑏 1 𝑠 Input values by enlarging the primary optic. Somewhere between a beamer that is very powerful with a tiny primary and one 1 :2 Vary 𝐷 such that𝐶 is minimized + + ± ± Given 𝛽 ,𝐷 , 𝑃 , 𝐷 𝑏 1 𝑠 that is very weak with a gigantic primary, there exists a (golden section search, GSS) happy medium that minimizes capex. Benford's results af- :1 ± Vary 𝑃 such that𝐶 is minimized + + Given 𝛽 , 𝐷 , 𝑃 , 𝐷 rm a rule of thumb used by microwave system designers 𝑠 𝑏 1 (golden section search, GSS) for rough estimates: Minimum capex is achieved when the 3a + Vary𝐷 such that𝛽 = 𝛽 + + cost is equally divided between antenna gain and radiated Given 𝛽 , 𝐷 , 𝑃 , 𝐷 𝑏 1 (bisection) power. Similar results have been obtained in the eld of a4 Integrate sailcraft trajectory to minimum beamed energy launch vehicles [12, 19]. Hereafter, we refer + Given 𝛽 , 𝐷 , 𝑃 , 𝐷 𝑏 1 𝑠 acceleration threshold (integration, RK45) to such analyses as \traditional beam-aperture tradeos". The \traditional beam-aperture tradeo " optimum sits 3b Vary 𝐷 such that 𝐶 is minimized on the axis between lasers and optics in g. 6. However, Given 𝛽 , 𝐷 , 𝑃 , 𝐷 1 𝑠 (golden section search, GSS) this is not a global minimum capex. This paper incorpo- b4 rates the key realization that beam duration is an indepen- Integrate sailcraft trajectory to 𝛽 = 𝛽 Given 𝛽 , 𝐷 , 𝑃 , 𝐷 𝑏 1 (integration, RK45) dent variable that can be traded vs. power and diameter to lower the capex even further. Each of the three costs shown in g. 6 can be traded for increases in the others: Dominant laser cost is mitigated by reducing power in fa- :1 :1 :1 𝑄 𝐶 vor of longer beam duration (increases losses) and larger :2 :2 :2 :2 beamer diameter (decreases losses). Dominant optics cost 𝐷 𝑃 𝑄 𝐶 𝑏 1 is mitigated by reducing beamer diameter (increases loses) in favor of higher power. Dominant energy storage cost is :3 :3 :3 :3 :3 Outputs 𝐷 𝑃 𝐷 𝑄 𝐶 𝑏 1 𝑠 mitigated by improving the eciency of energy transfer Figure 7: Solution procedure capital expenditure 6 Starting at the top of g. 7, the desired cruise veloc- that C=A is minimized (instead of P varied such than C b 1 ity is speci ed along with upper and lower bounds is minimized). Also, would be optimized instead of D . for the beamer power, P , and sailcraft diameter D . Given cost-optimal values for C=A , I , and I dt, a b 1 1 1 s Only a maximum beamer diameter D is speci ed. These separate optimization could then vary D to minimize C . values are brought into the outermost iteration (iteration This optimization would iterate the Goubau and cost mod- 1), a golden section search varying D between its chosen els without recalculating the computationally-expensive :2 bounds D such that beamer capex C is minimized. The trajectory. This amounts to a `traditional' power vs. aper- `:2' in the superscript means that C is already twice opti- ture tradeo as shown in g. 6. Thus, it is inferred that mized (with respect to P and D ). Shown at the bottom storage is an axis of the tradespace only because payload 1 b of the gure, iteration 1 returns the cost-optimal sailcraft mass is not a constant fraction of sail mass. Also, con- :3 :3 diameter D together with the corresponding capex C , stant costs that do not scale with beamer area C=A in now optimal over three dimensions. Before these optimal cost eq. (26) would again unfold storage as an axis of the values are returned, all the internal optimizations must tradespace. rst run, and these draw on the optimizer-chosen value of D , which is passed inwards to iteration 2. 6. Point Designs Iteration 2 is a golden section search varying P be- :1 tween its chosen bounds P such that capex C is mini- 1 6.1. 0:2 c Mission mized. The `:1' in the superscript means that it is already The 0:2 c point design embodies key elements of mid- optimized with respect to D only. This iteration returns 21st century Starshot missions to nearby stars, including :2 :2 P . It also returns the corresponding capex C to it- those to the Centauri System. eration 1, which uses its value to guide the optimizer in its choice of D for the next iteration. The newly-de ned 6.1.1. Inputs value of P is passed inwards to iterations 3a and 3b. The inputs to the mission point design are summa- Iteration 3a is a bisection solver varying D between + rized in table 1. A 1:06 μm wavelength is consistent with zero and its upper bound of D such that sailcraft velocity ytterbium-doped ber ampli ers. An initial sailcraft dis- at the end of trajectory integration, , is only just equal placement of 60 000 km is consistent with a low-thrust non- to its desired value of . By setting D in this way, the Keplarian orbit [20, 21] that keeps the sailcraft (and as- solution procedure assures that beam cuto always occurs sociated spacecraft) stationary in the sky relative to the as the sailcraft reaches exactly its desired velocity . It target star. remains to minimize the capex. Iteration a4 is a trajectory integration (implemented using the RK45 algorithm) returning the sailcraft velocity Table 1: System model inputs for 0:2 c mission at the point that its acceleration falls below a minimum 0:2 c target speed threshold, which terminates integration. 1:06 μm wavelength Iteration 3b is a golden section search varying D be- 60 000 km initial sail displacement from laser source tween its chosen upper bound D and its calculated lower bound D such that capex C is minimized. As described 1 g payload by cost eq. (26), C is an implicit function of the sailcraft 0:2 g m areal density trajectory, so this trajectory must be recalculated on each 10 spectral normal absorptance at 1:06 μm :1 iteration. Iteration 3b returns D and its corresponding 70% spectral normal re ectance at 1:06 μm :1 :1 capex C to iteration 2, which uses C to guide the opti- 625 K maximum temperature mizer in its choice of P for the next iteration. 0.01 total hemispherical emittance (2-sided, 625 K) Iteration b4 is a trajectory integration (implemented using the RK45 algorithm) that returns the pulse energy $0:01=W laser cost Q because it is needed in cost eq. (26) to compute energy $500=m optics cost storage cost. $50=kWh storage cost Finally, the result of all these iterations are values for 50% wallplug to laser eciency :3 :3 :3 three dimensions; D , P , and D ; that produce a min- 70% of beam power emerging from top of atmosphere 1 s :3 :3 imal capex C and cost-optimal values for Q . All these values are used in conjunction with auxiliary equations to A 1 g payload is bookkept separately from sail mass and quantify point designs. reserved for scienti c instrumentation and associated sup- port systems. Sail mass is calculated by the system model Simpli ed method for constant payload fraction. If pay- using the value of D chosen by the optimizer combined load were a constant fraction of sail mass, the solution with the input areal density. If the sail were an optimal- procedure could be simpli ed by omitting the outer itera- thickness layer of silicon dioxide, for example, it would tion on D . Beamer intensity I would then be varied such s b 7 2 have an areal density closer to 0:3 g m , a room tempera- Table 2: System model outputs for 0:2 c mission ture absorptance of less than 10 [22], and a normal re- $8.0B beamer capex comprised of: ectance of 12%. This point design uses a lower areal den- $2.0B lasers (200 GW max. transmitted power) 2 8 sity of 0:2 g m , higher absorptance of 10 , and a higher $2.8B optics (2:7 km primary eective diameter) re ectance of 70%, consistent with a hot two-dimensional $3.1B storage (63 GWh stored energy) nanohole photonic crystal that is somewhat tuned for high re ectance [23]. Because the system model does not yet $6M energy cost per Starshot (63 GWh @$0:1=kWh) incorporate a photonic crystal model, the strati ed layer 2.9% system energy eciency optical model is turned o and absorptance and re ectance remain constant throughout trajectory integration. 4:1 m sail diameter A maximum temperature of 625 K is placed on the sail 3:6 g sailcraft mass (includes payload mass) to prevent thermal runaway and to prevent damage to non- operating electronic or photonic components that are part 8 min (480 s) beam transmit duration of the sail or its payload. An assumed total hemispherical 9 min (550 s) sailcraft acceleration duration emittance of 0.01 accounts for the energy that is radiated by the sail in both the forward and backward directions at 40 Pa temperature-limited photon pressure 625 K. This emittance is much lower than that of typical 520 N temperature-limited force materials because the sailcraft is so thin and because its 14 900 g s temperature-limited acceleration absorptance is nearly zero at 1:06 μm. In comparison, a 2500 g s nal acceleration (0:13 au), 67 ls from source blackbody at 625 K emits most strongly at 4:6 μm. This means that to approach the radiative performance of a 37 kW m beamer maximum beam radiant exitance blackbody, the sail must switch from virtually invisible to 8:7 GW m sailcraft temperature-limited irradiance strongly absorbing/emitting within as small a wavelength range as possible. Cost factors are chosen such that beamer capex is less only one fortieth of the 2:5 TWh annual market for elec- than $10B. In particular, the laser cost is $0:01=W, four or- trical energy storage that is projected for 2040 [24], with ders of magnitude lower than the present cost, consistent electric vehicles being the dominant source of demand. with nearly automated production and high yields. Au- Though it would be impractical to convene two million tomated production lines for microwave oven magnetrons electric vehicle owners and drain their batteries for each long ago reached this cost. The optics cost is $500=m , starshot, it may be possible to use second-hand or donated three orders of magnitude lower than the present cost for battery packs that no longer perform well enough for elec- diraction-limited optics, again consistent with nearly au- tric vehicles. tomated production and high yields. This value is compa- The remainder of the beamer capex is incurred by the rable with the retail cost of computer monitors. $50=kWh optics; a lled array of telescope elements that form a energy storage is greater than current materials costs for 2:7 km eective diameter. This eective diameter is the some, but not all, energy storage technologies, and is more same as that of the Crescent Dunes solar energy concen- optimistic than current experience curve projections of fu- trator, located in the California desert. Of course, the ture electrical energy storage costs [24]. array modules for the beamer will be very dierent from solar concentrator mirrors. When operating at maximum 6.1.2. Results 2 power, the beamer's radiant exitance is 37 kW m (a spa- Upon running the system model using the inputs in tial average obtained by diving the power output by the table 1, the optimizers converge to the values given in ta- eective area of the primary optic). This exitance is forty ble 2. The sail diameter that minimizes beamer capex is times that of the solar concentrator operating at its peak, found to be 4:1 m, corresponding to a sailcraft mass of but three orders of magnitude lower than that of mili- 3:6 g. At 0:2 c, this mass has a relativistic kinetic energy tary laser beam directors. In the plane of the sail, the of 1:9 GWh, whereas the beamer expends 63 GWh, corre- beam converges to a much higher irradiance. By limit- sponding to 2.9% system energy eciency. This energy ing the beamer's power output, the system model reduces costs only $6M at a price of $0:1=kWh, making the en- the irradiance at the sail to its temperature-limited value ergy cost three orders of magnitude lower than the $8.0B of 8:7 GW m , which is three orders of magnitude lower beamer capex. than the ux of a 1 kW laser in a 10 μm ber, and ve The $8.0B beamer capex is comprised of three unequal orders of magnitude lower than the non-thermal ablation expenses with, surprisingly, laser cost being the smallest at threshold [25]. $2.0B. This is for a 200 GW maximum transmitted power, 2 The 8:7 GW m sail irradiance produces only 40 Pa which is far below the petawatt-class lasers now in exis- photon pressure, equivalent to a moderate breeze. But tence but many orders of magnitude greater in pulse en- the sail is very thin and the breeze moves at the speed ergy (and duration). The stored pulse energy of 63 GWh 0 of light, resulting in 14 900 g s initial acceleration. Such contributes the most to the beamer capex at $3.1B, but is 8 sailcraft 0 -1 10 1 10 Mars minimum distance 0.8 -1 -2 10 10 0.6 -2 -3 10 10 0.4 D86 beam diameter Lunar orbit -3 -4 10 10 0.2 Geostationary orbit -4 -5 10 0 10 0.2 1 10 0.95 0.15 10 0.9 0.85 once 0.1 10 0.8 0.75 2 0.05 10 Laniakea supercluster 0.7 twice escape velocity (0.01c) 0 0.65 10 16000 45 14000 P sailcraft rest frame beamer b 15 rest frame s 10 2000 s 0 0 10 650 600 2 kt TNT 8 600 6 550 4 500 2 450 0 400 0 0 1 2 3 4 5 6 7 8 9 10 0 1 2 3 4 5 6 7 8 9 10 0 1 2 3 4 5 6 7 8 9 10 Time in Earth frame [minutes] Time in Earth frame [minutes] Time in Earth frame [minutes] Figure 8: 0:2 c trajectory and related quantities acceleration is experienced by bullets and artillery shells, douse the beam if needed to protect unexpected aircraft but for fractions of a second. Even as the sailcraft reaches or ocks of birds. If the interruption were short enough, 0:2 c after 9 min, it is still accelerating at 2500 g s. Such is the trajectory could resume. the cost-optimum truncation point for the trajectory. The 1-D sailcraft trajectory and associated quantities During the sailcraft's acceleration, a satellite that tran- are plotted as functions of time in g. 8. This trajectory sits the beam in low Earth orbit would see a ash of lit- corresponds to the cost-optimal point design summarized tle more than 37 kW m for a fraction of second. This in table 2. Starting at the top left plot, the sailcraft begins is because the high- ux part of the beam is the focus, accelerating from a distance that would be just past the and the focus is initially past geostationary orbit and only radius of geostationary orbit. If the destination star were gets further away as the sail moves o. The beam and in the plane of the solar system, then the sailcraft would the ultra-high-acceleration sailcraft, behaving as if a speck pass the Moon's orbit within the rst minute, and could of dust in optical tweezers, may be able to dodge satel- be a third of the way to Mars by the end of acceleration lites in medium Earth orbits or supersynchronous orbits. around the 9 min mark. Alpha Centauri is at 60 dec- As is current practice, beam-satellite conjunction analyses lination, so its trajectory would be downwards out of the would be performed to help schedule times at which las- ing is allowable. Also, there would be interlocks to dim or at its minimum distance from Earth Relativistic kinetic energy [TJ] Acceleration [g’s] Speed, β [c] Distance from Earth [au] Temperature [K] Power [GW] Doppler shift [-] Unspilled fraction of beam, η [-] Force [N] Photon pressure [Pa] Fresnel number [-] Angular extent [arcsec] cost optimal minimum for 0.2 c at sailcraft at beamer 11 4 10 1 10 10 3 10 10 9 2 10 0.1 10 0.05 1000 1000 0.04 1,max 0.03 100 100 0.02 0.01 1,initial 0 10 10 10 100 650 temperature power limited limited regime regime 1 1 400 5 11 10 1 10 max. at sailcraft 4 8 10 10 3 5 10 10 2 2 10 0.1 10 1 10 1 10 1 10 Sail diameter [m] Sail diameter [m] Sail diameter [m] Figure 9: 0:2 c mission variation with sail diameter plane of the solar system. dently, it minimizes beamer capex to oversize the beamer The sailcraft speed increases almost linearly for the power. At some point, the beamer can no longer compen- rst half of the acceleration time. Thereafter, the speed sate for the increasing beam spillage, and the sail cools. By rises more gradually, until cuto occurs precisely as 0:2 c the end of the trajectory, less than 20% of the transmitted is reached. power reaches the sail. Essentially none of the beam is spilled until the trajec- Note that P as plotted in g. 8 does not include - tory reaches the 3 min mark. Until this time, the beamer nite photon travel time between the beamer and sailcraft, power is throttled back to prevent the sail from overheat- hence its time index is distorted. The undistorted plot ing. In the sail frame, this means that the incident and would be shifted and compressed left, then drop to zero at absorbed power, force, photon pressure and acceleration 8 min. To avoid the need to post-calculate the undistorted are all constant. As seen from the beamer frame, relativis- P , pulse energy Q is obtained by integrating (1 ) P 1 0 1 tic eects cause the apparent acceleration to pull slightly along with the trajectory. downward as speed increases. The sail's angular extent as seen from the beamer is After the 3 min mark, the beam begins to spill around initially within the pointing stability of telescopes such the sail as it moves further away. The beamer power grad- as the Hubble Space Telescope, but exceeds this as the ually ramps up to compensate for increasing losses. Evi- sail gets further away. Even at the speed of light, the max. at beamer at cutoff initial at beamer optics maximum lasers at cutoff total storage Acceleration [g’s] Beamer diameter [km] System energy efficiency [-] Capex [$] Range [au] Sailcraft mass [g] Stored energy [GWh] Relative acceleration [-] Mean radiant flux density [W/m²] Temperature [K] Power [GW] Beam duration [sec] round trip time between the beamer and sail varies from Intuitively, one might expect that, relative to the base- 0:4 s at the beginning of the trajectory to 134 s at the end. line point design in table 2, money is saved by using a Given also the high sail accelerations, it is obvious that the larger sail that is slowly accelerated by low power coming beamer cannot actively point to follow the sail. Instead, from a larger beamer over a longer period of time. This the sailcraft must be beam riding, seeking the axis of the limit corresponds to the right-hand side of each plot in beam by active or preferably passive stabilization schemes. g. 9. For sails that exceed 8 m diameter, the optimum The Doppler shift plot shows the 1:06 μm beam red- beamer diameter does indeed increase with sail diameter, shifting as seen from the sailcraft frame, eventually reach- as does beam duration and range. Consistent with expec- ing 1:30 μm wavelength. The re ected light again red- tations, acceleration decreases in this limit. However, the shifts, eventually reaching 1:60 μm wavelength. If the sail- laser power that is needed increases instead of decreases in craft were to transmit at the 0:85 μm wavelength com- the large sail diameter limit, and so capex also increases monly used by VCSEL laser diodes, it would be received due to the laser cost and the increased cost of energy stor- on Earth at 1:04 μm. age. Thus, economics does not favor the slow acceleration The sailcraft's relativistic kinetic energy reaches 6:7 TJ school of thought. (1:9 GWh). Per unit mass, this is 2 PJ kg . In compar- ison, the heat produced by Pu-238 alpha decay is three 6.1.4. Variation with respect to technology gures of merit orders of magnitude lower at 2 TJ kg . One way to tap The 0:2 c point design assumes particular cost perfor- the sailcraft's kinetic energy could be via the interstellar mance gures of merit for the laser, optics and storage. medium: The Local Interstellar Cloud is primarily com- These are referred to as k , k , and k in the simple cost l a e posed of 0:3 atoms=cm of partially ionized hydrogen [26]. model of eq. (26). But what happens if the lasers' cost From the sailcraft's perspective traveling at cruise velocity, performance, for example, is lower than expected? Also, this manifests as a monochromatic hydrogen beam that is what happens if the material absorbs more, or less than incident from the direction of travel, having a combined assumed, or the payload is heavier or lighter? kinetic energy of 55 W m , or 0:7 kW over the sail's area if it faces the direction of travel. Table 3: Relative impact of technology over/under performance In the highly unlikely event that a sailcraft were to col- Capex D P D b 1 s lide with a planetary atmosphere, then the energy released $0:1=W laser (10x) 2.4 1.5 0.46 0.82 would be equivalent to nearly 2 kt of TNT. On average, one $5000=m optics (10x) 2.4 0.52 2.5 1.4 asteroid per year enters the Earth's atmosphere with this $500=kWh storage (10x) 3.5 1.6 0.83 0.77 energy, though asteroids are orders of magnitude slower 10 absorptance (10x) 4.1 2.3 1.5 1.8 and heavier. The sailcraft would vaporize before it got 10 g payload (10x) 2.7 1.3 2.3 1.7 nearly as low into the atmosphere as asteroids do. $0:001=W laser (0.1x) 0.70 0.82 2.0 1.0 $50=m optics (0.1x) 0.48 2.0 0.45 0.74 6.1.3. Variation with respect to sail diameter $5=kWh storage (0.1x) 0.56 0.80 1.0 1.1 Sail diameter is a variable that the system model varies 10 absorptance (0.1x) 0.62 0.67 1.2 0.45 to minimize beamer capex. As shown in g. 7, it is varied 0:1 g payload (0.1x) 0.63 0.89 0.68 0.58 by the GSS algorithm of the outermost iteration (itera- 0.25 re ectance (0.36x) 3.3 1.9 2.1 1.2 tion 1). By turning o this iteration, the design space as 0.90 re ectance (1.3x) 0.76 0.87 0.84 0.96 seen by the optimizer is plotted in g. 9. Starting at the 0.99 re ectance (1.4x) 0.69 0.82 0.79 0.94 top left plot, the beamer (total) capex is minimum at the 0:1 c cruise speed (0.5x) 0.29 0.49 0.35 0.73 expected sail diameter of 4:1 m. On every plot, this diam- 0:4 c cruise speed (2x) 4.7 2.3 4.0 1.3 eter is marked by a vertical line. The line shows that the 0:9 c cruise speed (4.5x) 14 13 120 2.2 cost optimum does not exactly correspond to any other 0:99 c cruise speed (5x) 1400 44 1600 2.7 extrema; the optimum is a true tradeo between storage, 6 2 10 absorptance (10 x) 30 6.4 2.2 3.3 lasers and optics, as depicted in g. 6. 5 3 10 absorptance (10 x) 280 20 2.5 6.0 A striking feature in most of the plots in g. 9 is that the mission parameters vary in a qualitatively dierent Table 3 summarizes point designs, each of which varies way as sailcraft exceed 8 m diameter. This is because a technology gure of merit relative to the baseline. The smaller sailcraft operate in a temperature-limited regime, table lists the change in beamer capex, beamer diameter, whereas larger sailcraft operate in a beamer ux-limited beamer maximum transmit power, and sail diameter, rel- regime. On each trajectory integrator timestep, the sys- ative to the baseline of table 2. In the rst entry, the tem model calculates the sail temperature resulting from laser cost per watt is increased by 10 times. The result- maximum beamer power. If this temperature exceeds the ing capex increases only by 2.4 times because the cost- maximum, this maximum becomes the boundary condi- optimum power P halves and the cost-optimum beamer tion, and the corresponding beamer power (less than its diameter D increases (increases in transmission eciency) maximum) is calculated. to compensate. Similarly, increasing the optics cost per 11 init. at beamer max. at beamer 13 -1 4 10 10 10 -2 -3 6 -4 2 10 10 10 5 14 10 10 4 13 10 10 1,max 3 12 10 10 2 11 10 10 1,initial 1 10 10 10 0 9 10 10 -1 8 1 10 10 10 650 10 550 10 450 10 350 300 power temperature limited limited 10 250 regime regime 10 1 150 5 2 10 10 max. at sailcraft -1 -2 1 -3 0 10 10 10 0.001 0.01 0.1 0.001 0.01 0.1 0.001 0.01 0.1 Sail β at cutoff [c] Sail β at cutoff [c] Sail β at cutoff [c] Figure 10: Mission variation with cruise velocity unit area by 10 times causes the capex to increase by only subdominant and further reductions do not signi cantly 2.4 times. In this case, the cost-optimum solution halves aect the leading costs of storage and optics. For both the beamer diameter and dramatically increases the power optics and storage, reducing their costs to a tenth of the to compensate. Increasing the storage $/kWh by 10 times baseline has the eect of halving the capex. causes the capex to nearly quadruple, the biggest increase If the sail material absorbs 10 times more energy than so far because storage is the dominant cost, as seen in the expected , then the capex is more than quadruple that of top left plot of g. 9. The cost-optimum point design has the baseline, the largest increase of all. In comparison, a an increased system energy eciency, achieved by increas- tenfold increase in payload mass does not quite triple the ing the product of beamer and sail diameters, as described capex. The three order of magnitude absorptance range by the Goubau beam propagation eqs. (1) and (2). considered here is smaller than the ve order of magnitude In the happy event that a technology is even cheaper than needed, the capex is reduced less than might be in- Reducing absorptance is a good way to reduce trajectory dura- tion given that the sailcraft is mostly temperature-limited, as seen tuitively expected. Reducing laser cost to a tenth of its in g. 8. Shorter durations might be needed if the beam cannot baseline has the least eect, reducing capex to 70% of its be interrupted by passing satellites. Absorbing 0.1x increases the baseline. This is because, as seen in g. 9, laser cost is temperature-limited irradiance on the sailcraft by 10x, acceleration increases by 10x, and the duration falls from 9 min to less than 3 min. cost optimal absorbed by sail storage minimum to attain β optics total at cutoff at beamer lasers maximum at sailcraft Acceleration [g’s] Beamer diameter [m] Sail diameter [m] Capex [$] Range [au] Sailcraft mass [g] Stored energy [GWh] System energy efficiency [-] Mean radiant flux density [W/m²] Temperature [K] Power [W] Beam duration [sec] absorptance range of the current sail substrate material precursor mission point design embodies key elements of candidates [23]. For this reason, table 3 lists results for mid-21st century missions to probe the inner solar sys- higher absorptances at the end. tem through to the Oort cloud. A 0:01 c cruise velocity In all, a tenfold increase in any gure of merit does not enables the sailcraft to reach Mars in a day, Saturn in a increase the capex by nearly as much. Nor does a tenfold week, the Kuiper belt in a month, or the minimum solar decrease in a technology cost decrease the capex by nearly gravitational focus distance of 550 au in a year. A success- as much. This is because there are three constituent costs ful precursor mission proves that technologies, people and that are of comparable magnitudes after cost minimiza- processes are ready to scale up to the full 0:2 c mission. tion. The minimization process trades away large increases in one cost for smaller increases in others or all. There- 6.2.1. Inputs fore, cost minimization dampens the solution's sensitivity The inputs to the precursor mission point design are to unexpected changes in technology gures of merit. summarized in table 4. The 1:06 μm wavelength is con- Sail re ectance variations are also listed in table 3. sistent with ytterbium-doped ber ampli ers. An initial They show that if the baseline mission had used a single- sailcraft displacement of 300 km is consistent with a low layer dielectric sail with a re ectance of 25%, as shown in Earth orbit from which the sailcraft is entrained by a low- g. 5, then the capex would triple. This $18B dierence is power beam. why the baseline mission uses a photonic crystal sail with a higher re ectance of 70%. It is well known that nanostruc- Table 4: System model inputs for 0:01 c precursor mission turing the sail can improve the re ectance [23], and with 0:01 c target speed such a large nancial return, it is inconceivable that this 1:06 μm wavelength would not be done. Though almost perfect re ectors are 300 km initial sail displacement from laser source possible in principle, table 3 shows diminishing nancial returns much above 90% re ectance. Such performance 1 mg payload is attainable at a single wavelength, but the trajectory- 0:2 g m areal density averaged value will likely be lower as the Doppler shift 10 spectral normal absorptance at 1:06 μm lengthens the beam wavelength by 20% over the course 40% spectral normal re ectance at 1:06 μm of the sail acceleration. Also, the sail material's optical 625 K maximum temperature properties may be leveraged for communications, sensing, 0.01 total hemispherical emittance (2-sided, 625 K) or processing, and these competing uses would likely be worth a modest decrease in re ectance and/or increase in $1=W laser cost absorptance during sail acceleration. $10k/m optics cost $100=kWh storage cost 6.1.5. Variation with respect to cruise velocity 50% wallplug to laser eciency Cruise velocity is speci ed to be 0:2 c in the Break- 70% of beam power emerging from top of atmosphere through Starshot objectives. Nevertheless, it is interesting to vary this velocity to see how the cost-optimal mission A 1 mg payload (non-sail mass) is reserved for one or changes. The resulting plots are shown in g. 10. Each two sensors and associated support systems. Sail mass is quantity is plotted from 0:001 c to 0:99 c. The beamer (to- calculated by the system model based on the value of D tal) capex is shown in the top left plot. Also shown are chosen by the optimizer combined with the areal density constituent costs of storage, optics, and lasers. Choosing given as an input. Similar to the 0:2 c point design, this de- a mission that has half the cruise velocity, 0:1 c, decreases sign assumes a photonic crystal sail material with the same the beamer capex to a quarter of its baseline, whereas dou- thermal, mass and optical properties, except for a less am- bling it to 0:4 c quintuples it. That is to say, it costs $29B bitious re ectance of 40%. Again, the strati ed layer op- extra to halve the trip time or saves $6B to double it. In tical model is turned o, and absorptance and re ectance the case of Alpha Centauri, this shortens the 4:37 ly trip remain constant throughout trajectory integration. time from 22 years to 11 years, or lengthens it to 44 years. Cost factors are chosen to lie between present values Figure 10 also shows that cost-optimal missions are and those of the 0:2 c point design. The laser cost is $1=W, limited by sail temperature if the cruise velocity is faster two orders of magnitude lower than the present cost, con- than 0:03 c. This value of cruise velocity corresponds to sistent with requirements for laser-powered launch vehi- minimum beam duration. The minimum beamer diameter cles, yet two orders of magnitude greater than microwave reaches 10 km at 0:7 c and 100 km at greater than 0:99 c. oven magnetrons. The optics cost is $10k/m , two orders of magnitude lower than the present cost for diraction- 6.2. 0:01 c Precursor limited optics, consistent with signi cant production line Precursor missions demonstrate the key technologies automation. This value is comparable with the cost of ra- needed by Breakthrough Starshot, albeit at smaller scale, dio telescope aperture. $100=kWh energy storage is above lower speed, and lower cost than the 0:2 c missions. The current material cost oors for many energy storage tech- 13 3 nologies and in line with experience curve projections for of 9 protons=cm owing radially outward from the Sun electric vehicle battery packs circa 2035 [24]. at 0:001 c [27]. From the perspective of a sailcraft cruis- ing away from the Sun, this manifests as a 0:009 c proton 6.2.2. Results beam that is incident from the direction of travel, having a combined kinetic energy of 0:15 W m , or only 4:2 mW Upon running the system model using the inputs in over the sail's area if it faces the direction of travel. table 4, the optimizers converge to the values given in table 5. The cost-optimum sail diameter is found to be 6.3. Vacuum Tunnel 19 cm, corresponding to a mass of 6:6 mg. Forward [4] recognized that sails can be levitated and tested using medium power lasers in the 1 g gravity eld Table 5: System model outputs for 0:01 c mission of the Earth. Experimentally, this is a desirable con gu- $517M beamer capex comprised of: ration because the instrumentation stares at a stationary $285M lasers (285 MW max. transmitted power) sail within a vacuum chamber. After the sail beamrides at $224M optics (169 m primary eective diameter) 1 g, the dynamics then need to be tested and understood $8M storage (78 MWh stored energy) at ever-increasing accelerations. This can be accomplished horizontally in a vacuum tunnel. Gravity will try to pull $8k energy cost per mission (78 MWh @$0:1=kWh) the sail from the axis of the beam, and the sail should 0.01% system energy eciency counteract this and other perturbations. A key milestone is when the acceleration achieved by the sail equals the 19 cm sail diameter 14 900 g s acceleration needed by the 0:2 c mission. After 6:6 mg sailcraft mass (includes payload mass) that, the most visible milestone is reached when the sail becomes the fastest human-made craft. The current record 6 min (346 s) beam transmit duration is held by the Helios-B probe, which achieved a heliocen- 6 min (349 s) sail acceleration duration tric speed of 70:22 km s during its closest pass of the Sun in April 1976. It is expected that this record will be 23 Pa temperature-limited photon pressure broken by the 200 km s speeds of the upcoming Parker 0:64 N temperature-limited force Solar Probe. But is it practical for Breakthrough Starshot 10 000 g s temperature-limited acceleration to attempt the speed milestone in a ground-based vacuum 7 g s nal acceleration (0:007 au), 3:3 ls from source tunnel? 13 kW m beamer maximum beam radiant exitance 2 6.3.1. Inputs 8:7 GW m sailcraft temperature-limited irradiance At 0:01 c, the sailcraft has a relativistic kinetic energy Table 6: System model inputs for the vacuum tunnel of 8 kWh, whereas the beamer uses 78 MWh. This yields 1:06 μm wavelength 0.01% system energy eciency; two orders of magnitude 1 m initial sail displacement from laser source lower than the 0:2 c mission. However, this energy costs only $8k at a price of $0:1=kWh, making it ve orders of 10 ng payload magnitude cheaper than the beamer capex of $517M. 0:25 g m areal density At 169 m eective diameter, the 0:01 c beamer is 16 10 spectral normal absorptance at 1:06 μm times smaller than the 0:2 c beamer. Also, it has one 35% spectral normal re ectance at 1:06 μm third the maximum radiant exitance at 13 kW m vs. 625 K maximum temperature 37 kW m . The temperature-limited irradiance in the 0.01 total hemispherical emittance (2-sided, 625 K) sailcraft frame is the same as that of the 0:2 c mission because the sailcraft's absorptance, emittance, and tem- $100=W laser cost perature limit are the same. The initial sailcraft accelera- $1M/m optics cost tion is two thirds that of the 0:2 c mission because of the $500=kWh storage cost lower re ectance, but the cost-optimal trajectory ends at $10k/m vacuum tunnel cost only 7 g s, squeezing almost everything it can from the di- 50% wallplug to laser eciency minishing beam. Even then, the energy storage costs only $8M because storage is cheap relative to lasers and optics, The inputs to the vacuum tunnel system model are and because there is a limit to the extent that laser and detailed in table 6. Relative to the inputs for the 0:01 c optics costs can be reduced by increasing the pulse length. mission, the vacuum tunnel has a much shorter initial sail The sailcraft's relativistic kinetic energy reaches 30 MJ displacement of 1 m, consistent with a sail positioned close (8 kWh). Per unit mass, this is 4 TJ kg . In comparison, to an optic at the start of the tunnel. Also, the payload the heat produced by Pu-238 alpha decay is half as much at is reduced to a token 10 ng, the mass of a few cells, and 2 TJ kg . At 1 AU, the solar wind is primarily composed the re ectance is reduced to a near-term value of 35% per 14 10 10 1 -5 10 0.1 -6 10 0.01 4 7 10 10 10 1 10 10 4 1,max -1 -2 -3 2 0.1 10 10 -5 10 10 650 -6 -7 0.1 -8 10 400 5 6 12 10 10 10 maximum maximum at sail at cutoff 4 8 10 10 maximum at beamer 3 1 4 10 10 10 1 10 100 1 10 100 1 10 100 Sail speed at cutoff [km/s] Sail speed at cutoff [km/s] Sail speed at cutoff [km/s] Figure 11: Tunnel characteristics vs. sail speed at cuto the strati ed layer re ectance predictions in g. 5. Such a lower, then the re ectance would need to be increased to re ectance is consistent with commercially-available Si N obtain sucient acceleration. 3 4 membrane x-ray windows, for example. Again, the strati- The cost factors are chosen to be consistent with cur- ed layer optical model is turned o, and the absorptance rent market values. In addition to laser, optics, and stor- and re ectance remain constant throughout trajectory in- age costs, there is an additional vacuum tunnel cost of tegration. $10k/m. This tunnel cost factor is consistent with the To equal the acceleration of the 0:2 c mission, the sail cost of the LIGO vacuum tunnels [28] and includes the ideally has the same areal density, re ectance, absorp- beam tunnel and its enclosure, as well as vacuum equip- tance, and temperature limit as the 0:2 c sail. For the ment. To account for dierences from the 1 m LIGO tunnel vacuum tunnel, point designs assume higher areal density diameter, the tunnel cost is multiplied by the relative tun- and lower re ectance to be consistent with the idea that nel diameter. The tunnel cost does not include items that for the time being, a sail is a 1 mm diameter high-purity do not scale with length such as R&D, detectors, project lm that is not a photonic crystal. Correspondingly, the management, laboratory construction, and operations. absorptance is assumed to be an order of magnitude lower than that of the 0:2 c sail. If the absorptance were not optics cost optimal minimum to attain β tunnel total lasers Acceleration [g’s] Beamer diameter [m] Sail diameter [mm] Capex [$] Pipeline length [m] Sailcraft mass [g] Stored energy [Wh] System energy efficiency [-] Mean radiant flux density [W/m²] Temperature [K] Power [W] Beam duration [sec] 6.3.2. Results precursor mission. Also, a family of solutions is computed for a ground-based vacuum tunnel in which beam-riding Upon running the system model using the inputs in ta- and other aspects of the sail can be tested. All assume ble 6, the model converges to the values plotted in g. 11. the case of a circular dielectric sail that is accelerated by This gure is comprised of a family of point designs span- photon pressure from a 1:06 μm wavelength beam. Earlier ning the range of 1-500 km s sail speed at cuto. Slower investigators were led astray by pursuing high re ectance than 20 km s , laser cost is of primary importance and alone, leading them to metallic sails and/or heavy mul- tunnel cost is secondary. Faster than 30 km s , tunnel tilayer dielectric sails. This paper shows that re ectance cost is of primary importance and laser cost is secondary. divided by sailcraft areal density (or thickness, if the ma- In all cases, optics cost is tertiary. Storage cost turns out terial is held constant) is a gure of merit that leads to to be trivial, not even reaching $1000 at 500 km s , so the substantially improved system performance relative to ear- cost plot is not scaled to show it. lier sailcraft concepts. A strati ed layer model shows that If the sail had no payload, it would have a temperature- the sail acceleration per unit power is maximized by a limited acceleration exceeding 80 000 g s regardless of tun- single-layer dielectric that is =7 thick. An ideal technol- nel length; this acceleration is more than ve times that ogy would improve re ectance by removing, as opposed to of the 0:2 c mission. When payload mass is taken into ac- adding, mass. For this reason, two-dimensional nanohole count, the temperature-limited acceleration exceeds that photonic crystals are of interest as a future sail material. of the 0:2 c mission only in tunnels longer than 50 m. The 0:2 c point design minimizes capital cost by accel- For 20 km s sail speed at cuto, sucient to escape erating a 4:1 m diameter sailcraft for 9 min. In minimizing the solar system starting from Earth, the model infers: the cost, it is surprising that laser costs are secondary to A 0:4 km tunnel, 22 kW of lasers, and a 0:6 m diameter optics and storage costs, because this implies that there is telescope, costing a total of $5M. a limit to the amount that cheap laser power can compen- For 200 km s sail speed at cuto, which equals or sate for expensive optics or storage. The point design as- exceeds the fastest human-made craft, the model infers: sumes $0:01=W lasers and $500=m optics to achieve $8.0B A 28 km tunnel, 4:6 MW of lasers, and a 3:1 m diameter capital cost for the ground-based beamer. In contrast, the telescope, costing a total of $1.3B. energy needed to accelerate each sailcraft is a thousand There is a small but noticeable bump toward the upper times cheaper, making the 2.9% system energy eciency right corner of the sail diameter and other plots in g. 11, unimportant. With large xed costs and low incremental near the 50 km s mark. This bump is accompanied by costs, why not use the beamer to propel sailcraft to every a subtle change in gradient of sail diameter vs. speed. reachable star as often as possible? Therefore, a rational Decreasing the integration step size and convergence tol- outcome of Breakthrough Starshot might be to pave the erances does not remove the bump, nor does the model hit way for multi-lightyear pipelines of sailcraft that y past any kind of assumed limit on an internal variable. Also, each target star every few weeks ad in nitum. For the plotting costs vs. sail diameter for the speed at which the 0:2 c mission, it is clear that the sailcraft must be beam bump occurs shows that the minimum (found by iteration riding; the speed of light is too slow to allow the beam 1 of the solution procedure shown in g. 7) is a global min- to follow the sailcraft. If the beamer slews the beam dur- imum and that there are no competing local minima that ing acceleration, it will be to dodge satellites or ne-tune could cause the solution to jump from one value to another. the sailcraft's destination. It is also clear that nuclear It turns out that the source of this bump is the Goubau batteries are dead weight, having a thousand times lower beam model where the two eciencies described by eq. (2) speci c energy than the sailcraft's kinetic energy. From are spliced together. Though this splice preserves the con- the sailcraft's perspective cruising at 0:2 c, the interstellar tinuity of the function, it introduces a discontinuity in its medium manifests as a 0:7 kW monoenergetic hydrogen gradient. Smoothing this function in the vicinity of the beam that is incident from the direction of travel. A key splice has the eect of removing the bump in g. 11. question for future research is, what fraction of this power can be harvested? 7. Conclusions For an extra $29B, the cruise velocity can be doubled to 0:4 c, which halves the cruise time to Alpha Centauri to In this paper, a system model is formulated to describe 11 years. Following a successful 0:2 c mission, there would a beam-driven sailcraft. It minimizes beamer capital cost be high con dence in incrementally upgrading the beamer by trading o the relative expenses of lasers, optics, and to support 0:4 c missions. Hence, it is not unreasonable to storage. The system model employs nested numerical op- expect that private and government investments over sev- timizers and trajectory integration, whereas earlier models eral years could amount to $29B or more. However, with were based on closed-form approximations. The outcome increasing cruise velocity comes increasing beamer diame- is that the solution is cheaper and generates more accurate ter. There is a practical and desired limit to beamer diam- requirements, but it also exhibits more complex behaviors. eter, but how large is it? Cities are perhaps the best guide The system model is used to compute point designs for because they are human-engineered surfaces of the largest a 0:2 c Alpha Centauri mission and a 0:01 c solar system diameter. Greater London has a city area of 1572 km . If 16 it were circular, it would have a diameter of 45 km. Fig- assume less favorable values for laser and optics costs, and ure 10 predicts that a beamer the size of London would be the beam tunnel calculations assume present-day values. capable of propelling a sailcraft to greater than 0:9 c. Based in part on requirements derived from the system The 0:01 c point design diers from the 0:2 c point de- model, the array elements that comprise the beamer will sign in that it assumes nearer-term cost factors of $1=W become better de ned in future. Once designs are avail- and $10k/m to achieve $517M capital cost for its beamer able for the beamer array elements, it makes sense to up- and $8k energy cost per 19 cm diameter sailcraft. This date the cost inputs for the 0:2 c mission based on calcu- sailcraft is 20 times smaller and 600 times lighter than the lated laser and optics material cost oors combined with 0:2 c sailcraft, so its various subsystems need to be inte- an analysis of the array element production chain. grated into a much smaller mass and area. But, for such a The results presented here are based on uncertain input small photonic crystal sailcraft, what functionality is theo- values, and table 3 shows that the cost-optimal solutions retically possible and what sail area will it take? The 0:2 c are surprisingly resilient to them. Of the inputs, laser and mission will be bombarded by dust and radiation over two optics costs are the most uncertain. Their orders of mag- decades, so it may need many duplicates of each subsystem nitude uncertainty can be reduced by studies of materials to reach its mission objective. Hence, the 0:01 c sailcraft cost oors and production automation. The sail material could be developed as a single unit cell of the 0:2 c sail- is also an important source of uncertainty, in particular craft. If low mass and low energy cost translate into low its absorptance, and to lesser extents its emittance, re- incremental costs, then it makes sense to use many pre- ectance, areal density and maximum temperature. Can- cursors to prove the sailcraft technologies. Precursors can didate sail materials are wide band-gap dielectrics because probe magnetic elds and dust particle uxes ahead of the their absorptances are very low and rise the least as tem- 0:2 c missions. Indeed, the whole heliosphere is a testbed perature rises. This dependence of absorptance on tem- for sailcraft communication and sensor technologies. At perature is important and needs to be captured in future 0:01 c, the Kuiper belt is only a month away. Unlike the versions of the system model. Because of the large num- 0:2 c sailcraft, the 0:01 c sailcraft has very little power: Its ber of sail design parameters to explore, near-term work speci c kinetic energy is comparable with the speci c heat should de ne perhaps three promising reference materi- of Pu-238 alpha decay. From the perspective of a sailcraft als and characterize their temperature-dependent proper- cruising at 0:01 c, the solar wind manifests as only a 4 mW ties and other properties as needed by the system model. proton beam incident from the direction of travel. In com- Given current uncertainties, table 3 shows that the sail- parison, the solar power available at Kuiper belt distances craft material properties make more than a billion-dollar of 30 au is an order of magnitude higher. With solar power dierence to the beamer cost. Hence, it is worth consider- only, inner solar system missions are of course easier. able R&D investment to characterize and optimize mate- The ground-based vacuum tunnel assumes present-day rials speci cally for laser-driven sailing. cost factors of $100=W lasers, $1M/m optics, and $10k/m Finally, the modeling work presented here takes place vacuum tunnel for its family of point designs spanning the within the context of a wider systems engineering eort. range of 1-500 km s maximum sail speed. Primarily, a In future, system model development will be driven by the tunnel would exist to test sail beam-riding dynamics. Sec- increasing need to infer requirements and margins from an ondarily, the tunnel could be used to demonstrate high up-to-date and self-consistent understanding of the prob- sail velocities, and doing so would prove mastery of beam- lem. The system model can be extended to handle prob- riding dynamics transverse to the optical axis and also lon- abilistic inputs in order to estimate parameter sensitivi- gitudinally as the beam focus accelerates. Under the ma- ties, and to bound requirements and margins. It can also terial properties assumed here, speed milestones exceed- be extended to incorporate domain models. Such mod- ing 100 km s might be achieved at lowest cost by using els usually add computational complexity, so they need to space-based missions as opposed to ground-based tunnels. be blended with the existing simplistic models to converge However, tunnels may still be cheapest if the incremental system point designs in an acceptable time frame. Domain cost of upgrading an existing facility is low enough. Also, models would add greater detail and accuracy in the areas tunnel length and cost depend on sail material properties, of beamer elements and phasing; transatmospheric beam so this conclusion should be revisited as sail materials be- propagation; communications; cost modeling; and sailcraft come better characterized. thermal, structural, and electromagnetic properties. A key thesis of Starshot is that the cost of lasers and optics can and will fall. For 0:2 c missions, the laser cost 8. Acknowledgments per unit power and optics cost per unit area are chosen such that the beamer costs less than $10B. Consequently, This work was supported by the Breakthrough Prize $0:01=W laser cost and $500=m optics cost should be un- Foundation. I acknowledge helpful discussions with Pete derstood as requirements that derive from the $10B cost Klupar, Bruce Draine, James Benford, Ognjen Ilic, and cap. These requirements are to be achieved through tech- Harry Atwater. nology development and production line automation. Un- constrained by a cost cap, nearer-term precursor missions 17 References [24] O. Schmidt, A. Hawkes, A. Gambhir, I. Staell, The future cost of electrical energy storage based on experience rates, Nature Energy 2 (8) (2017) 17110. References [25] B. Stuart, M. Feit, A. Rubenchik, B. Shore, M. Perry, Laser- induced damage in dielectrics with nanosecond to subpicosecond [1] T. H. Maiman, Stimulated optical radiation in ruby, nature pulses, Physical Review Letters 74 (12) (1995) 2248. 187 (4736) (1960) 493{494. [26] P. C. Frisch, S. Red eld, J. D. Slavin, The interstellar medium [2] R. Forward, Pluto the gateway to the stars. missiles and rockets surrounding the sun, Annual Review of Astronomy and Astro- 10, 26{28; reprinted as pluto: Last stop before the stars, Science physics 49 (2011) 237{279. Digest 52 (1962) 70{75. [27] R. Schwenn, Solar wind: Global properties, Encyclopedia of As- [3] R. L. Forward, Roundtrip interstellar travel using laser-pushed tronomy and Astrophysics (Nature Publishing Group, Bristol, lightsails, Journal of Spacecraft and Rockets 21 (2) (1984) 187{ 2001). [28] P. Lindquist, LIGO nancial status, Tech. Rep. LIGO-G020261- [4] R. L. Forward, Laser weapon target practice with gee-whiz 00-P, California Institute of Technology, presented at The targets, in: Proceedings SDIO/DARPA Workshop on Laser Twelfth Program Advisory Committee (PAC12), Massachusetts Propulsion, Vol. 2, 1986, pp. 41{44. Institute of Technology, June 27 -28, 2002 (2002). [5] G. A. Landis, Optics and materials considerations for a laser- propelled lightsail, paper IAF-89-664 presented at the 40th In- ternational Astronautical Congress, Malaga, Spain, Oct. 7-13, [6] J. Guillochon, A. Loeb, SETI via leakage from light sails in ex- oplanetary systems, The Astrophysical Journal Letters 811 (2) (2015) L20. [7] P. Lubin, A roadmap to interstellar ight, Journal of the British Interplanetary Society 69 (2016) 40{72. [8] G. Goubau, F. Schwering, On the guided propagation of elec- tromagnetic wave beams, IRE Transactions on Antennas and Propagation 9 (3) (1961) 248{256. [9] G. Goubau, F. Schwering, 3.5 free space beam transmission, Microwave Power Engineering: Generation, transmission, recti- cation 1 (1968) 241. [10] R. Hansen, J. McSpadden, J. N. Benford, A universal power transfer curve, IEEE Microwave and Wireless Components Let- ters 15 (5) (2005) 369{371. [11] W. C. Brown, E. E. Eves, Beamed microwave power transmis- sion and its application to space, IEEE Transactions on Mi- crowave Theory and Techniques 40 (6) (1992) 1239{1250. [12] K. L. G. Parkin, T. Lambot, Microwave thermal propulsion nal report, Technical Publication NASA/TP{2017{219555, NASA (2015). [13] N. Kulkarni, P. M. Lubin, Q. Zhang, Relativistic solutions to directed energy, in: Proc. SPIE 9981, Planetary Defense and Space Environment Applications, 2016. [14] A. Einstein, On electrodynamics of moving bodies, Annals of Physics 322 (10) (1905) 891{921. [15] H. Robertson, Dynamical eects of radiation in the solar system, Monthly Notices of the Royal Astronomical Society 97 (1937) [16] N. Kulkarni, P. Lubin, Q. Zhang, Relativistic spacecraft pro- pelled by directed energy, The Astronomical Journal 155 (4) (2018) 155. [17] H. A. Macleod, Thin- lm optical lters, 4th Edition, CRC press, [18] J. Benford, Starship sails propelled by cost-optimized directed energy, Journal of the British Interplanetary Society 66 (2013) [19] J. T. Kare, K. L. G. Parkin, A comparison of laser and mi- crowave approaches to CW beamed energy launch, in: AIP Conference Proceedings, Vol. 830, AIP, 2006, pp. 388{399. [20] R. L. Forward, Light-levitated geostationary cylindrical orbits using perforated light sails, Journal of the Astronautical Sci- ences 32 (2) (1984) 221{226. [21] S. Baig, C. R. McInnes, Light-levitated geostationary cylindri- cal orbits are feasible, Journal of Guidance, Control, and Dy- namics 33 (3) (2010) 782{793. [22] M. Wandel, Attenuation in silica-based optical bers, Ph.D. the- sis, Department of Communications, Optics & Materials, Tech- nical University of Denmark (2005). [23] H. A. Atwater, A. R. Davoyan, O. Ilic, D. Jariwala, M. C. Sher- rott, C. M. Went, W. S. Whitney, J. Wong, Materials challenges for the starshot lightsail, Nature Materials (2018) 1.
Astrophysics – arXiv (Cornell University)
Published: Apr 24, 2018
You can share this free article with as many people as you like with the url below! We hope you enjoy this feature!
Read and print from thousands of top scholarly journals.
Already have an account? Log in
Bookmark this article. You can see your Bookmarks on your DeepDyve Library.
To save an article, log in first, or sign up for a DeepDyve account if you don’t already have one.
Copy and paste the desired citation format or use the link below to download a file formatted for EndNote
Access the full text.
Sign up today, get DeepDyve free for 14 days.
All DeepDyve websites use cookies to improve your online experience. They were placed on your computer when you launched this website. You can change your cookie settings through your browser.