Strong lenses are extremely useful probes of the distribution of matter on galaxy and cluster scales at cosmological distances, but are rare and difficult to find. The number of currently known lenses is on the order of 1,000. We wish to use crowdsourcing to carry out a lens search targeting massive galaxies selected from over 442 square degrees of photometric data from the Hyper Suprime-Cam (HSC) survey. We selected a sample of $\sim300,000$ galaxies with photometric redshifts in the range $0.2 < z_{phot} < 1.2$ and photometrically inferred stellar masses $\log{M_*} > 11.2$. We crowdsourced lens finding on this sample of galaxies on the Zooniverse platform, as part of the Space Warps project. The sample was complemented by a large set of simulated lenses and visually selected non-lenses, for training purposes. Nearly 6,000 citizen volunteers participated in the experiment. In parallel, we used YattaLens, an automated lens finding algorithm, to look for lenses in the same sample of galaxies. Based on a statistical analysis of classification data from the volunteers, we selected a sample of the most promising $\sim1,500$ candidates which we then visually inspected: half of them turned out to be possible (grade C) lenses or better. Including lenses found by YattaLens or serendipitously noticed in the discussion section of the Space Warps website, we were able to find 14 definite lenses, 129 probable lenses and 581 possible lenses. YattaLens found half the number of lenses discovered via crowdsourcing. Crowdsourcing is able to produce samples of lens candidates with high completeness and purity, compared to currently available automated algorithms. A hybrid approach, in which the visual inspection of samples of lens candidates pre-selected by discovery algorithms and/or coupled to machine learning is crowdsourced, will be a viable option for lens finding in the 2020s.

The major semiaxes of the planets in our Solar System obey a simple geometric progression known as the Titius-Bode Relation (TBR), whose physical origin remains disputed. It has been shown that the exoplanetary systems follow a similar (but not identical) progression of the form a_n= a_0 e^(bn), where a_0, b are constants to be determined for each system. Since its formulation, the Titius-Bode Relation has proved to be highly predictive in our Solar System. Using data from 27 exoplanetary systems with 5 or more planets and applying a proposed method, we conclude that reliable TB-like fits can be obtained for systems with at least 4 planets and that the precision of the TBR is 78%. By means of a statistical test we show that the periods of planets in real exoplanetary systems are not consistent with a random distribution. Rather, they show signs that their configuration is shaped by their mutual interactions.

As an illustration of the value of historical measures we present some examples of the dynamical evolution of multiple systems resembling the Orion Trapezium. We constructed models by combining carefully selected historical measures of the separations among components of young massive stellar systems with modern observations. By computing large numbers of fictitious systems resembling real trapezia we were able to simulate the dynamical evolution of such systems. Our results on the dynamical fate of the Orion Trapezium and of ten additional young clusters resembling the Orion Trapezium show extremely short dynamical lifetimes for these systems.

Under Newtonian dynamics, the relative motion of the components of a binary star should follow a Keplerian scaling with separation. Once orientation effects and a distribution of ellipticities are accounted for, dynamical evolution can be modelled to include the effects of Galactic tides and stellar mass perturbers, over the lifetime of the solar neighbourhood. This furnishes a prediction for the relative velocity between the components of a binary and their projected separation. Taking a carefully selected small sample of 81 solar neighbourhood wide binaries from the {\it Hipparcos} catalogue, we identify these same stars in the recent Gaia DR2, to test the prediction mentioned using the latest and most accurate astrometry available. The results are consistent with the Newtonian prediction for projected separations below 7000 AU, but inconsistent with it at larger separations, where accelerations are expected to be lower than the critical $a_{0}=1.2 \times 10^{-10} $ { m s$^{-2}$} value of MONDian gravity. This result challenges Newtonian gravity at low accelerations and shows clearly the appearance of gravitational anomalies of the type usually attributed to dark matter at galactic scales, now at much smaller stellar scales.

We present new CCD $VI$ photometry of the distant globular cluster Pal 13. Fourier decomposition of the light curves of the three cluster member RRab stars lead to estimations of [Fe/H]=-1.65, and a distance of 23.67$\pm$0.57 kpc. Light and colour near minimum phases for RRab stars leads to an estimate of $E(B-V)$=0.104 $\pm$ 0.001. A $V/(V-I)$ colour-magnitude diagram, built exclusively with likely star members, shows consistency with the above parameters and an age of 12 Gyrs. A search of variable stars in the field of view of our images revealed the variability of a red giant cluster member and of three probably non-member stars; two RRab stars and one W Virginis star or CW. The GAIA proper motions of member stars in Pal 13 show a significant scatter, consistent with the scenario of the cluster being tidally stripped.

We have selected archival observational data for several O and B trapezia in the Milky Way. For each of the main components of the trapezia we obtained transverse velocities from the historical separation data. With this information, and with the stellar masses of the main components, we studied the dynamical evolution of ensembles of multiple systems mimicking each one of the trapezia. For this purpose we conducted numerical $N-$body integrations using the best available values for the masses, the observed positions and {\bf transverse} velocities, randomly generated radial velocities, and random line-of-sight ($z$) positions for all components. Random perturbations were assigned to the observed quantities, compatible with the observational errors. A large fraction of the simulated systems ({\bf between 70 and 90 percent}) turned out to be unbound. The properties of the evolving systems are studied at different values of the evolution time. We find that the dynamical lifetimes of both the bound and unbound seems to be quite short, of less than 10~000 yr for the unbound systems, and of 10~000 to 20~000 yr for the bound systems. The end result of the simulations is usually a binary {\bf with semiaxes of a few hundred AU} sometimes a triple of hierarchical on non-hierarchical type. Non-hierarchical triples formed during the integrations (which are dynamically unstable) were found to have much longer lifetimes, of 250~000 to 500~000 yr. The frequency distributions of the major semi-axes and eccentricities of the resulting binaries are discussed and compared with observational properties of binary systems from the literature.

K2-138 is a moderately bright (V = 12.2, K = 10.3) main sequence K-star observed in Campaign 12 of the NASA K2 mission. It hosts five small (1.6-3.3R_Earth) transiting planets in a compact architecture. The periods of the five planets are 2.35 d, 3.56 d, 5.40 d, 8.26 d, and 12.76 d, forming an unbroken chain of near 3:2 resonances. Although we do not detect the predicted 2-5 minute transit timing variations with the K2 timing precision, they may be observable by higher cadence observations with, for example, Spitzer or CHEOPS. The planets are amenable to mass measurement by precision radial velocity measurements, and therefore K2-138 could represent a new benchmark systems for comparing radial velocity and TTV masses. K2-138 is the first exoplanet discovery by citizen scientists participating in the Exoplanet Explorers project on the Zooniverse platform.

We present the results of a systematic mapping of seasonally sculpted terrains on the South Polar region of Mars with the Planet Four: Terrains (P4T) online citizen science project. P4T enlists members of the general public to visually identify features in the publicly released Mars Reconnaissance Orbiter CTX images. In particular, P4T volunteers are asked to identify: 1) araneiforms (including features with a central pit and radiating channels known as 'spiders'); 2) erosional depressions, troughs, mesas, ridges, and quasi-circular pits characteristic of the South Polar Residual Cap (SPRC) which we collectively refer to as 'Swiss cheese terrain', and 3) craters. In this work we present the distributions of our high confidence classic spider araneiforms and Swiss cheese terrain identifications. We find no locations within our high confidence spider sample that also have confident Swiss cheese terrain identifications. Previously spiders were reported as being confined to the South Polar Layered Deposits (SPLD). Our work has provided the first identification of spiders at locations outside of the SPLD, confirmed with high resolution HiRISE imaging. We find araneiforms on the Amazonian and Hesperian polar units and the Early Noachian highland units, with 75% of the identified araneiform locations in our high confidence sample residing on the SPLD. With our current coverage, we cannot confirm whether these are the only geologic units conducive to araneiform formation on the Martian South Polar region. Our results are consistent with the current CO2 jet formation scenario with the process exploiting weaknesses in the surface below the seasonal CO2 ice sheet to carve araneiform channels into the regolith over many seasons. These new regions serve as additional probes of the conditions required for channel creation in the CO2 jet process. (Abridged)

We present very-high-energy $\gamma$-ray observations of the BL Lac object 1ES 2344+514 taken by the Very Energetic Radiation Imaging Telescope Array System (VERITAS) between 2007 and 2015. 1ES 2344+514 is detected with a statistical significance above background of $20.8\sigma$ in $47.2$ hours (livetime) of observations, making this the most comprehensive very-high-energy study of 1ES 2344+514 to date. Using these observations the temporal properties of 1ES 2344+514 are studied on short and long times scales. We fit a constant flux model to nightly- and seasonally-binned light curves and apply a fractional variability test, to determine the stability of the source on different timescales. We reject the constant-flux model for the 2007-2008 and 2014-2015 nightly-binned light curves and for the long-term seasonally-binned light curve at the $> 3\sigma$ level. The spectra of the time-averaged emission before and after correction for attenuation by the extragalactic background light are obtained. The observed time-averaged spectrum above 200 GeV is satisfactorily fitted (${\chi^2/NDF = 7.89/6}$) by a power-law function with index $\Gamma = 2.46 \pm 0.06_{stat} \pm 0.20_{sys} $ and extends to at least 8 TeV. The extragalactic-background-light-deabsorbed spectrum is adequately fit (${\chi^2/NDF = 6.73/6}$) by a power-law function with index $\Gamma = 2.15 \pm 0.06_{stat} \pm 0.20_{sys} $ while an F-test indicates that the power-law with exponential cutoff function provides a marginally-better fit ($\chi^2/NDF $ = $2.56 / 5 $) at the 2.1$\sigma$ level. The source location is found to be consistent with the published radio location and its spatial extent is consistent with a point source.

Large modern surveys require efficient review of data in order to find transient sources such as supernovae, and to distinguish such sources from artefacts and noise. Much effort has been put into the development of automatic algorithms, but surveys still rely on human review of targets. This paper presents an integrated system for the identification of supernovae in data from Pan-STARRS1, combining classifications from volunteers participating in a citizen science project with those from a convolutional neural network. The unique aspect of this work is the deployment, in combination, of both human and machine classifications for near real-time discovery in an astronomical project. We show that the combination of the two methods outperforms either one used individually. This result has important implications for the future development of transient searches, especially in the era of LSST and other large-throughput surveys.

We describe an experiment, located in south-east Colorado, USA, that measured aerosol optical depth profiles using two Lidar techniques. Two independent detectors measured scattered light from a vertical UV laser beam. One detector, located at the laser site, measured light via the inelastic Raman backscattering process. This is a common method used in atmospheric science for measuring aerosol optical depth profiles. The other detector, located approximately 40km distant, viewed the laser beam from the side. This detector featured a 3.5m2 mirror and measured elastically scattered light in a bistatic Lidar configuration following the method used at the Pierre Auger cosmic ray observatory. The goal of this experiment was to assess and improve methods to measure atmospheric clarity, specifically aerosol optical depth profiles, for cosmic ray UV fluorescence detectors that use the atmosphere as a giant calorimeter. The experiment collected data from September 2010 to July 2011 under varying conditions of aerosol loading. We describe the instruments and techniques and compare the aerosol optical depth profiles measured by the Raman and bistatic Lidar detectors.

We discuss recent observational data on the transverse and radial velocities, as well as on the masses of the main components of the Orion Trapezium. Based on the most reliable values of these quantities we study the dynamical evolution of ensembles of multiple systems mimicking the Orion Trapezium. To this end we conduct numerical $N-$body integrations using the observed masses, planar positions and velocities, radial velocities, and random line-of-sight ($z$) positions for all components. We include perturbations in these quantities compatible with the observational errors. We find the dynamical lifetimes of such systems to be quite short, of the order of 10 to 50 thousand years. The end result of the simulations is usually a tight binary, or sometimes a hierarchical triple. The properties of the evolved systems are studied at different values of the crossing times. The frequency distributions of the major semiaxes and eccentricities of the resulting binaries are discussed and compared with observations.

Employing a 40-kW radio-frequency transmitter just west of Delta, UT, and operating at 54.1 MHz, the TARA (Telescope Array RAdar) experiment seeks radar detection of extensive air showers (EAS) initiated by ultra-high energy cosmic rays (UHECR). For UHECR with energies in excess of $10^{19}$ eV, the Doppler-shifted "chirps" resulting from EAS shower core radar reflections should be observable above background (dominantly galactic) at distances of tens of km from the TARA transmitter. In order to stereoscopically reconstruct cosmic ray chirps, two remote, autonomous self-powered receiver stations have been deployed. Each remote station (RS) combines both low power consumption as well as low cost. Triggering logic, the powering and communication systems, and some specific details of hardware components are discussed.

Construction was completed during summer 2013 on the Telescope Array RAdar (TARA) bi-static radar observatory for Ultra-High Energy Cosmic Rays (UHECR). TARA is co-located with the Telescope Array, the largest "conventional" cosmic ray detector in the Northern Hemisphere, in radio-quiet Western Utah. TARA employs an 8 MW Effective Radiated Power (ERP) VHF transmitter and smart receiver system based on a 250 MS/s data acquisition system in an effort to detect the scatter of sounding radiation by UHECR-induced atmospheric ionization. TARA seeks to demonstrate bi-static radar as a useful new remote sensing technique for UHECRs, extending their detection aperture far beyond what is accessible by conventional means. In this report, we describe the design and performance of the TARA transmitter and receiver systems.

We present an improved catalog of halo wide binaries, compiled from an extensive literature search. Most of our binaries stem from the common proper motion binary catalogs by Allen et al. (2004), and Chanam\'e \& Gould. (2004) but we have also included binaries from the lists of Ryan (1992) and Zapatero-Osorio \& Martin (2004). All binaries were carefully checked and their distances and systemic radial velocities are included, when available. Probable membership to the halo population was tested by means of reduced proper motion diagrams for 251 candidate halo binaries. After eliminating obvious disk binaries we ended up with 211 probable halo binaries, for 150 of which radial velocities are available. We compute galactic orbits for these 150 binaries and calculate the time they spend within the galactic disk. Considering the full sample of 251 candidate halo binaries as well as several subsamples, we find that the distribution of angular separations (or expected major semiaxes) follows a power law $f(a) \sim a^{-1}$ (Oepik's relation) up to different limits. For the 50 most disk-like binaries, those that spend their entire lives within $z = \pm 500$~pc, this limit is found to be 19,000 AU (0.09 pc), while for the 50 most halo-like binaries, those that spend on average only 18\% of their lives within $z = \pm 500$~pc, the limit is 63,000 AU (0.31 pc). In a companion paper we employ this catalog to establish limits on the masses of the halo massive perturbers (MACHOs).

In order to determine an upper bound for the mass of the massive compact halo objets (MACHOs) we use the halo binaries contained in a recent catalog (Allen \& Monroy-Rodr\'{\i}guez 2013). To dynamically model their interactions with massive perturbers a Monte Carlo simulation is conducted, using an impulsive approximation method and assuming a galactic halo constituted by massive particles of a characteristic mass. The results of such simulations are compared with several subsamples of our improved catalog of candidate halo wide binaries. In accordance with Quinn et al. (2009) we also find our results to be very sensitive to the widest binaries. However, our larger sample, together with the fact that we can obtain galactic orbits for 150 of our systems, allows a more reliable estimate of the maximum MACHO mass than that obtained previously. If we employ the entire sample of 211 candidate halo stars we obtain an upper limit of $112 M_\sun$. However, using the 150 binaries in our catalog with computed galactic orbits we are able to refine our fitting criteria. Thus, for the 100 most halo-like binaries we obtain a maximum MACHO mass of $21-68 M_\sun$. Furthermore, we can estimate the dynamical effects of the galactic disk using binary samples that spend progressively shorter times within the disk. By extrapolating the limits obtained for our most reliable -albeit smallest- sample we find that as the time spent within the disk tends to zero the upper bound of the MACHO mass tends to less than $5 M_\sun$. The non-uniform density of the halo has also been taken into account, but the limit obtained, less than $5 M_\sun$, does not differ much from the previous one. Together with microlensing studies that provide lower limits on the MACHO mass, our results essentially exclude the existence of such objects in the galactic halo.

Recent observations for three types of astrophysical systems severely challenge the GR plus dark matter scenario, showing a phenomenology which is what modified gravity theories predict. Stellar kinematics in the outskirts of globular clusters show the appearance of MOND type dynamics on crossing the $a_{0}$ threshold. Analysis shows a ``Tully-Fisher'' relation in these systems, a scaling of dispersion velocities with the fourth root of their masses. Secondly, an anomaly has been found at the unexpected scales of wide binaries in the solar neighbourhood. Binary orbital velocities cease to fall along Keplerian expectations, and settle at a constant value, exactly on crossing the $a_{0}$ threshold. Finally, the inferred infall velocity of the bullet cluster is inconsistent with the standard cosmological scenario, where much smaller limit encounter velocities appear. This stems from the escape velocity limit present in standard gravity; the ``bullet'' should not hit the ``target'' at more than the escape velocity of the joint system, as it very clearly did. These results are consistent with extended gravity, but would require rather contrived explanations under GR, each. Thus, observations now put us in a situation where modifications to gravity at low acceleration scales cease to be a matter of choice, to now become inevitable.

Using the novel Diffracto-Astrometry technique, we analyze 44 Hubble Space Telescope Wide Field Planetary Camera 2 images of the Orion Trapezium (OT) taken over a span of 12 yr (1995-2007). We measure the relative positions of the six brighter OT components (A-F) and supplement these results with measurements of the relative separations and position angles taken from the literature, thus extending our analysis time base to ~200 yr. For every pair of components we find the relative rate of separation as well as the temporal rate of change of their position angles, which enable us to determine the relative kinematics of the system. Component E shows a velocity larger than the OT's escape velocity, thus confirming that it is escaping from the gravitational pull of this system.

That the stellar halo of the Milky Way has a density profile which to first approximation satisfies $\rho \propto r^{-3}$ has been known for a long time. More recently, it has become clear that M31 also has such an extended stellar halo, which approximately follows the same radial scaling. Studies of distant galaxies have revealed the same phenomenology. Also, we now know that the density profiles of the globular cluster systems of our Galaxy and Andromeda to first approximation follow $\rho \propto r^{-3}$, $\Sigma \propto R^{-2}$ in projection. Recently, diffuse populations of stars have been detected spherically surrounding a number of Galactic globular clusters, extending much beyond the Newtonian tidal radii, often without showing any evidence of tidal features. Within the standard Newtonian and GR scenario, numerous and diverse particular explanations have been suggested, individually tailored to each of the different classes of systems described above. Here we show that in a MONDian gravity scenario, any isothermal tenuous halo of tracer particles forming a small perturbation surrounding a spherically symmetric mass distribution, will have an equilibrium configuration which to first approximation satisfies a $\rho \propto r^{-3}$ scaling.

Over the past couple of years, a number of observational studies have confirmed the flattening of the radial velocity dispersion profiles for stars in various nearby globular clusters. As the projected radial coordinate is increased, a radius appears beyond which, the measured velocity dispersion ceases to drop and settles at a fixed value, $\sigma_{\infty}$. Under Newtonian gravity, this is explained by invoking tidal heating from the overall Milky Way potential on the outer, more loosely bound stars, of the globular clusters in question. From the point of view of modified gravity theories, such an outer flattening is expected on crossing the critical acceleration threshold $a_{0}$, beyond which, a transition to MONDian dynamics is expected, were equilibrium velocities cease to be a function of distance. In this paper we attempt to sort out between the above competing explanations, by looking at their plausibility in terms of an strictly empirical approach. We determine Newtonian tidal radii using masses accurately calculated through stellar population modelling, and hence independent of any dynamical assumptions, distances, size and orbital determinations for a sample of 16 globular clusters. We show that their Newtonian tidal radii at perigalacticon are generally larger that the radii at which the flattening in the velocity dispersion profiles occurs, by large factors of 4, on average. While this point makes the Newtonian tidal explanation suspect, it is found that the radii at which the flattening is observed on average correlate with the radii where the $a_{0}$ threshold is crossed, and that $\sigma_{\infty}$ values scale with the fourth root of the total masses, all features predicted under modified gravity theories.

Assuming Newton's gravity and GR to be valid at all scales leads to the dark matter hypothesis as a requirement demanded by the observed dynamics and measured baryonic content at galactic and extragalactic scales. Alternatively, modified gravity scenarios where a change of regime appears at acceleration scales $a<a_{0}$ have been proposed. This modified regime at $a<a_{0}$ will generically be characterised by equilibrium velocities which become independent of distance. Here we identify a critical test in this debate and we propose its application to samples of wide binary stars. Since for $1 M_{\odot}$ systems the acceleration drops below $a_{0}$ at scales of around 7000 AU, a statistical survey of wide binaries with relative velocities and separations reaching $10^{4}$ AU and beyond should prove useful to the above debate. We apply the proposed test to the best currently available data. Results show a constant upper limit to the relative velocities in wide binaries which is independent of separation for over three orders of magnitude, in analogy with galactic flat rotation curves in the same $a<a_{0}$ acceleration regime. Our results are suggestive of a breakdown of Kepler's third law beyond $a \approx a_{0}$ scales, in accordance with generic predictions of modified gravity theories designed not to require any dark matter at galactic scales and beyond.

Using dimensional analysis techniques we present an extension of Newton's gravitational theory built under the assumption that Milgrom's acceleration constant is a fundamental quantity of nature. The gravitational force converges to Newton's gravity and to a MOND-like description in two different mass and length regimes. It is shown that a modification on the force sector (and not in the dynamical one as MOND does) is more convenient and can reproduce and predict different phenomena usually ascribed to dark matter at the non-relativistic level.

Modified gravity scenarios where a change of regime appears at acceleration scales $a<a_{0}$ have been proposed. Since for $1 M_{\odot}$ systems the acceleration drops below $a_{0}$ at scales of around 7000 AU, a statistical survey of wide binaries with relative velocities and separations reaching $10^{4}$ AU and beyond should prove useful to the above debate. We apply the proposed test to the best currently available data. Results show a constant upper limit to the relative velocities in wide binaries which is independent of separation for over three orders of magnitude, in analogy with galactic flat rotation curves in the same $a<a_{0}$ acceleration regime. Our results are suggestive of a breakdown of Kepler's third law beyond $a \approx a_{0}$ scales, in accordance with generic predictions of modified gravity theories designed not to require any dark matter at galactic scales and beyond.

We present a new method for confirming transiting planets based on the combination of transit timingn variations (TTVs) and dynamical stability. Correlated TTVs provide evidence that the pair of bodies are in the same physical system. Orbital stability provides upper limits for the masses of the transiting companions that are in the planetary regime. This paper describes a non-parametric technique for quantifying the statistical significance of TTVs based on the correlation of two TTV data sets. We apply this method to an analysis of the transit timing variations of two stars with multiple transiting planet candidates identified by Kepler. We confirm four transiting planets in two multiple planet systems based on their TTVs and the constraints imposed by dynamical stability. An additional three candidates in these same systems are not confirmed as planets, but are likely to be validated as real planets once further observations and analyses are possible. If all were confirmed, these systems would be near 4:6:9 and 2:4:6:9 period commensurabilities. Our results demonstrate that TTVs provide a powerful tool for confirming transiting planets, including low-mass planets and planets around faint stars for which Doppler follow-up is not practical with existing facilities. Continued Kepler observations will dramatically improve the constraints on the planet masses and orbits and provide sensitivity for detecting additional non-transiting planets. If Kepler observations were extended to eight years, then a similar analysis could likely confirm systems with multiple closely spaced, small transiting planets in or near the habitable zone of solar-type stars.

Using the most recent proper-motion determination of the old, Solar-metallicity, Galactic open cluster M 67, in orbital computations in a non-axisymmetric model of the Milky Way, including a bar and 3D spiral arms, we explore the possibility that the Sun once belonged to this cluster. We have performed Monte Carlo numerical simulations to generate the present-day orbital conditions of the Sun and M 67, and all the parameters in the Galactic model. We compute 3.5 \times 10^5 pairs of orbits Sun-M 67 looking for close encounters in the past with a minimum distance approach within the tidal radius of M 67. In these encounters we find that the relative velocity between the Sun and M 67 is larger than 20 km/s. If the Sun had been ejected from M 67 with this high velocity by means of a three-body encounter, this interaction would destroy an initial circumstellar disk around the Sun, or disperse its already formed planets. We also find a very low probability, much less than 10^-7, that the Sun was ejected from M 67 by an encounter of this cluster with a giant molecular cloud. This study also excludes the possibility that the Sun and M 67 were born in the same molecular cloud. Our dynamical results convincingly demonstrate that M67 could not have been the birth cluster of our Solar System.

We present a theoretical analysis of the measuring algorithm we use when applying the Diffracto-Astrometry technique to Hubble Space Telescope Wide Field Planetary Camera 2 (WFPC2) saturated stellar images. Theoretical Point Spread Functions (PSFs) were generated using the Tiny Tim software for the four CCDs in the WFPC2 and for some of the available filters. These images were then measured with our Diffracto-Astrometry measuring algorithm using only their diffraction pattern, and positions for the simulated PSFs on each generated CCD-frame were obtained. The measuring algorithm recovers the original positions reasonably well ($\pm 0.1\ \rm{to} \ \pm 0.4$ pixels). However, slight deviations from the original values are observed. These also vary with position over the entire surface of the CCD. We adjust the difference between the real and the measured position with a quadratic function of the coordinates. The transformation coefficients also present a slight correlation with the filter effective wavelength. Application of these transformation coefficients allows us to determine the position of a stellar image with a precision of a few hundredths of a pixel.

Delta Scuti (delta Sct) stars are opacity-driven pulsators with masses of 1.5-2.5M$_{\odot}$, their pulsations resulting from the varying ionization of helium. In less massive stars such as the Sun, convection transports mass and energy through the outer 30 per cent of the star and excites a rich spectrum of resonant acoustic modes. Based on the solar example, with no firm theoretical basis, models predict that the convective envelope in delta Sct stars extends only about 1 per cent of the radius, but with sufficient energy to excite solar-like oscillations. This was not observed before the Kepler mission, so the presence of a convective envelope in the models has been questioned. Here we report the detection of solar-like oscillations in the delta Sct star HD 187547, implying that surface convection operates efficiently in stars about twice as massive as the Sun, as the ad hoc models predicted.

We calculate precise stellar radii and surface gravities from the asteroseismic analysis of over 500 solar-type pulsating stars observed by the Kepler space telescope. These physical stellar properties are compared with those given in the Kepler Input Catalog (KIC), determined from ground-based multi-color photometry. For the stars in our sample, we find general agreement but we detect an average overestimation bias of 0.23 dex in the KIC determination of log (g) for stars with log (g)_KIC > 4.0 dex, and a resultant underestimation bias of up to 50% in the KIC radii estimates for stars with R_KIC < 2 R sun. Part of the difference may arise from selection bias in the asteroseismic sample; nevertheless, this result implies there may be fewer stars characterized in the KIC with R ~ 1 R sun than is suggested by the physical properties in the KIC. Furthermore, if the radius estimates are taken from the KIC for these affected stars and then used to calculate the size of transiting planets, a similar underestimation bias may be applied to the planetary radii.

Kepler observations of the star HD 187091 (KID 8112039, hereafter KOI-54) revealed a remarkable light curve exhibiting sharp periodic brightening events every 41.8 days with a superimposed set of oscillations forming a beating pattern in phase with the brightenings. Spectroscopic observations revealed that this is a binary star with a highly eccentric orbit, e=0.83. We are able to match the Kepler light curve and radial velocities with a nearly face-on (i=5.5 degree) binary star model in which the brightening events are caused by tidal distortion and irradiation of nearly identical A stars during their close periastron passage. The two dominant oscillations in the light curve, responsible for the beating pattern, have frequencies that are the 91st and 90th harmonic of the orbital frequency. The power spectrum of the light curve, after removing the binary star brightening component, reveals a large number of pulsations, 30 of which have a signal-to-noise ratio > 7. Nearly all of these pulsations have frequencies that are either integer multiples of the orbital frequency or are tidally-split multiples of the orbital frequency. This pattern of frequencies unambiguously establishes the pulsations as resonances between the dynamic tides at periastron and the free oscillation modes of one of the stars. KOI-54 is only the 4th star to show such a phenomenon, and is by far the richest in terms of excited modes.

The Kepler Mission has recently announced the discovery of Kepler-10 b, the smallest exoplanet discovered to date and the first rocky planet found by the spacecraft. A second, 45-day period transit-like signal present in the photometry from the first eight months of data could not be confirmed as being caused by a planet at the time of that announcement. Here we apply the light-curve modeling technique known as BLENDER to explore the possibility that the signal might be due to an astrophysical false positive (blend). To aid in this analysis we report the observation of two transits with the Spitzer Space Telescope at 4.5 {\mu}m. When combined they yield a transit depth of 344 \pm 85 ppm that is consistent with the depth in the Kepler passband (376 \pm 9 ppm, ignoring limb darkening), which rules out blends with an eclipsing binary of a significantly different color than the target. Using these observations along with other constraints from high resolution imaging and spectroscopy we are able to exclude the vast majority of possible false positives. We assess the likelihood of the remaining blends, and arrive conservatively at a false alarm rate of 1.6 \times 10-5 that is small enough to validate the candidate as a planet (designated Kepler-10 c) with a very high level of confidence. The radius of this object is measured to be Rp = 2.227+0.052 -0.057 Earth radii. Kepler-10 c represents another example (with Kepler-9 d and Kepler-11 g) of statistical "validation" of a transiting exoplanet, as opposed to the usual "confirmation" that can take place when the Doppler signal is detected or transit timing variations are measured. It is anticipated that many of Kepler's smaller candidates will receive a similar treatment since dynamical confirmation may be difficult or impractical with the sensitivity of current instrumentation.

We investigate the possibility of nearly-equally spaced periods in 13 hot subdwarf B (sdB) stars observed with the Kepler spacecraft and one observed with CoRoT. Asymptotic limits for gravity (g-)mode pulsations provide relationships between equal period spacings of modes with differing degrees and relationships between periods of the same radial order but differing degrees. Period transforms, Kolmogorov-Smirnov tests, and linear least-squares fits have been used to detect and determine the significance of equal period spacings. We have also used Monte Carlo simulations to estimate the likelihood that the detected spacings could be produced randomly. Period transforms for nine of the Kepler stars indicate ell=1 period spacings, with five also showing peaks for ell=2 modes. 12 stars indicate ell=1 modes using the Kolmogorov-Smirnov test while another shows solely ell=2 modes. Monte Carlo results indicate that equal period spacings are significant in 10 stars above 99% confidence and 13 of the 14 are above 94% confidence. For 12 stars, the various methods find consistent regular period spacing values to within the errors, two others show some inconsistencies, likely caused by binarity, and the last has significant detections but the mode assignment disagrees between methods. We find a common ell=1 period spacing spanning a range from 231 to 272 s allowing us to correlate pulsation modes with 222 periodicities and that the ell=2 period spacings are related to the ell=1 spacings by the asymptotic relationship $1/\sqrt{3}$. We briefly discuss the impact of equal period spacings which indicate low-degree modes with a lack of significant mode trappings.

On 1 February 2011 the Kepler Mission released data for 156,453 stars observed from the beginning of the science observations on 2 May through 16 September 2009. There are 1235 planetary candidates with transit like signatures detected in this period. These are associated with 997 host stars. Distributions of the characteristics of the planetary candidates are separated into five class-sizes; 68 candidates of approximately Earth-size (radius < 1.25 Earth radii), 288 super-Earth size (1.25 Earth radii < radius < 2 Earth radii), 662 Neptune-size (2 Earth radii < radius < 6 Earth radii), 165 Jupiter-size (6 Earth radii < radius < 15 Earth radii), and 19 up to twice the size of Jupiter (15 Earth radii < radius < 22 Earth radii). In the temperature range appropriate for the habitable zone, 54 candidates are found with sizes ranging from Earth-size to larger than that of Jupiter. Five are less than twice the size of the Earth. Over 74% of the planetary candidates are smaller than Neptune. The observed number versus size distribution of planetary candidates increases to a peak at two to three times Earth-size and then declines inversely proportional to area of the candidate. Our current best estimates of the intrinsic frequencies of planetary candidates, after correcting for geometric and sensitivity biases, are 6% for Earth-size candidates, 7% for super-Earth size candidates, 17% for Neptune-size candidates, and 4% for Jupiter-size candidates. Multi-candidate, transiting systems are frequent; 17% of the host stars have multi-candidate systems, and 33.9% of all the candidates are part of multi-candidate systems.

We present results from the final six months of a survey to search for pulsations in white dwarfs and hot subdwarf stars with the Kepler spacecraft. Spectroscopic observations are used to separate the objects into accurate classes, and we explore the physical parameters of the subdwarf B (sdB) stars and white dwarfs in the sample. From the Kepler photometry and our spectroscopic data, we find that the sample contains 5 new pulsators of the V1093 Her type, one AM CVn type cataclysmic variable, and a number of other binary systems. This completes the survey for compact pulsators with Kepler. No V361 Hya type of short-period pulsating sdB stars were found in this half, leaving us with a total of one single multiperiodic V361 Hya and 13 V1093 Her pulsators for the full survey. Except for the sdB pulsators, no other clearly pulsating hot subdwarfs or white dwarfs were found, although a few low-amplitude candidates still remain. The most interesting targets discovered in this survey will be observed throughout the remainder of the Kepler Mission, providing the most long-term photometric datasets ever made on such compact, evolved stars. Asteroseismic investigations of these datasets will be invaluable in revealing the interior structure of these stars, and will boost our understanding of their evolutionary history.

We report results of initial work done on selected candidate Cepheids to be observed with the Kepler space telescope. Prior to the launch 40 candidates were selected from previous surveys and databases. The analysis of the first 322 days of Kepler photometry, and recent ground-based follow-up multicolour photometry and spectroscopy allowed us to confirm that one of these stars, V1154 Cyg (KIC 7548061), is indeed a 4.9-d Cepheid. Using the phase lag method we show that this star pulsates in the fundamental mode. New radial velocity data are consistent with previous measurements, suggesting that a long-period binary component is unlikely. No evidence is seen in the ultra-precise, nearly uninterrupted Kepler photometry for nonradial or stochastically excited modes at the micromagnitude level. The other candidates are not Cepheids but an interesting mix of possible spotted stars, eclipsing systems and flare stars.

The large number of stars for which uninterrupted high-precision photometric timeseries data are being collected with \textit{Kepler} and CoRoT initiated the development of automated methods to analyse the stochastically excited oscillations in main-sequence, subgiant and red-giant stars. Aims: We investigate the differences in results for global oscillation parameters of G and K red-giant stars due to different methods and definitions. We also investigate uncertainties originating from the stochastic nature of the oscillations. Methods: For this investigation we use Kepler data obtained during the first four months of operation. These data have been analysed by different groups using already published methods and the results are compared. We also performed simulations to investigate the uncertainty on the resulting parameters due to different realizations of the stochastic signal. Results: We obtain results for the frequency of maximum oscillation power (nu_max) and the mean large separation (<delta nu>) from different methods for over one thousand red-giant stars. The results for these parameters agree within a few percent and seem therefore robust to the different analysis methods and definitions used here. The uncertainties for nu_max and <delta nu> due to differences in realization noise are not negligible and should be taken into account when using these results for stellar modelling.

We derived the absolute proper motion (PM) of the old, solar-metallicity Galactic open cluster M67 using observations collected with CFHT (1997) and with LBT (2007). About 50 galaxies with relatively sharp nuclei allow us to determine the absolute PM of the cluster. We find (mu_alpha cos(delta),mu_delta)_J2000.0 = (-9.6+/-1.1,-3.7+/-0.8) mas/yr. By adopting a line-of-sight velocity of 33.8+/-0.2 km/s, and assuming a distance of 815+/-50 pc, we explore the influence of the Galactic potential, with and without the bar and/or spiral arms, on the galactic orbit of the cluster.

The Kepler Mission, launched on Mar 6, 2009 was designed with the explicit capability to detect Earth-size planets in the habitable zone of solar-like stars using the transit photometry method. Results from just forty-three days of data along with ground-based follow-up observations have identified five new transiting planets with measurements of their masses, radii, and orbital periods. Many aspects of stellar astrophysics also benefit from the unique, precise, extended and nearly continuous data set for a large number and variety of stars. Early results for classical variables and eclipsing stars show great promise. To fully understand the methodology, processes and eventually the results from the mission, we present the underlying rationale that ultimately led to the flight and ground system designs used to achieve the exquisite photometric performance. As an example of the initial photometric results, we present variability measurements that can be used to distinguish dwarf stars from red giants.

We report the discovery and the Rossiter-McLaughlin effect of Kepler-8b, a transiting planet identified by the NASA Kepler Mission. Kepler photometry and Keck-HIRES radial velocities yield the radius and mass of the planet around this F8IV subgiant host star. The planet has a radius RP = 1.419 RJ and a mass, MP = 0.60 MJ, yielding a density of 0.26 g cm^-3, among the lowest density planets known. The orbital period is P = 3.523 days and orbital semima jor axis is 0.0483+0.0006/-0.0012 AU. The star has a large rotational v sin i of 10.5 +/- 0.7 km s^-1 and is relatively faint (V = 13.89 mag), both properties deleterious to precise Doppler measurements. The velocities are indeed noisy, with scatter of 30 m s^-1, but exhibit a period and phase consistent with the planet implied by the photometry. We securely detect the Rossiter-McLaughlin effect, confirming the planet's existence and establishing its orbit as prograde. We measure an inclination between the projected planetary orbital axis and the projected stellar rotation axis of lambda = -26.9 +/- 4.6 deg, indicating a moderate inclination of the planetary orbit. Rossiter-McLaughlin measurements of a large sample of transiting planets from Kepler will provide a statistically robust measure of the true distribution of spin-orbit orientations for hot jupiters in general.

The Kepler Mission Science Operations Center (SOC) performs several critical functions including managing the ~156,000 target stars, associated target tables, science data compression tables and parameters, as well as processing the raw photometric data downlinked from the spacecraft each month. The raw data are first calibrated at the pixel level to correct for bias, smear induced by a shutterless readout, and other detector and electronic effects. A background sky flux is estimated from ~4500 pixels on each of the 84 CCD readout channels, and simple aperture photometry is performed on an optimal aperture for each star. Ancillary engineering data and diagnostic information extracted from the science data are used to remove systematic errors in the flux time series that are correlated with these data prior to searching for signatures of transiting planets with a wavelet-based, adaptive matched filter. Stars with signatures exceeding 7.1 sigma are subjected to a suite of statistical tests including an examination of each star's centroid motion to reject false positives caused by background eclipsing binaries. Physical parameters for each planetary candidate are fitted to the transit signature, and signatures of additional transiting planets are sought in the residual light curve. The pipeline is operational, finding planetary signatures and providing robust eliminations of false positives.

The Kepler Mission seeks to detect Earth-size planets transiting solar-like stars in its ~115 deg^2 field of view over the course of its 3.5 year primary mission by monitoring the brightness of each of ~156,000 Long Cadence stellar targets with a time resolution of 29.4 minutes. We discuss the photometric precision achieved on timescales relevant to transit detection for data obtained in the 33.5-day long Quarter 1 (Q1) observations that ended 2009 June 15. The lower envelope of the photometric precision obtained at various timescales is consistent with expected random noise sources, indicating that Kepler has the capability to fulfill its mission. The Kepler light curves exhibit high precision over a large dynamic range, which will surely permit their use for a large variety of investigations in addition to finding and characterizing planets. We discuss the temporal characteristics of both the raw flux time series and the systematic error-corrected flux time series produced by the Kepler Science Pipeline, and give examples illustrating Kepler's large dynamic range and the variety of light curves obtained from the Q1 observations.

We have made a search for common proper motion (CPM) companions to the wide binaries in the solar vicinity. We found that the binary GJ 282AB has a very distant CPM companion (NLTT 18149) at a separation $s=1.09 \arcdeg$. Improved spectral types and radial velocities are obtained, and ages determined for the three components. The Hipparcos trigonometric parallaxes and the new radial velocities and ages turn out to be very similar for the three stars, and provide strong evidence that they form a physical system. At a projected separation of 55733AU from GJ 282AB, NLTT 18149 ranks among the widest physical companions known.

We present absolute astrometry of four radio sources in the Becklin-Neugebauer/Kleinman-Low (BN/KL) region, derived from archival data (taken in 1991, 1995, and 2000) as well as from new observations (taken in 2006). All data consist of 3.6 cm continuum emission and were taken with the Very Large Array in its highest angular resolution A configuration. We confirm the large proper motions of the BN object, the radio source I (GMR I) and the radio counterpart of the infrared source n (Orion-n), with values from 15 to 26 km/s. The three sources are receding from a point between them from where they seem to have been ejected about 500 years ago, probably via the disintegration of a multiple stellar system. We present simulations of very compact stellar groups that provide a plausible dynamical scenario for the observations. The radio source Orion-n appeared as a double in the first three epochs, but as single in 2006. We discuss this morphological change. The fourth source in the region, GMR D, shows no statistically significant proper motions. We also present new, accurate relative astrometry between BN and radio source I that restrict possible dynamical scenarios for the region. During the 2006 observations, the radio source GMR A, located about 1' to the NW of the BN/KL region, exhibited an increase in its flux density of a factor of ~3.5 over a timescale of one hour. This rapid variability at cm wavelengths is similar to that previously found during a flare at millimeter wavelengths that took place in 2003.

Absolute proper motions for six new globular clusters have recently been determined. This motivated us to obtain the Galactic orbits of these six clusters both in an axisymmetric Galactic potential and in a barred potential, such as the one of our Galaxy. Orbits are also obtained for a Galactic potential that includes spiral arms. The orbital characteristics are compared and discussed for these three cases. Tidal radii and destruction rates are also computed and discussed.

The frequency distribution f(a) of semi-major axis of double and multiple systems, as well as their eccentricities and mass ratios, contain valuable fossil information about the process of star formation and the dynamical history of the systems. In order to advance in the understanding of these questions, we have made an extensive analysis of the frequency distribution f (a) for wide binaries (a>25 AU) in the various published catalogues, as well as in our own (Poveda et al., 1994; Allen et al., 2000; Poveda & Hernandez, 2003). Based upon all these studies we have established that the frequency f(a) is function of the age of the system and follows Oepik's distribution f(a) ~ 1/a in the range of 100 AU < a < a[c](t), where a[c](t) is a critical semi-major axis beyond which binaries have dissociated by encounters with massive objects. We argue that the physics behind the distribution f(a) ~ 1/a is a process of energy relaxation, analogous to that present in stellar clusters (secular relaxation) or in spherical galaxies (violent relaxation). The frequency distribution of mass ratios in triple systems as well as the existence of runaway stars, indicate that both types of relaxation are important in the process of binary and multiple star formation.

Wide or fragile pairs are sensitive probes of the galactic potential, and they have been used to provide information about the galactic tidal field, the density of GMC and the masses of dark matter perturbers present in both the disk and the halo. Halo wide binaries and moving clusters, since they are likely to be the remains of past mergers or of dissolved clusters, can provide information on the dynamical and merger history of our Galaxy. Such remnants should continue to show similar motions over times of the order of their ages. We have looked for phase space groupings among the low-metallicity stars of Schuster et al. (2006) and have identified a number of candidate moving clusters. In several of the moving clusters we found a wide CPM binary already identified in our catalogue of wide binaries among high-velocity and metal-poor stars (Allen et al 2000a). Spectroscopic follow-up studies of these stars would confirm the physical reality of the groups, as well as allow us to distinguish whether their progenitors are dissolved clusters or accreted extragalactic systems

The effect of a barred potential (such as the one of the Milky Way) on the galactic orbits of forty-eight globular clusters for which absolute proper motions are known is studied. The orbital characteristics are compared with those obtained for the case of an axisymmetric galactic potential. Tidal radii are computed and discussed for both the better known axisymmetric case and that including a bar. The destruction rates due to bulge and disk shocking are calculated and compared in both galactic potentials.

A recent article has employed the determination from groundbased images of high proper motions in the Orion Nebula Cluster to argue that JW 349, JW 355, and JW 451 are high velocity (38 km/s, 89 km/s, and 69 km/s, respectively) low mass runaway stars. We report on measurement of the proper motions of these stars using images made by the Hubble Space Telescope's WFPC2 imager and find that there is no evidence for motions above 6.2 km/s for JW 349 and 7.9 km/s for JW 355, while the motion of 5.5 km/s for JW 451 is only slightly larger than the measurement uncertainty of 3.9 km/s. We conclude that there is no observational support for these stars being high velocity runaway stars.

We present absolute astrometry of 35 radio sources in the Orion Trapezium and Becklin-Neugebauer/Kleinman-Low regions, obtained from Very Large Array archival observations collected over a period of 15 years. By averaging the results for all the sources, we estimate the mean absolute proper motion of Orion to be --in Galactic coordinates-- $\mu_\ell \cos b$ = +2.1 $\pm$ 0.2 mas yr$^{-1}$; $\mu_b$ = $-$0.1 $\pm$ 0.2 mas yr$^{-1}$. These values agree remarkably well with those expected from the differential rotation of the Milky Way. Subtraction of this mean motion from the individual measurements allows us to register all proper motions to the rest frame of the Orion nebula, and identify radio sources with large residual velocities. In the KL region, we find three sources in this situation: the BN object, the radio source I, and the radio counterpart of the infrared source n. All three objects appear to be moving away from a common point where they must all have been located about 500 years ago. This suggests that all three sources were originally part of a multiple massive stellar system that recently disintegrated as a result of a close dynamical

In the course of a search for common proper motion binaries, in the Jones & Walker (JW) catalogue of proper motions in the Orion Nebula Cluster, we came across several faint stars with proper motions larger than one arcsecond per century and probabilities of membership P larger than 0.90. Such stars are interesting because they could be low-mass runaway stars recently accelerated by n-body interactions in compact multiple systems. Of particular interest among these stars is JW 451, which has a P = 0.98, the largest transverse velocity among all the stars with P larger or equal than 0.5 (69 km/s), and a proper motion vector which suggests that it was accelerated by the component C of the Trapezium some 1000 years ago. A closer examination of those JW stars with proper motion larger than one arcsec per century revealed that two other stars, JW 349 and JW 355 (with transverse velocities of 38 and 90 km/s respectively), in spite of being listed with P = 0 by JW, should also be considered part of the cluster, because these objects are also externally ionized proplyds. In fact, Hillenbrand (1997) assigns to them probabilities of membership of 0.99. Moreover, the proper motion errors of these two stars are relatively small, and so they are good candidates to be runaway stars recently accelerated in the Orion Nebula Cluster.

We present absolute astrometry of the core of the Orion molecular cloud, made with Very Large Array archive data taken over the last two decades. Our analysis reveals that both the BN object and the radio source I have proper motions: the BN object has a proper motion of $12.6 \pm 0.6$ mas yr$^{-1}$ (corresponding to a velocity of $27 \pm 1$ km s$^{-1}$ at an adopted distance of 450 pc) to the northwest, while the radio source I has a proper motion of $5.6 \pm 0.7$ mas yr$^{-1}$ (corresponding to a velocity of $12 \pm 2$ km s$^{-1}$) to the southeast. The motion of the two sources is nearly antiparallel, diverging from a point in between them, where they were located about 500 years ago. These results suggest that the BN object and the radio source I were part of a multiple young stellar system that disintegrated in the recent past.

We have collected from the literature a list of early-type stars, situated at large distances from the galactic plane, for which evidence of youth seems convincing. We discuss two possible formation mechanisms for these stars: ejection from the plane by dynamical interactions within small clusters, and formation away from the plane, via induced shocks created by spiral density waves. We identify the stars that could be explained by each mechanism. We conclude that the ejection mechanism can account for about two thirds of the stars, while a combination of star formation at z = 500-800 pc from the plane and ejection, can account for 90 percent of these stars. Neither mechanism, nor both together, can explain the most extreme examples.

RICE is designed to detect ultrahigh energy (>100 TeV) neutrinos from astrophysical sources. It will consist of an array of compact radio (100 to 1000 MHz) receivers buried in ice at the South Pole. During the 1995-96 and 1996-97 austral summers, several receivers transmitters were deployed in bore holes drilled for the AMANDA project, at depths of 141 to 260 m. This was the first in situ test of radio receivers in deep ice for neutrino astronomy.