Comparison of the physical acoustic channel response of a line array of thin rectangular bars to an equivalent model of thin vibrating rectangular pistons.

The resolution of an array is determined by the number and spatial distance of apertures (channels) within the array and the geometry of each aperture. The accurate design of acoustic sensing arrays relies on an a prioiri estimate of the expected far field radiation pattern of reciprocally behaved elements chosen for each aperture which is difficult to calculate under damped and loaded conditions. The estimated response of one channel of a vertical line array, when modeled as a series of rectangular vibrating pistons on a rigid baffle, is compared to the measured response of one channel of a line array comprised of a series of thin rectangular bars under load and operating off resonance. Although simple modeling can predict the 3dB main lobe width of the channel with some accuracy, loading and damping effects will alter the individual element response and hence the sensitivity of the array and side lobe magnitudes when off axis steering. This is important to note when estimating array gain and noise contributions from sidelobes under steered conditions

The (In)Stability of Planetary Systems

We present results of numerical simulations which examine the dynamical stability of known planetary systems, a star with two or more planets. First we vary the initial conditions of each system based on observational data. We then determine regions of phase space which produce stable planetary configurations. For each system we perform 1000 ~1 million year integrations. We examine upsilon And, HD83443, GJ876, HD82943, 47UMa, HD168443, and the solar system (SS). We find that the resonant systems, 2 planets in a first order mean motion resonance, (HD82943 and GJ876) have very narrow zones of stability. The interacting systems, not in first order resonance, but able to perturb each other (upsilon And, 47UMa, and SS) have broad regions of stability. The separated systems, 2 planets beyond 10:1 resonance, (we only examine HD83443 and HD168443) are fully stable. Furthermore we find that the best fits to the interacting and resonant systems place them very close to unstable regions. The boundary in phase space between stability and instability depends strongly on the eccentricities, and (if applicable) the proximity of the system to perfect resonance. In addition to million year integrations, we also examined stability on ~100 million year timescales. For each system we ran ~10 long term simulations, and find that the Keplerian fits to these systems all contain configurations which may be regular on this timescale.Comment: 37 pages, 49 figures, 13 tables, submitted to Ap

QYMSYM: A GPU-Accelerated Hybrid Symplectic Integrator That Permits Close Encounters

We describe a parallel hybrid symplectic integrator for planetary system integration that runs on a graphics processing unit (GPU). The integrator identifies close approaches between particles and switches from symplectic to Hermite algorithms for particles that require higher resolution integrations. The integrator is approximately as accurate as other hybrid symplectic integrators but is GPU accelerated.Comment: 17 pages, 2 figure

A Lagrangian Integrator for Planetary Accretion and Dynamics (LIPAD)

We presented the first particle based, Lagrangian code that can follow the collisional/accretional/dynamical evolution of a large number of km-sized planetesimals through the entire growth process to become planets. We refer to it as the 'Lagrangian Integrator for Planetary Accretion and Dynamics' or LIPAD. LIPAD is built on top of SyMBA, which is a symplectic $N$-body integrator. In order to handle the very large number of planetesimals required by planet formation simulations, we introduce the concept of a `tracer' particle. Each tracer is intended to represent a large number of disk particles on roughly the same orbit and size as one another, and is characterized by three numbers: the physical radius, the bulk density, and the total mass of the disk particles represented by the tracer. We developed statistical algorithms that follow the dynamical and collisional evolution of the tracers due to the presence of one another. The tracers mainly dynamically interact with the larger objects (`planetary embryos') in the normal N-body way. LIPAD's greatest strength is that it can accurately model the wholesale redistribution of planetesimals due to gravitational interaction with the embryos, which has recently been shown to significantly affect the growth rate of planetary embryos . We verify the code via a comprehensive set of tests which compare our results with those of Eulerian and/or direct N-body codes.Comment: Accepted to the Astronomical Journal. See http://www.boulder.swri.edu/~hal/LIPAD.html for more detail including animation

Drawing Planar Graphs with a Prescribed Inner Face

Given a plane graph $G$ (i.e., a planar graph with a fixed planar embedding) and a simple cycle $C$ in $G$ whose vertices are mapped to a convex polygon, we consider the question whether this drawing can be extended to a planar straight-line drawing of $G$. We characterize when this is possible in terms of simple necessary conditions, which we prove to be sufficient. This also leads to a linear-time testing algorithm. If a drawing extension exists, it can be computed in the same running time

A Symplectic Integrator for Hill's Equations

Hill's equations are an approximation that is useful in a number of areas of astrophysics including planetary rings and planetesimal disks. We derive a symplectic method for integrating Hill's equations based on a generalized leapfrog. This method is implemented in the parallel N-body code, PKDGRAV and tested on some simple orbits. The method demonstrates a lack of secular changes in orbital elements, making it a very useful technique for integrating Hill's equations over many dynamical times. Furthermore, the method allows for efficient collision searching using linear extrapolation of particle positions.Comment: 15 pages, 2 figures; minor revisions; accepted for publication in the Astronomical Journa

On the formation of hot Neptunes and super-Earths

The discovery of short-period Neptune-mass objects, now including the remarkable system HD69830 (Lovis et al. 2006) with three Neptune analogues, raises difficult questions about current formation models which may require a global treatment of the protoplanetary disc. Several formation scenarios have been proposed, where most combine the canonical oligarchic picture of core accretion with type I migration (e.g. Terquem & Papaloizou 2007) and planetary atmosphere physics (e.g. Alibert et al. 2006). To date, published studies have considered only a small number of progenitors at late times. This leaves unaddressed important questions about the global viability of the models. We seek to determine whether the most natural model -- namely, taking the canonical oligarchic picture of core accretion and introducing type I migration -- can succeed in forming objects of 10 Earth masses and more in the innermost parts of the disc. This problem is investigated using both traditional semianalytic methods for modelling oligarchic growth as well as a new parallel multi-zone N-body code designed specifically for treating planetary formation problems with large dynamic range (McNeil & Nelson 2009). We find that it is extremely difficult for oligarchic tidal migration models to reproduce the observed distribution. Even under many variations of the typical parameters, we form no objects of mass greater than 8 Earth masses. By comparison, it is relatively straightforward to form icy super-Earths. We conclude that either the initial conditions of the protoplanetary discs in short-period Neptune systems were substantially different from the standard disc models we used, or there is important physics yet to be understood.Comment: 19 pages, 18 figures. Final version accepted to MNRAS 30 September 200

From planetesimals to terrestrial planets: N-body simulations including the effects of nebular gas and giant planets

We present results from a suite of N-body simulations that follow the accretion history of the terrestrial planets using a new parallel treecode that we have developed. We initially place 2000 equal size planetesimals between 0.5--4.0 AU and the collisional growth is followed until the completion of planetary accretion (> 100 Myr). All the important effect of gas in laminar disks are taken into account: the aerodynamic gas drag, the disk-planet interaction including Type I migration, and the global disk potential which causes inward migration of secular resonances as the gas dissipates. We vary the initial total mass and spatial distribution of the planetesimals, the time scale of dissipation of nebular gas, and orbits of Jupiter and Saturn. We end up with one to five planets in the terrestrial region. In order to maintain sufficient mass in this region in the presence of Type I migration, the time scale of gas dissipation needs to be 1-2 Myr. The final configurations and collisional histories strongly depend on the orbital eccentricity of Jupiter. If today's eccentricity of Jupiter is used, then most of bodies in the asteroidal region are swept up within the terrestrial region owing to the inward migration of the secular resonance, and giant impacts between protoplanets occur most commonly around 10 Myr. If the orbital eccentricity of Jupiter is close to zero, as suggested in the Nice model, the effect of the secular resonance is negligible and a large amount of mass stays for a long period of time in the asteroidal region. With a circular orbit for Jupiter, giant impacts usually occur around 100 Myr, consistent with the accretion time scale indicated from isotope records. However, we inevitably have an Earth size planet at around 2 AU in this case. It is very difficult to obtain spatially concentrated terrestrial planets together with very late giant impacts.Comment: 51 pages, 19 figures, 2 tables, published in Icaru