Strong quasar-galaxy lensing provides a powerful tool to probe the inter-stellar medium (ISM) of the lens galaxy using radiation from the background quasar. Using the Cosmic Origins Spectrograph (COS) on board the Hubble Space Telescope, we study the cold ISM properties of the lens galaxy in B1152+199 at a redshift of z=0.4377. Since existing optical extinction and X-ray absorption measurements of the lens have revealed a large amount of cold ISM, we expected to detect a damped Lya absorption (DLA) system in the near ultraviolet spectrum; however, our upper limit on the HI column density is several orders of magnitude below the expectation. We also marginally detect OI and CII absorption lines associated with the lens galaxy in the COS spectrum. Thus, the lens galaxy is identified as a ghostly DLA system, and further investigations of these ghostly DLA systems would be important to characterize the biases of using DLAs to probe the matter density of the universe. Although preliminary, the most likely explanation of the non-detection of the DLA is because of the Lya emission of the lens galaxy that fills in the absorption trough, with a Lya luminosity of 4e42 erg/s.

Here we confirm the short periodic variability of a subparsec supermassive binary black hole (SMBBH) candidate Mrk 231 in the extended optical photometric data set collected by the Catalina Real-Time Transient Survey (CRTS) and All-Sky Automated Survey for Supernovae (ASAS-SN). Using the Lomb-Scargle periodogram and 2DHybrid method, we detected the significant periodicity of ~ 1.1 yr beyond a damped random walk model in the CRTS+ASAS-SN optical data set. Mrk 231 has been previously proposed as a SMBBH candidate with a highly unequal mass ratio (q~ 0.03), very tight mutual separation of ~590 AU, and an orbital period of ~1.2 yr. Hence, our result further supports, even though not prove, the intriguing hypothesis that SMBBHs with low mass ratios may be more common than close-equal mass SMBBHs. This result, however, was obtained from the contribution of CRTS data with limited sampling cadence and photometric accuracy, and further monitoring of Mrk 231 is crucial to confirm the periodicity.

We show in detail that the recently derived expression for evaluating the integrated Sachs-Wolfe (ISW) temperature shift in the cosmic microwave background (CMB) caused by individual embedded (compensated) lenses is equivalent to the conventional approach for flat background cosmologies. The conventional approach requires evaluating an integral of the time derivative of the lensing potential, whereas the new Fermat potential approach is simpler and only requires taking a derivative of the potential part of the time delay.

Most galaxies in clusters have supermassive black holes at their center, and a fraction of those supermassive black holes show strong activity. These active galactic nuclei(AGNs) are an important probe of environmental dependence of galaxy evolution, intra-cluster medium, and cluster-scale feedback. We investigated AGN fraction in one of the largest samples of X-ray selected clusters from the ROSAT and their immediate surrounding field regions below z < 0.5. We found lower average AGN fraction in clusters, (2.37+-0.39)% than for the fields (5.12+-0.16)%. The lower AGN fractions in clusters were measured, after dividing the clusters into five redshift intervals between 0.0 and 0.5, in each redshift interval, and we found an increase in the fraction for both cluster and field galaxies with redshift below z < 0.5, which clearly indicates an environment and redshift dependence. We further divided the clusters into low-mass and high-mass objects using a mass cut at log(M500/Msun) = 13.5, finding comparable AGN fractions for both classifications, while a significantly higher AGN fraction in field. We also measured increasing AGN fractions with clustercentric distance for all redshift bins, further confirming the environmental dependence of AGN activities. In addition, we did not find an obvious trend between AGN fraction and SDSS-R absolute magnitudes among different redshift bins. We conclude that the lower AGN fraction in clusters relative to fields indicate that factors, such as inefficient galaxy mergers and ram pressure stripping cause a deficit of cold gas available in high density regions to fuel the central super-massive black hole. Clusters and fields in present universe have lost more gas relative to their high redshift counterparts resulting in a lower AGN fraction observed today.

We construct a sample of 10,680 wall galaxies and 3,064 void galaxies with $M_R \lesssim -20$ by cross referencing a void catalog from literature with Baryon Oscillation Spectroscopic Survey (BOSS) CMASS and WiggleZ galaxies, where the CMASS survey targets redder galaxies and the WiggleZ survey targets bluer galaxies. Comparing the density profiles of the red and blue galaxies as a function of the void radius, we find that the number ratio of red-to-blue galaxies increases with distances from the void centers, suggesting a deficit of luminous and normal red galaxies in voids. We find a mean (g -- r) magnitude color of 1.298 and 1.210 for the wall and void galaxies, respectively, when considering the combined red and blue sample, which is found to be a significant difference. However, when considering the blue and red samples separately, we find no significant color difference. We conclude that the constituents galaxies of each population, rather than intrinsic color difference, is the main driver in the apparent average color difference of galaxies in voids and walls, indicating a deficit of luminous and normal red galaxies in voids. Our analysis suggests that the primary environmental-dependence effect on galaxy evolution for normal and luminous galaxies between void and wall regions is manifested in the number of red galaxies, which depends on the environmental-dependent merger history. Using a semi-analytic simulation model, we can successfully reproduce the apparent color difference between the void and wall galaxies.

We present a new method to reconstruct three dimensional electron density of Coronal Mass Ejection (CME) based on genetic algorithm, namely genetic reconstruction method (GRM). GRM is firstly applied to the model CMEs with different orientations and shapes. A set of analytic GL98 model CMEs are employed to produce synthetic CME images for GRM reconstruction. Model CMEs with longitude of 0, 45, 90, 135, 180 degrees and latitude of 0, 15, 30, 45 degrees are used to test performance of GRM. We also obscure the CME with simulated occulter of coronagraph to find out the influence of incompleteness of CME brightness. GRM reconstruction is carried out using data from STEREO A and B with separation angle of 90 degrees and from STEREO A and SOHO with separation angle of 73 degrees, respectively. Pearson correlation coefficient (pcc) is calculated to analyze similarity between model and reconstructed CMEs for brightness and electron density. Comparisons based on the similarity analysis under various conditions stated above give us valuable insights of advantages and limitations of GRM reconstruction. Then the method is applied to real coronagraph data from STEREO-A, B and SOHO on September 30th, 2013.

Quasar microlensing serves as a unique probe of discrete objects within galaxies and galaxy clusters. Recent advancement of the technique shows that it can constrain planet-scale objects beyond our native galaxy by studying their induced microlensing signatures, the energy shift of emission lines originated in the vicinity of the black hole of high redshift background quasars. We employ this technique to exert effective constraints on the planet-mass object distribution within two additional lens systems, Q J0158$-$4325 ($z_l = 0.317$) and SDSS J1004+4112 ($z_l = 0.68$) using Chandra observations of the two gravitationally-lensed quasars. The observed variations of the emission line peak energy can be explained as microlensing of the FeK$\alpha$ emission region induced by planet-mass microlenses. To corroborate this, we perform microlensing simulations to determine the probability of a caustic transiting the source region and compare this with the observed line shift rates. Our analysis yields constraints on the sub-stellar population, with masses ranging from Moon ($10^{-8} M_{\odot}$) to Jupiter ($10^{-3} M_{\odot}$) sized bodies, within these galaxy or cluster scale structures, with total mass fractions of $\sim 3\times10^{-4}$ and $\sim 1\times10^{-4}$ with respect to halo mass for Q J0158$-$4325 and SDSS J1004+4112, respectively. Our analysis suggests that unbound planet-mass objects are universal in galaxies, and we surmise the objects to be either free-floating planets or primordial black holes. We present the first-ever constraints on the sub-stellar mass distribution in the intra-cluster light of a galaxy cluster. Our analysis yields the most stringent limit for primordial black holes at the mass range.

In this work, we calculated the sizes of unresolved X-ray emission regions in three gravitationally lensed radio-loud quasars, B\,1422+231, MG\,J0414+0534 and Q\,0957+561, using a combination of imaging and spectral analysis on the X-ray data taken from the \textit{Chandra X-Ray Observatory}. We tentatively detected FeK$\alpha$ emission lines in MG\,J0414+0534 and Q\,0957+561 with over 95\% significance, whereas, we did not significantly detect FeK$\alpha$ emission in B\,1422+231. We constructed differential microlensing light curves from absorption corrected count rates. We subsequently performed a microlensing analysis on the X-ray microlensing light curves to measure the X-ray source sizes in soft (0.83--3.6 keV), hard (3.6--21.8 keV), and full (0.83--21.8 keV) bands, based on either Bayesian or maximum likelihood probabilities. For B\,1422+231, sizes from the two methods are consistent with each other, e.g. $R_X^{hard}/R_G = 6.17\pm5.48 \text{ (Bayesian), } 11.81\pm3.75 \text{ (maximum likelihood)}$, where $R_G=GM_{BH}/c^2)$. However, for MG\,J0414+0534 and Q\,0957+561, the two methods yield completely different results suggesting that more frequently sampled data with better signal-to-noise ratio are needed to measure the source size for these two objects. Comparing the acquired size values with the radio-quiet sample in the literature we found that our results are consistent with X-ray source size scaling approximately as $R_X \propto M_{BH}$ with the mass of the central supermassive black hole. Our results also indicate that radio-loud quasars tend to have larger unresolved X-ray emission sizes compared to the radio-quiet ones.

We present an analysis of Chandra spectra of five gravitationally lensed active galactic nuclei. We confirm the previous detections of FeK$\alpha$ emission lines in most images of these objects with high significance. The line energies range from 5.8 to 6.8 keV with widths from unresolved to 0.6 keV, consistent with emission close to spinning black holes viewed at different inclination angles. We also confirm the positive offset from the Iwasawa-Taniguchi effect, the inverse correlation between the FeK$\alpha$ equivalent width and the X-ray luminosity in AGN, where our measured equivalent widths are larger in lensed quasars. We attribute this effect to microlensing, and perform a microlensing likelihood analysis to constrain the emission size of the relativistic reflection region and the spin of supermassive black holes, assuming that the X-ray corona and the reflection region, responsible for the iron emission line, both follow power-law emissivity profiles. The microlensing analysis yields strong constraints on the spin and emissivity index of the reflection component for Q2237+0305, with $a > 0.92$ and $n > 5.4$. For the remaining four targets, we jointly constrain the two parameters, yielding $a=0.8\pm0.16$ and an emissivity index of $n=4.0\pm 0.8$, suggesting that the relativistic X-ray reflection region is ultra-compact and very close to the innermost stable circular orbits of black holes, which are spinning at close to the maximal value. We successfully constrain the half light radius of the emission region to $< 2.4$ $r_g$ ($r_g = GM/c^2$) for Q2237+0305 and in the range 5.9--7.4 $r_g$ for the joint sample.

Microlensing by stars within distant galaxies acting as strong gravitational lenses of multiply-imaged quasars, provides a unique and direct measurement of the internal structure of the lensed quasar on nano-arcsecond scales. The measurement relies on the temporal variation of high-magnification caustic crossings which vary on timescales of days to years. Multiwavelength observations provide information from distinct emission regions in the quasar. Through monitoring of these strong gravitational lenses, a full tomographic view can emerge with Astronomical-Unit scale resolution. Work to date has demonstrated the potential of this technique in about a dozen systems. In the 2020s there will be orders of magnitude more systems to work with. Monitoring of lens systems for caustic-crossing events to enable triggering of multi-platform, multi-wavelength observations in the 2020s will fulfill the potential of quasar microlensing as a unique and comprehensive probe of active black hole structure and dynamics.

Optical and X-ray observations of strongly gravitationally lensed quasars (especially when four separate images of the quasar are produced) determine not only the amount of matter in the lensing galaxy but also how much is in a smooth component and how much is composed of compact masses (e.g., stars, stellar remnants, primordial black holes, CDM sub-halos, and planets). Future optical surveys will discover hundreds to thousands of quadruply lensed quasars, and sensitive X-ray observations will unambiguously determine the ratio of smooth to clumpy matter at specific locations in the lensing galaxies and calibrate the stellar mass fundamental plane, providing a determination of the stellar $M/L$. A modest observing program with a sensitive, sub-arcsecond X-ray imager, combined with the planned optical observations, can make those determinations for a large number (hundreds) of the lensing galaxies, which will span a redshift range of $\sim$$0.25<z<1.5$

We are starting a project to find $\gamma$-ray millisecond pulsars (MSPs) among the unidentified sources detected by the Large Area Telescope (LAT) onboard the {Fermi Gamma-Ray Space Telescope (Fermi)}, by radio observations. The selection of good candidates from analysis of the LAT data is an important part of the project. Given that there is more than 10 years worth of LAT data and the advent of the newly released LAT 8-year point source list (FL8Y), we have conducted a selection analysis, on the basis of our previous analysis, and report the results here. Setting the requirements for the unidentified sources in FL8Y of Galactic latitudes $|b|>$ 5$^\circ$ and curvature significances $>$3$\sigma$, there are 202 sources with detection signficances $>$6$\sigma$. We select 57 relatively bright ones (detection significances $>$15$\sigma$) and analyze their 10.2 years of LAT data. Their variability is checked to exclude variable sources (likely~blazars), test statistic maps are constructed to avoid contaminated sources, and curvature significances are re-obtained and compared to their $\gamma$-ray spectra to exclude non-significant sources. In the end, 48 candidates are found. Based on the available information, mostly from multi-wavelength studies, we discuss the possible nature of several of the candidates. Most of these candidates are currently being observed with the 65-meter Shanghai Tian Ma Radio Telescope.

Many analyses have concluded that the accretion disc sizes measured from the microlensing variability of quasars are larger than the expectations from the standard thin disc theory by a factor of $\sim4$. We propose a simply model by invoking a strong wind from the disc to flatten its radial temperature profile, which can then reconcile the size discrepancy problem. This wind model has been successfully applied to several microlensed quasars with a wind strength $s\lesssim1.3$ by only considering the inward decreasing of the mass accretion rate (where $s$ is defined through $\dot{M}(R)\propto({R}/{R_{0}})^{s}$ ). After further incorporating the angular momentum transferred by the wind, our model can resolve the disc size problem with an even lower wind parameter. The corrected disc sizes under the wind model are correlated with black hole masses with a slope in agreement with our modified thin disc model.

We analyze 168 Swift monitoring observations of the nearest broad absorption line quasar Mrk 231 in the UV and X-ray bands, where we detect significant variability in the UV ($\sim$2246\AA) light curve with a null probability of $4.3\times10^{-10}$ for a constant model. Separately, from an archival sample of Swift observed active galactic nuclei (AGN), we measure the relation between UV excess variance and luminosity, finding that the normalized UV excess variance decreases with luminosity. Comparing to this mean relation, the normalized UV excess variance of Mrk 231 is smaller, however within the scatter characterising the full population. The upper limit of the X-ray excess variance is consistent with other AGN. The power spectrum density of the UV light curve can be well fit by a power law model with a slope of $1.82\pm0.14$ between $10^{-7.5}$ and $10^{-6}$ Hz, consistent with those for typical AGN, with no obvious quasi-periodical oscillation peaks. The UV variability and its power spectrum suggest that a significant amount of the UV emission of Mrk 231 is from the accretion disk. The consistencies in the normalized UV variability and the shape of the power spectrum density between Mrk 231 and other normal AGN suggest that the origin of UV variability of broad absorption line quasars is similar to other AGN, and dust scattering at large scales such as the torus is not a dominating process for the UV emission of Mrk 231. Significant scattering, if present, is constrained to smaller than $\sim$10 light days. We perform lagged correlation analysis between the UV and X-ray light curves and find the correlation insignificant within the present data.

We define a second-order moment of the observational differential microlensing curves that can be used to impose constraints on physical properties of lensed quasars. We show that this quantity is sensitive both to variations in the source size and the deflector mass. We formulize a methodology to recover the source size from the observational measurements when the mass spectrum is fixed. As a case study, we test it with a sample of four quadruple lenses, both in simulated scenarios and with real data from the $\textit{Chandra X-ray Observatory}$. In our simulations with a uniform stellar population the method works best to detect sources around $0.1$ Einstein radii, giving correct upper/lower limits for much smaller/bigger sizes without requiring a big leap in additional computational effort as compared to a single-epoch approach, yet taking advantage of multi-epoch information. We apply the method to a small sample of X-ray data from four objects assuming a range of star masses, and obtain a degeneracy relation between the source size and deflector mass. Combined with previous estimates for the size of the X-ray corona, the degeneracy relation suggests that X-ray microlensing is mainly induced by planetary mass objects.

We summarize and reanalyze observations bearing upon missing galactic baryons, where we propose a consistent picture for halo gas in L >~ L* galaxies. The hot X-ray emitting halos are detected to 50-70 kpc, where typically, M_hot(<50 kpc) ~ 5E9 Msun, and with density n \propto r^-3/2. When extrapolated to R200, the gas mass is comparable to the stellar mass, but about half of the baryons are still missing from the hot phase. If extrapolated to 1.9-3 R200, the baryon to dark matter ratio approaches the cosmic value. Significantly flatter density profiles are unlikely for R < 50 kpc and they are disfavored but not ruled out for R > 50 kpc. For the Milky Way, the hot halo metallicity lies in the range 0.3-1 solar for R < 50 kpc. Planck measurements of the thermal Sunyaev-Zeldovich effect toward stacked luminous galaxies (primarily early-type) indicate that most of their baryons are hot, near the virial temperature, and extend beyond R200. This stacked SZ signal is nearly an order of magnitude larger than that inferred from the X-ray observations of individual (mostly spiral) galaxies with M_* > 10^11.3 Msun. This difference suggests that the hot halo properties are distinct for early and late type galaxies, possibly due to different evolutionary histories. For the cooler gas detected in UV absorption line studies, we argue that there are two absorption populations: extended halos; and disks extending to ~50 kpc, containing most of this gas, and with masses a few times lower than the stellar masses. Such extended disks are also seen in 21 cm HI observations and in simulations.

Previously, planets have been detected only in the Milky Way galaxy. Here, we show that quasar microlensing provides a means to probe extragalactic planets in the lens galaxy, by studying the microlensing properties of emission close to the event horizon of the supermassive black hole of the background quasar, using the current generation telescopes. We show that a population of unbound planets between stars with masses ranging from Moon to Jupiter masses is needed to explain the frequent Fek line energy shifts observed in the gravitationally lensed quasar RXJ1131-1231 at a lens redshift of $z=0.295$ or 3.8 billion light-years away. We constrain the planet mass fraction to be larger than 0.0001 of the halo mass, which is equivalent to 2,000 objects ranging from Moon to Jupiter mass per main sequence star.

We present Hubble Space Telescope imaging confirming the optical disappearance of the failed supernova (SN) candidate identified by Gerke et al. (2015). This $\sim 25~M_{\odot}$ red supergiant experienced a weak $\sim 10^{6}~L_{\odot}$ optical outburst in 2009 and is now at least 5 magnitudes fainter than the progenitor in the optical. The mid-IR flux has slowly decreased to the lowest levels since the first measurements in 2004. There is faint ($2000-3000~L_{\odot}$) near-IR emission likely associated with the source. We find the late-time evolution of the source to be inconsistent with obscuration from an ejected, dusty shell. Models of the spectral energy distribution indicate that the remaining bolometric luminosity is $>6$ times fainter than that of the progenitor and is decreasing as $\sim t^{-4/3}$. We conclude that the transient is unlikely to be a SN impostor or stellar merger. The event is consistent with the ejection of the envelope of a red supergiant in a failed SN and the late-time emission could be powered by fallback accretion onto a newly-formed black hole. Future IR and X-ray observations are needed to confirm this interpretation of the fate for the star.

We combine new Lucky Imaging astrometry from NTT/AstraLux Sur with already published astrometry from the AstraLux Large M-dwarf Multiplicity Survey to compute orbital elements and individual masses of the 2MASS J10364483+1521394 triple system belonging to the Ursa-Major moving group. The system consists of one primary low-mass M-dwarf orbited by two less massive companions, for which we determine a combined dynamical mass of $M_{\rm{B}+\rm{C}}= 0.48 \pm 0.14\ M_\odot$. We show from the companions relative motions that they are of equal mass (with a mass ratio of $1.00 \pm 0.03$), thus $0.24 \pm 0.07\ M_\odot$ individually, with a separation of $3.2 \pm 0.3\ $AU and conclude that these masses are significantly higher ($30\%$) than what is predicted by theoretical stellar evolutionary models. The biggest uncertainty remains the distance to the system, here adopted as $20.1 \pm 2.0$ pc based on trigonometric parallax, whose ambiguity has a major impact on the result. With the new observational data we are able to conclude that the orbital period of the BC pair is $8.41^{+0.04}_{-0.02}\ $years.

We present $Suzaku$ off-center observations of two poor galaxy groups, NGC 3402 and NGC 5129, with temperatures below 1 keV. Through spectral decomposition, we measure their surface brightnesses and temperatures out to 330 and 680 times the critical density of the universe for NGC 3402 and NGC 5129, respectively. These quantities are consistent with extrapolations from existing inner measurements of the two groups. With the refined X-ray luminosities, both groups prefer $L_X-T$ relations without a break in the group regime. Furthermore, we measure the electron number densities and hydrostatic masses at these radii. We find that the electron number density profiles require three $\beta$ model components, with nearly flat slopes in the 3$^{rd}$ $\beta$ component for both groups. However, we find the effective slope in the outskirts to be $\beta_{out}$ = 0.59 and 0.49 for NGC 3402 and NGC 5129, respectively. Adding the gas mass measured from the X-ray data and stellar mass from group galaxy members, we measure baryon fractions of $f_b$ = 0.113 $\pm$ 0.013 and 0.091 $\pm$ 0.006 for NGC 3402 and NGC 5129, respectively. Combining other poor groups with well measured X-ray emission to the outskirts, we find an average baryon fraction of $f_{b,ave}$ = 0.100 $\pm$ 0.004 for X-ray bright groups with temperatures between 0.8$-$1.3 keV, extending existing constraints to lower mass systems.

Following our work presented in Dai et al. (2016), we report our detailed data analysis for another 38 Fermi gamma-ray un-associated sources. These sources are selected from the Fermi Large Area Telescope (LAT) third source catalog on the basis of the properties of known gamma-ray millisecond pulsars (MSPs) and for the purpose of finding likely candidate MSPs. From our analysis of the LAT data, we identify that among the 38 sources, 28 of them are single point-like sources with clean background and their spectra show significant curvature. We also conduct analysis of archival X-ray data available for 24 of the 28 sources. In the fields of 10 sources, there are at least one X-ray object, and in those of the other 14 sources, no X-ray object is detected but probably due to the X-ray observations being short. We discuss the possible MSP nature for these sources. Six of them(J0514.6-4406, J1035.7-6720, J1624.2-4041, J1744.1-7619, J1946.4-5403, and J2039.6-5618) are most likely associated with pulsars because of multi-wavelength identifications including direct radio or gamma-ray detection of pulsations. To firmly establish the associations or verify the MSP nature for other sources, deep X-ray and/or optical observations are needed.

Recently Graphics Processing Units (GPUs) have been used to speed up very CPU-intensive gravitational microlensing simulations. In this work, we use the Xeon Phi coprocessor to accelerate such simulations and compare its performance on a microlensing code with that of NVIDIA's GPUs. For the selected set of parameters evaluated in our experiment, we find that the speedup by Intel's Knights Corner coprocessor is comparable to that by NVIDIA's Fermi family of GPUs with compute capability 2.0, but less significant than GPUs with higher compute capabilities such as the Kepler. However, the very recently released second generation Xeon Phi, Knights Landing, is about 5.8 times faster than the Knights Corner, and about 2.9 times faster than the Kepler GPU used in our simulations. We conclude that the Xeon Phi is a very promising alternative to GPUs for modern high performance microlensing simulations.

The PICASSO dark matter search experiment operated an array of 32 superheated droplet detectors containing 3.0 kg of C$_{4}$F$_{10}$ and collected an exposure of 231.4 kgd at SNOLAB between March 2012 and January 2014. We report on the final results of this experiment which includes for the first time the complete data set and improved analysis techniques including \mbox{acoustic} localization to allow fiducialization and removal of higher activity regions within the detectors. No signal consistent with dark matter was observed. We set limits for spin-dependent interactions on protons of $\sigma_p^{SD}$~=~1.32~$\times$~10$^{-2}$~pb (90\%~C.L.) at a WIMP mass of 20 GeV/c$^{2}$. In the spin-independent sector we exclude cross sections larger than $\sigma_p^{SI}$~=~4.86~$\times$~10$^{-5 }$~pb~(90\% C.L.) in the region around 7 GeV/c$^{2}$. The pioneering efforts of the PICASSO experiment have paved the way forward for a next generation detector incorporating much of this technology and experience into larger mass bubble chambers.

We present a promising new technique, the g-distribution method, for measuring the inclination angle (i), the innermost stable circular orbit (ISCO), and the spin of a supermassive black hole. The g-distribution method uses measurements of the energy shifts in the relativistic iron line emitted by the accretion disk of a supermassive black hole due to microlensing by stars in a foreground galaxy relative to the g-distribution shifts predicted from microlensing caustic calculations. We apply the method to the gravitationally lensed quasars RX J1131-1231 (z_s=0.658, z_l=0.295), QJ 0158-4325 (z_s=1.294, z_l=0.317), and SDSS 1004+4112 (z_s=1.734, z_l=0.68). For RX J1131-1231 our initial results indicate that r_ISCO<8.5 gravitational radii (r_g) and i > 76 degrees. We detect two shifted Fe lines, in several observations, as predicted in our numerical simulations of caustic crossings. The current DeltaE-distribution of RX J1131-1231 is sparsely sampled but further X-ray monitoring of RX J1131-1231 and other lensed quasars will provide improved constraints on the inclination angles, ISCO radii and spins of the black holes of distant quasars.

We present an X-ray photometric analysis of six gravitationally lensed quasars, with observation campaigns spanning from 5 to 14 years, measuring the total (0.83 - 21.8 keV restframe), soft (0.83 - 3.6 keV), and hard (3.6 - 21.8 keV) band image flux ratios or each epoch. Using the ratios of the model-predicted macro-magnifications as baselines, we build differential microlensing light curves and obtain joint likelihood functions for the average X-ray emission region sizes. Our analysis yields a Probability Distribution Function for the average half-light radius of the X-Ray emission region in the sample that peaks slightly above 1 gravitational radius and with nearly indistinguishable 68% confidence (one-sided) upper limits of 17.8 and 18.9 gravitational radii for the soft and hard X-ray emitting regions, assuming a mean stellar mass of 0.3 solar masses. We see hints of energy dependent microlensing between the soft and hard bands in two of the objects. In a separate analysis on the root-mean-square of the microlensing variability, we find significant differences between the soft and hard bands but the sign of the difference is not consistent across the sample. This suggests the existence of some kind of spatial structure to the X-ray emission in an otherwise extremely compact source. We also discover a correlation between the RMS microlensing variability and the average microlensing amplitude.

The solar galactic neighbourhood contains a number of young co-moving associations of stars (so-called `young moving groups') with ages of ~10--150 Myr, which are prime targets for a range of scientific studies, including direct imaging planet searches. The late-type stellar population of such groups still remain in their pre-main sequence phase, and are thus well suited for purposes such as isochronal dating. Close binaries are particularly useful in this regard, since they allow for a model-independent dynamical mass determination. Here we present a dedicated effort to identify new close binaries in nearby young moving groups, through high-resolution imaging with the AstraLux Sur Lucky Imaging camera. We surveyed 181 targets, resulting in the detection of 61 companions or candidates, of which 38 are new discoveries. An interesting example of such a case is 2MASS J00302572-6236015 AB, which is a high-probability member of the Tucana-Horologium moving group, and has an estimated orbital period of less than 10 years. Among the previously known objects is a serendipitous detection of the deuterium burning boundary circumbinary companion 2MASS J01033563-5515561 (AB)b in the z'-band, thereby extending the spectral coverage for this object down to near-visible wavelengths.

The 2015 March 15 coronal mass ejection as one of the two that together drove the largest geomagnetic storm of solar cycle 24 so far was associated with sympathetic filament eruptions. We investigate the relations between the different filaments involved in the eruption. A surge-like small-scale filament motion is confirmed as the trigger that initiated the erupting filament with multi-wavelength observations and using a forced magnetic field extrapolation method. When the erupting filament moved to an open magnetic field region, it experienced an obvious acceleration process and was accompanied by a C-class flare and the rise of another larger filament that eventually failed to erupt. We measure the decay index of the background magnetic field, which presents a critical height of 118 Mm. Combining with a potential field source surface extrapolation method, we analyze the distributions of the large-scale magnetic field, which indicates that the open magnetic field region may provide a favorable condition for F2 rapid acceleration and have some relation with the largest solar storm. The comparison between the successful and failed filament eruptions suggests that the confining magnetic field plays an important role in the preconditions for an eruption.

Star-forming galaxies produce gamma-rays primarily via pion production, resulting from inelastic collisions between cosmic ray protons and the interstellar medium (ISM). The dense ISM and high star formation rates of luminous and ultra-luminous infrared galaxies (LIRGs and ULIRGs) imply that they should be strong gamma-ray emitters, but so far only two LIRGs have been detected. Theoretical models for their emission depend on the unknown fraction of cosmic ray protons that escape these galaxies before interacting. We analyze Fermi-LAT data for 82 of the brightest IRAS LIRGs and ULIRGs. We examine each system individually and carry out a stacking analysis to constrain their gamma-ray fluxes. We report the detection of the nearest ULIRG Arp 220 (~4.6sigma). We observe a gamma-ray flux (0.8--100 GeV) of 2.4e-10 phot cm^-2 s^-1 with photon index of 2.23 (8.2e10^41 ergs s^-1 at 77 Mpc) We also derive upper limits for the stacked LIRGs and ULIRGs. The gamma-ray luminosity of Arp~220 and the stacked upper limits agree with calorimetric predictions for dense star-forming galaxies. With the detection of Arp 220, we extend the gamma-ray--IR luminosity correlation to the high luminosity regime with log(L_[0.1-100 GeV]) = 1.25log(L_[8-1000]) + 26.7 as well as the gamma-ray--radio continuum luminosity correlation with log(L_[0.1-100 GeV]) = 1.22log(L_[1.4 GHz]) + 13.3. The current survey of Fermi-LAT is on the verge of detecting more LIRGs/ULIRGs in the local universe, and we expect even more detections with deeper Fermi-LAT observations or the next generation of gamma-ray detectors.

We report our detailed data analysis for 39 $\gamma$-ray sources selected from the 992 unassociated sources in the \textit{Fermi} Large Area Telescope (LAT) third source catalog. The selection criteria, which were set for finding candidate millisecond pulsars (MSPs), are non-variables with curved spectra and $>$5$^{\circ}$ Galactic latitudes. From our analysis, 24 sources were found to be point-like sources not contaminated by background or nearby unknown sources. Three of them, J1544.6$-$1125, J1625.1$-$0021, and J1653.6$-$0158, have been previously studied, indicating that they are likely MSPs. The spectra of J0318.1+0252 and J2053.9+2922 do not have properties similar to that of known $\gamma$-ray MSPs, and we thus suggest that they are not MSPs. Analysis of archival X-ray data for most of the 24 sources were also conducted. Four sources were found with X-ray objects in their error circles, and 16 with no detection. The ratios between the $\gamma$-ray fluxes and X-ray fluxes or flux upper limits are generally lower than those of the known $\gamma$-ray MSPs, suggesting that if the $\gamma$-ray sources are MSPs, none of the X-ray objects are the counterparts. Deep X-ray or radio observations of these sources are needed in order to identify their MSP nature.

We study 203 (of 442) Swift AGN and Cluster Survey extended X-ray sources located in the SDSS DR8 footprint to search for galaxy over-densities in three dimensional space using SDSS galaxy photometric redshifts and positions near the Swift cluster candidates. We find 104 Swift clusters with a >3sigma galaxy over-density. The remaining targets are potentially located at higher redshifts and require deeper optical follow-up observations for confirmation as galaxy clusters. We present a series of cluster properties including the redshift, BCG magnitude, BCG-to-X-ray center offset, optical richness, and X-ray luminosity. We also detect red sequences in ~85% of the 104 confirmed clusters. The X-ray luminosity and optical richness for the SDSS confirmed Swift clusters are correlated and follow previously established relations. The distribution of the separations between the X-ray centroids and the most likely BCG is also consistent with expectation. We compare the observed redshift distribution of the sample with a theoretical model, and find that our sample is complete for z <~ 0.3 and is still 80% complete up to z ~= 0.4, consistent with the SDSS survey depth. These analysis results suggest that our Swift cluster selection algorithm has yielded a statistically well-defined cluster sample for further studying cluster evolution and cosmology. We also match our SDSS confirmed Swift clusters to existing cluster catalogs, and find 42, 23 and 1 matches in optical, X-ray and SZ catalogs, respectively, so the majority of these clusters are new detections.

We analyzed a large sample of radio-loud and radio-quiet quasar spectra at redshift 1.0 < z < 1.2 to compare the inferred underlying quasar continuum slopes (after removal of the host galaxy contribution) with accretion disk models. The latter predict redder (decreasing) alpha_3000 continuum slopes (L_\nu~\nu^alpha at 3000Ang) with increasing black hole mass, bluer alpha_3000 with increasing luminosity at 3000Ang, and bluer alpha_3000 with increasing spin of the black hole, when all other parameters are held fixed. We find no clear evidence for any of these predictions in the data. In particular we find that: (i) alpha_3000 shows no significant dependence on black hole mass or luminosity. Dedicated Monte Carlo tests suggest that the substantial observational uncertainties in the black hole virial masses can effectively erase any intrinsic dependence of alpha_3000 on black hole mass, in line with some previous studies. (ii) The mean slope alpha_3000 of radio-loud sources, thought to be produced by rapidly spinning black holes, is comparable to, or even redder than, that of radio-quiet quasars. Indeed, although quasars appear to become more radio loud with decreasing luminosity, we still do not detect any significant dependence of alpha_3000 on radio loudness. The predicted mean alpha_3000 slopes tend to be bluer than in the data. Disk models with high inclinations and dust extinction tend to produce redder slopes closer to empirical estimates. Our mean alpha_3000 values are close to the ones independently inferred at z<0.5 suggesting weak evolution with redshift, at least for moderately luminous quasars.

We present late-time Hubble and Spitzer Space Telescope imaging of SN 2008S and NGC 300 2008OT-1, the prototypes of a common class of stellar transients whose true nature is debated. Both objects are still fading and are now >15 times fainter than the progenitors in the mid-IR and are undetected in the optical and near-IR. Data from the Large Binocular Telescope and Magellan show that neither source has been variable in the optical since fading in 2010. We present models of surviving sources obscured by dusty shells or winds and find that extreme dust models are needed for surviving stars to be successfully hidden by dust. Explaining these transients as supernovae explosions, such as the electron capture supernovae believed to be associated with extreme AGB stars, seems an equally viable solution. Though SN 2008S is not detected in Chandra X-Ray Observatory data taken in 2012, the flux limits allow the fading IR source to be powered solely by the shock interaction of ejecta with the circumstellar medium if the shock velocity at the time of the observation was >20% slower than estimated from emission line widths while the transient was still optically bright. Continued SST monitoring and 10-20 micron observations with JWST can resolve any remaining ambiguities.

We study the role of the coronal magnetic field configuration of an active region in determining the propagation direction of a coronal mass ejection (CME). The CME occurred in the active region 11944 (S09W01) near the disk center on 2014 January 7 and was associated with an X1.2 flare. A new CME reconstruction procedure based on a polarimetric technique is adopted, which shows that the CME changed its propagation direction by around 28$^\circ$ in latitude within 2.5 R$_\odot$ and 43$^\circ$ in longitude within 6.5 R$_\odot$ with respect to the CME source region. This significant non-radial motion is consistent with the finding of M$\ddot{o}$stl et al. (2015). We use nonlinear force-free field (NLFFF) and potential field source surface (PFSS) extrapolation methods to determine the configurations of the coronal magnetic field. We also calculate the magnetic energy density distributions at different heights based on the extrapolations. Our results show that the active region coronal magnetic field has a strong influence on the CME propagation direction. This is consistent with the "channelling" by the active region coronal magnetic field itself, rather than deflection by nearby structures. These results indicate that the active region coronal magnetic field configuration has to be taken into account in order to determine CME propagation direction correctly.

We review results from our monitoring observations of several lensed quasars performed in the optical, UV, and X-ray bands. Modeling of the multi-wavelength light curves provides constraints on the extent of the optical, UV, and X-ray emission regions. One of the important results of our analysis is that the optical sizes as inferred from the microlensing analysis are significantly larger than those predicted by the theoretical-thin-disk estimate. In a few cases we also constrain the slope of the size-wavelength relation. Our size constraints of the soft and hard X-ray emission regions of quasars indicate that in some objects of our sample the hard X-ray emission region is more compact than the soft and in others the soft emission region is smaller. This difference may be the result of the relative strengths of the disk-reflected (harder and extended) versus corona-direct (softer and compact) components in the quasars of our sample. Finally, we present the analysis of several strong microlensing events where we detect an evolution of the relativistic Fe line profile as the magnification caustic traverses the accretion disk. These caustic crossings are used to provide constraints on the innermost stable circular orbit (ISCO) radius and the accretion disk inclination angle of the black hole in quasar RX J1131-1231.

Supermassive binary black holes (BBHs) are unavoidable products of galaxy mergers and are expected to exist in the cores of many quasars. Great effort has been made during the past several decades to search for BBHs among quasars; however, observational evidence for BBHs remains elusive and ambiguous, which is difficult to reconcile with theoretical expectations. In this paper, we show that the distinct optical-to-UV spectrum of Mrk 231 can be well interpreted as emission from accretion flows onto a BBH, with a semimajor axis of ~590AU and an orbital period of ~1.2 year. The flat optical and UV continua are mainly emitted from a circumbinary disk and a mini-disk around the secondary black hole (BH), respectively; and the observed sharp drop off and flux deficit at wavelength lambda ~ 4000-2500 Angstrom is due to a gap (or hole) opened by the secondary BH migrating within the circumbinary disk. If confirmed by future observations, this BBH will provide a unique laboratory to study the interplay between BBHs and accretion flows onto them. Our result also demonstrates a new method to find sub-parsec scale BBHs by searching for deficits in the optical-to-UV continuum among the spectra of quasars.

We analyze the optical, UV, and X-ray microlensing variability of the lensed quasar SDSS J0924+0219 using six epochs of Chandra data in two energy bands (spanning 0.4-8.0 keV, or 1-20 keV in the quasar rest frame), 10 epochs of F275W (rest-frame 1089A) Hubble Space Telescope data, and high-cadence R-band (rest-frame 2770A) monitoring spanning eleven years. Our joint analysis provides robust constraints on the extent of the X-ray continuum emission region and the projected area of the accretion disk. The best-fit half-light radius of the soft X-ray continuum emission region is between 5x10^13 and 10^15 cm, and we find an upper limit of 10^15 cm for the hard X-rays. The best-fit soft-band size is about 13 times smaller than the optical size, and roughly 7 GM_BH/c^2 for a 2.8x10^8 M_sol black hole, similar to the results for other systems. We find that the UV emitting region falls in between the optical and X-ray emitting regions at 10^14 cm < r_1/2,UV < 3x10^15 cm. Finally, the optical size is significantly larger, by 1.5*sigma, than the theoretical thin-disk estimate based on the observed, magnification-corrected I-band flux, suggesting a shallower temperature profile than expected for a standard disk.

Active galactic nuclei (AGNs) and quasars are important astrophysical objects to understand. Recently, microlensing observations have constrained the size of the quasar X-ray emission region to be of the order of 10 gravitational radii of the central supermassive black hole. For distances within a few gravitational radii, light paths are strongly bent by the strong gravity field of the central black hole. If the central black hole has nonzero angular momentum (spin), a photon's polarization plane will be rotated by the gravitational Faraday effect. The observed X-ray flux and polarization will then be influenced significantly by the strong gravity field near the source. Consequently, linear gravitational lensing theory is inadequate for such extreme circumstances. We present simple algorithms computing strong lensing effects of Kerr black holes, including effects on polarization. Our algorithms are realized in a program "KERTAP" in two versions: MATLAB and Python. The key ingredients of KERTAP are: a graphic user interface, a {\it backward} ray-tracing algorithm, a polarization propagator dealing with gravitational Faraday rotation, and algorithms computing observables such as flux magnification and polarization angles. Our algorithms can be easily realized in other programming languages such as FORTRAN, C, and C++. The MATLAB version of KERTAP is parallelized using the MATLAB Parallel Computing Toolbox and the Distributed Computing Server. The Python code was sped up using Cython and supports full implementation of MPI using 'mpi4py' package. As an example, we investigate the inclination angle dependence of the observed polarization and the strong lensing magnification of AGN X-ray emission. We conclude that it is possible to perform complex numerical-relativity-related computations using interpreted languages such as MATLAB and Python.

We present a simple gravitational lens model to illustrate the ease of using the embedded lensing theory when studying cosmic voids. It confirms the previously used repulsive lensing models for deep voids. We start by estimating magnitude fluctuations and weak lensing shears of background sources lensed by large voids. We find that sources behind large ($\sim$$90\,\rm Mpc$) and deep voids (density contrast about $-0.9$) can be magnified or demagnified with magnitude fluctuations of up to $\sim$$0.05\,\rm mag$ and that the weak-lensing shear can be up to the $\sim$$10^{-2}$ level in the outer regions of large voids. Smaller or shallower voids produce proportionally smaller effects. We investigate the "wiggling" of the primary cosmic microwave background (CMB) temperature anisotropies caused by intervening cosmic voids. The void-wiggling of primary CMB temperature gradients is of the opposite sign to that caused by galaxy clusters. Only extremely large and deep voids can produce wiggling amplitudes similar to galaxy clusters, $\sim$$\rm 15\,\mu K$ by a large void of radius $\sim$$4^\circ$ and central density contrast $-0.9$ at redshift 0.5 assuming a CMB background gradient of $\sim$$\rm10\,\mu K\, arcmin^{-1}$. The dipole signal is spread over the entire void area, and not concentrated at the lens' center as it is for clusters. Finally we use our model to simulate CMB sky maps lensed by large cosmic voids. Our embedded theory can easily be applied to more complicated void models and used to study gravitational lensing of the CMB, to probe dark-matter profiles, to reduce the lensing-induced systematics in supernova Hubble diagrams, as well as study the integrated Sachs-Wolfe effect.

We describe a method using the integrated Sachs-Wolfe (ISW) effect caused by individual inhomogeneities to determine the cosmological parameters, $H_0$, $\Omega_{\rm m}$, and $\Omega_\Lambda$, etc. This ISW-redshift test requires detailed knowledge of the internal kinematics of a set of individual density perturbations, e.g., galaxy clusters and/or cosmic voids, in particular their density and velocity profiles, and their mass accretion rates. It assumes the density perturbations are isolated and imbedded (equivalently compensated) and makes use of the newly found relation between the ISW temperature perturbation of the CMB and the Fermat potential of the lens. Given measurements of the amplitudes of the temperature variations in the CMB caused by such clusters or voids at various redshifts and estimates of their angular sizes or masses, one can constrain the cosmological parameters. More realistically, the converse is more likely, i.e., if the background cosmology is sufficiently constrained, measurement of ISW profiles of clusters and voids (e.g., hot and cold spots and rings) can constrain dynamical properties of the dark matter, including accretion, associated with such lenses and thus constrain the evolution of these objects with redshift.

We derive the Fermat potential for a spherically symmetric lens embedded in an FLRW cosmology and use it to investigate the late-time integrated Sachs-Wolfe (ISW) effect, i.e., secondary temperature fluctuations in the cosmic microwave background (CMB) caused by individual large scale clusters and voids. We present a simple analytical expression for the temperature fluctuation in the CMB across such a lens as a derivative of the lens' Fermat potential. This formalism is applicable to both linear and nonlinear density evolution scenarios, to arbitrarily large density contrasts, and to all open and closed background cosmologies. It is much simpler to use and makes the same predictions as conventional approaches. In this approach the total temperature fluctuation can be split into a time-delay part and an evolutionary part. Both parts must be included for cosmic structures that evolve and both can be equally important. We present very simple ISW models for cosmic voids and galaxy clusters to illustrate the ease of use of our formalism. We use the Fermat potentials of simple cosmic void models to compare predicted ISW effects with those recently extracted from WMAP and \emph{Planck} data by stacking large cosmic voids using the aperture photometry method. If voids in the local universe with large density contrasts are no longer evolving we find that the time delay contribution alone predicts values consistent with the measurements. However, we find that for voids still evolving linearly, the evolutionary contribution cancels a significant part of the time delay contribution and results in predicted signals that are much smaller than recently observed.

The Swift AGN and Cluster Survey (SACS) uses 125 deg^2 of Swift XRT serendipitous fields with variable depths surrounding gamma-ray bursts to provide a medium depth (4e-15 erg/s/cm^2) and area survey filling the gap between deep, narrow Chandra/XMM-Newton surveys and wide, shallow ROSAT surveys. Here we present a catalog of 22,563 point sources and 442 extended sources and examine the number counts of the AGN and galaxy cluster populations. SACS provides excellent constraints on the AGN number counts at the bright end with negligible uncertainties due to cosmic variance, and these constraints are consistent with previous measurements. We use Wise mid-infrared (MIR) colors to classify the sources. For AGN we can roughly separate the point sources into MIR-red and MIR-blue AGN, finding roughly equal numbers of each type in the soft X-ray band (0.5-2 keV), but fewer MIR-blue sources in the hard X-ray band (2-8 keV). The cluster number counts, with 5% uncertainties from cosmic variance, are also consistent with previous surveys but span a much larger continuous flux range. Deep optical or IR follow-up observations of this cluster sample will significantly increase the number of higher redshift (z > 0.5) X-ray-selected clusters.

We examine a sample of $\sim 250 000$ "locally brightest galaxies" selected from the Sloan Digital Sky Survey to be central galaxies within their dark matter halos. These were previously stacked by the Planck Collaboration to measure the Sunyaev-Zel'dovich signal as a function of central galaxy stellar mass. Here, we stack the X-ray emission from these halos using data from the ROSAT All-Sky Survey. We detect emission across almost our entire sample, including emission which we attribute to hot gas around galaxies spanning a range of 1.2 dex in stellar mass (corresponding to nearly two orders of magnitude in halo mass) down to $M* = 10^{10.8} M_{\odot}$ ($M_{500} \approx 10^{12.6} M_{\odot}$). Over this range, the X-ray luminosity can be fit by a power-law, either of stellar mass or of halo mass. A single unified scaling relation between mass and $L_X$ applies for galaxies, groups, and clusters. This relation has a steeper slope than expected for self-similarity, in contrast to the $Y_{SZ}$-$M_{500}$ relation, showing the importance of non-gravitational heating. If this heating is predominantly due to AGN feedback, the lack of a break in our relation suggests that AGN feedback is tightly self-regulated and fairly gentle, in agreement with recent high-resolution simulations. Our results are consistent with earlier measurements of the $L_X$-$L_K$ relation for elliptical galaxies and of the $L_X$-$M_{500}$ relation for optically-selected galaxy clusters. However, our $L_X$-$M_{500}$ relation lies more than a factor of two below most previous relations based on X-ray-selected cluster samples. We argue that optical selection offers a less biased view of the $L_X$-$M_{500}$ relation.

Galaxy clusters are predicted to produce gamma-rays through cosmic ray interactions and/or dark matter annihilation, potentially detectable by the Fermi Large Area Telescope (Fermi-LAT). We present a new, independent stacking analysis of Fermi-LAT photon count maps using the 78 richest nearby clusters (z<0.12) from the Two Micron All-Sky Survey (2MASS) cluster catalog. We obtain the lowest limit on the photon flux to date, 2.3e-11 ph/s/cm^2 (95% confidence) per cluster in the 0.8-100 GeV band, which corresponds to a luminosity limit of 3.5e44 ph/s. We also constrain the emission limits in a range of narrower energy bands. Scaling to recent cosmic ray acceleration and gamma-ray emission models, we find that cosmic rays represent a negligible contribution to the intra-cluster energy density and gas pressure.

Microlensing has proven an effective probe of the structure of the innermost regions of quasars, and an important test of accretion disk models. We present light curves of the lensed quasar HE 0435-1223 in the R band and in the ultraviolet, and consider them together with X-ray light curves in two energy bands that are presented in a companion paper. Using a Bayesian Monte Carlo method, we constrain the size of the accretion disk in the rest-frame near- and far-UV, and constrain for the first time the size of the X-ray emission regions in two X-ray energy bands. The R-band scale size of the accretion disk is about 10^15.23 cm (~23 r_g), slightly smaller than previous estimates, but larger than would be predicted from the quasar flux. In the UV, the source size is weakly constrained, with a strong prior dependence. The UV to R-band size ratio is consistent with the thin disk model prediction, with large error bars. In soft and hard X-rays, the source size is smaller than ~10^14.8 cm (~10 r_g) at 95% confidence. We do not find evidence of structure in the X-ray emission region, as the most likely value for the ratio of the hard X-ray size to the soft X-ray size is unity. Finally, we find that the most likely value for the mean mass of stars in the lens galaxy is ~0.3 M_sun, consistent with other studies.

There is growing evidence of a higher intrinsic fraction of broad absorption line quasars (BALQSOs) than that obtained in optical surveys, on which most previous X-ray studies of BALQSOs have focused. Here we present Chandra observations of 18 BALQSOs at $z\sim2$, selected from a near-infrared (2MASS) sample, where the BALQSO fraction is likely to be close to the intrinsic fraction. We measure photon indices using the stacked spectra of the optically-faint ($i-K_s\geq 2.3$ mag) and optically-bright ($i-K_s < 2.3$ mag) samples to be $\Gamma \simeq 1.5$--$2.1$. We constrain their intrinsic column density by modelling the X-ray fractional hardness ratio, finding a mean column density of $3.5\times10^{22}$ \cmsq\ assuming neutral absorption. We incorporate SDSS optical measurements (rest frame UV) to study the broadband spectral index between the X-ray and UV bands, and compare this to a large sample of normal quasars. We estimate that the optically-faint BALQSOs are X-ray weaker than the optically-bright ones, and the entire sample of BALQSOs are intrinsically X-ray weak when compared to normal AGN. Correcting for magnification of X-ray emission via gravitational lensing by the central black hole viewed at large inclination angles makes these BALQSOs even more intrinsically X-ray weak. Finally, we estimate AGN kinetic feedback efficiencies of a few percent for an X-ray wind of $0.3c$ in high-ionization BALQSOs. Combined with energy carried by low-ionization BALQSOs and UV winds, the total kinetic energy in BALQSOs can be sufficient to provide AGN kinetic feedback required to explain the co-evolution between black holes and host galaxies.

We present a simplified version of the lowest-order embedded point mass gravitational lens theory and then make the extension of this theory to any embedded transparent lens. Embedding a lens effectively reduces the gravitational potential's range, i.e., partially shields the lensing potential because the lens mass is made a contributor to the mean mass density of the universe and not simply superimposed upon it. We give the time-delay function for the embedded point mass lens from which we can derive the simplified lens equation by applying Fermat's least-time principle. Even though rigorous derivations are only made for the point mass in a flat background, the generalization of the lens equation to lowest-order for any distributed lens in any homogeneous background is obvious. We find from this simplified theory that embedding can introduce corrections above the few percent level in weak lensing shears caused by large clusters but only at large impacts. The potential part of the time delay is also affected in strong lensing at the few percent level. Additionally we again confirm that the presence of a cosmological constant alters the gravitational deflection of passing photons.

The Sudbury Neutrino Observatory (SNO) has confirmed the standard solar model and neutrino oscillations through the observation of neutrinos from the solar core. In this paper we present a search for neutrinos associated with sources other than the solar core, such as gamma-ray bursters and solar flares. We present a new method for looking for temporal coincidences between neutrino events and astrophysical bursts of widely varying intensity. No correlations were found between neutrinos detected in SNO and such astrophysical sources.

We homogeneously analyzed the \chandra\ X-ray observations of 10 gravitational lenses, HE 0047-1756, QJ 0158-4325, SDSS 0246-0805, HE 0435-1223, SDSS 0924+0219, SDSS 1004+4112, HE 1104-1805, PG 1115+080, Q 1355-2257, and Q 2237+0305, to measure the differential X-ray absorption between images, the metallicity, and the dust-to-gas ratio of the lens galaxies. We detected differential absorption in all lenses except SDSS 0924+0219 and HE 1104-1805. This doubles the sample of dust-to-gas ratio measurements in cosmologically distant lens galaxies. We successfully measured the gas phase metallicity of three lenses, Q 2237+0305, SDSS 1004+4112, and B 1152+199 from the X-ray spectra. Our results suggest a linear correlation between metallicity and dust-to-gas ratio (i.e., a constant metal-to-dust ratio), consistent with what is found for nearby galaxies. We obtain an average dust-to-gas ratio $E(B-V)/N_H=1.17^{+0.41}_{-0.31} \times 10^{-22}\rm mag\,cm^2\,atom^{-1}$ in the lens galaxies, with an intrinsic scatter of $\rm0.3\,dex$. Combining these results with data from GRB afterglows and quasar foreground absorbers, we found a mean dust-to-gas ratio $\mdtg,$ now significantly lower than the average Galactic value, $1.7\,\times 10^{-22}\,\rm mag\, cm^{2}\, atoms^{-1}.$ This suggests evolution of dust-to-gas ratios with redshift and lower average metallicities for the higher redshift galaxies, consistent with current metal and dust evolution models of interstellar medium. The slow evolution in the metal-to-dust ratio with redshift implies very rapid dust formation in high redshift ($z>2$) galaxies.

The gravitationally lensed quasar HE 1104-1805 has been observed at a variety of wavelengths ranging from the mid-infrared to X-ray for nearly 20 years. We combine flux ratios from the literature, including recent Chandra data, with new observations from the SMARTS telescope and HST, and use them to investigate the spatial structure of the central regions using a Bayesian Monte Carlo analysis of the microlensing variability. The wide wavelength coverage allows us to constrain not only the accretion disk half-light radius r_1/2, but the power-law slope \xi\ of the size-wavelength relation r_1/2 ~ \lambda^\xi. With a logarithmic prior on the source size, the (observed-frame) R-band half-light radius log(r_1/2/cm) is 16.0+0.3-0.4, and the slope \xi\ is 1.0+0.30-0.56. We put upper limits on the source size in soft (0.4-1.2 keV) and hard (1.2-8 keV) X-ray bands, finding 95% upper limits on log (r_1/2/cm) of 15.33 in both bands. A linear prior yields somewhat larger sizes, particularly in the X-ray bands. For comparison, the gravitational radius, using a black hole mass estimated using the H\beta\ line, is log(r_g/cm) = 13.94. We find that the accretion disk is probably close to face-on, with cos i = 1.0 being four times more likely than cos i = 0.5. We also find probability distributions for the mean mass of the stars in the foreground lensing galaxy, the direction of the transverse peculiar velocity of the lens, and the position angle of the projected accretion disk's major axis (if not face-on).

We present 279 epochs of optical monitoring data spanning 5.4 years from 2007 January to 2012 June for the largest image separation (22.6 arcsec) gravitationally lensed quasar, SDSS J1029+2623. We find that image A leads the images B and C by dt_AB = (744+-10) days (90% confidence); the uncertainty includes both statistical uncertainties and systematic differences due to the choice of models. With only a ~1% fractional error, the interpretation of the delay is limited primarily by cosmic variance due to fluctuations in the mean line-of-sight density. We cannot separate the fainter image C from image B, but since image C trails image B by only 2-3 days in all models, the estimate of the time delay between image A and B is little affected by combining the fluxes of images B and C. There is weak evidence for a low level of microlensing, perhaps created by the small galaxy responsible for the flux ratio anomaly in this system. Interpreting the delay depends on better constraining the shape of the gravitational potential using the lensed host galaxy, other lensed arcs and the structure of the X-ray emission.

We use gravitational microlensing to determine the size of the X-ray and optical emission regions of the quadruple lens system Q 2237+0305. The optical half-light radius, log(R_{1/2,V}/cm)=16.41\pm0.18 (at lambda_{rest}=2018 \AA), is significantly larger than the observed soft, log(R_{1/2,soft}/cm)=15.76^{+0.41}_{-0.34} (1.1-3.5 keV in the rest frame), and hard, log(R_{1/2,hard}/cm)=15.46^{+0.34}_{-0.29} (3.5-21.5 keV in the rest frame), band X-ray emission. There is a weak evidence that the hard component is more compact than the soft, with log(R_{1/2,soft}/R_{1/2,hard}) \sim 0.30^{+0.53}_{-0.45}. This wavelength-dependent structure agrees with recent results found in other lens systems using microlensing techniques, and favors geometries in which the corona is concentrated near the inner edge of the accretion disk. While the available measurements are limited, the size of the X-ray emission region appears to be roughly proportional to the mass of the central black hole.

We place general constraints on the luminosity and mass of hot X-ray emitting gas residing in extended "hot halos" around nearby massive galaxies. We examine stacked images of 2165 galaxies from the 2MASS Very Isolated Galaxy Catalog (2MVIG), as well as subsets of this sample based on galaxy morphology and K-band luminosity. We detect X-ray emission at high confidence (ranging up to nearly 10\sigma) for each subsample of galaxies. The average L_X within 50 kpc is 1.0\pm0.1 (statistical) \pm0.2 (systematic) x10^40 erg/s, although the early-type galaxies are more than twice as luminous as the late-type galaxies. Using a spatial analysis, we also find evidence for extended emission around five out of seven subsamples (the full sample, the luminous galaxies, early-type galaxies, luminous late-type galaxies, and luminous early-type galaxies) at 92.7%, 99.3%, 89.3%, 98.7%, and 92.1% confidence, respectively. Several additional lines of evidence also support this conclusion and suggest that about 1/2 of the total emission is extended, and about 1/3 of the extended emission comes from hot gas. For the sample of luminous galaxies, which has the strongest evidence for extended emission, the average hot gas mass is 4x10^9 Msun within 50 kpc and the implied accretion rate is 0.4 Msun/yr.

Recent quasar microlensing observations have constrained the sizes of X-ray emission regions to be within about 10 gravitational radii of the central supermassive black hole. Therefore, the X-ray emission from lensed quasars is first strongly lensed by the black hole before it is lensed by the foreground galaxy and star fields. We present a scheme that combines the initial strong lensing of a Kerr black hole with standard linearized microlensing by intervening stars. We find that X-ray microlensed light curves incorporating Kerr strong gravity can differ significantly from standard curves. The amplitude of the fluctuations in the light curves can increase or decrease by ~0.65-0.75 mag by including Kerr strong gravity. Larger inclination angles give larger amplitude fluctuations in the microlensing light curves. Consequently, current X-ray microlensing observations might have under or overestimated the sizes of the X-ray emission regions. We estimate this bias using a simple metric based on the amplitude of magnitude fluctuations. The half light radius of the X-ray emission region can be underestimated up to ~50% or overestimated up to ~20%. Underestimates are found in most situations we have investigated. The only exception is for a disk with large spin, radially flat emission profile, and observed nearly face on, where an overestimate is found. Thus, more accurate microlensing size constraints should be obtainable by including Kerr lensing. The caustic crossing time can differ by months after including Kerr strong gravity. A simultaneous monitoring of gravitational lensed quasars in both X-ray and optical bands with densely sampled X-ray light curves might reveal this feature. We conclude that it should be possible to constrain important parameters such as inclination angles and black hole spins from combined Kerr and microlensing effects.

Recent quasar microlensing observations have constrained the X-ray emission sizes of quasars to be about 10 gravitational radii, one order of magnitude smaller than the optical emission sizes. Using a new ray-tracing code for the Kerr space-time, we find that the observed X-ray flux is strongly influenced by the gravity field of the central black hole, even for observers at moderate inclination angles. We calculate inclination-dependent flux profiles of active galactic nuclei in the optical and X-ray bands by combining the Kerr lensing and projection effects for future references. We further study the dependence of the X-ray-to-optical flux ratio on the inclination angle caused by differential lensing distortion of the X-ray and optical emission, assuming several corona geometries. The strong lensing X-ray-to-optical magnification ratio can change by a factor of ~10 for normal quasars in some cases, and another factor of ~10 for broad absorption line quasars (BALs) and obscured quasars. Comparing our results with the observed distributions in normal and broad absorption line quasars, we find that the inclination angle dependence of the magnification ratios can change the X-ray-to-optical flux ratio distributions significantly. In particular, the mean value of the spectrum slope parameter $\alpha_{ox},$ $0.3838\log F_{2 keV}/F_{2500 {\AA}}$, can differ by ~0.1-0.2 between normal and broad absorption line quasars, depending on corona geometries, suggesting larger intrinsic absorptions in BALs.

We present Hubble Space Telescope (HST) ACS and WFC3 observations of SDSS J1029+2623, a three-image quasar lens system produced by a foreground cluster at z=0.584. Our strong lensing analysis reveals 6 additional multiply imaged galaxies. We confirm the complex nature of the mass distribution of the lensing cluster, with a bimodal distribution which deviates from the Chandra X-ray SB distribution. The Einstein radius is estimated to be \theta_E=15.2" \pm 0.5" for the quasar redshift of z=2.197. We derive a radial mass distribution from the combination of strong lensing, HST/ACS weak lensing, and Subaru/Suprime-cam weak lensing analysis results, finding a best-fit virial mass of M_vir=(1.5+0.40-0.35) \times 10^14 h^-1 M_sun and a concentration parameter of c_vir=25.7+14.1-7.5. The lensing mass estimate at the outer radius is smaller than the X-ray mass estimate by a factor of ~2. We ascribe this large mass discrepancy to shock heating of the intracluster gas during a merger, which is also suggested by the complex mass and gas distributions and the high value of the concentration parameter. In the HST image, we also identify a probable galaxy, GX, in the vicinity of the faintest quasar image C. The inclusion of GX explains the anomalous flux ratios between the quasar images. The morphology of the highly elongated quasar host galaxy is also well reproduced. The best-fit model suggests large total magnifications of 30 for the quasar and 35 for the quasar host galaxy, and has an AB time delay consistent with the measured value. (Abridged)

We use the progenitor of SN2012aw to illustrate the consequences of modeling circumstellar dust using Galactic (interstellar) extinction laws that (1) ignore dust emission in the near-IR and beyond; (2) average over dust compositions, and (3) mis-characterize the optical/UV absorption by assuming that scattered photons are lost to the observer. The primary consequences for the progenitor of SN2012aw are that both the luminosity and the absorption are significantly over-estimated. In particular, the stellar luminosity is most likely in the range 10^4.8 < L/Lsun < 10^5.0 and the star was not extremely massive for a Type IIP progenitor, with M < 15Msun. Given the properties of the circumstellar dust and the early X-ray/radio detections of SN2012aw, the star was probably obscured by an on-going wind with Mdot ~ 10^-5.5 to 10^-5.0 Msun/year at the time of the explosion, roughly consistent with the expected mass loss rates for a star of its temperature (T_* ~ 3600K) and luminosity. In the spirit of Galactic extinction laws, we supply simple interpolation formulas for circumstellar extinction by dusty graphitic and silicate shells as a function of wavelength (>0.3 micron) and total (absorption plus scattering) V-band optical depth (tau < 20). These do not include the contributions of dust emission, but provide a simple, physical alternative to incorrectly using interstellar extinction laws.

We present an analysis of selection biases in the M-sigma relation using Monte- Carlo simulations including the sphere of influence resolution selection bias and a selection bias in the velocity dispersion distribution. We find that the sphere of influence selection bias has a significant effect on the measured slope of the M-sigma relation, modeled as \beta_intrinsic = -4.69 + 2.22\beta_measured, where the measured slope is shallower than the model slope in the parameter range of \beta > 4, with larger corrections for steeper model slopes. Therefore, when the sphere of influence is used as a criterion to exclude unreliable measurements, it also in- troduces a selection bias that needs to be modeled to restore the intrinsic slope of the relation. We find that the selection effect due to the velocity dispersion distribution of the sample, which might not follow the overall distribution of the population, is not important for slopes of \beta ~ 4-6 of a logarithmically linear M-sigma relation, which could impact some studies that measure low (e.g., \beta < 4) slopes. Combining the selection biases in velocity dispersions and the sphere of influence cut, we find the uncertainty of the slope is larger than the value without modeling these effects, and estimate an intrinsic slope of \beta = 5.28^{+0.84}_{-0.55}.

We present results from Chandra observations of the cluster lens SDSS J1029+2623 at z_l=0.58, which is a gravitationally lensed quasar with the largest known image separation. We clearly detect X-ray emission both from the lensing cluster and the three lensed quasar images. The cluster has an X-ray temperature of kT = 8.1 (+2.0, -1.2) keV and bolometric luminosity of L_X = 9.6e44 erg s^-1. Its surface brightness is centered near one of the brightest cluster galaxies, and it is elongated East-West. We identify a subpeak North-West of the main peak, which is suggestive of an ongoing merger. Even so, the X-ray mass inferred from the hydrostatic equilibrium assumption appears to be consistent with the lensing mass from the Einstein radius of the system. We find significant absorption in the soft X-ray spectrum of the faintest quasar image, which can be caused by an intervening material at either the lens or source redshift. The X-ray flux ratios between the quasar images (after correcting for absorption) are in reasonable agreement with those at optical and radio wavelengths, and all the flux ratios are inconsistent with those predicted by simple mass models. This implies that microlensing effect is not significant for this system and dark matter substructure is mainly responsible for the anomalous flux ratios.

We observed SN2002bu in the near-IR with the Hubble Space Telescope, the mid-IR with the Spitzer Space Telescope and in X-rays with Swift 10 years after the explosion. If the faint L_H\sim100 Lsun HST near-IR source at the transient position is the near-IR counterpart of SN2002bu, then the source has dramatically faded between 2004 and 2012, from L\sim10^6.0 Lsun to L\sim10^4.5 Lsun. It is still heavily obscured, tau_V\sim5 in graphitic dust models, with almost all the energy radiated in the mid-IR. The radius of the dust emission is increasing as R\simt^(0.7+/-0.4) and the optical depth is dropping as tau_V\simt^(-1.3+/-0.4). The evolution expected for an expanding shell of material, tau_V\sim1/t^2, is ruled out at approximately 2 sigma while the tau_V\simt^(-0.8) to t^(-1) optical depth scaling for a shock passing through a pre-existing wind is consistent with the data. If the near-IR source is a chance superposition, the present day source can be moderately more luminous, significantly more obscured and evolving more slowly. While we failed to detect X-ray emission, the X-ray flux limits are consistent with the present day emissions being powered by an expanding shock wave. SN2002bu is clearly a member of the SN2008S class of transients, but continued monitoring of the evolution of the spectral energy distribution is needed to conclusively determine the nature of the transient.

We give analytic expressions for image properties of objects seen around point mass lenses embedded in a flat $\Lambda$CDM universe. An embedded lens in an otherwise homogeneous universe offers a more realistic representation of the lens's gravity field and its associated deflection properties than does the conventional linear superposition theory. Embedding reduces the range of the gravitational force acting on passing light beams thus altering all quantities such as deflection angles, amplifications, shears and Einstein ring sizes. Embedding also exhibits the explicit effect of the cosmological constant on these same lensing quantities. In this paper we present these new results and demonstrate how they can be used. The effects of embedding on image properties, although small i.e., usually less than a fraction of a percent, have a more pronounced effect on image distortions in weak lensing where the effects can be larger than 10%. Embedding also introduces a negative surface mass density for both weak and strong lensing, a quantity altogether absent in conventional Schwarzschild lensing. In strong lensing we find only one additional quantity, the potential part of the time delay, which differs from conventional lensing by as much as 4%, in agreement with our previous numerical estimates.

We present four new seasons of optical monitoring data and six epochs of X-ray photometry for the doubly-imaged lensed quasar Q J0158-4325. The high-amplitude, short-period microlensing variability for which this system is known has historically precluded a time delay measurement by conventional methods. We attempt to circumvent this limitation by application of a Monte Carlo microlensing analysis technique, but we are only able to prove that the delay must have the expected sign (image A leads image B). Despite our failure to robustly measure the time delay, we successfully model the microlensing at optical and X-ray wavelengths to find a half light radius for soft X-ray emission log(r_{1/2,X,soft}/cm) = 14.3^{+0.4}_{-0.5}, an upper limit on the half-light radius for hard X-ray emission log(r_{1/2,X,hard}/cm) <= 14.6 and a refined estimate of the inclination-corrected scale radius of the optical R-band (rest frame 3100 Angstrom) continuum emission region of log(r_s/cm) = 15.6+-0.3.

We present results from monitoring observations of the gravitationally lensed quasar RX J1131-1231 performed with the Chandra X-ray Observatory. The X-ray observations were planned with relatively long exposures that allowed a search for energy-dependent microlensing in the soft (0.2-2 keV) and hard (2-10 keV) light curves of the images of RX J1131-1231. We detect significant microlensing in the X-ray light-curves of images A and D, and energy-dependent microlensing of image D. The magnification of the soft band appears to be larger than that in the hard band by a factor of ~ 1.3 when image D becomes more magnified. This can be explained by the difference between a compact, softer-spectrum corona that is producing a more extended, harder spectrum reflection component off the disk. This is supported by the evolution of the fluorescent iron line in image D over three consecutive time-averaged phases of the light curve. In the first period, a Fe line at E = 6.36(-0.16,+0.13) keV is detected (at > 99% confidence). In the second period, two Fe lines are detected, one at E = 5.47(-0.08,+0.06) keV (detected at > 99% confidence) and another at E = 6.02(-0.07,+0.09) keV (marginally detected at > 90% confidence), and in the third period, a broadened Fe line at 6.42(-0.15,+0.19) keV is detected (at > 99% confidence). This evolution of the Fe line profile during the microlensing event is consistent with the line distortion expected when a caustic passes over the inner disk where the shape of the fluorescent Fe line is distorted by General Relativistic and Doppler effects.

We present \emph{Chandra} monitoring data for six gravitationally lensed quasars: QJ 0158$-$4325, HE 0435$-$1223, HE 1104$-$1805, SDSS 0924+0219, SDSS 1004+4112, and Q 2237+0305. We detect X-ray microlensing variability in all six lenses with high confidence. We detect energy dependent microlensing in HE 0435$-$1223, SDSS 1004+4112, SDSS 0924+0219 and Q 2237+0305. We present a detailed spectral analysis for each lens, and find that simple power-law models plus Gaussian emission lines give good fits to the spectra. We detect intrinsic spectral variability in two epochs of Q 2237+0305. We detect differential absorption between images in four lenses. We also detect the \feka\ emission line in all six lenses, and the Ni XXVII K$\alpha$ line in two images of Q 2237+0305. The rest frame equivalent widths of the \feka\ lines are measured to be 0.4--1.2 keV, significantly higher than those measured in typical active galactic nuclei of similar X-ray luminosities. This suggests that the \feka\ emission region is more compact or centrally concentrated than the continuum emission region.

We present our XMM-Newton observation of the fastest rotating spiral galaxy UGC 12591. We detect hot gas halo emission out to 110 kpc from the galaxy center, and constrain the halo gas mass to be smaller than 3.5e11 solar masses. We also measure the temperature of the hot gas as T=0.64\pm0.03 keV. Combining our X-ray constraints and the near-infrared and radio measurements in the literature, we find a baryon mass fraction of 0.03-0.04 in UGC 12591, suggesting a missing baryon mass of 75% compared with the cosmological mean value. Combined with another recent measurement in NGC 1961, the result strongly argues that the majority of missing baryons in spiral galaxies does not reside in their hot halos. We also find that UGC 12591 lies significantly below the baryonic Tully-Fisher relationship. Finally, we find that the baryon fractions of massive spiral galaxies are similar to those of galaxy groups with similar masses, indicating that the baryon loss is ultimately controlled by the gravitational potential well. The cooling radius of this gas halo is small, similar to NGC 1961, which argues that the majority of stellar mass of this galaxy is not assembled as a result of cooling of this gas halo.

We report results from a combined analysis of solar neutrino data from all phases of the Sudbury Neutrino Observatory. By exploiting particle identification information obtained from the proportional counters installed during the third phase, this analysis improved background rejection in that phase of the experiment. The combined analysis resulted in a total flux of active neutrino flavors from 8B decays in the Sun of (5.25 \pm 0.16(stat.)+0.11-0.13(syst.))\times10^6 cm^{-2}s^{-1}. A two-flavor neutrino oscillation analysis yielded \Deltam^2_{21} = (5.6^{+1.9}_{-1.4})\times10^{-5} eV^2 and tan^2{\theta}_{12}= 0.427^{+0.033}_{-0.029}. A three-flavor neutrino oscillation analysis combining this result with results of all other solar neutrino experiments and the KamLAND experiment yielded \Deltam^2_{21} = (7.41^{+0.21}_{-0.19})\times10^{-5} eV^2, tan^2{\theta}_{12} = 0.446^{+0.030}_{-0.029}, and sin^2{\theta}_{13} = (2.5^{+1.8}_{-1.5})\times10^{-2}. This implied an upper bound of sin^2{\theta}_{13} < 0.053 at the 95% confidence level (C.L.).

This paper details the solar neutrino analysis of the 385.17-day Phase-III data set acquired by the Sudbury Neutrino Observatory (SNO). An array of $^3$He proportional counters was installed in the heavy-water target to measure precisely the rate of neutrino-deuteron neutral-current interactions. This technique to determine the total active $^8$B solar neutrino flux was largely independent of the methods employed in previous phases. The total flux of active neutrinos was measured to be $5.54^{+0.33}_{-0.31}(stat.)^{+0.36}_{-0.34}(syst.)\times 10^{6}$ cm$^{-2}$ s$^{-1}$, consistent with previous measurements and standard solar models. A global analysis of solar and reactor neutrino mixing parameters yielded the best-fit values of $\Delta m^2 = 7.59^{+0.19}_{-0.21}\times 10^{-5}{eV}^2$ and $\theta = 34.4^{+1.3}_{-1.2}$ degrees.

We present our long term Chandra X-ray monitoring data for the gravitationally lensed quasar Q2237+0305 with 20 epochs spanning 10 years. We easily detect microlensing variability between the images in the full (0.2--8 keV), soft (0.2--2 keV), and hard (2--8 keV) bands at very high confidence. We also detect, for the first time, chromatic microlensing differences between the soft and hard X-ray bands. The hard X-ray band is more strongly microlensed than the soft band, suggesting that the corona above the accretion disk thought to generate the X-rays has a non-uniform electron distribution, in which the hotter and more energetic electrons occupy more compact regions surrounding the black holes. Both the hard and soft X-ray bands are more strongly microlensed than the optical (restframe UV) emission, indicating that the X-ray emission is more compact than the optical, confirming the microlensing results from other lenses.

We give the lens equation for light deflections caused by point mass condensations in an otherwise spatially homogeneous and flat universe. We assume the signal from a distant source is deflected by a single condensation before it reaches the observer. We call this deflector an embedded lens because the deflecting mass is part of the mean density. The embedded lens equation differs from the conventional lens equation because the deflector mass is not simply an addition to the cosmic mean. We prescribe an iteration scheme to solve this new lens equation and use it to compare our results with standard linear lensing theory. We also compute analytic expressions for the lowest order corrections to image amplifications and distortions caused by incorporating the lensing mass into the mean.

We present Suzaku observations of three iron low-ionization broad absorption line quasars (FeLoBALs). We detect J1723+5553 (3\sigma) in the observed 2--10 keV band, and constrain its intrinsic nh column density to nh > 6\times10^{23} \cmsq by modeling its X-ray hardness ratio. We study the broadband spectral index, aox, between the X-ray and UV bands by combining the X-ray measurements and the UV flux extrapolated from 2MASS magnitudes, assuming a range of intrinsic column densities, and then comparing the aox values for the three FeLoBALs with those from a large sample of normal quasars. We find that the FeLoBALs are consistent with the spectral energy distribution (SED) of normal quasars if the intrinsic nh column densities are nh > 7\times10^{23} \cmsq for J0943+5417, nh > 2\times10^{24} \cmsq for J1352+4293, and 6\times10^{23} < \nh < 3\times10^{24} \cmsq for J1723+5553. At these large intrinsic column densities, the optical depth from Thomson scattering can reach \sim6, which will significantly modulate the UV flux. Our results suggest that the X-ray absorbing material is located at a different place from the UV absorbing wind, likely between the X-ray and UV emitting regions. We find a significant kinetic feedback efficiency for FeLoBALs, indicating that the outflows are an important feedback mechanism in quasars.

Low-ionization (MgII, FeII, FeIII) broad absorption line quasars (LoBALs) probe a relatively obscured quasar population, and could be at an early evolutionary stage for quasars. We study the intrinsic fractions of LoBALs using the SDSS, 2MASS, and FIRST surveys. We find that the LoBAL fractions of the near infra-red (NIR) and radio samples are approximately 5--7 times higher than those measured in the optical sample. This suggests that the fractions measured in the NIR and radio bands are closer to the intrinsic fractions of the populations, and that the optical fractions are significantly biased due to obscuration effects, similar to high-ionization broad absorption line quasars (HiBALs). We also find that the LoBAL fractions decrease with increasing radio luminosities, again, similar to HiBALs. In addition, we find tentative evidence for high fractions of LoBALs at high NIR luminosities, especially for FeLoBALs with a fraction of ~18 per cent at M_K_s < -31 mag. This population of NIR luminous LoBALs may be at an early evolutionary stage of quasar evolution. We use a two-component model of LoBALs including a pure geometric component and a luminosity dependent component at high NIR luminosities, and obtain better fits than those from a pure geometric model. Therefore, the LoBAL population can be modelled as a hybrid of both the geometric and evolutionary models, where the geometric component constitutes 3.4\pm0.3, 5.8\pm0.4, and 1.5\pm0.3 per cent of the quasar population for BI-LoBALs, AI-LoBALs, and FeLoBALs, respectively. Considering a population of obscured quasars that do not enter the SDSS survey, which could have a much higher LoBAL fraction, we expect that intrinsic fraction of LoBALs could be even higher.

We report new results obtained in calibrations of superheated liquid droplet detectors used in dark matter searches with different radiation sources (n,$\alpha$,$\gamma$). In particular, detectors were spiked with alpha-emitters located inside and outside the droplets. It is shown that the responses are different, depending on whether alpha particles or recoil nuclei create the signals. The energy thresholds for $\alpha$-emitters are compared with test beam measurements using mono-energetic neutrons, as well as with theoretical predictions. Finally a model is presented which describes how the observed intensities of particle induced acoustic signals can be related to the dynamics of bubble growth in superheated liquids. An improved understanding of the bubble dynamics is an important first step in obtaining better discrimination between particle types interacting in detectors of this kind.

We present the baryon fractions of 2MASS groups and clusters as a function of cluster richness using total and gas masses measured from stacked ROSAT X-ray data and stellar masses estimated from the infrared galaxy catalogs. We detect X-ray emission even in the outskirts of clusters, beyond r_200 for richness classes with X-ray temperatures above 1 keV. This enables us to more accurately determine the total gas mass in these groups and clusters. We find that the optically selected groups and clusters have flatter temperature profiles and higher stellar-to-gas mass ratios than the individually studied, X-ray bright clusters. We also find that the stellar mass in poor groups with temperatures below 1 keV is comparable to the gas mass in these systems. Combining these results with individual measurements for clusters, groups, and galaxies from the literature, we find a break in the baryon fraction at ~1 keV. Above this temperature, the baryon fraction scales with temperature as f_b \propto T^0.20\pm0.03. We see significantly smaller baryon fractions below this temperature, and the baryon fraction of poor groups joins smoothly onto that of systems with still shallower potential wells such as normal and dwarf galaxies where the baryon fraction scales with the inferred velocity dispersion as f_b \propto \sigma^1.6. The small scatter in the baryon fraction at any given potential well depth favors a universal baryon loss mechanism and a preheating model for the baryon loss. The scatter is, however, larger for less massive systems. Finally, we note that although the broken power-law relation can be inferred from data points in the literature alone, the consistency between the baryon fractions for poor groups and massive galaxies inspires us to fit the two categories of objects (galaxies and clusters) with one relation.

We compute time delays for gravitational lensing in a flat LambdaCDM Swiss cheese universe. We assume a primary and secondary pair of light rays are deflected by a single point mass condensation described by a Kottler metric (Schwarzschild with Lambda) embedded in an otherwise homogeneous cosmology. We find that the cosmological constant's effect on the difference in arrival times is non-linear and at most around 0.002% for a large cluster lens; however, we find differences from time delays predicted by conventional linear lensing theory that can reach ~4% for these large lenses. The differences in predicted delay times are due to the failure of conventional lensing to incorporate the lensing mass into the mean mass density of the universe.

We compute the deflection angle to order (m/r_0)^2 and m/r_0*Lambda r_0^2 for a light ray traveling in a flat LambdaCDM cosmology which encounters a completely condensed mass region. We use a Swiss cheese model for the inhomogeneities and find that the most significant correction to the Einstein angle occurs not because of the non-linear terms but instead occurs because the condensed mass is embedded in a background cosmology.The Swiss cheese model predicts a decrease in the deflection angle of ~2% for weakly lensed galaxies behind the rich cluster A1689, and that the reduction can be as large as ~5% for similar rich clusters at z\approx 1. Weak lensing deflection angles caused by galaxies can likewise be reduced by as much as ~4%. We show that the lowest order correction in which Lambda appears is proportional to m/r_0*\sqrt{Lambda r_0^2} and could cause as much as a ~0.02% increase in the deflection angle for light that passes through a rich cluster. The lowest order non-linear correction in the mass is proportional to m/r_0*\sqrt{m/r_0} and can increase the deflection angle by ~0.005% for weak lensing by galaxies.

The metallicity in galaxy clusters is expected to originate from the stars in galaxies, with a population dominated by high mass stars likely being the most important stellar component, especially in rich clusters. We examine the relationship between the metallicity and the prominence of galaxies as measured by the star to baryon ratio, M$_*$/M$_{bary}$. Counter to expectations, we rule out a metallicity that is proportional to M$_*$/M$_{bary}$, where the best fit has the gas phase metallicity decreasing with M$_*$/M$_{bary}$, or the metallicity of the gas plus the stars being independent of M$_*$/M$_{bary}$. This implies that the population of stars responsible for the metals is largely proportional to the total baryonic mass of the cluster, not to the galaxy mass within the cluster. If generally applicable, most of the heavy elements in the universe were not produced within galaxies.

We discuss how the effective radius Phi(Re) function (ERF) recently worked out by Bernardi et al. (2009) represents a new testbed to improve the current understanding of Semi-analytic Models of Galaxy formation. In particular, we here show that a detailed hierarchical model of structure formation can broadly reproduce the correct peak in the size distribution of local early-type galaxies, although it significantly overpredicts the number of very compact and very large galaxies. This in turn is reflected in the predicted size-mass relation, much flatter than the observed one, due to too large (~3 kpc) low-mass galaxies (<10^11 \msun), and to a non-negligible fraction of compact (< 0.5-1 kpc) and massive galaxies (> 10^11 \msun). We also find that the latter discrepancy is smaller than previously claimed, and limited to only ultracompact (Re < 0.5 kpc) galaxies when considering elliptical-dominated samples. We explore several causes behind these effects. We conclude that the former problem might be linked to the initial conditions, given that large and low-mass galaxies are present at all epochs in the model. The survival of compact and massive galaxies might instead be linked to their very old ages and peculiar merger histories. Overall, knowledge of the galactic stellar mass {\em and} size distributions allows a better understanding of where and how to improve models.

In a sample of about 45,700 early-type galaxies extracted from SDSS, we find that the shape, normalization, and dispersion around the mean size-stellar mass relation is the same for young and old systems, provided the stellar mass is greater than 3*10^10 Msun. This is difficult to reproduce in pure passive evolution models, which generically predict older galaxies to be much more compact than younger ones of the same stellar mass. However, this aspect of our measurements is well reproduced by hierarchical models of galaxy formation. Whereas the models predict more compact galaxies at high redshifts, subsequent minor, dry mergers increase the sizes of the more massive objects, resulting in a flat size-age relation at the present time. At lower masses, the models predict that mergers are less frequent, so that the expected anti-correlation between age and size is not completely erased. This is in good agreement with our data: below 3*10^10 Msun, the effective radius R_e is a factor of ~2 lower for older galaxies. These successes of the models are offset by the fact that the predicted sizes have other serious problems, which we discuss.

We cross-correlate the SDSS DR3 quasar sample with FIRST and the Vestergaard et al. black hole (BH) mass sample to compare the mean accretion histories of optical and radio quasars. We find significant statistical evidence that radio quasars have a higher mean Eddington ratio Lambda at z > 2 with respect to optical quasars, while the situation is clearly reverse at z < 1. At z > 2 radio quasars happen to be less massive than optical quasars; however, as redshift decreases radio quasars appear in increasingly more massive BHs with respect to optical quasars. These two trends imply that radio sources are not a mere random subsample of optical quasars. No clear correlation between radio activity and BH mass and/or accretion rate is evident from our data, pointing to other BH properties, possibly the spin, as the driver of radio activity. We have checked that our main results do not depend on any evident bias. We perform detailed modelling of reasonable accretion histories for optical and radio quasars, finding that radio quasars grow by a factor of a few, at the most, since z ~ 4. The comparison between the predicted mass function of active radio quasars and the observed optical luminosity function of radio quasars, implies a significantly lower probability for lower mass BHs to be radio loud at all epochs, in agreement with what is observed in the local universe.

We use gravitational microlensing of the four images of the z=0.658 quasar RXJ1131-1231 to measure the sizes of the optical and X-ray emission regions of the quasar. The (face-on) scale length of the optical disk at rest frame 400 nm is 1.3 10^15cm, while the half-light radius of the rest frame 0.3-17 keV X-ray emission is 2.3 10^14cm. The formal uncertainties are factors of 1.6 and 2.0, respectively. With the exception of the lower limit on the X-ray size, the results are very stable against any changes in the priors used in the analysis. Based on the Hbeta line-width, we estimate that the black hole mass is ~10^8 Msun, which corresponds to a gravitational radius of r_g~2 10^13 cm. Thus, the X-ray emission is emerging on scales of ~10r_g and the 400 nm emission on scales of ~70 r_g. A standard thin disk of this size should be significantly brighter than observed. Possible solutions are to have a flatter temperature profile or to scatter a large fraction of the optical flux on larger scales after it is emitted. While our calculations were not optimized to constrain the dark matter fraction in the lens galaxy, dark matter dominated models are favored. With well-sampled optical and X-ray light curves over a broad range of frequencies there will be no difficulty in extending our analysis to completely map the structure of the accretion disk as a function of wavelength.

Using the sample of long Gamma-ray bursts (GRBs) detected by Swift-BAT before June 2007, we measure the cumulative distribution of the peak photon fluxes (log N - log P) of the Swift bursts. Compared with the BATSE sample, we find that the two distributions are consistent after correcting the band pass difference, suggesting that the two instruments are sampling the same population of bursts. We also compare the log N - log P distributions for sub-samples of the Swift bursts, and find evidence for a deficit (99.75% confident) of dark bursts without optical counterparts at high peak flux levels, suggesting different redshift or Gamma-ray luminosity distributions for these bursts. The consistency between the log N - log P distributions for the optically detected bursts with and without redshift measurements indicates that the current sample of the Swift bursts with redshift measurements, although selected heterogeneously, represents a fair sample of the non-dark bursts. We calculate the luminosity functions of this sample in two redshift bins (z<1 and z>1), and find a broken power-law is needed to fit the low redshift bin, where dN/dL \propto L^{-1.27\pm0.06} for high luminosities (L_{peak} > 5e48 erg/s) and dN/dL \propto L^{-2.3\pm0.3} at for low luminosities, confirming the results of several studies for a population of low luminosity GRBs.

We analyzed Chandra observations of three gravitational lenses, SBS0909+523, FBQS0951+2635, and B1152+199, to measure the differential X-ray absorption and the dust-to-gas ratio of the lens galaxies. We successfully detected the differential X-ray absorption in SBS0909+523 and B1152+199, and failed to detect it in FBQS0951+2635 due to the dramatic drop in its flux from the ROSAT epoch. These measurements significantly increase the sample of dust-to-gas ratio measurements in cosmologically-distant, normal galaxies. Using the larger sample, we obtain an average dust-to-gas ratio of E(B-V)/NH = (1.5\pm0.5)e-22 mag cm^2/atoms with an estimated intrinsic dispersion in the ratio of \simeq 40%. This average dust-to-gas ratio is consistent with our previous measurement, and the average Galactic value of 1.7e-22 mag cm^2/atoms and the estimated intrinsic dispersion is also consistent with the 30% observed in the Galaxy. A larger sample size is still needed to improve the measurements and to begin studying the evolution in the ratio with cosmic time. We also detected X-ray microlensing in SBS0909+523 and significant X-ray variability in FBQS0951+2635.

Using the 2.4m MDM and 8.4m Large Binocular Telescope, we observed nine GRB afterglows to systematically probe the late time behaviors of afterglows including jet breaks, flares, and supernova bumps. In particular, the LBT observations have typical flux limits of 25-26 mag in the Sloan r' band, which allows us to extend the temporal baseline for measuring jet breaks by another decade in time scale. We detected four jet breaks (including a "textbook" jet break in GRB070125) and a fifth candidate, all of which are not detectable without deep, late time optical observations. In the other four cases, we do not detect the jet breaks either because of contamination from the host galaxy light, the presence of a supernova bump, or the intrinsic faintness of the optical afterglow. This suggests that the basic picture that GRBs are collimated is still valid and that the apparent lack of Swift jet breaks is due to poorly sampled afterglow light curves, particularly at late times. Besides the jet breaks, we also detected late time flares, which could attribute to late central engine activities, and two supernova bumps.

We present X-ray spectral analysis of five Chandra and XMM-Newton observations of the gravitationally-lensed blazar PKS1830-211 from 2000 to 2004. We show that the X-ray absorption toward PKS1830-211 is variable, and the variable absorption is most likely to be intrinsic with amplitudes of about 2-30e22 cm^-2 depending on whether or not the absorber is partially covering the X-ray source. Our results confirm the variable absorption observed previously, although interpreted differently, in a sequence of ASCA observations. This large variation in the absorption column density can be interpreted as outflows from the central engine in the polar direction, consistent with recent numerical models of inflow/outflows in AGNs. In addition, it could possibly be caused by the interaction between the blazar jet and its environment, or the variation from the geometric configuration of the jet. While the spectra can also be fitted with a variable absorption at the lens redshift, we show that this model is unlikely. We also rule out the simple microlensing interpretation of variability which was previously suggested.

Using the 8.4m Large Binocular Telescope, we observed six GRB afterglows from 2.8 hours to 30.8 days after the burst triggers to systematically probe the late time behaviors of afterglows including jet breaks, flares, and supernova bumps. We detected five afterglows with Sloan r' magnitudes ranging from 23.0-26.3 mag. The depth of our observations allows us to extend the temporal baseline for measuring jet breaks by another decade in time scale. We detected two jet breaks and a third candidate, all of which are not detectable without deep, late time optical observations. In the other three cases, we do not detect the jet breaks either because of contamination from the host galaxy light, the presence of a supernova bump, or the intrinsic faintness of the optical afterglow. This suggests that the basic picture that GRBs are collimated is still valid and that the apparent lack of Swift jet breaks is due to poorly sampled afterglow light curves, particularly at late times.

We present galaxy luminosity functions at 3.6, 4.5, 5.8, and 8.0 micron measured by combining photometry from the IRAC Shallow Survey with redshifts from the AGN and Galaxy Evolution Survey of the NOAO Deep Wide-Field Survey Bootes field. The well-defined IRAC samples contain 3800-5800 galaxies for the 3.6-8.0 micron bands with spectroscopic redshifts and z < 0.6. We obtained relatively complete luminosity functions in the local redshift bin of z < 0.2 for all four IRAC channels that are well fit by Schechter functions. We found significant evolution in the luminosity functions for all four IRAC channels that can be fit as an evolution in M* with redshift, \Delta M* = Qz. While we measured Q=1.2\pm0.4 and 1.1\pm0.4 in the 3.6 and 4.5 micron bands consistent with the predictions from a passively evolving population, we obtained Q=1.8\pm1.1 in the 8.0 micron band consistent with other evolving star formation rate estimates. We compared our LFs with the predictions of semi-analytical galaxy formation and found the best agreement at 3.6 and 4.5 micron, rough agreement at 8.0 micron, and a large mismatch at 5.8 micron. These models also predicted a comparable Q value to our luminosity functions at 8.0 micron, but predicted smaller values at 3.6 and 4.5 micron. We also measured the luminosity functions separately for early and late-type galaxies. While the luminosity functions of late-type galaxies resemble those for the total population, the luminosity functions of early-type galaxies in the 3.6 and 4.5 micron bands indicate deviations from the passive evolution model, especially from the measured flat luminosity density evolution. Combining our estimates with other measurements in the literature, we found (53\pm18)% of the present stellar mass of early-type galaxies has been assembled at z=0.7.

We present results from a monitoring campaign performed with the Chandra X-ray Observatory of the gravitationally lensed quasars RX J1131-1231 and HE 1104-1805. We detect significant X-ray variability in all images of both quasars. The flux variability detected in image A of RX J1131-1231 is of particular interest because of its high amplitude (a factor of ~ 20). We interpret it as arising from microlensing since the variability is uncorrelated with that of the other images and the X-ray flux ratios show larger changes than the optical as we would expect for microlensing of the more compact X-ray emission regions. The differences between the X-ray and optical flux ratios of HE 1104-1805 are less dramatic, but there is no significant soft X-ray or dust absorption, implying the presence of X-ray microlensing in this system as well. Combining the X-ray data with the optical light curves we find that the X-ray emitting region of HE 1104-1805 is compact with a half-light radius ~ 6 r_g, where the gravitational radius is r_g = 3.6 x 10^14 cm, thus placing significant constraints on AGN corona models. We also find that the microlensing in HE 1104-1805 favors mass models for the lens galaxy that are dominated by dark matter. Finally, we better characterize the massive foreground cluster near RX J1131-1231, set limits on other sources of extended X-ray emission, and limit the fluxes of any central odd images to be 30-50 (3 sigma) times fainter than the observed images.

We find that the fraction of classical Broad Absorption Line quasars (BALQSOs) among the FIRST radio sources in the Sloan Data Release 3, is 20.5^{+7.3}_{-5.9}% at the faintest radio powers detected (L_{\rm 1.4 GHz}~10^{32} erg/s), and rapidly drops to <8% at L_{\rm 1.4 GHz}~3*10^{33} erg/s. Similarly, adopting the broader Absorption Index (AI) definition of Trump et al. (2006) we find the fraction of radio BALQSOs to be 44^{+8.1}_{-7.8}% reducing to 23.1^{+7.3}_{-6.1}% at high luminosities. While the high fraction at low radio power is consistent with the recent near-IR estimates by Dai et al. (2008), the lower fraction at high radio powers is intriguing and confirms previous claims based on smaller samples. The trend is independent of the redshift range, the optical and radio flux selection limits, or the exact definition of a radio match. We also find that at fixed optical magnitude, the highest bins of radio luminosity are preferentially populated by non-BALQSOs, consistent with the overall trend. We do find, however, that those quasars identified as AI-BALQSOs but \emph{not} under the classical definition, do not show a significant drop in their fraction as a function of radio power, further supporting independent claims for which these sources, characterized by lower equivalent width, may represent an independent class with respect to the classical BALQSOs. We find the balnicity index, a measure of the absorption trough in BALQSOs, and the mean maximum wind velocity to be roughly constant at all radio powers. We discuss several plausible physical models which may explain the observed fast drop in the fraction of the classical BALQSOs with increasing radio power, \emph{although no one is entirely satisfactory}. (abridged).

We report on multi-wavelength observations, ranging from the X-ray to radio wave bands, of the IPN-localized gamma-ray burst GRB 070125. Spectroscopic observations reveal the presence of absorption lines due to O I, Si II, and C IV, implying a likely redshift of z = 1.547. The well-sampled light curves, in particular from 0.5 to 4 days after the burst, suggest a jet break at 3.7 days, corresponding to a jet opening angle of ~7.0 degrees, and implying an intrinsic GRB energy in the 1 - 10,000 keV band of around E = (6.3 - 6.9)x 10^(51) erg (based on the fluences measured by the gamma-ray detectors of the IPN network). GRB 070125 is among the brightest afterglows observed to date. The spectral energy distribution implies a host extinction of Av < 0.9 mag. Two rebrightening episodes are observed, one with excellent time coverage, showing an increase in flux of 56% in ~8000 seconds. The evolution of the afterglow light curve is achromatic at all times. Late-time observations of the afterglow do not show evidence for emission from an underlying host galaxy or supernova. Any host galaxy would be subluminous, consistent with current GRB host-galaxy samples. Evidence for strong Mg II absorption features is not found, which is perhaps surprising in view of the relatively high redshift of this burst and the high likelihood for such features along GRB-selected lines of sight.

We analyzed the microlensing of the X-ray and optical emission of the lensed quasar PG 1115+080. We find that the effective radius of the X-ray emission is 1.3(+1.1 -0.5) dex smaller than that of the optical emission. Viewed as a thin disk observed at inclination angle i, the optical accretion disk has a scale length, defined by the point where the disk temperature matches the rest frame energy of the monitoring band (kT=hc/lambda_rest with lambda_rest=0.3 micron), of log[(r_{s,opt}/cm)(cos(i) / 0.5)^{1/2}] = 16.6 \pm 0.4. The X-ray emission region (1.4-21.8 keV in the rest frame) has an effective half-light radius of log[r_{1/2,X}/cm] = 15.6 (+0.6-0.9}. Given an estimated black hole mass of 1.2 * 10^9 M_sun, corresponding to a gravitational radius of log[r_g/cm] = 14.3, the X-ray emission is generated near the inner edge of the disk while the optical emission comes from scales slightly larger than those expected for an Eddington-limited thin disk. We find a weak trend supporting models with low stellar mass fractions near the lensed images, in mild contradiction to inferences from the stellar velocity dispersion and the time delays.

We present results from a statistical analysis of 173 bright radio-quiet AGNs selected from the Chandra Deep Field-North and Chandra Deep Field-South surveys (hereafter, CDFs) in the redshift range of 0.1 < z < 4. We find that the X-ray power-law photon index (Gamma) of radio-quiet AGNs is correlated with their 2-10 keV rest-frame X-ray luminosity (L_X) at the > 99.5 percent confidence level in two redshift bins, 0.3 < z < 0.96, and 1.5 < z < 3.3 and is slightly less significant in the redshift bin 0.96 < z < 1.5. We investigate the redshift evolution of the correlation between the power-law photon index and the hard X-ray luminosity and find that the slope and offset of a linear fit to the correlation change significantly (at the > 99.9 percent confidence level) between redshift bins of 0.3 < z < 0.96 and 1.5 < z < 3.3. We explore physical scenarios explaining the origin of this correlation and its possible evolution with redshift in the context of steady corona models focusing on its dependency on variations of the properties of the hot corona with redshift.

The intrinsic fraction of broad absorption line quasars (BALQSOs) is important in constraining geometric and evolutionary models of quasars. We present the fraction of BALQSOs in 2MASS detected quasars within the SDSS DR3 sample in the redshift range of 1.7 < z < 4.38. The fraction of BALQSOs is 40.4^{+3.4}_{-3.3}% in the 2MASS 99% database K_s band completeness sample, and 38.5^{+1.7}_{-1.7}% in the larger 2MASS sample extending below the completeness limit. These fractions are significantly higher than the 26% reported in the optical bands for the same parent sample. We also present the fraction of BALQSOs as functions of apparent magnitudes, absolute magnitudes, and redshift in the 2MASS and SDSS bands. The 2MASS fractions are consistently higher than the SDSS fractions in every comparison, and the BALQSO fractions steadily increase with wavelength from the SDSS u to the 2MASS K_s bands. Furthermore, the i - K_s color distributions of BALQSOs and non-BALQSOs indicate that BALQSOs are redder than non-BALQSOs, with a K-S test probability of 2e-12. These results are consistent with the spectral difference between BALQSOs and non-BALQSOs including both the absorption troughs and dust extinction in BALQSOs, which leads to significant selection biases against BALQSOs in the optical bands. Using a simple simulation incorporating the luminosity function of quasars and the amount of obscuration for BALQSOs, we simultaneously fit the BALQSO fractions in the SDSS and 2MASS bands. We obtain a true BALQSO fraction of 43\pm2% for luminous quasars (M_{K_s} \lesssim -30.1 mag).

We present results from Chandra and XMM-Newton observations of the low-ionization broad absorption line (LoBAL) quasar H 1413+117. Our spatial and spectral analysis of a recent deep Chandra observation confirms a microlensing event in a previous Chandra observation performed about 5 years earlier. We present constraints on the structure of the accretion flow in H 1413+117 based on the time-scale of this microlensing event. Our analysis of the combined spectrum of all the images indicates the presence of two emission peaks at rest-frame energies of 5.35 keV and 6.32 keV detected at the > 98% and > 99% confidence levels, respectively. The double peaked Fe emission line is fit well with an accretion-disk line model, however, the best-fitting model parameters are neither well constrained nor unique. Additional observations are required to constrain the model parameters better and to confirm the relativistic interpretation of the double peaked Fe Kalpha line. Another possible interpretation of the Fe emission is fluorescent Fe emission from the back-side of the wind. The spectra of images C and D show significant high-energy broad absorption features that extend up to rest-frame energies of 9 keV and 15 keV respectively. We propose that a likely cause of these differences is significant variability of the outflow on time-scales that are shorter than the time-delays between the images. The Chandra observation of H 1413+117 has made possible for the first time the detection of the inner regions of the accretion disk and/or wind and the high ionization component of the outflowing wind of a LoBAL quasar.

We obtained 98 R-band and 18 B, r', i' images of the optical afterglow of GRB 060526 (z=3.21) with the MDM 1.3m, 2.4m, and the PROMPT telescopes in Cerro Tololo over the 5 nights following the burst trigger. Combining these data with other optical observations reported in GCN and the Swift-XRT observations, we compare the optical and X-ray afterglow light curves of GRB 060526. Both the optical and X-ray afterglow light curves show rich features, such as flares and breaks. The densely sampled optical observations provide very good coverage at T>1.e4 sec. We observed a break at 2.4e5 sec in the optical afterglow light curve. Compared with the X-ray afterglow light curve, the break is consistent with an achromatic break supporting the beaming models of GRBs. However, the pre-break and post-break temporal decay slopes are difficult to explain in simple afterglow models. We estimated a jet angle of \theta_j ~ 7 degrees and a prompt emission size of R_{prompt} ~ 2e14 cm. In addition, we detected several optical flares with amplitudes of \Delta m ~ 0.2, 0.6, and 0.2 mag. The X-ray afterglows detected by Swift have shown complicated decay patterns. Recently, many well-sampled optical afterglows also show decays with flares and multiple breaks. GRB 060526 provides an additional case of such a complex, well observed optical afterglow. The accumulated well-sampled afterglows indicate that most of the optical afterglows are complex.

We stacked the X-ray data from the ROSAT All Sky Survey for over 4,000 clusters selected from the 2MASS catalog and divided into five richness classes. We detected excess X-ray emission over background at the center of the stacked images in all five richness bins. The interrelationships between the mass, X-ray temperature and X-ray luminosity of the stacked clusters agree well with those derived from catalogs of X-ray clusters. Poisson variance in the number of galaxies occupying halos of a given mass leads to significant differences between the average richness at fixed mass and the average mass at fixed richness that we can model relatively easily using a simple model of the halo occupation distribution. These statistical effects probably explain recent results in which optically-selected clusters lie on the same X-ray luminosity-temperature relations as local clusters but have lower optical richnesses than observed for local clusters with the same X-ray properties. When we further binned the clusters by redshift, we did not find significant redshift-dependent biases in the sense that the X-ray luminosities for massive clusters of fixed optical richness show little dependence on redshift beyond that expected from the effects of Poisson fluctuations. Our results demonstrate that stacking of RASS data from optically selected clusters can be a powerful test for biases in cluster selection algorithms.

We report on two recent z~4 gamma-ray bursts (GRBs), GRB 060206 and GRB 060210, for which we have obtained well-sampled optical light curves. Our data, combined with early optical data reported in the literature, shows unusual behavior for both afterglows. In R-band GRB 060206 (z=4.045) experienced a slow early decay, followed by a rapid increase in brightness by factor ~2.5 about 1 hour after the burst. Its afterglow then faded in a broken power-law fashion, with a smooth break at t_b=0.6 days, but with additional, less dramatic (~10%) ``bumps and wiggles'', well detected in the densely sampled light curve. The R-band afterglow of GRB 060210 (z=3.91) is also unusual: the light curves was more or less flat between 60 and 300 sec after the burst, followed by ~70% increase at ~600 sec after the burst, after which the light curve declined as a t^{-1.3} power-law. Despite earlier reports to the contrary, we find that for GRB 060206 X-rays also more or less follow the optical decay, but with significant variations on short timescales. However, the X-ray afterglow is contaminated by a nearby, variable source, which especially at late times obscures the behavior of the X-ray afterglow. We argue that ``anomalous'' optical afterglows are likely to be the norm, and that the rapid variations often seen in Swift-XRT data would also be seen in the optical light curves, given good enough sampling. As a result, some of the often employed procedures, such as deriving the jet opening angle using a smooth broken power-law fit to the optical light curves, in many cases might have a poor statistical significance. We argue that the early increase in brighness for both bursts might be due to the turn-on of the external shock. Existence of such features could provide valuable additional information about the burst. (Abridged)

Microlensing of gravitationally lensed quasars by the stars in the foreground lens galaxy can be used to probe the nature of dark matter, to determine the mean stellar mass in the lens galaxy, and to measure the internal structure of quasar accretion disks. Until recently, little progress has been made toward using microlensing for these purposes because of the difficulty in obtaining the necessary data and the lack of good analysis methods. In the last few years, both problems have been solved. In particular, Bayesian analysis methods provide a general approach to measuring quantities of physical interest and their uncertainties from microlensing light curves. We discuss the data and the analysis methods and show preliminary results for all three astrophysical applications.

We report the discovery of a cluster-scale lensed quasar, SDSS J1029+2623, selected from the Sloan Digital Sky Survey. The lens system exhibits two lensed images of a quasar at z_s=2.197. The image separation of 22.5" makes it the largest separation lensed quasar discovered to date. The similarity of the optical spectra and the radio loudnesses of the two components support the lensing hypothesis. Images of the field show a cluster of galaxies at z_l~0.55 that is responsible for the large image separation. The lensed images and the cluster light center are not collinear, which implies that the lensing cluster has a complex structure.

We report the detection of differential gas column densities in three gravitational lenses, MG0414+0534, HE1104-1805, and PKS1830-211. Combined with the previous differential column density measurements in B1600+434 and Q2237+0305 and the differential extinction measurements of these lenses, we probe the dust-to-gas ratio of a small sample of cosmologically distant normal galaxies. We obtain an average dust-to-gas ratio of E(B-V)/NH =(1.4\pm0.5) e-22 mag cm^2/atoms with an estimated intrinsic dispersion in the ratio of ~40%. This average dust-to-gas ratio is consistent with the average Galactic value of 1.7e-22 mag cm^2/atoms and the estimated intrinsic dispersion is also consistent with the 30% observed in the Galaxy.

We have measured the three time delays for the quadruple gravitational lens RX J1131-1231 using two seasons of monitoring data. The short delays between cusp images are A-B = 11.98 +1.52,-1.27 days and A-C = 9.61 +1.97,-1.57 days. The long A-D delay for the counter image is not as precisely determined because of the season gaps, but the data suggest a delay of -87 +- 8 days. The short delays are difficult to explain using standard isothermal halo models of the lensing potential, which instead prefer A-B and A-C delays of ~1 day for reasonable values of the Hubble constant. Matching the cusp delays is possible by adding a significant (~5 x 10^10 M_sun) amount of matter nearly coincident (~0.05" South-East) with the A image. Adding such a satellite also helps improve the quasar and lens astrometry of the model, reduces the velocity dispersion of the main lens and shifts it closer to the Fundamental Plane. This is suggestive of a satellite galaxy to the primary lens, but its expected luminosity and proximity to both image A and the system's bright Einstein ring make visual identification impossible even with the existing HST data. We also find evidence for significant structure along the line of sight toward the lens. Archival Chandra observations show two nearby regions of extended X-ray emission, each with bolometric X-ray luminosities of 2-3 x 10^43 ergs/s. The brighter region is located 153" from the lens and centered on a z=0.1 foreground cD galaxy, and the fainter and presumably more distant region is 4-5 times closer (in angular separation) to the lens and likely corresponds to the weaker of two galaxy red sequences (which includes the lens galaxy) previously detected at optical wavelengths.

We observed B1600+434 and B1608+656 with CXO/ACIS, detecting both quasar images in B1600+434 and three of four images in B1608+656. We did not detect significant X-ray emission from nearby galaxy groups or clusters associated with each lens galaxy. The upper limits on the X-ray luminosity of any cluster within 4 arcmin of each lens and at each lens redshift are of ~2\times10^{42} and ~6\times10^{42} erg/s for B1600+434 and B1608+656, respectively. The radio-loud source quasars have power-law photon indices of \Gamma=1.9\pm0.2 and \Gamma=1.4\pm0.3 and X-ray luminosities of 1.4^{+0.2}_{-0.1}\times10^{45} and 2.9^{+0.7}_{-0.4}\times10^{44} erg/s for B1600+434 and B1608+656, respectively before correcting for the magnification. We detected a differential absorption column density of \Delta N_H ~ 3\times10^{21} cm^{-2} between the two images of B1600+434, roughly consistent with expectations from differential extinction estimates of \Delta E(B-V)=0.1 mag and a standard dust-to-gas ratio. The differential absorption observed in gravitational lenses may serve as an important probe to study the gas content in high redshift galaxies since it can separate the absorbing column originating from the lens galaxy and those intrinsic to quasars. We also detected 157 serendipitous X-ray sources in the two Chandra fields and identified the brighter optical counterparts using the SDSS and DPOSS surveys.

The possibility that long gamma-ray burst (GRB) jets are structured receives growing attention recently, and we have suggested that most GRBs and their softer, less energetic fraternity, X-ray flashes (XRFs), can be understood within a quasi-universal structured jet picture, given that the jet structure for each individual burst is a Gaussian or similar function. Here we perform a global test on such a quasi-universal Gaussian-like structured jet by comparing Monte-Carlo simulation results with a broad spectrum of observational data. Using the same set of input parameters as in the previous work (Zhang et al. 2004), we confront the model with more observational constraints. These constraints include the burst redshift distribution, jet break angle distribution, two-dimensional redshift vs. jet break angle distribution, luminosity function, and Log N - Log P distribution. The results indicate that the model is generally compatible with the data. This conclusion, together with our previous tests with the observed jet break angle vs. isotropic energy and observed peak energy vs. fluence relations, suggests that current long GRB and XRF data are generally consistent with such a quasi-standard-energy and quasi-standard-angle jet picture. With future homogeneous burst samples (such as the one to be retrieved from the Swift mission), the refined GRB jet structure can be further constrained through a global comparison between various observed and predicted burst property distributions and relations.

We present results from a mini-survey of relatively high redshift (1.7 < z < 4) gravitationally lensed radio-quiet quasars observed with the Chandra X-ray Observatory and with XMM-Newton. The lensing magnification effect allows us to search for changes in quasar spectroscopic and flux variability properties with redshift over three orders of magnitude in intrinsic X-ray luminosity. It extends the study of quasar properties to unlensed X-ray flux levels as low as a few times 10^{-15} erg cm^{-2} s^{-1} in the observed 0.4-8 keV band. We find a possible correlation between the X-ray powerlaw photon index and X-ray luminosity of the gravitationally lensed radio-quiet quasar sample. The X-ray spectral slope steepens as the X-ray luminosity increases. This correlation is still significant when we combine our data with other samples of radio-quiet quasars with z > 1.5, especially in the low luminosity range between 10^{43}-10^{45.5} erg s^{-1}. This result is surprising considering that such a correlation is not found for quasars with redshifts below 1.5. We suggest that this correlation can be understood in the context of the hot-corona model for X-ray emission from quasar accretion disks, under the hypothesis that the quasars in our sample accrete very close to their Eddington limits and the observed luminosity range is set by the range of black hole masses (this hypothesis is consistent with recent predictions of semi-analytic models for quasar evolution). The upper limits of X-ray variability of our relatively high redshift sample of lensed quasars are consistent with the known correlation between variability and luminosity observed in Seyfert 1s when this correlation is extrapolated to the larger luminosities of our sample.

Two separate topics are discussed. (1) We describe the classifications of the long GRB early afterglow lightcurves within the framework of the fireball shock model, focusing on the interplay between the reverse and forward shock emission components. We will also provide evidence that the central engine of at least two bursts are entrained with strong magnetic fields, and discuss the implications of this result for our understanding of the GRB phenomenon; (2) We argue that the current gamma-ray burst (GRB) and X-ray flash (XRF) data are consistent with a picture that all GRB-XRF jets are structured and quasi-universal, with a typical Gaussian-like jet structure.

An observed correlation Ep \propto (E_{iso})^{1/2} extending from Gamma-ray Bursts (GRBs) to X-ray flashes (XRFs) poses problems both for a power-law universal jet model where the energy per solid angle decreases as the inverse square of the angle respect to the jet axis, and for a conical jet model with a uniform energy density within the jet beam and a sharp energy cut-off at the jet edge. Here we show that the current GRB-XRF prompt emission/afterglow data can be understood in terms of a picture in which the GRB-XRF jets are quasi-universal and structured, with a Gaussian-like or similar structure, i.e., one where the jet has a characteristic angle, with a mild variation of energy inside and a rapid (e.g. exponential) decrease of energy outside of it. A Monte Carlo simulation shows that the current data is compatible with such a quasi-universal Gaussian jet with a typical opening angle of 5.7(+3.4)(-2.1) degrees, and with a standard jet energy of about log (E_j/1 erg)= 51.1+-0.3. According to this model, the true-to-observed number ratio of the whole GRB-XRF population is about 14 with the current instrumental sensitivity.

The idea that GRBs originate from uniform jets has been used to explain numerous observations of breaks in the GRB afterglow lightcurves. We explore the possibility that GRBs instead originate from a structured jet that may be quasi-universal, where the variation in the observed properties of GRBs is due to the variation in the observer viewing angle. We test how various models reproduce the jet data of Bloom, Frail, & Kulkarni (2003), which show a negative correlation between the isotropic energy output and the inferred jet opening angle (in a uniform jet configuration). We find, consistent with previous studies, that a power-law structure for the jet energy as a function of angle gives a good description. However, a Gaussian jet structure can also reproduce the data well, particularly if the parameters of the Gaussian are allowed some scatter. We place limits on the scatter of the parameters in both the Gaussian and power-law models needed to reproduce the data, and discuss how future observations will better distinguish between these models for the GRB jet structure. In particular, the Gaussian model predicts a turnover at small opening angles and in some cases a sharp cutoff at large angles, the former of which may already have been observed. We also discuss the predictions each model makes for the observed luminosity function of GRBs and compare these predictions with the existing data.

We present the observations of the gravitationally lensed system QSO 2237+0305 performed with the Chandra X-ray Observatory on 2000 Sept. 6, and on 2001 Dec. 8 for 30.3 ks and 9.5 ks, respectively. Imaging analysis resolves the four X-ray images of the Einstein Cross. A possible fifth image is detected; however, this detection less certain. Fits to the combined spectrum of all images of the Einstein Cross assuming a simple power law with Galactic and intervening absorption at the lensing galaxy yield a photon index of 1.90(+0.05,-0.05). For the first observation, this spectral model yields a 0.4-8.0 keV X-ray flux of 4.6e-13 erg cm-2 s-1 and a 0.4-8.0 keV lensed luminosity of 1.0e46 erg s-1. The source exhibits variability both over long and short time scales. The X-ray flux has dropped by 20% between the two observations, and the K-S test showed that image A is variable at the 97% confidence level within the first observation. Furthermore, a possible time-delay of 2.7(+0.5,-0.9) hours between images A and B with image A leading is detected in the first observation. The X-ray flux ratios of the images are consistent with the optical flux ratios which are affected by microlensing suggesting that the X-ray emission is also microlensed. A comparison between our measured column densities and those inferred from extinction measurements suggests a higher dust-to-gas ratio in the lensing galaxy than the average value of our Galaxy. Finally, we report the detection at the 99.99% confidence level of a broad emission feature near the redshifted energy of the Fe K\alpha line in only the spectrum of image A. The rest frame energy, width, and equivalent width of this feature are E = 5.7(+0.2,-0.3) keV, sigma = 0.87(+0.30,-0.15) keV, and EW = 1200(+300,-200) eV, respectively.

The mini Broad Absorption Line (BAL) quasar RX J0911.4+0551 was observed with the Advanced CCD Imaging Spectrometer (ACIS) of the Chandra X-ray Observatory for ~ 29 ks as part of a gravitational lens (GL) survey aimed at measuring time-delays. Timing analysis of the light-curve of the lensed image A2 shows a rapid flux variation with a duration of about 2000s. A Kolmogorov-Smirnov test shows that the probability that a constant-intensity source would produce the observed variability is less than ~ 0.2 percent. We discuss possible origins for the observed short-term X-ray variability. Our gravitational lens models for the RX J0911.4+0551 GL system predict a time-delay of less than a day between images A1 and A2. The rapid variability combined with the predicted short-time delay make RX J0911.4+0551 an ideal system to apply the GL method for estimating the Hubble constant. We describe the prospects of measuring H_0 within single X-ray observations of GL systems with relatively short time delays. Modeling of the spectrum of the mini-BAL quasar RX J0911.4+0551 suggests the presence of an intrinsic absorber. Partial covering models are slightly preferred over models that contain absorption due to intrinsic ionized or neutral gas.