We study the late-time evolution of the central regions of two Milky Way-like simulations of galaxies formed in a cosmological context, one hosting a fast bar and the other a slow one. We find that bar length, R_b, measurements fluctuate on a dynamical timescale by up to 100%, depending on the spiral structure strength and measurement threshold. The bar amplitude oscillates by about 15%, correlating with R_b. The Tremaine-Weinberg-method estimates of the bars' instantaneous pattern speeds show variations around the mean of up to ~20%, typically anti-correlating with the bar length and strength. Through power spectrum analyses, we establish that these bar pulsations, with a period in the range ~60-200 Myr, result from its interaction with multiple spiral modes, which are coupled with the bar. Because of the presence of odd spiral modes, the two bar halves typically do not connect at exactly the same time to a spiral arm, and their individual lengths can be significantly offset. We estimated that in about 50% of bar measurements in Milky Way-mass external galaxies, the bar lengths of SBab type galaxies are overestimated by ~15% and those of SBbc types by ~55%. Consequently, bars longer than their corotation radius reported in the literature, dubbed "ultra-fast bars", may simply correspond to the largest biases. Given that the Scutum-Centaurus arm is likely connected to the near half of the Milky Way bar, recent direct measurements may be overestimating its length by 1-1.5 kpc, while its present pattern speed may be 5-10 km/s/kpc smaller than its time-averaged value.

Spin off events and impacts can eject boulders from an asteroid surface and rubble pile asteroids can accumulate from debris following a collision between large asteroids. These processes produce a population of gravitational bound objects in orbit that can impact an asteroid surface at low velocity and with a distribution of impact angles. We present laboratory experiments of low velocity spherical projectiles into a fine granular medium, sand. We delineate velocity and impact angles giving ricochets, those giving projectiles that roll-out from the impact crater and those that stop within their impact crater. With high speed camera images and fluorescent markers on the projectiles we track spin and projectile trajectories during impact. We find that the projectile only reaches a rolling without slipping condition well after the marble has reached peak penetration depth. The required friction coefficient during the penetration phase of impact is 4-5 times lower than that of the sand suggesting that the sand is fluidized near the projectile surface during penetration. We find that the critical grazing impact critical angle dividing ricochets from roll-outs, increases with increasing impact velocity. The critical angles for ricochet and for roll-out as a function of velocity can be matched by an empirical model during the rebound phase that balances a lift force against gravity. We estimate constraints on projectile radius, velocity and impact angle that would allow projectiles on asteroids to ricochet or roll away from impact, finally coming to rest distant from their initial impact sites.

With backwards orbit integration we estimate birth locations of young stellar associations and moving groups identified in the solar neighborhood that are younger than 70 Myr. The birth locations of most of these stellar associations are at smaller galactocentric radius than the Sun, implying that their stars moved radially outwards after birth. Exceptions to this rule include are the Argus and Octans associations which formed outside the Sun's Galactocentric radius. Variations in birth heights of the stellar associations suggest that they were born in a corrugated disk of molecular clouds, similar to that inferred from the current filamentary molecular cloud distribution and dust extinction maps. Multiple spiral arm features with different but near corotation pattern speeds and at different heights could account for the stellar association birth sites. We find that the young stellar associations are located in between peaks in the UV stellar velocity distribution for stars in the solar neighborhood. This would be expected if they were born in a spiral arm which perturbs stellar orbits that cross it. In contrast, stellar associations seem to be located near peaks in the vertical phase space distribution, suggesting that the gas in which stellar associations are born moves vertically together with the low velocity dispersion disk stars.

Mass-spring model simulations are used to investigate past spin states of a viscoelastic Phobos and Deimos. From an initially tidally locked state, we find crossing of a spin-orbit resonance with Mars or a mean motion resonance with each other does not excite tumbling in Phobos or Deimos. However, once tumbling our simulations show that these moons can remain so for an extended period and during this time their orbital eccentricity can be substantially reduced. We attribute the tendency for simulations of an initially tumbling viscoelastic body to drop into spin-synchronous state at very low eccentricity to the insensitivity of the tumbling chaotic zone volume to eccentricity. After a tumbling body enters the spin synchronous resonance, it can exhibit long lived non-principal axis rotation and this too can prolong the period of time with enhanced tidally generated energy dissipation. The low orbital eccentricities of Phobos and Deimos could in part be due to spin excitation by nearly catastrophic impacts rather than tidal evolution following orbital resonance excitation.

We present ALMA CO(1-0) observations toward the dust lane of the nearest elliptical and radio galaxy, NGC 5128 (Centaurus A), with high angular resolution ($\sim$ 1 arcsec, or 18 pc), including information from large to small spatial scales and total flux. We find a total molecular gas mass of 1.6$\times$10$^9$ $M_\odot$ and we reveal the presence of filamentary components more extended than previously seen, up to a radius of 4 kpc. We find that the global star formation rate is $\sim$1 \Msol yr$^{-1}$, which yields a star formation efficiency (SFE) of 0.6 Gyr$^{-1}$ (depletion time $\tau =$1.5 Gyr), similar to those in disk galaxies. We show the most detailed view to date (40\,pc resolution) of the relation between molecular gas and star formation within the stellar component of an elliptical galaxy, from several kpc scale to the circumnuclear region close to the powerful radio jet. Although on average the SFEs are similar to those of spiral galaxies, the circumnuclear disk (CND) presents SFEs of 0.3 Gyr$^{-1}$, lower by a factor of 4 than the outer disk. The low SFE in the CND is in contrast to the high SFEs found in the literature for the circumnuclear regions of some nearby disk galaxies with nuclear activity, probably as a result of larger shear motions and longer AGN feedback. The higher SFEs in the outer disk suggests that only central molecular gas or filaments with sufficient density and strong shear motions will remain in $\sim$1 Gyr, which will later result in the compact molecular distributions and low SFEs usually seen in other giant ellipticals with cold gas.

We consider how an impact generated seismic pulse affects the surface of an asteroid distant from the impact site. With laboratory experiments on dry polydisperse gravel mixtures, we track the trajectories of particles ejected from the surface by a single strong upward propagating pressure pulse. High speed video images show that ejecta trajectories are independent of particle size, and collisions primarily take place upon landing. When they land particles are ballistically sorted, as proposed by Shinbrot et al. (2017), leaving larger particles on the surface and smaller particles more widely dispersed. A single strong pulse can leave previously buried boulders stranded on the surface. Boulder stranding due to an impact excited seismic pulse is an additional mechanism that could leave large boulders present on the surface of rubble asteroids such as 162173 Ryugu, 101955 Bennu and 25143 Itokawa.

We show that Smoothed Particle Hydrodynamics (SPH) simulations of dwarf galaxies interacting with a Milky Way-like disk produce moving groups in the simulated stellar disk. We compare the observed disk moving groups in the Solar neighborhood with velocity perturbations produced in SPH simulations with and without dwarf galaxy interactions. Our simulations include three different cases: one that includes three satellites -- with a Sagittarius-like dwarf galaxy interaction, as well as the Large and Small Magellanic Clouds; another that orbits in the plane of the Milky Way, and a null case, i.e., one with no dwarf galaxy interaction. The velocity sub-structures are identified using an algorithm to identify moving groups based on velocity over-densities in the plane of the disk. The properties of the identified moving groups change as the interacting simulations evolve. There are clear differences in the number and distribution of moving groups formed in the interacting simulations relative to the isolated simulation that we have analyzed. Our analysis suggests that some of the moving groups in the Milky Way may have formed due to dynamical interactions with perturbing dwarf satellites.

Using viscoelastic mass spring model simulations to track heat distribution inside a tidally perturbed body, we measure the asymmetry of heating in the crust of a spin synchronous Moon in eccentric orbit about the Earth. With the Moon within a few Earth radii of the Earth, we find that tidal heating per unit area in a lunar crustal shell is asymmetric due to the octupole order moment in the Earth's tidal field and is 10 to 20\% higher on its near side than on its far side. Tidal heating reduces the crustal basal heat flux and the rate of magma ocean crystallization. Assuming that the local crustal growth rate depends on the local basal heat flux and the distribution of tidal heating in latitude and longitude, a heat conductivity model illustrates that a moderately asymmetric and growing lunar crust could maintain its asymmetry. Near/far side asymmetric tidal heating could serve as a trigger for uneven lunar crustal growth.

Using a damped mass-spring model, we simulate wobble of spinning homogeneous viscoelastic ellipsoids undergoing non-principal axis rotation. Energy damping rates are measured for oblate and prolate bodies with different spin rates, spin states, viscoelastic relaxation timescales, axis ratios, and strengths. Analytical models using a quality factor by Breiter et al. (2012) and for the Maxwell rheology by Frouard & Efroimsky (2018) match our numerical measurements of the energy dissipation rate after we modify their predictions for the numerically simulated Kelvin-Voigt rheology. Simulations of nearly spherical but wobbling bodies with hard and soft cores show that the energy dissipation rate is more sensitive to the material properties in the core than near the surface. The sensitivity to viscoelastic model implies that inferred statistics of tumbling lifetimes in asteroids might be interpreted in terms of differences in their material properties.

We use the second data releases of the ESA Gaia astrometric survey and the high-resolution GALAH spectroscopic survey to analyse the structure of our Galaxy's disc components. With GALAH, we separate the alpha-rich and alpha-poor discs (with respect to Fe), which are superposed in both position and velocity space, and examine their distributions in action space. We study the distribution of stars in the zV_z phase plane, for both V_phi and V_R, and recover the remarkable "phase spiral" discovered by Gaia. We identify the anticipated quadrupole signature in zV_z of a tilted velocity ellipsoid for stars above and below the Galactic plane. By connecting our work with earlier studies, we show that the phase spiral is likely to extend well beyond the narrow solar neighbourhood cylinder in which it was found. The phase spiral is a signature of corrugated waves that propagate through the disc, and the associated non-equilibrium phase mixing. The radially asymmetric distribution of stars involved in the phase spiral reveals that the corrugation, which is mostly confined to the alpha-poor disc, grows in z-amplitude with increasing radius. We present new simulations of tidal disturbance of the Galactic disc by the Sagittarius (Sgr) dwarf. The effect on the zV_z phase plane lasts >2 Gyr but a subsequent disc crossing wipes out the coherent structure. We find that the phase spiral was excited < 0.5 Gyr ago by an object like Sgr with total mass 3 x 10^10 Msun (stripped down from 5 x 10^10 Msun when it first entered the halo) passing through the plane.

We consider an impact on an asteroid that is energetic enough to cause resurfacing by seismic reverberation and just below the catastrophic disruption threshold, assuming that seismic waves are not rapidly attenuated. In asteroids with diameter less than 1 km we identify a regime where rare energetic impactors can excite seismic waves with frequencies near those of the asteroid's slowest normal modes. In this regime, the distribution of seismic reverberation is not evenly distributed across the body surface. With mass-spring model elastic simulations, we model impact excitation of seismic waves with a force pulse exerted on the surface and using three different asteroid shape models. The simulations exhibit antipodal focusing and normal mode excitation. If the impulse excited vibrational energy is long lasting, vibrations are highest at impact point, its antipode and at high surface elevations such as an equatorial ridge. A near equatorial impact launches a seismic impulse on a non-spherical body that can be focused on two additional points on an the equatorial ridge. We explore simple flow models for the morphology of vibration induced surface slumping. We find that the initial seismic pulse is unlikely to cause large shape changes. Long lasting seismic reverberation on Bennu caused by a near equatorial impact could have raised the height of its equatorial ridge by a few meters and raised two peaks on it, one near impact site and the other near its antipode.

The solar neighborhood contains disc stars that have recently crossed spiral arms in the Galaxy. We propose that boundaries in local velocity distributions separate stars that have recently crossed or been perturbed by a particular arm from those that haven't. Ridges in the stellar velocity distributions constructed from the second Gaia data release trace orbits that could have touched nearby spiral arms at apocentre or pericentre. The multiple ridges and arcs seen in local velocity distributions are consistent with the presence of multiple spiral features and different pattern speeds and imply that the outer Galaxy is flocculent rather than grand design.

Before the publication of the Gaia DR2 we confirmed with RAVE and TGAS an observation recently made with the GALAH survey by Quillen ey al. concerning the Coma Berenices moving group in the Solar neighbourhood, namely that it is only present at negative Galactic latitudes. This allowed us to show that it is coherent in vertical velocity, providing a first evidence for incomplete vertical phase-mixing. We estimated for the first time from dynamical arguments that the moving group must have formed at most ~ 1.5 Gyr ago, and related this to a pericentric passage of the Sagittarius dwarf satellite galaxy. The present note is a rewritten version of the original arXiv post on this result now also including a confirmation of our finding with Gaia DR2.

Using GALAH survey data of nearby stars, we look at how structure in the planar (u,v) velocity distribution depends on metallicity and on viewing direction within the Galaxy. In nearby stars, with distance d < 1 kpc, the Hercules stream is most strongly seen in higher metallicity stars [Fe/H] > 0.2. The Hercules stream peak v value depends on viewed galactic longitude, which we interpret as due to the gap between the stellar stream and more circular orbits being associated with a specific angular momentum value of about 1640 km/s kpc. The association of the gap with a particular angular momentum value supports a bar resonant model for the Hercules stream. Moving groups previously identified in Hipparcos observations are easiest to see in stars nearer than 250 pc, and their visibility and peak velocities in the velocity distributions depends on both viewing direction (galactic longitude and hemisphere) and metallicity. We infer that there is fine structure in local velocity distributions that varies over distances of a few hundred pc in the Galaxy.

Using tracer particles embedded in self-gravitating shearing sheet N-body simulations, we investigate the distance in guiding centre radius that stars or star clusters can migrate in a few orbital periods. The standard deviations of guiding centre distributions and maximum migration distances depend on the Toomre or critical wavelength and the contrast in mass surface density caused by spiral structure. Comparison between our simulations and estimated guiding radii for a few young super-solar metallicity open clusters, including NGC 6583, suggests that the contrast in mass surface density in the solar neighbourhood has standard deviation (in the surface density distribution) divided by mean of about 1/4 and larger than measured using COBE data by Drimmel and Spergel. Our estimate is consistent with a standard deviation of $\sim$0.07 dex in the metallicities measured from high-quality spectroscopic data for 38 young open clusters (<1 Gyr) with mean galactocentric radius 7-9 kpc.

A Hamiltonian model is constructed for the spin axis of a planet perturbed by a nearby planet with both planets in orbit about star. We expand the planet-planet gravitational potential perturbation to first order in orbital inclinations and eccentricities, finding terms describing spin resonances involving the spin precession rate and the two planetary mean motions. Convergent planetary migration allows the spinning planet to be captured into spin resonance. With initial obliquity near zero, the spin resonance can lift the planet's obliquity to near 90 or 180 degrees depending upon whether the spin resonance is first or zero-th order in inclination. Past capture of Uranus into such a spin resonance could give an alternative non-collisional scenario accounting for Uranus's high obliquity. However we find that the time spent in spin resonance must be so long that this scenario cannot be responsible for Uranus's high obliquity. Our model can be used to study spin resonance in satellite systems. Our Hamiltonian model explains how Styx and Nix can be tilted to high obliquity via outward migration of Charon, a phenomenon previously seen in numerical simulations.

We present the distribution and kinematics of the molecular gas in the circumnuclear disk (CND, 400 pc x 200 pc) of Centaurus A with resolutions of ~5 pc (0.3 arcsec) and shed light onto the mechanism feeding the Active Galactic Nucleus (AGN) using CO(3-2), HCO+(4-3), HCN(4-3), and CO(6-5) observations obtained with ALMA. Multiple filaments or streamers of tens to a hundred parsec scale exist within the CND, which form a ring-like structure with an unprojected diameter of 9 x 6 arcsec (162pc x 108pc) and a position angle PA = 155deg. Inside the nuclear ring, there are two leading and straight filamentary structures with lengths of about 30-60pc at PA = 120deg on opposite sides of the AGN, with a rotational symmetry of 180deg and steeper position-velocity diagrams, which are interpreted as nuclear shocks due to non-circular motions. Along the filaments, and unlike other nearby AGNs, several dense molecular clumps present low HCN/HCO+(4-3) ratios (~0.5). The filaments abruptly end in the probed transitions at r = 20pc from the AGN, but previous near-IR H2 (J=1-0) S(1) maps show that they continue in an even ~1000 K), winding up in the form of nuclear spirals, and forming an inner ring structure with another set of symmetric filaments along the N-S direction and within r = 10pc. The molecular gas is governed primarily by non-circular motions, being the successive shock fronts at different scales where loss of angular momentum occurs, a mechanism which may feed efficiently powerful radio galaxies down to parsec scales.

New Horizons mission observations show that the small satellites Styx, Nix, Kerberos and Hydra, of the Pluto-Charon system, have not tidally spun-down to near synchronous spin states and have high obliquities with respect to their orbit about the Pluto-Charon binary (Weaver et al. 2016). We use a damped mass-spring model within an N-body simulation to study spin and obliquity evolution for single spinning non-round bodies in circumbinary orbit. Simulations with tidal dissipation alone do not show strong obliquity variations from tidally induced spin-orbit resonance crossing and this we attribute to the high satellite spin rates and low orbital eccentricities. However, a tidally evolving Styx exhibits intermittent obliquity variations and episodes of tumbling. During a previous epoch where Charon migrated away from Pluto, the minor satellites could have been trapped in orbital mean motion inclination resonances. An outward migrating Charon induces large variations in Nix and Styx's obliquities. The cause is a commensurability between the mean motion resonance frequency and the spin precession rate of the spinning body. As the minor satellites are near mean motion resonances, this mechanism could have lifted the obliquities of all four minor satellites. The high obliquities need not be primordial if the minor satellites were at one time captured into mean motion resonances.

Until 2013, only the giant planets were known to host ring systems. In June 2013, a stellar occulation revealed the presence of narrow and dense rings around Chariklo, a small Centaur object that orbits between Saturn and Uranus. Meanwhile, the Cassini spacecraft revealed evidence for the possible past presence of rings around the Saturnian satellites Rhea and Iapetus. Mars and Pluto are expected to have tenuous dusty rings, though they have so far evaded detection. More remotely, transit events observed around a star in 2007 may have revealed for the first time exoplanetary rings around a giant planet orbiting that star. So, evidence is building to show that rings are more common features in the universe than previously thought. Several interesting issues arising from the discovery (or suspicion) of new ring systems are described in this chapter.

We revisit the nature of large dips in flux from extinction by dusty circumstellar material that is observed by Kepler for many young stars in the Upper Sco and $\rho$ Oph star formation regions. These young, low-mass "dipper" stars are known to have low accretion rates and primarily host moderately evolved dusty circumstellar disks. Young low mass stars often exhibit rotating star spots that cause quasi-periodic photometric variations. We found no evidence for periods associated with the dips that are different from the star spot rotation period in spectrograms constructed from the light curves. The material causing the dips in most of these light curves must be approximately corotating with the star. We find that disk temperatures computed at the disk corotation radius are cool enough that dust should not sublime. Crude estimates for stellar magnetic field strengths and accretion rates are consistent with magnetospheric truncation near the corotation radius. Magnetospheric truncation models can explain why the dips are associated with material near corotation and how dusty material is lifted out of the midplane to obscure the star which would account for the large fraction of young low mass stars that are dippers. We propose that variations in disk orientation angle, stellar magnetic field dipole tilt axis, and disk accretion rate are underlying parameters accounting for differences in the dipper light curves.

We use numerical simulations to measure the sensitivity of the tidal spin down rate of a homogeneous triaxial ellipsoid to its axis ratios by comparing the drift rate in orbital semi-major axis to that of a spherical body with the same mass, volume and simulated rheology. We use a mass-spring model approximating a viscoelastic body spinning around its shortest body axis, with spin aligned with orbital spin axis, and in circular orbit about a point mass. The torque or drift rate can be estimated from that predicted for a sphere with equivalent volume if multiplied by $0.5 (1 + b^4/a^4)(b/a)^{-4/3} (c/a)^{-\alpha_c}$ where $b/a$ and $c/a$ are the body axis ratios and index $\alpha_c \approx 1.05$ is consistent with the random lattice mass spring model simulations but $\alpha_c = 4/3$ suggested by scaling estimates. A homogeneous body with axis ratios 0.5 and and 0.8, like Haumea, has orbital semi-major axis drift rate about twice as fast as a spherical body with the same mass, volume and material properties. A simulation approximating a mostly rocky body but with 20\% of its mass as ice concentrated at its ends has a drift rate 10 times faster than the equivalent homogeneous rocky sphere. However, this increase in drift rate is not enough to allow Haumea's satellite, Hi'iaka, to have tidally drifted away from Haumea to its current orbital semi-major axis.

Supermassive black holes in galaxy centres can grow by the accretion of gas, liberating energy that might regulate star formation on galaxy-wide scales. The nature of the gaseous fuel reservoirs that power black hole growth is nevertheless largely unconstrained by observations, and is instead routinely simplified as a smooth, spherical inflow of very hot gas. Recent theory and simulations instead predict that accretion can be dominated by a stochastic, clumpy distribution of very cold molecular clouds - a departure from the "hot mode" accretion model - although unambiguous observational support for this prediction remains elusive. Here we report observations that reveal a cold, clumpy accretion flow towards a supermassive black hole fuel reservoir in the nucleus of the Abell 2597 Brightest Cluster Galaxy (BCG), a nearby (redshift z=0.0821) giant elliptical galaxy surrounded by a dense halo of hot plasma. Under the right conditions, thermal instabilities can precipitate from this hot gas, producing a rain of cold clouds that fall toward the galaxy's centre, sustaining star formation amid a kiloparsec-scale molecular nebula that inhabits its core. The observations show that these cold clouds also fuel black hole accretion, revealing "shadows" cast by the molecular clouds as they move inward at about 300 kilometres per second towards the active supermassive black hole in the galaxy centre, which serves as a bright backlight. Corroborating evidence from prior observations of warmer atomic gas at extremely high spatial resolution, along with simple arguments based on geometry and probability, indicate that these clouds are within the innermost hundred parsecs of the black hole, and falling closer towards it.

We present new rest frame UV and visible observations of 22 high-redshift (1 < z < 2.5) 3C radio galaxies and QSOs obtained with the Hubble Space Telescope's (HST) Wide Field Camera 3 (WFC3) instrument. Using a custom data reduction strategy in order to assure the removal of cosmic rays, persistence signal, and other data artifacts, we have produced high-quality science-ready images of the targets and their local environments. We observe targets with regions of UV emission suggestive of active star formation. In addition, several targets exhibit highly distorted host galaxy morphologies in the rest frame visible images. Photometric analyses reveals that brighter QSOs tend to be generally redder than their dimmer counterparts. Using emission line fluxes from the literature, we estimate that emission line contamination is relatively small in the rest frame UV images for the QSOs. Using archival VLA data, we have also created radio map overlays for each of our targets, allowing for analysis of the optical and radio axes alignment.

Signatures of "evaporative" winds from exo-planets on short (hot) orbits around their host star have been observed in a number of systems. In this paper we present global AMR simulations that track the launching of the winds, their expansion through the circumstellar environment, and their interaction with a stellar wind. We focus on purely hydrodynamic flows including the anisotropy of the wind launching and explore the orbital/fluid dynamics of the resulting flows in detail. In particular we find that a combination of the tidal and Coriolis forces strongly distorts the planetary "Parker" wind creating "up-orbit" and "down-orbit" streams. We characterize the flows in terms of their orbital elements which change depending on their launch position on the planet. We find that the anisotropy in the atmospheric temperature leads to significant backflow on to the planet. The planetary wind interacts strongly with the stellar wind creating instabilities that cause eventual deposition of planetary gas onto the star. We present synthetic observations of both transit and absorption line-structure for our simulations. For our initial conditions, we find that the orbiting wind material produces absorption signatures at significant distances from the planet and substantial orbit to orbit variability. Ly-{\alpha} absorption shows red and blueshifted features out to 70 km/s. Finally, using semi-analytic models we constrain the effect of radiation pressure, given the approximation of uniform stellar absorption.

Close tidal encounters among large planetesimals and moons should have been more common than grazing or normal impacts. Using a mass spring model within an N-body simulation, we simulate the deformation of the surface of an elastic spherical body caused by a close parabolic tidal encounter with a body that has similar mass as that of the primary body. Such an encounter can induce sufficient stress on the surface to cause brittle failure of an icy crust and simulated fractures can extend a large fraction of the radius of body. Strong tidal encounters may be responsible for the formation of long graben complexes and chasmata in ancient terrain of icy moons such as Dione, Tethys, Ariel and Charon.

We use a set of high-resolution N-body simulations of the Galactic disk to study its interactions with the population of satellites predicted cosmologically. One simulation illustrates that multiple passages of massive satellites with different velocities through the disk generate a wobble, having the appearance of rings in face-on projections of the stellar disk. They also produce flares in the disk outer parts and gradually heat the disk through bending waves. A different numerical experiment shows that an individual satellite as massive as the Sagittarius dwarf galaxy passing through the disk will drive coupled horizontal and vertical oscillations of stars in underdense regions, with small significant associated heating. This experiment shows that vertical excursions of stars in these low-density regions can exceed 1 kpc in the Solar neighborhood, resembling the coherent vertical oscillations recently detected locally. They can also induce non-zero vertical streaming motions as large as 10-20 km s$^{-1}$, consistent with recent observations in the Galactic disk. This phenomenon appears as a local ring, with modest associated disk heating.

We investigate the plausibility of a cometary source of the unusual transits observed in the KIC 8462852 light curve. A single comet of similar size to those in our solar system produces a transit depth of the order of $10^{-3}$ lasting less than a day which is much smaller and shorter than the largest dip observed ($\sim20\%$ for $\sim3$ days), but a large, closely traveling cluster of comets can fit the observed depths and durations. We find that a series of large comet swarms, with all but one on the same orbit, provides a good fit for the KIC 8462852 data during Quarters 16 and 17, but does not explain the large dip observed during Quarter 8. However, the transit dips only loosely constrain the orbits and can be fit by swarms with periastrons differing by a factor of 10. To reach a transit depth of $\sim0.2$, the comets need to be in a close group of $\sim30$, if they are $\sim100$ km in radius or in a group of $\sim300$ if they are $\sim10$ km in radius. The total number of comets required to fit all of the dips is $\sim70$ $\sim$100 km or $\sim700$ $\sim10$ km comets. A single comet family from a tidally disrupted Ceres-sized progenitor or the start of a Late Heavy Bombardment period explains the last $\sim60$ days of the unusual KIC 8462852 light curve.

We use a damped mass-spring model within an N-body code to simulate the tidal evolution of the spin and orbit of a self-gravitating viscoelastic spherical body moving around a point-mass perturber. The damped mass-spring model represents a Kelvin-Voigt viscoelastic solid. We measure the tidal quality function (the dynamical Love number $\,k_2\,$ divided by the tidal quality factor $\,Q\,$) from the numerically computed tidal drift of the semimajor axis of the binary. The shape of $\,k_2/Q\,$, as a function of the principal tidal frequency, reproduces the kink shape predicted by Efroimsky (2012a; CeMDA 112$\,:\,$283) for the tidal response of near-spherical homogeneous viscoelastic rotators. We demonstrate that we can directly simulate the tidal evolution of spinning viscoelastic objects. In future, the mass-spring N-body model can be generalised to inhomogeneous and/or non-spherical bodies.

A dynamical characterization of dark matter dominated dwarf galaxies from their observed effects on galactic disks (i.e. Galactoseismology) has remained an elusive goal. Here, we present preliminary results from spectroscopic observations of three clustered Cepheid candidates identified from $K$-band light curves towards Norma. The average heliocentric radial velocity of these stars is $\sim$ 156 km/s, which is large and distinct from that of the Galaxy's stellar disk. These objects at $l \sim 333 ^\circ$ and $b \sim -1 ^\circ$ are therefore halo stars; using the $3.6~\micron$ period-luminosity relation of Type I Cepheids, they are at $\sim$ 73 kpc. Our ongoing $I$-band photometry indicates variability on the same time scale as the period determined from the $K_{s}$-band light curve. Distances determined from the $K$-band period-luminosity relation and the 3.6 $\micron$ period-luminosity relation are comparable. The observed radial velocity of these stars agrees roughly with predictions from dynamical models. If these stars are indeed members of the predicted dwarf galaxy that perturbed the outer HI disk of the Milky Way, this would mark the first application of Galactoseismology.

We report the identification of extended tidal debris potentially associated with the globular cluster NGC 3201, using the RAVE catalogue. We find the debris stars are located at a distance range of 1-7 kpc based on the forthcoming RAVE distance estimates. The derived space velocities and integrals of motion show interesting connections to NGC 3201, modulo uncertainties in the proper motions. Three stars, which are among the 4 most likely candidates for NGC 3201 tidal debris, are separated by 80 degrees on the sky yet are well matched by the 12 Gyr, [Fe/H] = -1.5 isochrone appropriate for the cluster. This is the first time tidal debris around this cluster has been reported over such a large spatial extent, with implications for the cluster$'$s origin and dynamical evolution.

We use test-particle integrations to show that epicyclic motions excited by a pericentre passage of a dwarf galaxy could account for bulk vertical velocity streaming motions recently observed in the Galactic stellar disc near the Sun. We use fixed potential test-particle integrations to isolate the role of phase wrapping of epicyclic perturbations from bending and breathing waves or modes, which require self-gravity to oscillate. Perturbations from a fairly massive Sagittarius dwarf galaxy, $M_d \sim 2.5 \times 10^{10} M_\odot$, are required to account for the size of the observed streaming motions from its orbital pericentre approximately a Gyr ago. A previous passage of the dwarf through the Galactic disc approximately 2.2 Gyr ago (with a then more massive dwarf galaxy) is less effective. If phase wrapping of epicyclic perturbations is responsible for stellar streaming motions in the Galactic disc, then there should be variations in velocity gradients on scales of a few kpc in the vicinity of the Sun.

We present initial results of a new campaign of simulations focusing on the interaction of planetary winds with stellar environments using Adaptive Mesh Refinement methods. We have confirmed the results of Stone & Proga 2009 that an azimuthal flow structure is created in the planetary wind due to day/night temperatures differences. We show that a backflow towards the planet will occur with a strength that depends on the escape parameter. When a stellar outflow is included, we see unstable bow waves forming through the outflow's interaction with the planetary wind.

The temperature of a circumbinary disc edge should undulate due to variations in illumination as a function of binary orbital phase. We explore circumbinary disc temperature variations as a source of broad-band infrared light curve variability. Approximating the wall of a circumbinary disc edge as a wide optically thick cylinder with surface temperature dependent on its illumination, we find that a binary comprised of 1 M$_\odot$ and 0.5 M$_\odot$ pre-main sequence stars in a $\sim$15.5 day period, would exhibit the largest amplitude variations of $\sim$9% at 3.77 and 4.68 {\mu}m as seen by a distant observer. The amplitude of variations and shape of the light curve is sensitive to the luminosity and mass ratios of the stars in the binary, the radius of the circumbinary disc clearing, the binary separation, and the orbital inclination. The light curve variations are smooth and very red with a non-sinusoidal shape for most of the parameter space explored. Possible morphologies include a single peak with a flat region, two peaks of different heights or a single dip.

We present a multiwavelength morphological analysis of star forming clouds and filaments in the central ($< 50$ kpc) regions of 16 low redshift ($z<0.3$) cool core brightest cluster galaxies (BCGs). New Hubble Space Telescope (HST) imaging of far ultraviolet continuum emission from young ($\sim 10$ Myr), massive ($> 5$ \Msol) stars reveals filamentary and clumpy morphologies, which we quantify by means of structural indices. The FUV data are compared with X-ray, Ly$\alpha$, narrowband H$\alpha$, broadband optical/IR, and radio maps, providing a high spatial resolution atlas of star formation locales relative to the ambient hot ($\sim10^{7-8}$ K) and warm ionised ($\sim 10^4$ K) gas phases, as well as the old stellar population and radio-bright AGN outflows. Nearly half of the sample possesses kpc-scale filaments that, in projection, extend toward and around radio lobes and/or X-ray cavities. These filaments may have been uplifted by the propagating jet or buoyant X-ray bubble, or may have formed {\it in situ} by cloud collapse at the interface of a radio lobe or rapid cooling in a cavity's compressed shell. The morphological diversity of nearly the entire FUV sample is reproduced by recent hydrodynamical simulations in which the AGN powers a self-regulating rain of thermally unstable star forming clouds that precipitate from the hot atmosphere. In this model, precipitation triggers where the cooling-to- freefall time ratio is $t_{\mathrm{cool}}/t_{\mathrm{ff}}\sim 10$. This condition is roughly met at the maxmial projected FUV radius for more than half of our sample, and clustering about this ratio is stronger for sources with higher star formation rates.

This is a brief rebuttal to arXiv:1502.03821, which claims to provide the first observational proof of dark matter interior to the solar circle. We point out that this result is not new, and can be traced back at least a quarter century.

We estimate the size and distribution of the parent populations for the 6 largest (at least 20 stars in the Solar neighborhood) chemical groups identified in the Chemical Tagging experiment by Mitschang et al.~2014. Stars in the abundance groups tend to lie near a boundary in angular momentum versus eccentricity space where the probability is highest for a star to be found in the Solar neighborhood and where orbits have apocenter approximately equal to the Sun's galactocentric radius. Assuming that the parent populations are uniformly distributed at all azimuthal angles in the Galaxy, we estimate that the parent populations of these abundance groups contain at least 200,000 members. The spread in angular momentum of the groups implies that the assumption of a uniform azimuthal distribution only fails for the two youngest groups and only for the highest angular momentum stars in them. The parent populations of three thin disk groups have narrow angular momentum distributions, but tails in the eccentricity and angular momentum distributions suggest that only a small fraction of stars have migrated and increased in eccentricity. In contrast, the parent populations of the thick disk groups exhibit both wide angular momentum and eccentricity distributions implying that both heating and radial migration has taken place.

Winds from short-period Earth and Neptune mass exoplanets, driven by high energy radiation from a young star, may evaporate a significant fraction of a planet's mass. If the momentum flux from the evaporative wind is not aligned with the planet/star axis, then it can exert a torque on the planet's orbit. Using steady-state one-dimensional evaporative wind models we estimate this torque using a lag angle that depends on the product of the speed of the planet's upper atmosphere and a flow timescale for the wind to reach its sonic radius. We also estimate the momentum flux from time-dependent one-dimensional hydrodynamical simulations. We find that only in a very narrow regime in planet radius, mass and stellar radiation flux is a wind capable of exerting a significant torque on the planet's orbit. Similar to the Yarkovsky effect, the wind causes the planet to drift outward if atmospheric circulation is prograde (super-rotating) and in the opposite direction if the circulation is retrograde. A close-in super Earth mass planet that loses a large fraction of its mass in a wind could drift a few percent of its semi-major axis. While this change is small, it places constraints on the evolution of resonant pairs such as Kepler 36 b and c.

Distant regions close to the plane of our Galaxy are largely unexplored by optical surveys as they are hidden by dust. We have used near-infrared data (that minimizes dust obscuration) from the ESO Public survey VISTA Variables of the Via Lactea (VVV) (Minniti et al. 2011; Saito et al. 2012; henceforth S12) to search for distant stars at low latitudes. We have discovered four Cepheid variables within an angular extent of one degree centered at Galactic longitude of $l = -27.4^\circ$ and Galactic latitude of $b = -1.08 ^\circ$. We use the tightly constrained period-luminosity relationship that these pulsating stars obey (Persson et al. 2004; Matsunaga et al. 2011) to derive distances. We infer an average distance to these Cepheid variables of 90 kpc. The Cepheid variables are highly clustered in angle (within one degree) and in distance (the standard deviation of the distances is 12 kpc). They are at an average distance of $\sim 2~\rm kpc$ from the plane and their maximum projected separation is $\sim 1~ \rm kpc$. These young ($\sim$ 100 Myr old), pulsating stars (Bono et al. 2005) are unexpected at such large distances from the Galactic disk, which terminates at $\sim$ 15 kpc (Minniti et al. 2011). The highly clustered nature in distance and angle of the Cepheid variables suggests that the stars may be associated with a dwarf galaxy, one that was earlier predicted by a dynamical analysis (Chakrabarti \& Blitz 2009).

We investigate the nature of the unusual eclipsing star OGLE LMC-ECL-11893 (OGLE J05172127-6900558) in the Large Magellanic Cloud recently reported by Dong et al. 2014. The eclipse period for this star is 468 days, and the eclipses exhibit a minimum of ~1.4 mag, preceded by a plateau of ~0.8 mag. Spectra and optical/IR photometry are consistent with the eclipsed star being a lightly reddened B9III star of inferred age ~150 Myr and mass of ~4 solar masses. The disk appears to have an outer radius of ~0.2 AU with predicted temperatures of ~1100-1400 K. We model the eclipses as being due to either a transiting geometrically thin dust disk or gaseous accretion disk around a secondary object; the debris disk produces a better fit. We speculate on the origin of such a dense circumstellar dust disk structure orbiting a relatively old low-mass companion, and on the similarities of this system to the previously discovered EE Cep.

Numerical integrations of the closely-packed inner Uranian satellite system show that variations in semi-major axes can take place simultaneously between three or four consecutive satellites. We find that the three-body Laplace angle values are distributed unevenly and have histograms showing structure, if the angle is associated with a resonant chain, with both pairs of bodies near first-order two-body resonances. Estimated three-body resonance libration frequencies can be only an order of magnitude lower than those of first-order resonances. Their strength arises from a small divisor from the distance to the first-order resonances and insensitivity to eccentricity, which make up for their dependence on moon mass. Three-body resonances associated with low-integer Laplace angles can also be comparatively strong due to the many multiples of the angle contributed from Fourier components of the interaction terms. We attribute small coupled variations in semi-major axis, seen throughout the simulation, to ubiquitous and weak three-body resonant couplings. We show that a system with two pairs of bodies in first-order mean-motion resonance can be transformed to resemble the well-studied periodically-forced pendulum with the frequency of a Laplace angle serving as a perturbation frequency. We identify trios of bodies and overlapping pairs of two-body resonances in each trio that have particularly short estimated Lyapunov timescales.

We search for systems hosting eclipsing discs using a complete sample of eclipsing binaries (EBs); those previously identified in the third phase of the Optical Gravitational Lensing Experiment (OGLE-III). Within a subsample of 2,823 high-cadence, high-photometric precision and large eclipsing depth detached EBs previously identified in the Large Magellanic Cloud (LMC), we find that the skewness and kurtosis of the light curves magnitude distribution within the primary eclipse can distinguish EBs hosting a disc from those without. Two systems with previously identified eclipsing discs (OGLE-LMC-ECL-11893 and OGLE-LMC-ECL-17782) are identified with near zero skewness ($|S|<0.5$) and positive kurtosis. No additional eclipsing disc systems were found in the OGLE-III LMC, Small Magellanic Cloud (SMC) or Galactic Disc (GD) EB light curves. We estimate that the fraction of detached near main-sequence LMC EBs (which have a primary with an $I$-band magnitude brighter than $\simeq 19\,\rm{mag}$) that host a disc is approximately 1/1000. As circumstellar disc lifetimes are short, we expected to primarily find eclipsing discs around young stars. In addition, as there is more room for a disc in a widely separated binary and because a disk close to a luminous star would be above the dust sublimation temperature, we expected to primarily find eclipsing discs in long period binaries. However, OGLE-LMC-ECL-17782 is a 13.3 day period B star system with a transient and hot ($\sim 6000\,\rm{K}$, $\sim 0.1\,\rm{AU}$ radius) disc and Scott et al. (in prep.) estimate an age of $150\,\rm{Myr}$ for OGLE-LMC-ECL-11893. Both discs are unexpected in the EB sample and impel explanation.

We infer the past orbit of the Sagittarius (Sgr) dwarf galaxy in the Milky Way halo by integrating backwards from its observed position and proper motions, including the effects of dynamical friction. Given measured proper motions, we show that there is a relation between the eccentricity ($e$) of Sgr's orbit and the mass of the Milky Way ($M_{T}$) in the limit of no dynamical friction. That relation can be fit by a power-law of the form: $e \approx 0.49 \left(M_{T}/10^{12} M_{\odot}\right)^{-0.88}$. At a fixed Milky Way mass, the dynamical friction term increases the mean eccentricity of the orbit and lowers the spread in eccentricities in proportion to the mass of the Sgr dwarf. We explore the implications of various observational constraints on Sgr's apocenter on the $e-M$ relation; Sgr masses outside the range $10^{9} M_{\odot} \la M_{\rm Sgr} \la 5 \times 10^{10} M_{\odot}$ are precluded, for Milky Way masses $\sim 1 - 2.5 \times 10^{12} M_{\odot}$. If Belokurov et al.'s (2014) observations represent the farthest point of Sgr's stream, then Milky Way masses in excess of $2 \times 10^{12} M_{\odot}$ are excluded for $M_{\rm Sgr} \la 10^{10} M_{\odot}$. Deeper observations of Sgr's tidal debris, from upcoming surveys such as GAIA, will allow better measurement of the Milky Way mass and of the Sgr dwarf.

We searched the light curves of over 40000 stars in the 2MASS calibration database, spanning approximately 4 years, for objects that have significant day long dimming events. We also searched the multi-color light curves for red dimming events that could be due to transient extinction. In the color independent sigma-limited search, we found 46 previously unknown eclipsing binaries, 6 periodic variable stars likely to be intrinsic pulsators and 21 young stellar objects in the $\rho$ Ophiucus star formation region previously studied by Parks et al. 2013. An additional 11 objects exhibited dimming events and most of these are unclassified. The search for red dimming events primarily reveals a population of low luminosity active galaxies that become bluer when they are brighter, and variable young stellar objects exhibiting high cross-correlation coefficients between color and brightness. The young stellar objects exhibit brightness and color variations in the direction of interstellar extinction whereas the active galaxies can have a bowed distribution in color and magnitude with reduced variation in color when the object is brightest. Among the objects that are usually quiescent (not strongly variable), we failed to find any dimming events deeper than 0.2 magnitude and lasting longer than a day. Two of the least embedded young stellar objects, however, dimmed by 0.2 mag for longer than a day without strong color variation. Having failed to find new exotic objects, we explored ways to eliminate commonly found objects so that a larger number of objects may be searched. We find that all but one eclipsing binary are excluded by requiring moderate color variation during a dimming event and most of the active galaxies and young stellar objects are excluded by placing a limit on the standard deviation of the magnitude distribution in the light curve.

The velocity dispersions of stars near the Sun are known to increase with stellar age, but age can be difficult to determine so a proxy like the abundance of alpha elements (e.g., Mg) with respect to iron, [alpha/Fe], is used. Here we report an unexpected behavior found in the velocity dispersion of a sample of giant stars from the RAdial Velocity Experiment (RAVE) survey with high quality chemical and kinematical information, in that it decreases strongly for stars with [Mg/Fe] > 0.4 dex (i.e., those that formed in the first Gyr of the Galaxy's life). These findings can be explained by perturbations from massive mergers in the early Universe, which have affected more strongly the outer parts of the disc, and the subsequent radial migration of stars with cooler kinematics from the inner disc. Similar reversed trends in velocity dispersion are also found for different metallicity subpopulations. Our results suggest that the Milky Way disc merger history can be recovered by relating the observed chemo-kinematic relations to the properties of past merger events.

Convergent migration allows pairs of planet to become trapped into mean motion resonances. Once in resonance, the planets' eccentricities grow to an equilibrium value that depends on the ratio of migration time scale to the eccentricity damping timescale, $K=\tau_a/\tau_e$, with higher values of equilibrium eccentricity for lower values of $K$. For low equilibrium eccentricities, $e_{eq}\propto K^{-1/2}$. The stability of a planet pair depends on eccentricity so the system can become unstable before it reaches its equilibrium eccentricity. Using a resonant overlap criterion that takes into account the role of first and second order resonances and depends on eccentricity, we find a function $K_{min}(\mu_p, j)$ that defines the lowest value for $K$, as a function of the ratio of total planet mass to stellar mass ($\mu_p$) and the period ratio of the resonance defined as $P_1/P_2=j/(j+k)$, that allows two convergently migrating planets to remain stable in resonance at their equilibrium eccentricities. We scaled the functions $K_{min}$ for each resonance of the same order into a single function $K_c$. The function $K_{c}$ for planet pairs in first order resonances is linear with increasing planet mass and quadratic for pairs in second order resonances with a coefficient depending on the relative migration rate and strongly on the planet to planet mass ratio. The linear relation continues until the mass approaches a critical mass defined by the 2/7 resonance overlap instability law and $K_c \to \infty$. We compared our analytic boundary with an observed sample of resonant two planet systems. All but one of the first order resonant planet pair systems found by radial velocity measurements are well inside the stability region estimated by this model. We calculated $K_c$ for Kepler systems without well-constrained eccentricities and found only weak constraints on $K$.

Interfering two- and three-arm spiral patterns have previously been inferred to exist in many galaxies and also in numerical simulations, and invoked to explain important dynamical properties, such as lack of symmetry, kinks in spiral arms, and star formation in armlets. The non-axisymmetric galactic mean-field dynamo model of Chamandy et al. is generalized to allow for such multiple co-existing spiral patterns in the kinetic alpha_k effect, leading to the existence of magnetic spiral arms in the large-scale magnetic field with several new properties. The large-scale magnetic field produced by an evolving superposition of two- and three-arm (or two- and four-arm) patterns evolves with time along with the superposition. Magnetic arms can be stronger and more extended in radius and in azimuth when produced by two interfering patterns rather than by one pattern acting alone. Transient morphological features arise in the magnetic arms, including bifurcations, disconnected armlets, and temporal and spatial variation in arm strength and pitch angles. Pitch angles of the large-scale magnetic field and magnetic arm structures (ridges) are smaller than those typically inferred from observations of spiral galaxies for model parameters of Chamandy et al., but can become comparable to typically inferred values for certain (still realistic) parameters. The magnetic field is sometimes strongest in between the alpha_k-arms, unlike in standard models with a single pattern, where it is strongest within the alpha_k-arms. Moreover, for models with a two- and three-arm pattern, some amount of m=1 azimuthal symmetry is found to be present in the magnetic field, which is generally not the case for forcing by single two- or three-arm patterns. (abridged)

We explore a second order Hamiltonian vertical resonance model for X-shaped or peanut-shaped galactic bulges. The X-shape is caused by the 2:1 vertical Lindblad resonance with the bar, with two vertical oscillation periods per orbital period in the bar frame. We examine N-body simulations and find that due to the bar slowing down and disk thickening during bar buckling, the resonance and associated peanut-shape moves outward. The peanut-shape is consistent with the location of the vertical resonance, independent of whether the bar buckled or not. We estimate the resonance width from the potential m=4 Fourier component and find that the resonance is narrow, affecting orbits over a narrow range in the angular momentum distribution, dL/L ~ 0.05. As the resonance moves outward, stars originally in the mid plane are forced out of the mid plane into orbits just within the resonance separatrix. The height of the separatrix orbits, estimated from the Hamiltonian model, is approximately consistent with the peanut-shape height. The X-shape is comprised of stars in the vicinity of the resonance separatrix. The velocity distributions from the simulations illustrate that low inclination orbits are depleted within resonance. Within resonance, the vertical velocity distribution is broad, consistent with resonant heating caused by the passage of the resonance through the disk. In the Milky Way bulge we relate the azimuthally averaged mid-plane mass density near the vertical resonance to the rotation curve and bar pattern speed. At an estimated vertical resonance galactocentric radius of ~1.3 kpc, we confirm a mid-plane density of ~5x10^8 Msol/kpc^3, consistent with recently estimated mass distributions. We find that the rotation curve, bar pattern speed, 2:1 vertical resonance location, X-shape tips, and mid-plane mass density, are all self-consistent in the Milky Way galaxy bulge.

We explore scenarios for the origin of two different density planets in the Kepler 36 system in adjacent orbits near the 7:6 mean motion resonance. We find that fine tuning is required in the stochastic forcing amplitude, the migration rate and planet eccentricities to allow two convergently migrating planets to bypass mean motion resonances such as the 4:3, 5:4 and 6:5, and yet allow capture into the 7:6 resonance. Stochastic forcing can eject the system from resonance causing a collision between the planets, unless the disk inducing migration and stochastic forcing is depleted soon after resonance capture. We explore a scenario with approximately Mars mass embryos originating exterior to the two planets and migrating inwards toward two planets. We find that gravitational interactions with embryos can nudge the system out of resonances. Numerical integrations with about a half dozen embryos can leave the two planets in the 7:6 resonance. Collisions between planets and embryos have a wide distribution of impact angles and velocities ranging from accretionary to disruptive. We find that impacts can occur at sufficiently high impact angle and velocity that the envelope of a planet could have been stripped, leaving behind a dense core. Some of our integrations show the two planets exchanging locations, allowing the outer planet that had experienced multiple collisions with embryos to become the innermost planet. A scenario involving gravitational interactions and collisions with embryos may account for both the proximity of the Kepler 36 planets and their large density contrast.

HR 8799 is a four planet system that also hosts a debris disk. By numerically integrating both planets and a planetesimal disk, we find interactions between an exterior planetesimal disk and the planets can influence the lifetime of the system. We first consider resonant planetary configurations that remained stable for at least 7 Myrs sans debris disk. An exterior debris disk with only 1 per cent the mass of the outermost planet (approximately a Neptune mass) was sufficiently large enough to pull the system out of resonance after 2 to 6 Myrs. Secondly, we consider configurations which are unstable in less than a few hundred thousand years. We find that these can be stabilized by a debris disk with a mass of more than 10 per cent that of the outermost planet. Our two sets of simulations suggest that estimates of the long term stability of a planetary system should take into account the role of the debris disk.

Non-axisymmetric components, such as spirals and central bars, play a major role in shaping galactic discs. An important aspect of the disc secular evolution driven by these perturbers is the radial migration of stars. It has been suggested recently that migration can populate a thick-disc component from inner-disc stars with high vertical energies. Since this has never been demonstrated in simulations, we study in detail the effect of radial migration on the disc velocity dispersion and disc thickness, by separating simulated stars into migrators and non-migrators. We apply this method to three isolated barred Tree-SPH N-body galaxies with strong radial migration. Contrary to expectations, we find that as stellar samples migrate, on the average, their velocity dispersion change (by as much as 50%) in such a way as to approximately match the non-migrating population at the radius at which they arrive. We show that, in fact, migrators suppress heating in parts of the disc. To confirm the validity of our findings, we also apply our technique to three cosmological re-simulations, which use a completely different simulation scheme and, remarkably, find very similar results. We believe the inability of migration to thicken discs is a fundamental property of internal disc evolution, irrespective of the migration mechanism at work. We explain this with the approximate conservation of the (average) vertical and radial actions rather than the energy. This "action mixing" can be used to constrain the migration rate in the Milky Way: estimates of the average vertical action in observations for different populations of stars should reveal flattening with radius for older groups of stars.

We investigate the evolution of galactic disks in N-body Tree-SPH simulations. We find that disks, initially truncated at three scale-lengths, can triple their radial extent, solely driven by secular evolution. Both Type I (single exponential) and Type II (down-turning) observed disk surface-brightness profiles can be explained by our findings. We relate these results to the strong angular momentum outward transfer, resulting from torques and radial migration associated with multiple patterns, such as central bars and spiral waves of different multiplicity. We show that even for stars ending up on cold orbits, the changes in angular momentum exhibit complex structure as a function of radius, unlike the expected effect of transient spirals alone. Focussing on one of our models, we find evidence for non-linear coupling among m=1, 2, 3 and 4 density waves, where m is the pattern multiplicity. We suggest that the naturally occurring larger resonance widths at galactic radii beyond four scale-lengths may have profound consequences on the formation and location of breaks in disk density profiles, provided spirals are present at such large distances. We also consider the effect of gas inflow and show that when in-plane smooth gas accretion of ~5 M_sun/yr is included, the outer disks become more unstable, leading to a strong increase in the stellar velocity dispersion. This, in turn, causes the formation of a Type III (up-turning) profile in the old stellar population. We propose that observations of Type III surface brightness profiles, combined with an up-turn in the stellar velocity dispersions beyond the disk break, could be a signature of ongoing gas-accretion. The results of this study suggest that disk outskirts comprised of stars migrated from the inner disk would have relatively large radial velocity dispersions, and significant thickness when seen edge-on. [Abridged]

A fraction of multiple planet candidate systems discovered from transits by the Kepler mission contain pairs of planet candidates that are in orbital resonance or are spaced slightly too far apart to be in resonance. We focus here on the four planet system, KOI 730, that has planet periods satisfying the ratios 8:6:4:3. By numerically integrating four planets initially in this resonant configuration in proximity to an initially exterior cold planetesimal disk, we find that of the order of a Mars mass of planet-orbit-crossing planetesimals is sufficient to pull this system out of resonance. Approximately one Earth mass of planet-orbit-crossing planetesimals increases the interplanetary spacings sufficiently to resemble the multiple planet candidate Kepler systems that lie just outside of resonance. This suggests that the closely spaced multiple planet Kepler systems, host only low mass debris disks or their debris disks have been extremely stable. We find that the planetary inclinations increase as a function of the mass in planetesimals that have crossed the orbits of the planets. If systems are left at zero inclination and in resonant chains after depletion of the gas disk then we would expect a correlation between distance to resonance and mutual planetary inclinations. This may make it possible to differentiate between dynamical mechanisms that account for the fraction of multiple planet systems just outside of resonance.

We identify migrating stars in an N-body simulation of a Milky-Way-like disk. Outward migration can occur when a star in a low eccentricity orbit lags a short-lived local spiral arm density peak. We interpret short lived local density peaks, that appear and fade on approximately an orbital period, as arising from positive interference between spiral density wave patterns. Stars near such a peak can migrate over a significant distance in galactocentric radius during the peak lifetime, providing that the peak is sufficiently dense. Using a Gaussian bar model for the potential perturbation associated with a narrow transient spiral feature, estimates of the migration rate, angular offset between particle and spiral feature, and maximum eccentricity for migrators roughly agrees with the values measured in our simulation. When multiple spiral density waves are present, local density peaks can appear and disappear on timescales faster than the timescale estimated for growth and decay of individual waves and the peak surface density can be larger than for any individual wave. Consequently, migration induced by transient density peaks may be more pervasive than that mediated by the growth and decay of individual patterns and occurring at their corotation resonance. We discuss interpretation of transient-like behavior in terms of interfering patterns, including estimating a coherence time for features that appear due to constructive interference, their effective angular rotation rates and the speed and direction that a density maximum would move across a galaxy inducing a localized and traveling burst of star formation.

New Chandra X-ray and Herschel FIR observations enable a multiwavelength study of active galactic nucleus (AGN) heating and intracluster medium (ICM) cooling in the brightest cluster galaxy of Abell 2597. The new Chandra observations reveal the central < 30 kiloparsec X-ray cavity network to be more extensive than previously thought, and associated with enough enthalpy to theoretically inhibit the inferred classical cooling flow. Nevertheless, we present new evidence, consistent with previous results, that a moderately strong residual cooling flow is persisting at 4%-8% of the classically predicted rates in a spatially structured manner amid the feedback-driven excavation of the X-ray cavity network. New Herschel observations are used to estimate warm and cold dust masses, a lower-limit gas-to-dust ratio, and a star formation rate consistent with previous measurements. The cooling time profile of the ambient X-ray atmosphere is used to map the locations of the observational star formation entropy threshold as well as the theoretical thermal instability threshold. Both lie just outside the < 30 kpc central region permeated by X-ray cavities, and star formation as well as ionized and molecular gas lie interior to both. The young stars are distributed in an elongated region that is aligned with the radio lobes, and their estimated ages are both younger and older than the X-ray cavity network, suggesting both jet-triggered as well as persistent star formation over the current AGN feedback episode. Bright X-ray knots that are coincident with extended Ly-alpha and FUV continuum filaments motivate a discussion of structured cooling from the ambient hot atmosphere along a projected axis that is perpendicular to X-ray cavity and radio axis. We conclude that the cooling ICM is the dominant contributor of the cold gas reservoir fueling star formation and AGN activity in the Abell 2597 BCG.

We present new Chandra X-ray observations of the brightest cluster galaxy (BCG) in the cool core cluster Abell 2597. The data reveal an extensive kpc-scale X-ray cavity network as well as a 15 kpc filament of soft-excess gas exhibiting strong spatial correlation with archival VLA radio data. In addition to several possible scenarios, multiwavelength evidence may suggest that the filament is associated with multiphase (10^3 - 10^7 K) gas that has been entrained and dredged-up by the propagating radio source. Stemming from a full spectral analysis, we also present profiles and 2D spectral maps of modeled X-ray temperature, entropy, pressure, and metal abundance. The maps reveal an arc of hot gas which in projection borders the inner edge of a large X-ray cavity. Although limited by strong caveats, we suggest that the hot arc may be (a) due to a compressed rim of cold gas pushed outward by the radio bubble or (b) morphologically and energetically consistent with cavity-driven active galactic nucleus (AGN) heating models invoked to quench cooling flows, in which the enthalpy of a buoyant X-ray cavity is locally thermalized as ambient gas rushes to refill its wake. If confirmed, this would be the first observational evidence for this model.

<abridged> The large relative sizes of circumstellar and circumplanetary disks imply that they might be seen in eclipse in stellar light curves. We estimate that a survey of ~10^4 young (~10 Myr old) post-accretion pre-MS stars monitored for ~10 years should yield at least a few deep eclipses from circumplanetary disks and disks surrounding low mass companion stars. We present photometric and spectroscopic data for a pre-MS K5 star (1SWASP J140747.93-394542.6), a newly discovered ~0.9 Msun member of the ~16 Myr-old Upper Cen-Lup subgroup of Sco-Cen at a kinematic distance of 128 pc. SuperWASP and ASAS light curves for this star show a remarkably long, deep, and complex eclipse event centered on 29 April 2007. At least 5 multi-day dimming events of >0.5 mag are identified, with a >3.3 mag deep eclipse bracketed by two pairs of ~1 mag eclipses symmetrically occurring +-12 days and +-26 days before and after. Hence, significant dimming of the star was taking place on and off over at least a ~54 day period in 2007, and a strong >1 mag dimming event occurred over a ~12 day span. We place a firm lower limit on the period of 850 days (i.e. the orbital radius of the eclipser must be >1.7 AU and orbital velocity must be <22 km/s). The shape of the light curve is similar to the lop-sided eclipses of the Be star EE Cep. We suspect that this new star is being eclipsed by a low-mass object orbited by a dense inner disk, girded by at least 3 dusty rings of lower optical depth. Between these rings are at least two annuli of near-zero optical depth (i.e. gaps), possibly cleared out by planets or moons, depending on the nature of the secondary. For possible periods in the range 2.33-200 yr, the estimated total ring mass is ~8-0.4 Mmoon (if the rings have optical opacity similar to Saturn's rings), and the edge of the outermost detected ring has orbital radius ~0.4-0.09 AU.

By logging encounters between planetesimals and planets we compute the distribution of encounters in a numerically integrated two planet system that is migrating due to interactions with an exterior planetesimal belt. Capture of an irregular satellite in orbit about a planet through an exchange reaction with a binary planetesimal is only likely when the binary planetesimal undergoes a slow and close encounter with the planet. In our simulations we find that close and slow encounters between planetesimals and a planet primarily occur with the outermost and not innermost planet. Taking care to consider where a planet orbit crossing binary planetesimal would first be tidally disrupted, we estimate the probability of both tidal disruption and irregular satellite capture. We estimate that the probability that the secondary of a binary planetesimal is captured and becomes an irregular satellite about a Neptune mass outer planet is about 1/100 for binaries with masses and separations similar to transneptunian planetesimal binaries. If young exoplanetary debris disks host a binary planetesimal population then outwards migrating outer planets should host captured irregular satellite populations. We discuss interpretation of emission associated with the exoplanet Fomalhaut b in terms of collisional evolution of a captured irregular satellite population that is replenished due to planetary migration.

Disc non-axisymmetrc components, such as spirals and central bars, are nowadays known to play an important role in shaping galactic discs. Here we use Tree-SPH N-body simulations to examine the effect of these perturbers on two aspects: the occurrence of multiple patterns in discs and the effects of radial migration on disc thickening. We find that, in addition to a central bar, multiple spiral patterns and lopsided modes develop in all models. Interaction among these asymmetric features results in a large scale stellar migration. However, we show that, despite the strong radial mixing, discs cannot be thickened sufficiently to match observed thick discs. We relate this to the adiabatic cooling as stars migrate radially outwards. We also find that the bulge contribution to a thick-disc component for an Sa-type galaxy at ~2.5 disc scale-lengths is less than 1% and zero in the case of a Milky Way-like, Sb-type. Our findings cast doubt on the plausibility of thick disc formation via stellar radial migration.

We compute the strengths of zero-th order (in eccentricity) three-body resonances for a co-planar and low eccentricity multiple planet system. In a numerical integration we illustrate that slowly moving Laplace angles are matched by variations in semi-major axes among three bodies with the outer two bodies moving in the same direction and the inner one moving in the opposite direction, as would be expected from the two quantities that are conserved in the three-body resonance. A resonance overlap criterion is derived for the closely and equally spaced, equal mass system with three-body resonances overlapping when interplanetary separation is less than an order unity factor times the planet mass to the one quarter power. We find that three-body resonances are sufficiently dense to account for wander in semi-major axis seen in numerical integrations of closely spaced systems and they are likely the cause of instability of these systems. For interplanetary separations outside the overlap region, stability timescales significantly increase. Crudely estimated diffusion coefficients in eccentricity and semi-major axis depend on a high power of planet mass and interplanetary spacing. An exponential dependence previously fit to stability or crossing timescales is likely due to the limited range of parameters and times possible in integration and the strong power law dependence of the diffusion rates on these quantities.

An N-body hybrid simulation, integrating both massive and tracer particles, of a Galactic disk is used to study the stellar phase space distribution or velocity distributions in different local neighborhoods. Pattern speeds identified in Fourier spectrograms suggest that two-armed and three-armed spiral density waves, a bar and a lopsided motion are coupled in this simulation, with resonances of one pattern lying near resonances of other patterns. We construct radial and tangential (uv) velocity distributions from particles in different local neighborhoods. More than one clump is common in these local velocity distributions regardless of the position in the disk. Features in the velocity distribution observed at one galactic radius are also seen in nearby neighborhoods (at larger and smaller radii) but with shifted mean v values. This is expected if the v velocity component of a clump sets the mean orbital galactic radius of its stars. We find that gaps in the velocity distribution are associated with the radii of kinks or discontinuities in the spiral arms. These gaps also seem to be associated with Lindblad resonances with spiral density waves and so denote boundaries between different dominant patterns in the disk. We discuss implications for interpretations of the Milky Way disk based on local velocity distributions. Velocity distributions created from regions just outside the bar's Outer Lindblad resonance and with the bar oriented at 45 degrees from the Sun-Galactic center line more closely resemble that seen in the solar neighborhood (triangular in shape at lower uv and with a Hercules like stream) when there is a strong nearby spiral arm, consistent with the observed Centaurus Arm tangent, just interior to the solar neighborhood.

Using HST NICMOS 2 observations we have measured 1.6-micron near infrared nuclear luminosities of 100 3CR radio galaxies with z<0.3, by modeling and subtracting the extended emission from the host galaxy. We performed a multi-wavelength statistical analysis (including optical and radio data) of the properties of the nuclei following classification of the objects into FRI and FRII, and LIG (low-ionization galaxies), HIG (high-ionization galaxies) and BLO (broad-lined objects) using the radio morphology and optical spectra, respectively. The correlations among near infrared, optical, and radio nuclear luminosity support the idea that the near infrared nuclear emission of FRIs has a non-thermal origin. Despite the difference in radio morphology, the multi-wavelength properties of FRII LIG nuclei are statistically indistinguishable from those of FRIs, an indication of a common structure of the central engine. All BLOs show an unresolved near infrared nucleus and a large near infrared excess with respect to FRII LIGs and FRIs of equal radio core luminosity. This requires the presence of an additional (and dominant) component other than the non-thermal light. Considering the shape of their spectral energy distribution, we ascribe the origin of their near infrared light to hot circumnuclear dust. A near infrared excess is also found in HIGs, but their nuclei are substantially fainter than those of BLO. This result indicates that substantial obscuration along the line-of-sight to the nuclei is still present at 1.6 micron. Nonetheless, HIGs nuclei cannot simply be explained in terms of dust obscuration: a significant contribution from light reflected in a circumnuclear scattering region is needed to account for their multiwavelength properties.

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.

Quillen et al. and O'Dea et al. carried out a Spitzer study of a sample of 62 brightest cluster galaxies (BCGs) from the ROSAT Brightest Cluster Sample chosen based on their elevated H-alpha flux. We present Hubble Space Telescope Advanced Camera for Surveys (ACS) far ultraviolet (FUV) images of the Ly-alpha and continuum emission of the luminous emission-line nebulae in 7 BCGs found to have an infrared excess. We confirm that the BCGs are actively forming stars suggesting that the IR excess seen in these BCGs is indeed associated with star formation. The FUV continuum emission extends over a region of ~7-28 kpc (largest linear size) and even larger in Ly-alpha. The young stellar population required by the FUV observations would produce a significant fraction of the ionizing photons required to power the emission line nebulae. Star formation rates estimated from the FUV continuum range from ~3 to ~14 times lower than those estimated from the IR, however both the Balmer decrement in the central few arcseconds and detection of CO in most of these galaxies imply that there are regions of high extinction that could have absorbed much of the FUV continuum. Analysis of archival VLA observations reveals compact radio sources in all 7 BCGs and kpc scale jets in Abell 1835 and RXJ 2129+00. The four galaxies with archival deep Chandra observations exhibit asymmetric X-ray emission, the peaks of which are offset from the center of the BCGs by ~10 kpc on average. A low feedback state for the AGN could allow increased condensation of the hot gas into the center of the galaxy and the feeding of star formation.

We use a hybrid test particle/N-body simulation to integrate 4 million massless test particle trajectories within a fully self-consistent 10^5 particle N-body simulation. The number of massless particles allows us to resolve fine structure in the spatial distribution and phase space of a dwarf galaxy as it is disrupted in the tidal field of a Milky Way type galaxy. The tidal tails exhibit nearly periodic clumping or a smoke-like appearance. By running simulations with different satellite particle mass, halo particle mass, number of massive and massless particles and with and without a galaxy disk, we have determined that the instabilities are not due to numerical noise, amplification of structure in the halo, or shocking as the satellite passes through the disk of the Galaxy. We measure Jeans wavelengths and growth timescales in the tidal tail and show that the Jeans instability is a viable explanation for the clumps. We find that the instability causes velocity perturbations of order 10 km/s. Clumps in tidal tails present in the Milky Way could be seen in stellar radial velocity surveys as well as number counts. We find that the unstable wavelength growth is sensitive to the simulated mass of dark matter halo particles. A simulation with a smoother halo exhibits colder and thinner tidal tails with more closely spaced clumps than a simulation with more massive dark matter halo particles. Heating by the halo particles increases the Jeans wavelength in the tidal tail affecting substructure development, suggesting an intricate connection between tidal tails and dark matter halo substructure.

The dust destruction timescales in the cores of clusters of galaxies are relatively short given their high central gas densities. However, substantial mid-infrared and sub-mm emission has been detected in many brightest cluster galaxies. In this letter we present Herschel PACS and SPIRE photometry of the brightest cluster galaxy in three strong cooling flow clusters, A1068, A2597 and Zw3146. This photometry indicates that a substantial mass of cold dust is present (>3 x 10^7 Mo) at temperatures significantly lower (20-28K) than previously thought based on limited MIR and/or sub-mm results. The mass and temperature of the dust appear to match those of the cold gas traced by CO with a gas-to-dust ratio of 80-120.

The question of how much gas cools in the cores of clusters of galaxies has been the focus of many, multiwavelength studies in the past 30 years. In this letter we present the first detections of the strongest atomic cooling lines, [C II], [O I] and [N I] in two strong cooling flow clusters, A1068 and A2597, using Herschel PACS. These spectra indicate that the substantial mass of cold molecular gas (>10^9 Mo) known to be present in these systems is being irradiated by intense UV radiation, most probably from young stars. The line widths of these FIR lines indicate that they share dynamics similar but not identical to other ionised and molecular gas traced by optical, near-infrared and CO lines. The relative brightness of the FIR lines compared to CO and FIR luminosity is consistent with other star-forming galaxies indicating that the properties of the molecular gas clouds in cluster cores and the stars they form are not unusual. These results provide additional evidence for a reservoir of cold gas that is fed by the cooling of gas in the cores of the most compact clusters and provide important diagnostics of the temperature and density of the dense clouds this gas resides in.

Tidal tails composed of stars should be unstable to the Jeans instability and this can cause them to look like beads on a string. The Jeans wavelength and tail diameter determine the wavelength and growth rate of the fastest growing unstable mode. Consequently the distance along the tail to the first clump and spacing between clumps can be used to estimate the mass density in the tail and its longitudinal velocity dispersion. Clumps in the tidal tails of the globular cluster Palomar 5 could be due to Jeans instability. We find that their spacing is consistent with the fastest growing mode if the velocity dispersion in the tail is similar to that in the cluster itself. While all tidal tails should exhibit gravitational instability, we find that clusters or galaxies with low concentration parameters are most likely to exhibit short wavelength rapidly growing Jeans modes in their tidal tails.

Previous work has related the Galactic bar to structure in the local stellar velocity distribution. Here we show that the bar also influences the spatial gradients of the velocity vector via the Oort constants. By numerical integration of test-particles we simulate measurements of the Oort C value in a gravitational potential including the Galactic bar. We account for the observed trend that C is increasingly negative for stars with higher velocity dispersion. By comparing measurements of C with our simulations we improve on previous models of the bar, estimating that the bar pattern speed is Omega_b/Omega_0=1.87\pm0.04, where Omega_0 is the local circular frequency, and the bar angle lies within 20[deg]<phi_0<45[deg]. We find that the Galactic bar affects measurements of the Oort constants A and B less than ~2 km/s/kpc for the hot stars.

We compile position and inclination angles for tilted ring fits to the warped dusty and gaseous disk spanning radius 1.8 to 6500 pc from recent observations. For radii exterior to a kpc, tilted ring orientations lie on an arc on a polar inclination versus position angle plot, suggesting that precession following a merger can account for the ring morphology. Three kinks in the ring orientations are seen on the polar plot, the one at radius of about 1.3 kpc we suspect corresponds to the location where self-gravity in the disk affects the ring precession rate. Another at a radius of about 600 pc may be associated with a gap in the gas distribution. A third kink is seen at a radius of 100 pc. A constant inclination tilted disk precessing about the jet axis may describe the disk between 100 and 20 pc but not interior to this. A model with disk orientation matching the molecular circumnuclear disk at 100 pc that decays at smaller radii to an inner flat disk perpendicular to the jet may account for disk orientations within 100 pc. Neither model would account for the cusps or changes in disk orientation at 100 or 600 pc.

We present 19 nearby (z<0.3) 3CR radio galaxies imaged at low- and high-excitation as part of a Cycle 15 Hubble Space Telescope snapshot survey with the Advanced Camera for Surveys. These images consist of exposures of the H-alpha (6563 \AA, plus [NII] contamination) and [OIII] 5007 \AA emission lines using narrow-band linear ramp filters adjusted according to the redshift of the target. To facilitate continuum subtraction, a single-pointing 60 s line-free exposure was taken with a medium-band filter appropriate for the target's redshift. We discuss the steps taken to reduce these images independently of the automated recalibration pipeline so as to use more recent ACS flat-field data as well as to better reject cosmic rays. We describe the method used to produce continuum-free (pure line-emission) images, and present these images along with qualitative descriptions of the narrow-line region morphologies we observe. We present H-alpha+[NII] and [OIII] line fluxes from aperture photometry, finding the values to fall expectedly on the redshift-luminosity trend from a past HST/WFPC2 emission line study of a larger, generally higher redshift subset of the 3CR. We also find expected trends between emission line luminosity and total radio power, as well as a positive correlation between the size of the emission line region and redshift. We discuss the associated interpretation of these results, and conclude with a summary of future work enabled by this dataset.

Using test particle simulations we examine the structure of the outer Galactic disk as it is perturbed by a satellite in a tight eccentric orbit about the Galaxy. A satellite of mass a few times 10^9 Msol can heat the outer Galactic disk, excite spiral structure and a warp and induce streams in the velocity distribution. We examine particle eccentricity versus the change in mean radius between initial and current orbits. Correlations between these quantities are reduced after a few satellite pericenter passages. Stars born in the outer galaxy can be moved in radius from their birth positions and be placed in low eccentricity orbits inside their birth radii. We propose that mergers and perturbations from satellite galaxies and subhalos can induce radial mixing in the stellar metallicity distribution.

We perform numerical simulations of a stellar galactic disk with initial conditions chosen to represent an unrelaxed population which might have been left following a merger. Stars are unevenly distributed in radial action angle, though the disk is axisymmetric. The velocity distribution in the simulated Solar neighborhood exhibits waves traveling in the direction of positive v, where u,v are the radial and tangential velocity components. As the system relaxes and structure wraps in phase space, the features seen in the uv-plane move closer together. We show that these results can be obtained also by a semi-analytical method. We propose that this model could provide an explanation for the high velocity streams seen in the Solar neighborhood at approximate v in km/s, of -60 (HR 1614), -80 (Arifyanto and Fuchs 2006), -100 (Arcturus), and -160 (Klement et al. 2008). In addition, we predict four new features at v ~ -140, -120, 40 and 60 km/s. By matching the number and positions of the observed streams, we estimate that the Milky Way disk was strongly perturbed ~1.9 Gyr ago. This event could have been associated with Galactic bar formation.

Abridged: The golden age of astrophysics is upon us with both grand discoveries (extra-solar planets, dark matter, dark energy) and precision cosmology. Fundamental understanding of the working of stars and galaxies is within reach, thanks to newly available precision measurements. We highlight the importance of distances and model independent distances and masses. Distances are fundamental in astrophysics and their knowledge can change our perception of phenomena dramatically: e.g., in antiquity, the Heliocentric model was rejected because the predicted stellar parallaxes were not observed. Distance measurements are directly related to the history & fate of the uni- verse as they provide 2 of 3 methods available to date the universe. The 1st method is based on the ages of stars, which can be ascertained if their lumi- nosities (distances) are accurately known. The 2nd method relies on cosmolo- gical methods. To 1st order, the age of the universe is the inverse H_0. As stressed by the previous decadal report, "the fundamental goal of ... astrophysics is to understand how the universe ... galaxies [and] stars ... formed, how they evolved, and what their destiny will be." These questions can be answered partly by micro-arcsecond astrometry: 1) Galactic archeology: a detailed reconstruction of the formation history of the Milky Way and other Local Group galaxies, 2) the very oldest stars in the Milky Way and the age of the Universe, and 3) H_0 and concordance cosmology. These goals are achievable by combining muas-arcsecond astrometry from the proposed SIM-Lite mission supplemented with ground-based spectroscopy. The results of our proposed project will force the biggest reassessment of stellar astrophysics in 50 years, which will affect most branches of astrophysics.

Radial velocity searches for extrasolar planets have discovered stars undergoing slow accelerations. The accelerations are likely due to one or more gas giant planets or brown dwarfs orbiting with period longer than the total time span of observations. The stellar acceleration is proportional to the mass of the companion divided by the square of its radius from the star. In this paper we predict the distribution of accelerations using a Monte Carlo method and assuming a semi-major axis and mass distribution for the companions. Radial velocity surveys find that ~20 and 10% of stars surveyed exhibit accelerations above ~10 and 25 m s^-1 yr^-1, respectively. We find that an extrapolation of the size and period distribution found by radial velocity surveys within ~5 AU or companion imaging surveys predicts too few systems with detectable accelerations. The fractions of stars with accelerations above these values can be matched if every star hosts a gas giant planet or brown dwarf with semi-major axes between 4 and 20 AU or if a significant fraction (~70%) of stars host 3 or more massive planets in this range. If brown dwarfs are not discovered in imaging surveys then we would infer that there is a large population of gas giant planets present with semi-major axes in the 5 to 30 AU range.

We simulate planet migration caused by interactions between planets and a planetesimal disk. We use an N-body integrator optimized for near-Keplerian motion that runs in parallel on a video graphics card, and that computes all pair-wise gravitational interactions. We find that the fraction of planetesimals found in mean motion resonances is reduced and planetary migration rates are on average about 50% slower when gravitational interactions between the planetesimals are computed than when planetesimal self-gravity is neglected. This is likely due to gravitational stirring of the planetesimal disk that is not present when self-gravity is neglected that reduces their capture efficiency because of the increased particle eccentricity dispersion. We find that migration is more stochastic when the disk is self-gravitating or comprised of more massive bodies. Previous studies have found that if the planetesimal disk density is below a critical level, migration is "damped" and the migration rate decays exponentially, otherwise it is "forced" and the planet's migration rate could accelerate exponentially. Migration rates measured from our undamped simulations suggest that the migration rate saturates at a level proportional to disk density and subsequently is approximately power law in form with time.

The morphology of the outer rings of early-type spiral galaxies is compared to integrations of massless collisionless particles initially in nearly circular orbits. Particles are perturbed by a quadrupolar gravitational potential corresponding to a growing and secularly evolving bar. We find that outer rings with R1R2 morphology and pseudorings are exhibited by the simulations even though they lack gaseous dissipation. Simulations with stronger bars form pseudorings earlier and more quickly than those with weaker bars. We find that the R1 ring, perpendicular to the bar, is fragile and dissolves after a few bar rotation periods if the bar pattern speed increases by more than ~ 8%, bar strength increases (by >~ 140%) after bar growth, or the bar is too strong (Q_T>0.3). If the bar slows down after formation, pseudoring morphology persists and the R2 ring perpendicular to the bar is populated due to resonance capture. The R2 ring remains misaligned with the bar and increases in ellipticity as the bar slows down. The R2 ring becomes scalloped and does not resemble any ringed galaxies if the bar slows down more than 3.5% suggesting that bars decrease in strength before they slow down this much. We compare the morphology of our simulations to B-band images of 9 ringed galaxies from the Ohio State University Bright Spiral Galaxy Survey, and we find a reasonable match in morphologies to R1R2' pseudorings seen within a few bar rotation periods of bar formation. Some of the features previously interpreted in terms of dissipative models may be due to transient structure associated with recent bar growth and evolution.

We present the second part of an H-band (1.6 microns) atlas of z<0.3 3CR radio galaxies, using the Hubble Space Telescope Near Infrared Camera and Multi-Object Spectrometer (HST NICMOS2). We present new imaging for 21 recently acquired sources, and host galaxy modeling for the full sample of 101 (including 11 archival) -- an 87% completion rate. Two different modeling techniques are applied, following those adopted by the galaxy morphology and the quasar host galaxy communities. Results are compared, and found to be in excellent agreement, although the former breaks down in the case of strongly nucleated sources. Companion sources are tabulated, and the presence of mergers, tidal features, dust disks and jets are catalogued. The tables form a catalogue for those interested in the structural and morphological dust-free host galaxy properties of the 3CR sample, and for comparison with morphological studies of quiescent galaxies and quasar host galaxies. Host galaxy masses are estimated, and found to typically lie at around 2*10^11 solar masses. In general, the population is found to be consistent with the local population of quiescent elliptical galaxies, but with a longer tail to low Sersic index, mainly consisting of low-redshift (z<0.1) and low-radio-power (FR I) sources. A few unusually disky FR II host galaxies are picked out for further discussion. Nearby external sources are identified in the majority of our images, many of which we argue are likely to be companion galaxies or merger remnants. The reduced NICMOS data are now publicly available from our website (http://archive.stsci.edu/prepds/3cr/)

We present the second data release of the Radial Velocity Experiment (RAVE), an ambitious spectroscopic survey to measure radial velocities (RVs) and stellar atmosphere parameters of up to one million stars using the 6dF multi-object spectrograph on the 1.2-m UK Schmidt Telescope of the Anglo-Australian Observatory (AAO). It is obtaining medium resolution spectra (median R=7,500) in the Ca-triplet region (8,410--8,795 \AA) for southern hemisphere stars in the magnitude range 9<I<12. Following the first data release (Steinmetz et al. 2006) the current release doubles the sample of published RVs, now containing 51,829 RVs for 49,327 individual stars observed on 141 nights between April 11 2003 and March 31 2005. Comparison with external data sets shows that the new data collected since April 3 2004 show a standard deviation of 1.3 km/s, about twice better than for the first data release. For the first time this data release contains values of stellar parameters from 22,407 spectra of 21,121 individual stars. They were derived by a penalized \chi^2 method using an extensive grid of synthetic spectra calculated from the latest version of Kurucz models. From comparison with external data sets, our conservative estimates of errors of the stellar parameters (for a spectrum with S/N=40) are 400 K in temperature, 0.5 dex in gravity, and 0.2 dex in metallicity. We note however that the internal errors estimated from repeat RAVE observations of 822 stars are at least a factor 2 smaller. We demonstrate that the results show no systematic offsets if compared to values derived from photometry or complementary spectroscopic analyses. The data release includes proper motion and photometric measurements. It can be accessed via the RAVE webpage: http://www.rave-survey.org and through CDS.

The link between turbulence in star formatting environments and protostellar jets remains controversial. To explore issues of turbulence and fossil cavities driven by young stellar outflows we present a series of numerical simulations tracking the evolution of transient protostellar jets driven into a turbulent medium. Our simulations show both the effect of turbulence on outflow structures and, conversely, the effect of outflows on the ambient turbulence. We demonstrate how turbulence will lead to strong modifications in jet morphology. More importantly, we demonstrate that individual transient outflows have the capacity to re-energize decaying turbulence. Our simulations support a scenario in which the directed energy/momentum associated with cavities is randomized as the cavities are disrupted by dynamical instabilities seeded by the ambient turbulence. Consideration of the energy power spectra of the simulations reveals that the disruption of the cavities powers an energy cascade consistent with Burgers'-type turbulence and produces a driving scale-length associated with the cavity propagation length. We conclude that fossil cavities interacting either with a turbulent medium or with other cavities have the capacity to sustain or create turbulent flows in star forming environments. In the last section we contrast our work and its conclusions with previous studies which claim that jets can not be the source of turbulence.

In this paper we explore the relationship between protostellar outflows and turbulence in molecular clouds. Using 3-D numerical simulations we focus on the hydrodynamics of multiple outflows interacting within a parsec scale volume. We explore the extent to which transient outflows injecting directed energy and momentum into a sub-volume of a molecular cloud can be converted into random turbulent motions. We show that turbulence can readily be sustained by these interactions and show that it is possible to broadly characterize an effective driving scale of the outflows. We compare the velocity spectrum obtained in our studies to that of isotropically forced hydrodynamic turbulence finding that in outflow driven turbulence a power law is indeed achieved. However we find a steeper spectrum (beta ~ 3) is obtained in outflow driven turbulence models than in isotropically forced simulations (beta ~ 2). We discuss possible physical mechanisms responsible for these results as well and their implications for turbulence in molecular clouds where outflows will act in concert with other processes such as gravitational collapse.

Quillen et al.(2007) presented an imaging survey with the {\it Spitzer Space Telescope} of 62 brightest cluster galaxies with optical line emission located in the cores of X-ray luminous clusters. They found that at least half of these sources have signs of excess infrared emission. Here we discuss the nature of the IR emission and its implications for cool core clusters. The strength of the mid-IR excess emission correlates with the luminosity of the optical emission lines. Excluding the four systems dominated by an AGN, the excess mid-infrared emission in the remaining brightest cluster galaxies is likely related to star formation. The mass of molecular gas (estimated from CO observations) is correlated with the IR luminosity as found for normal star forming galaxies. The gas depletion time scale is about 1 Gyr. The physical extent of the infrared excess is consistent with that of the optical emission line nebulae. This supports the hypothesis that the star formation occurs in molecular gas associated with the emission line nebulae and with evidence that the emission line nebulae are mainly powered by ongoing star formation. We find a correlation between mass deposition rates (${\dot M}_X$) estimated from the X-ray emission and the star formation rate estimated from the infrared luminosity. The star formation rates are 1/10 to 1/100 of the mass deposition rates suggesting that the re-heating of the ICM is generally very effective in reducing the amount of mass cooling from the hot phase but not eliminating it completely.

We introduce a differential equation for star formation in galaxies that incorporates negative feedback with a delay. When the feedback is instantaneous, solutions approach a self-limiting equilibrium state. When there is a delay, even though the feedback is negative, the solutions can exhibit cyclic and episodic solutions. We find that periodic or episodic star formation only occurs when two conditions are satisfied. Firstly the delay timescale must exceed a cloud consumption timescale. Secondly the feedback must be strong. This statement is quantitatively equivalent to requiring that the timescale to approach equilibrium be greater than approximately twice the cloud consumption timescale. The period of oscillations predicted is approximately 4 times the delay timescale. The amplitude of the oscillations increases with both feedback strength and delay time. We discuss applications of the delay differential equation (DDE) model to star formation in galaxies using the cloud density as a variable. The DDE model is most applicable to systems that recycle gas and only slowly remove gas from the system. We propose likely delay mechanisms based on the requirement that the delay time is related to the observationally estimated time between episodic events. The proposed delay timescale accounting for episodic star formation in galaxy centers on periods similar to P 10 Myrs, irregular galaxies with P 100 Myrs, and the Milky Way disk with P~ 2Gyr, could be that for exciting turbulence following creation of massive stars, that for gas pushed into the halo to return and interact with the disk and that for spiral density wave evolution, respectively.

Giant planets embedded in circumstellar discs are expected to open gaps in these discs. We examine the vertical structure of the gap edges. We find that the planet excites spiral arms with significant (Mach number of a half) vertical motion of the gas, and discuss the implications of these motions. In particular, the spiral arms will induce strong vertical stirring of the dust, making the edge appeared `puffed up' relative to the bulk of the disc. Infra-red observations (sensitive to dust) would be dominated by the light from the thick inner edge of the disc. Sub-millimetre observations (sensitive to gas velocities) would appear to be hot in `turbulent' motions (actually the ordered motion caused by the passage of the spiral arms), but cold in chemistry. Resolved sub-millimetre maps of circumstellar discs might even be able to detect the spiral arms directly.

Previous work has related the Galactic Bar to structure in the local stellar velocity distribution. Here we show that the Bar also influences the spatial gradients of the velocity vector via the Oort constants. By numerical integration of test-particles we simulate measurements of the Oort C value in a gravitational potential including the Galactic Bar. We account for the observed trend that C is increasingly negative for stars with higher velocity dispersion. By comparing measurements of C with our simulations we improve on previous models of the Bar, estimating that the Bar pattern speed is Omega_b/Omega_0=1.87\pm0.04, where Omega_0 is the local circular frequency, and the Bar angle lies within 20<phi_0<45 deg. We find that the Galactic Bar affects measurements of the Oort constants A and B less than ~2 km/s/kpc for the hot stars.

We report on an imaging survey with the Spitzer Space Telescope of 62 brightest cluster galaxies with optical line emission. These galaxies are located in the cores of X-ray luminous clusters selected from the ROSAT All-Sky Survey. We find that about half of these sources have a sign of excess infrared emission; 22 objects out of 62 are detected at 70 microns, 18 have 8 to 5.8 micron flux ratios above 1.0 and 28 have 24 to 8 micron flux ratios above 1.0. Altogether 35 of 62 objects in our survey exhibit at least one of these signs of infrared excess. Four galaxies with infrared excesses have a 4.5/3.6 micron flux ratio indicating the presence of hot dust, and/or an unresolved nucleus at 8 microns. Three of these have high measured [OIII](5007A)/Hbeta flux ratios suggesting that these four, Abell 1068, Abell 2146, and Zwicky 2089, and R0821+07, host dusty active galactic nuclei (AGNs). 9 objects (including the four hosting dusty AGNs) have infrared luminosities greater than 10^11 L_sol and so can be classified as luminous infrared galaxies (LIRGs). Excluding the four systems hosting dusty AGNs, the excess mid-infrared emission in the remaining brightest cluster galaxies is likely related to star formation.

We report on the results of infrared spectroscopic mapping observations carried out in the nuclear region of Centaurus A (NGC5128). The 500 pc bipolar dust shell discovered by Quillen et al.(2006) is even more clearly seen in the 11.3 micron dust emission feature than previous broad band imaging. The pure rotational lines of molecular hydrogen other than the S(0) line are detected above the dusty disk and associated with the oval dust shell. The molecular hydrogen transitions indicate the presence of warm gas at temperatures 250--720K. The ratio of the dust emission features at 7.7 and 11.3 micron and the ratio of the [NeII](12.8) and 11.3 dust emission feature are lower in the 500 pc dust shell than in the star forming disk. The clearer shell morphology at 11.3 micron, warm molecular hydrogen emission in the shell, and variation in line ratios in the shell compared to those in the disk, confirm spectroscopically that this shell is a separate coherent entity and is unlikely to be a chance superposition of dust filaments. The physical conditions in the shell are most similar to Galactic supernova remnants where blast waves encounter molecular clouds. The lines requiring the highest level of ionization, [NeV] and [OIV], are detected 20--25 arcsec north-east and south-west of the nucleus and at position angles near the radio jet axis. Fine structure line ratios and limits from this region suggest that the medium is low density and illuminated by a hard radiation field at low ionization parameter. These higher S molecular hydrogen pure rotational transitions are also particularly bright in the same region as the [OIV] and [NeV] emission. This suggests that the gas associated with the dust shell has been excited near the jet axis and is part of an ionization cone.

We investigate the effect of spiral structure on the Galactic disk as viewed by pencil beams centered on the Sun, relevant to upcoming surveys such as ARGOS, SEGUE, and GAIA. We create synthetic Galactic maps which we call Pencil Beam Maps (PBMs) of the following observables: line-of-sight velocities, the corresponding velocity dispersion, and the stellar number density that are functions of distance from the observer. We show that such maps can be used to infer spiral structure parameters, such as pattern speed, solar phase angle, and number of arms. The mean line-of-sight velocity and velocity dispersion are affected by up to ~35 km/s which is well within the detectable limit for forthcoming radial velocity surveys. One can measure the pattern speed by searching for imprints of resonances. In the case of a two-armed spiral structure it can be inferred from the radius of a high velocity dispersion ring situated at the 2:1 ILR. This information, however, must be combined with information related to the velocities and stellar number density in order to distinguish from a four-armed structure. If the pattern speed is such that the 2:1 ILR is hidden inside the Galactic bulge the 2:1 OLR will be present in the outer Galaxy and thus can equivalently be used to estimate the pattern speed. Once the pattern speed is known the solar angle can be estimated from the line-of-sight velocities and the number density PBMs. Forthcoming radial velocity surveys are likely to provide powerful constraints of the structure of the Milky Way disk.

We consider the minimum mass planet, as a function of radius, that is capable of opening a gap in an $\alpha$-accretion disk. We estimate that a half Jupiter mass planet can open a gap in a disk with accretion rate $\dot M \la 10^{-8} M_\odot$/yr for viscosity parameter $\alpha =0.01$, and Solar mass and luminosity. The minimum mass is approximately proportional to $\dot M^{0.48} \alpha^{0.8} M_*^{0.42} L_*^{-0.08}$. This estimate can be used to rule out the presence of massive planets in gapless accretion disks. We identify two radii at which an inwardly migrating planet may become able to open a gap and so slow its migration; the radius at which the heating from viscous dissipation is similar to that from stellar radiation in a flared disk, and the radius at which the disk has optical depth order 1 in a self-shadowed disk. If a gap opening planet cuts off disk accretion allowing the formation of a central hole or clearing in the disk then we would estimate that the clearing radius would approximately be proportional to the stellar mass.

We present computational results showing eccentricity growth in the inner portions of a protoplanetary disc. We attribute this to the evolving surface density of the disc. The planet creates a gap, which adjusts the balance between the 3:1 (eccentricity exciting) and 2:1 (eccentricity damping) resonances. The eccentricity of the inner disc can rise as high as 0.3, which is sufficient to cause it to be accreted onto the star. This offers an alternative mechanism for producing the large holes observed in the discs of CoKu Tau/4, GM Aur and DM Tau.

Infrared spectra from the Spitzer Space Telescope (SSC) of many debris disks are well fit with a single black body temperature which suggest clearings within the disk. We assume that inside the clearing orbital instability due to planets removes dust generating planetesimal belts and dust generated by the outer disk that is scattered or drifts into the clearing. From numerical integrations we estimate a minimum planet spacing required for orbital instability (and so planetesimal and dust removal) as a function of system age and planet mass. We estimate that a 10^8 year old debris disk with a dust disk edge at a radius of 50 AU hosted by an A star must contain approximately 5 Neptune mass planets between the clearing radius and the iceline in order to remove all primordial objects within it. We infer that known debris disk systems contain at least a fifth of a Jupiter mass in massive planets. The number of planets and spacing required is insensitive to the assumed planet mass. However an order of magnitude higher total mass in planets could reside in these systems if the planets are more massive.

We investigate isophotal properties and dust morphology in the nuclear regions of 84 radio galaxies, imaged in the optical and near-infrared as part of Hubble Space Telescope snapshot surveys. We present a sample-wide trend between host galaxy isophotal structure and the inclination of dusty circumnuclear disks at the centers of 13 of these objects. We find that galaxies containing edge-on disks are invariably seen to possess boxy isophotes, while round, face-on disks are seen exclusively in objects with round or elliptical isophotes. Dust-rich sources with disky isophotes are observed only to possess dust in the form of extended filamentary lanes, and not in settled distributions like disks. As we do not expect that edge-on and face-on disks reside in different populations of galaxies, we conclude that perceived isophotal boxiness is dependent upon the angle at which the observer views the host galaxy's axis of symmetry. We discuss our results in the context of dissipative merger scenarios, and infer that dusty disks primarily reside in old, boxy remnants of gas-poor galaxy mergers, whereas filamentary dust lanes reside in younger disky remnants of gas-rich mergers.

We consider constraints on the planetesimal population residing in the disks of AU Microscopii, Beta Pictoris and Fomalhaut taking into account their observed thicknesses and normal disk opacities. We estimate that bodies of radius 5, 180 and 70 km are responsible for initiating the collisional cascade accounting for the dust production for AU-Mic, Beta-Pic and Fomalhaut's disks, respectively, at break radii from the star where their surface brightness profiles change slope. Larger bodies, of radius 1000km and with surface density of order 0.01 g/cm^2, are required to explain the thickness of these disks assuming that they are heated by gravitational stirring. A comparison between the densities of the two sizes suggests the size distribution in the largest bodies is flatter than that observed in the Kuiper belt. AU Mic's disk requires the shallowest size distribution for bodies with radius greater than 10km suggesting that the disk contains planetary embryos experiencing a stage of runaway growth.

We perform test-particle simulations in a 2D, differentially rotating stellar disk, subjected to a two-armed steady state spiral density wave perturbation in order to estimate the influence of spiral structure on the local velocity field. By using Levenberg-Marquardt least-squares fit we decompose the local velocity field (as a result of our simulations) into Fourier components to fourth order. Thus we obtain simulated measurements of the Oort constants, A, B, C, and K. We get relations between the Fourier coefficients and some galactic parameters, such as the phase angle of the Solar neighborhood and the spiral pattern speed. We show that systematic errors due to the presence of spiral structure are likely to affect the measurements of the Oort constants. Moderate strength spiral structure causes errors of order 5 km/s/kpc in A and B. We find variations of the Fourier coefficients with velocity dispersion, pattern speed, and sample depth. For a sample at an average heliocentric distance of 0.8 kpc we can summarize our findings as follows:(i) if our location in the Galaxy is near corotation then we expect a vanishing value for C for all phase angles;(ii) for a hot disk, spiral structure induced errors for all Oort constants vanish at, and just inward of the corotation radius; (iii) as one approaches the 4:1 indblad resonances |C| increases and so does its variation with galactic azimuth; (iv) for all simulations |C|, on average, is larger for lower stellar velocity dispersions, contrary to recent measurements.

Two dimensional particle disks in proximity to a planet are numerically integrated to determine when a planet in a circular orbit can truncate a particle disk. Collisions are treated by giving each particle a series of velocity perturbations during the integration. We estimate the mass of a planet required to truncate a particle disk as a function of collision rate, related to the disk optical depth, and velocity perturbation size, related to the disk velocity dispersion. We find that for particle disks in the regime estimated for debris disks, a Neptune mass planet is sufficiently massive to truncate the disk. If both the velocity dispersion and the disk optical depth are low (dispersion less than approximately 0.02 in units of circular motion, and optical depth less than 10^-4) then an Earth mass planet suffices. We find that the disk is smooth and axisymmetric unless the velocity perturbation is small and the planet mass is of order or greater than a Neptune mass in which case azimuthal structure is seen near prominent mean motion resonances.

We investigate the evolution of fossil cavities produced by extinct YSO jets and wide angle outflows. Fossil cavities are ellipsoidal or cylindrical shells of swept-up of ambient (molecular cloud) material moving at low velocities. The cavities form when the momentum in a YSO jet or wide angle outflow decays in time allowing the bowshock or swept-up shell to decelerate to velocities near the turbulent speed in the cloud. It has been suggested in previous studies that cavities provide efficient coupling between the jets/outflows and the cloud and, as such, are the agents by which cloud turbulence is can be re-energized. In this paper we carry forward a series of numerical simulations of jets and outflows whose momentum flux decrease in time. We compare simulations with decaying momentum fluxes to those with constant flux. We show that decaying flux models exhibit deceleration of the outflow head and backfilling via expansion off of the cavity walls. They also have lower density contrast, are longer lived and wider than their continuously driven counterparts. The simulations recover the basic properties of observed fossil cavities. In addition, we provide synthetic observations in terms of P-V diagrams which demonstrate that fossil cavities form both jets and wide angle outflows are characterized by linear "Hubble-law" expansions patterns superimposed on "spur" patterns indicative of the head of a bow shock.

We consider particles with low free or proper eccentricity that are orbiting near planets on eccentric orbits. Via collisionless particle integration we numerically find the location of the boundary of the chaotic zone in the planet's corotation region. We find that the distance in semi-major axis between the planet and boundary depends on the planet mass to the 2/7 power and is independent of the planet eccentricity, at least for planet eccentricities below 0.3. Our integrations reveal a similarity between the dynamics of particles at zero eccentricity near a planet in a circular orbit and with zero free eccentricity particles near an eccentric planet. The 2/7 law has been previously explained by estimating the semi-major at which the first order mean motion resonances are large enough to overlap. Orbital dynamics near an eccentric planet could differ due to first order corotation resonances that have strength proportional to the planet's eccentricity. However, we find the corotation resonance width at low free eccentricity is small. Also the first order resonance width at zero free eccentricity is the same as that for a zero eccentricity particle near a planet in a circular orbit. This accounts for insensitivity of the chaotic zone width to planet eccentricity. Particles at zero free eccentricity near an eccentric planet have similar dynamics to those at zero eccentricity near a planet in a circular orbit.

We investigate the effect of a planet on an eccentric orbit on a two dimensional low mass gaseous disk. At a planet eccentricity above the planet's Hill radius divided by its semi-major axis, we find that the disk morphology differs from that exhibited by a disk containing a planet in a circular orbit. An eccentric gap is created with eccentricity that can exceed the planet's eccentricity and precesses with respect to the planet's orbit. We find that a more massive planet is required to open a gap when the planet is on an eccentric orbit. We attribute this behavior to spiral density waves excited at corotation resonances that increase the disk's eccentricity and exert a torque opposite in sign to that exerted by the Lindblad resonances. The reduced torque makes it more difficult for waves driven by the planet to overcome viscous inflow in the disk. Spectral energy distributions of transitional disks may reveal the presence of an eccentric planet if they are matched by an inner edge comprising a range of temperatures set by the semi-major axis and eccentricity of the hole.

We propose that the eccentricity and sharpness of the edge of Fomalhaut's disk are due to a planet just interior to the ring edge. The collision timescale consistent with the disk opacity is long enough that spiral density waves cannot be driven near the planet. The ring edge is likely to be located at the boundary of a chaotic zone in the corotation region of the planet. We find that this zone can open a gap in a particle disk as long as the collision timescale exceeds the removal or ejection timescale in the zone. We use the slope measured from the ring edge surface brightness profile to place an upper limit on the planet mass. The removal timescale in the chaotic zone is used to estimate a lower limit. The ring edge has eccentricity caused by secular perturbations from the planet. These arguments imply that the planet has a mass between that of Neptune and that of Saturn, a semi-major axis of approximately 119 AU and longitude of periastron and eccentricity, 0.1, the same as that of the ring edge.

We present the first data release of the Radial Velocity Experiment (RAVE), an ambitious spectroscopic survey to measure radial velocities and stellar atmosphere parameters (temperature, metallicity, surface gravity) of up to one million stars using the 6dF multi-object spectrograph on the 1.2-m UK Schmidt Telescope of the Anglo-Australian Observatory (AAO). The RAVE program started in 2003, obtaining medium resolution spectra (median R=7,500) in the Ca-triplet region ($\lambda\lambda$ 8,410--8,795 \AA) for southern hemisphere stars drawn from the Tycho-2 and SuperCOSMOS catalogs, in the magnitude range 9<I<12. The first data release is described in this paper and contains radial velocities for 24,748 individual stars (25,274 measurements when including re-observations). Those data were obtained on 67 nights between 11 April 2003 to 03 April 2004. The total sky coverage within this data release is $\sim$4,760 square degrees. The average signal to noise ratio of the observed spectra is 29.5, and 80% of the radial velocities have uncertainties better than 3.4 km/s. Combining internal errors and zero-point errors, the mode is found to be 2 km/s. Repeat observations are used to assess the stability of our radial velocity solution, resulting in a variance of 2.8 km/s. We demonstrate that the radial velocities derived for the first data set do not show any systematic trend with color or signal to noise. The RAVE radial velocities are complemented in the data release with proper motions from Starnet 2.0, Tycho-2 and SuperCOSMOS, in addition to photometric data from the major optical and infrared catalogs (Tycho-2, USNO-B, DENIS and 2MASS). The data release can be accessed via the RAVE webpage: http://www.rave-survey.org.

We present newly acquired images of the near-infrared counterpart of 3CR radio sources. All the sources were selected to have a redshift of less than 0.3 to allow us to obtain the highest spatial resolution. The observations were carried out as a snapshot program using the Near-Infrared Camera and Multiobject Spectrograph (NICMOS) on-board the Hubble Space Telescope (HST). In this paper we describe 69 radio galaxies observed for the first time with NICMOS during HST cycle 13. All the objects presented here are elliptical galaxies. However, each of them has unique characteristics such as close companions, dust lanes, unresolved nuclei, arc-like features, globular clusters and jets clearly visible from the images or with basic galaxy subtraction.

We report the discovery of a new optical-IR synchrotron jet in the radio galaxy 3C133 from our HST/NICMOS snapshot survey. The jet and eastern hotspot are well resolved, and visible at both optical and IR wavelengths. The IR jet follows the morphology of the inner part of the radio jet, with three distinct knots identified with features in the radio. The radio-IR SED's of the knots are examined, along with those of two more distant hotspots at the eastern extreme of the radio feature. The detected emission appears to be synchrotron, with peaks in the NIR for all except one case, which exhibits a power-law spectrum throughout.

Among radio galaxies containing nuclear dust disks, the bipolar jet axis is generally observed to be perpendicular to the disk major axis. The FR I radio source 3C 449 is an outlier to this statistical majority, as it possesses a nearly parallel jet/disk orientation on the sky. We examine the 600 pc dusty disk in this galaxy with images from the Hubble Space Telescope. We find that a colormap of the disk exhibits a twist in its isocolor contours (isochromes). We model the colormap by integrating galactic starlight through an absorptive disk, and find that the anomalous twist in the isochromes can be reproduced in the model with a vertically thin, warped disk. The model predicts that the disk is nearly perpendicular to the jet axis within 100 pc of the nucleus. We discuss physical mechanisms capable of causing such a warp. We show that precessional models or a torque on the disk arising from a possible binary black hole in the AGN causes precession on a timescale that is too long to account for the predicted disk morphology. However, we estimate that the pressure in the X-ray emitting interstellar medium is large enough to perturb the disk, and argue that jet-driven anisotropy in the excited ISM may be the cause of the warp. In this way, the warped disk in 3C 449 may be a new manifestation of feedback from an active galactic nucleus.

Spitzer mid-infrared images of the dusty warped disk in the galaxy Centaurus A show a parallelogram-shaped structure. We successfully model the observed mid-infrared morphology by integrating the light from an emitting, thin, and warped disk, similar to that inferred from previous kinematic studies. The models with the best match to the morphology lack dust emission within the inner 0.1 to 0.8 kpc, suggesting that energetic processes near the nucleus have disturbed the inner molecular disk, creating a gap in the molecular gas distribution.

Spitzer Space Telescope mid-infrared images of the radio galaxy Centaurus A reveal a shell-like, bipolar, structure 500 pc to the north and south of the nucleus. This shell is seen in 5.8, 8.0 and 24 micron broad-band images. Such a remarkable shell has not been previously detected in a radio galaxy and is the first extragalactic nuclear shell detected at mid-infrared wavelengths. We estimate that the shell is a few million years old and has a mass of order million solar masses. A conservative estimate for the mechanical energy in the wind driven bubble is 10^53 erg. The shell could have created by a small few thousand solar mass nuclear burst of star formation. Alternatively, the bolometric luminosity of the active nucleus is sufficiently large that it could power the shell. Constraints on the shell's velocity are lacking. However, if the shell is moving at 1000 km/s then the required mechanical energy would be 100 times larger.

A reverberation-mapping program on NGC 4395, the least-luminous known Seyfert 1 galaxy, undertaken with the Space Telescope Imaging Spectrograph on the Hubble Space Telescope, yields a measurement of the mass of the central black hole of 360,000 solar masses. The observations consist of two visits of 5 orbits each, in 2004 April and July. During each of these visits, the UV continuum varied by at least 10% (rms) and only C IV 1549 showed corresponding variations large enough to reliably determine the emission-line lag, which was measured to be of order one hour for both visits. The size of the C IV-emitting region is about a factor of three smaller than expected if the slope of the broad-line region radius-luminosity relationship is identical to that for the H-beta emission line. NGC 4395 is underluminous even for its small black hole mass; the Eddington ratio of 0.0012 is lower than that of any other active galactic nucleus for which a black hole mass measurement has been made by emission-line reverberation.

A migrating planet can capture planetesimals into mean motion resonances. However, resonant trapping can be prevented when the drift or migration rate is sufficiently high. Using a simple Hamiltonian system for first and second order resonances, we explore how the capture probability depends on the order of the resonance, drift rate and initial particle eccentricity. We present scaling factors as a function of the planet mass and resonance strength to estimate the planetary migration rate above which the capture probability drops to less than 1/2. Applying our framework to multiple extra solar planetary systems that have two planets locked in resonance, we estimate lower limits for the outer planet's migration rate allowing resonance capture of the inner planet. Mean motion resonances are comprised of multiple resonant subterms. We find that the corotation subterm can reduce the probability of capture when the planet eccentricity is above a critical value. We present factors that can be used to estimate this critical planet eccentricity. Applying our framework to the migration of Neptune, we find that Neptune's eccentricity is near the critical value that would make its 2:1 resonance fail to capture twotinos. The capture probability is affected by the separation between resonant subterms and so is also a function of the precession rates of the longitudes of periapse of both planet and particle near resonance.

We consider a differentially rotating, 2D stellar disk perturbed by two steady state spiral density waves moving at different patterns speeds. Our investigation is based on direct numerical integration of initially circular test-particle orbits. We examine a range of spiral strengths and spiral speeds and show that stars in this time dependent gravitational field can be heated (their random motions increased).This is particularly noticeable in the simultaneous propagation of a 2-armed spiral density wave near the corotation resonance (CR), and a weak 4-armed one near the inner and outer 4:1 Lindblad resonances. In simulations with 2 spiral waves moving at different pattern speeds we find: (1) the variance of the radial velocity, sigma_R^2, exceeds the sum of the variances measured from simulations with each individual pattern; (2) sigma_R^2 can grow with time throughout the entire simulation; (3) sigma_R^2 is increased over a wider range of radii compared to that seen with one spiral pattern; (4) particles diffuse radially in real space whereas they don't when only one spiral density wave is present. Near the CR with the stronger, 2-armed pattern, test particles are observed to migrate radially. These effects take place at or near resonances of both spirals so we interpret them as the result of stochastic motions. This provides a possible new mechanism for increasing the stellar velocity dispersion in galactic disks. If multiple spiral patterns are present in the Galaxy we predict that there should be large variations in the stellar velocity dispersion as a function of radius.

We propose that the two armed spiral features seen in visible Hubble Space Telescope images of scattered light in HD100546's circumstellar disk are caused by the illumination of a warped outer disk. A tilt of 6-15 degrees from the symmetry plane can cause the observed surface brightness variations providing the disk is very twisted (highly warped) at radii greater than 200 AU where the spiral features are seen. Dust lanes are due in part to shadowing in the equatorial plane from the inner disk within a radius of 100 AU. HD100546's outer disk, if viewed edge-on, would appear similar to that of Beta Pictorus. A disk initially misaligned with a planetary system, becomes warped due to precession induced by planetesimal bodies and planets. However, the twistedness of HD100546's disk cannot be explained by precession during the lifetime of the system induced by a single Jovian mass planet within the clearing at ~13 AU. One possible explanation for the corrugated disk is that precession was induced by massive of bodies embedded in the disk at larger radius. This would require approximately a Jupiter mass of bodies well outside the central clearing at 13 AU and within the location of the spiral features or at radii approximately between 50-200 AU.

Recent spectral observations by the Spitzer Space Telescope (SST) reveal that some discs around young ($\sim {\rm few} \times 10^6$ yr old) stars have remarkably sharp transitions to a low density inner region in which much of the material has been cleared away. It has been recognized that the most plausible mechanism for the sharp transition at a specific radius is the gravitational influence of a massive planet. This raises the question of whether the planet can also account for the hole extending all the way to the star. Using high resolution numerical simulations, we show that Jupiter-mass planets drive spiral waves which create holes on time scales $\sim 10$ times shorter than viscous or planet migration times. We find that the theory of spiral-wave driven accretion in viscous flows by Takeuchi et al. (1996) can be used to provide a consistent interpretation of the simulations. In addition, although the hole surface densities are low, they are finite, allowing mass accretion toward the star. Our results therefore imply that massive planets can form extended, sharply bounded spectral holes which can still accommodate substantial mass accretion rates. The results also imply that holes are more likely than gaps for Jupiter mass planets around solar mass stars.

We study the effects of an annular gap induced by an embedded proto-planet on disk scattered light images and the infrared spectral energy distribution. We find that the outer edge of a gap is brighter in the scattered light images than a similar location in a gap-free disk. The stellar radiation that would have been scattered by material within in the gap is instead scattered by the disk wall at the outer edge of the gap, producing a bright ring surrounding the dark gap in the images. Given sufficient resolution, such gaps can be detected by the presence of this bright ring in scattered light images. A gap in a disk also changes the shape of the SED. Radiation that would have been absorbed by material in the gap is instead reprocessed by the outer gap wall. This leads to a decrease in the SED at wavelengths corresponding to the temperature at the radius of the missing gap material and a corresponding flux increase at longer wavelengths corresponding to the temperature of the outer wall. We note that, unlike an inner hole in the disk, the presence of an annular gap does not change the bolometric IR flux; it simply redistributes the radiation, previously produced by material within the gap, to longer wavelengths. This implies that the changes in the SED generally will be smaller for gaps than holes. Although it will be difficult on the basis of the SED alone to distinguish between the presence of a gap and other physical effects, the level of changes can be sufficiently large to be measurable with current instruments (e.g., Spitzer).

With 2D hydrodynamical simulations of disks perturbed externally by stars, brown dwarfs or planets we investigate possible scenarios that can account for the spiral structure in circumstellar disks. We consider two scenarios, spiral structure driven by an external bound planet or low mass star and that excited by a previous stellar close encounter or flyby. We find that both scenarios produce morphology similar to that observed in the outer disks of HD 141569A and HD 100546; moderately open 2-armed outer spiral structure. The outer two-armed spiral structure observed in the disk of HD 141569A is qualitatively reproduced with tidal perturbations from its companion binary HD 141569B,C on a prograde orbit near periapse. Our simulation accounts for the outer spiral arms, but is less successful than the secular model of Augereau and Papaloizou at matching the lopsidedness or asymmetry of the disk edge at 300AU. The disk has been previously truncated by the tidal force from the binary. A bound object (stellar or planetary) is unlikely to explain the spiral structure in HD 100546. A co-eval planet or brown dwarf in the disk of sufficient mass to account for the amplitude of the spiral structure would be detectable in NICMOS and STIS images, however existing images reveal no such object. A previous encounter could explain the observed structure, provided that the encounter occurred less than a few thousand year ago. The object responsible for causing the spiral structure must then be within a few arcminutes of the star. However, the USNO-B proper motion survey reveals no candidate object. Moreover, the probability that a field star encountered HD 100546 in the past few thousand years is very low.

We present submillimeter continuum observations at 870 microns of the cores of low redshift 3CRR radio galaxies, observed at the Heinrich Hertz Submillimeter Telescope. The cores are nearly flat spectrum between the radio and submillimeter which implies that the submillimeter continuum is likely to be synchrotron emission and not thermal emission from dust. The emitted power from nuclei detected at optical wavelengths and in the X-rays is similar in the submillimeter, optical and X-rays. The submillimeter to optical and X-ray power ratios suggest that most of these sources resemble misdirected BL Lac type objects with synchrotron emission peaking at low energies. However we find three exceptions, the FR I galaxy 3C264 and the FR II galaxies 3C390.3 and 3C338 with high X-ray to submillimeter luminosity ratios. These three objects are candidate high or intermediate energy peaked BL Lac type objects. With additional infrared observations and from archival data, we compile spectral energy distributions (SEDs) for a subset of these objects. The steep dips observed near the optical wavelengths in many of these objects suggest that extinction inhibits the detection and reduces the flux of optical continuum core counterparts. High resolution near or mid-infrared imaging may provide better measurements of the underlying synchrotron emission peak.

Clumps in the solar neighborhood's stellar velocity distribution could be caused by spiral density waves. In the solar neighborhood, stellar velocities corresponding to orbits that are nearly closed in the frame rotating with a spiral pattern represent likely regions for stellar concentrations. Via particle integration, we show that orbits can intersect the solar neighborhood when they are excited by Lindblad resonances with a spiral pattern. We find that a two-armed spiral density wave with pattern speed placing the Sun near the 4:1 Inner Lindblad Resonance (ILR) can cause two families of nearly closed orbits in the solar neighborhood. One family corresponds to square shaped orbits aligned so their peaks lie on top of, and support, the two dominant stellar arms. The second family correspond to orbits 45 degrees out of phase with the other family. Such a spiral density pattern could account for two major clumps in the solar neighborhood's velocity distribution. The Pleiades/Hyades moving group corresponds to the first family of orbits and the Coma Berenices moving group corresponds to the second family. This model requires a spiral pattern speed of approximately 0.66 +- 0.03 times the angular rotation rate of the Sun or 18.1 +- 0.8 km/s/kpc.

We investigate gap formation in gaseous protostellar disks by a planet in a circular orbit in the limit of low disk viscosity. This regime may be appropriate to an aging disk after the epoch of planet formation. We find that the distance of planet to the gap outer boundary can be between the location of the $m=2$ and $m=1$ outer Lindblad resonances. This distance is weakly dependent upon both the planet's mass and disk viscosity. We find that the evolution of the disk edge takes place on two timescales. The first timescale is set by the spiral density waves driven by the planet. The second timescale depends on the viscosity of the disk. The disk approaches a state where the outward angular momentum flux caused by the disk viscosity is balanced by the dissipation of spiral density waves which are driven at the Lindblad resonances. This occurs inefficiently however because of the extremely low gas density near the planet. We find that the distance between the planet and the peak density at the disk outer edge is only weakly dependent on the viscosity and planet mass, however the ratio of the gas density near the planet to that in the disk (or the slope of density along the disk edge) is strongly dependent upon both quantities. We find that the disk density profile along the edge scales approximately with disk viscosity divided by the square of the planet mass. We account for this behavior with a simple scenario in which the dissipation of angular momentum from the spiral density waves is balanced against diffusion in the steep edge of the disk.

Outflows from young stellar objects have been identified as a possible source of turbulence in molecular clouds. To investigate the relationship between outflows, cloud dynamics and turbulence, we compare the kinematics of the molecular gas associated with NGC 1333, traced in 13CO(1-0), with the distribution of young stellar objects (YSOs) within. We find a velocity dispersion of ~ 1-1.6 km/s in 13CO that does not significantly vary across the cloud, and is uncorrelated with the number of nearby young stellar outflows identified from optical and submillimeter observations. However, from velocity channel maps we identify about 20 cavities or depressions in the 13CO intensity of scales > 0.1-0.2 pc and velocity widths 1-3 km/s. The cavities exhibit limb brightened rims in both individual velocity channel maps and position velocity diagrams, suggesting that they are slowly expanding. We interpret these cavities to be remnants of past YSO outflow activity: If these cavities are presently empty, they would fill in on time scales of a million years. This can exceed the lifetime of a YSO outflow phase, or the transit time of the central star through the cavity, explaining the the absence of any clear correlation between the cavities and YSO outflows. We find that the momentum and energy deposition associated with the expansion of the cavities is sufficient to power the turbulence in the cloud. In this way we conclude that the cavities are an important intermediary step between the conversion of YSO outflow energy and momentum into cloud turbulent motions.

We investigate the relation between Forbush cosmic ray decrease recovery time and coronal mass ejection transit time between the Sun and Earth. We identify 17 Forbush decreases from ground based neutron count rates between 1978 and 2003 that occur during the same phase in the solar cycle and can be associated with single coronal mass ejections (CMEs) in the SOHO LASCO CME Catalog or previously published reports, and with specific interplanetary coronal mass ejections (ICMEs) crossing the vicinity of Earth. We find an anti-correlation between Forbush recovery times and CME transit time that contradicts the predictions of simple cosmic ray diffusive barrier models. The anti-correlation suggests that the decay rate of ICMEs is anti-correlated with their travel speed. Forbush recovery times range from seven times the transit time for the fastest disturbance to a fifth the Sun-Earth transit time for the slowest. To account for the large range of measured recovery times we infer that the slowest disturbances must decay rapidly with radius whereas the fastest ones must remain strong. The longest recovery times suggest that the fastest disturbances in our sample decayed less rapidly with radius than the ambient solar wind magnetic field strength.

Spitzer observations of the young star CoKuTau/4 reveal a disk with a 10 AU hole that is most likely caused by a newly formed planet. Assuming that the planet opened a gap in the viscous disk, we estimate that the planet mass is greater than 0.1 Jupiter masses. This estimate depends on a lower limit to the disk viscosity derived from the time scale needed to accrete the inner disk, creating the now detectable hole. The planet migration time scale must at least modestly exceed the time for the spectrally inferred hole to clear. The proximity of the planet to the disk edge implied by our limits suggests that the latter is perturbed by the nearby planet and may exhibit a spiral pattern rotating with the planet. This pattern might be resolved with current ground based mid-infrared cameras and optical cameras on the Hubble Space Telescope. The required sub-Myr planet formation may challenge core accretion formation models. However, we find that only if the planet mass is larger than about 10 Jupiter masses, allowing for a high enough surface density without inducing migration, would formation by direct gravitational instability be possible.

We numerically investigate the possibility that a close stellar encounter could account for the high inclinations of the Kuiper belt, as originally proposed by Ida, Larwood and Burkert, however we consider encounters with pericenters within those explored by previous works. A star of mass 0.2 Msol entering the solar system on a nearly parabolic, low inclination, retrograde orbit with perihelion of 50 AU can scatter ~30% of the Kuiper belt into moderate inclination (i>10 degrees) and eccentricity orbits, while leaving the rest of the belt at low eccentricities and inclinations. This scenario provides a possible explanation for the dichotomy of the Kuiper belt, accounting for the scattered and classical populations, including objects with high eccentricities and inclinations. High eccentricity objects that were not placed into resonance with Neptune are more likely to be removed from the system by subsequent long timescale evolution. Consequently this scenario suggests that some Plutinos could have originally been objects in circular orbits which were scattered into resonant high eccentricity and high inclination orbits by the stellar encounter. Objects could have been placed into stable resonant regions with Neptune by the encounter. However long lived regions of dynamical stability should differ from those selected by resonance capture. This scenario predicts the presence of objects in resonances with Neptune that are unlikely to have been captured via Neptune's migration.

We compare the soft diffuse X-ray emission from Chandra images of 12 nearby intermediate inclination spiral galaxies to the morphology seen in Halpha, molecular gas, and mid-infrared emission. We find that diffuse X-ray emission is often located along spiral arms in the outer parts of spiral galaxies but tends to be distributed in a rounder morphology in the center. The X-ray morphology in the spiral arms matches that seen in the mid-infrared or Halpha and so implies that the X-ray emission is associated with recent active star formation. We see no strong evidence for X-ray emission trailing the location of high mass star formation in spiral arms. However, population synthesis models predict a high mechanical energy output rate from supernovae for a time period that is about 10 times longer than the lifetime of massive ionizing stars, conflicting with the narrow appearance of the arms in X-rays. The fraction of supernova energy that goes into heating the ISM must depend on environment and is probably higher near sites of active star formation. The X-ray estimated emission measures suggest that the volume filling factors and scale heights are high in the galaxy centers but low in the outer parts of these galaxies. The differences between the X-ray properties and morphology in the centers and outer parts of these galaxies suggest that galactic fountains operate in outer galaxy disks but that winds are primarily driven from galaxy centers.

In the inner 3kpc of M51 we find that logarithmic spirals provide good fits to the peak intensities in molecular gas observed by BIMA in the CO (J=1-0) emission line along the spiral arms. However, we measure significant asymmetries between the location and density of the arms observed on one side of the galaxy compared to those on the opposite side. Between a radius of 1 and 2.2 kpc the gas distribution traced in CO is lopsided with densities twice as large in one arm compared to the opposite one. In the same region, the spiral arms are offset by 20 degrees from the position of the arm on the opposite side of the galaxy after a rotation of 180 degrees. We use the ratio of CO emission to that in the Pa alpha hydrogen recombination line to estimate the efficiency of star formation along each arm. Except for a narrow region at about 2kpc where star formation is surprisingly inefficient, the gas depletion time is approximately 200 million years despite large variations in the molecular gas density in the arms. We account for the deviations from two armed bi-symmetric structure with a simple model that includes an additional three-armed spiral density wave. This model accounts for the angular offset between the arm on one side compared to its opposite, the lopsided distribution in molecular gas, and inter-arm star formation to the north east of the galaxy nucleus. Since the star formation efficiency is unaffected by the variations in gas density, and the variations in gas density can be accounted for by the presence of an additional three-armed density wave, we suggest that the star formation rate is variable, and is highest where and when the maxima of the two spiral density waves coincide or constructively add.

We present nuclear spectral energy distributions (SEDs) in the range 0.4-16micron for an expanded CfA sample of Seyfert galaxies. The spectral indices from 1-16micron range from alpha_IR 0.9 to 3.8. The shapes of the spectra are correlated with Seyfert type in the sense that steeper nuclear SEDs (nu*f_nu increasing with increasing wavelength) tend to be found in Seyfert 2s and flatter SEDs (nu *f_nu constant) in Seyfert 1-1.5s. The galaxies optically classified as Seyferts 1.8s and 1.9s display values of alpha_IR as in type 1 objects, or values intermediate between those of Seyfert 1s and Seyfert 2s. The intermediate SEDs of many Seyfert 1.8-1.9s may be consistent with the presence of a pure Seyfert 1 viewed through a moderate amount (A_V <5mag) of foreground galaxy extinction. We find, however, that between 10 and 20% of galaxies with broad optical line components have steep infrared SEDs. Torus models usually adopt high equatorial opacities to reproduce the infrared properties of Seyfert 1s and 2s, resulting in a dichotomy of infrared SEDs (flat for type 1s, and steep for type 2s). Such a dichotomy, however, is not observed in our sample. The wide range of spectral indices observed in the type 2 objects, the lack of extremely steep SEDs, and the large numbers of objects with intermediate spectral indices cannot be reconciled with predictions from existing optically thick torus models. We discuss possible modifications to improve torus models, including low optical depth tori, clumpy dusty tori, and high-optical-depth tori with an extended optically thin component.

The star S0-2, which is orbiting Sgr A* with a 15-yr period, almost certainly did not form in situ. We propose that it was injected into this close orbit by the tidal disruption of a massive-star binary, whose primary was more massive than S0-2 and at least 60 Msun. From numerical integrations we find that 1-2% of incoming binaries with closest approach equal to 130 AU leave the secondary in an orbit with eccentricity within 0.01 of that of SO-2. If additional stars are found orbiting Sgr A* with relatively short periods, they could be used to probe the formation of massive stars in the Galactocentric region, even though the massive stars themselves have long since perished.

We consider the possibility that an infalling bulge core or stellar cluster could form an eccentric disk, following tidal disruption by a massive black hole in the center of a galaxy. As a function of central black hole mass, we constrain the core radii and central densities of cluster progenitors capable of becoming nearly Keplerian disks which can support lopsided slow modes. We find that progenitor stellar clusters with core radii less than a pc and densities above a few times 10^5 Msol/pc^3 are likely eccentric disk progenitors near a massive black hole of mass 10^7 to 10^8 Msol. Lower density and larger progenitor cores are capable of causing eccentric stellar disks near more massive black holes. Our constraints on the progenitor cores are consistent with existing N-body simulations, which in one case has produced an eccentric disk. For M31 and NGC 4486B, the estimated progenitor cluster cores are dense and compact compared to Galactic globular clusters, however the cores of nearby galaxies such as M33, M32 and M31 itself are in the right regime. If galaxy mergers can create eccentric disks, then they would be a natural consequence of hierarchical galaxy formation.

The morphology of the epsilon Eridani dust ring is reproduced by a numerical simulation of dust particles captured into the 5:3 and 3:2 exterior mean-motion resonances with a 0.3 eccentricity 10^-4 solar mass planet at periastron at a semi-major axis of 40 AU. The morphology will differ when the planet is at aphelion, in about 140 years. Moderate eccentricity planets in outer extra-solar systems will cause observable variations in the morphology of associated dusty rings.

We present a simple resonant Hamiltonian model for the vertical response of a stellar disk to the growth of a bar perturbation. As a bar perturbation grows stars become trapped in vertical Inner Lindblad resonances and are lifted into higher amplitude orbits. The vertical structure of a boxy and peanut shaped bulge as a function of radius and azimuthal angle in the galaxy plane can be predicted from the strength and speed of the bar perturbation and the derivatives of the gravitational potential. This model predicts that stars on the outer side of the resonance are lifted higher than stars on the inner side, offering an explanation for the sharp outer edge of the boxy/peanut.

We announce the initial release of data from the Ohio State University Bright Spiral Galaxy Survey, a BVRJHK imaging survey of a well-defined sample of 205 bright, nearby spiral galaxies. We present H-band morphological classification on the Hubble sequence for the OSU Survey sample. We compare the H-band classification to B-band classification from our own images and from standard galaxy catalogs. Our B-band classifications match well with those of the standard catalogs. On average, galaxies with optical classifications from Sa through Scd appear about one T-type earlier in the H-band than in the B-band, but with large scatter. This result does not support recent claims made in the literature that the optical and near-IR morphologies of spiral galaxies are uncorrelated. We present detailed descriptions of the H-band morphologies of our entire sample, as well as B- and H-band images for a set of 17 galaxies chosen as type examples, and BRH color-composite images of six galaxies chosen to demonstrate the range in morphological variation as a function of wavelength. Data from the survey are accessible at http://www.astronomy.ohio-state.edu/~survey/

We consider the nature of orbits near the solar neighborhood which are perturbed by local spiral arms and the Milky Way bar. We present a simplified Hamiltonian model which includes resonant terms from both types of perturbations and is similar to the forced pendulum. Via numerical integration of this model we construct Poincare maps to illustrate the nature and stability of the phase space. We find that resonance overlap is most likely to cause widespread chaos when the pattern of the spiral structure puts the solar neighborhood near the 2:1 inner Lindblad resonance (ILR) in the case of a 2-armed pattern, or near the 4:1 ILR in the case of a 4-armed pattern. When this happens both the quasiperiodic orbits which support the spiral structure and those that oscillate with the bar are disrupted near the bar's 2:1 outer Lindblad resonance (OLR). Consequently the pattern speed of spiral structure which passes through the OLR must be faster than 0.45 times the solar neighborhood angular rotation rate if it is 2-armed or faster than 0.75 times this value if it is 4-armed. Alternatively the OLR may form a boundary between spiral modes at different pattern speeds. In all cases we find that spiral structure is disrupted by the OLR over a narrow range of radius and the extent of the orbits aligned perpendicular to the bar at the OLR is limited by the spiral perturbations.

Using star counts in both color and magnitude from the Two Micron All Sky Survey (2MASS) Second Incremental Release Point Source Catalog we search for evidence of non-uniform extinction and stellar population density changes in the Galactic Plane. Extinction causes the entire main sequence to shift toward redder colors on a color magnitude diagram. A local increase in the stellar density causes an increase in the star counts along a line parallel to the main sequence. We find streaks in star count color magnitude contour plots along the angle of the main sequence which are likely to be caused by distant gas clouds and stellar density variations. The distance of a gas cloud or stellar density change can be estimated from the location of the shift in the star count contours. We identify features in these diagrams which are coherent across at least 10 degrees in Galactic longitude. A series of features is seen at the plausible distance of the expected Perseus spiral arm at a distance of 2 to 4 kpc from the sun. However other features as prominent are found at both at larger and smaller distances. These structures are over 300 pc in size and so likely to be associated with large scale coherent structures in the gas distribution such as weak spiral arms. The presence of multiple and weak spiral arms, and lack of strong ones suggests that the outer Milky Way is flocculent in its morphology.

Hipparcos parallaxes and proper motions have made it possible to construct Hertzsprung-Russell (HR) diagrams of nearby clusters with unprecedented accuracy. The standard deviation of high fidelity non-binary non-variable stars about a model stellar evolution isochrone in the Hyades cluster is about 0.04 magnitudes. We use this deviation to estimate an upper limit on the scatter in metallicities in stars in this cluster. From the gradient of the isochrones evolution in the HR diagram we estimate an upper limit for the scatter of metallicities Delta [Fe/H] <~ 0.03 dex, a smaller limit than has been measured previously spectroscopically. This suggests that stars in open clusters are formed from gas that is nearly homogeneous in its metallicity. We consider the hypothesis that processes associated with planet formation can pollute the convection zone of stars. If the position on the HR diagram is insensitive to the metallicity of the convection zone and atmosphere, then stars which have very polluted convection zones can be identified from a comparison between their metallicity and position on the HR diagram. Alternatively if the pollution of the star by metals results in a large change in the position of the star on the HR diagram in a direction perpendicular to the isochrone, then the low scatter of stars in the Hyades can be used to place constraints on quantity of high-Z material that could have polluted the stars.

A wind passing over a surface may cause an instability in the surface such as the flapping seen when wind blows across a flag or waves when wind blows across water. We show that when a radially outflowing wind blows across a dense thin rotating disk, an initially flat disk is unstable to warping. When the wind is subsonic, the growth rate is dependent on the lift generated by the wind and the phase lag between the pressure perturbation and the vertical displacement in the disk caused by drag. When the wind is supersonic, the grow rate is primarily dependent on the form drag caused by the surface. While the radiative warping instability proposed by Pringle is promising for generating warps near luminous accreting objects, we expect the wind driven instability introduced here would dominate in objects which generate energetic outflows.

Because a chaotic zone can reduce the long timescale capture probabilities and cause catastrophic events such as close encounters with a planet or star during temporary capture, the dynamics of migrating planets is likely to be strongly dependent on the widths of the resonance chaotic zones. Previous theoretical work on the resonant capture of particles into mean-motion resonances by orbital migration has been restricted to the study of integrable models. By exploring toy drifting Hamiltonian models we illustrate how the structure in phase space of a resonance can be used to generalize this integrable capture theory to include the richer phenomenology of chaotic resonances. We show that particles are temporarily captured into the chaotic zone of a resonance with fixed shape for a time that is given by the width of the chaotic zone divided by the resonance drift rate. Particles can be permanently captured into a drifting chaotic resonance only if they are captured into a growing integrable region. Therefore resonances containing wide chaotic zones have lower permanent capture probabilities than those lacking chaotic zones. We expect large deviations from the predictions of integrable capture theories when the chaotic zone is large and the migration rate is sufficiently long that many Lyapunov times pass while particles are temporarily captured in the resonance. The 2:1 mean-motion resonance with Neptune in the Kuiper Belt contains a chaotic zone, even when integrated on short timescales such as a million years. Capture probabilities are lower than estimated previously from drifting integrable models. This may offer an explanation for low eccentricity Kuiper Belt objects between 45-47 AU which should have been previously captured in the 2:1 resonance by Neptune's migration.

We measure reddenings to HII regions in NGC 2903, NGC 1512, M51, NGC 4449 and NGC 6946 from Hubble Space Telescope Pa alpha and H alpha images. Extinctions range from A_V ~ 5 - 0 depending upon the galaxy. For the galaxies with HST images in both lines, NGC 2903, NGC 1512 and M51, the Pa alpha and H alpha emission are almost identical in morphology which implies that little emission from bright HII regions is hidden from view by regions of comparatively high extinction. The scatter in the measured extinctions is only +- 0.5 mag. We compare the reddenings we measure in five galaxies using the Pa alpha to H alpha ratios to those measured previously from the Balmer decrement in the LMC and as a function of radius in M101 and M51. We find that luminosity weighted mean extinctions of these ensembles of HI regions are correlated with gas surface density and metallicity. The correlation is consistent with the mean extinction depending on dust density where the dust to gas mass ratio scales with the metallicity. This trend is expected if HII regions tend to be located near the mid-plane of a gas disk and emerge from their parent molecular clouds soon after birth. In environments with gas densities below a few hundred Msol/pc^2 star formation rates estimated from integrated line fluxes and mean extinctions are likely to be fairly accurate.

JHKL'M (1-5 micron) imaging of a sample of Seyfert (Sy) 2 galaxies is presented. We confirm that the 1-2.2 micron nuclear continuum of some Sy 2s is dominated by stellar emission, whereas the continuum emission at longer wavelengths (3-5 micron) is almost entirely non-stellar in origin. The IR non-stellar spectral energy distributions (SED) (up to 15 micron) of Sy galaxies show a variety of shapes, and they are well reproduced with the tapered disk models of Efstathiou & Rowan-Robinson (1995). We have used two models, one including an optically thin cone component and a coneless model. Although our modelling of the SEDs does not allow us to favor either model to account for all the observed SEDs, we find that the viewing angle towards the central source is well constrained by both models. We have also investigated non-stellar color-color diagrams. The colors of the Sy galaxies with viewing angles theta < 30 degree are better reproduced with the cone model. These diagrams provide a good means to separate Sy 2s with moderate obscuration (A_V < 20 mag, from hard X-ray observations) from those with high obscuration. The ground-based 4.8 micron and ISO 9.6 micron luminosities are well correlated with the hard X-ray luminosities of Sy 1s and 2s. These continuum emissions appear as a good indicator of the AGN luminosity, at least in the cases of hard X-ray Compton-thin Sy galaxies. We finally stress the finding that some Compton thick galaxies show bright non-stellar emission at 5 micron. This suggests that the near-IR emission in Sys is produced in an extended component illuminated by the central source, that is more visible from all viewing angles. We discuss possible implications of mid-IR surveys for the search of counterparts of highly obscured hard X-ray sources. (Abridged)

The modest extinction and reasonably face-on viewing geometry make the luminous IR galaxy NGC1614 an ideal laboratory for study of a powerful starburst. HST/NICMOS observations show: 1.) deep CO stellar absorption, tracing a starburst nucleus about 45 pc in diameter; 2.) surrounded by a ~600 pc diameter ring of supergiant HII regions revealed in Pa$\alpha$ line emission; 3.) lying within a molecular ring indicated by its extinction shadow in H-K; 4.) all at the center of a disturbed spiral galaxy. The luminosities of the giant HII regions in the ring are extremely high, an order of magnitude brighter than 30 Doradus. The relation of deep stellar CO bands to surrounding ionized gas ring to molecular gas indicates that the luminous starburst started in the nucleus and is propagating outward into the surrounding molecular ring. This hypothesis is supported by evolutionary starburst modeling that shows that the properties of NGC1614 can be fitted with two short-lived bursts of star formation separated by 5 Myr. The total dynamical mass of the starburst region of 1.3x10^9 Msol is mostly accounted for by the old pre-star burst stellar population. Although our starburst models use a modified Salpeter IMF, the tight mass budget suggests that the IMF may contain relatively more 10-30 Msol stars and fewer low mass stars than the Salpeter function. The dynamical mass is nearly 4 times smaller than the mass of molecular gas estimated from the standard ratio of 12CO(1-0) to H2. A number of arguments place the mass of gas in the starburst region at ~25% of the dynamical mass, nominally about 1/15 and with an upper limit of 1/10 of the amount estimated from 12CO and the standard ratio. (Abridged)

We examine 112 Seyfert galaxies observed by the Hubble Space Telescope (HST) at 1.6 microns. We find that ~50% of the Seyfert 2.0 galaxies which are part of the Revised Shapeley-Ames (RSA) Catalog or the CfA redshift sample contain unresolved continuum sources at 1.6 microns. All but a couple of the Seyfert 1.0-1.9 galaxies display unresolved continuum sources. The unresolved sources have fluxes of order a mJy, near-infrared luminosities of order 10^41 erg/s and absolute magnitudes M_H ~-16. Comparison non-Seyfert galaxies from the RSA Catalog display significantly fewer (~20%), somewhat lower luminosity nuclear sources, which could be due to compact star clusters. We find that the luminosities of the unresolved Seyfert 1.0-1.9 sources at 1.6 microns are correlated with [OIII] 5007A and hard X-ray luminosities, implying that these sources are non-stellar. Assuming a spectral energy distribution similar to that of a Seyfert 2 galaxy, we estimate that a few percent of local spiral galaxies contain black holes emitting as Seyferts at a moderate fraction, 10^-1 to 10^-4, of their Eddington luminosities. With increasing Seyfert type the fraction of unresolved sources detected at 1.6 microns and the ratio of 1.6 microns to [OIII] fluxes tend to decrease. These trends are consistent with the unification model for Seyfert 1 and 2 galaxies.

We re-examine the age-velocity dispersion relation in the solar neighborhood using improved stellar age estimates that are based on Hipparcos parallaxes and recent stellar evolution calculations. The resulting relation shows that the Milky Way stellar disk was relatively quiescent, suffering little heating or dispersion increase between 3 and 9 Gyrs. However, at an age of 9 Gyrs there is an abrupt increase in the stellar velocity dispersions. To explain the abrupt increase (by almost a factor of 2) in dispersions we propose that the Milky Way suffered a minor merger 9 Gyrs ago that created the thick disk. The quiescent phase is consistent with inefficient heating caused by scattering from tightly wound transient spiral structure. It is possible that many galaxy disks have radial velocity dispersions at the point where heating from spiral structure becomes inefficient. Surface brightness profiles and dispersion measurements can be used to predict spiral arm winding numbers at this saturation point. A comparison between the observed and predicted winding numbers could test this hypothesis and also search for evidence of merger induced disk heating in nearby galaxies.

We present a study of Seyfert 1.5-2.0 galaxies observed at two epochs with the Hubble Space Telescope (HST) at 1.6 microns. We find that unresolved nuclear emission from 9 of 14 nuclei varies at the level of 10-40% on timescales of 0.7-14 months, depending upon the galaxy. A control sample of Seyfert galaxies lacking unresolved sources and galaxies lacking Seyfert nuclei show less than 3% instrumental variation in equivalent aperture measurements. This proves that the unresolved sources are non-stellar and associated with the central pc of active galactic nuclei. Unresolved sources in Seyfert 1.8 and 1.9 galaxies are not usually detected in HST optical surveys, however high angular resolution infrared observations will provide a way to measure time delays in these galaxies.

We have obtained new Hubble Space Telescope (HST) observations of M84, a nearby massive elliptical galaxy whose nucleus contains a 1.5 X 10^9 Msun dark compact object, which presumably is a supermassive black hole. Our Space Telescope Imaging Spectrograph (STIS) spectrum provides the first clear detection of emission lines in the blue (e.g., [O II] 3727, Hbeta, and [O III] 4959, 5007), which arise from a compact region 0.28" across centered on the nucleus. Our Near Infrared Camera and Multi-Object Spectrometer (NICMOS) images exhibit the best view through the prominent dust lanes evident at optical wavelengths and provide a more accurate correction for the internal extinction. The relative fluxes of the emission lines we have detected in the blue together with those detected in the wavelength range 6295 - 6867 AA by Bower et al. (1998, ApJ, 492, L111) indicate that the gas at the nucleus is photoionized by a nonstellar process, instead of hot stars. Stellar absorption features from cool stars at the nucleus are very weak. We update the spectral energy distribution of the nuclear point source and find that although it is roughly flat in most bands, the optical to UV continuum is very red, similar to the spectral energy distribution of BL Lac. Thus, the nuclear point source seen in high-resolution optical images (Bower et al. 1997, ApJ, 483, L33) is not a star cluster but is instead a nonstellar source. Assuming isotropic emission from this source, we estimate that the ratio of bolometric luminosity to Eddington luminosity is 5 X 10^(-7). However, this could be underestimated if this source is a misaligned BL Lac object, which is a possibility suggested by the spectral energy distribution and the evidence of optical variability we describe.

We have determined the fraction of barred galaxies in the H-band for a statistically well-defined sample of 186 spirals drawn from the Ohio State University Bright Spiral Galaxy survey. We find 56% of our sample to be strongly barred at H, while another 16% is weakly barred. Only 27% of our sample is unbarred in the near-infrared. The RC3 and the Carnegie Atlas of Galaxies both classify only about 30% of our sample as strongly barred. Thus strong bars are nearly twice as prevalent in the near-infrared as in the optical. The frequency of genuine optically hidden bars is significant, but lower than many claims in the literature: 40% of the galaxies in our sample that are classified as unbarred in the RC3 show evidence for a bar in the H-band, while for the Carnegie Atlas this fraction is 66%. Our data reveal no significant trend in bar fraction as a function of morphology in either the optical or H-band. Optical surveys of high redshift galaxies may be strongly biased against finding bars, as bars are increasingly difficult to detect at bluer rest wavelengths.

We report the discovery of gas inflow motions towards the active nucleus of the Circinus galaxy caused by the non-axisymmetric potential of a nuclear gas bar. Evidence for dust associated with the bar comes from the HST/NICMOS H-K color map, whereas the streaming motions along the gas bar are seen in the velocity field of the H2 S(1)(1-0) emission line. The gas bar is about 100 pc long with a visual extinction in excess of 10 mag. Indication for the gaseous nature of this bar comes from the lack of a stellar counterpart even in the K band where the extinction is greatly reduced. We also use the NICMOS emission line images (Pa-alpha, [SiVI], and [FeII]) to study the innermost region of the ionization cones and the nuclear star forming activity. We discuss the possible relationship of these components with the gaseous bar.

During orbital migration of a giant extrasolar planet via ejection of planetesimals (Murray et al.~1998), inner mean motion resonances can be strong enough to cause planetesimals to graze or impact the star. We integrate numerically the motions of particles which pass through the 3:1 or 4:1 mean motion resonances of a migrating Jupiter mass planet. We find that many particles can be trapped in the 3:1 or 4:1 resonances and pumped to high enough eccentricities that they impact the star. This implies that for a planet migrating a substantial fraction of its semi major axis, a significant fraction of its mass in planetesimals could impact the star. This process may be capable of enriching the metallicity of the star, and at a time when the star is no longer fully convective. Upon close approaches to the star the surfaces of these planetesimals will be sublimated. Orbital migration should cause continuing production of evaporating bodies, suggesting that this process should be detectable with searches for transient absorption lines in young stars. The remainder of the particles will not impact the star but can be subsequently ejected by the planet as it migrates further inwards. This allows the planet to migrate a substantial fraction of its initial semi-major axis via ejection of planetesimals.

We present NICMOS imaging of broad band and molecular hydrogen emission in Seyfert galaxies. In 6 of 10 Seyferts we detect resolved or extended emission in the 1-0 S(1) 2.121 or 1-0 S(3) 1.9570 micron molecular hydrogen lines. We did not detect emission in the most distant galaxy or in the 2 Seyfert 1 galaxies in our sample because of the luminosity of the nuclear point sources. In NGC 5643, NGC 2110 and MKN 1066, molecular hydrogen emission is detected in the extended narrow line region on scales of a few hundred pc from the nucleus. Emission is coincident with [OIII] and H alpha+[NII] line emission. This emission is also near dust lanes observed in the visible to near-infrared color maps suggesting that a multiphase medium exists near the ionization cones and that the morphology of the line emission is dependent on the density of the ambient media. The high 1-0 S(1) or S(3) H2 to H alpha flux ratio suggests that shock excitation of molecular hydrogen (rather than UV fluorescence) is the dominant excitation process in these extended features. In NGC 2992 and NGC 3227 the molecular hydrogen emission is from 800 and 100 pc diameter `disks' (respectively) which are not directly associated with [OIII] emission and are near high levels of extinction (AV > 10). In NGC 4945 the molecular hydrogen emission appears to be from the edge of a 100 pc superbubble. In these 3 galaxies the molecular gas could be excited by processes associated with local star formation. We confirm previous spectroscopic studies finding that no single mechanism is likely to be responsible for the molecular hydrogen excitation in Seyfert galaxies.

We present a NICMOS 1.6 micron imaging isophotal study of 27 early-type galaxies. Our selected sample is not as biased against dusty galaxies as the visible WFPC samples observed with HST. Core galaxies have reduced ellipticity and boxiness near and within their core or break radius. This supports a core formation mechanism which mixes or scatters stars such as scattering caused by a binary black hole. We find no correlation between core properties and dust mass or X-ray luminosity, suggesting that processes determining the current gas content (e.g., such as minor mergers and cooling flows) are unrelated to processes occurring during core formation. Core galaxies exist in a variety of environments ranging from poor groups to large clusters. A combined sample suggests that galaxy groups may harbor more luminous power law galaxies than clusters such as Virgo and Fornax.

Upper limits for the disk mass-to-light ratios for the low surface brightness galaxies UGC 6614 and F568-6 (Malin 2) are estimated by considering the minimum velocity perturbations in the HI velocity field that should result from the spiral structure observed in the R band images. The weak observed response in the $\phi$ velocity component limits the mass-to-light ratios of the disk within a scale length to M/L <~ 3 and 6 for UGC 6614 for F568-6 respectively (in solar units) based upon azimuthal variations observed in the R band images. These limits are sufficiently strong to require a significant dark matter component even in the central regions of these galaxies. Our limits furthermore imply that this dark matter component cannot be in the form of a cold disk since a cold disk would necessarily be involved in the spiral structure. However, a more massive disk could be consistent with the observations because of a non-linear gas response or if the gas is driven by bar-like distortions instead of spiral structure. To produce the large observed arm/interarm HI density variations it is likely that the spiral arm potential perturbation is sufficiently strong to produce shocks in the gas. For a forcing that is greater than 2% of the axisymmetric force, M/L >~ 1 is required in both galaxies in the outer regions. These lower limits imply that the stellar surface density is at least of the same order as the gas surface density. This is consistent with the large scale morphology of the spiral structure, and the stability of the gas disk, both which suggest that a moderate stellar component is required to produce the observed spiral structure. (Shortened abstract)

In radio galaxies such as M84 dust features tend to be nearly perpendicular to radio jets yet are not aligned with the galaxy isophotes. The timescale for precession in the galaxy is short suggesting that an alternative mechanism causes the gas disk to be misaligned with the galaxy. In M84 we estimate the pressure on the disk required to overcome the torque from the galaxy and find that it is small compared to the thermal pressure in the hot ambient ISM estimated from the X-ray emission. We therefore propose that pressure gradients in a jet associated hot interstellar medium exert a torque on the gas disk in M84 causing it to be misaligned with the galaxy isophotal major axis. This mechanism could explain the connection between gas disk angular momentum and jet axes in nearby radio galaxies. By integrating the light of the galaxy through a warped gas and dust disk we find that the geometry of gas disk in M84 is likely to differ that predicted from a simple precession model. We find that the morphology of the gas disk in M84 is consistent with a warped geometry where precession is caused by a combination of a galactic torque and a larger torque due to pressure gradients in the ambient X-ray emitting gas. Precession occurs at an axis between the jet and galaxy major axis, but nearer to the jet axis implying that the pressure torque is 2-4 times larger than the galactic torque. A better model to the morphology of the disk is achieved when precession takes place about an elliptical rather than circular path. This suggests that the isobars in the hot medium are strongly dependent on angle from the jet axis.

We construct near-infrared to visible broad band NICMOS/WFPC color maps for 4 early-type Seyfert galaxies with S-shaped or one-sided ionization cones. We find that dust lanes are near or connected to many of the features seen in the [OIII] and H$\alpha$+[NII] line emission maps. This suggests that much of the structure of line emission in these ionization cones is determined by the distribution of ambient dense galactic gas. Spiral arms, dust lanes caused by bars, or gaseous warps provide dense gas which when illuminated by a conical beam of ultraviolet photons can result in the complicated line emission morphologies observed.

We have used ISO to observe at 12$\mu$m seven E+A galaxies plus an additional emission line galaxy, all in the Coma cluster. E+A galaxies lacking narrow emission lines have 2.2$\mu$m to 12$\mu$m flux density ratios or limits similar to old stellar populations (typical of early-type galaxies). Only galaxies with emission lines have enhanced 12$\mu$m flux density. Excess 12$\mu$m emission is therefore correlated with the presence of on-going star formation or an active galactic nucleus (AGN). By comparing the current star formation rates with previous rates estimated from the Balmer absorption features, we divide the galaxies into two groups: those for which star formation has declined significantly following a dramatic peak $\sim$ 1 Gyr ago; and those with a significant level of ongoing star formation or/and an AGN. There is no strong difference in the spatial distribution on the sky between these two groups. However, the first group has systemic velocities above the mean cluster value and the second group below that value. This suggests that the two groups differ kinematically. Based on surveys of the Coma cluster in the radio, the IRAS sources, and galaxies detected in H$\alpha$ emission, we sum the far infrared luminosity function of galaxies in the cluster. We find that star formation in late type galaxies is probably the dominant component of the Coma cluster far infrared luminosity. The presence of significant emission from intracluster dust is not yet firmly established. The member galaxies also account for most of the far infrared output from nearby rich clusters in general.

I summarize a series of studies where the underlying theme is to exploit non-axisymmetry disk structures such as bars and spiral arms. By measuring the strength of gravitational forces from these structures we can measure the mass in them. The advantage of this approach is that mass to light ratios can be constrained in a way independent of the dark matter profile. By examining a variety of types of galaxies it appears that a correlation between morphology and dark matter fraction emerges.

We discuss the possibility that gaseous giant planets drive strong outflows during early phases of their formation. We consider the range of parameters appropriate for magneto-centrifugally driven stellar and disk outflow models and find that if the proto-Jovian planet or accretion disk had a magnetic field of >~ 10 Gauss and moderate mass inflow rates through the disk of less than 10^-7 M_J/yr that it is possible to drive an outflow. Estimates based both on scaling from empirical laws observed in proto-stellar outflows and the magneto-centrigugal disk and stellar+disk wind models suggest that winds with mass outflow rates of 10^-8 M_J/yr and velocities of order ~ 20 km/s could be driven from proto-Jovian planets. Prospects for detection and some implications for the formation of the solar system are briefly discussed.

NGC 2915 is a remarkable blue compact dwarf galaxy which contains an inner HI bar and strong spiral structure in its HI disk extending well beyond its optical disk. We propose here that NGC 2915 has a dark or faint (not yet observed) massive disk component. By massive we mean of the same order as the HI gas mass which ranges from 0.5 -- 10 Msol/pc^2. We present three different dynamical arguments based on bar and spiral structure that an additional disk mass component is required. The additional mass requirements from these three dynamical arguments depend upon the distance of the galaxy in the same way. If the galaxy were at a distance twice that estimated then the additional mass required would be negligible. Virgocentric models show that then its radial velocity would be 250-500 km/s below its local bulk flow value. Local flow models suggest that this is unlikely, however they also do not exclude it at high confidence levels. Detection of a quiescent or faint stellar disk population would be an exciting prospect since it would be strong evidence for a previous epoch of star formation in a previously low surface brightness galaxy. Then NGC 2915 would be an example of a galaxy which has faded and so could have been one of the faint blue galaxy population observed at moderate redshift.

We consider seven distant galaxies with clearly evident spiral structure from HST images. Three of these were chosen from Vogt et al. (1996) (VFP) and have measured rotational velocities. Five were chosen from the Medium Deep Survey and are studied in Sarajedini et al. 1996 (SGGR), and one galaxy is found in both papers. We place upper limits on their mass-to-light ratios (M/L) by computing M/L_B for a maximal disk. We find that these galaxies have maximal disk mass-to-light ratios M/L_B = 1.5 - 3.5 M_sol/L_Bsol at the low end, but within the range seen in nearby galaxies. The mass-to-light ratios are low enough to suggest that the galaxies contain a young, rapidly formed stellar population. By using a Toomre stability criterion for formation of spiral structure, we place constraints on the ratio of M/L to the stellar velocity dispersion. If these galaxies have maximal disks they would have to be nearly unstable so as to have small enough velocity dispersions that their disks are not unrealistically thick. This suggests that there is a substantial amount of dark matter present in the luminous regions of the galaxy.

In M84 dust features are not aligned with the galaxy isophotes but are perpendicular to the radio jets. We estimate the timescale that the gas disk should precess in the galaxy. Since this timescale is short ($\sim 2\times 10^7$ years at 100 pc) we consider mechanisms that could cause the gas disk to be misaligned with the galaxy. One possibility is that a cooling flow replenishes the disk on this timescale and at angles that vary with radius. Another possibility is that there is a pressure on the disk that is large enough to overcome the torque from the galaxy. We estimate this pressure and find that it can be provided by ram pressure from an energetic interstellar medium that is consistent with velocity dispersions observed in ionized gas and densities estimated from the X-ray emission. We therefore propose that an AGN associated energetic interstellar medium is responsible for causing the gas disk in M84 to be misaligned with the galaxy isophotal major axis for $r \lta 600$ pc. We propose that AGN associated outflows or kinetic motion in low density media could be responsible for jet/disk alignments observed on 100 pc scales in nearby radio galaxies.

We consider 9 fields of low ecliptic latitude observed at 12 microns centered on QSOs. In these fields we detect 7 additional background sources at the level of 0.4-4mJy in 70 square arcminutes. 4 of these sources correspond to galaxies seen in the Digitized Sky Survey. 3 of these sources are not observed in this survey. The 4 sources with optical counterparts are low redshift objects with low infrared to optical flux ratios and so low inferred infrared luminosities L_FIR <~ 10^9 L_solar. Statistical arguments suggest that the unidentified sources are more likely to be distant L_FIR > 10^10 L_solar starbursting galaxies, rather than asteroids or late-type stars. 2 of our QSO fields are likely to have starbursting galaxies at moderately high redshift near the QSO (in angular scale). Since the lensing optical depth towards these QSOs is expected to peak at moderate redshift there is some possibility that these galaxies are in clusters that lense the QSOs. If none of our sources are asteroids then the number of small (km) sized main belt asteroids cannot be larger than 500,000 per absolute magnitude bin.

Outer rings are located at the greatest distance from the galaxy center of any feature resonant with a bar. Because of their large scale, their morphology is sensitive to the distribution of the dark matter in the galaxy. We introduce here how study of these rings can constrain the mass-to-light ratio of the bar, and so the percentage of dark matter in the center of these galaxies. We compare periodic orbits integrated in the ringed galaxy NGC 6782 near the outer Lindblad resonance to the shape of the outer ring. The non-axisymmetric component of the potential resulting from the bar is derived from a near-infrared image of the galaxy. The axisymmetric component is derived assuming a flat rotation curve. We find that the pinched non-self-intersecting periodic orbits are more elongated for higher bar mass-to-light ratios and faster bars. The inferred mass-to-light ratio of the bar depends on the assumed inclination of the galaxy. With an assumed galaxy inclination of i=41 degrees, for the orbits to be consistent with the observed ring morphology the mass-to-light ratio of the bar must be high, greater than 70% of a maximal disk value. For i=45 degrees, the mass-to-light ratio of the bar is $75\pm 15%$ of the maximal disk value. Since the velocity field of these rings can be used to constrain the galaxy inclination as well as which periodic orbit is represented in the ring, further study will yield tighter constraints on the mass-to-light ratio of the bar. If a near maximal disk value for the bar is required, then either there would be little dark matter within the bar, or the dark matter contained in the disk of the galaxy would be non-axisymmetric and would rotate with the bar.

We present near-Infrared images of low surface brightness galaxies: H band images of UGC 6614 and F568-6 (Malin 2). The H band images show spiral structure so we are confident that we are seeing the galaxy disks as well as the central bulges. The optical-IR colors of these galaxies, $R-H = 2.2\pm 0.2$, $B-H = 3.5-4.2$, are extremely red and are similar to those of many S0 and elliptical galaxies. This represents strong evidence for a component of old stars. Recent studies find that low surface brightness galaxies span a wide range of morphologies and colors. In this context UGC 6614 and F568-6 could be examples of low surface brightness galaxies which have undergone a past epoch of more vigorous star formation.

We measure the extinction law in a galaxy's spiral arm and interarm regions using a visual and infrared (BVRJHK) imaging study of the interacting galaxies NGC 2207 and IC 2163. This is an overlapping spiral galaxy pair in which NGC 2207 partially occults IC 2163. This geometry enables us to directly measure the extinction of light from the background galaxy as it passes through the disk of the foreground galaxy. We measure the extinction as a function of wavelength, and find that there is less extinction in the optical bands than expected from a normal Galactic extinction law. This deviation is significantly larger in the interarm region than in the spiral arm. The extinction curve in the spiral arm resembles a Milky Way $R_V=5.0$ dust model and the interarm extinction curve is flatter (``greyer'') still. We examine the effect of scattering of background galaxy light into the line of sight and find that it is negligible. We also examine the effect of an unresolved patchy dust distribution using a simple two-component dust model as well as the clumpy dust model of Witt & Gordon (1996). Both models clearly demonstrate that an unresolved patchy dust distribution can flatten the extinction curve significantly. When fit to the data, both models suggest that the observed difference between the arm and interarm extinction curves is caused by the interarm region of NGC 2207 having a higher degree of dust patchiness (density ratio between high-density and low-density phases) than the spiral arm region. We note that an unresolved patchy dust distribution will cause us to underestimate the average column depth of gas in a galaxy if based solely on the visual extinction. It is much better to use the infrared extinction for this purpose.

We report on the discovery of a boxy/peanut shaped bulge in the highly inclined barred Seyfert 2 galaxy NGC~7582. The peanut shape is clearly evident in near infrared $JHK$ images but obscured by extinction from dust in visible $BVR$ images. This suggests that near infrared imaging surveys will discover a larger number of boxy/peanut morphologies than visible surveys, particularly in galaxies with heavy extinction such as NGC~7582. The bulge in NGC~7582 exhibits strong boxiness compared to other boxy/peanut shaped bulges. If the starburst was mediated by the bar, then it is likely that the bar formed in less than a few bar rotation periods or a few $\times 10^8$ years ago. If the bar also caused the peanut, then the peanut would have formed quickly; on a timescale of a few bar rotation periods.

We compare stellar orbits in an early (NGC 4314) and a late-type barred galaxy (NGC 1073). We find that these bars are cold in the sense that the majority of stars can be described as being nearby to periodic orbits. We place limits of (65km/s)^2 and (50km/s)^2 on the components of the diagonalized velocity dispersion ellipsoid matrix for stars in the bars of NGC 4314 and NGC 1073 respectively. Both bars end near the inner 4:1 Lindblad (ultraharmonic) resonance. We conjecture that a large class of bars end near this resonance. The morphology of the bars depends on the strength of the 4:1 resonance, which is strong in the early-type barred galaxy and weak in the other. This results in a flat bar major axis surface brightness profile for the early-type bar and an exponential profile for the late-type bar.

We describe a technique that uses radial color gradients in disk galaxies to detect the presence of bulk non-circular motion or elliptical orbits. In a disk galaxy with both a radial color gradient and non-circular motion, isochromes or iso-color contours should follow the shape of closed stellar orbits, and the ellipticity of the isophotes should vary as a function of wavelength. A difference in the ellipticity of isochromes and the isophotes can be used to detect the presence of non-circular motion. A model galaxy is constructed which demonstrates this phenomenon. The difference between isochrome and isophote ellipticity is directly related to the ellipticity of the potential. This provides a new way to measure the ellipticity of the dark matter in the outer parts of galaxies. As an example, we apply this technique to two dwarf galaxies NGC 1800 and NGC 7764. We detect a bar in NGC 1800 which has only previously been suggested from the HI velocity field. In NGC 7764 there is no color gradient along its bar so we cannot detect non-circular motion in this region; however ellipticities observed in a star forming ring at the end of the bar are consistent with this ring being located near the corotation resonance.

We present $B,V,R,J,H,K$ broad--band images of the barred galaxy NGC 1097. The optical and infrared colors maps trace the location of the major dust features. The dust lanes are relatively free of star formation and have low opacity. The depth of the dust lanes decreases as a function of wavelength so that they are deepest at $B$ and shallower at $H$ and $K$ where they are difficult to detect. We find that the dust lane on the north-west side of the galaxy has colors consistent with those of a dust screen with opacity derived from the galactic extinction law. Thus, it must be located on the near side of the galaxy. The colors in the dust lane on the south-east side of the galaxy, on the other hand, suggest that there are stars both in front of and behind the dust, consistent with this dust lane being located on the far side of the galaxy. Near the inner spiral ring (at $\sim 17'' = 1.4$kpc from the nucleus), the dust lane becomes double peaked with peak to valley extinction ratio similar to that observed in radio continuum by Ondrechen and van der Hulst (1983). This suggests that if the cold diffuse component of the interstellar medium (traced by the dust) moves with the cosmic ray electron component then the magnetic field pressure is not a significant force in the shocks. Alternatively if the two phases do not move together a more detailed comparison may show differences in the shock shapes. From the $K$ images, we estimate that the total stellar mass of the star forming spiral ring at a radius of $\sim$ 1kpc from the nucleus is of the same order as the molecular gas mass in the ring. The infrared images show a short bar inside the inner spiral ring. In the principal plane of the galaxy, the short bar is not perpendicular to the prominent outer bar. This suggests that there are torques between the inner bar and the spiral

We present visible and near-infrared images of the amorphous dwarf galaxy NGC 1705 in the $B$, $V$, $R$, $J$, $H$, and $K$ bands. Optical and near infrared colors in the galaxy are consistent with a composite population with colors similar to those of the bluest Sc nuclei. The contribution of young stars to the disk population of the galaxy decreases outwards from the center. In the near infrared images we detect three unresolved point sources, two of which have colors and luminosities comparable to intermediate age LMC clusters.

In this paper, by considering a simple fluid model, we investigate the role of phase transitions in the ISM on the galaxy scale gas dynamics. Cooling and heating timescales in the ISM are typically shorter than typical galactic rotation timescales, so the individual phases in the ISM can be assumed to be in temperature equilibrium with the radiation field. Using this approximation we can construct an equation of state which depends upon the average density and mass fractions in the individual phases. Previous studies suggest that there is an equilibrium phase fraction as a function of pressure. We incorporate evolution towards this equilibrium state as a relaxation term with a time to obtain equilibrium $\tau$. We derive a condition in terms of a critical Mach number when one dimensional shocks should be continuous. For small values of the relaxation time $\tau$ we show that the relaxation term acts like a viscosity. We show with one dimensional simulations that increasing $\tau$ causes shocks to become smoother. This work suggests that phase changes can strongly effect the gas dynamics of the ISM across spiral arms and bars.

We present near-infrared images of NGC 4314 in the $J$,$H$, and $K$ bands. Colors in these bands are constant across the bar of NGC 4314 indicating that the mass-to-light ratio is constant over much of the galaxy. We have developed an efficient procedure for mapping the gravitational potential of a non-axisymmetric galaxy based on its appearance in the near-infrared using a Fourier Transform method. We account for the thickness of the galaxy's stellar disk by convolving its $K$ image with a function that depends upon the vertical stellar scale height. We apply this procedure to NGC 4314 and find that prograde stellar orbits integrated in the potential are consistent with the shape of the bar for a disk vertical scale height of $h = 350 \pm 100$ pc and a corotation radius of $r=3.5 \pm 0.5$ kpc. We find that the inner $m=4$ Lindblad resonance is approximately at the location of knob-like features seen near the end of the bar. The shape of nearly closed orbits at this radius are lozenge or diamond shaped and have low speeds along the major axis of the bar. The low speeds imply an increase in density along the major axis of the bar. The angle of the major axis of the $m=4$ component of the surface density begins to twist at the location of its resonance, and the angle of the $m=2$ component begins to twist at a larger radius near the corotation radius.

We present images of the barred galaxy NGC 7479 in the optical and near infrared broad--bands $B,V,R,J,H,K$, and images in H$\alpha$+[N II] and CO(J=1$\rightarrow$0) emission. The H$\alpha$ and CO emission in the bar are coincident and confined to narrow linear features that are offset from the center of the bar as observed in the near infrared. We estimate the gravitational potential from the $K$ image, which provides an estimate of the torque on the gas at the position of the CO emission in the bar. We find that the implied gas inflow velocity derived from the torque is $10-20$ km s$^{-1}$. Our inflow velocity is independent of the large streaming motions which can be observed in CO and HI.