Alien suns reversing in exoplanet skies

Earth’s rapid spin, modest tilt, and nearly circular orbit ensure that the sun always appears to move forward, rising in the east and setting in the west. However, for some exoplanets, solar motion can reverse causing alien suns to apparently move backward. Indeed, this dramatic motion marginally occurs for Mercury in our own solar system. For exoplanetary observers, we study the scope of solar motion as a function of eccentricity, spin–orbit ratio, obliquity, and nodal longitude, and we visualize the motion in spatial and spacetime plots. For zero obliquity, reversals occur when a planet’s spin angular speed is between its maximum and minimum orbital angular speeds, and we derive exact nonlinear equations for eccentricity and spin–orbit to bound reversing and non-reversing motion. We generalize the notion of solar day to gracefully handle the most common reversals.

The recent surge in exoplanet discoveries expands our inventory of solar systems even as it challenges our understanding of them. Many exoplanets have unusual orbits, including those around multiple star systems 1,2 . Such bizarre solar dances imply exotic sunrises and sunsets, as when a planet orbits one star that orbits another, so that when the planet is opposite the two stars it experiences two sunsets and sunrises per day, but when it is between the two stars it experiences constant daylight.
In this article, we investigate the range of apparent solar motion for planets orbiting single stars in simple models using theory and computation. In particular, we explore the ranges of orbital eccentricity, spin-orbit ratio, and axial tilt that permit apparent solar reversals, an especially striking motion: not only can suns stand still in the sky, they can sometimes move backward. Planets in plausible pseudosynchronized rotations or spin-orbit resonances should experience solar reversals. Although the dynamics of most exoplanets are not yet fully known, we begin locating some exoplanets in the abstract parameter space of solar reversals.
We first focus on the important special case of zero tilt and derive exact nonlinear equations delimiting the range of eccentricity and spin-orbit ratio that allows solar reversals. We note that Mercury is just inside this range. Spacetime "coil" plots illustrate important examples. We next extend these results numerically to planetary tilt, parameterized by obliquity and nodal angle. Spatial "ball" plots illustrate representative cases. We then suggest a reversal-aware generalization of a solar day and consider implications of reversals for civilizations that witness them. Finally, we qualitatively discuss chaotic and multiple star systems and locate the current work in a broader context.

Without tilt
A planet of mass m orbits a fixed star of mass M in an ellipse of radii r − ≤ r ≤ r + and speeds v + ≤ v ≤ v − , where subscripts ± indicate periapsis and apoapsis values. Its constant angular momentum magnitude and energy The difference in energies between the apopapsis and periapsis is Eliminate the speeds with the corresponding angular momenta to get Assuming constant angular spin ω s and zero obliquity, apparent solar reversals occur when the spin angular speed is between the minimum and maximum orbital angular speeds, or where the spin-orbit ratio ρ = ω s /ω o . The dimensionless function quantifies the reversal, where R < 0 means a sun always appears to move one-way in the sky, and R > 0 means its movement sometimes appears to reverse. Figure 1 summarizes the possibilities. Solid black curves are plots of orbital eccentricity e (ordinate) versus critically large and small spin-orbit ratios ρ ± = ω ± /ω o (abscissa). Colors code reversals R , with red negative, blue positive, and saturation decreasing with magnitude. Dots indicate Venus and its unusual retrograde spin 3,4 , Mercury and its spectacular 3-2 resonance 5 , and β Pictoris b 6 , the first exoplanet with a measured spin-orbit ratio. Mercury is just inside the blue reversal region. Other solar system planets would appear near β Pictoris b. Horizontal dashed line indicates isolated Kepler planets with typical eccentricities e ∼ 0.3 , which experience reversals for spin-orbit ratio ρ ∼ 1 (while Kepler multiples tend to have nearly circular orbits) 7,8 . Dot-dashed curve indicates pseudosynchronized planets 9 , such as the eccentric exoplanet HD 80606 b, all of which would experience reversals.
We visualize the apparent solar motion using Mathematica 10 to numerically integrate the nonlinear polar equations of motion starting at the apoapsis r = r + , θ = 0 , ṙ = 0 , θ = ω + and compute the angle www.nature.com/scientificreports/ and distance r of the sun with respect to an equatorial observer (neglecting the planet's radius), as in the Fig. 2 spacetime "coil" plots. A 2 + 1 dimensional spacetime extrusion (red) represents the observer on the local horizon (green). A tube represents the apparent solar motion, both forward (yellow) and backward (cyan). Inset (topright) represents orbital motion in the solar system's reference frame, where spin (red arrow) is perpendicular to the orbits and dot (green) is the sun. Figure 2a illustrates an Earth-like planet whose sun apparently rises in the east and sets in the west with no reversals; parameters are {ρ, e} = {2, 0.1} . Figure 2b illustrates a non-Earth-like planet with a more eccentric orbit whose sun rises in the east and sets in the west but reverses at every noon and every perihelion (although equatorial observer at the antipodes would not see these solar reversals); parameters are {ρ, e} = {2, 0.5}.  www.nature.com/scientificreports/ We demonstrate the range of possibilities by drawing coil plots perpendicular to a grid of spin-orbit ratio and eccentricity {ρ, e} , as in Fig. 3. Red-white-blue background colors again code the reversal R . Red extrusions indicate observers, with time increasing outward for 8 orbital periods. For tidally-locked planets in circular orbits, {ρ, e} = {1, 0} and their suns appear to hang motionless in their skies; as eccentricity increases and the orbital speed varies, their suns appear to move backward and forward, even rising and setting in the east followed by rising and setting in the west near {ρ, e} = {1, 0.9} . For a given orbital eccentricity e, only a range of spin-orbit ratios between ρ ± = ω ± /ω o produce reversals (cyan and yellow).
Mercury is just inside the reversal region, and Fig. 4 focuses on its coil plot (neglecting perturbations due to asphericity, other planets, general relativity, and so on). On Mercury, one (solar) day lasts 2 years, and once a year an equatorial observer witnesses a brief solar reversal, surely a special day for any future inhabitants.

With tilt
More generally, consider a planet whose spin is tilted by an obliquity angle ǫ with respect to the orbital angular momentum vector � L ∝ẑ and whose equatorial and orbital planes intersect at a nodal angle with respect to the orbital eccentricity vector (or long axis of the elliptical orbit) � ǫ ∝x . At time t, an equatorial observer has a local north-up-east coordinate system where R φû is a rotation matrix through an angle φ about a direction û . More explicitly, the local unit vectors are The 3D Fig. 5a plot shows the region of reversals as a function of orbital eccentricity ǫ , spin-orbit ratio ρ = ω s /ω o , and obliquity ǫ , for nodal angle = 0 . For planets or satellites that maintain significant obliquities, dynamical evolution and forced nodal precession can lead to capture into long-term or secular spin-orbit resonance known as Cassini states 11 . Such planets maintain obliquity despite tidal dissipation, and the equilibrium rotation rate is a function of eccentricity and obliquity 12 , as indicated by the Fig. 5a checkered surface, all of which undergo reversals. The 2D Fig. 5b plot shows that small obliquities do not significantly change the reversals, where the dashed line represents Earth's ǫ = 23.44 • ∼ 0.4 obliquity, although large obliquities expand the range. Dot-dashed curve indicates Cassini states. Similarly, small changes in the nodal angle have a small effect.
We visualize the apparent solar motion by the Fig. 6 spatial "ball" plots. A dot (red) represents the observer on the local horizon (green). A tube represents the apparent solar motion, both forward (yellow) and backward (cyan). Inset (top-right) represents orbital motion in the solar system's reference frame, where spin (red arrow) is tilted and dot (green) is the sun. Figure 6a illustrates an Earth-like planet whose sun apparently rises in the east and sets in the west with altitude varying with season and with no reversals; parameters are {ρ, e, ǫ, �} = {8, 0.1, 0.4, 3.8} . Figure 6b illustrates a non-Earth-like planet in a more eccentric orbit with a more oblique tilt whose sun rises in the east and sets in the west but reverses twice a year and not near periapsis; parameters are {8, 0.6, 1.5, 3.8} . The Fig. 7 nested plots of {ǫ, �} at each {ρ, e} summarize the scope of solar motion. Again, moderate tilt is secondary to orbital eccentricity and spin-orbit ratio in determining the passage of the sun in the sky, including its reversals.
Nonequatorial observers can also experience solar reversals. Figure 8a simulates the view from an Earth-like planet whose sun appears to rise in the east and set in the west every day for a year; rainbow hues code time and the sun is drawn at equal time intervals. Figure 8b simulates a Mercury-like exoplanet, nicknamed planet "Z", whose sun appears to reverse its motion just once a year at noon and at periapsis where the sun appears (17a) N = {sin ǫ sin �, sin ǫ sin �, cos ǫ}, (17b) U = {cos � sin ω s t + sin � cos ǫ cos ω s t, sin � sin ω s t − cos � cos ǫ cos ω s t, − sin ǫ cos ω s t}, (17c) E = {cos � cos ω s t − sin � cos ǫ sin ω s t, sin � cos ω s t + cos � cos ǫ sin ω s t, sin ǫ sin ω s t}.

Apoday
For small obliquity, which may be common and includes the exoplanet β Pictoris b 13 , and where the planar nature of the geometry even admits convenient analysis using complex numbers 14 , the without-tilt results suffice and reversals occur only near periapsis, as in the Fig. 8b strobe plot. A sun might even appear to rise in the east and set in the west, and then rise in the west for a "reverse" or "anti" day before setting in the east and repeating.  www.nature.com/scientificreports/ For planets in eccentric orbits with small obliquity and constant spin angular speed between their extreme orbital angular speeds, many "normal" days, when the planet orbits slowly near its apopapsis, will be punctuated once a year by a reversal day, when the planet orbits quickly near its periapsis, and the sun appears to crisscross the local meridian three times, as in Fig. 8b. For such cases, the notion of a solar day as the time between successive noons might need to be revised or discarded. We suggest defining the (variable) apoday to be the time between successive solar meridian crossings in the dominant, apoapsical direction. We expect solar reversals to have enormous cultural implications for any intelligent life that witnesses them.

Discussion
The motion of alien suns in exoplanet skies, long the subject of science fiction 15 , is now becoming a focus of science fact. For example, observations 16 by the recently launched James Webb Space Telescope of the periapsis heating of the high-eccentricity exoplanet HD 80606 b should elucidate its potential for both pseudosynchronization and solar reversals.
We do not yet know what fraction of exoplanets experience solar reversals and for how long. In our solar system, only Mercury experiences reversals, and just barely. The striking spin-orbit resonance and large eccentricity that make this possible may reflect a history (and future) of chaotic dynamics 5,17 . Also in our solar system, observers on Hyperion see Saturn and the sun rise and set chaotically 18 , while observers on Nix and Hydra see Pluto and the sun rise and set chaotically 19 . Similarly, some exoplanets may experience chaotic sunsets, at least for some of their lifetimes (as eccentricities, spins, and tilts can wander and evolve).
Exoplanets may be stable in multiple star systems 20 , where multiple reversals may be common. The recently discovered sextuply-eclipsing sextuple star system 21 suggests dramatic possibilities. Like the growing interest in planetary habitability 22 and alien skies 23 , knowledge of alien sunsets, sunrises, and solar motion enlarges our conception of possible worlds. Study of potential and other solar systems provides context for study of our own.

Data availability
A Mathematica notebook that generates most of the results in this article is available from the corresponding author on reasonable request.