We present the results of a detailed investigation into the physical conditions in interstellar material interacting with the supernova remnant IC 443. Our analysis is based on a comprehensive examination of high-resolution far-ultraviolet spectra obtained with the Space Telescope Imaging Spectrograph onboard the Hubble Space Telescope of two stars behind IC 443. One of our targets (HD 43582) probes gas along the entire line of sight through the supernova remnant, while the other (HD 254755) samples material located ahead of the primary supernova shock front. We identify low velocity quiescent gas in both directions and find that the densities and temperatures in these components are typical of diffuse atomic and molecular clouds. Numerous high velocity components are observed in the absorption profiles of neutral and singly-ionized atomic species toward HD 43582. These components exhibit a combination of greatly enhanced thermal pressures and significantly reduced dust-grain depletions. We interpret this material as cooling gas in a recombination zone far downstream from shocks driven into neutral gas clumps. The pressures derived for a group of ionized gas components at high positive velocity toward HD 43582 are lower than those of the other shocked components, pointing to pressure inhomogeneities across the remnant. A strong very high velocity component near -620 km/s is seen in the absorption profiles of highly-ionized species toward HD 43582. The velocity of this material is consistent with the range of shock velocities implied by observations of soft thermal X-ray emission from IC 443. Moderately high-velocity gas toward HD 254755 may represent shocked material from a separate foreground supernova remnant.

The observed internal plateau of X-ray emission in some short GRBs suggests the formation of a remnant supra-massive magnetar following a double neutron star (NS) merger. In this paper, we assume that the rotational energy is lost mainly via gravitational wave (GW) radiation instead of magnetic dipole (MD) radiation, and present further constraints on the NS nuclear equation of state (EoS) via mass quadrupole deformation and r-mode fluid oscillations of the magnetar. We present two short GRBs with measured redshifts, 101219A and 160821B, whose X-ray light curves exhibit an internal plateau. This suggests that a supra-massive NS may survive as the central engine. By considering twelve NS EoSs, within the mass quadrupole deformation scenario we find that the GM1, DD2, and DDME2 models give an $M_{\rm p}$ band falling within the 2$\sigma$ region of the proto-magnetar mass distribution for $\varepsilon=0.01$. This is consistent with the constraints from the MD radiation dominated model of rotational energy loss. However, for an r-mode fluid oscillation model with $\alpha=0.1$ the data suggest that the NS EOS is close to the Shen and APR models, which is obviously different from the MD radiation dominated and mass quadrupole deformation cases.

Short GRBs with extended emission (EE) that are composed initially of a short-hard spike and followed by a long-lasting EE, are thought to be classified as a subsection of short GRBs. The narrow energy band available during the {\em Swift} era combined with a lack of spectral information prevented discovery of the intrinsic properties of those events. In this paper, we performed a systematic search of short GRBs with EE by using all available {\em Fermi}/GBM data. The search identified 26 GBM-detected short GRBs with EE that are similar to GRB 060614 observed by {\em Swift}/BAT. We focus on investigating the spectral and temporal properties for both the hard spike and the EE components of all 26 GRBs, and explore differences and possible correlations between them. We find that while the peak energy ($E_{\rm p}$) of the hard spikes is a little bit harder than that of the EE, but their fluences are comparable. The harder $E_{\rm p}$ seems to correspond to a larger fluence and peak flux with a large scatter for both the hard spike and EE components. Moreover, the $E_{\rm p}$ of both the hard spikes and EE are compared to other short GRBs. Finally, we also compare the properties of GRB 170817A with those short GRBs with EE and find no significant statistical differences between them. We find that GRB 170817A has the lowest $E_{\rm p}$, likely because it was off-axis.

In response to the need for the Astro2020 Decadal Survey to explicitly engage early career astronomers, the National Academies of Sciences, Engineering, and Medicine hosted the Early Career Astronomer and Astrophysicist Focus Session (ECFS) on October 8-9, 2018 under the auspices of Committee of Astronomy and Astrophysics. The meeting was attended by fifty six pre-tenure faculty, research scientists, postdoctoral scholars, and senior graduate students, as well as eight former decadal survey committee members, who acted as facilitators. The event was designed to educate early career astronomers about the decadal survey process, to solicit their feedback on the role that early career astronomers should play in Astro2020, and to provide a forum for the discussion of a wide range of topics regarding the astrophysics career path. This white paper presents highlights and themes that emerged during two days of discussion. In Section 1, we discuss concerns that emerged regarding the coming decade and the astrophysics career path, as well as specific recommendations from participants regarding how to address them. We have organized these concerns and suggestions into five broad themes. These include (sequentially): (1) adequately training astronomers in the statistical and computational techniques necessary in an era of "big data", (2) responses to the growth of collaborations and telescopes, (3) concerns about the adequacy of graduate and postdoctoral training, (4) the need for improvements in equity and inclusion in astronomy, and (5) smoothing and facilitating transitions between early career stages. Section 2 is focused on ideas regarding the decadal survey itself, including: incorporating early career voices, ensuring diverse input from a variety of stakeholders, and successfully and broadly disseminating the results of the survey.

The event GW170817/GRB 170817A, discovered via the successful joint observation of its gravitational wave radiation and its multi-wavelength electromagnetic counterparts, was the first definite "smoking-gun" from the merger of two neutron stars (NSs). However, the remnant of the merger remains unknown. Piro et al. recently claimed that a low-significance X-ray variability in GRB 170817A. By systematically comparing the properties of variability in the afterglow of GRB 170817A and X-ray flares in GRB afterglows, we find that this X-ray variability seems to share similar statistical correlations with X-ray flares in GRB afterglows. We further investigate several possible merger product scenarios to see whether they can produce the observed X-ray variability in GRB 170817A. The first scenario invokes a stable magnetar as the central engine producing the later X-ray variability via differential rotation or fall-back accretion onto the NS. The second scenario invokes a black hole as the central engine with a fall-back accretion process. The final scenario is a central engine with a long-lived supra-massive NS. We find that the first two scenarios have difficulty producing the later X-ray variability, which requires either an impractical NS magnetic field or an extraordinarily large stellar envelope and an extremely long accretion timescale. However, the third scenario seems to be consistent with observations, and the later X-ray variability can be produced by the magnetosphere which is expelled following the collapse of the NS with a $B_p\in(3.6, 13.5)\times10^{13}$ G.

Two-episode emission components separated by quiescent gaps in the prompt emission of gamma-ray bursts (GRBs) have been observed in the {\em Swift} era, but there is a lack of spectral information due to the narrow energy band of the {\em Swift}/Burst Alert Telescope. In this paper, a systematic analysis of the spectral and temporal properties of the prompt emission of 101 {\em Fermi}/Gamma-ray Burst Monitor detected long GRBs show the existence of two-episode emission components in the light curves, with quiescent times of up to hundreds of seconds. We focus on investigating the differences of those two emission episodes. We find that the light curves of the two emission components exhibit different behavior, e.g., a soft emission component that either precedes or follows the main prompt emission or that the intensity of the two emission episodes are comparable with each other. No statistically significant correlation in the duration of the two emission episodes can be claimed. We define a new parameter $\varepsilon$ as the ratio of the peak flux of the first and second emission episodes and find that a higher $\varepsilon$ corresponds to a larger fluence. The preferred spectral model in our analysis is a cutoff power-law model for most GRBs. The distribution of $E_p$ for episodes I and II range from tens of keV to 1000 keV with a lognormal fit and there are no significant differences between them. Moreover, we do not find significant relationships between $\varepsilon$ and $E_p$ for the two emission episodes. Those results suggest that these two-episode emission components likely share the same physical origin.

The calibration hardware system of the Large Synoptic Survey Telescope (LSST) is designed to measure two quantities: a telescope's instrumental response and atmospheric transmission, both as a function of wavelength. First of all, a "collimated beam projector" is designed to measure the instrumental response function by projecting monochromatic light through a mask and a collimating optic onto the telescope. During the measurement, the light level is monitored with a NIST-traceable photodiode. This method does not suffer from stray light effects or the reflections (known as ghosting) present when using a flat-field screen illumination, which has a systematic source of uncertainty from uncontrolled reflections. It allows for an independent measurement of the throughput of the telescope's optical train as well as each filter's transmission as a function of position on the primary mirror. Second, CALSPEC stars can be used as calibrated light sources to illuminate the atmosphere and measure its transmission. To measure the atmosphere's transfer function, we use the telescope's imager with a Ronchi grating in place of a filter to configure it as a low resolution slitless spectrograph. In this paper, we describe this calibration strategy, focusing on results from a prototype system at the Cerro Tololo Inter-American Observatory (CTIO) 0.9 meter telescope. We compare the instrumental throughput measurements to nominal values measured using a laboratory spectrophotometer, and we describe measurements of the atmosphere made via CALSPEC standard stars during the same run.

Using visible, radio, microwave, and sub-mm data, we study several lines of sight toward stars generally closer than 1 kpc on a component-by-component basis. We derive the component structure seen in absorption at visible wavelengths from Ca II, Ca I, K I, CH, CH$^{+}\!,$ and CN and compare it to emission from H I, CO and its isotopologues, and C$^{+}$ from the GOT C+ survey. The correspondence between components in emission and absorption help create a more unified picture of diffuse atomic and molecular gas in the interstellar medium. We also discuss how these tracers are related to the CO-dark H$_{2}$ gas probed by C$^{+}$ emission and discuss the kinematic connections among the species observed.

The cosmological evolution of primordial black holes (PBHs) is considered. A comprehensive view of the accretion and evaporation histories of PBHs across the entire cosmic history is presented, with focus on the critical mass holes. The critical mass of a PBH for current era evaporation is $M_{cr}\sim 5.1\times10^{14}$ g. Across cosmic time such a black hole will not accrete radiation or matter in sufficient quantity to hasten the inevitable evaporation, if the black hole remains within an average volume of the universe. The accretion rate onto PBHs is most sensitive to the mass of the hole, the sound speed in the cosmological fluid, and the energy density of the accreted components. It is not easy for a PBH to accrete the average cosmological fluid to reach $30M_\odot$ by $z\sim0.1$, the approximate mass and redshift of the merging BHs that were the sources of the gravitational wave events GW150914 and GW151226. A PBH located in an overdense region can undergo enhanced accretion leading to the possibility of growing by many orders of magnitude across cosmic history. Thus, two merging PBHs are a plausible source for the observed gravitational wave events. However, it is difficult for isolated PBHs to grow to supermassive black holes (SMBHs) at high redshift with masses large enough to fit observational constraints.

GRB 160821B is a short gamma-ray burst (GRB) at redshift $z=0.16$, with a duration less than 1 second and without detection of any "extended emission" up to more than 100 seconds in both {\em Swift}/BAT and {\em Fermi}/GBM bands. An X-ray plateau with a sharp drop 180 seconds after the BAT trigger was observed with {\em Swift}/XRT. No supernova or kilo-nova signature was detected. Assuming the central engine of this SGRB is a recently born supra-massive magnetar, we can explain the SGRB as jet radiation and its X-ray plateau as the internal energy dissipation of the pulsar wind as it spins down. We constrain its surface magnetic field as $B_{\rm p}<3.12\times 10^{16}$ G and initial spin period as $P_0< 8.5\times 10^{-3}$ seconds. Its equation of state is consistent with the GM1 model with $M_{\rm TOV} \sim 2.37 M_\odot$ and ellipticity $\epsilon<0.07$. Its gravitational wave (GW) radiation may be detectable with the future Einstein Telescope, but is much weaker than the current detectability limit of advanced-LIGO. The GW radiation of such an event would be detectable by advanced-LIGO if it occurred at a distance of 100 Mpc ($z=0.023$).

The rebrightening phenomenon is an interesting feature in some X-ray, optical, and radio afterglows of gamma-ray bursts (GRBs). Here, we propose a possible energy-supply assumption to explain the rebrightenings of radio afterglows, in which the central engine with multiple active phases can supply at least two GRB pulses in a typical GRB duration time. Considering the case of double pulses supplied by the central engine, the double pulses have separate physical parameters, except for the number density of the surrounding interstellar medium (ISM). Their independent radio afterglows are integrated by the ground detectors to form the rebrightening phenomenon. In this work, we firstly simulate diverse rebrightening light curves under consideration of different and independent physical parameters. Using this assumption, we also give our best fit to the radio afterglow of GRB 970508 at three frequencies of 1.43, 4.86, and 8.46 GHz. We suggest that the central engine may be active continuously at a timescale longer than that of a typical GRB duration time as many authors have suggested, and that it may supply enough energy to cause the long-lasting rebrightenings observed in some GRB afterglows.

The Taiwanese-American Occultation Survey (TAOS) aims to detect serendipitous occultations of stars by small (about 1 km diameter) objects in the Kuiper Belt and beyond. Such events are very rare (<0.001 events per star per year) and short in duration (about 200 ms), so many stars must be monitored at a high readout cadence. TAOS monitors typically around 500 stars simultaneously at a 5 Hz readout cadence with four telescopes located at Lulin Observatory in central Taiwan. In this paper, we report the results of the search for small Kuiper Belt Objects (KBOs) in seven years of data. No occultation events were found, resulting in a 95% c.l. upper limit on the slope of the faint end of the KBO size distribution of q = 3.34 to 3.82, depending on the surface density at the break in the size distribution at a diameter of about 90 km.

We present a detailed temperature and magnetic investigation of the T Tauri star V410 Tau by means of a simultaneous Doppler- and Zeeman-Doppler Imaging. Moreover we introduce a new line profile reconstruction method based on a singular value decomposition (SVD) to extract the weak polarized line profiles. One of the key features of the line profile reconstruction is that the SVD line profiles are amenable to radiative transfer modeling within our Zeeman-Doppler Imaging code iMap. The code also utilizes a new iterative regularization scheme which is independent of any additional surface constraints. To provide more stability a vital part of our inversion strategy is the inversion of both Stokes I and Stokes V profiles to simultaneously reconstruct the temperature and magnetic field surface distribution of V410 Tau. A new image-shear analysis is also implemented to allow the search for image and line profile distortions induced by a differential rotation of the star. The magnetic field structure we obtain for V410 Tau shows a good spatial correlation with the surface temperature and is dominated by a strong field within the cool polar spot. The Zeeman-Doppler maps exhibit a large-scale organization of both polarities around the polar cap in the form of a twisted bipolar structure. The magnetic field reaches a value of almost 2 kG within the polar region but smaller fields are also present down to lower latitudes. The pronounced non-axisymmetric field structure and the non-detection of a differential rotation for V410 Tau supports the idea of an underlying $\alpha^2$-type dynamo, which is predicted for weak-lined T Tauri stars.

Modern time-domain surveys continuously monitor large swaths of the sky to look for astronomical variability. Astrophysical discovery in such data sets is complicated by the fact that detections of real transient and variable sources are highly outnumbered by bogus detections caused by imperfect subtractions, atmospheric effects and detector artefacts. In this work we present a machine learning (ML) framework for discovery of variability in time-domain imaging surveys. Our ML methods provide probabilistic statements, in near real time, about the degree to which each newly observed source is astrophysically relevant source of variable brightness. We provide details about each of the analysis steps involved, including compilation of the training and testing sets, construction of descriptive image-based and contextual features, and optimization of the feature subset and model tuning parameters. Using a validation set of nearly 30,000 objects from the Palomar Transient Factory, we demonstrate a missed detection rate of at most 7.7% at our chosen false-positive rate of 1% for an optimized ML classifier of 23 features, selected to avoid feature correlation and over-fitting from an initial library of 42 attributes. Importantly, we show that our classification methodology is insensitive to mis-labelled training data up to a contamination of nearly 10%, making it easier to compile sufficient training sets for accurate performance in future surveys. This ML framework, if so adopted, should enable the maximization of scientific gain from future synoptic survey and enable fast follow-up decisions on the vast amounts of streaming data produced by such experiments.

Efficient and automated classification of periodic variable stars is becoming increasingly important as the scale of astronomical surveys grows. Several recent papers have used methods from machine learning and statistics to construct classifiers on databases of labeled, multi--epoch sources with the intention of using these classifiers to automatically infer the classes of unlabeled sources from new surveys. However, the same source observed with two different synoptic surveys will generally yield different derived metrics (features) from the light curve. Since such features are used in classifiers, this survey-dependent mismatch in feature space will typically lead to degraded classifier performance. In this paper we show how and why feature distributions change using OGLE and \textit{Hipparcos} light curves. To overcome survey systematics, we apply a method, \textit{noisification}, which attempts to empirically match distributions of features between the labeled sources used to construct the classifier and the unlabeled sources we wish to classify. Results from simulated and real--world light curves show that noisification can significantly improve classifier performance. In a three--class problem using light curves from \textit{Hipparcos} and OGLE, noisification reduces the classifier error rate from 27.0% to 7.0%. We recommend that noisification be used for upcoming surveys such as Gaia and LSST and describe some of the promises and challenges of applying noisification to these surveys.

We present a first and detailed study of the bright and active K0IV-III star HD 123351. The star is found to be a single-lined spectroscopic binary with a period of 147.8919+-0.0003 days and a large eccentricity of e=0.8086+-0.0001. The rms of the orbital solution is just 47 m/s, making it the most precise orbit ever obtained for an active binary system. The rotation period is constrained from long-term photometry to be 58.32+-0.01 days. It shows that HD 123351 is a very asynchronous rotator, rotating five times slower than the expected pseudo-synchronous value. Two spotted regions persisted throughout the 12 years of our observations. Four years of Halpha, CaII H&K and HeI D3 monitoring identifies the same main periodicity as the photometry but dynamic spectra also indicate that there is an intermittent dependence on the orbital period, in particular for Ca ii H&K in 2008. Line-profile inversions of a pair of Zeeman sensitive/insensitive iron lines yield an average surface magnetic-flux density of 542+-72 G. The time series for 2008 is modulated by the stellar rotation as well as the orbital motion, such that the magnetic flux is generally weaker during times of periastron and that the chromospheric emissions vary in anti-phase with the magnetic flux. We also identify a broad and asymmetric lithium line profile and measure an abundance of log n(Li) = 1.70+-0.05. The star's position in the H-R diagram indicates a mass of 1.2+-0.1 Msun and an age of 6-7 Gyr. We interpret the anti-phase relation of the magnetic flux with the chromospheric emissions as evidence that there are two magnetic fields present at the same time, a localized surface magnetic field associated with spots and a global field that is oriented towards the (low-mass) secondary component.

Despite the great promise of machine-learning algorithms to classify and predict astrophysical parameters for the vast numbers of astrophysical sources and transients observed in large-scale surveys, the peculiarities of the training data often manifest as strongly biased predictions on the data of interest. Typically, training sets are derived from historical surveys of brighter, more nearby objects than those from more extensive, deeper surveys (testing data). This sample selection bias can cause catastrophic errors in predictions on the testing data because a) standard assumptions for machine-learned model selection procedures break down and b) dense regions of testing space might be completely devoid of training data. We explore possible remedies to sample selection bias, including importance weighting (IW), co-training (CT), and active learning (AL). We argue that AL---where the data whose inclusion in the training set would most improve predictions on the testing set are queried for manual follow-up---is an effective approach and is appropriate for many astronomical applications. For a variable star classification problem on a well-studied set of stars from Hipparcos and OGLE, AL is the optimal method in terms of error rate on the testing data, beating the off-the-shelf classifier by 3.4% and the other proposed methods by at least 3.0%. To aid with manual labeling of variable stars, we developed a web interface which allows for easy light curve visualization and querying of external databases. Finally, we apply active learning to classify variable stars in the ASAS survey, finding dramatic improvement in our agreement with the ACVS catalog, from 65.5% to 79.5%, and a significant increase in the classifier's average confidence for the testing set, from 14.6% to 42.9%, after a few AL iterations.

The helium-peculiar star a Cen exhibits line profile variations of elements such as iron, nitrogen and oxygen in addition to its well-known extreme helium variability. New high S/N, high-resolution spectra are used to perform a quantitative measurement of the abundances of the star and determine the relation of the concentrations of the heavier elements on the surface of the star to the helium concentration and the magnetic field orientation. Doppler images have been created using programs described in earlier papers by Rice and others. An alternative surface abundance mapping code has been used to model the helium line variations after our Doppler imaging of certain individual helium lines produced mediocre results. We confirm the long-known existence of helium-rich and helium-poor hemispheres on a Cen and we measure a difference of more than two orders of magnitude in helium abundance from one side of the star to the other. Helium is overabundant by a factor of about 5 over much of the helium-rich hemisphere. Of particular note is our discovery that the helium-poor hemisphere has a very high abundance of helium-3, approximately equal to the helium-4 abundance. a Cen is therefore a new member of the small group of helium-3 stars and the first well-established magnetic member of the class. For the three metals investigated here, there are two strong concentrations of abundance near the equator consistent with the positive magnetic maximum and two somewhat weaker concentrations of abundance where the helium concentration is centered and roughly where the negative peak of the magnetic field would be found. Another strong concentration is found near the equator and this is not explainable in terms of any simple symmetry with the helium abundance or the apparent magnetic field main polar locations.

The Taiwanese-American Occultation Survey (TAOS) project has collected more than a billion photometric measurements since 2005 January. These sky survey data-covering timescales from a fraction of a second to a few hundred days-are a useful source to study stellar variability. A total of 167 star fields, mostly along the ecliptic plane, have been selected for photometric monitoring with the TAOS telescopes. This paper presents our initial analysis of a search for periodic variable stars from the time-series TAOS data on one particular TAOS field, No. 151 (RA = 17$^{\rm h}30^{\rm m}6\fs$67, Dec = 27\degr17\arcmin 30\arcsec, J2000), which had been observed over 47 epochs in 2005. A total of 81 candidate variables are identified in the 3 square degree field, with magnitudes in the range 8 < R < 16. On the basis of the periodicity and shape of the lightcurves, 29 variables, 15 of which were previously unknown, are classified as RR Lyrae, Cepheid, delta Scuti, SX Phonencis, semi-regular and eclipsing binaries.

The Taiwanese-American Occultation Survey (TAOS) monitors fields of up to ~1000 stars at 5 Hz simultaneously with four small telescopes to detect occultation events from small (~1 km) Kuiper Belt Objects (KBOs). The survey presents a number of challenges, in particular the fact that the occultation events we are searching for are extremely rare and are typically manifested as slight flux drops for only one or two consecutive time series measurements. We have developed a statistical analysis technique to search the multi-telescope data set for simultaneous flux drops which provides a robust false positive rejection and calculation of event significance. In this paper, we describe in detail this statistical technique and its application to the TAOS data set.

We have analyzed the first 3.75 years of data from TAOS, the Taiwanese American Occultation Survey. TAOS monitors bright stars to search for occultations by Kuiper Belt Objects (KBOs). This dataset comprises 5e5 star-hours of multi-telescope photometric data taken at 4 or 5 Hz. No events consistent with KBO occultations were found in this dataset. We compute the number of events expected for the Kuiper Belt formation and evolution models of Pan & Sari (2005), Kenyon & Bromley (2004), Benavidez & Campo Bagatin (2009), and Fraser (2009). A comparison with the upper limits we derive from our data constrains the parameter space of these models. This is the first detailed comparison of models of the KBO size distribution with data from an occultation survey. Our results suggest that the KBO population is comprised of objects with low internal strength and that planetary migration played a role in the shaping of the size distribution.

We analyzed data accumulated during 2005 and 2006 by the Taiwan-American Occultation Survey (TAOS) in order to detect short-period variable stars (periods of <~ 1 hour) such as delta Scuti. TAOS is designed for the detection of stellar occultation by small-size Kuiper Belt Objects (KBOs) and is operating four 50cm telescopes at an effective cadence of 5Hz. The four telescopes simultaneously monitor the same patch of the sky in order to reduce false positives. To detect short-period variables, we used the Fast Fourier Transform algorithm (FFT) inasmuch as the data points in TAOS light-curves are evenly spaced. Using FFT, we found 41 short-period variables with amplitudes smaller than a few hundredths of a magnitude and periods of about an hour, which suggest that they are low-amplitude delta Scuti stars (LADS). The light-curves of TAOS delta Scuti stars are accessible online at the Time Series Center website (http://timemachine.iic.harvard.edu)

We present the results of a search for occultation events by objects at distances between 100 and 1000 AU in lightcurves from the Taiwanese-American Occultation Survey (TAOS). We searched for consecutive, shallow flux reductions in the stellar lightcurves obtained by our survey between 7 February 2005 and 31 December 2006 with a total of $\sim4.5\times10^{9}$ three-telescope simultaneous photometric measurements. No events were detected, allowing us to set upper limits on the number density as a function of size and distance of objects in Sedna-like orbits, using simple models.

The star BD+29 1748 was resolved to be a close binary from its occultation by the asteroid 87 Sylvia on 2006 December 18 UT. Four telescopes were used to observe this event at two sites separated by some 80 km apart. Two flux drops were observed at one site, whereas only one flux drop was detected at the other. From the long-term variation of Sylvia, we inferred the probable shape of the shadow during the occultation, and this in turn constrains the binary parameters: the two components of BD+29 1748 have a projected separation of 0.097" to 0.110" on the sky with a position angle 104 deg to 107 deg. The asteroid was clearly resolved with a size scale ranging from 130 to 290 km, as projected onto the occultation direction. No occultation was detected for either of the two known moonlets of 87 Sylvia.

Results from the first two years of data from the Taiwanese-American Occultation Survey (TAOS) are presented. Stars have been monitored photometrically at 4 Hz or 5 Hz to search for occultations by small (~3 km) Kuiper Belt Objects (KBOs). No statistically significant events were found, allowing us to present an upper bound to the size distribution of KBOs with diameters 0.5 km < D < 28 km.

This paper treats the problem of detecting periodicity in a sequence of photon arrival times, which occurs, for example, in attempting to detect gamma-ray pulsars. A particular focus is on how auxiliary information, typically source intensity, background intensity, and incidence angles and energies associated with each photon arrival should be used to maximize the detection power. We construct a class of likelihood-based tests, score tests, which give rise to event weighting in a principled and natural way, and derive expressions quantifying the power of the tests. These results can be used to compare the efficacies of different weight functions, including cuts in energy and incidence angle. The test is targeted toward a template for the periodic lightcurve, and we quantify how deviation from that template affects the power of detection.

General relativistic kinematics and the cosmological principle alone imply a monotonicity constraint in the Hubble diagram, which we confront to present-day supernova data. We use the running gradient method of statistical inference by Hall & Heckman (2000). We find no significant departure from monotonicity. The method seems well adapted and we recommend its use with future data.

The Taiwan-America Occultation Survey (TAOS) aims to determine the number of small icy bodies in the outer reach of the Solar System by means of stellar occultation. An array of 4 robotic small (D=0.5 m), wide-field (f/1.9) telescopes have been installed at Lulin Observatory in Taiwan to simultaneously monitor some thousand of stars for such rare occultation events. Because a typical occultation event by a TNO a few km across will last for only a fraction of a second, fast photometry is necessary. A special CCD readout scheme has been devised to allow for stellar photometry taken a few times per second. Effective analysis pipelines have been developed to process stellar light curves and to correlate any possible flux changes among all telescopes. A few billion photometric measurements have been collected since the routine survey began in early 2005. Our preliminary result of a very low detection rate suggests a deficit of small TNOs down to a few km size, consistent with the extrapolation of some recent studies of larger (30--100 km) TNOs.

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.