We present an exquisite, 30-min cadence Kepler (K2) light curve of the Type Ia supernova (SN Ia) 2018oh (ASASSN-18bt), starting weeks before explosion, covering the moment of explosion and the subsequent rise, and continuing past peak brightness. These data are supplemented by multi-color Pan-STARRS1 and CTIO 4-m DECam observations obtained within hours of explosion. The K2 light curve has an unusual two-component shape, where the flux rises with a steep linear gradient for the first few days, followed by a quadratic rise as seen for typical SNe Ia. This "flux excess" relative to canonical SN Ia behavior is confirmed in our $i$-band light curve, and furthermore, SN 2018oh is especially blue during the early epochs. The flux excess peaks 2.14$\pm0.04$ days after explosion, has a FWHM of 3.12$\pm0.04$ days, a blackbody temperature of $T=17,500^{+11,500}_{-9,000}$ K, a peak luminosity of $4.3\pm0.2\times10^{37}\,{\rm erg\,s^{-1}}$, and a total integrated energy of $1.27\pm0.01\times10^{43}\,{\rm erg}$. We compare SN 2018oh to several models that may provide additional heating at early times, including collision with a companion and a shallow concentration of radioactive nickel. While all of these models generally reproduce the early K2 light curve shape, we slightly favor a companion interaction, at a distance of $\sim$$2\times10^{12}\,{\rm cm}$ based on our early color measurements, although the exact distance depends on the uncertain viewing angle. Additional confirmation of a companion interaction in future modeling and observations of SN 2018oh would provide strong support for a single-degenerate progenitor system.

For decades optical time-domain searches have been tuned to find ordinary supernovae, which rise and fall in brightness over a period of weeks. Recently, supernova searches have improved their cadences and a handful of fast-evolving luminous transients (FELTs) have been identified. FELTs have peak luminosities comparable to Type Ia supernovae, but rise to maximum in $<10$ days and fade from view in $<$month. Here we present the most extreme example of this class thus far, KSN2015K, with a rise time of only 2.2 days and a time above half-maximum ($t_{1/2}$) of only 6.8 days. Here we show that, unlike Type Ia supernovae, the light curve of KSN2015K was not powered by the decay of radioactive elements. We further argue that it is unlikely that it was powered by continuing energy deposition from a central remnant (a magnetar or black hole). Using numerical radiation hydrodynamical models, we show that the light curve of KSN2015K is well fit by a model where the supernova runs into external material presumably expelled in a pre-supernova mass loss episode. The rapid rise of KSN2015K therefore probes the venting of photons when a hypersonic shock wave breaks out of a dense extended medium.

We discovered two transient events in the Kepler field with light curves that strongly suggest they are type II-P supernovae. Using the fast cadence of the Kepler observations we precisely estimate the rise time to maximum for KSN2011a and KSN2011d as 10.5$\pm 0.4$ and 13.3$\pm 0.4$ rest-frame days respectively. Based on fits to idealized analytic models, we find the progenitor radius of KSN2011a (280$\pm 20$ R$_\odot$) to be significantly smaller than that for KSN2011d (490$\pm 20$ R$_\odot$) but both have similar explosion energies of 2.0$\pm 0.3\times 10^{51}$ erg. The rising light curve of KSN2011d is an excellent match to that predicted by simple models of exploding red supergiants (RSG). However, the early rise of KSN2011a is faster than the models predict possibly due to the supernova shockwave moving into pre-existing wind or mass-loss from the RSG. A mass loss rate of $10^{-4}$ M$_\odot$ yr$^{-1}$ from the RSG can explain the fast rise without impacting the optical flux at maximum light or the shape of the post-maximum light curve. No shock breakout emission is seen in KSN2011a, but this is likely due to the circumstellar interaction suspected in the fast rising light curve. The early light curve of KSN2011d does show excess emission consistent with model predictions of a shock breakout. This is the first optical detection of a shock breakout from a type II-P supernova.

We report on candidate active galactic nuclei (AGN) discovered during the monitoring of $\sim$500 bright (r < 18 mag) galaxies over several years with the Kepler Mission. Most of the targets were sampled every 30 minutes nearly continuously for a year or more. Variations of 0.001 mag and often less could be detected reliably. About 4.0% (19) of our random sample continuously fluctuated with amplitudes increasing with longer timescales, but the majority are close to the limits of detectability with Kepler. We discuss our techniques to mitigate the long term instrumental trends in Kepler light curves and our resulting structure function curves. The amplitudes of variability over four month periods, as seen in the structure functions and PSDs, can dramatically change for many of these AGN candidates. Four of the candidates have features in their Structure Functions that may indicate quasi-periodic behavior, although other possibilities are discussed.

We develop Bayesian statistical methods for discovering and assigning probabilities to non-random (e.g., physical) stellar companions. These companions are either presently bound or were previously bound. The probabilities depend on similarities in corrected proper motion parallel and perpendicular to the brighter component's motion, parallax, and the local phase-space density of field stars. Control experiments are conducted to understand the behavior of false positives. The technique is applied to the Hipparcos Catalogue within 100 pc. This is the first all-sky survey to locate escaped companions still drifting along with each other. In the <100 pc distance range, ~220 high probability companions with separations between 0.01 - 1 pc are found. The first evidence for a population (~300) of companions separated by 1 - 8 pc is found. We find these previously unnoticed naked-eye companions (both with V<6th mag): Capella & 50 Per, delta Vel & HIP 43797, Alioth (epsilon UMa), Megrez (delta UMa) & Alcor (80 UMA), gamma & tau Cen, phi Eri & eta Hor, 62 & 63 Cnc, gamma & tau Per, zeta & delta Hya, beta01, beta02 & beta03 Tuc, N Vel & HIP 47479, HIP 98174 & HIP 97646, 44 & 58 Oph, s Eri & HIP 14913, and pi & rho Cep. High probability fainter companions (V > 6th mag) of primaries with V < 4th mag are found for: Fomalhaut (alpha PsA), gamma UMa, alpha Lib, Alvahet (iota Cephi), delta Ara, beta Ser, iota Peg, beta Pic, kappa Phe and gamma Tuc.

Astrometry can detect rocky planets in a broad range of masses and orbital distances and measure their masses and three-dimensional orbital parameters, including eccentricity and inclination, to provide the properties of terrestrial planets. The masses of both the new planets and the known gas giants can be measured unambiguously, allowing a direct calculation of the gravitational interactions, both past and future. Such dynamical interactions inform theories of the formation and evolution of planetary systems, including Earth-like planets. Astrometry is the only technique technologically ready to detect planets of Earth mass in the habitable zone (HZ) around solar-type stars within 20 pc. These Earth analogs are close enough for follow-up observations to characterize the planets by infrared imaging and spectroscopy with planned future missions such as the James Webb Space Telescope (JWST) and the Terrestrial Planet Finder/Darwin. Employing a demonstrated astrometric precision of 1 microarcsecond and a noise floor under 0.1 micro-arcseconds, SIM Lite can make multiple astrometric measurements of the nearest 60 F-, G-, and K-type stars during a five-year mission. SIM Lite directly tests theories of rocky planet formation and evolution around Sun-like stars and identifies the nearest potentially habitable planets for later spaceborne imaging, e.g., with Terrestrial Planet Finder and Darwin. SIM was endorsed by the two recent Decadal Surveys and it meets the highest-priority goal of the 2008 AAAC Exoplanet Task Force.

Stellar mass plays a central role in our understanding of star formation and aging. Stellar astronomy is largely based on two maps, both dependent on mass, either indirectly or directly: the Hertzprung-Russell Diagram (HRD) and the Mass-Luminosity Relation (MLR). The extremes of both maps, while not terra incognita, are characterized by large uncertainties. A precise HRD requires precise distance obtained by direct measurement of parallax. A precise MLR requires precise measurement of binary orbital parameters, with the ultimate goal the critical test of theoretical stellar models. Such tests require mass accuracies of ~1%. Substantial improvement in both maps requires astrometry with microsecond of arc measurement precision. Why? First, the 'tops' of both stellar maps contain relatively rare objects, for which large populations are not found until the observing horizon reaches hundreds or thousands of parsecs. Second, the 'bottoms' and 'sides' of both maps contain stars, either intrinsically faint, or whose rarity guarantees great distance, hence apparent faintness. With an extensive collection of high accuracy masses that can only be provided by astrometry with microsecond of arc measurement precision, astronomers will be able to stress test theoretical models of stars at any mass and at every stage in their aging processes.

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.

The fact that dark matter (DM), thus far, has revealed itself only on scales of galaxies and larger, again thrusts onto astrophysics the opportunity and the responsibility to confront the age old mystery "What is the nature of matter?" By deriving basic data on the nature of DM - e.g., mass of its particle(s), present mean temperature, distribution in galaxies and other structures in the universe, and capacity for dissipational collapse - we will be uncovering the properties of the dominant species of matter in the universe and significantly extending the standard models of particle physics. Determining the mass of the DM particle to an order of magnitude would help to sort out the particle family to which it (or they) belongs. Beyond mass, there are issues of stability. The DM particle may be unstable with a measurable half-life, or it may become unstable after absorbing a certain amount of energy from collisions. In both cases it would contribute to the present hot dark matter component. Some key parameters of DM can most accurately be measured in the very nearby universe because DM dominates the mass in the outer Milky Way (MW), in other galaxies in the Local Group, and in the Local Group in its entirety. The presence and amount of DM can be quantified by study of dynamical processes observable in fine detail within these entities. Precise measurements of 3-D velocities for stars, coherent star streams, and stars in satellite stellar systems out to the edge of the Galaxy can reveal "what is the shape, orientation, density law, and lumpiness of the dark matter halo" as well as "what is the total mass of the Galaxy?"

Precision astrometry at microarcsecond accuracy has application to a wide range of astrophysical problems. This paper is a study of the science questions that can be addressed using an instrument that delivers parallaxes at about 4 microarcsec on targets as faint as V = 20, differential accuracy of 0.6 microarcsec on bright targets, and with flexible scheduling. The science topics are drawn primarily from the Team Key Projects, selected in 2000, for the Space Interferometry Mission PlanetQuest (SIM PlanetQuest). We use the capabilities of this mission to illustrate the importance of the next level of astrometric precision in modern astrophysics. SIM PlanetQuest is currently in the detailed design phase, having completed all of the enabling technologies needed for the flight instrument in 2005. It will be the first space-based long baseline Michelson interferometer designed for precision astrometry. SIM will contribute strongly to many astronomical fields including stellar and galactic astrophysics, planetary systems around nearby stars, and the study of quasar and AGN nuclei. SIM will search for planets with masses as small as an Earth orbiting in the `habitable zone' around the nearest stars using differential astrometry, and could discover many dozen if Earth-like planets are common. It will be the most capable instrument for detecting planets around young stars, thereby providing insights into how planetary systems are born and how they evolve with time. SIM will observe significant numbers of very high- and low-mass stars, providing stellar masses to 1%, the accuracy needed to challenge physical models. Using precision proper motion measurements, SIM will probe the galactic mass distribution and the formation and evolution of the Galactic halo. (abridged)

We investigate how the uncertainty on the Hubble constant (H_0) affects the uncertainty in the Equation of State (EOS) of Dark Energy and the total density of the Universe (Omega_tot). We use the approximate relations between the cosmological parameters [Spergel etal (2007)] and use error-propagation to estimate the effects of improving the CMB parameters and H_0 on the EOS of Dark Energy (DE). First we assume that the additional data does not improve significantly, but decrease the error on H_0 by a factor <~10. Second, we allow improved additional data but current H_0 errors (i.e., the DE Task Force case). In the 1st scenario, improvements of the CMB parameters hardly change the accuracy of the EOS and Omega_tot, unless H_0 can be measured with an accuracy of a few %. We find that a combination of moderate improvements for both H_0 and other data significantly constrains the evolution of dark energy, but at a reduced cost. We review several methods (and their strengths and weaknesses) that might yield extra-galactic distances with errors of about 1%. We review: the Velocity Field method, two Maser methods, two Light Echo techniques, the Binary Star method, and the Rotational Parallax (RP) technique. Because these methods substantially rely on geometry rather than astrophysics or cosmology, their results are quite robust. We focus on the advantages of the RP technique which can provide single-step, bias-free distances to nearby spirals. These distances can be used to improve the zero-point for other methods which in turn allow for a much improved H_0 errors. Achieving an accuracy of ~2% in the distances to M31, M33 and the LMC by the RP method requires proper motions from future astrometric missions (SIM, GAIA and OBSS, or the SKA).

We propose a MIDEX-class space mission with the goal to find and characterize roughly 10,000 transiting planets. When transits occur, a much more detailed characterization of the planet is possible and so a large data base of transiting planets will provide planets with a large range in periods and radii for follow-up studies. Our survey will be all-sky and focused on stars brighter than V=14.8. Down to V=12, LEAVITT will be able to detect Neptune-sized objects. Because of it's high cadence, LEAVITT is about 100 times more sensitive at detecting transits than GAIA, while it will find more than 20 times as many transits as KEPLER. LEAVITT has multi-band photometric capability implemented via a low-resolution dispersive element which can obtain 0.2% (2 mmag) photometry down to V=14.8. LEAVITT's high multi-band photometric accuracy reduces the number of false-positives significantly.

The astrometric signature imposed by a planet on its primary increases substantially towards longer periods (proportinal to P^2/3), so that long-period planets can be more easily detected, in principle. For example, a one Solar-mass (M_Sun) star would be pulled by roughly 1 mas by a one Jupiter-mass (M_J) planet with a period of one-hundred years at a distance of 20 pc. Such position accuracies can now be obtained with both ground-based and space-based telescopes. The difficulty was that it often takes many decades before a detectable position shift will occur. However, by the time the next generation of astrometric missions such as SIM will be taking data, several decades will have past since the first astrometric mission, HIPPARCOS. Here we propose to use a new astrometric method that employs a future, highly accurate SIM Quick-Look survey and HIPPARCOS data taken twenty years prior. Using position errors for SIM of 4 muas, this method enables the detection and characterization of Solar-system analogs (SOSAs) with periods up to 240 (500) years for 1 (10) M_J companions. Because many tens of thousands nearby stars can be surveyed this way for a modest expenditure of SIM time and SOSAs may be quite abundant, we expect to find many hundreds of extra-solar planets with long-period orbits. Such a data set would nicely complement the short-period systems found by the radial-velocity method. Brown dwarfs and low-mass stellar companions can be found and characterized if their periods are shorter than about 500 years. This data set will provide invaluable constraints on models of planet formation, as well as a database for systems where the location of the giant planets allow for the formation of low-mass planets in the habitable zone. [Abridged]

We have designed and constructed a ``dispersed Fourier Transform Spectrometer'' (dFTS), consisting of a conventional FTS followed by a grating spectrometer. By combining these two devices, we negate a substantial fraction of the sensitivity disadvantage of a conventional FTS for high resolution, broadband, optical spectroscopy, while preserving many of the advantages inherent to interferometric spectrometers. In addition, we have implemented a simple and inexpensive laser metrology system, which enables very precise calibration of the interferometer wavelength scale. The fusion of interferometric and dispersive technologies with a laser metrology system yields an instrument well-suited to stellar spectroscopy, velocimetry, and extrasolar planet detection, which is competitive with existing high-resolution, high accuracy stellar spectrometers. In this paper, we describe the design of our prototype dFTS, explain the algorithm we use to efficiently reconstruct a broadband spectrum from a sequence of narrowband interferograms, and present initial observations and resulting velocimetry of stellar targets.

I analyze several catalogs of known visual and spectroscopic binaries and conclude that a large number of binaries is missing in current catalogs. Samples of the best studied (nearby and bright) stars indicate that the true binary fraction may be as high as 95%. A preliminary analysis indicates that these binaries can a ect the astrometry signicantly.

The Oort constants describe the local variations of the stellar streaming field. Classically, they are determined from stellar proper motions. We discuss problems arising in this procedure. A large, hitherto overlooked, source of systematic error arises from longitudinal variations of the mean stellar parallax. Together with the solar reflex motion, these variations create contributions to the longitudinal proper motions mu_l*(l) that are indistinguishable from the Oort Constants at the 20% level. Fortunately, we can correct for this "mode mixing" using the latitudinal proper motions mu_b(l). We use ~10^6 stars from the ACT/Tycho-2 catalogs with proper motion errors of ~ 3 mas/yr. We find significant deviations from expectations based on a smooth axisymmetric equilibrium disk, in particular non-zero C for old red giant stars. We also find variations of the Oort constants with the asymmetric drift of the sub-sample. These correlations are different in nature than those expected for an axisymmetric Galaxy. The most reliable tracers for the ``true'' Oort constants are red giants, which are old enough to be in equilibrium and distant enough to be unaffected by possible local anomalies. For these stars we find, A ~ 16, B ~- 17, A - B ~ 33, and C ~- 10 km/s/kpc with internal errors of about 1-2 and external error of perhaps the same order. These values are consistent with our knowledge of the Milky Way (flat rotation curve and Omega=A - B ~ 28 +/- 2). Based on observations made with the ESA Hipparcos astrometry satellite. (Abridged)

Three planned astrometry survey satellites, FAME, DIVA, and GAIA, all aim at observing magnitude-limited samples. We argue that substantial additional scientific opportunities are within the reach of these mission if they devote a modest fraction of their catalogs to selected targets that are fainter than their magnitude limits. We show that the addition of ~10^6 faint (R>15) targets to the 40 10^6 object FAME catalog can improve the precision of the reference frame by a factor 2.5, to 7 micro-as/yr, increase Galactocentric distance at which halo rotation can be precisely (2 km/s) measured by a factor 4, to 25 kpc, and increase the number of late M dwarfs, L dwarfs, and white dwarfs with good parallaxes by an order of magnitude. In most cases, the candidate quasars, horizontal branch stars, and dim dwarfs that should be observed to achieve these aims are not yet known. We present various methods to identify candidates from these classes, and assess the efficiencies of these methods. The analysis presented here can be applied to DIVA with modest modifications. Application to GAIA should be deferred until the characteristics of potential targets are better constrained.

(Abridged) Axisymmetric models of the Milky Way exhibit strong interrelations between the Galactic constants (R_0 and T_0), the stellar columndensity (S_*) and the shape of the dark matter (DM) halo. Here we present analytical relations that can be used to investigate the effects of the uncertain gaseous velocity dispersion on the HI flaring constraints. The contribution of cosmic rays and magnetic fields to the pressure gradients is small. A significantly flattened dark matter halo is only possible if R_0 <~ 6.8 kpc. If R_0 is larger than ~7 kpc, or T_0 >~ 170 km/s, we can rule out two DM candidates that require a highly flattened DM halo: 1) decaying massive neutrinos; and 2) a disk of cold molecular hydrogen. It is only possible to construct self-consistent models of the Galaxy based on the IAU-recommended values for the Galactic constants in the unlikely case that the the stellar columndensity is smaller than ~18 M_sun/pc^2. If we assume that the halo is oblate and S_* = 35 +/- 5 M_sun/pc^2, R_0 <~ 8 kpc and T_0 <~ 200 km/s. Combining the best kinematical and star-count estimates of S_*, we conclude that: 25 <~ S_* <~ 45 M_sun/pc^2. Kuijken & Gilmore's (1991) determination of the columndensity of matter with |z|<=1.1 kpc is robust and valid over a wide range of Galactic constants. Our mass models show that the DM density in the Galactic centre is uncertain by a factor 1000. In the Solar neighbourhood we find: rho_DM ~0.42 GeV/c^2/cm^3 or (11 +/- 5) mM_sun/pc^3 -- roughly 15% of rho_tot.

We have obtained optical long-slit spectroscopy of the nucleus of M32 using the Space Telescope Imaging Spectrograph aboard the Hubble Space Telescope. The stellar rotation velocity and velocity dispersion, as well as the full line-of-sight velocity distribution (LOSVD), were determined as a function of position along the slit using two independent spectral deconvolution algorithms. We see three clear kinematical signatures of the nuclear black hole: a sudden upturn, at ~0.3 arc seconds from the center, in the stellar velocity dispersions; a flat or rising rotation curve into the center; and strong, non-Gaussian wings on the central LOSVD. The central velocity dispersion is ~130 km/s (Gaussian fit) or ~175 km/s (corrected for the wings). Both the velocity dispersion spike and the shape of the central LOSVD are consistent with the presence of a supermassive compact object in M32 with a mass in the range 2-5 x 10^6 solar masses. These data are a significant improvement on previous stellar kinematical data, making M32 the first galaxy for which the imprint of the black hole's gravitation on the stellar velocities has been observed with a resolution comparable to that of gas-dynamical studies.

Accurate distances to stars and galaxies allow for significant tests of stellar evolution, galaxy formation and evolution and cosmology. NASA's Space Interferometry Mission (SIM) will obtain precision [1-30 microarcsec] astrometry for 10-20k objects brighter than V=20 mag. In this paper we discuss a method to determine the distance of nearby spiral galaxies using the technique of rotational parallaxes, and show that it is possible to achieve distance errors of only a few percent. With such distances at hand, it becomes possible to determine an accurate zero-point for the Tully-Fisher relation, one of the tools to measure distances throughout the universe. The rotational parallax method employs the common orbital motions of a of stars within a galaxy to determine their distance. For stars on circular orbits, the two proper motions and the radial velocity suffice to determine the inclination of the orbit, the rotational velocity and the distance. Several factors complicate the application of this simple technique to real galaxies: 1) the target galaxy may have substantial space-motion, 2) stars in real galaxies are not on circular orbits (e.g. due to spiral-arm streaming motions or induced by random motions), 3) stars in the galaxy are at significantly different distances from the Sun (the near side of M 31 is ~5% closer than the far side). Notwithstanding these complication, distances can be determined to within several percent for the nearest spirals. (abridged)

In order to test the reliability of determinations of the shapes of galaxies' dark matter halos, we have made such measurements for the Milky Way by two independent methods, which make use of the stellar kinematics in the solar neighbourhood and the observed flaring of the Galactic HI layer to estimate the flattening of the Galactic dark halo. These techniques are found to produce a consistent estimate for the halo shape, with a shortest-to-longest axis ratio of q ~ 0.8, but only if one adopts somewhat non-standard values for the distance to the Galactic centre, R_0, and the local Galactic rotation speed, Theta_0. For consistency, one requires values of R_0 < 7.6 kpc and Theta_0 < 190 km/s. Although differing significantly from the current IAU-sanctioned values, these upper limits are consistent with all existing observational constraints. If future measurements confirm these lower values for the Galactic constants, then the validity of the gas layer flaring method will be confirmed. Further, dark matter candidates such as cold molecular gas and massive decaying neutrinos, which predict very flat dark halos with q < 0.2, will be ruled out. Conversely, if the Galactic constants were found to be close to the more conventional values, then there would have to be some systematic error in the methods for measuring dark halo shapes, so the existing modeling techniques would have to be viewed with some scepticism.

The local stellar kinematics of the Milky Way offer a useful tool for studying the rotation curve of the Galaxy. These kinematics -- usually parameterized by the Oort constants A and B -- depend on the local gradient of the rotation curve as well as its absolute value (Theta_0), and the Sun's distance to the Galactic center (R_0). The density of interstellar gas in the Milky Way is shown to vary non-monotonically with radius, and so contributes significantly to the local gradient of the rotation curve. We have therefore calculated mass models for the Milky Way that include this component, and have derived the corresponding radial variation in the Oort constants. Between 0.9R_0 and 1.2R_0, the Oort functions, A(R) and B(R), differ significantly from the general Theta/R dependence. Various previously-inexplicable observations are shown to be consistent with these predictions. These models can explain the 40% difference between the values for 2 A R_0 derived from radial velocity data originating in the inner and outer Galaxy. They also go some way toward explaining the different shapes of the velocity ellipsoids of giant and dwarf stars in the solar neighbourhood. However, a consistent picture only emerges if one adopts small values of R_0 = 7.1 +/- 0.4 kpc and Theta_0 = 184 +/- 8 km/s. With these Galactic constants, the Milky Way's rotation curve declines slowly in the outer Galaxy; V_rot(20 kpc) = 166 kms. Our low value for R_0 agrees well with the only direct determination (7.2 +/- 0.7 kpc, Reid 1993). Using these Galactic constants, we find that the proper motion of Sgr A^* is consistent with the observational constraints. The radial velocities and proper motions of our best fit model are entirely consistent with the radial velocities of Cepheids and the Hipparcos measurements of their proper motions.

The local stellar kinematics of the Milky Way, parameterized by the Oort constants A and B, depend on the local gradient of the rotation curve, its absolute value (Theta_0), and the distance to the Galactic center (R_0). The surface density of interstellar gas in the Milky Way varies non-monotonically with radius, and so contributes significantly to the local gradient of the rotation curve, and the Oort constants. Because of this, the Oort functions A(R) and B(R) differ significantly from the dominant ~Theta_0/R dependence, in the Solar neighborhood and other locations in the Galaxy. These models may explain the ~40% difference between the values for 2 A R_0 derived from radial velocity data originating in the inner and outer Galaxy (Merrifield 1992). Incorporating these local non-linearities explains the significant differences between the Oort constants derived from nearby stars (d <= 1 kpc; Hanson 1987) and distant Cepheids (d=0.5-6 kpc; Feast & Whitelock 1997). However, a consistent picture only emerges if one adopts small values for the Galactic constants: R_0 = 7.1 +/- 0.4 kpc, and Theta_0 = 184 +/- 8 km/s. These values are consistent with most kinematical methods of determining R_0, including the proper motion of Sgr A* (Backer 1996), the direct determination of R_0 using water masers (7.2 +/- 0.7 kpc; Reid 1993), and constraints set by the shape of the Milky Way's dark halo (Olling & Merrifield 1997b).

The distribution of shapes of galaxies' dark halos provides a basic test for models of galaxy formation. To-date, few dark halo shapes have been measured, and the results of different methods appear contradictory. Here, we add to the sample of measured shapes by calculating the flattening of the Milky Way's dark halo based on the manner in which the gas layer in the Galaxy flares with radius. We also test the validity of this technique -- which has already been applied to several other galaxies -- by comparing the inferred halo flattening to that obtained from a stellar-kinematic analysis, which can only be applied to the Milky Way. Both methods return consistent values for the shape of the Milky Way's halo, with a shortest-to-longest axis ratio for the dark matter distribution of q = 0.75 +/- 0.25. However, this consistency is only achieved if we adopt a value of R_0 = 7 +/- 1 kpc for the Sun's distance to the Galactic center. Although this value is smaller than the IAU-sanctioned R_0 = 8.5 kpc, it is quite consistent with current observations. Whatever value of R_0 is adopted, neither method returns halo parameters consistent with a disk-like mass distribution, for which q ~ 0.2. This finding rules out cold molecular gas and decaying massive neutrinos as dark matter candidates.

In a previous paper (Paper I) a method was developed to determine the shapes of dark matter halos of spiral galaxies from the flaring and velocity dispersion of the gas layer. Here I present the results for the almost edge-on Scd galaxy NGC 4244 (Paper II Paper III) and preliminary results for the Milky Way. NGC 4244's dark matter halo is found to be highly flattened with a shortest-to-longest axis ratio of 0.2 +0.3-0.1. If the dark matter is disk-like, the vertical velocity dispersion of the dark matter must be ~20% larger than the measured tangential dispersion in the HI. The flaring of the Milky Way's gas layer, the local column of identified stars and the total column within 1.1 kpc from the plane are consistent with a moderately flattened dark halo (E7-E0) and galactic constants of \RSUN = 7.1 kpc, \VSUN = 180 \kms, while the stellar disk is ~80% of maximal.

In a previous paper (Olling 1995, astro-ph/9505002) a method was developed to determine the shapes of dark matter halos of spiral galaxies from an accurate determination of the rotation curve, the flaring of the gas layer and the velocity dispersion in the HI. Here I report the results for the almost edge-on Scd galaxy NGC 4244 (Olling 1996a,astro-ph/9605110; 1996b, astro-ph/9605111). The observed flaring of the HI beyond the optical disk puts significant constraints on the shape of the dark matter halo, which are almost independent of the stellar mass-to-light ratio. NGC 4244's dark matter halo is found to be highly flattened with a shortest-to-longest axis ratio of 0.2 (-0.1,+0.3). If the dark matter is disk-like, the data presented in this paper imply that the vertical velocity dispersion of the dark matter must be 10% - 30% larger than the measured tangential dispersion in the HI.

I present sensitive high resolution HI VLA B+C+D array observations of the edge-on Scd galaxy NGC 4244. The gas layer of NGC 4244 is rather symmetric in all respects, i.e. the surface density distribution, flaring and warping, so that the rotation curve (RC) is determined accurately, despite the fact that i~90^o. The RC rises slowly in the inner 6 kpc, is ~ constant at 100 km/s out to 10kpc, and decreases in Keplerian fashion by 15% at the last measured point at 14 kpc. The RC constrains the stellar mass-to-light ratio to lie between 50 and 100% of the ``maximum-disk'' value. A new technique is presented to determine simultaneously the inclination and the thickness of the gas layer from HI observations which uses the apparent widths at many channels and works for i >= 60^o. The inclination of the un- warped disk is about 84.5 degrees: beyond D25/2 the inclination changes gra- dually till 82.5^o, while at large radii the disk seems to warp back to the plane defined by the inner disk. The gaseous velocity dispersion is roughly constant within the optical disk (8.5+/-1 km/s) and increases slightly beyond. The thickness of the gas layer increases gradually from 0.4 to 1.5 kpc. In an accompanying paper (AJ, Aug. 1996; astro-ph/9605111). I use the measurements presented in this paper to infer that the dark matter halo of NGC 4244 is highly flattened.

In a previous paper (Olling 1995, \aj, 110, 591; astro-ph/9505002) a method was developed to determine the shapes of dark matter halos of spiral galaxies from an accurate determination of the rotation curve, the flaring of the gas layer and the velocity dispersion in the HI. Here this method is applied to the almost edge-on Scd galaxy NGC 4244 for which the necessary parameters are determined in the accompanying paper (AJ, Aug. 1996; astro-ph/9605110). The observed flaring of the HI beyond the optical disk puts significant constraints on the shape of the dark matter halo, which are almost independent of the stellar mass-to-light ratio. NGC 4244's dark matter halo is found to be highly flattened with a shortest-to-longest axis ratio of 0.2 (-0.1)(+0.3). If the dark matter is disk-like, the data presented in this paper imply that the vertical velocity dispersion of the dark matter must be 10% - 30% larger than the measured tangential dispersion in the HI. Alternatively, the measured flaring curve is consistent with a round halo if the gaseous velocity dispersion ellipsoid is anisotropic. In that case the vertical dispersion of the gas is 50 - 70% of the measured tangential velocity dispersion.

I present a new method of deriving the shape of the dark matter (DM) halos of spiral galaxies. The method relies on the comparison of model predictions with high spectral and spatial resolution HI observations of the gas layer. The potential arising from the {\em total} mass distribution of the galaxy is used in the calculation of the vertical distribution of the gas. I developed a new algorithm to calculate the force field of an arbitrary, azimuthally symmetric, density distribution. This algorithm is used to calculate the forces due to the radially truncated stellar disk as well as of the flaring gas layer. I use a simple two-parameter family of disk-halo models which have essentially the same observed equatorial rotation curve but different vertical forces. This mass model is composed of a stellar disk with constant M/L, and a DM-halo with a given axial ratio. I approximate the radial force due to the gaseous disk, and iteratively determine the vertical force due to the global distribution of the gas. The thickness of the gaseous disk is sensitive to both the flattening of the DM-halo and the self-gravity of the gas, but not to the particular choice of disk-halo decomposition. I show that the determination of the thickness of the gas layer is not restricted to edge-on galaxies, but can be measured for moderately inclined systems as well.