Search

We present an in-depth analysis of a newly proposed correlation function in visibility space, between the E and B modes of the linear polarization, hereafter the EB-correlation, for a set of time-averaged GRMHD simulations compared with the phase map from different semi-analytic models as well as the Event Horizon Telescope (EHT) 2017 data for M87* source. We demonstrate that the phase map of the time-averaged EB-correlation contains novel information that might be linked to the BH spin, accretion state and the electron temperature. A detailed comparison with a semi-analytic approach with different azimuthal expansion modes shows that to recover the morphology of the real/imaginary part of the correlation function and its phase, we require higher orders of these azimuthal modes. To extract the phase features, we propose to use the Zernike polynomial reconstruction developing an empirical metric to break degeneracies between models with different BH spins that are qualitatively similar. We use a set of different geometrical ring models with various magnetic and velocity field morphologies and show that both the image space and visibility based EB-correlation morphologies in MAD simulations can be explained with simple fluid and magnetic field geometries as used in ring models. SANEs by contrast are harder to model, demonstrating that the simple fluid and magnetic field geometries of ring models are not sufficient to describe them owing to higher Faraday Rotation depths. A qualitative comparison with the EHT data demonstrates that some of the features in the phase of EB-correlation might be well explained by the current models for BH spins as well as electron temperatures, while others may require a larger theoretical surveys., Comment: 26 pages, 21 Figures

We explore the plasma matter content in the innermost accretion disk/jet in M87* as relevant for an enthusiastic search for the signatures of anti-matter in the next generation of the Event Horizon Telescope (ngEHT). We model the impact of non-zero positron-to-electron ratio using different emission models including a constant electron to magnetic pressure (constant $\beta_e$ model) with a population of non-thermal electrons as well as a R-beta model populated with thermal electrons. In the former case, we pick a semi-analytic fit to the force-free region of a general relativistic magnetohydrodynamic (GRMHD) simulation, while in the latter case, we analyze the GRMHD simulations directly. In both cases, positrons are being added at the post-processing level. We generate polarized images and spectra for some of these models and find out that at the radio frequencies, both of the linear and the circular polarizations get enhanced per adding pairs. On the contrary, we show that at higher frequencies a substantial positron fraction washes out the circular polarization. We report strong degeneracies between different emission models and the positron fraction, though our non-thermal models show more sensitivities to the pair fraction than the thermal models. We conclude that a large theoretical image library is indeed required to fully understand the trends probed in this study, and to place them in the context of large set of parameters which also affect polarimetric images, such as magnetic field strength, black hole spin, and detailed aspects of the electron temperature and the distribution function., Comment: 14 pages, 9 figures

We propose to trace the dynamical motion of a shearing hot spot near the SgrA* source through a dynamical image reconstruction algorithm, StarWarps. Such a hot spot may form as the exhaust of magnetic reconnection in a current sheet near the black hole horizon. A hot spot that is ejected from the current sheet into an orbit in the accretion disk may shear and diffuse due to instabilities at its boundary during its orbit, resulting in a distinct signature. We subdivide the motion to two distinct phases; the first phase refers to the appearance of the hot spot modelled as a bright blob, followed by a subsequent shearing phase simulated as a stretched ellipse. We employ different observational arrays, including EHT(2017,2022) and the next generation event horizon telescope (ngEHTp1, ngEHT) arrays, in which few new additional sites are added to the observational array. We make dynamical image reconstructions for each of these arrays. Subsequently, we infer the hot spot phase in the first phase followed by the axes ratio and the ellipse area in the second phase. We focus on the direct observability of the orbiting hot spot in the sub-mm wavelength. Our analysis demonstrates that newly added dishes may easily trace the first phase as well as part of the second phase, before the flux is reduced substantially. The algorithm used in this work can be extended to any other types of the dynamical motion. Consequently, we conclude that the ngEHT is a key to directly observe the dynamical motions near variable sources, such as SgrA*., Comment: 16 pages, 8 pages

We investigate general relativistic magnetohydrodynamic simulations (GRMHD) to determine the physical origin of the twisty patterns of linear polarization seen in spatially resolved black hole images and explain their morphological dependence on black hole spin. By characterising the observed emission with a simple analytic ring model, we find that the twisty morphology is determined by the magnetic field structure in the emitting region. Moreover, the dependence of this twisty pattern on spin can be attributed to changes in the magnetic field geometry that occur due to the frame dragging. By studying an analytic ring model, we find that the roles of Doppler boosting and lensing are subdominant. Faraday rotation may cause a systematic shift in the linear polarization pattern, but we find that its impact is subdominant for models with strong magnetic fields and modest ion-to-electron temperature ratios. Models with weaker magnetic fields are much more strongly affected by Faraday rotation and have more complicated emission geometries than can be captured by a ring model. However, these models are currently disfavoured by the recent EHT observations of M87*. Our results suggest that linear polarization maps can provide a probe of the underlying magnetic field structure around a black hole, which may then be usable to indirectly infer black hole spins. The generality of these results should be tested with alternative codes, initial conditions, and plasma physics prescriptions., Comment: 26 pages, 20 figures, accepted for publication in ApJ

We present an in-depth analysis of gas morphologies for a sample of 25 Milky Way-like galaxies from the IllustrisTNG TNG50 simulation. We constrain the morphology of cold, warm, hot gas, and gas particles as a whole using a Local Shell Iterative Method (LSIM) and explore its observational implications by computing the hard-to-soft X-ray ratio, which ranges between $10^{-3}$-$10^{-2}$ in the inner $\sim 50 \rm kpc$ of the distribution and $10^{-5}$-$10^{-4}$ at the outer portion of the hot gas distribution. We group galaxies into three main categories: simple, stretched, and twisted. These categories are based on the radial reorientation of the principal axes of the reduced inertia tensor. We find that a vast majority ($77\%$) of the galaxies in our sample exhibit twisting patterns in their radial profiles. Additionally, we present detailed comparisons between 1) the gaseous distributions belonging to individual temperature regimes, 2) the cold gas distributions and stellar distributions, and 3) the gaseous distributions and dark matter (DM) halos. We find a strong correlation between the morphological properties of the cold gas and stellar distributions. Furthermore, we find a correlation between gaseous distributions with DM halo that increases with gas temperature, implying that we may use the warm-hot gaseous morphology as a tracer to probe the DM morphology. Finally, we show gaseous distributions exhibit significantly more prolate morphologies than the stellar distributions and DM halos, which we hypothesize is due to stellar and AGN feedback., Comment: 26 pages, 12 figures

We analyze the velocity anisotropy of stars in real and energy space for a sample of Milky Way-like galaxies in the TNG50 simulation. We employ different selection criteria, including spatial, kinematic and metallicity cuts, and make three halo classes ($\mathcal{A}$-$\mathcal{C}$) which show mild-to-strong sensitivity to different selections. The above classes cover 48%, 16% and 36% of halos, respectively. We analyze the $\beta$ radial profiles and divide them into either monotonically increasing radial profiles or ones with peaks and troughs. We demonstrate that halos with monotonically increasing $\beta$ profiles are mostly from class $\mathcal{A}$, whilst those with peaks/troughs are part of classes $\mathcal{B}$-$\mathcal{C}$. This means that care must be taken as the observationally reported peaks/troughs might be a consequence of different selection criteria. We infer the anisotropy parameter $\beta$ energy space and compare that against the $\beta$ radial profile. It is seen that 65% of halos with very mild sensitivity to different selections in real space, are those for which the $\beta$ radial and energy profiles are closely related. Consequently, we propose that comparing the $\beta$ radial and energy profiles might be a novel way to examine the sensitivity to different selection criteria and thus examining the robustness of the anisotropy parameter in tracing stellar kinematics. We compare simulated $\beta$ radial profiles against various observations and demonstrate that, in most cases, the model diversity is comparable with the error bars from different observations, meaning that the TNG50 models are in good overall agreement with observations., Comment: 21 pages, 16 figures, Accepted for publication in ApJ

We investigate the morphology of the stellar distribution in a sample of Milky Way (MW) like galaxies in the TNG50 simulation. Using a local in shell iterative method (LSIM) as the main approach, we explicitly show evidence of twisting (in about 52% of halos) and stretching (in 48% of them) in the real space. This is matched with the re-orientation observed in the eigenvectors of the inertia tensor and gives us a clear picture of having a re-oriented stellar distribution. We make a comparison between the shape profile of dark matter (DM) halo and stellar distribution and quite remarkably see that their radial profiles are fairly close, especially at small galactocentric radii where the stellar disk is located. This implies that the DM halo is somewhat aligned with stars in response to the baryonic potential. The level of alignment mostly decreases away from the center. We study the impact of substructures in the orbital circularity parameter. It is demonstrated that in some cases, far away substructures are counter-rotating compared with the central stars and may flip the sign of total angular momentum and thus the orbital circularity parameter. Truncating them above 150 kpc, however, retains the disky structure of the galaxy as per initial selection. Including the impact of substructures in the shape of stars, we explicitly show that their contribution is subdominant. Overlaying our theoretical results to the observational constraints from previous literature, we establish fair agreement., Comment: 29 pages, 21 figures, Accepted for publication in ApJ

We present a comprehensive analysis of the shape of dark matter (DM) halos in a sample of 25 Milky Way-like galaxies in TNG50 simulation. Using an Enclosed Volume Iterative Method (EVIM), we infer an oblate-to-triaxial shape for the DM halo with the median $T \simeq 0.24 $. We group DM halos in 3 different categories. Simple halos (32% of population) establish principal axes whose ordering in magnitude does not change with radius and whose orientations are almost fixed throughout the halo. Twisted halos (32% of population), experience levels of gradual rotations throughout their radial profiles. Finally, stretched halos (36% of population) demonstrate a stretching in their principal axes lengths where the ordering of different eigenvalues change with radius. Subsequently, the halo experiences a "rotation" of $\sim$90 deg where the stretching occurs. Visualizing the 3D ellipsoid of each halo, for the first time, we report signs of re-orienting ellipsoid in twisted and stretched halos. We examine the impact of baryonic physics on DM halo shape through a comparison to dark matter only (DMO) simulations. This suggests a triaxial (prolate) halo. We analyze the impact of substructure on DM halo shape in both hydro and DMO simulations and confirm that their impacts are subdominant. We study the distribution of satellites in our sample. In simple and twisted halos, the angle of satellites' angular momentum with galaxy's angular momentum grows with radius. However, stretched halos show a flat distribution of angles. Overlaying our theoretical outcome on the observational results presented in the literature establishes a fair agreement., Comment: 32 pages, 28 plots, Accepted for publication in ApJ

We present a new measurement of the systemic proper motion of the Small Magellanic Cloud (SMC), based on an expanded set of 30 fields containing background quasars and spanning a $\sim$3 year baseline, using the \textit{Hubble Space Telescope} (\textit{HST}) Wide Field Camera 3. Combining this data with our previous 5 \textit{HST} fields, and an additional 8 measurements from the \textit{Gaia}-Tycho Astrometric Solution Catalog, brings us to a total of 43 SMC fields. We measure a systemic motion of $\mu_{W}$ = $-0.82$ $\pm$ 0.02 (random) $\pm$ 0.10 (systematic) mas yr$^{-1}$ and $\mu_{N}$ = $-1.21$ $\pm$ 0.01 (random) $\pm$ 0.03 (systematic) mas yr$^{-1}$. After subtraction of the systemic motion, we find little evidence for rotation, but find an ordered mean motion radially away from the SMC in the outer regions of the galaxy, indicating that the SMC is in the process of tidal disruption. We model the past interactions of the Clouds with each other based on the measured present-day relative velocity between them of $103 \pm 26$ km s$^{-1}$. We find that in 97\% of our considered cases, the Clouds experienced a direct collision $147 \pm 33$ Myr ago, with a mean impact parameter of $7.5 \pm 2.5$ kpc., Comment: 13 pages, 12 figures, 3 tables, accepted to ApJ

How black holes accrete surrounding matter is a fundamental, yet unsolved question in astrophysics. It is generally believed that matter is absorbed into black holes via accretion disks, the state of which depends primarily on the mass-accretion rate. When this rate approaches the critical rate (the Eddington limit), thermal instability is supposed to occur in the inner disc, causing repetitive patterns of large-amplitude X-ray variability (oscillations) on timescales of minutes to hours. In fact, such oscillations have been observed only in sources with a high mass accretion rate, such as GRS 1915+105. These large-amplitude, relatively slow timescale, phenomena are thought to have physical origins distinct from X-ray or optical variations with small amplitudes and fast ($\lesssim$10 sec) timescales often observed in other black hole binaries (e.g., XTE J1118+480 and GX 339-4). Here we report an extensive multi-colour optical photometric data set of V404 Cygni, an X-ray transient source containing a black hole of nine solar masses (and a conpanion star) at a distance of 2.4 kiloparsecs. Our data show that optical oscillations on timescales of 100 seconds to 2.5 hours can occur at mass-accretion rates more than ten times lower than previously thought. This suggests that the accretion rate is not the critical parameter for inducing inner-disc instabilities. Instead, we propose that a long orbital period is a key condition for these large-amplitude oscillations, because the outer part of the large disc in binaries with long orbital periods will have surface densities too low to maintain sustained mass accretion to the inner part of the disc. The lack of sustained accretion -- not the actual rate -- would then be the critical factor causing large-amplitude oscillations in long-period systems., Comment: Published in Nature on January 7th, 2016

We present proper motions for the Large & Small Magellanic Clouds (LMC & SMC) based on three epochs of \textit{Hubble Space Telescope} data, spanning a $\sim 7$ yr baseline, and centered on fields with background QSOs. The first two epochs, the subject of past analyses, were obtained with ACS/HRC, and have been reanalyzed here. The new third epoch with WFC3/UVIS increases the time baseline and provides better control of systematics. The three-epoch data yield proper motion random errors of only 1-2% per field. For the LMC this is sufficient to constrain the internal proper motion dynamics, as will be discussed in a separate paper. Here we focus on the implied center-of-mass proper motions: mu_W(LMC) = -1.910 +/- 0.020 mas/yr, mu_N(LMC) = 0.229 +/- 0.047 mas/yr, and mu_W(SMC) = -0.772 +/- 0.063 mas/yr, mu_N(SMC) = -1.117 +/- 0.061 mas/yr. We combine the results with a revised understanding of the solar motion in the Milky Way to derive Galactocentric velocities: v_{tot,LMC} = 321 +/- 24 km/s and v_{tot,SMC} = 217 +/- 26 km/s. Our proper motion uncertainties are now dominated by limitations in our understanding of the internal kinematics and geometry of the Clouds, and our velocity uncertainties are dominated by distance errors. Orbit calculations for the Clouds around the Milky Way allow a range of orbital periods, depending on the uncertain masses of the Milky Way and LMC. Periods $\lesssim 4$ Gyr are ruled out, which poses a challenge for traditional Magellanic Stream models. First-infall orbits are preferred (as supported by other arguments as well) if one imposes the requirement that the LMC and SMC must have been a bound pair for at least several Gyr., Comment: 27 pages, 15 figures, accepted to ApJ

We present 663 QSO candidates in the Large Magellanic Cloud (LMC) selected using multiple diagnostics. We started with a set of 2,566 QSO candidates from our previous work selected using time variability of the MACHO LMC lightcurves. We then obtained additional information for the candidates by crossmatching them with the Spitzer SAGE, the MACHO UBVI, the 2MASS, the Chandra and the XMM catalogs. Using this information, we specified six diagnostic features based on mid-IR colors, photometric redshifts using SED template fitting, and X-ray luminosities in order to further discriminate high confidence QSO candidates in the absence of spectra information. We then trained a one-class SVM (Support Vector Machine) model using the diagnostics features of the confirmed 58 MACHO QSOs. We applied the trained model to the original candidates and finally selected 663 high confidence QSO candidates. Furthermore, we crossmatched these 663 QSO candidates with the newly confirmed 144 QSOs and 275 non-QSOs in the LMC fields. On the basis of the counterpart analysis, we found that the false positive rate is less than 1%., Comment: 13 pages, 17 figures. accepted for publication in ApJ

We present a new QSO selection algorithm using a Support Vector Machine (SVM), a supervised classification method, on a set of extracted times series features including period, amplitude, color, and autocorrelation value. We train a model that separates QSOs from variable stars, non-variable stars and microlensing events using 58 known QSOs, 1,629 variable stars and 4,288 non-variables using the MAssive Compact Halo Object (MACHO) database as a training set. To estimate the efficiency and the accuracy of the model, we perform a cross-validation test using the training set. The test shows that the model correctly identifies ~80% of known QSOs with a 25% false positive rate. The majority of the false positives are Be stars. We applied the trained model to the MACHO Large Magellanic Cloud (LMC) dataset, which consists of 40 million lightcurves, and found 1,620 QSO candidates. During the selection none of the 33,242 known MACHO variables were misclassified as QSO candidates. In order to estimate the true false positive rate, we crossmatched the candidates with astronomical catalogs including the Spitzer Surveying the Agents of a Galaxy's Evolution (SAGE) LMC catalog and a few X-ray catalogs. The results further suggest that the majority of the candidates, more than 70%, are QSOs., Comment: 17 pages, 11 figures; accepted for the publication in ApJ

We study the dynamics of the Magellanic Clouds in a model for the Local Group whose mass is constrained using the timing argument/two-body limit of the action principle. The goal is to evaluate the role of M31 in generating the high angular momentum orbit of the Clouds, a puzzle that has only been exacerbated by the latest $HST$ proper motion measurements. We study the effects of varying the total Local Group mass, the relative mass of the Milky Way and M31, the proper motion of M31, and the proper motion of the LMC on this problem. Over a large part of this parameter-space we find that tides from M31 are insignificant. For a range of LMC proper motions approximately $3\sigma$ higher than the mean and total Local Group mass $> 3.5\times 10^{12} M_\odot$, M31 can provide a significant torque to the LMC orbit. However, if the LMC is bound to the MW, then M31 is found to have negligible effect on its motion and the origin of the high angular momentum of the system remains a puzzle. Finally, we use the timing argument to calculate the total mass of the MW-LMC system based on the assumption that they are encountering each other for the first time, their previous perigalacticon being a Hubble time ago, obtaining $M_{\rm MW} + M_{\rm LMC} = (8.7 \pm 0.8) \times 10^{11} M_\odot$., Comment: 10 pages, 5 figures, accepted to ApJ

Catalogs of periodic variable stars contain large numbers of periodic light-curves (photometric time series data from the astrophysics domain). Separating anomalous objects from well-known classes is an important step towards the discovery of new classes of astronomical objects. Most anomaly detection methods for time series data assume either a single continuous time series or a set of time series whose periods are aligned. Light-curve data precludes the use of these methods as the periods of any given pair of light-curves may be out of sync. One may use an existing anomaly detection method if, prior to similarity calculation, one performs the costly act of aligning two light-curves, an operation that scales poorly to massive data sets. This paper presents PCAD, an unsupervised anomaly detection method for large sets of unsynchronized periodic time-series data, that outputs a ranked list of both global and local anomalies. It calculates its anomaly score for each light-curve in relation to a set of centroids produced by a modified k-means clustering algorithm. Our method is able to scale to large data sets through the use of sampling. We validate our method on both light-curve data and other time series data sets. We demonstrate its effectiveness at finding known anomalies, and discuss the effect of sample size and number of centroids on our results. We compare our method to naive solutions and existing time series anomaly detection methods for unphased data, and show that PCAD's reported anomalies are comparable to or better than all other methods. Finally, astrophysicists on our team have verified that PCAD finds true anomalies that might be indicative of novel astrophysical phenomena.

We present a new approach to the analysis of time symmetry in light curves, such as those in the x-ray at the center of the Scorpius X-1 occultation debate. Our method uses a new parameterization for such events (the bilogistic event profile) and provides a clear, physically relevant characterization of each event's key features. We also demonstrate a Markov Chain Monte Carlo algorithm to carry out this analysis, including a novel independence chain configuration for the estimation of each event's location in the light curve. These tools are applied to the Scorpius X-1 light curves presented in Chang et al. (2007), providing additional evidence based on the time series that the events detected thus far are most likely not occultations by TNOs., Comment: 24 pages, 18 figures. Preprint typeset using LaTeX style emulateapj v. 04/20/08

We conducted a search for occultations of bright stars by Kuiper Belt Objects (KBOs) to estimate the density of sub-km KBOs in the sky. We report here the first results of this occultation survey of the outer solar system conducted in June 2007 and June/July 2008 at the MMT Observatory using Megacam, the large MMT optical imager. We used Megacam in a novel shutterless continuous--readout mode to achieve high precision photometry at 200 Hz. We present an analysis of 220 star hours at signal-to-noise ratio of 25 or greater. The survey efficiency is greater than 10% for occultations by KBOs of diameter d>=0.7 km, and we report no detections in our dataset. We set a new 95% confidence level upper limit for the surface density \Sigma_N(d) of KBOs larger than 1 km: \Sigma_N(d>=1 km) <= 2.0e8 deg^-2, and for KBOs larger than 0.7 km \Sigma_N(d>= 0.7 km) <= 4.8e8 deg^-2., Comment: 13 pages, 12 figures, submitted to AJ, modified fig 11 that did now display properly

We present a de-trending algorithm for the removal of trends in time series. Trends in time series could be caused by various systematic and random noise sources such as cloud passages, changes of airmass, telescope vibration or CCD noise. Those trends undermine the intrinsic signals of stars and should be removed. We determine the trends from subsets of stars that are highly correlated among themselves. These subsets are selected based on a hierarchical tree clustering algorithm. A bottom-up merging algorithm based on the departure from normal distribution in the correlation is developed to identify subsets, which we call clusters. After identification of clusters, we determine a trend per cluster by weighted sum of normalized light-curves. We then use quadratic programming to de-trend all individual light-curves based on these determined trends. Experimental results with synthetic light-curves containing artificial trends and events are presented. Results from other de-trending methods are also compared. The developed algorithm can be applied to time series for trend removal in both narrow and wide field astronomy., Comment: Revised version according to the referee's second review

The HST proper motion (PM) measurements of the Clouds have severe implications for their interaction history with the Milky Way (MW) and with each other. The Clouds are likely on their first passage about the MW and the SMC's orbit about the LMC is better described as quasi-periodic rather than circular. Binary L/SMC orbits that satisfy observational constraints on their mutual interaction history (e.g. the formation of the Magellanic Bridge during a collision between the Clouds ~300 Myr ago) can be located within 1 sigma of the mean PMs. However, these binary orbits are not co-located with the Magellanic Stream (MS) when projected on the plane of the sky and the line-of-sight velocity gradient along the LMC's orbit is significantly steeper than that along the MS. These combined results ultimately rule out a purely tidal origin for the MS: tides are ineffective without multiple pericentric passages and can neither decrease the velocity gradient nor explain the offset stream in a polar orbit configuration. Alternatively, ram pressure stripping of an extended gaseous disk may naturally explain the deviation. The offset also suggests that observations of the little-explored region between RA 21h and 23h are crucial for characterizing the full extent of the MS., Comment: 6 pages, 5 figures, to be published in proceedings of IAU Symposium No. 256: The Magellanic System: Stars, Gas, and Galaxies, Jacco Th. van Loon & Joana M. Oliveira, eds. 2009 (Cambridge: Cambridge University Press)

In HST Cycles 11 and 13 we obtained two epochs of ACS/HRC data for fields in the Magellanic Clouds centered on background quasars. We used these data to determine the proper motions of the LMC and SMC to better than 5% and 15% respectively. The results had a number of unexpected implications for the Milky Way-LMC-SMC system. The implied three-dimensional velocities were larger than previously believed and close to the escape velocity in a standard 10^12 solar mass Milky Way dark halo, implying that the Clouds may be on their first passage. Also, the relative velocity between the LMC and SMC was larger than expected, leaving open the possibility that the Clouds may not be bound to each other. To further verify and refine our results we requested an additional epoch of data in Cycle 16 which is being executed with WFPC2/PC due to the failure of ACS. We present the results of an ongoing analysis of these WFPC2 data which indicate good consistency with the two-epoch results., Comment: 6 pages, 4 figures. To be published in Proceedings of IAU Symposium 256, The Magellanic System: Stars, Gas, and Galaxies, Jacco Th. van Loon & Joana M. Oliveira, eds. 2009 (Cambridge: Cambridge University Press)

We present results from a study of the orbits of eclipsing binary stars (EBs) in the Magellanic Clouds. The samples comprise 4510 EBs found in the Large Magellanic Cloud (LMC) by the MACHO project, 2474 LMC EBs found by the OGLE-II project (of which 1182 are also in the MACHO sample), 1380 in the Small Magellanic Cloud (SMC) found by the MACHO project, and 1317 SMC EBs found by the OGLE-II project (of which 677 are also in the MACHO sample); we also consider the EROS sample of 79 EBs in the bar of the LMC. Statistics of the phase differences between primary and secondary minima allow us to infer the statistics of orbital eccentricities within these samples. We confirm the well-known absence of eccentric orbit in close binary stars. We also find evidence for rapid circularization in longer period systems when one member evolves beyond the main sequence, as also found by previous studies., Comment: 37 pages, 16 figures, accepted for publication in ApJ. Added a new reference and updated information on on line material

We present a new sample of 4634 eclipsing binary stars in the Large Magellanic Cloud (LMC), expanding on a previous sample of 611 objects and a new sample of 1509 eclipsing binary stars in the Small Magellanic Cloud (SMC), that were identified in the light curve database of the MACHO project. We perform a cross correlation with the OGLE-II LMC sample, finding 1236 matches. A cross correlation with the OGLE-II SMC sample finds 698 matches. We then compare the LMC subsamples corresponding to center and the periphery of the LMC and find only minor differences between the two populations. These samples are sufficiently large and complete that statistical studies of the binary star populations are possible., Comment: 67 pages, 40 figures

Recent proper motion measurements of the Large and Small Magellanic Clouds (LMC and SMC, respectively) by Kallivayalil et al (2006a,b) suggest that the 3D velocities of the Clouds are substantially higher (~100 km/s) than previously estimated and now approach the escape velocity of the Milky Way (MW). Previous studies have also assumed that the Milky Way can be adequately modeled as an isothermal sphere to large distances. Here we re-examine the orbital history of the Clouds using the new velocities and a LCDM-motivated MW model with virial mass Mvir = 1e12 Msun (e.g. Klypin et al 2002). We conclude that the Clouds are either currently on their first passage about the MW or, if the MW can be accurately modeled by an isothermal sphere to distances >200 kpc (i.e., Mvir > 2e12 Msun), that their orbital period and apogalacticon distance must be a factor of two larger than previously estimated, increasing to 3 Gyr and 200 kpc, respectively. A first passage scenario is consistent with the fact that the LMC and SMC appear to be outliers when compared to other satellite galaxies of the MW: they are irregular in appearance and are moving faster. We discuss the implications of this orbital analysis for our understanding of the star formation history, the nature of the warp in the MW disk and the origin of the Magellanic Stream (MS), a band of HI gas trailing the LMC and SMC that extends roughly 100 degrees across the sky. Specifically, as a consequence of the new orbital history of the Clouds, the origin of the MS may not be explainable by current tidal and ram pressure stripping models.

We report on the identification of the lens responsible for microlensing event MACHO-LMC-20. As part of a \textit{Spitzer}/IRAC program conducting mid-infrared follow-up of the MACHO Large Magellanic Cloud microlensing fields, we discovered a significant flux excess at the position of the source star for this event. These data, in combination with high resolution near-infrared \textit{Magellan}/PANIC data has allowed us to classify the lens as an early M dwarf in the thick disk of the Milky Way, at a distance of $\sim 2$ kpc. This is only the second microlens to have been identified, the first also being a M dwarf star in the disk. Together, these two events are still consistent with the expected frequency of nearby stars in the Milky Way thin and thick disks acting as lenses., Comment: 6 pages, 4 figures, submitted to ApJ Letters

We present a measurement of the systemic proper motion of the Small Magellanic Cloud (SMC) made using the Advanced Camera for Surveys (ACS) on the \textit{Hubble Space Telescope} (\textit{HST}). We tracked the SMC's motion relative to 4 background QSOs over a baseline of approximately 2 years. The measured proper motion is : $\mu_W = -1.16 \pm 0.18 \masyr, \mu_N = -1.17 \pm 0.18 \masyr$. This is the best measurement yet of the SMC's proper motion. We combine this new result with our prior estimate of the proper motion of the Large Magellanic Cloud (LMC) from the same observing program to investigate the orbital evolution of both Clouds over the past 9 Gyr. The current relative velocity between the Clouds is $105 \pm 42 \kms$. Our investigations of the past orbital motions of the Clouds in a simple model for the dark halo of the Milky Way imply that the Clouds could be unbound from each other. However, our data are also consistent with orbits in which the Clouds have been bound to each other for approximately a Hubble time. Smaller proper motion errors and better understanding of the LMC and SMC masses would be required to constrain their past orbital history and their bound vs. unbound nature unambiguously. The new proper motion measurements should be sufficient to allow the construction of improved models for the origin and properties of the Magellanic Stream. In turn, this will provide new constraints on the properties of the Milky Way dark halo., Comment: Accepted for publication in ApJ; 41 pages, 15 figures

We present a measurement of the systemic proper motion of the Large Magellanic Cloud (LMC) from astrometry with the High Resolution Camera (HRC) of the Advanced Camera for Surveys (ACS) on the Hubble Space Telescope (HST). We observed LMC fields centered on 21 background QSOs that were discovered from their optical variability in the MACHO database. The QSOs are distributed homogeneously behind the central few degrees of the LMC. With 2 epochs of HRC data and a ~2 year baseline we determine the proper motion of the LMC to better than 5% accuracy: mu_W = -2.03 +/- 0.08 mas/yr; mu_N = 0.44 +/- 0.05 mas/yr. This is the most accurate proper motion measurement for any Milky Way satellite thus far. When combined with HI data from the Magellanic Stream this should provide new constraints on both the mass distribution of the Galactic Halo and models of the Stream., Comment: 40 pages, 15 figures, submitted to ApJ

The light curve of a blue variable in the MACHO LMC database (FTS ID 78.5979.72) appeared nearly unvarying for ~4 years (quasi-flat segment) but then rapidly changed to become periodic with noisy minima for the remaining 4 years (periodic segment); there are no antecedent indications of a gradual approach to this change. Lomb Periodogram analyses indicate the presence of two distinct periods of ~61 days and 8 days in both the quasi-flat and the periodic segments. Minima of the periodic segment cover at least 50% of the orbital period and contain spikes of light with the 8-day period; maxima do not show this short period. The system typically shows maxima to be redder than minima. The most recent OGLE-III light curve shows only a 30-day periodicity. The variable's V and R magnitudes and color are those of a Be star, and recent sets of near infrared spectra four days apart, secured during the time of the OGLE-III data, show H-alpha emission near and at a maximum, confirming its Be star characteristics. The model that best fits the photometric behavior consists of a thin ring-like circumstellar disk of low mass with four obscuring sectors orbiting the central B star in unison at the 61-day period. The central star peers through the three equi- spaced separations between the four sectors producing the 8-day period. The remainder of the disk contains hydrogen in emission making maxima appear redder. A companion star of lower mass in an inclined and highly eccentric orbit produces an impulsive perturbation near its periastron to change the disk's orientation, changing eclipses from partial to complete within ~ 10 days., Comment: 42 pages, 14 figures, and 2 tables Submitted to AJ v3: Title changed, figures added, model modified

We present an algorithm that allows fast and efficient detection of transits, including planetary transits, from light-curves. The method is based on building an ensemble of fiducial models and compressing the data using the MOPED algorithm. We describe the method and demonstrate its efficiency by finding planet-like transits in simulated Pan-STARRS light-curves. We show that that our method is independent of the size of the search space of transit parameters. In large sets of light-curves, we achieve speed up factors of order of $10^{8}$ times over the full $\chi2$ search. We discuss how the algorithm can be used in forthcoming large surveys like Pan-STARRS and LSST and how it may be optimized for future space missions like Kepler and COROT where most of the processing must be done on board., Comment: 9 pages, 9 figures

We have carried out photometry of the microlensing event MACHO-LMC-5 with Spitzer's IRAC ten years after the magnification of the LMC source star was recorded. This event is unique in the annals of gravitational microlensing: the lensing star itself has been observed using HST (once with WFPC2 and twice with ACS/HRC). Since the separation between the source and lens at the epoch of the Spitzer observations was $\sim0.24''$, the two stars cannot be resolved in the Spitzer images. However, the IRAC photometry clearly establishes that the lens is a M5 dwarf star from its infrared excess, which in turn yields a mass of $\sim0.2 M_{\odot}$. This demonstrates the potential of Spitzer to detect the lenses in other gravitational microlensing events., Comment: Accepted for publication, ApJS Spitzer Special Issue

The Taiwanese-American Occultation Survey (TAOS) will detect objects in the Kuiper Belt, by measuring the rate of occultations of stars by these objects, using an array of three to four 50cm wide-field robotic telescopes. Thousands of stars will be monitored, resulting in hundreds of millions of photometric measurements per night. To optimize the success of TAOS, we have investigated various methods of gathering and processing the data and developed statistical methods for detecting occultations. In this paper we discuss these methods. The resulting estimated detection efficiencies will be used to guide the choice of various operational parameters determining the mode of actual observation when the telescopes come on line and begin routine observations. In particular we show how real-time detection algorithms may be constructed, taking advantage of having multiple telescopes. We also discuss a retrospective method for estimating the rate at which occultations occur., Comment: 14

We have performed a search for halo white dwarfs as high proper motion objects in a second epoch WFPC2 image of the Groth-Westphal strip. We identify 24 high proper motion objects with mu > 0.014 ''/yr. Five of these high proper motion objects are identified as strong white dwarf candidates on the basis of their position in a reduced proper motion diagram. We create a model of the Milky Way thin disk, thick disk and stellar halo and find that this sample of white dwarfs is clearly an excess above the < 2 detections expected from these known stellar populations. The origin of the excess signal is less clear. Possibly, the excess cannot be explained without invoking a fourth galactic component: a white dwarf dark halo. We present a statistical separation of our sample into the four components and estimate the corresponding local white dwarf densities using only the directly observable variables, V, V-I, and mu. For all Galactic models explored, our sample separates into about 3 disk white dwarfs and 2 halo white dwarfs. However, the further subdivision into the thin and thick disk and the stellar and dark halo, and the subsequent calculation of the local densities are sensitive to the input parameters of our model for each Galactic component. Using the lowest mean mass model for the dark halo we find a 7% white dwarf halo and six times the canonical value for the thin disk white dwarf density (at marginal statistical significance), but possible systematic errors due to uncertainty in the model parameters likely dominate these statistical error bars. The white dwarf halo can be reduced to around 1.5% of the halo dark matter by changing the initial mass function slightly. The local thin disk white dwarf density in our solution can be made consistent with the canonical value by assuming a larger thin disk scaleheight of 500 pc., Comment: revised version, accepted by ApJ, results unchanged, discussion expanded

We entertain the idea that robust theoretical expectations can become a tool in removing hidden observational or data-reduction biases. We illustrate this approach for a specific problem associated with gravitational microlensing. Using the fact that a group is more than just a collection of individuals, we derive formulae for correcting the distribution of the dimensionless impact parameters of events, u_min. We refer to the case when undetected biases in the u_min distribution can be alleviated by multiplication of impact parameters of all events by a common constant factor. We show that in this case the general maximum likelihood problem of solving an infinite number of equations reduces to two constraints, and we find an analytic solution. Under the above assumptions, this solution represents a state in which the ``entropy'' of a microlensing ensemble is at its maximum, that is, the distribution of u_min resembles a specific, theoretically expected, box-like distribution to the highest possible extent. We also show that this technique does not allow one to correct the parameters of individual events on the event by event basis independently from each other., Comment: 16 pages, version accepted by ApJ, results unchanged, additional discussion regarding conditions suitable for application of the presented method

The accretion of comets onto DA white dwarfs can produce observable metal absorption lines. We show here that comet systems around the progenitor main sequence star are vulnerable to being lost during asymptotic giant branch mass loss, if the mass loss is sufficiently asymmetric to impart modest linear momentum to the white dwarf. This may have bearing on the frequency of observation of heavy elements in white dwarf stars, and on inferences regarding the frequency of comet systems, if the imparted linear velocities of white dwarfs can be estimated., Comment: 19 pages, 8 figures, AASTeX-4.0, to be published in ApJ (7/1/98); replaced with slight modifications to intro text and references

We present a community-led assessment of the solar system investigations achievable with NASA's next-generation space telescope, the Wide Field Infrared Survey Telescope (WFIRST). WFIRST will provide imaging, spectroscopic, and coronagraphic capabilities from 0.43 to 2.0 m and will be a potential contemporary and eventual successor to the James Webb Space Telescope (JWST). Surveys of irregular satellites and minor bodies are where WFIRST will excel with its 0.28 deg2 field-of-view Wide Field Instrument. Potential ground- breaking discoveries from WFIRST could include detection of the first minor bodies orbiting in the inner Oort Cloud, identification of additional Earth Trojan asteroids, and the discovery and characterization of asteroid binary systems similar to Ida/Dactyl. Additional investigations into asteroids, giant planet satellites, Trojan asteroids, Centaurs, Kuiper belt objects, and comets are presented. Previous use of astrophysics assets for solar system science and synergies between WFIRST, Large Synoptic Survey Telescope, JWST, and the proposed Near-Earth Object Camera mission is discussed. We also present the case for implementation of moving target tracking, a feature that will benefit from the heritage of JWST and enable a broader range of solar system observations.