We analyzed the young (2.8-Myr-old) binary system FS Tau A using near-infrared (H-band) high-contrast polarimetry data from Subaru/HiCIAO and sub-millimeter CO (J=2-1) line emission data from ALMA. Both the near-infrared and sub-millimeter observations reveal several clear structures extending to $\sim$240 AU from the stars. Based on these observations at different wavelengths, we report the following discoveries. One arm-like structure detected in the near-infrared band initially extends from the south of the binary with a subsequent turn to the northeast, corresponding to two bar-like structures detected in ALMA observations with an LSRK velocity of 1.19-5.64 km/s. Another feature detected in the near-infrared band extends initially from the north of the binary, relating to an arm-like structure detected in ALMA observations with an LSRK velocity of 8.17-16.43 km/s. From their shapes and velocities, we suggest that these structures can mostly be explained by two streamers that connect the outer circumbinary disk and the central binary components. These discoveries will be helpful for understanding the evolution of streamers and circumstellar disks in young binary systems.

The SR24 multi-star system hosts both circumprimary and circumsecondary disks, which are strongly misaligned from each other. The circumsecondary disk is circumbinary in nature. Interestingly, both disks are interacting, and they possibly rotate in opposite directions. To investigate the nature of this unique twin disk system, we present 0.''1 resolution near-infrared polarized intensity images of the circumstellar structures around SR24, obtained with HiCIAO mounted on the Subaru 8.2 m telescope. Both the circumprimary disk and the circumsecondary disk are resolved and have elongated features. While the position angle of the major axis and radius of the NIR polarization disk around SR24S are 55$^{\circ}$ and 137 au, respectively, those around SR24N are 110$^{\circ}$ and 34 au, respectively. With regard to overall morphology, the circumprimary disk around SR24S shows strong asymmetry, whereas the circumsecondary disk around SR24N shows relatively strong symmetry. Our NIR observations confirm the previous claim that the circumprimary and circumsecondary disks are misaligned from each other. Both the circumprimary and circumsecondary disks show similar structures in $^{12}$CO observations in terms of its size and elongation direction. This consistency is because both NIR and $^{12}$CO are tracing surface layers of the flared disks. As the radius of the polarization disk around SR24N is roughly consistent with the size of the outer Roche lobe, it is natural to interpret the polarization disk around SR24N as a circumbinary disk surrounding the SR24Nb-Nc system.

We present H-band polarized scattered light imagery and JHK high-contrast spectroscopy of the protoplanetary disk around HD 163296 observed with the HiCIAO and SCExAO/CHARIS instruments at Subaru Observatory. The polarimetric imagery resolve a broken ring structure surrounding HD 163296 that peaks at a distance along the major axis of 0.65 (66 AU) and extends out to 0.98 (100 AU) along the major axis. Our 2011 H-band data exhibit clear axisymmetry, with the NW- and SE- side of the disk exhibiting similar intensities. Our data are clearly different than 2016 epoch H-band observations from VLT/SPHERE that found a strong 2.7x asymmetry between the NW- and SE-side of the disk. Collectively, these results indicate the presence of time variable, non-azimuthally symmetric illumination of the outer disk. Based on our 3D-MCRT modeling of contemporaneous IR spectroscopic and H-band polarized intensity imagery of the system, we suggest that while the system could plausibly host an inclined inner disk component, such a component is unlikely to be responsible for producing the observed time-dependent azimuthal variations in the outer scattered light disk of the system. While our SCExAO/CHARIS data are sensitive enough to recover the planet candidate identified from NIRC2 in the thermal IR, we fail to detect an object with a corresponding JHK brightness estimated from the atmospheric models of Baraffe et al. 2003. This suggests that the candidate is either fainter in JHK bands than model predictions, possibly due to extinction from the disk or atmospheric dust/clouds, or that it is an artifact of the dataset/data processing. Our SCExAO/CHARIS data lower the IR mass limits for planets inferred at larger stellocentric separations; however, these ALMA-predicted protoplanet candidates are currently still consistent with direct imaging constraints.

We report an orbital characterization of GJ1108Aab that is a low-mass binary system in pre-main-sequence phase. Via the combination of astrometry using adaptive optics and radial velocity measurements, an eccentric orbital solution of $e$=0.63 is obtained, which might be induced by the Kozai-Lidov mechanism with a widely separated GJ1108B system. Combined with several observed properties, we confirm the system is indeed young. Columba is the most probable moving group, to which the GJ1108A system belongs, although its membership to the group has not been established. If the age of Columba is assumed for GJ1108A, the dynamical masses of both GJ1108Aa and GJ1108Ab ($M_{\rm dynamical,GJ1108Aa}=0.72\pm0.04 M_{\odot}$ and $M_{\rm dynamical,GJ1108Ab}=0.30\pm0.03 M_{\odot}$) are more massive than what an evolutionary model predicts based on the age and luminosities. We consider the discrepancy in mass comparison can attribute to an age uncertainty; the system is likely older than stars in Columba, and effects that are not implemented in classical models such as accretion history and magnetic activity are not preferred to explain the mass discrepancy. We also discuss the performance of the evolutionary model by compiling similar low-mass objects in evolutionary state based on the literature. Consequently, it is suggested that the current model on average reproduces the mass of resolved low-mass binaries without any significant offsets.

We present $H$- and $K_{\rm s}$-bands observations of the LkH$\alpha$ 330 disk with a multi-band detection of the large gap and spiral-like structures. The morphology of the outer disk ($r\sim$$0\farcs3$) at PA=0--45$^\circ$ and PA=180--290$^\circ$ are likely density wave-induced spirals and comparison between our observational results and simulations suggests a planet formation. We have also investigated the azimuthal profiles at the ring and the outer-disk regions as well as radial profiles in the directions of the spiral-like structures and semi-major axis. Azimuthal analysis shows a large variety in wavelength and implies that the disk has non-axisymmetric dust distributions. The radial profiles in the major-axis direction (PA=$271^\circ$) suggest that the outer region ($r\geq0\farcs25$) may be influenced by shadows of the inner region of the disk. The spiral-like directions (PA=10$^\circ$ and 230$^\circ$) show different radial profiles, which suggests that the surfaces of the spiral-like structures are highly flared and/or have different dust properties. Finally, a color-map of the disk shows a lack of an outer eastern region in the $H$-band disk, which may hint the presence of an inner object that casts a directional shadow onto the disk.

We conducted high-contrast polarimetry observations of T Tau in the H-band, using the HiCIAO instrument mounted on the Subaru Telescope, revealing structures as near as 0.$\arcsec$1 from the stars T Tau N and T Tau S. The whole T Tau system is found to be surrounded by nebula-like envelopes, and several outflow-related structures are detected in these envelopes. We analyzed the detailed polarization patterns of the circumstellar structures near each component of this triple young star system and determined constraints on the circumstellar disks and outflow structures. We suggest that the nearly face-on circumstellar disk of T Tau N is no larger than 0.$\arcsec$8, or 117 AU, in the northwest, based on the existence of a hole in this direction, and no larger than 0.$\arcsec$27, or 40 AU, in the south. A new structure "N5" extends to about 0.$\arcsec$42, or 59 AU, on the southwest of the star, believed to be part of the disk. We suggest that T Tau S is surrounded by a highly inclined circumbinary disk with a radius of about 0.$\arcsec$3, or 44 AU, with a position angle of about 30$^\circ$, that is misaligned with the orbit of the T Tau S binary. After analyzing the positions and polarization vector patterns of the outflow-related structures, we suggest that T Tau S should trigger the well-known E-W outflow, and is also likely to be responsible for a southwest precessing outflow "coil" and a possible south outflow.

Hyper Suprime-Cam (HSC) is a wide-field imaging camera on the prime focus of the 8.2m Subaru telescope on the summit of Maunakea in Hawaii. A team of scientists from Japan, Taiwan and Princeton University is using HSC to carry out a 300-night multi-band imaging survey of the high-latitude sky. The survey includes three layers: the Wide layer will cover 1400 deg$^2$ in five broad bands ($grizy$), with a $5\,\sigma$ point-source depth of $r \approx 26$. The Deep layer covers a total of 26~deg$^2$ in four fields, going roughly a magnitude fainter, while the UltraDeep layer goes almost a magnitude fainter still in two pointings of HSC (a total of 3.5 deg$^2$). Here we describe the instrument, the science goals of the survey, and the survey strategy and data processing. This paper serves as an introduction to a special issue of the Publications of the Astronomical Society of Japan, which includes a large number of technical and scientific papers describing results from the early phases of this survey.

Large exoplanet surveys have successfully detected thousands of exoplanets to-date. Utilizing these detections and non-detections to constrain our understanding of the formation and evolution of planetary systems also requires a detailed understanding of the basic properties of their host stars. We have determined the basic stellar properties of F, K, and G stars in the Strategic Exploration of Exoplanets and Disks with Subaru (SEEDS) survey from echelle spectra taken at the Apache Point Observatory's 3.5m telescope. Using ROBOSPECT to extract line equivalent widths and TGVIT to calculate the fundamental parameters, we have computed Teff, log(g), vt, [Fe/H], chromospheric activity, and the age for our sample. Our methodology was calibrated against previously published results for a portion of our sample. The distribution of [Fe/H] in our sample is consistent with that typical of the Solar neighborhood. Additionally, we find the ages of most of our sample are $< 500 Myrs$, but note that we cannot determine robust ages from significantly older stars via chromospheric activity age indicators. The future meta-analysis of the frequency of wide stellar and sub-stellar companions imaged via the SEEDS survey will utilize our results to constrain the occurrence of detected co-moving companions with the properties of their host stars.

The Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) is a three-layered imaging survey aimed at addressing some of the most outstanding questions in astronomy today, including the nature of dark matter and dark energy. The survey has been awarded 300 nights of observing time at the Subaru Telescope and it started in March 2014. This paper presents the first public data release of HSC-SSP. This release includes data taken in the first 1.7 years of observations (61.5 nights) and each of the Wide, Deep, and UltraDeep layers covers about 108, 26, and 4 square degrees down to depths of i~26.4, ~26.5, and ~27.0 mag, respectively (5sigma for point sources). All the layers are observed in five broad bands (grizy), and the Deep and UltraDeep layers are observed in narrow bands as well. We achieve an impressive image quality of 0.6 arcsec in the i-band in the Wide layer. We show that we achieve 1-2 per cent PSF photometry (rms) both internally and externally (against Pan-STARRS1), and ~10 mas and 40 mas internal and external astrometric accuracy, respectively. Both the calibrated images and catalogs are made available to the community through dedicated user interfaces and database servers. In addition to the pipeline products, we also provide value-added products such as photometric redshifts and a collection of public spectroscopic redshifts. Detailed descriptions of all the data can be found online. The data release website is https://hsc-release.mtk.nao.ac.jp/.

With dozens of Jovian and super-Jovian exoplanets known to orbit their host stars in or near the stellar habitable zones, it has recently been suggested that moons the size of Mars could offer abundant surface habitats beyond the solar system. Several searches for such exomoons are now underway, and the exquisite astronomical data quality of upcoming space missions and ground-based extremely large telescopes could make the detection and characterization of exomoons possible in the near future. Here we explore the effects of tidal heating on the potential of Mars- to Earth-sized satellites to host liquid surface water, and we compare the tidal heating rates predicted by tidal equilibrium model and a viscoelastic model. In addition to tidal heating, we consider stellar radiation, planetary illumination and thermal heat from the planet. However, the effects of a possible moon atmosphere are neglected. We map the circumplanetary habitable zone for different stellar distances in specific star-planet-satellite configurations, and determine those regions where tidal heating dominates over stellar radiation. We find that the `thermostat effect' of the viscoelastic model is significant not just at large distances from the star, but also in the stellar habitable zone, where stellar radiation is prevalent. We also find that tidal heating of Mars-sized moons with eccentricities between 0.001 and 0.01 is the dominant energy source beyond 3--5 AU from a Sun-like star and beyond 0.4--0.6 AU from an M3 dwarf star. The latter would be easier to detect (if they exist), but their orbital stability might be under jeopardy due to the gravitational perturbations from the star.

We present high-contrast observations of 68 young stellar objects (YSOs) explored as part of the SEEDS survey on the Subaru telescope. Our targets are very young ($<$10 Myr) stars, which often harbor protoplanetary disks where planets may be forming. We achieve a typical contrast of $\sim$$10^{-4}$--$10^{-5.5}$ at an angular distance of 1\arcsec\ from the central star, corresponding to typical mass sensitivities (assuming hot-start evolutionary models) of $\sim$10 ${\rm M_J}$ at 70 AU and $\sim$6 ${\rm M_J}$ at 140 AU. We detected a new stellar companion to HIP 79462 and confirmed the substellar objects GQ Lup b and ROXs 42B b. An additional six companion candidates await follow-up observations to check for common proper motion. Our SEEDS YSO observations probe the population of planets and brown dwarfs at the very youngest ages; these may be compared to the results of surveys targeting somewhat older stars. Our sample and the associated observational results will help enable detailed statistical analyses of giant planet formation.

By performing non-masked polarization imaging with Subaru/HiCIAO, polarized scattered light from the inner region of the disk around the GG Tau A system was successfully detected in the $H$ band with a spatial resolution of approximately 0.07$\arcsec$, revealing the complicated inner disk structures around this young binary. This paper reports the observation of an arc-like structure to the north of GG Tau Ab and part of a circumstellar structure that is noticeable around GG Tau Aa extending to a distance of approximately 28 AU from the primary star. The speckle noise around GG Tau Ab constrains its disk radius to <13 AU. Based on the size of the circumbinary ring and the circumstellar disk around GG Tau Aa, the semi-major axis of the binary's orbit is likely to be 62 AU. A comparison of the present observations with previous ALMA and near-infrared (NIR) H$_2$ emission observations suggests that the north arc could be part of a large streamer flowing from the circumbinary ring to sustain the circumstellar disks. According to the previous studies, the circumstellar disk around GG Tau Aa has enough mass and can sustain itself for a duration sufficient for planet formation; thus, our study indicates that planets can form within close (separation $\lesssim$ 100 AU) young binary systems.

We present high-contrast H-band polarized intensity (PI) images of the transitional disk around the young solar-like star GM Aur. The near-infrared direct imaging of the disk was derived by polarimetric differential imaging using the Subaru 8.2-m Telescope and HiCIAO. An angular resolution and an inner working angle of 0."07 and r~0."05, respectively, were obtained. We clearly resolved a large inner cavity, with a measured radius of 18+/-2 au, which is smaller than that of a submillimeter interferometric image (28 au). This discrepancy in the cavity radii at near-infrared and submillimeter wavelengths may be caused by a 3-4M_Jup planet about 20 au away from the star, near the edge of the cavity. The presence of a near-infrared inner is a strong constraint on hypotheses for inner cavity formation in a transitional disk. A dust filtration mechanism has been proposed to explain the large cavity in the submillimeter image, but our results suggest that this mechanism must be combined with an additional process. We found that the PI slope of the outer disk is significantly different from the intensity slope obtained from HST/NICMOS, and this difference may indicate the grain growth process in the disk.

We present high signal-to-noise ratio, precise $YJH$ photometry and $Y$ band (\gpiwave~$\mu$m) spectroscopy of HD 1160 B, a young substellar companion discovered from the Gemini NICI Planet Finding Campaign, using the Subaru Coronagraphic Extreme Adaptive Optics instrument and the Gemini Planet Imager. HD 1160 B has typical mid-M dwarf-like infrared colors and a spectral type of M5.5$^{+1.0}_{-0.5}$, where the blue edge of our $Y$ band spectrum rules out earlier spectral types. Atmospheric modeling suggests HD 1160 B having an effective temperature of 3000--3100 $K$, a surface gravity of log $g$ = 4--4.5, a radius of~\bestfitradius~$R_{\rm J}$, and a luminosity of log $L$/$L_{\odot} = -2.76 \pm 0.05$. Neither the primary's Hertzspring-Russell diagram position nor atmospheric modeling of HD 1160 B show evidence for a sub-solar metallicity. The interpretation of the HD 1160 B depends on which stellar system components are used to estimate an age. Considering HD 1160 A, B and C jointly, we derive an age of 80--125 Myr, implying that HD 1160 B straddles the hydrogen-burning limit (70--90 $M_{\rm J}$). If we consider HD 1160 A alone, younger ages (20--125 Myr) and a brown dwarf-like mass (35--90 $M_{\rm J}$) are possible. Interferometric measurements of the primary, a precise GAIA parallax, and moderate resolution spectroscopy can better constrain the system's age and how HD 1160 B fits within the context of (sub)stellar evolution.

We present H-band (1.6 {\mu}m) scattered light observations of the transitional disk RX J1615.3-3255, located in the ~1 Myr old Lupus association. From a polarized intensity image, taken with the HiCIAO instrument of the Subaru Telescope, we deduce the position angle and the inclination angle of the disk. The disk is found to extend out to 68 $\pm$ 12 AU in scattered light and no clear structure is observed. Our inner working angle of 24 AU does not allow us to detect a central decrease in intensity similar to that seen at 30 AU in the 880 {\mu}m continuum observations. We compare the observations with multiple disk models based on the Spectral Energy Distribution (SED) and submm interferometry and find that an inner rim of the outer disk at 30 AU containing small silicate grains produces a polarized intensity signal which is an order of magnitude larger than observed. We show that a model in which the small dust grains extend smoothly into the cavity found for large grains is closer to the actual H-band observations. A comparison of models with different dust size distributions suggests that the dust in the disk might have undergone significant processing compared to the interstellar medium.

We report the direct imaging detection of a low-mass companion to a young, moderately active star V450 And, that was previously identified with the radial velocity method. The companion was found in high-contrast images obtained with the Subaru Telescope equipped with the HiCIAO camera and AO188 adaptive optics system. From the public ELODIE and SOPHIE archives we extracted available high-resolution spectra and radial velocity (RV) measurements, along with RVs from the Lick planet search program. We combined our multi-epoch astrometry with these archival, partially unpublished RVs, and found that the companion is a low-mass star, not a brown dwarf, as previously suggested. We found the best-fitting dynamical masses to be $m_1=1.141_{-0.091}^{+0.037}$ and $m_2=0.279^{+0.023}_{-0.020}$ M$_\odot$. We also performed spectral analysis of the SOPHIE spectra with the iSpec code. The Hipparcos time-series photometry shows a periodicity of $P=5.743$ d, which is also seen in SOPHIE spectra as an RV modulation of the star A. We interpret it as being caused by spots on the stellar surface, and the star to be rotating with the given period. From the rotation and level of activity, we found that the system is $380^{+220}_{-100}$ Myr old, consistent with an isochrone analysis ($220^{+2120}_{-90}$ Myr). This work may serve as a test case for future studies of low-mass stars, brown dwarfs and exoplanets by combination of RV and direct imaging data.

We present high-contrast H-band polarized intensity images of the transitional disk around the young solar-like star LkCa 15. By utilizing Subaru/HiCIAO for polarimetric differential imaging, both the angular resolution and the inner working angle reach 0.07" and r=0.1", respectively. We obtained a clearly resolved gap (width <~ 27 AU) at ~ 48 AU from the central star. This gap is consistent with images reported in previous studies. We also confirmed the existence of a bright inner disk with a misaligned position angle of 13+/-4 degree with respect to that of the outer disk, i.e., the inner disk is possibly warped. The large gap and the warped inner disk both point to the existence of a multiple planetary system with a mass of <~1Mjup.

We find a new substellar companion to the Pleiades member star, Pleiades HII 3441, using the Subaru telescope with adaptive optics. The discovery is made as part of the high-contrast imaging survey to search for planetary-mass and substellar companions in the Pleiades and young moving groups. The companion has a projected separation of 0".49 +/- 0".02 (66 +/- 2 AU) and a mass of 68 +/- 5 M_J based on three observations in the J-, H-, and K_S-band. The spectral type is estimated to be M7 (~2700 K), and thus no methane absorption is detected in the H band. Our Pleiades observations result in the detection of two substellar companions including one previously reported among 20 observed Pleiades stars, and indicate that the fraction of substellar companions in the Pleiades is about 10.0 +26.1/-8.8 %. This is consistent with multiplicity studies of both the Pleiades stars and other open clusters.

Dust trapping accelerates the coagulation of dust particles, and thus it represents an initial step toward the formation of planetesimals. We report $H$-band (1.6 um) linear polarimetric observations and 0.87 mm interferometric continuum observations toward a transitional disk around LkHa 330. As results, a pair of spiral arms were detected in the $H$-band emission and an asymmetric (potentially arm-like) structure was detected in the 0.87 mm continuum emission. We discuss the origin of the spiral arm and the asymmetric structure, and suggest that a massive unseen planet is the most plausible explanation. The possibility of dust trapping and grain growth causing the asymmetric structure was also investigated through the opacity index (beta) by plotting the observed SED slope between 0.87 mm from our SMA observation and 1.3 mm from literature. The results imply that grains are indistinguishable from ISM-like dust in the east side ($beta = 2.0 pm 0.5$), but much smaller in the west side $beta = 0.7^{+0.5}_{-0.4}$, indicating differential dust size distribution between the two sides of the disk. Combining the results of near-infrared and submillimeter observations, we conjecture that the spiral arms exist at the upper surface and an asymmetric structure resides in the disk interior. Future observations at centimeter wavelengths and differential polarization imaging in other bands (Y to K) with extreme AO imagers are required to understand how large dust grains form and to further explore the dust distribution in the disk.

We present new analysis of multi-epoch, H-band, scattered light images of the AB Aur system. We used a Monte Carlo, radiative transfer code to simultaneously model the system's SED and H-band polarized intensity imagery. We find that a disk-dominated model, as opposed to one that is envelope dominated, can plausibly reproduce AB Aur's SED and near-IR imagery. This is consistent with previous modeling attempts presented in the literature and supports the idea that at least a subset of AB Aur's spirals originate within the disk. In light of this, we also analyzed the movement of spiral structures in multi-epoch H-band total light and polarized intensity imagery of the disk. We detect no significant rotation or change in spatial location of the spiral structures in these data, which span a 5.8 year baseline. If such structures are caused by disk-planet interactions, the lack of observed rotation constrains the location of the orbit of planetary perturbers to be >47 AU.

We present high-contrast angular differential imaging (ADI) observations of the debris disk around HD 32297 in H-band, as well as the first polarimetric images for this system in polarized differential imaging (PDI) mode with Subaru/HICIAO. In ADI, we detect the nearly edge-on disk at >5sigma levels from ~0.45 arcsec to ~1.7 arcsec (50-192 AU) from the star and recover the spine deviation from the midplane already found in previous works. We also find for the first time imaging and surface brightness (SB) indications for the presence of a gapped structure on both sides of the disk at distances of ~0.75 arcsec (NE side) and ~0.65 arcsec (SW side). Global forward-modelling work delivers a best-fit model disk and well-fitting parameter intervals that essentially match previous results, with high-forward scattering grains and a ring located at 110 AU. However, this single ring model cannot account for the gapped structure seen in our SB profiles. We create simple double ring models and achieve a satisfactory fit with two rings located at 60 and 95 AU, respectively, low-forward scattering grains and very sharp inner slopes. In polarized light we retrieve the disk extending from ~0.25-1.6 arcsec, although the central region is quite noisy and high S/N are only found in the range ~0.75-1.2 arcsec. The disk is polarized in the azimuthal direction, as expected, and the departure from the midplane is also clearly observed. Evidence for a gapped scenario is not found in the PDI data. We obtain a linear polarization degree of the grains that increases from ~10% at 0.55 arcsec to ~25% at 1.6 arcsec. The maximum is found at scattering angles of ~90degrees, either from the main components of the disk or from dust grains blown out to larger radii.

We present the first near-infrared scattered-light detection of the transitional disk around V1247 Ori, which was obtained using high-resolution polarimetric differential imaging observations with Subaru/HiCIAO. Our imaging in the H band reveals the disk morphology at separations of ~0.14"-0.86" (54-330 au) from the central star. The polarized intensity (PI) image shows a remarkable arc-like structure toward the southeast of the star, whereas the fainter northwest region does not exhibit any notable features. The shape of the arm is consistent with an arc of 0.28" $\pm$ 0.09" in radius (108 au from the star), although the possibility of a spiral arm with a small pitch angle cannot be excluded. V1247 Ori features an exceptionally large azimuthal contrast in scattered, polarized light; the radial peak of the southeastern arc is about three times brighter than the northwestern disk measured at the same distance from the star. Combined with the previous indication of an inhomogeneous density distribution in the gap at $\lesssim$46 au, the notable asymmetry in the outer disk suggests the presence of unseen companions and/or planet-forming processes ongoing in the arc.

A radial velocity (RV) survey for intermediate-mass giants has been operated for over a decade at Okayama Astrophysical Observatory (OAO). The OAO survey has revealed that some giants show long-term linear RV accelerations (RV trends), indicating the presence of outer companions. Direct imaging observations can help clarify what objects generate these RV trends. We present the results of high-contrast imaging observations or six intermediate-mass giants with long-term RV trends using the Subaru Telescope and HiCIAO camera. We detected co-moving companions to $\gamma$ Hya B ($0.61^{+0.12}_{-0.14} M_\odot$), HD 5608 B ($0.10 \pm 0.01 M_\odot$), and HD 109272 B ($0.28 \pm 0.06 M_\odot$). For the remaining targets($\iota$ Dra, 18 Del, and HD 14067) we exclude companions more massive than 30-60 $M_\mathrm{Jup}$ at projected separations of 1arcsec-7arcsec. We examine whether these directly imaged companions or unidentified long-period companions can account for the RV trends observed around the six giants. We find that the Kozai mechanism can explain the high eccentricity of the inner planets $\iota$ Dra b, HD 5608 b, and HD 14067 b.

A high angular resolution near-infrared polarized-intensity image of the GG Tau A binary system was obtained with the Subaru Telescope. The image shows the circumbinary disk scattering the light from the central binary. The azimuthal profile of the polarized intensity of the circumbinary disk is roughly reproduced by a simple disk model with the Henyey-Greenstein function and the Rayleigh function, indicating small dust grains at the surface of the disk. Combined with a previous observation of the circumbinary disk, our image indicates that the gap structure in the circumbinary disk orbits anti-clockwise, while material in the disk orbit clockwise. We propose a shadow of material located between the central binary and the circumbinary disk. The separations and position angles of the stellar components of the binary in the past 20 years are consistent with the binary orbit with a = 33.4 AU and e = 0.34.

We present the first spatially resolved polarized scattered light H-band detection of the DoAr 28 transitional disk. Our two epochs of imagery detect the scattered light disk from our effective inner working angle of 0.10" (13 AU) out to 0.50" (65 AU). This inner working angle is interior to the location of the system's gap inferred by previous studies using SED modeling (15 AU). We detected a candidate point source companion 1.08" northwest of the system; however, our second epoch of imagery strongly suggests that this object is a background star. We constructed a grid of Monte Carlo Radiative Transfer models of the system, and our best fit models utilize a modestly inclined (50 deg), 0.01 Msun disk that has a partially depleted inner gap from the dust sublimation radius out to ~8 AU. Subtracting this best fit, axi-symmetric model from our polarized intensity data reveals evidence for two small asymmetries in the disk, which could be attributable to variety of mechanisms.

We present the photometric properties of a sample of infrared (IR) bright dust obscured galaxies (DOGs). Combining wide and deep optical images obtained with the Hyper Suprime-Cam (HSC) on the Subaru Telescope and all-sky mid-IR (MIR) images taken with Wide-Field Infrared Survey Explorer (WISE), we discovered 48 DOGs with $i - K_\mathrm{s} > 1.2$ and $i - [22] > 7.0$, where $i$, $K_\mathrm{s}$, and [22] represent AB magnitude in the $i$-band, $K_\mathrm{s}$-band, and 22 $\mu$m, respectively, in the GAMA 14hr field ($\sim$ 9 deg$^2$). Among these objects, 31 ($\sim$ 65 %) show power-law spectral energy distributions (SEDs) in the near-IR (NIR) and MIR regime, while the remainder show a NIR bump in their SEDs. Assuming that the redshift distribution for our DOGs sample is Gaussian, with mean and sigma $z$ = 1.99 $\pm$ 0.45, we calculated their total IR luminosity using an empirical relation between 22 $\mu$m luminosity and total IR luminosity. The average value of the total IR luminosity is (3.5 $\pm$ 1.1) $\times$ $10^{13}$ L$_{\odot}$, which classifies them as hyper-luminous infrared galaxies (HyLIRGs). We also derived the total IR luminosity function (LF) and IR luminosity density (LD) for a flux-limited subsample of 18 DOGs with 22 $\mu$m flux greater than 3.0 mJy and with $i$-band magnitude brighter than 24 AB magnitude. The derived space density for this subsample is log $\phi$ = -6.59 $\pm$ 0.11 [Mpc$^{-3}$]. The IR LF for DOGs including data obtained from the literature is well fitted by a double-power law. The derived lower limit for the IR LD for our sample is $\rho_{\mathrm{IR}}$ $\sim$ 3.8 $\times$ 10$^7$ [L$_{\odot}$ Mpc$^{-3}$] and its contributions to the total IR LD, IR LD of all ultra-luminous infrared galaxies (ULIRGs), and that of all DOGs are $>$ 3 %, $>$ 9 %, and $>$ 15 %, respectively.

Coronagraphic imagery of the circumstellar disk around HD 169142 in H-band polarized intensity (PI) with Subaru/HiCIAO is presented. The emission scattered by dust particles at the disk surface in 0.2" <= r <= 1.2", or 29 <= r <= 174 AU, is successfully detected. The azimuthally-averaged radial profile of the PI shows a double power-law distribution, in which the PIs in r=29-52 AU and r=81.2-145 AU respectively show r^{-3}-dependence. These two power-law regions are connected smoothly with a transition zone (TZ), exhibiting an apparent gap in r=40-70 AU. The PI in the inner power-law region shows a deep minimum whose location seems to coincide with the point source at \lambda = 7 mm. This can be regarded as another sign of a protoplanet in TZ. The observed radial profile of the PI is reproduced by a minimally flaring disk with an irregular surface density distribution or with an irregular temperature distribution or with the combination of both. The depletion factor of surface density in the inner power-law region (r< 50 AU) is derived to be <= 0.16 from a simple model calculation. The obtained PI image also shows small scale asymmetries in the outer power-law region. Possible origins for these asymmetries include corrugation of the scattering surface in the outer region, and shadowing effect by a puffed up structure in the inner power-law region.

We present new high-resolution ($\sim$0\farcs09) $H$-band imaging observations of the circumstellar disk around the T Tauri star SU Aur. Our observations with Subaru-HiCIAO have revealed the presence of scattered light as close as 0\farcs15 ($\sim$20 AU) to the star. Within our image, we identify bright emission associated with a disk with a minimum radius of $\sim$90 AU, an inclination of $\sim$35\degr from the plane of the sky, and an approximate P.A. of 15\degr for the major axis. We find a brightness asymmetry between the northern and southern sides of the disk due to a non-axisymmetric disk structure. We also identify a pair of asymmetric tail structures extending east and west from the disk. The western tail extends at least 2\farcs5 (350 AU) from the star, and is probably associated with a reflection nebula previously observed at optical and near-IR wavelengths. The eastern tail extends at least 1\arcsec (140 AU) at the present signal-to-noise. These tails are likely due to an encounter with an unseen brown dwarf, but our results do not exclude the explanation that these tails are outflow cavities or jets.

We present a new Subaru/HiCIAO high-contrast H-band polarized intensity (PI) image of a nearby transitional disk associated with TW Hydrae. The scattered light from the disk was detected from 0.2" to 1.5" (11 - 81 AU) and the PI image shows a clear axisymmetric depression in polarized intensity at ~ 0.4" (~ 20 AU) from the central star, similar to the ~ 80 AU gap previously reported from HST images. Azimuthal polarized intensity profile also shows the disk beyond 0.2" is almost axisymmetric. We discuss two possible scenarios explaining the origin of the polarized intensity depression: 1) a gap structure may exist at ~ 20 AU from the central star because of shallow slope seen in the polarized intensity profile, and 2) grain growth may be occurring in the inner region of the disk. Multi-band observations at NIR and millimeter/sub-millimeter wavelengths play a complementary role in investigating dust opacity and may help reveal the origin of the gap more precisely.

Tidal heating of exomoons may play a key role in their habitability, since the elevated temperature can melt the ice on the body even without significant solar radiation. The possibility of life is intensely studied on Solar System moons such as Europa or Enceladus, where the surface ice layer covers tidally heated water ocean. Tidal forces may be even stronger in extrasolar systems, depending on the properties of the moon and its orbit. For studying the tidally heated surface temperature of exomoons, we used a viscoelastic model for the first time. This model is more realistic than the widely used, so-called fixed Q models, because it takes into account the temperature dependency of the tidal heat flux, and the melting of the inner material. With the use of this model we introduced the circumplanetary Tidal Temperate Zone (TTZ), that strongly depends on the orbital period of the moon, and less on its radius. We compared the results with the fixed Q model and investigated the statistical volume of the TTZ using both models. We have found that the viscoelastic model predicts 2.8 times more exomoons in the TTZ with orbital periods between 0.1 and 3.5 days than the fixed Q model for plausible distributions of physical and orbital parameters. The viscoelastic model gives more promising results in terms of habitability, because the inner melting of the body moderates the surface temperature, acting like a thermostat.

The known population of exoplanets exhibits a much wider range of orbital eccentricities than Solar System planets and has a much higher average eccentricity. These facts have been widely interpreted to indicate that the Solar System is an atypical member of the overall population of planetary systems. We report here on a strong anti-correlation of orbital eccentricity with multiplicity (number of planets in the system) among catalogued radial velocity (RV) systems. The mean, median and rough distribution of eccentricities of Solar System planets fits an extrapolation of this anti-correlation to the eight planet case rather precisely despite the fact that no more than 2 Solar System planets would be detectable with RV data comparable to that in the exoplanet sample. Moreover, even if regarded as a single or double planetary system, the Solar System lies in a reasonably heavily populated region of eccentricity-multiplicity space. Thus, the Solar System is not anomalous among known exoplanetary systems with respect to eccentricities when its multiplicity is taken into account. Specifically, as the multiplicity of a system increases the eccentricity decreases roughly as a power law of index -1.20. A simple and plausible but ad hoc and model-dependent interpretation of this relationship implies that approximately 80% of the one planet and 25% of the two planet systems in our sample have additional, as yet undiscovered, members but that systems of higher observed multiplicity are largely complete (i.e., relatively rarely contain additional undiscovered planets). If low eccentricities indeed favor high multiplicities, habitability may be more common in systems with a larger number of planets.

The formation scenario of a gapped disk, i.e., transitional disk, and its asymmetry is still under debate. Proposed scenarios such as disk-planet interaction, photoevaporation, grain growth, anticyclonic vortex, eccentricity, and their combinations would result in different radial distributions of the gas and the small (sub-$\mu$m size) and large (millimeter size) dust grains as well as asymmetric structures in a disk. Optical/near-infrared (NIR) imaging observations and (sub-)millimeter interferometry can trace small and large dust grains, respectively; therefore multi-wavelength observations could help elucidate the origin of complicated structures of a disk. Here we report SMA observations of the dust continuum at 1.3~mm and $^{12}$CO~$J=2\rightarrow1$ line emission of the pre-transitional protoplanetary disk around the solar-mass star PDS~70. PDS~70, a weak-lined T Tauri star, exhibits a gap in the scattered light from its disk with a radius of $\sim$65~AU at NIR wavelengths. However, we found a larger gap in the disk with a radius of $\sim$80~AU at 1.3~mm. Emission from all three disk components (the gas and the small and large dust grains) in images exhibits a deficit in brightness in the central region of the disk, in particular, the dust-disk in small and large dust grains has asymmetric brightness. The contrast ratio of the flux density in the dust continuum between the peak position to the opposite side of the disk reaches 1.4. We suggest the asymmetries and different gap-radii of the disk around PDS~70 are potentially formed by several (unseen) accreting planets inducing dust filtration.

We present the first resolved near infrared imagery of the transition disk Oph IRS 48 (WLY 2-48), which was recently observed with ALMA to have a strongly asymmetric sub-millimeter flux distribution. H-band polarized intensity images show a $\sim$60AU radius scattered light cavity with two pronounced arcs of emission, one from Northeast to Southeast and one smaller, fainter and more distant arc in the Northwest. K-band scattered light imagery reveals a similar morphology, but with a clear third arc along the Southwestern rim of the disk cavity. This arc meets the Northwestern arc at nearly a right angle, revealing the presence of a spiral arm or local surface brightness deficit in the disk, and explaining the East-West brightness asymmetry in the H-band data. We also present 0.8-5.4$\mu$m IRTF SpeX spectra of this object, which allow us to constrain the spectral class to A0$\pm$1 and measure a low mass accretion rate of 10$^{-8.5}$M$_{\odot}$/yr, both consistent with previous estimates. We investigate a variety of reddening laws in order to fit the mutliwavelength SED of Oph IRS 48 and find a best fit consistent with a younger, higher luminosity star than previous estimates.

We compared the number of faint stars detected in deep survey fields with the current stellar distribution model of the Galaxy and found that the detected number in the H band is significantly smaller than the predicted number. This indicates that M-dwarfs, the major component, are fewer in the halo and the thick disk. We used archived data of several surveys in both the north and south field of GOODS (Great Observatories Origins Deep Survey), MODS in GOODS-N, and ERS and CANDELS in GOODS-S. The number density of M-dwarfs in the halo has to be 20+/-13% relative to that in the solar vicinity, in order for the detected number of stars fainter than 20.5 mag in the H band to match with the predicted value from the model. In the thick disk, the number density of M-dwarfs must be reduced (52+/-13%) or the scale height must be decreased (~600 pc). Alternatively, overall fractions of the halo and thick disks can be significantly reduced to achieve the same effect, because our sample mainly consists of faint M-dwarfs. Our results imply that the M-dwarf population in regions distant from the Galactic plane is significantly smaller than previously thought. We then discussed the implications this has on the suitability of the model predictions for the prediction of non-companion faint stars in direct imaging extrasolar planet surveys by using the best-fit number densities.

We conduct a statistical analysis of a combined sample of direct imaging data, totalling nearly 250 stars. The stars cover a wide range of ages and spectral types, and include five detections ($\kappa$ And b, two $\sim$60 M$_{\rm J}$ brown dwarf companions in the Pleiades, PZ Tel B, and CD$-$35 2722B). For some analyses we add a currently unpublished set of SEEDS observations, including the detections GJ 504b and GJ 758B. We conduct a uniform, Bayesian analysis of all stellar ages using both membership in a kinematic moving group and activity/rotation age indicators. We then present a new statistical method for computing the likelihood of a substellar distribution function. By performing most of the integrals analytically, we achieve an enormous speedup over brute-force Monte Carlo. We use this method to place upper limits on the maximum semimajor axis of the distribution function derived from radial-velocity planets, finding model-dependent values of $\sim$30--100 AU. Finally, we model the entire substellar sample, from massive brown dwarfs to a theoretically motivated cutoff at $\sim$5 M$_{\rm Jup}$, with a single power law distribution. We find that $p(M, a) \propto M^{-0.65\pm0.60} a^{-0.85\pm0.39}$ (1$\sigma$ errors) provides an adequate fit to our data, with 1.0--3.1\% (68\% confidence) of stars hosting 5--70 $M_{\rm Jup}$ companions between 10 and 100 AU. This suggests that many of the directly imaged exoplanets known, including most (if not all) of the low-mass companions in our sample, formed by fragmentation in a cloud or disk, and represent the low-mass tail of the brown dwarfs.

We describe the expected scientific capabilities of CHARIS, a high-contrast integral-field spectrograph (IFS) currently under construction for the Subaru telescope. CHARIS is part of a new generation of instruments, enabled by extreme adaptive optics (AO) systems (including SCExAO at Subaru), that promise greatly improved contrasts at small angular separation thanks to their ability to use spectral information to distinguish planets from quasistatic speckles in the stellar point-spread function (PSF). CHARIS is similar in concept to GPI and SPHERE, on Gemini South and the Very Large Telescope, respectively, but will be unique in its ability to simultaneously cover the entire near-infrared $J$, $H$, and $K$ bands with a low-resolution mode. This extraordinarily broad wavelength coverage will enable spectral differential imaging down to angular separations of a few $\lambda/D$, corresponding to $\sim$$0.\!\!''1$. SCExAO will also offer contrast approaching $10^{-5}$ at similar separations, $\sim$$0.\!\!''1$--$0.\!\!''2$. The discovery yield of a CHARIS survey will depend on the exoplanet distribution function at around 10 AU. If the distribution of planets discovered by radial velocity surveys extends unchanged to $\sim$20 AU, observations of $\sim$200 mostly young, nearby stars targeted by existing high-contrast instruments might find $\sim$1--3 planets. Carefully optimizing the target sample could improve this yield by a factor of a few, while an upturn in frequency at a few AU could also increase the number of detections. CHARIS, with a higher spectral resolution mode of $R \sim 75$, will also be among the best instruments to characterize planets and brown dwarfs like HR 8799 cde and $\kappa$ And b.

We present a new method of analysis for determining the surface geometry of five protoplanetary disks observed with near-infrared imaging polarimetry using Subaru-HiCIAO. Using as inputs the observed distribution of polarized intensity (PI), disk inclination, assumed properties for dust scattering, and other reasonable approximations, we calculate a differential equation to derive the surface geometry. This equation is numerically integrated along the distance from the star at a given position angle. We show that, using these approximations, the local maxima in the PI distribution of spiral arms (SAO 206462, MWC 758) and rings (2MASS J16042165-2130284, PDS 70) are associated with local concave-up structures on the disk surface. We also show that the observed presence of an inner gap in scattered light still allows the possibility of a disk surface that is parallel to the light path from the star, or a disk that is shadowed by structures in the inner radii. Our analysis for rings does not show the presence of a vertical inner wall as often assumed in studies of disks with an inner gap. Finally, we summarize the implications of spiral and ring structures as potential signatures of ongoing planet formation.

The diversity and quantity of moons in the Solar System suggest a manifold population of natural satellites exist around extrasolar planets. Of peculiar interest from an astrobiological perspective, the number of sizable moons in the stellar habitable zones may outnumber planets in these circumstellar regions. With technological and theoretical methods now allowing for the detection of sub-Earth-sized extrasolar planets, the first detection of an extrasolar moon appears feasible. In this review, we summarize formation channels of massive exomoons that are potentially detectable with current or near-future instruments. We discuss the orbital effects that govern exomoon evolution, we present a framework to characterize an exomoon's stellar plus planetary illumination as well as its tidal heating, and we address the techniques that have been proposed to search for exomoons. Most notably, we show that natural satellites in the range of 0.1 - 0.5 Earth mass (i) are potentially habitable, (ii) can form within the circumplanetary debris and gas disk or via capture from a binary, and (iii) are detectable with current technology.

We present results from the first three years of observations of moving group targets in the SEEDS high-contrast imaging survey of exoplanets and disks using the Subaru telescope. We achieve typical contrasts of ~10^5 at 1" and ~10^6 beyond 2" around 63 proposed members of nearby kinematic moving groups. We review each of the kinematic associations to which our targets belong, concluding that five, \beta Pictoris (~20 Myr), AB Doradus (~100 Myr), Columba (~30 Myr), Tucana-Horogium (~30 Myr), TW Hydrae (~10 Myr), are sufficiently well-defined to constrain the ages of individual targets. Somewhat less than half of our targets are high-probability members of one of these moving groups. For all of our targets, we combine proposed moving group membership with other age indicators where available, including Ca II HK emission, X-ray activity, and rotation period, to produce a posterior probability distribution of age. SEEDS observations discovered a substellar companion to one of our targets, \kappa And, a late B star. We do not detect any other substellar companions, but do find seven new close binary systems, of which one still needs to be confirmed. A detailed analysis of the statistics of this sample, and of the companion mass constraints given our age probability distributions and exoplanet cooling models, will be presented in a forthcoming paper.

To reveal the structures of a transition disk around a young stellar object in Lupus, Sz 91, we have performed aperture synthesis 345 GHz continuum and CO(3--2) observations with the Submillimeter Array ($\sim1\arcsec$--3$\arcsec$ resolution), and high-resolution imaging of polarized intensity at the $K_s$-band by using the HiCIAO instrument on the Subaru Telescope ($0\farcs25$ resolution). Our observations successfully resolved the inner and outer radii of the dust disk to be 65 AU and 170 AU, respectively, which indicates that Sz 91 is a transition disk source with one of the largest known inner holes. The model fitting analysis of the spectral energy distribution reveals an H$_2$ mass of $2.4\times10^{-3}$ $M_\sun$ in the cold ($T<$30 K) outer part at $65<r<170$ AU by assuming a canonical gas-to-dust mass ratio of 100, although a small amount ($>3\times10^{-9}$ $M_\sun$) of hot ($T\sim$180 K) dust possibly remains inside the inner hole of the disk. The structure of the hot component could be interpreted as either an unresolved self-luminous companion body (not directly detected in our observations) or a narrow ring inside the inner hole. Significant CO(3--2) emission with a velocity gradient along the major axis of the dust disk is concentrated on the Sz 91 position, suggesting a rotating gas disk with a radius of 420 AU. The Sz 91 disk is possibly a rare disk in an evolutionary stage immediately after the formation of protoplanets because of the large inner hole and the lower disk mass than other transition disks studied thus far.

Most exoplanets detected by direct imaging so far have been characterized by relatively hot (> ~1000 K) and cloudy atmospheres. A surprising feature in some of their atmospheres has been a distinct lack of methane, possibly implying non-equilibrium chemistry. Recently, we reported the discovery of a planetary companion to the Sun-like star GJ 504 using Subaru/HiCIAO within the SEEDS survey. The planet is substantially colder (<600 K) than previously imaged planets, and has indications of fewer clouds, which implies that it represents a new class of planetary atmospheres with expected similarities to late T-type brown dwarfs in the same temperature range. If so, one might also expect the presence of significant methane absorption, which is characteristic of such objects. Here, we report the detection of deep methane absorption in the atmosphere of GJ 504 b, using the Spectral Differential Imaging mode of HiCIAO to distinguish the absorption feature around 1.6 um. We also report updated JHK photometry based on new Ks-band data and a re-analysis of the existing data. The results support the notion that GJ 504 b has atmospheric properties distinct from other imaged exoplanets, and will become a useful reference object for future planets in the same temperature range.

We report a discovery of a companion candidate around one of {\it Kepler} Objects of Interest (KOIs), KOI-94, and results of our quantitative investigation of the possibility that planetary candidates around KOI-94 are false positives. KOI-94 has a planetary system in which four planetary detections have been reported by {\it Kepler}, suggesting that this system is intriguing to study the dynamical evolutions of planets. However, while two of those detections (KOI-94.01 and 03) have been made robust by previous observations, the others (KOI-94.02 and 04) are marginal detections, for which future confirmations with various techniques are required. We have conducted high-contrast direct imaging observations with Subaru/HiCIAO in $H$ band and detected a faint object located at a separation of $\sim0.6''$ from KOI-94. The object has a contrast of $\sim 1\times 10^{-3}$ in $H$ band, and corresponds to an M type star on the assumption that the object is at the same distance of KOI-94. Based on our analysis, KOI-94.02 is likely to be a real planet because of its transit depth, while KOI-94.04 can be a false positive due to the companion candidate. The success in detecting the companion candidate suggests that high-contrast direct imaging observations are important keys to examine false positives of KOIs. On the other hand, our transit light curve reanalyses lead to a better period estimate of KOI-94.04 than that on the KOI catalogue and show that the planetary candidate has the same limb darkening parameter value as the other planetary candidates in the KOI-94 system, suggesting that KOI-94.04 is also a real planet in the system.

We previously reported the direct detection of a low mass companion at a projected separation of 55+-2 AU around the B9 type star {\kappa} Andromedae. The properties of the system (mass ratio, separation) make it a benchmark for the understanding of the formation and evolution of gas giant planets and brown dwarfs on wide-orbits. We present new angular differential imaging (ADI) images of the Kappa Andromedae system at 2.146 (Ks), 3.776 (L'), 4.052 (NB 4.05) and 4.78 {\mu}m (M') obtained with Keck/NIRC2 and LBTI/LMIRCam, as well as more accurate near-infrared photometry of the star with the MIMIR instrument. We derive a more accurate J = 15.86 +- 0.21, H = 14.95 +- 0.13, Ks = 14.32 +- 0.09 mag for {\kappa} And b. We redetect the companion in all our high contrast observations. We confirm previous contrasts obtained at Ks and L' band. We derive NB 4.05 = 13.0 +- 0.2 and M' = 13.3 +- 0.3 mag and estimate Log10(L/Lsun) = -3.76 +- 0.06. We build the 1-5 microns spectral energy distribution of the companion and compare it to seven PHOENIX-based atmospheric models in order to derive Teff = 1900+100-200 K. Models do not set constrains on the surface gravity. ``Hot-start" evolutionary models predict masses of 14+25-2 MJup based on the luminosity and temperature estimates, and considering a conservative age range for the system (30+120-10 Myr). ``warm-start" evolutionary tracks constrain the mass to M >= 11 MJup. Therefore, the mass of {\kappa} Andromedae b mostly falls in the brown-dwarf regime, due to remaining uncertainties in age and mass-luminosity models. According to the formation models, disk instability in a primordial disk could account for the position and a wide range of plausible masses of {\kappa} And b.

Several exoplanets have recently been imaged at wide separations of >10 AU from their parent stars. These span a limited range of ages (<50 Myr) and atmospheric properties, with temperatures of 800--1800 K and very red colors (J - H > 0.5 mag), implying thick cloud covers. Furthermore, substantial model uncertainties exist at these young ages due to the unknown initial conditions at formation, which can lead to an order of magnitude of uncertainty in the modeled planet mass. Here, we report the direct imaging discovery of a Jovian exoplanet around the Sun-like star GJ 504, detected as part of the SEEDS survey. The system is older than all other known directly-imaged planets; as a result, its estimated mass remains in the planetary regime independent of uncertainties related to choices of initial conditions in the exoplanet modeling. Using the most common exoplanet cooling model, and given the system age of 160 [+350, -60] Myr, GJ 504 b has an estimated mass of 4 [+4.5, -1.0] Jupiter masses, among the lowest of directly imaged planets. Its projected separation of 43.5 AU exceeds the typical outer boundary of ~30 AU predicted for the core accretion mechanism. GJ 504 b is also significantly cooler (510 [+30, -20] K) and has a bluer color (J-H = -0.23 mag) than previously imaged exoplanets, suggesting a largely cloud-free atmosphere accessible to spectroscopic characterization. Thus, it has the potential of providing novel insights into the origins of giant planets, as well as their atmospheric properties.

We present near-infrared coronagraphic imaging polarimetry of RY Tau. The scattered light in the circumstellar environment was imaged at H-band at a high resolution (~0".05) for the first time, using Subaru-HiCIAO. The observed polarized intensity (PI) distribution shows a butterfly-like distribution of bright emission with an angular scale similar to the disk observed at millimeter wavelengths. This distribution is offset toward the blueshifted jet, indicating the presence of a geometrically thick disk or a remnant envelope, and therefore the earliest stage of the Class II evolutionary phase. We perform comparisons between the observed PI distribution and disk models with: (1) full radiative transfer code, using the spectral energy distribution (SED) to constrain the disk parameters; and (2) monochromatic simulations of scattered light which explore a wide range of parameters space to constrain the disk and dust parameters. We show that these models cannot consistently explain the observed PI distribution, SED, and the viewing angle inferred by millimeter interferometry. We suggest that the scattered light in the near-infrared is associated with an optically thin and geometrically thick layer above the disk surface, with the surface responsible for the infrared SED. Half of the scattered light and thermal radiation in this layer illuminates the disk surface, and this process may significantly affect the thermal structure of the disk.

We carried out an imaging survey for extrasolar planets around stars in the Pleiades (125 Myr, 135 pc) in the $H$ and $K_{S}$ bands using HiCIAO combined with the adaptive optics, AO188, on the Subaru telescope. We found 13 companion candidates fainter than 14.5 mag in the $H$ band around 9 stars. Five of these 13 were confirmed to be background stars by measurement of their proper motion. One was not found in the second epoch observation, and thus was not a background or companion object. One had multi-epoch image, but the precision of its proper motion was not sufficient to conclude whether it was background object. Four other candidates are waiting for second epoch observations to determine their proper motion. Finally, the remaining 2 were confirmed to be 60 $M_{J}$ brown dwarf companions orbiting around HD 23514 (G0) and HII 1348 (K5) respectively, as had been reported in previous studies. In our observations, the average detection limit for a point source was 20.3 mag in the $H$ band beyond 1''.5 from the central star. On the basis of this detection limit, we calculated the detection efficiency to be 90% for a planet with 6 to 12 Jovian masses and a semi-major axis of 50--1000 AU. For this we extrapolated the distribution of planet mass and semi-major axis derived from RV observations and adopted the planet evolution model of Baraffe et al. (2003). As there was no detection of a planet, we estimated the frequency of such planets to be less than 17.9% ($2\sigma$) around one star of the Pleiades cluster.

Debris disks around young main-sequence stars often have gaps and cavities which for a long time have been interpreted as possibly being caused by planets. In recent years, several giant planet discoveries have been made in systems hosting disks of precisely this nature, further implying that interactions with planets could be a common cause of such disk structures. As part of the SEEDS high-contrast imaging survey, we are surveying a population of debris disk-hosting stars with gaps and cavities implied by their spectral energy distributions, in order to attempt to spatially resolve the disk as well as to detect any planets that may be responsible for the disk structure. Here we report on intermediate results from this survey. Five debris disks have been spatially resolved, and a number of faint point sources have been discovered, most of which have been tested for common proper motion, which in each case has excluded physical companionship with the target stars. From the detection limits of the 50 targets that have been observed, we find that beta Pic b-like planets (~10 Mjup planets around G--A-type stars) near the gap edges are less frequent than 15--30%, implying that if giant planets are the dominant cause of these wide (27 AU on average) gaps, they are generally less massive than beta Pic b.

We demonstrate the ability of existing and planned telescopes, on the ground and in space, to directly image tidally heated exomoons orbiting gas-giant exoplanets. Tidally heated exomoons can plausibly be far more luminous than their host exoplanet and as much as 0.1% as bright as the system's stellar primary if it is a low mass star. Because emission from exomoons can be powered by tidal forces, they can shine brightly at arbitrarily large separations from the system's stellar primary with temperatures of several hundreds degrees Kelvin or even higher in extreme cases. Furthermore, these high temperatures can occur in systems that are billions of years old. Tidally heated exomoons may thus be far easier targets for direct imaging studies than giant exoplanets which must be both young and at a large projected separation from their primary to be accessible to current generation direct imaging studies. For example, the (warm) Spitzer Space Telescope and the next generation of ground based instruments could detect an exomoon roughly the size of the Earth at a temperature ~600K and a distance ~5 parsecs. Future mid-infrared space telescopes, such as JWST and SPICA, will be capable of directly imaging tidally heated exomoons around the nearest two dozen stars with brightness temperature >300K and R>1R_earth orbiting at >12AU from the primary star at a 5sigma confidence level in a 10,000 second integration. In addition it is possible that some of the exoplanets which have already been directly imaged are actually tidally heated exomoons or blends of such objects with hot young planets. If such exomoons exist and are sufficiently common (i.e. nearby), it may well be far easier to directly image an exomoon with surface conditions that allow the existence of liquid water than it will be to resolve an Earth-like planet in the classical Habitable Zone of its primary.

We present the first near infrared (NIR) spatially resolved images of the circumstellar transitional disk around SR21. These images were obtained with the Subaru HiCIAO camera, adaptive optics and the polarized differential imaging (PDI) technique. We resolve the disk in scattered light at H-band for stellocentric 0.1"<r<0.6" (12<r<75AU). We compare our results with previously published spatially-resolved 880 micron continuum Submillimeter Array (SMA) images that show an inner r<36AU cavity in SR21. Radiative transfer models reveal that the large disk depletion factor invoked to explain SR21's sub-mm cavity cannot be "universal" for all grain sizes. Even significantly more moderate depletions (delta=0.1, 0.01 relative to an undepleted disk) than those that reproduce the sub-mm cavity (delta~10^-6) are inconsistent with our H-band images when they are assumed to carry over to small grains, suggesting that surface grains scattering in the NIR either survive or are generated by whatever mechanism is clearing the disk midplane. In fact, the radial polarized intensity profile of our H-band observations is smooth and steeply inwardly-increasing (r^-3), with no evidence of a break at the 36AU sub-mm cavity wall. We hypothesize that this profile is dominated by an optically thin disk envelope or atmosphere component. We also discuss the compatibility of our data with the previously postulated existence of a sub-stellar companion to SR21 at r~10-20AU, and find that we can neither exclude nor verify this scenario. This study demonstrates the power of multiwavelength imaging of transitional disks to inform modeling efforts, including the debate over precisely what physical mechanism is responsible for clearing these disks of their large midplane grains.

We describe Algorithms for Calibration, Optimized Registration, and Nulling the Star in Angular Differential Imaging (ACORNS-ADI), a new, parallelized software package to reduce high-contrast imaging data, and its application to data from the SEEDS survey. We implement several new algorithms, including a method to register saturated images, a trimmed mean for combining an image sequence that reduces noise by up to ~20%, and a robust and computationally fast method to compute the sensitivity of a high-contrast observation everywhere on the field-of-view without introducing artificial sources. We also include a description of image processing steps to remove electronic artifacts specific to Hawaii2-RG detectors like the one used for SEEDS, and a detailed analysis of the Locally Optimized Combination of Images (LOCI) algorithm commonly used to reduce high-contrast imaging data. ACORNS-ADI is written in python. It is efficient and open-source, and includes several optional features which may improve performance on data from other instruments. ACORNS-ADI requires minimal modification to reduce data from instruments other than HiCIAO. It is freely available for download at www.github.com/t-brandt/acorns-adi under a BSD license.

We present Subaru/HiCIAO H-band high-contrast images of the debris disk around HIP 79977, whose pres- ence was recently inferred from an infrared excess. Our images resolve the disk for the first time, allowing characterization of its shape, size, and dust grain properties. We use angular differential imaging (ADI) to reveal the disk geometry in unpolarized light out to a radius of ~2", as well as polarized differential imaging (PDI) to measure the degree of scattering polarization out to ~1.5". In order to strike a favorable balance between suppression of the stellar halo and conservation of disk flux, we explore the application of principal component analysis (PCA) to both ADI and reference star subtraction. This allows accurate forward modeling of the effects of data reduction on simulated disk images, and thus direct comparison with the imaged disk. The resulting best-fit values and well-fitting intervals for the model parameters are a surface brightness power-law slope of S_out = -3.2 [-3.6,-2.9], an inclination of i = 84{\deg} [81{\deg},86{\deg}], a high Henyey-Greenstein forward-scattering parameter of g = 0.45 [0.35, 0.60], and a non-significant disk-star offset of u = 3.0 [-1.5, 7.5] AU = 24 [-13, 61] mas along the line of nodes. Furthermore, the tangential linear polarization along the disk rises from ~10% at 0.5" to ~45% at 1.5". These measurements paint a consistent picture of a disk of dust grains produced by collisional cascades and blown out to larger radii by stellar radiation pressure.

We study the variability of major atmospheric absorption features in the disk-integrated spectra of the Earth with future application to Earth-analogs in mind, concentrating on the diurnal timescale. We first analyze observations of the Earth provided by the EPOXI mission, and find 5-20% fractional variation of the absorption depths of H2O and O2 bands, two molecules that have major signatures in the observed range. From a correlation analysis with the cloud map data from the Earth Observing Satellite (EOS), we find that their variation pattern is primarily due to the uneven cloud cover distribution. In order to account for the observed variation quantitatively, we consider a simple opaque cloud model, which assumes that the clouds totally block the spectral influence of the atmosphere below the cloud layer, equivalent to assuming that the incident light is completely scattered at the cloud top level. The model is reasonably successful, and reproduces the EPOXI data from the pixel-level EOS cloud/water vapor data. A difference in the diurnal variability patterns of H2O and O2 bands is ascribed to the differing vertical and horizontal distribution of those molecular species in the atmosphere. On the Earth, the inhomogeneous distribution of atmospheric water vapor is due to the existence of its exchange with liquid and solid phases of H2O on the planet's surface on a timescale short compared to atmospheric mixing times. If such differences in variability patterns were detected in spectra of Earth-analogs, it would provide the information on the inhomogeneous composition of their atmospheres.

We present the first near-IR scattered light detection of the transitional disk associated with the Herbig Ae star MWC 758 using data obtained as part of the Strategic Exploration of Exoplanets and Disks with Subaru, and 1.1 micron HST/NICMOS data. While sub-millimeter studies suggested there is a dust-depleted cavity with r=0.35, we find scattered light as close as 0.1 (20-28 AU) from the star, with no visible cavity at H, K', or Ks. We find two small-scaled spiral structures which asymmetrically shadow the outer disk. We model one of the spirals using spiral density wave theory, and derive a disk aspect ratio of h ~ 0.18, indicating a dynamically warm disk. If the spiral pattern is excited by a perturber, we estimate its mass to be 5+3,-4 Mj, in the range where planet filtration models predict accretion continuing onto the star. Using a combination of non-redundant aperture masking data at L' and angular differential imaging with Locally Optimized Combination of Images at K' and Ks, we exclude stellar or massive brown dwarf companions within 300 mas of the Herbig Ae star, and all but planetary mass companions exterior to 0.5. We reach 5-sigma contrasts limiting companions to planetary masses, 3-4 MJ at 1.0 and 2 MJ at 1.55 using the COND models. Collectively, these data strengthen the case for MWC 758 already being a young planetary system.

We present the direct imaging discovery of an extrasolar planet, or possible low-mass brown dwarf, at a projected separation of 55 +/- 2 AU (1.058 +/- 0.007 arcsec) from the B9-type star Kappa And. The planet was detected with Subaru/HiCIAO during the SEEDS survey, and confirmed as a bound companion via common proper motion measurements. Observed near-infrared magnitudes of J = 16.3 +/- 0.3, H = 15.2 +/- 0.2, Ks = 14.6 +/- 0.4, and L' = 13.12 +/- 0.09 indicate a temperature of ~1700 K. The galactic kinematics of the host star are consistent with membership in the Columba association, implying a corresponding age of 30 +20 -10 Myr. The system age, combined with the companion photometry, points to a model-dependent companion mass ~12.8 MJup. The host star's estimated mass of 2.4-2.5 Msun places it among the most massive stars ever known to harbor an extrasolar planet or low-mass brown dwarf. While the mass of the companion is close to the deuterium burning limit, its mass ratio, orbital separation, and likely planet-like formation scenario imply that it may be best defined as a `Super-Jupiter' with properties similar to other recently discovered companions to massive stars.

We report high-resolution (0.07 arcsec) near-infrared polarized intensity images of the circumstellar disk around the star 2MASS J16042165-2130284 obtained with HiCIAO mounted on the Subaru 8.2 m telescope. We present our $H$-band data, which clearly exhibits a resolved, face-on disk with a large inner hole for the first time at infrared wavelengths. We detect the centrosymmetric polarization pattern in the circumstellar material as has been observed in other disks. Elliptical fitting gives the semimajor axis, semiminor axis, and position angle (P.A.) of the disk as 63 AU, 62 AU, and -14 $^{\circ}$, respectively. The disk is asymmetric, with one dip located at P.A.s of $\sim85^{\circ}$. Our observed disk size agrees well with a previous study of dust and CO emission at submillimeter wavelength with Submillimeter Array. Hence, the near-infrared light is interpreted as scattered light reflected from the inner edge of the disk. Our observations also detect an elongated arc (50 AU) extending over the disk inner hole. It emanates at the inner edge of the western side of the disk, extending inward first, then curving to the northeast. We discuss the possibility that the inner hole, the dip, and the arc that we have observed may be related to the existence of unseen bodies within the disk.

Through detailed radiative transfer modeling, we present a disk+cavity model to simultaneously explain both the SED and Subaru H-band polarized light imaging for the pre-transitional protoplanetary disk PDS 70. Particularly, we are able to match not only the radial dependence, but also the absolute scale, of the surface brightness of the scattered light. Our disk model has a cavity 65 AU in radius, which is heavily depleted of sub-micron-sized dust grains, and a small residual inner disk which produces a weak but still optically thick NIR excess in the SED. To explain the contrast of the cavity edge in the Subaru image, a factor of ~1000 depletion for the sub-micron-sized dust inside the cavity is required. The total dust mass of the disk may be on the order of 1e-4 M_sun, only weakly constrained due to the lack of long wavelength observations and the uncertainties in the dust model. The scale height of the sub-micron-sized dust is ~6 AU at the cavity edge, and the cavity wall is optically thick in the vertical direction at H-band. PDS 70 is not a member of the class of (pre-)transitional disks identified by Dong et al. (2012), whose members only show evidence of the cavity in the millimeter-sized dust but not the sub-micron-sized dust in resolved images. The two classes of (pre-)transitional disks may form through different mechanisms, or they may just be at different evolution stages in the disk clearing process.

We report that HAT-P-7 has a common proper motion stellar companion. The companion is located at $\sim3.9$ arcsec to the east and estimated as an M5.5V dwarf based on its colors. We also confirm the presence of the third companion, which was first reported by Winn et al. (2009), based on long-term radial velocity measurements. We revisit the migration mechanism of HAT-P-7b given the presence of those companions, and propose sequential Kozai migration as a likely scenario in this system. This scenario may explain the reason for an outlier in the discussion of the spin-orbit alignment timescale for HAT-P-7b by Albrecht et al. (2012).

We present high resolution H-band polarized intensity (PI; FWHM = 0."1: 14 AU) and L'-band imaging data (FWHM = 0."11: 15 AU) of the circumstellar disk around the weak-lined T Tauri star PDS 70 in Centaurus at a radial distance of 28 AU (0."2) up to 210 AU (1."5). In both images, a giant inner gap is clearly resolved for the first time, and the radius of the gap is ~70 AU. Our data show that the geometric center of the disk shifts by ~6 AU toward the minor axis. We confirm that the brown dwarf companion candidate to the north of PDS 70 is a background star based on its proper motion. As a result of SED fitting by Monte Carlo radiative transfer modeling, we infer the existence of an optically thick inner disk at a few AU. Combining our observations and modeling, we classify the disk of PDS 70 as a pre-transitional disk. Furthermore, based on the analysis of L'-band imaging data, we put an upper limit mass of companions at ~30 to ~50MJ within the gap. Taking account of the presence of the large and sharp gap, we suggest that the gap could be formed by dynamical interactions of sub-stellar companions or multiple unseen giant planets in the gap.

Ground-based telescopes equipped with adaptive-optics (AO) systems and specialized science cameras are now capable of directly detecting extrasolar planets. We present the expected scientific capabilities of CHARIS, the Coronagraphic High Angular Resolution Imaging Spectrograph, which is being built for the Subaru 8.2 m telescope of the National Astronomical Observatory of Japan. CHARIS will be implemented behind the new extreme adaptive optics system at Subaru, SCExAO, and the existing 188-actuator system AO188. CHARIS will offer three observing modes over near-infrared wavelengths from 0.9 to 2.4 microns (the y-, J-, H-, and K-bands), including a low-spectral-resolution mode covering this entire wavelength range and a high-resolution mode within a single band. With these capabilities, CHARIS will offer exceptional sensitivity for discovering giant exoplanets, and will enable detailed characterization of their atmospheres. CHARIS, the only planned high-contrast integral field spectrograph on an 8m-class telescope in the Northern Hemisphere, will complement the similar instruments such as Project 1640 at Palomar, and GPI and SPHERE in Chile.

We present H-band polarimetric imagery of UX Tau A taken with HiCIAO/AO188 on the Subaru Telescope. UX Tau A has been classified as a pre-transitional disk object, with a gap structure separating its inner and outer disks. Our imagery taken with the 0.15 (21 AU) radius coronagraphic mask has revealed a strongly polarized circumstellar disk surrounding UX Tau A which extends to 120 AU, at a spatial resolution of 0.1 (14 AU). It is inclined by 46 \pm 2 degree as the west side is nearest. Although SED modeling and sub-millimeter imagery suggested the presence of a gap in the disk, with the inner edge of the outer disk estimated to be located at 25 - 30 AU, we detect no evidence of a gap at the limit of our inner working angle (23 AU) at the near-infrared wavelength. We attribute the observed strong polarization (up to 66 %) to light scattering by dust grains in the disk. However, neither polarization models of the circumstellar disk based on Rayleigh scattering nor Mie scattering approximations were consistent with the observed azimuthal profile of the polarization degrees of the disk. Instead, a geometric optics model of the disk with nonspherical grains with the radii of 30 micron meter is consistent with the observed profile. We suggest that the dust grains have experienced frequent collisional coagulations and have grown in the circumstellar disk of UX Tau A.

One of the key predictions of modeling from the IR excess of Herbig Ae stars is that for protoplanetary disks, where significant grain growth and settling has occurred, the dust disk has flattened to the point that it can be partially or largely shadowed by the innermost material at or near the dust sublimation radius. When the self-shadowing has already started, the outer disk is expected to be detected in scattered light only in the exceptional cases that the scale height of the dust disk at the sublimation radius is smaller than usual. High-contrast imaging combined with the IR spectral energy distribution allow us to measure the degree of flattening of the disk, as well as to determine the properties of the outer disk. We present polarimetric differential imaging in $H$ band obtained with Subaru/HiCIAO of one such system, MWC 480. The HiCIAO data were obtained at a historic minimum of the NIR excess. The disk is detected in scattered light from 0\farcs2-1\farcs0 (27.4-137AU). Together with the marginal detection of the disk from 1998 February 24 by HST/NICMOS, our data constrain the opening half angle for the disk to lie between 1.3$\leq\theta\leq 2.2^\circ$. When compared with similar measures in CO for the gas disk from the literature, the dust disk subtends only $\sim$30% of the gas disk scale height (H/R$\sim$0.03). Such a dust disk is a factor of 5-7 flatter than transitional disks, which have structural signatures that giant planets have formed.

Life arose on Earth sometime in the first few hundred million years after the young planet had cooled to the point that it could support water-based organisms on its surface. The early emergence of life on Earth has been taken as evidence that the probability of abiogenesis is high, if starting from young-Earth-like conditions. We revisit this argument quantitatively in a Bayesian statistical framework. By constructing a simple model of the probability of abiogenesis, we calculate a Bayesian estimate of its posterior probability, given the data that life emerged fairly early in Earth's history and that, billions of years later, curious creatures noted this fact and considered its implications. We find that, given only this very limited empirical information, the choice of Bayesian prior for the abiogenesis probability parameter has a dominant influence on the computed posterior probability. Although terrestrial life's early emergence provides evidence that life might be common in the Universe if early-Earth-like conditions are, the evidence is inconclusive and indeed is consistent with an arbitrarily low intrinsic probability of abiogenesis for plausible uninformative priors. Finding a single case of life arising independently of our lineage (on Earth, elsewhere in the Solar System, or on an extrasolar planet) would provide much stronger evidence that abiogenesis is not extremely rare in the Universe.

Existing and planned optical telescopes and surveys can detect artificially-illuminated objects comparable in total brightness to a major terrestrial city out to the outskirts of the Solar System. Orbital parameters of Kuiper belt objects (KBOs) are routinely measured to exquisite precisions of <10^{-3}. Here we propose to measure the variation of the observed flux F from such objects as a function of their changing orbital distances D. Sunlight-illuminated objects will show a logarithmic slope alpha=(dlogF/dlog D)=-4 whereas artificially-illuminated objects should exhibit alpha=-2. Planned surveys using the proposed LSST will provide superb data that would allow measurement of alpha for thousands of KBOs. If objects with alpha=-2 are found, follow-up observations can measure their spectra to determine if they are illuminated by artificial lighting. The search can be extended beyond the Solar System with future generations of telescopes on the ground and in space, which would be capable of detecting phase modulation due to very strong artificial illumination on the night-side of planets as they orbit their parent stars.

We present high-resolution, H-band, imaging observations, collected with Subaru/HiCIAO, of the scattered light from the transitional disk around SAO 206462 (HD 135344B). Although previous sub-mm imagery suggested the existence of the dust-depleted cavity at r~46AU, our observations reveal the presence of scattered light components as close as 0.2" (~28AU) from the star. Moreover, we have discovered two small-scale spiral structures lying within 0.5" (~70AU). We present models for the spiral structures using the spiral density wave theory, and derive a disk aspect ratio of h~0.1, which is consistent with previous sub-mm observations. This model can potentially give estimates of the temperature and rotation profiles of the disk based on dynamical processes, independently from sub-mm observations. It also predicts the evolution of the spiral structures, which can be observable on timescales of 10-20 years, providing conclusive tests of the model. While we cannot uniquely identify the origin of these spirals, planets embedded in the disk may be capable of exciting the observed morphology. Assuming that this is the case, we can make predictions on the locations and, possibly, the masses of the unseen planets. Such planets may be detected by future multi-wavelengths observations.

Transitional circumstellar disks around young stellar objects have a distinctive infrared deficit around 10 microns in their Spectral Energy Distributions (SED), recently measured by the Spitzer Infrared Spectrograph (IRS), suggesting dust depletion in the inner regions. These disks have been confirmed to have giant central cavities by imaging of the submillimeter (sub-mm) continuum emission using the Submillimeter Array (SMA). However, the polarized near-infrared scattered light images for most objects in a systematic IRS/SMA cross sample, obtained by HiCIAO on the Subaru telescope, show no evidence for the cavity, in clear contrast with SMA and Spitzer observations. Radiative transfer modeling indicates that many of these scattered light images are consistent with a smooth spatial distribution for micron-sized grains, with little discontinuity in the surface density of the micron-sized grains at the cavity edge. Here we present a generic disk model that can simultaneously account for the general features in IRS, SMA, and Subaru observations. Particularly, the scattered light images for this model are computed, which agree with the general trend seen in Subaru data. Decoupling between the spatial distributions of the micron-sized dust and mm-sized dust inside the cavity is suggested by the model, which, if confirmed, necessitates a mechanism, such as dust filtration, for differentiating the small and big dust in the cavity clearing process. Our model also suggests an inwardly increasing gas-to-dust-ratio in the inner disk, and different spatial distributions for the small dust inside and outside the cavity, echoing the predictions in grain coagulation and growth models.

We present high-contrast images of HR 4796 A taken with Subaru/HiCIAO in H-band, resolving the debris disk in scattered light. The application of specialized angular differential imaging methods (ADI) allows us to trace the inner edge of the disk with high precision, and reveals a pair of "streamers" extending radially outwards from the ansae. Using a simple disk model with a power-law surface brightness profile, we demonstrate that the observed streamers can be understood as part of the smoothly tapered outer boundary of the debris disk, which is most visible at the ansae. Our observations are consistent with the expected result of a narrow planetesimal ring being ground up in a collisional cascade, yielding dust with a wide range of grain sizes. Radiation forces leave large grains in the ring and push smaller grains onto elliptical, or even hyperbolic trajectories. We measure and characterize the disk's surface brightness profile, and confirm the previously suspected offset of the disk's center from the star's position along the ring's major axis. Furthermore, we present first evidence for an offset along the minor axis. Such offsets are commonly viewed as signposts for the presence of unseen planets within a disk's cavity. Our images also offer new constraints on the presence of companions down to the planetary mass regime (~9 Jupiter masses at 0.5", ~3 Jupiter masses at 1").

Radial velocity (RV) observations of an exoplanet system giving a value of M_T sin(i) condition (i.e. give information about) not only the planet's true mass M_T but also the value of sin(i) (where i is the orbital inclination angle). Thus the value of sin(i) for a system with any particular observed value of M_T sin(i) cannot be assumed to be drawn randomly from a uniform distribution between zero and unity (corresponding to an isotropic i distribution). The actual distribution from which it is drawn depends on the intrinsic distribution of M_T for the exoplanet population being studied. We give a simple Bayesian derivation of this relationship and apply it to several "toy models" for the (currently unknown) intrinsic distribution of M_T. The results show that the effect can be an important one. For example, even for simple power-law distributions of M_T, the median value of sin(i) in an observed RV sample can vary between 0.25 and 0.71 (as compared to the 0.5 value for an isotropic i distribution) for indices of the power-law in the quite plausible range between -2 and -0.5, respectively. Over the same range of indicies, the 95% confidence upper bound on M_T ranges from 4.5 to 400 times larger than M_T sin(i), respectively, due to sin(i) uncertainty alone. More complex, but still simple and plausible, distributions of M_T yield still more complicated and less intuitive $sin(i)$ distributions. In particular, if the M_T distribution contains any characteristic mass scale M_c, the sin(i) distribution will depend on the ratio of M_T sin(i) to $M_c$, often in a non-trivial way. Our qualitative conclusion is that RV studies of exoplanets, both individual objects and statistical samples, should regard the sin(i) factor as more than a "numerical constant of order unity" with simple and well understood statistical properties. (abridged)

As a test-bed for future investigations of directly imaged terrestrial exoplanets, we present the recovery of the surface components of the Earth from multi-band diurnal light curves obtained with the EPOXI spacecraft. We find that the presence and longitudinal distribution of ocean, soil and vegetation are reasonably well reproduced by fitting the observed color variations with a simplified model composed of a priori known albedo spectra of ocean, soil, vegetation, snow and clouds. The effect of atmosphere, including clouds, on light scattered from surface components is modeled using a radiative transfer code. The required noise levels for future observations of exoplanets are also determined. Our model-dependent approach allows us to infer the presence of major elements of the planet (in the case of the Earth, clouds and ocean) with observations having S/N $\gtrsim 10$ in most cases and with high confidence if S/N $\gtrsim 20$. In addition, S/N $\gtrsim 100$ enables us to detect the presence of components other than ocean and clouds in a fairly model-independent way. Degradation of our inversion procedure produced by cloud cover is also quantified. While cloud cover significantly dilutes the magnitude of color variations compared to the cloudless case, the pattern of color changes remains. Therefore, the possibility of investigating surface features through light curve fitting remains even for exoplanets with cloud cover similar to the Earth's.

The measurement of the light scattered from extrasolar planets informs atmospheric and formation models. With the discovery of many hot Jupiter planets orbiting nearby stars, this motivates the development of robust methods of characterisation from follow up observations. In this paper we discuss two methods for determining the planetary albedo in transiting systems. First, the most widely used method for measuring the light scattered by hot Jupiters (Collier Cameron et al.) is investigated for application for typical echelle spectra of a transiting planet system, showing that detection requires high signal-to-noise ratio data of bright planets. Secondly a new Fourier analysis method is also presented, which is model-independent and utilises the benefits of the reduced number of unknown parameters in transiting systems. This approach involves solving for the planet and stellar spectra in Fourier space by least-squares. The sensitivities of the methods are determined via Monte Carlo simulations for a range of planet-to-star fluxes. We find the Fourier analysis method to be better suited to the ideal case of typical observations of a well constrained transiting system than the Collier Cameron et al. method. We apply the Fourier analysis method for extracting the light scattered by transiting hot Jupiters from high resolution spectra to echelle spectra of HD 209458 and HD 189733. Unfortunately we are unable to improve on the previous upper limit of the planet-to-star flux for HD 209458b set by space-based observations. A 1{\sigma}upper limit on the planet-to-star flux of HD 189733b is measured in the wavelength range of 558.83-599.56 nm yielding {\epsilon} < 4.5 \times 10-4. Improvement in the measurement of the upper limit of the planet-to-star flux of this system, with ground-based capabilities, requires data with a higher signal-to-noise ratio, and increased stability of the telescope.

GJ 758 B is a cold (~600K) companion to a Sun-like star at 29 AU projected separation, which was recently detected with high-contrast imaging. Here we present photometry of the companion in seven photometric bands from Subaru/HiCIAO, Gemini/NIRI and Keck/NIRC2, providing a rich sampling of the spectral energy distribution in the 1-5 micron wavelength range. A clear detection at 1.58 micron combined with an upper limit at 1.69 micron shows methane absorption in the atmosphere of the companion. The mass of the companion remains uncertain, but an updated age estimate indicates that the most likely mass range is ~30-40 Mjup. In addition, we present an updated astrometric analysis that imposes tighter constraints on GJ 758 B's orbit and identifies the proposed second candidate companion, "GJ 758 C", as a background star.

We report high-resolution 1.6 $\micron$ polarized intensity ($PI$) images of the circumstellar disk around the Herbig Ae star AB Aur at a radial distance of 22 AU ($0."15$) up to 554 AU (3.$"$85), which have been obtained by the high-contrast instrument HiCIAO with the dual-beam polarimetry. We revealed complicated and asymmetrical structures in the inner part ($\lesssim$140 AU) of the disk, while confirming the previously reported outer ($r$ $\gtrsim$200 AU) spiral structure. We have imaged a double ring structure at $\sim$40 and $\sim$100 AU and a ring-like gap between the two. We found a significant discrepancy of inclination angles between two rings, which may indicate that the disk of AB Aur is warped. Furthermore, we found seven dips (the typical size is $\sim$45 AU or less) within two rings as well as three prominent $PI$ peaks at $\sim$40 AU. The observed structures, including a bumpy double ring, a ring-like gap, and a warped disk in the innermost regions, provide essential information for understanding the formation mechanism of recently detected wide-orbit ($r$ $>$20 AU) planets.

We present H- and Ks-band imaging data resolving the gap in the transitional disk around LkCa 15, revealing the surrounding nebulosity. We detect sharp elliptical contours delimiting the nebulosity on the inside as well as the outside, consistent with the shape, size, ellipticity, and orientation of starlight reflected from the far-side disk wall, whereas the near-side wall is shielded from view by the disk's optically thick bulk. We note that forward-scattering of starlight on the near-side disk surface could provide an alternate interpretation of the nebulosity. In either case, this discovery provides confirmation of the disk geometry that has been proposed to explain the spectral energy distributions (SED) of such systems, comprising an optically thick outer disk with an inner truncation radius of ~46 AU enclosing a largely evacuated gap. Our data show an offset of the nebulosity contours along the major axis, likely corresponding to a physical pericenter offset of the disk gap. This reinforces the leading theory that dynamical clearing by at least one orbiting body is the cause of the gap. Based on evolutionary models, our high-contrast imagery imposes an upper limit of 21 Jupiter masses on companions at separations outside of 0.1" and of 13 Jupiter masses outside of 0.2". Thus, we find that a planetary system around LkCa 15 is the most likely explanation for the disk architecture.

Radial velocity (RV) observations of an exoplanet system giving a value of M_T sin(i) condition (ie. give information about) not only the planet's true mass M_T but also the value of sin(i) for that system (where i is the orbital inclination angle). Thus the value of sin(i) for a system with any particular observed value of M_T sin(i) cannot be assumed to be drawn randomly from a distribution corresponding to an isotropic i distribution, i.e. the presumptive prior distribution . Rather, the posterior distribution from which it is drawn depends on the intrinsic distribution of M_T for the exoplanet population being studied. We give a simple Bayesian derivation of this relationship and apply it to several "toy models" for the (currently unknown) intrinsic distribution of M_T. The results show that the effect can be an important one. For example, even for simple power-law distributions of M_T, the median value of sin(i) in an observed RV sample can vary between 0.860 and 0.023 (as compared to the 0.866 value for an isotropic i distribution) for indices (alpha) of the power-law in the range between -2 and +1, respectively. Over the same range of indicies, the 95% confidence interval on M_T varies from 1.002-4.566 (alpha = -2) to 1.13-94.34 (alpha = +1) times larger than M_T sin(i) due to sin(i) uncertainty alone. Our qualitative conclusion is that RV studies of exoplanets, both individual objects and statistical samples, should regard the sin(i) factor as more than a "numerical constant of order unity" with simple and well understood statistical properties. We argue that reports of M_T sin(i) determinations should be accompanied by a statement of the corresponding confidence bounds on M_T at, say, the 95% level based on an explicitly stated assumed form of the true M_T distribution in order to more accurately reflect the mass uncertainties associated with RV studies.

Characterizing the surfaces of rocky exoplanets via the scattered light will be an essential challenge to investigate the existence of life on habitable exoplanets. We present a simple reconstruction method for fractional areas of different surface types from photometric variations, or colors, of a second Earth. We create mock light curves for Earth without clouds using empirical data. Then these light curves are fitted to the isotropic scattering model consisting of 4 surface types: ocean, soil, snow and vegetation. In an idealized situation where the photometric errors are only photon shot noise, we are able to reproduce the fractional areas of those components fairly well. We may be even able to detect a signature of vegetation from the distinct feature of photosynthesis on the Earth, known as the red edge. In our reconstruction method, Rayleigh scattering due to the atmosphere has an important effect, and for terrestrial exoplanets with atmosphere similar to our Earth, it is possible to estimate the presence of oceans and an atmosphere simultaneously.

We present results of direct imaging observations for HAT-P-7 taken with the Subaru HiCIAO and the Calar Alto AstraLux. Since the close-in transiting planet HAT-P-7b was reported to have a highly tilted orbit, massive bodies such as giant planets, brown dwarfs, or a binary star are expected to exist in the outer region of this system. We show that there are indeed two candidates for distant faint stellar companions around HAT-P-7. We discuss possible roles played by such companions on the orbital evolution of HAT-P-7b. We conclude that as there is a third body in the system as reported by Winn et al. (2009, ApJL, 763, L99), the Kozai migration is less likely while planet-planet scattering is possible.

Planetesimal formation is a common by-product of the star formation process. Taking the dynamical history of the Solar system as a guideline -- in which the planetesimal belts were heavily depleted due to graviational perburbation with the giant planets -- and assuming similar processes have take place in other planetary systems, one would expect the interestellar space to be filled with extra-solar planetesimals. However, not a single one of these objects has been detected so far entering the Solar system, even though it would clearly be distinguishable from a Solar system comet due to its highly hyperbolic orbit. The Large Synoptic Survey Telescope (LSST) will provide wide coverage maps of the sky to a very high sensitivity, ideal to detect moving objects like comets, both active and inactive. In anticipation of these observations, we estimate how many inactive "interstellar comets" might be detected during the duration of the survey. The calculation takes into account estimates (from observations and models) of the number density of stars, the amount of solids available to form planetesimals, the frequency of planet and planetesimal formation, the efficiency of planetesimal ejection, and the possible size distribution of these small bodies.

We present an improved measurement of the Rossiter-McLaughlin effect for the exoplanetary system HD 17156, based on radial-velocity data gathered with the Subaru 8.2m telescope throughout the planetary transit of UT 2008 November 7. The data allow for a precise and independent determination of the projected spin-orbit angle of this system: $\lambda = 10.0^{\circ} \pm 5.1^{\circ}$. This result supersedes the previous claim of $\lambda = 62^{\circ} \pm 25^{\circ}$ by Narita et al., which was based on lower-precision data with poor statistics. Thus the stellar spin and planetary orbital axes of the HD 17156 system are likely to be well-aligned, despite the planet's large orbital eccentricity suggesting a history of strong dynamical interactions.

We present photometric and spectroscopic observations of the 2009 Feb. 2 transit of the exoplanet XO-3b. The new data show that the planetary orbital axis and stellar rotation axis are misaligned, as reported earlier by Hebrard and coworkers. We find the angle between the sky projections of the two axes to be 37.3 +/- 3.7 deg, as compared to the previously reported value of 70 +/- 15 deg. The significance of this discrepancy is unclear because there are indications of systematic effects. XO-3b is the first exoplanet known to have a highly inclined orbit relative to the equatorial plane of its parent star, and as such it may fulfill the predictions of some scenarios for the migration of massive planets into close-in orbits. We revisit the statistical analysis of spin-orbit alignment in hot-Jupiter systems. Assuming the stellar obliquities to be drawn from a single Rayleigh distribution, we find the mode of the distribution to be 13^{+5}_{-2} deg. However, it remains the case that a model representing two different migration channels--in which some planets are drawn from a perfectly-aligned distribution and the rest are drawn from an isotropic distribution--is favored over a single Rayleigh distribution.

The identification of an extrasolar planet as Earth-like will depend on the detection of atmospheric signatures or surface non-uniformities. In this paper we present spatially unresolved flux light curves of Earth for the purpose of studying a prototype extrasolar terrestrial planet. Our monitoring of the photometric variability of earthshine revealed changes of up to 23 % per hour in the brightness of Earth's scattered light at around 600 nm, due to the removal of specular reflection from the view of the Moon. This variability is accompanied by reddening of the spectrum, and results from a change in surface properties across the continental boundary between the Indian Ocean and Africa's east coast. Our results based on earthshine monitoring indicate that specular reflection should provide a useful tool in determining the presence of liquid water on extrasolar planets via photometric observations.

We present new spectroscopic and photometric observations of the HAT-P-1 planetary system. Spectra obtained during three transits exhibit the Rossiter-McLaughlin effect, allowing us to measure the angle between the sky projections of the stellar spin axis and orbit normal, \lambda = 3.7 +/- 2.1 degrees. The small value of \lambda for this and other systems suggests that the dominant planet migration mechanism preserves spin-orbit alignment. Using two new transit light curves, we refine the transit ephemeris and reduce the uncertainty in the orbital period by an order of magnitude. We find a upper limit on the orbital eccentricity of 0.067, with 99% confidence, by combining our new radial-velocity measurements with those obtained previously.

We investigate the estimation of orbital parameters by least-$\chi^2$ Keplerian fits to radial velocity (RV) data using synthetic data sets. We find that while the fitted period is fairly accurate, the best-fit eccentricity and $M_p\sin i$ are systematically biased upward from the true values for low signal-to-noise ratio $K/\sigma\lesssim 3$ and moderate number of observations $N_{\rm obs}\lesssim 60$, leading to a suppression of the number of nearly circular orbits. Assuming intrinsic distributions of orbital parameters, we generate a large number of mock RV data sets and study the selection effect on the eccentricity distribution. We find the overall detection efficiency only mildly decreases with eccentricity. This is because although high eccentricity orbits are more difficult to sample, they also have larger RV amplitudes for fixed planet mass and orbital semi-major axis. Thus the primary source of uncertainties in the eccentricity distribution comes from biases in Keplerian fits to detections with low-amplitude and/or small $N_{\rm obs}$, rather than from selection effects. Our results suggest that the abundance of low-eccentricity exoplanets may be underestimated in the current sample and we urge caution in interpreting the eccentricity distributions of low-amplitude detections in future RV samples.

We monitored the Doppler shift of the G0V star TrES-2 throughout a transit of its giant planet. The anomalous Doppler shift due to stellar rotation (the Rossiter-McLaughlin effect) is discernible in the data, with a signal-to-noise ratio of 2.9, even though the star is a slow rotator. By modeling this effect we find that the planet's trajectory across the face of the star is tilted by -9 +/- 12 degrees relative to the projected stellar equator. With 98% confidence, the orbit is prograde.

We use the high magnification event seen in the 1999 OGLE campaign light curve of image C of the quadruply imaged gravitational lens Q2237+0305 to study the structure of the quasar engine. We have obtained g'- and r'-band photometry at the Apache Point Observatory 3.5m telescope where we find that the event has a smaller amplitude in the r'-band than in the g'- and OGLE V-bands. By comparing the light curves with microlensing simulations we obtain constraints on the sizes of the quasar regions contributing to the g'- and r'-band flux. Assuming that most of the surface mass density in the central kiloparsec of the lensing galaxy is due to stars and by modeling the source with a Gaussian profile, we obtain for the Gaussian width 1.20 x 10^15 sqrt(M/0.1M_sun)cm < sigma_g' < 7.96 x 10^15 sqrt(M/0.1Msun) cm, where M is the mean microlensing mass, and a ratio sigma_r'/sigma_g'=1.25^{+0.45}_{-0.15}. With the limits on the velocity of the lensing galaxy from Gil-Merino et al. (2005) as our only prior, we obtain 0.60 x 10^15 sqrt(M/0.1Msun) cm < sigma_g' < 1.57 x 10^15 sqrt(M/0.1Msun) cm and a ratio sigma_r'/sigma_g'=1.45^{+0.90}_{-0.25} (all values at 68 percent confidence). Additionally, from our microlensing simulations we find that, during the chromatic microlensing event observed, the continuum emitting region of the quasar crossed a caustic at >72 percent confidence.

We present a new technique for detecting scattered starlight from transiting, close-orbiting extrasolar giant planets (CEGPs) that has the virtues of simplicity, robustness, linearity, and model-independence. Given a series of stellar spectra obtained over various phases of the planetary orbit, the goal is to measure the strength of the component scattered by the planet relative to the component coming directly from the star. We use two complementary strategies, both of which rely on the predictable Doppler shifts of both components and on combining the results from many spectral lines and many exposures. In the first strategy, we identify segments of the stellar spectrum that are free of direct absorption lines and add them after Doppler-shifting into the planetary frame. In the second strategy, we compare the distribution of equivalent-width ratios of the scattered and direct components. Both strategies are calibrated with a ``null test'' in which scrambled Doppler shifts are applied to the spectral segments. As an illustrative test case, we apply our technique to spectra of HD 209458 taken when the planet was near opposition (with orbital phases ranging from 11 to 34$\arcdeg$, where 0$\arcdeg$ is at opposition), finding that the planet-to-star flux ratio is $(1.4 \pm 2.9)\times10^{-4}$ in the wavelength range 554$-$681 nm. This corresponds to a geometric albedo of $0.8 \pm 1.6$, assuming the phase function of a Lambert sphere. Although the result is not statistically significant, the achieved sensitivity and relatively small volume of data upon which it is based are very encouraging for future ground-based spectroscopic studies of scattered light from transiting CEGP systems.

We present cosmological results from the statistics of lensed quasars in the Sloan Digital Sky Survey (SDSS) Quasar Lens Search. By taking proper account of the selection function, we compute the expected number of quasars lensed by early-type galaxies and their image separation distribution assuming a flat universe, which is then compared with 7 lenses found in the SDSS Data Release 3 to derive constraints on dark energy under strictly controlled criteria. For a cosmological constant model (w=-1) we obtain \Omega_\Lambda=0.74^{+0.11}_{-0.15}(stat.)^{+0.13}_{-0.06}(syst.). Allowing w to be a free parameter we find \Omega_M=0.26^{+0.07}_{-0.06}(stat.)^{+0.03}_{-0.05}(syst.) and w=-1.1\pm0.6(stat.)^{+0.3}_{-0.5}(syst.) when combined with the constraint from the measurement of baryon acoustic oscillations in the SDSS luminous red galaxy sample. Our results are in good agreement with earlier lensing constraints obtained using radio lenses, and provide additional confirmation of the presence of dark energy consistent with a cosmological constant, derived independently of type Ia supernovae.

We report the first results of our systematic search for strongly lensed quasars using the spectroscopically confirmed quasars in the Sloan Digital Sky Survey (SDSS). Among 46,420 quasars from the SDSS Data Release 3 (~4188 deg^2), we select a subsample of 22,683 quasars that are located at redshifts between 0.6 and 2.2 and are brighter than the Galactic extinction corrected i-band magnitude of 19.1. We identify 220 lens candidates from the quasar subsample, for which we conduct extensive and systematic follow-up observations in optical and near-infrared wavebands, in order to construct a complete lensed quasar sample at image separations between 1'' and 20'' and flux ratios of faint to bright lensed images larger than 10^{-0.5}. We construct a statistical sample of 11 lensed quasars. Ten of these are galaxy-scale lenses with small image separations (~1''-2'') and one is a large separation (15'') system which is produced by a massive cluster of galaxies, representing the first statistical sample of lensed quasars including both galaxy- and cluster-scale lenses. The Data Release 3 spectroscopic quasars contain an additional 11 lensed quasars outside the statistical sample.

We report a measurement of the Rossiter--McLaughlin effect in the transiting extrasolar planetary system TrES-1, via simultaneous spectroscopic and photometric observations with the Subaru and MAGNUM telescopes. By modeling the radial velocity anomaly that was observed during a transit, we determine the sky-projected angle between the stellar spin axis and the planetary orbital axis to be $\lambda = 30 \pm 21$ [deg]. This is the third case for which $\lambda$ has been measured in a transiting exoplanetary system, and the first demonstration that such measurements are possible for relatively faint host stars ($V \sim 12$, as compared to $V \sim 8$ for the other systems). We also derive a time of mid-transit, constraints on the eccentricity of the TrES-1b orbit ($e = 0.048 \pm 0.025$), and upper limits on the mass of the Trojan companions ($\lesssim$14 $M_{\oplus}$) at the 3$\sigma$ level.

The short-period exoplanet HD 147506b (also known as HAT-P-2b) has an eccentric orbit, raising the possibility that it migrated through planet-planet scattering or Kozai oscillations accompanied by tidal dissipation. Either of these scenarios could have significantly tilted the orbit relative to the host star's equatorial plane. Here we present spectroscopy of a transit of HD 147506b, and assess the spin-orbit alignment via the Rossiter-McLaughlin effect. We find the sky projections of the stellar spin axis and orbital axis to be aligned within 14 deg. Thus we find no corroborating evidence for scattering or Kozai migration, although these scenarios cannot be ruled out with the present data.

The chemical composition of Earth's atmosphere has undergone substantial evolution over the course of its history. It is possible, even likely, that terrestrial planets in other planetary systems have undergone similar changes; consequently, the age distribution of nearby stars is an important consideration in designing surveys for Earth-analogues. Valenti & Fischer (2005) provide age and metallicity estimates for 1039 FGK dwarfs in the Solar Neighbourhood. Using the Hipparcos catalogue as a reference to calibrate potential biases, we have extracted volume-limited samples of nearby stars from the Valenti-Fischer dataset. Unlike other recent investigations, our analysis shows clear evidence for an age-metallicity relation in the local disk, albeit with substantial dispersion at any epoch. The mean metallicity increases from -0.3 dex at a lookback time of ~10 Gyrs to +0.15 dex at the present day. Supplementing the Valenti-Fischer measurements with literature data to give a complete volume-limited sample, the age distribution of nearby FGK dwarfs is broadly consistent with a uniform star-formation rate over the history of the Galactic disk. In striking contrast, most stars known to have planetary companions are younger than 5 Gyrs; however, stars with planetary companions within 0.4 AU have a significantly flatter age distribution, indicating that those systems are stable on timescales of many Gyrs. Several of the older, lower metallicity host stars have enhanced [alpha/Fe] ratios, implying membership of the thick disk. If the frequency of terrestrial planets is also correlated with stellar metallicity, then the median age of such planetary system is likely to be ~3 Gyrs. We discuss the implications of this hypothesis in designing searches for Earth analogues among the nearby stars.

We present spectroscopy of a transit of the exoplanet HD 189733b. By modeling the Rossiter-McLaughlin effect (the anomalous Doppler shift due to the partial eclipse of the rotating stellar surface), we find the angle between the sky projections of the stellar spin axis and orbit normal to be lambda = -1.4 +/- 1.1 deg. This is the third case of a ``hot Jupiter'' for which lambda has been measured. In all three cases lambda is small, ruling out random orientations with 99.96% confidence, and suggesting that the inward migration of hot Jupiters generally preserves spin-orbit alignment.

We compare the most successful and widely used map of Galactic dust extinction, provided by Schlegel, Finkbeiner & Davis (1998; hereafter SFD), to the galaxy number counts in the Sloan Digital Sky Survey (SDSS) photometric/spectroscopic DR4 sample. We divide the SDSS survey area into 69 disjoint subregions according to the dust extinction provided by SFD and compare the surface number density of galaxies in each subregion. As expected, the galaxy surface number density decreases with increasing extinction but only for SFD extinction values above about 0.1 to 0.2 magnitudes (depending on the band). At lower values of the SFD extinction, we find that the sky surface density of galaxies increases with increasing extinction, precisely the opposite of the effect expected from Galactic dust. We suggest that the far infrared (FIR) brightness of the sky in regions of true low dust extinction is significantly ``contaminated'' by the FIR emission from background galaxies. We show that such an explanation is both qualitatively and quantitatively consistent with the available data. Based on this interpretation we conclude that systematic errors in the SFD extinction map due to extragalactic FIR emission are quite small, of order hundredths of a magnitude, but nevertheless statistically detectable. (Abridged)

We present an algorithm for selecting an uniform sample of gravitationally lensed quasar candidates from low-redshift (0.6<z<2.2) quasars brighter than i=19.1 that have been spectroscopically identified in the SDSS. Our algorithm uses morphological and color selections that are intended to identify small- and large-separation lenses, respectively. Our selection algorithm only relies on parameters that the SDSS standard image processing pipeline generates, allowing easy and fast selection of lens candidates. The algorithm has been tested against simulated SDSS images, which adopt distributions of field and quasar parameters taken from the real SDSS data as input. Furthermore, we take differential reddening into account. We find that our selection algorithm is almost complete down to separations of 1'' and flux ratios of 10^-0.5. The algorithm selects both double and quadruple lenses. At a separation of 2'', doubles and quads are selected with similar completeness, and above (below) 2'' the selection of quads is better (worse) than for doubles. Our morphological selection identifies a non-negligible fraction of single quasars: To remove these we fit images of candidates with a model of two point sources and reject those with unusually small image separations and/or large magnitude differences between the two point sources. We estimate the efficiency of our selection algorithm to be at least 8% at image separations smaller than 2'', comparable to that of radio surveys. The efficiency declines as the image separation increases, because of larger contamination from stars. We also present the magnification factor of lensed images as a function of the image separation, which is needed for accurate computation of magnification bias.

We present results from Chandra observations of SDSS J1004+4112, a strongly lensed quasar system with a maximum image separation of 15". All four bright images of the quasar, as well as resolved X-ray emission originating from the lensing cluster, are clearly detected. The emission from the lensing cluster extends out to approximately 1.5 arcmin. We measure the bolometric X-ray luminosity and temperature of the lensing cluster to be 4.7e44 erg s^-1 and 6.4 keV, consistent with the luminosity-temperature relation for distant clusters. The mass estimated from the X-ray observation shows excellent agreement with the mass derived from gravitational lensing. The X-ray flux ratios of the quasar images differ markedly from the optical flux ratios, and the combined X-ray spectrum of the images possesses an unusually strong Fe Kalpha emission line, both of which are indicative of microlensing.

We report the results of a prototype photometric search for transiting extrasolar planets using Subaru Suprime-Cam. Out of about 100,000 stars monitored around the Galactic plane (l=90 degree,b=0 degree), we find that 7,700 (27,000) stars show the photometric precision below 1% (3%) for 60 second exposures which is required to detect extrasolar planets by the transit method. Thus Suprime-Cam has the photometric stability and accuracy required for a transiting planet survey. During this observing run, we detected three transiting planetary candidates (i'-band magnitude around 18.5) which exhibit a single full transit-like light curve with a fractional depth of <5%. While future photometric and/or spectroscopic follow-ups remain to be done, the estimated parameters for the three systems are consistent with planetary companions around main-sequence stars. We also found two eclipsing binary candidates and eleven variable stars exhibiting W UMa-like light curves.

An extensive photometric monitoring of KH 15D, an enigmatic variable in the young star cluster NGC 2264, has been conducted. Simultaneous and accurate near-infrared (JHKs-bands) photometry is presented between 2003 December and 2005 March covering most of the variable phase. The infrared variability is characterized by large-amplitude and long-lasting eclipse, as observed at optical. The period of variability is 48.3 +/- 0.2 days, the maximum photometric amplitude of variability is ~4.2 mag, and the eclipse duration is \~0.5 in phase units. These are consistent with the most recent period, amplitude, and duration at optical. The blueing of the J-H color (~0.16 mag) during the eclipse, which has been suggested before, is unambiguously confirmed; a similar blueing at H-Ks is less clear but is probably present at a similar level. The overall shape of the JHKs light curves is very similar to the optical one, including a fair time-symmetry and a less stable flux during the eclipse with a slight hump near the zero phase. Most of these variability features of KH 15D observed at near-infrared wavelengths can be explained with the recent model employing an eclipse by the inclined, precessing disk and an outer scattering region around a pre-main-sequence binary.

We report follow-up observations of two gravitational lens candidates identified in the Sloan Digital Sky Survey (SDSS) dataset. We have confirmed that SDSS J102111.02+491330.4 is a previously unknown gravitationally lensed quasar. This lens system exhibits two images of a $z = 1.72$ quasar, with an image separation of $1{\farcs}14 \pm 0.04$. Optical and near-IR imaging of the system reveals the presence of the lensing galaxy between the two quasar images. Observations of SDSS J112012.12+671116.0 indicate that it is more likely a binary quasar than a gravitational lens. This system has two quasars at a redshift of $z = 1.49$, with an angular separation of $1{\farcs}49 \pm 0.02$. However, the two quasars have markedly different SEDs and no lens galaxy is apparent in optical and near-IR images of this system. We also present a list of 31 SDSS lens candidates which follow-up observations have confirmed are \textit{not} gravitational lenses.

We determine the stellar, planetary, and orbital properties of the transiting planetary system HD 209458, through a joint analysis of high-precision radial velocities, photometry, and timing of the secondary eclipse. Of primary interest is the strong detection of the Rossiter-McLaughlin effect, the alteration of photospheric line profiles that occurs because the planet occults part of the rotating surface of the star. We develop a new technique for modeling this effect, and use it to determine the inclination of the planetary orbit relative to the apparent stellar equator [lambda = (-4.4 +/- 1.4) degrees], and the line-of-sight rotation speed of the star [v*sin(I) = (4.70 +/- 0.16) km/s]. The uncertainty in these quantities has been reduced by an order of magnitude relative to the pioneering measurements by Queloz and collaborators. The small but nonzero misalignment is probably a relic of the planet formation epoch, because the expected timescale for tidal coplanarization is larger than the age of the star. Our determination of v*sin(I) is a rare case in which rotational line broadening has been isolated from other broadening mechanisms.

We have searched for absorption in several common atomic species due to the atmosphere or exosphere of the transiting extrasolar planet HD 209458b, using high precision optical spectra obtained with the Subaru High Dispersion Spectrograph (HDS). Previously we reported an upper limit on H alpha absorption of 0.1% (3 sigma) within a 5.1\AA band. Using the same procedure, we now report upper limits on absorption due to the optical transitions of Na D, Li, H alpha, H beta, H gamma, Fe, and Ca. The 3 sigma upper limit for each transition is approximately 1% within a 0.3\AA band (the core of the line), and a few tenths of a per cent within a 2\AA band (the full line width). The wide-band results are close to the expected limit due to photon-counting (Poisson) statistics, although in the narrow-band case we have encountered unexplained systematic errors at a few times the Poisson level. These results are consistent with all previously reported detections and upper limits, but are significantly more sensitive.

Earth's deciduous plants have a sharp order-of-magnitude increase in leaf reflectance between approximately 700 and 750 nm wavelength. This strong reflectance of Earth's vegetation suggests that surface biosignatures with sharp spectral features might be detectable in the spectrum of scattered light from a spatially unresolved extrasolar terrestrial planet. We assess the potential of Earth's step-function-like spectroscopic feature, referred to as the "red edge", as a tool for astrobiology. We review the basic characteristics and physical origin of the red edge and summarize its use in astronomy: early spectroscopic efforts to search for vegetation on Mars and recent reports of detection of the red edge in the spectrum of Earthshine (i.e., the spatially integrated scattered light spectrum of Earth). We present Earthshine observations from Apache Point Observatory to emphasize that time variability is key to detecting weak surface biosignatures such as the vegetation red edge. We briefly discuss the evolutionary advantages of vegetation's red edge reflectance, and speculate that while extraterrestrial "light harvesting organisms" have no compelling reason to display the exact same red edge feature as terrestrial vegetation, they might have similar spectroscopic features at different wavelengths than terrestrial vegetation. This implies that future terrestrial-planet-characterizing space missions should obtain data that allow time-varying, sharp spectral features at unknown wavelengths to be identified. We caution that some mineral reflectance edges are similar in slope and strength to vegetation's red edge (albeit at different wavelengths); if an extrasolar planet reflectance edge is detected care must be taken with its interpretation.

There is evidence that the transiting planet HD 209458b has a large exosphere of neutral hydrogen, based on a 15% decrement in Lyman-alpha flux that was observed by Vidal-Madjar et al. during transits. Here we report upper limits on H-alpha absorption by the exosphere. The results are based on optical spectra of the parent star obtained with the Subaru High Dispersion Spectrograph. Comparison of the spectra taken inside and outside of transit reveals no exospheric H-alpha signal greater than 0.1% within a 5.1A band (chosen to have the same Delta_lambda/lambda as the 15% Ly-alpha absorption). The corresponding limit on the column density of n=2 neutral hydrogen is N_2 <~ 10^9 cm^{-2}. This limit constrains proposed models involving a hot (~10^4 K) and hydrodynamically escaping exosphere.

We analyze the properties of quasar variability using repeated SDSS imaging data in five UV-to-far red photometric bands, accurate to 0.02 mag, for 13,000 spectroscopically confirmed quasars. The observed time lags span the range from 3 hours to over 3 years, and constrain the quasar variability for rest-frame time lags of up to two years, and at rest-frame wavelengths from 1000 Ang. to 6000 Ang. We demonstrate that 66,000 SDSS measurements of magnitude differences can be described within the measurement noise by a simple function of only three free parameters. The addition of POSS data constrains the long-term behavior of quasar variability and provides evidence for a turn-over in the structure function. This turn-over indicates that the characteristic time scale for optical variability of quasars is of the order 1 year.

We present seven epochs of spectroscopy on the quadruply imaged quasar SDSS J1004+4112, spanning observed-frame time delays from 1 to 322 days. The spectra reveal differences in the emission lines between the lensed images. Specifically, component A showed a strong enhancement in the blue wings of several high-ionization lines relative to component B, which lasted at least 28 days (observed frame) then faded. Since the predicted time delay between A and B is <30 days, our time coverage suggests that the event was not intrinsic to the quasar. We attribute these variations to microlensing of part of the broad emission line region of the quasar, apparently resolving structure in the source plane on a scale of ~10^{16} cm at z=1.734. In addition, we observed smaller differences in the emission line profiles between components A and B that persisted throughout the time span, which may also be due to microlensing or millilensing. Further spectroscopic monitoring of this system holds considerable promise for resolving the structure of the broad emission line region in quasars.

We study the recently discovered gravitational lens SDSS J1004+4112, the first quasar lensed by a cluster of galaxies. It consists of four images with a maximum separation of 14.62''. The system has been confirmed as a lensed quasar at z=1.734 on the basis of deep imaging and spectroscopic follow-up observations. We present color-magnitude relations for galaxies near the lens plus spectroscopy of three central cluster members, which unambiguously confirm that a cluster at z=0.68 is responsible for the large image separation. We find a wide range of lens models consistent with the data, but they suggest four general conclusions: (1) the brightest cluster galaxy and the center of the cluster potential well appear to be offset by several kpc; (2) the cluster mass distribution must be elongated in the North--South direction, which is consistent with the observed distribution of cluster galaxies; (3) the inference of a large tidal shear (~0.2) suggests significant substructure in the cluster; and (4) enormous uncertainty in the predicted time delays between the images means that measuring the delays would greatly improve constraints on the models. We also compute the probability of such large separation lensing in the SDSS quasar sample, on the basis of the CDM model. The lack of large separation lenses in previous surveys and the discovery of one in SDSS together imply a mass fluctuation normalization \sigma_8=1.0^{+0.4}_{-0.2} (95% CL), if cluster dark matter halos have an inner slope -1.5. Shallower profiles would require higher values of \sigma_8. Although the statistical conclusion might be somewhat dependent on the degree of the complexity of the lens potential, the discovery is consistent with the predictions of the abundance of cluster-scale halos in the CDM scenario. (Abridged)

We report discovery of actively star forming elliptical galaxies in a morphologically classified sample of bright galaxies at a low redshift obtained from the Sloan Digital Sky Survey. The emission lines of these galaxies do not show the characteristics of active galactic nuclei, and thus their strong H$\alpha$ emission is ascribed to star formation with a rate nearly as high as that is seen in typical late spiral galaxies. This is taken as evidence against the traditional view that all elliptical galaxies formed early and now evolve only passively. The frequency of such star forming elliptical galaxies is a few tenths of a percent in the sample, but increases to 3% if we include active S0 galaxies. We may identify these galaxies as probable progenitors of so-called E+A galaxies that show the strong Balmer absorption feature of A stars superimposed on an old star population. The approximate match of the abundance of active elliptical plus S0 galaxies with that of E+A galaxies indicates that the duration of such late star formation episodes is of the order of $\gsim 1$ Gyr. If we interpret these galaxies as new additions to the early-type galaxy population, and if we assume a power law for their number evolution, the abundance of early-type galaxies at $z=1$ is about 30% less than that at $z=0$.

Gravitational lensing is a powerful tool for the study of the distribution of dark matter in the Universe. The cold-dark-matter model of the formation of large-scale structures predicts the existence of quasars gravitationally lensed by concentrations of dark matter so massive that the quasar images would be split by over 7 arcsec. Numerous searches for large-separation lensed quasars have, however, been unsuccessful. All of the roughly 70 lensed quasars known, including the first lensed quasar discovered, have smaller separations that can be explained in terms of galaxy-scale concentrations of baryonic matter. Although gravitationally lensed galaxies with large separations are known, quasars are more useful cosmological probes because of the simplicity of the resulting lens systems. Here we report the discovery of a lensed quasar, SDSS J1004+4112, which has a maximum separation between the components of 14.62 arcsec. Such a large separation means that the lensing object must be dominated by dark matter. Our results are fully consistent with theoretical expectations based on the cold-dark-matter model.

We report the discovery of SDSSJ115517.35+634622.0, a previously unknown gravitationally lensed quasar. The lens system exhibits two images of a $z = 2.89$ quasar, with an image separation of $1{\farcs}832 \pm 0.007$ . Near-IR imaging of the system reveals the presence of the lensing galaxy between the two quasar images. Based on absorption features seen in the Sloan Digital Sky Survey (SDSS) spectrum, we determine a lens galaxy redshift of $z = 0.1756$. The lens is rather unusual in that one of the quasar images is only $0{\farcs}22\pm0{\farcs}07$ ($\sim 0.1 R_{\rm eff}$) from the center of the lens galaxy and photometric modeling indicates that this image is significantly brighter than predicted by a SIS model. This system was discovered in the course of an ongoing search for strongly lensed quasars in the dataset from the SDSS.

Recently, Chartas et al. (2001) detected a rapid X-ray flare in the gravitationally lensed, multiple image quasar RX J0911.4+0551. Dramatic events, such as rapid X-ray flares, are useful in providing high precision measurements of the time delays between multiple images. In this paper, we argue that there is a new possibility in measurements of time delays between multiple images of gravitationally lensed quasars; constrain the locations of putative flares that give rise to the intrinsic rapid variabilities of quasars. The realization, however, of these goals cannot be presently achieved due to the limited accuracy of the current measurements. We predict that timing flares with accuracies of the order of a few seconds will be needed to probe the location of the flares. Our proposing method will work with better instruments in near future, such as XEUS.

We present an algorithm for selecting gravitational lens candidates from amongst Sloan Digital Sky Survey (SDSS) quasars. In median Early Data Release (EDR) conditions, the algorithm allows for the recovery of pairs of equal flux point sources down to separations of $\sim 0{\farcs}7$ or with flux ratios up to $\sim$ 10:1 at a separation of $1\farcs5$. The algorithm also recovers a wide variety of plausible quad geometries. We also present a method for determining the selection function of this algorithm through the use of simulated SDSS images and introduce a method for calibrating our simulated images through truth-testing with real SDSS data. Finally, we apply our algorithm and selection function to SDSS quasars from the EDR to get an upper bound on the lensing fraction. We find 13 candidates among 5120 z $>$ 0.6 SDSS quasars, implying an observed lensing fraction of not more than 4 $\times 10^{-3}$. There is one likely lens system in our final sample, implying an observed lensing fraction of not less than $3 \times 10^{-5}$ (95% confidence levels).

For the first time in human history the possibility of detecting and studying Earth-like planets is on the horizon. Terrestrial Planet Finder (TPF), with a launch date in the 2015 timeframe, is being planned by NASA to find and characterize planets in the habitable zones of nearby stars. The mission Darwin from ESA has similar goals. The motivation for both of these space missions is the detection and spectroscopic characterization of extrasolar terrestrial planet atmospheres. Of special interest are atmospheric biomarkers--such as O2, O3, H2O, CO and CH4--which are either indicative of life as we know it, essential to life, or can provide clues to a planet's habitability. A mission capable of measuring these spectral features would also obtain sufficient signal-to-noise to characterize other terrestrial planet properties. For example, physical characteristics such as temperature and planetary radius can be constrained from low- resolution spectra. In addition, planet characteristics such as weather, rotation rate, presence of large oceans or surface ice, and existence of seasons could be derived from photometric measurements of the planet's variability. We will review the potential to characterize terrestrial planets beyond their spectral signatures. We will also discuss the possibility to detect strong surface biomarkers--such as Earth's vegetation red edge near 700 nm--that are different from any known atomic, molecular, or mineralogical signature.

The light scattered by an extrasolar Earth-like planet's surface and atmosphere will vary in intensity and color as the planet rotates; the resulting light curve will contain information about the planet's properties. Since most of the light comes from a small fraction of the planet's surface, the temporal flux variability can be quite significant, $\sim$ 10-100%. In addition, for cloudless Earth-like extrasolar planet models, qualitative changes to the surface (such as ocean fraction, ice cover) significantly affect the light curve. Clouds dominate the temporal variability of the Earth but can be coherent over several days. In contrast to Earth's temporal variability, a uniformly, heavily clouded planet (e.g. Venus), would show almost no flux variability. We present light curves for an unresolved Earth and for Earth-like model planets calculated by changing the surface features. This work suggests that meteorological variability and the rotation period of an Earth-like planet could be derived from photometric observations. The inverse problem of deriving surface properties from a given light curve is complex and will require much more investigation.

Strong gravitational lensing of distant supernovae (SNe), particularly Type Ia's, has some exploitable properties not available when other sorts of cosmologically distant sources are lensed. One such property is that the ``standard candle'' nature of the SN at peak brightness allows a direct determination of the lensing magnification factor for each well observed image. Another is that the duration of an SN event is of the same order as the differential time delays between the various lens images. The comparable time scales of the photometric event and the time delay invite a variety of applications, including high precision measurements of the delay and the targeting of especially interesting phases of the explosion for intensive observation when they appear in trailing images. As an initial exploration of these possibilities we present calculations of SN lensing statistics in a ``concordance cosmology''. We emphasize magnification and time delay effects. Plausible SN surveys, such as the proposed SNAP space mission, would discover several to some tens of strongly lensed SNe Ia per year, and at least a few of these will be at redshifts well beyond those that would be accessible via unlensed events. The total number of such anomalously high redshift SNe Ia will be a useful test of high redshift star formation models. SN surveys of finite duration will, of course, miss the appearance of some images; we quantify this selection bias. Finally, we investigate how well the appearance of trailing images can be predicted based on various amounts of available information on the lensing event. Knowledge of the magnification factor for the leading image makes it possible to predict the appearance of a trailing image relatively accurately if the lens redshift is also known. (Abridged)

We present a data set of images of the gravitationally lensed quasar Q2237+0305, that was obtained at the Apache Point Observatory (APO) between June 1995 and January 1998. Although the images were taken under variable, often poor seeing conditions and with coarse pixel sampling, photometry is possible for the two brighter quasar images A and B with the help of exact quasar image positions from HST observations. We obtain a light curve with 73 data points for each of the images A and B. There is evidence for a long (>~ 100 day) brightness peak in image A in 1996 with an amplitude of about 0.4 to 0.5 mag (relative to 1995), which indicates that microlensing has been taking place in the lensing galaxy. Image B does not vary much over the course of the observation period. The long, smooth variation of the light curve is similar to the results from the OGLE monitoring of the system (Wozniak et al. 2000a).

The fraction of quasar's and gamma-ray burst (GRB) afterglows that vary due to microlensing by the stellar populations of intervening elliptical/S0 galaxies is computed by combining the joint distribution of effective microlensing convergence (kappa) and shear (gamma) with microlensing magnification patterns. Microlensing is common in multiply imaged sources. We find that 1 in 3 multiply imaged quasars should vary by more than 0.5 magnitudes per decade due to microlensing, while 10% of macrolensed GRB afterglows should show a departure of more than 0.5 magnitudes from their intrinsic light-curve during the first 30 days. However microlensing by stars is rare in general, with only 1 source in ~500 varying by more than 0.5 magnitudes during the same periods. We find that most microlensing by stars will be observed in a regime where gamma>0.1. Thus point-mass lenses do not provide an adequate description for most microlensing events. If dark matter halos contain a large fraction of mass in compact objects, the fraction of microlensed (by 0.5 magnitudes) images rises significantly to ~1 in 10 for quasars and ~1 in 5 for GRB afterglows. Comparison of variability between macrolensed and normal quasar images, and a moderate number of well sampled GRB afterglow light-curves should therefore discover or refute the existence of stellar mass compact objects in galaxy halos. While microlensing results in departures of the distribution of magnifications from that of a smooth model, the effect on the macrolensing magnification bias for the discovery of lenses in quasar surveys is small. On the other hand, microlensing significantly broadens the distribution of macrolensed image flux ratios.

Published parametric models of the Einstein Cross gravitational lens demonstrate that the image geometry can be reproduced by families of models. In particular, the slope of the mass-profile for the lens galaxy is unconstrained. However, recent models predict a dependence of image flux ratios on the slope of the mass profile. We use this dependence to constrain the mass profile by calculating the likelihood of the slope using published mid-IR flux ratios (including microlensing variability). We find that the galaxy is likely to be flatter than isothermal, and therefore that the mass-to-light ratio is decreasing in the inner kpc.

The distribution of differential time delays \Delta t between images produced by strong gravitational lensing contains information on the mass distributions in the lensing objects as well as on cosmological parameters such as H_0. We derive an explicit expression for the conditional probability distribution function of time delays P(\Delta t|\theta), given an image separation between multiple images \theta, and related statistics. We consider lensing halos described by the singular isothermal sphere (SIS) approximation and by its generalization as proposed by Navarro, Frenk, & White (NFW) which has a density profile \rho \propto r^{-\alpha} in the innermost region. The time delay distribution is very sensitive to these profiles; steeper inner slopes tend to produce larger time delays. For example, if H_0=70 km/s/Mpc, a \Lambda-dominated cosmology and a source redshift z_S=1.27 are assumed, lenses with \theta=5'' produce a time delay of \Delta t[yr]=1.5^{+1.7}_{-0.9}, 0.39^{+0.37}_{-0.22}, 0.15^{+0.11}_{-0.09}, and 0.071^{+0.054}_{-0.038} (50% confidence interval), for SIS, generalized NFW with \alpha=1.5, \alpha=1.0, and \alpha=0.5, respectively. At a fixed image separation, the time delay is determined by the difference in the lensing potential between the position of the two images, which typically occur at different impact parameters. Although the values of \Delta t are proportional to the inverse of H_0, P(\Delta t|\theta) is rather insensitive to all other cosmological model parameters, source redshifts, magnification biases and so on. A knowledge of P(\Delta t|\theta) will also be useful in designing the observing program of future large scale synoptic variability surveys and for evaluating possible selection biases operating against large splitting lens systems.

It is conventional to calculate the probability of microlensing for a cosmologically distant source based on the Press-Gunn approximation that the lensing objects are uniformly and randomly distributed in the intervening space with a constant comoving density. Here we investigate more realistic cosmological microlensing statistics by considering the strong spatial clustering of likely lensing objects with each other in galaxies and their association with the clumps of dark matter that make up the massive halos of galaxies. The distribution of microlensing optical depth (kappa) along randomly chosen sight lines is calculated as is the conditional distribution of kappa along sight lines near one which is strongly microlensed. Our overall result is that the Press-Gunn approximation is a useful order-of-magnitude approximation if the massive halos of galaxies are made of dark compact objects but that it fails badly and can be qualitatively misleading in the more likely case in which only the ordinary stellar populations of galaxies are the dominant source of cosmological microlensing events. In particular, we find that microlensing by stars is limited to of order 1 percent of high redshift sources at any one time. Furthermore, even though only a small fraction of high redshift sources are multiply-imaged (by galaxies), it is these sources that are most likely to be microlensed by stars. Consequently, microlensing by stars is usually observed at kappa's near 1 where the simple isolated point mass lens approximation is not appropriate. However, if CDM halos are composed of condensed objects, then more than 10 percent of high redshift sources are microlensed at any given time. The vast majority of these sources are not multiply-imaged, and have kappa's smaller than 0.01.

Deep mid-infrared (MIR) observations could reveal a population of faint, high-redshift (z>3) dusty starburst galaxies that are the progenitors of present-day spheroids or bulges, and are beyond the reach of current instruments. We utilize a semi-analytic galaxy formation scheme to find an extreme model for the MIR galaxy counts, designed to maximize the number of detectable sources down to a flux level of a few nJy. The model incorporates the formation of heavily dust-enshrouded stellar populations at high redshift, and is consistent with existing observations, including faint counts at 1.6um in the NICMOS Hubble Deep Field, and the upper limit on the extragalactic MIR background from TeV gamma rays. Our models predict upto 0.5 galaxies/sq.arcsec at the threshold of 100 nJy at 6um, with a comparable or larger surface density at longer MIR wavelengths. We conclude that a significant new population of high-redshift galaxies could be detected by the Space Infrared Telescope Facility (SIRTF) and Next Generation Space Telescope (NGST). Such a population would constitute background noise for the Terrestrial Planet Finder (TPF), and could necessitate repeat observations: every {\it TPF} resolution element have a 10 percent chance of being contaminated by a background galaxy.

The detection of massive planets orbiting nearby stars has become almost routine, but current techniques are as yet unable to detect terrestrial planets with masses comparable to the Earth's. Future space-based observatories to detect Earth-like planets are being planned. Terrestrial planets orbiting in the habitable zones of stars-where planetary surface conditions are compatible with the presence of liquid water-are of enormous interest because they might have global environments similar to Earth's and even harbor life. The light scattered by such a planet will vary in intensity and colour as the planet rotates; the resulting light curve will contain information about the planet's properties. Here we report a model that predicts features that should be discernible in light curves obtained by low-precision photometry. For extrasolar planets similar to Earth we expect daily flux variations up to hundreds of percent, depending sensitively on ice and cloud cover. Qualitative changes in surface or climate generate significant changes in the predicted light curves. This work suggests that the meteorological variability and the rotation period of an Earth-like planet could be derived from photometric observations. Other properties such as the composition of the surface (e.g., ocean versus land fraction), climate indicators (for example ice and cloud cover), and perhaps even signatures of Earth-like plant life could be constrained or possibly, with further study, even uniquely determined.

If it is hypothesised that there is no dark matter then some alternative gravitational theory must take the place of general relativity (GR) on the largest scales. Dynamical measurements can be used to investigate the nature of such a theory, but only where there is visible matter. Gravitational lensing is potentially a more powerful probe as it can be used to measure deflections far from the deflector, and, for sufficiently large separations, allow it to be treated as a point-mass. Microlensing within the local group does not yet provide any interesting constraints, as only images formed close to the deflectors are appreciably magnified, but stacking of multiple light-curves and observations of microlensing on cosmological scales may be able to discriminate between GR and non-dark matter theories. Galaxy-galaxy lensing is likely to be an even more powerful probe of gravity, with the Sloan Digital Sky Survey (SDSS) commissioning data used here to constrain the deflection law of galaxies to be A(R) propto R^(0.1 +/- 0.1) for impact parameters in the range 50 kpc < R < 1 Mpc. Together with observations of flat rotation curves, these results imply that, in any gravitational theory, photons must experience (close to) twice the deflection of massive particles moving at the speed of light (at least on these physical scales). The full SDSS data-set will also be sensitive to asymmetry in the lensing signal and variation of the deflection law with galaxy type. A detection of either of these effects would represent an independent confirmation that galaxies are dark matter-dominated; conversely azimuthal symmetry of the shear signal would rule out the typically ellipsoidal haloes predicted by most simulations of structure formation.

Modified Newtonian dynamics (MOND) is an alternative theory of gravity that aims to explain large-scale dynamics without recourse to any form of dark matter. However the theory is incomplete, lacking a relativistic counterpart, and so makes no definite predictions about gravitational lensing. The most obvious form that MONDian lensing might take is that photons experience twice the deflection of massive particles moving at the speed of light, as in general relativity (GR). In such a theory there is no general thin-lens approximation (although one can be made for spherically-symmetric deflectors), but the three-dimensional acceleration of photons is in the same direction as the relativistic acceleration would be. In regimes where the deflector can reasonably be approximated as a single point-mass (specifically low-optical depth microlensing and weak galaxy-galaxy lensing), this naive formulation is consistent with observations. Forthcoming galaxy-galaxy lensing data and the possibility of cosmological microlensing have the potential to distinguish unambiguously between GR and MOND. Some tests can also be performed with extended deflectors, for example by using surface brightness measurements of lens galaxies to model quasar lenses, although the breakdown of the thin-lens approximation allows an extra degree of freedom. Nonetheless, it seems unlikely that simple ellipsoidal galaxies can explain both constraints. Further, the low-density universe implied by MOND must be completely dominated by the cosmological constant (to fit microwave background observations), and such models are at odds with the low frequency of quasar lenses. These conflicts might be resolved by a fully consistent relativistic extension to MOND; the alternative is that MOND is not an accurate description of the universe.

Strong lensing is a powerful probe of the distribution of matter in the cores of clusters of galaxies. Recent studies suggest that the cold dark matter model predicts cores that are denser than those observed in galaxies, groups and clusters. One possible resolution of the discrepancy is that the dark matter has strong interactions (SIDM), which leads to lower central densities. A generalized form of the Navarro, Frenk and White profile (Zhao profile) may be used to describe these halos. In this paper we examine gravitational lensing statistics for this class of model. The optical depth to multiple imaging is a very sensitive function of the profile parameters in the range of interest for SIDM halos around clusters of galaxies. Less concentrated profiles, which result from larger self-interaction cross-sections, can produce many fewer lensed pairs. Lensing statistics provide a powerful test for SIDM. More realistic and observationally oriented calculations remain to be done, however larger self-interaction cross-sections may well be ruled out by the very existence of strong lenses on galaxy cluster scales. The inclusion of centrally dominant cluster galaxies should boost the cross-section to multiple imaging. However our preliminary calculations suggest that the additional multiple imaging rate is small with respect to the differences in multiple imaging rate for different halo profiles. In future statistical studies, it will be important to properly account for the scatter among halo profiles since the optical depth to multiple imaging is dominated by the most concentrated members of a cluster population.

Gravitational lensing is most often used as a tool to investigate the distribution of (dark) matter in the universe, but, if the mass distribution is known a priori, it becomes, at least in principle, a powerful probe of gravity itself. Lensing observations are a more powerful tool than dynamical measurements because they allow measurements of the gravitational field far away from visible matter. For example, modified Newtonian dynamics (MOND) has no relativistic extension, and so makes no firm lensing predictions, but galaxy-galaxy lensing data can be used to empirically the deflection law of a point-mass. MONDian lensing is consistent with general relativity, in so far as the deflection experienced by a photon is twice that experienced by a massive particle moving at the speed of light. With the deflection law in place and no invisible matter, MOND can be tested wherever lensing is observed. The implications are that either MONDian lensing is completely non-linear or that MOND is not an accurate description of the universe.

If the halo of the lensing galaxy 0957+561 is made of massive compact objects (MACHOs), they must affect the lightcurves of the quasar images Q0957+561 A and B differently. We search for this microlensing effect in the double quasar by comparing monitoring data for the two images A and B - obtained with the 3.5m Apache Point Observatory from 1995 to 1998 - with intensive numerical simulations. This way we test whether the halo of the lensing galaxy can be made of MACHOs of various masses. We can exclude a halo entirely made out of MACHOs with masses between 10e(-6) Msun and 10e(-2) Msun for quasar sizes of less than 3*10e(14)cm (h = 0.6), hereby extending previous limits upwards by one order of magnitude.

We report on a wide-area (48' x 9') imaging survey of faint galaxies in R and I bands toward the 1338+27 field where an unusual concentration of five QSOs at z\sim1.1, embedded in a larger-scale clustering of 23 QSOs, is known to exist. Using a quite homogeneous galaxy catalog with a detection completeness limit of I~23.5, we detect a significant clustering signature of faint red galaxies with I > 21 and R-I > 1.2 over a scale extending to ~20h^{-1}_{50} Mpc. Close examination of the color-magnitude diagram, the luminosity function, and the angular correlation function indeed suggests that those galaxies are located at z~1.1 and trace the underlying large-scale structure at that epoch, together with the group of 5 QSOs. Since the whole extent of the cluster of 23 QSOs (~70h^{-1}_{50} Mpc) is roughly similar to the local ``Great Wall'', the area may contain a high-redshift counterpart of superclusters in the local universe.

Gravitational lensing of a gamma ray burst (GRB) by a single point mass will produce a second, delayed signal. Several authors have discussed using microlensed GRBs to probe a possible cosmological population of compact objects. We analyse a closely related phenomenon; the effect of microlensing by low to medium optical depth in compact objects on the averaged observed light-curve of a sample of GRBs. We discuss the cumulative measured flux as a function of time resulting from delays due to microlensing by cosmological compact objects. The time-scale and curvature of this function describe unique values for the compact object mass and optical depth. For GRBs with durations larger than the detector resolution, limits could be placed on the mass and optical depth of cosmological compact objects. The method does not rely on the separation of lensed bursts from those which are spatially coincident.

The first investigations of the response of the microlensing magnification pattern (at an optical depth of order unity) to the mass function of the microlenses found that the resulting statistics depend only on the mean microlens mass <m>. In particular the mean microlensing caustic crossing rate was found to be proportional to the squareroot of <m>. We show that while this is true in the limit of mass functions with a narrow range of mass, in general the magnification pattern shows structure that reflects the contribution to the optical depth of microlenses with different masses. We present a better approximation, relating the microlens mass function to light-curve statistics. We show that the variability statistics of quasar microlensing light-curves can (in principle) be inverted to obtain the mass function of the microlenses in the mass range over which the mass density remains comparable, ie. p(m)dm~C/m. A preliminary analysis of the structure function for Q2237+0305 suggests that there is not a significant contribution to the optical depth from very low mass objects (0.001 solar mass). However observations of multiple microlensed quasars for a period of ~20 years may in the future yield a detailed p(m)dm. In the mass range where the number density is comparable, ie. p(m)dm~const., the distribution of flux factors could be inverted to find the microlens mass function. This may be used as a probe of the abundance of planets with orbital radii >100 AU.

Microlensing in Q2237+0305 between 1985 and 1995 (eg. Irwin et al. 1989; Corrigan et al. 1991; Ostensen et al. 1996) has been interpreted in two different ways; as microlensing by stellar mass objects of a continuum source having dimensions significantly smaller than the microlens Einstein radius (ER) (eg. Wambsganss, Paczynski & Schneider 1990; Rauch & Blandford 1991), and as microlensing by very low mass objects of a source as large as 5 ER (Refsdal & Stabell 1993; Haugan 1996). In this paper we present evidence in favour of a small source. Limits on the source size (in units of ER) are obtained from the combination of limits on the number of microlens Einstein radii crossed by the source during the monitoring period with two separate light-curve features. Firstly, recently published monitoring data (Wozniak et al. 2000; OGLE web page) show large variations (~0.8-1.5 magnitudes) between image brightnesses over a period of 700 days or ~15% of the monitoring period. Secondly, the 1988 peak in the image A light-curve had a duration that is a small fraction (<0.02) of the monitoring period. Such rapid microlensing rises and short microlensing peaks only occur for small sources. We find that the observed large-rapid variation limits the source size to be <0.2 ER (95% confidence). The width of the light-curve peak provides a stronger constraint of <0.02 ER (99% confidence). The Einstein radius (projected into the source plane) of the average microlens mass (m) in Q2237+0305 is ER ~ 10^{17}\sqrt{m} cm. The interpretation that stars are responsible for microlensing in Q2237+0305 therefore results in limits on the continuum source size that are consistent with current accretion disc theory.

Spectrophotometric observation of the gravitationally microlensed quasar Q2237+0305 during a High Magnification Event (HME) is potentially a very powerful tool for probing the structure of the quasars accretion disc on scales of less than 10^-8 arc seconds. How often we can expect to observe a HME is dependent on the lens system parameters of galactic transverse velocity, mean microlens mass and the size of the magnified continuum source. We have previously used published microlensed light-curves to obtain expressions for the likely-hood of the values for these parameters (Wyithe, Webster & Turner 1999b,2000; Wyithe, Webster, Turner & Mortlock 2000). Here we use this information to investigate the expected rate of SHMEs. We find the average rate of SHMEs as well as the number that we can expect to observe over periods of a decade and of a single observing season. We find that the average SHME rate summed over all images in Q2237+0305 is 1.5\pm0.6 - 6.3\pm1.3 events per decade. During the period following a caustic crossing we find that the event rate in the corresponding image is enhanced by 50-100%, and therefore that the overall event rate may be higher during these periods. From the distribution of events expected during a 6 month period we find that there is 1 chance in 4-10 of observing a SHME per observing season. The systematic dependence in these values arises from the different assumptions for smooth matter content, orientation of the galactic transverse velocity and the size of photometric error in the monitoring data. The results support continued monitoring of Q2237+0305 with the aim of obtaining detailed spectroscopic and photometric observations of a SHME.

The central regions of the gravitationally lensed quasar Q2237+0305 can be indirectly resolved on nano-arcsecond scales if viewed spectrophotometricly during a microlensing high magnification event (HME). Q2237+0305 is currently being monitored from the ground (eg. OGLE collaboration, Apache Point Observatory), with the goal, among others, of triggering ground and spacecraft based target of opportunity (TOO) observations of an HME. In this work we investigate the rate of change (trigger) in image brightness that signals an imminent HME and importantly, the separation between the trigger and the event peak. In addition, we produce colour dependent model light-curves by combining high-resolution microlensing simulations with a realistic model for a thermal accretion disc source. We make hypothetical target of opportunity spectroscopic observations using our determination of the appropriate trigger as a guide. We find that if the source spectrum varies with source radius, a 3 observation TOO program should be able to observe a microlensing change in the continuum slope following a light-curve trigger with a success rate of >80%.

The four bright images of the gravitationally lensed quasar Q2237+0305 are being monitored from the ground (eg. OGLE collaboration, Apache Point Observatory) in the hope of observing a high magnification event (HME). Over the past three seasons (1997-1999) the OGLE collaboration has produced microlensing light-curves with unprecedented coverage. These demonstrate smooth, independent (therefore microlensing) variability between the images (Wozniak et al. 2000a,b; OGLE web page). We have retrospectively compared probability functions for high-magnification event parameters with several observed light-curve features. We conclude that the 1999 image C peak was due to the source having passed outside of a cusp rather than to a caustic crossing. In addition, we find that the image C light-curve shows evidence for a caustic crossing between the 1997 and 1998 observing seasons involving the appearance of new critical images. Our models predict that the next image C event is most likely to arrive 500 days following the 1999 peak, but with a large uncertainty (100-2000 days). Finally, given the image A light-curve derivative at the end of the 1999 observing season, our modelling suggests that a caustic crossing will occur between the 1999 and 2000 observing seasons, implying a minimum for the image A light-curve ~1-1.5 magnitudes fainter than the November 1999 level.

Following the detection of a gravitational microlensing high magnification event (HME) in Q2237+0305A attempts have been made to place limits on the dimensions of the quasar continuum source. The analyses have studied either the observed event magnitude or the event duration. The latter approach has been hampered by lack of knowledge about the transverse velocity of the lensing galaxy. We obtain both upper and lower statistical limits on the size of the continuum source from the observed HME using determinations of transverse velocity obtained from the published monitoring data. Our calculations take account of the caustic orientation as well as the component of the caustic velocity that results from stellar proper motions. Our determination of source size relies on an estimated duration of 52 days for the HME, and so will be refined when more HMEs are observed. We find that the upper and lower limits on the magnified region of the R-band continuum source are 6 times 10^15 and 2 times 10^13 cm respectively (99% confidence). Through consideration of the joint probability for source size and mean microlens mass we find that the mean mass lies between ~0.01 solar masses and ~1 solar mass (95% confidence).

Gravitational microlensing at cosmological distances is potentially a powerful tool for probing the mass functions of stars and compact objects in other galaxies. In the case of multiply-imaged quasars, microlensing data has been used to determine the average microlens mass. However the measurements have relied on an assumed transverse velocity for the lensing galaxy. Since the measured mass scales with the square of the transverse velocity, published mass limits are quite uncertain. In the case of Q2237+0305 we have properly constrained this uncertainty. The distribution of light curve derivatives allows quantitative treatment of the relative rates of microlensing due to proper motions of microlenses, the orbital stream motion of microlenses and the bulk galactic transverse velocity. By demanding that the microlensing rate due to the motions of microlenses is the minimum that should be observed we determine lower limits for the average mass of stars and compact objects in the bulge of Q2237+0305. If microlenses are assumed to move in an orbital stream the lower limit ranges between 0.005 and 0.023 solar masses where the the systematic dependence is due to the fraction of smooth matter and the size of photometric error assumed for published monitoring data. However, if the microlenses are assumed to move according to an isotropic velocity dispersion then a larger lower limit of 0.019-0.11 solar masses is obtained. A significant contribution of Jupiter mass compact objects to the mass distribution of the galactic bulge of Q2237+0305 is therefore unambiguously ruled out.

We present a method for computing the probability distribution of microlensed light curve derivatives both in the case of a static lens with a transverse velocity, and in the case of microlensing that is produced through stellar proper motions. The distributions are closely related in form, and can be considered equivalent after appropriate scaling of the input transverse velocity. The comparison of the distributions in this manner provides a consistent way to consider the relative contribution to microlensing (both large and small fluctuations) of the two classes of motion, a problem that is otherwise an extremely expensive computational exercise. We find that the relative contribution of stellar proper motions to the microlensing rate is independent of the mass function assumed for the microlenses, but is a function of optical depth and shear. We find that stellar proper motions produce a higher overall microlensing rate than a transverse velocity of the same magnitude. This effect becomes more pronounced at higher optical depth. With the introduction of shear, the relative rates of microlensing become dependent on the direction of the transverse velocity. This may have important consequences in the case of quadruply lensed quasars such as Q2237+0305, where the alignment of the shear vector with the source trajectory varies between images.

We derive a simple empirical photometric redshift estimator using a training set of galaxies with multiband photometry and measured redshifts in the Hubble Deep Field (HDF). This estimator is model-independent; it does not use spectral templates. The dispersion between the estimated redshifts and the spectroscopically measured ones is small; the dispersions range from $\sigma_z\simeq 0.03$ to 0.1 for $z\la 2$ galaxies, and from $\sigma_z\simeq 0.14$ to 0.25 for $z\ga 2$ galaxies. The predictions provided by our empirical redshift estimator agree well with recently measured galaxy redshifts. We illustrate how our empirical redshift estimator can be modified to include flat spectrum galaxies with $1.4\la z \la 2$.

Determination of microlensing parameters in the gravitationally lensed quasar Q2237+0305 from the statistics of high magnification events will require monitoring for more than 100 years (Wambsganss, Paczynski & Schneider 1990). However we show that the effective transverse velocity of the lensing galaxy can be determined on a more realistic time-scale through consideration of the distribution of light-curve derivatives. The 10 years of existing monitoring data for Q2237+0305 are analysed. These data display strong evidence for microlensing that is not associated with a high magnification event. An upper limit of v < 500 km/sec is obtained for the galactic transverse velocity which is smaller than previously assumed values. The analysis suggests that the observed microlensing variation may be predominantly due to stellar proper motions. The statistical significance of the results obtained from our method will be increased by the addition of data points from current and future monitoring campaigns. However reduced photometric errors will be more valuable than an increased sampling rate.

The optical-ultraviolet continuum from active galactic nuclei (AGN) seem to originate from optically thick and/or thin disks, and occasionally from associated circumnuclear starburst regions. These different possible origins can, in principle, be discriminated by observations of gravitational microlensing event. We performed numerical simulations of microlensing of an AGN disk by a single lensing star passing in front of the AGN. Calculated spectral variations and light curves show distinct behavior, depending on the nature of emitting region; the time variation of a few months with strong wavelength dependence is expected in the case of an optically thick disk ({\it standard disk}), while an optically thin disk ({\it advection-dominated disk}) will produce shorter, nearly wavelength-independent variation. In the case of a starburst region much slower variations (over a year) will be superposed on the shorter variations.

We analyze observations of the quadruple lensed quasar Q2237+0305, taken with the WFPC2 camera in the F336W and F300W bands. 25 exposures were performed within 15 hours real time on 3 November 1995. On a timescale of 3--4 hours, we observe no variation in component A of greater than 0.02 mag. The other components are constant over a period of 10 hours to within about 0.05 mag. In the final 5 hours there is some evidence (not conclusive) for variation of component D by about 0.1 mag. Component A was brighter than component B by about 0.3 mag. Components C and D were fainter than component A by about 1.3 and 1.4 mag. Any fifth (central) component was at least 6.5 mag fainter than component A. Using the PC chip, we measure the relative distances of the four components with high accuracy. Our values are systematically larger than other investigators' (by 0.1% to 2.0%). The F336W filter had been chosen for the observations because it could have allowed us to see extended Ly-alpha emission from the Broad-Line Region (BLR). However, the quasar components are consistent with a point source. We conclude that there cannot be a Ly-alpha feature in the image plane brighter than about 23.5 mag in F336W and further from the quasar core than 100 mas. According to a lensing model by Rix, Scheider & Bahcall (1992), this would preclude any such features in the source plane further than 20 mas (~ 100 h^{-1} pc, assuming q_0 = 0.5) from the quasar core and brighter than 25 mag before magnification.

We present an R-band and J-band photometry of an optical transient which is likely to be associated with the gamma-ray burst event GRB 971214. Our first measurement took place 13 hours after the gamma-ray event. The brightness decayed with a power-law exponent alpha = -1.20 +- 0.02, which is similar to those of GRB 970228 and GRB 970508 which had exponents of alpha = -1.10 +- 0.04 and alpha = -1.141 +- 0.014 respectively. The transient decayed monotonically during the first four days following the gamma-ray event in contrast with the optical transient associated with GRB 970508 which increased in brightness, peaking two days after the burst, before settling to a power-law decay.

We derive simple empirical color-redshift relations for $z\la 4$ galaxies in the Hubble Deep Field (HDF) using a linear function of three photometric colors ($U-B$, $B-V$, $V-I$). The dispersion between the estimated redshifts and the spectroscopically observed ones is small for relations derived in several separate color regimes; the dispersions range from $\sigma_z\simeq 0.03$ to 0.1 for $z\la 2$ galaxies, and from $\sigma_z\simeq 0.14$ to 0.25 for $z\ga 2 $ galaxies. We apply the color-redshift relations to the HDF photometric catalog and obtain estimated redshifts that are consistent with those derived from spectral template fitting methods. The advantage of these color-redshift relations is that they are simple and easy to use and do not depend on the assumption of any particular spectral templates; they provide model independent redshift estimates for $z\la 4$ galaxies using only multi-band photometry, and they apply to about 90% of all galaxies. We provide a color-based estimated redshift catalog of HDF galaxies to $z\la 4$. We use the estimated redshifts to investigate the redshift distribution of galaxies in the HDF; we find peaks in the redshift distribution that suggest large-scale clustering of galaxies to at least $z\sim 1$ and that are consistent with those identified in spectroscopic probes of the HDF.

The innermost regions of quasars can be resolved by a gravitational-lens {\lq}telescope{\rq} on scales down to a few AU. For the purpose, X-ray observations are most preferable, because X-rays originating from the innermost regions, can be selectively amplified by microlensing due to the so-called `caustic crossing'. If detected, X-ray variations will constrain the size of the X-ray emitting region down to a few AU. The maximum attainable resolution depends mainly on the monitoring intervals of lens events, which should be much shorter than the crossing time. On the basis of this idea, we performe numerical simulations of microlensing of an optically-thick, standard-type disk as well as an optically-thin, advection-dominated accretion flow (ADAF). Calculated spectral variations and light curves show distinct behaviors, depending on the photon energy. X-ray radiation which is produced in optically thin region, exhibits intensity variation over a few tens of days. In contrast, optical-UV fluxes, which are likely to come from optically thick region, exhibit more gradual light changes, which is consistent with the microlensing events so far observed in Q2237+0305. Currently, Q2237+0305 is being monitored in the optical range at Apache Point Observatory. Simultaneous multi-wavelength observations by X-ray sattelites (e.g., ASCA, AXAF, XMM) as well as HST at the moment of a microlens event enable us to reveal an AU scale structure of the central accretion disk around black hole.

Aperiodic optical variability is a common property of Active Galactic Nuclei (AGNs), though its physical origin is still open to question. To study the origin of the optical -- ultraviolet variability in AGN, we compare light curves of two models to observations of quasar 0957+561 in terms of a structure function analysis. In the starburst (SB) model, random superposition of supernovae in the nuclear starburst region produce aperiodic luminosity variations, while in the disk-instability (DI) model, variability is caused by instabilities in the accretion disk around a supermassive black hole. We calculate fluctuating light curves and structure functions, $V(\tau)$, by simple Monte-Carlo simulations on the basis of the two models. Each resultant $V(\tau)$ possesses a power-law portion, $[V(\tau)]^{1/2} \propto \tau^{\beta}$, at short time lags ($\tau$). The two models can be distinguished by the logarithmic slope, $\beta$; $\beta \sim$ 0.74--0.90 in the SB model and $\beta \sim$ 0.41--0.49 in the DI model, while the observed light curves exhibit $\beta \sim$ 0.35. Therefore, we conclude that the DI model is favored over the SB model to explain the slopes of the observational structure function, in the case of 0957+561, though this object is a radio-loud object and thus not really a fair test for the SB model. In addition, we examine the time-asymmetry of the light curves by calculating $V(\tau)$ separately for brightening and decaying phases. The two models exhibit opposite trends of time-asymmetry to some extent, although the present observation is not long enough to test this prediction.

We analyze bright-end (K = 10 -17) galaxy counts from a number of near-infrared galaxy surveys. All studies available as of mid-1997, considered individually or collectively, show that the observed near-infrared galaxy number counts at low magnitudes are inconsistent with a simple no-evolution model. We examine evolutionary effects and a local underdensity model as possible causes of this effect. We find that the data are fit by either a factor of 1.7 - 2.4 deficiency of galaxies out to a redshift of z = 0.10 - 0.23, depending on the k corrections and evolution (e-) corrections used and the adopted values of the Schechter luminosity function parameters, or by unexpectedly strong low redshift evolution in the K-band, leading to (e+k)-corrections at z= 0.5 that are as much as 60% larger than accepted values. The former possibility would imply that the local expansion rate on scales of several hundred Mpc exceeds the global value of H_0 by up to 30% and that the amplitude of very large scale density fluctuations is far larger than expected in any current cosmogonic scenario. The latter possibility would mean that even the apparently most secure aspects of our understanding of galaxy evolution and spectral energy distributions are seriously flawed.

Flux estimates for faint sources or transients are systematically biased high because there are far more truly faint sources than bright. Corrections which account for this effect are presented as a function of signal-to-noise ratio and the (true) slope of the faint-source number-flux relation. The corrections depend on the source being originally identified in the image in which it is being photometered. If a source has been identified in other data, the corrections are different; a prescription for calculating the corrections is presented. Implications of these corrections for analyses of surveys are discussed; the most important is that sources identified at signal-to-noise ratios of four or less are practically useless.

Low amplitude (linear regime) cosmic density fluctuations lead to spatial variations in the locally measurable value of $H_0$ (denoted as $H_L$), $\delta_H \equiv (H_L-H_0)/H_0$, which are of order 3-6% (95% confidence interval) in a sphere of 200$ h^{-1}$Mpc in diameter, and of order 1-2% in a sphere of 400$ h^{-1}$Mpc in diameter, for three currently viable structure formation models (tilted CDM, $\Lambda$CDM, and open CDM) as normalized by the 4 year COBE DMR data. However, the true matter distribution power spectrum may differ from any of the currently viable models. For example, it may contain sharp features which have escaped detection so far. The measured CMB dipole velocity (the Galaxy's peculiar velocity with respect to the CMB rest frame) provides additional constraints on the probability distribution of $\delta_H$ that supplement our limited knowledge of the power spectrum. For a matter power spectrum which consists of the smooth power spectrum of a viable cosmological model plus a delta-function bump, we find that given the CMB dipole velocity, the 95% CL upper limit of $|\delta_H|$ increases approximately by a factor of two, but the probability distribution of $\delta_H$ is non-Gaussian, with increased probability at small $\delta_H$ compared to Gaussian. Abandoning model power spectra entirely, we find that the observed CMB dipole velocity alone provides a very robust limit, $\sqrt{< \delta_H^2 >_R} <10.5 h^{-1}{Mpc}/R$ at 95% CL, in a sphere of radius $R$, for an arbitrary power spectrum. Thus, variations between currently available local measures of $H_0$ and its true global value of a few to several percent are to be expected and differences as large as 10% are possible based on our current knowledge of the CMB anisotropies.

We study strong gravitational lensing by spiral galaxies, modeling them as infinitely thin uniform disks embedded in singular isothermal spheres. We derive general properties of the critical curves and caustics analytically. The multiple-image cross section is a sensitive function of the inclination angle of the disk relative to the observer. We compute the inclination averaged cross section for several sets of lensing parameters. For realistic disk mass and size parameters, we find that the cross section for multiple imaging is increased by only a modest factor and no dramatic increase in the optical depth for strong lensing of QSOs would be expected. However, the cross section for high magnifications is significantly increased due to the inclusion of a disk, especially for nearly edge-on configurations; due to the strong observational selection effects favoring high magnifications, there might be significant consequences for lensing statistics.

Every astrophysical object (dark or not) is a gravitational lens, as well as a receiver/observer of the light from sources lensed by other objects in its neighborhood. For a given pair of source and lens, there is a thin on-axis tubelike volume behind the lens in which the radiation flux from the source is greatly increased due to gravitational lensing. Any objects which pass through such a thin tube or beam will experience strong bursts of radiation, i.e., Extreme Gravitational Lensing Events (EGLEs). We have studied the physics and statistics of EGLEs. EGLEs may have interesting astrophysical effects, such as the destruction of dust grains, ignition of masers, etc. Here we illustrate the possible astrophysical effects of EGLEs with one specific example, the destruction of dust grains in globular clusters.

Interferometric gravitational wave detectors could measure the frequency sweep of a binary inspiral [characterized by its chirp mass] to high accuracy. The observed chirp mass is the intrinsic chirp mass of the binary source multiplied by $(1+z)$, where $z$ is the redshift of the source. Assuming a non-zero cosmological constant, we compute the expected redshift distribution of observed events for an advanced LIGO detector. We find that the redshift distribution has a robust and sizable dependence on the cosmological constant; the data from advanced LIGO detectors could provide an independent measurement of the cosmological constant.

Photometric monitoring of the gravitational lens system 0957+561A,B in the g and r bands with the Apache Point Observatory (APO) 3.5 m telescope during 1996 shows a sharp g band event in the trailing (B) image light curve at the precise time predicted from the observation of an event during 1995 in the leading (A) image with a delay of 415 days. This success confirms the "short delay," and the lack of any feature at a delay near 540 days rejects the "long delay" for this system, resolving a long-standing controversy. A series of statistical analyses of our light curve data yield a best fit delay of 417 +/- 3 days (95% confidence interval). Recent improvements in the modeling of the lens system (consisting of a galaxy and cluster) allow us to derive a value of the global (at z = 0.36) value of Hubble's constant H_0 using Refsdal's method, a simple and direct distance determination based on securely understood physics and geometry. The result is H_0 = 63 +/- 12 km/s/Mpc (for Omega = 1) where this 95% confidence interval is dominated by remaining lens model uncertainties.

Foreground galaxies that amplify the light from background quasars may also dim that light if the galaxies contain enough dust. Extinction by dust in lenses could hide the large number of lensed systems predicted for a flat universe with a large value of the cosmological constant $\Lambda$. We look for one signature of dust, namely reddening, by examining optical-infrared colors of gravitationally lensed images of quasars. We find that the lensed systems identified in radio and infrared searches have redder optical-IR colors than optically selected ones. This could be due to a bias against selecting reddened (hence extincted) quasars in the optical surveys, or due to the differences in the intrinsic colors of optical and radio quasars. Comparison of the radio-selected lensed and unlensed quasars shows that the lensed ones have redder colors. We therefore conclude that at least part of the color difference between the two lens samples is due to dust. From the color difference between lensed and unlensed radio quasars (and assuming Galactic extinction law) we can reconcile a large cosmological constant ($\Lambda=0.9$) with the number of lensed systems observed in flux limited optical surveys. These results substantially weaken the strongest constraint on cosmological scenarios that invoke a non-zero cosmological constant to explain age discrepancy problems, satisfy predictions of inflationary models of the early universe and play a role in large scale structure formation models. They also raise the prospect of using gravitational lenses to study the interstellar medium in high redshift galaxies.

We consider an astrophysical system with a population of sources and a population of lenses. For each pair of source and lens, there is a thin on-axis tube-like volume behind the lens in which the radiation flux from the source is greatly increased due to gravitational lensing. Any objects (such as dust grains) which pass through such a thin tube will experience strong bursts of radiation, i.e., Extreme Gravitational Lensing Events (EGLE). We study the physics and statistics of EGLE for the case in which the shear is larger or comparable to the finite source size. We find that the presence of shear has only a small effect on the EGLE statistics.

For a given source and lens pair, there is a thin on-axis tube-like volume behind the lens in which the radiation flux from the source is greatly increased due to gravitational lensing. Any objects (such as dust grains) which pass through such a thin tube will experience strong bursts of radiation, i.e., Extreme Gravitational Lensing Events (EGLE). We study the physics and statistics of EGLE for the case in which finite source size is more important than shear. One of the several possible significant astrophysical effects is investigated with an illustrative calculation.

We discuss the gravitational lensing of gravitational waves from merging neutron star binaries, in the context of advanced LIGO type gravitational wave detectors. We consider properties of the expected observational data with cut on the signal-to-noise ratio $\rho$, i.e., $\rho>\rho_0$. An advanced LIGO should see unlensed inspiral events with a redshift distribution with cut-off at a redshift $z_{\rm max} < 1$ for $h \leq 0.8$. Any inspiral events detected at $z>z_{\rm max}$ should be lensed. We compute the expected total number of events which are present due to gravitational lensing and their redshift distribution for an advanced LIGO in a flat Universe. If the matter fraction in compact lenses is close to 10\%, an advanced LIGO should see a few strongly lensed events per year with $\rho >5$.

We propose a new method to measure the mass of large-scale filaments in galaxy redshift surveys. The method is based on the fact that the mass per unit length of isothermal filaments depends only on their transverse velocity dispersion. Filaments that lie perpendicular to the line of sight may therefore have their mass per unit length measured from their thickness in redshift space. We present preliminary tests of the method and find that it predicts the mass per unit length of filaments in an N-body simulation to an accuracy of ~35%. Applying the method to a select region of the Perseus-Pisces supercluster yields a mass-to-light ratio M/L_B around 460h in solar units to within a factor of two. The method measures the mass-to-light ratio on length scales of up to 50h^(-1) Mpc and could thereby yield new information on the behavior of the dark matter on mass scales well beyond that of clusters of galaxies.

A unique reconstruction of the image of a high redshift source galaxy responsible for multiple long arcs in the z = 0.4 cluster 0024+1654 is obtained by inverse lensing. Deep B and I imaging with the Hubble Space Telescope (Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science Institute, which is operated by AURA, under NASA contract NAS 5-26555) enables high resolution of the arcs due to strong gravitational lensing of the background source. The gravitational lens in the foreground cluster is thus used to obtain a magnified view of the distant source. Four strongly lensed images of the source lead to a unique reconstruction. Each of the long arcs, when unlensed, leads to the same reconstructed source image exhibiting a beaded ring-like morphology. The U luminosity of the ring alone is equivalent to a normal galaxy. This is likely a galaxy in formation.

CCD photometry of the gravitational lens system 0957+561A,B in the g and r bands was obtained on alternate nights, weather permitting, from December 1994 through May 1995 using the Double Imaging Spectrograph (DIS) on the Apache Point Observatory (APO) 3.5-meter telescope. The remote observing and fast instrument change capabilities of this facility allowed accumulation of light curves sampled frequently and consistently. The Honeycutt ensemble photometry algorithm was applied to the data set and yielded typical relative photometric errors of approximately 0.01 magnitudes. Image A exhibited a sharp drop of about 0.1 magnitudes in late December 1994; no other strong features were recorded in either image. This event displays none of the expected generic features of a microlensing-induced flux variation and is likely to be intrinsic to the quasar; if so, it should also be seen in the B image with the lensing differential time delay. We give the expected 1996 image B light curves based on two values of the time delay and brightness ratio which have been proposed and debated in the literature. Continued monitoring of the system in the first half of 1996 should easily detect the image B event and thus resolve the time-delay controversy.

In the now well established conventional view (see Rees [1] and references therein), quasi-stellar objects (QSOs) and related active galactic nuclei (AGN) phenomena are explained as the result of accretion of plasma onto giant black holes which are postulated to form via gravitational collapse of the high density regions in the centers of massive host galaxies. This model is supported by a wide variety of indirect evidence and seems quite likely to apply at least to some observed AGN phenomena. However, one surprising set of new Hubble Space Telescope (HST) observations [2-4] directly challenges the conventional model, and the well known evolution of the QSO population raises some additional, though not widely recognized, difficulties. We propose here an alternative possibility: the Universe contains a substantial independent population of super-massive black holes, and QSO's are a phenomenon that occurs due to their collisions with galaxies or gas clouds in the intergalactic medium (IGM). This hypothesis would naturally explain why the QSO population declines very rapidly towards low redshift, as well as the new HST data.

Most physical properties derived for quasars, as single entities or as a population, depend upon the cosmology assumed. In this paper, we calculate the quasar luminosity function and some related quantities for a flat universe dominated by a cosmological constant $\Lambda$ ($\Lambda=0.9, \Omega=0.1$) and compare them with those deduced for a flat universe with zero cosmological constant ($\Lambda=0,\Omega=1$). We use the AAT quasar survey data (Boyle et al. 1990) as input in both cases. The data are fit well by a pure luminosity evolution model for both the cosmologies, but with different evolutionary parameters. From the luminosity function, we predict (extrapolate) a greater number of quasars at faint apparent magnitudes (twice the number at B=24, $z < 2.2$) for the $\Lambda$ dominated universe. This population of faint quasars at high redshift would result in a higher incidence of gravitational lensing. The total luminosity of the quasar population and the total mass tied up in black hole remnants of quasars is not sensitive to the cosmology. However, for a $\Lambda$ cosmology this mass is tied up in fewer but more massive black holes.

Gravitational lensing provides a strict test of cosmogonic models. Detailed numerical propagation of light rays through a universe having a distribution of inhomogeneities derived from the standard CDM (cold dark matter) scenario with the aid of massive, fully nonlinear computer simulations is used to test the model. It predicts that we should have seen far more widely split quasar images than have been found. These and other inconsistencies rule out the COBE normalized CDM model with $\Omega=1$ and $H_0=50$km/s/Mpc; but variants of this model might be constructed, which could pass the stringent tests provided by strong gravitational lensing.

We investigate the impact of hierarchical galaxy merging on the statistics of gravitational lensing of distant sources. Since no definite theoretical predictions for the merging history of luminous galaxies exist, we adopt a parametrized prescription, which allows us to adjust the expected number of pieces comprising a typical present galaxy at z=0.65. The existence of global parameter relations for elliptical galaxies and constraints on the evolution of the phase space density in dissipationless mergers, allow us to limit the possible evolution of galaxy lens properties under merging. We draw two lessons from implementing this lens evolution into statistical lens calculations: (1) The total optical depth to multiple imaging (e.g. of quasars) is quite insensitive to merging. (2) Merging leads to a smaller mean separation of observed multiple images. Because merging does not reduce drastically the expected lensing frequency it cannot make $\lambda$-dominated cosmologies compatible with the existing lensing observations. A comparison with the data from the HST Snapshot Survey shows that models with little or no evolution of the lens population are statistically favored over strong merging scenarios. The specific merging scenario proposed by Toomre (1977) can be rejected (95\% level) by such a comparison. Some versions of the scenario proposed by Broadhurst, Ellis \& Glazebrook (1992), are statistically acceptable.

The evolution of nonlinear density fluctuations around the Jeans mass shortly after cosmological recombination is analyzed using a 3D hydrodynamics/dark--matter code. The Cosmic Background Radiation (CBR) exerts Compton friction on free electrons due to peculiar velocities. The dynamics therefore depends strongly on the gas ionization history. Under a variety of ionization conditions and in systems with or without non-baryonic components, the baryons lose angular momentum efficiently and collapse to form a compact optically--thick object which would probably quickly evolve into a massive black hole. Attention is concentrated on elucidating some of the novel physical effects in early cosmological collapses, but ways in which more realistic calculations might be made and in which the scenario could be incorporated into a more complete cosmogonic model are discussed.