In this paper we study the morphological properties of spiral galaxies, including measurements of spiral arm number and pitch angle. Using Galaxy Zoo 2, a stellar mass-complete sample of 6,222 SDSS spiral galaxies is selected. We use the machine vision algorithm SpArcFiRe to identify spiral arm features and measure their associated geometries. A support vector machine classifier is employed to identify reliable spiral features, with which we are able to estimate pitch angles for half of our sample. We use these machine measurements to calibrate visual estimates of arm tightness, and hence estimate pitch angles for our entire sample. The properties of spiral arms are compared with respect to various galaxy properties. The star formation properties of galaxies vary significantly with arm number, but not pitch angle. We find that galaxies hosting strong bars have spiral arms substantially ($4-6^\mathrm{o}$) looser than unbarred galaxies. Accounting for this, spiral arms associated with many-arm structures are looser (by 2$^\mathrm{o}$) than those in two-arm galaxies. In contrast to this average trend, galaxies with greater bulge-to-total stellar mass ratios display both fewer and looser spiral arms. This effect is primarily driven by the galaxy disc, such that galaxies with more massive discs contain more spiral arms with tighter pitch angles. This implies that galaxy central mass concentration is not the dominant cause of pitch angle and arm number variations between galaxies, which in turn suggests that not all spiral arms are governed by classical density waves or modal theories.

The Galaxy Zoo has provided morphological data on many galaxies. Several biases have been identified in the Galaxy Zoo data. Here we report on a newly discovered selection effect: astronomers interested in studying spiral galaxies may select a set of spiral galaxies based upon a threshold in spirality (the fraction of Galaxy Zoo humans who report seeing spiral structure). SpArcFiRe is an automated tool that decomposes a spiral galaxy into its constituent spiral arms, providing objective, quantitative data on their structure. SpArcFiRe measures the pitch angle of spiral arms. We have observed that when selecting a set of spiral galaxies based on a threshold on spirality, the pitch angle of spiral arms appear increase with redshift. We hypothesize that this is a selection effect: tightly-wound spiral arms become less visible as images degrade with increasing redshift, leading to fewer such galaxies being included in the sample at higher redshifts. We corroborate this hypothesis by artificially degrading images of nearby galaxies, then using a machine learning algorithm trained on Galaxy Zoo data to provide a spirality for each artificially degraded image. It correctly predicts that spirality decreases as image quality degrades. Thus, the mean pitch angle of those galaxies remaining above the spirality threshold is higher than those eliminated by the selection effect. This demonstrates that users who select samples of galaxies using a threshold of Galaxy Zoo votes must carefully consider the possibility of selection effects on morphological measures, even if the measure itself is believed to be objective and unbiased. Finally, we also perform an empirical sensitivity analysis to demonstrate that SpArcFiRe's output changes in a smooth and predictable fashion to changes in its internal algorithmic parameters.

The Galaxy Zoo 1 catalog displays a bias towards the S-wise winding direction in spiral galaxies which has yet to be explained. The lack of an explanation confounds our attempts to verify the Cosmological Principle, and has spurred some debate as to whether a bias exists in the real universe. The bias manifests not only in the obvious case of trying to decide if the universe as a whole has a winding bias, but also in the more insidious case of selecting which galaxies to include in a winding direction survey. While the former bias has been accounted for in a previous image-mirroring study, the latter has not. Furthermore, the bias has never been {\em corrected} in the GZ1 catalog, as only a small sample of the GZ1 catalog was re-examined during the mirror study. We show that the existing bias is a human {\em selection} effect rather than a human chirality bias. In effect, the excess S-wise votes are spuriously "stolen" from the elliptical and edge-on-disk categories, not the Z-wise category. Thus, when selecting a set of spiral galaxies by imposing a threshold $T$ so that $\max(P_S,P_Z) > T$ or $P_S+P_Z>T$, we spuriously select more S-wise than Z-wise galaxies. We show that when a provably unbiased machine selects which galaxies are spirals independent of their chirality, the S-wise surplus vanishes, even if humans are still used to determine the chirality. Thus, when viewed across the entire GZ1 sample (and by implication, the Sloan catalog), the winding direction of arms in spiral galaxies as viewed from Earth is consistent with the flip of a fair coin.

Given an approximately centered image of a spiral galaxy, we describe an entirely automated method that finds, centers, and sizes the galaxy and then automatically extracts structural information about the spiral arms. For each arm segment found, we list the pixels in that segment and perform a least-squares fit of a logarithmic spiral arc to the pixels in the segment. The algorithm takes about 1 minute per galaxy, and can easily be scaled using parallelism. We have run it on all ~644,000 Sloan objects classified as "galaxy" and large enough to observe some structure. Our algorithm is stable in the sense that the statistics across a large sample of galaxies vary smoothly based on algorithmic parameters, although results for individual galaxies can sometimes vary in a non-smooth but easily understood manner. We find a very good correlation between our quantitative description of spiral structure and the qualitative description provided by humans via Galaxy Zoo. In addition, we find that pitch angle often varies significantly segment-to-segment in a single spiral galaxy, making it difficult to define "the" pitch angle for a single galaxy. Finally, we point out how complex arm structure (even of long arms) can lead to ambiguity in defining what an "arm" is, leading us to prefer the term "arm segments".

We describe a system for the automatic quantification of structure in spiral galaxies. This enables translation of sky survey images into data needed to help address fundamental astrophysical questions such as the origin of spiral structure---a phenomenon that has eluded theoretical description despite 150 years of study (Sellwood 2010). The difficulty of automated measurement is underscored by the fact that, to date, only manual efforts (such as the citizen science project Galaxy Zoo) have been able to extract information about large samples of spiral galaxies. An automated approach will be needed to eliminate measurement subjectivity and handle the otherwise-overwhelming image quantities (up to billions of images) from near-future surveys. Our approach automatically describes spiral galaxy structure as a set of arcs, precisely describing spiral arm segment arrangement while retaining the flexibility needed to accommodate the observed wide variety of spiral galaxy structure. The largest existing quantitative measurements were manually-guided and encompassed fewer than 100 galaxies, while we have already applied our method to more than 29,000 galaxies. Our output matches previous information, both quantitatively over small existing samples, and qualitatively against human classifications from Galaxy Zoo.

We demonstrate that the system of Jovian planets (Sun+Jupiter+Saturn+Uranus+Neptune), integrated for 200 million years as an isolated 5-body system using many sets of initial conditions all within the uncertainty bounds of their currently known positions, can display both chaos and near-integrability. The conclusion is consistent across four different integrators, including several comparisons against integrations utilizing quadruple precision. We demonstrate that the Wisdom-Holman symplectic map using simple symplectic correctors as implemented in Mercury 6.2 (Chambers 1999) gives a reliable characterization of the existence of chaos for a particular initial condition only with timesteps less than about 10 days, corresponding to about 400 steps per orbit. We also integrate the canonical DE405 initial condition out to 5 Gy, and show that it has a Lyapunov Time of 200--400 My, opening the remote possibility of accurate prediction of the Jovian planetary positions for 5 Gy.

The existence of chaos in the system of Jovian planets has been in question for the past 15 years. Various investigators have found Lyapunov times ranging from about 5 millions years upwards to infinity, with no clear reason for the discrepancy. In this paper, we resolve the issue. The position of the outer planets is known to only a few parts in 10 million. We show that, within that observational uncertainty, there exist Lyapunov timescales in the full range listed above. Thus, the ``true'' Lyapunov timescale of the outer Solar System cannot be resolved using current observations.

A shadow of a numerical solution to a chaotic system is an_exact_ solution to the equations of motion that remains close to the numerical solution for a long time. In a collisionless n-body system, we know that particle motion is governed by the global potential rather than by inter-particle interactions. As a result, the trajectory of each individual particle in the system is independently shadowable. It is thus meaningful to measure the number of particles that have shadowable trajectories as a function of time. We find that the number of shadowable particles decays exponentially with time as exp(-mu t), and that for eps in [~0.2,1] (in units of the local mean inter-particle separation $\bar n$), there is an explicit relationship between the decay constant mu, the timestep h of the leapfrog integrator, the softening eps, and the number of particles N in the simulation. Thus, given N and eps, it is possible to pre-compute the timestep h necessary to acheive a desired fraction of shadowable particles after a given length of simulation time. We demonstrate that a large fraction of particles remain shadowable over ~100 crossing times even if particles travel up to about 1/3 of the softening length per timestep. However, a sharp decrease in the number of shadowable particles occurs if the timestep increases to allow particles to travel further than 1/3 the softening length in one timestep, or if the softening is decreased below ~0.2$\bar n$.

A {\it shadow} is an exact solution to a chaotic system of equations that remains close to a numerically computed solution for a long time, ending in a {\it glitch}. We study the distribution of shadow durations at low dimension and how shadow durations scale as dimension increases up to 300 in a slightly simplified gravitational n-body system. We find that ``softened'' systems are shadowable for many tens of crossing times even for large n, while in an ``unsoftened'' system each particle encounters glitches independently as a Poisson process, giving shadow durations that scale as 1/n.

Simple ``solar systems'' are generated with planetary orbital radii r distributed uniformly random in log(r) between 0.2 and 50 AU. A conservative stability criterion is imposed by requiring that adjacent planets are separated by a minimum distance of k Hill radii, for values of k ranging from 1 to 8. Least-squares fits of these systems to generalized Bode laws are performed, and compared to the fit of our own Solar System. We find that this stability criterion, and other ``radius-exclusion'' laws, generally produce approximately geometrically spaced planets that fit a Titius-Bode law about as well as our own Solar System. We then allow the random systems the same exceptions that have historically been applied to our own Solar System. Namely, one gap may be inserted, similar to the gap between Mars and Jupiter, and up to 3 planets may be ``ignored'', similar to how some forms of Bode's law ignore Mercury, Neptune, and Pluto. With these particular exceptions, we find that our Solar System fits significantly better than the random ones. However, we believe that this choice of exceptions, designed specifically to give our own Solar System a better fit, gives it an unfair advantage that would be lost if other exception rules were used. We conclude that the significance of Bode's law is simply that stable planetary systems tend to be regularly spaced; this conclusion could be strengthened by the use of more stringent methods of rejecting unstable solar systems, such as long-term orbit integrations.