• Accelerating gravitational microlensing simulations using the Xeon Phi coprocessor(1703.09707)

March 28, 2017 cs.DC, astro-ph.IM, astro-ph.HE
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.
• Algorithms And Programs For Strong Gravitational Lensing In Kerr Space-time Including Polarization(1505.02714)

May 11, 2015 astro-ph.HE
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.
• Evidence for AGN Feedback in the Broad Absorption Lines and Reddening of Mrk 231(1404.6795)

May 3, 2014 astro-ph.GA
We present the first J-band spectrum of Mrk 231, which reveals a large \ion{He}{1}*$\lambda 10830$ broad absorption line with a profile similar to that of the well-known \ion{Na}{1} broad absorption line. Combining this spectrum with optical and UV spectra from the literature, we show that the unusual reddening noted by \citet{veilleux13} is explained by a reddening curve like those previously used to explain low values of total-to-selective extinction in SNe Ia. The nuclear starburst may be the origin and location of the dust. Spatially-resolved emission in the broad absorption line trough suggests nearly full coverage of the continuum emission region. The broad absorption lines reveal higher velocities in the \ion{He}{1}* lines (produced in the quasar-photoionized \ion{H}{2} region) compared with the \ion{Na}{1} and \ion{Ca}{2} lines (produced in the corresponding partially-ionized zone). {\it Cloudy} simulations show that a density increase is required between the \ion{H}{2} and partially-ionized zones to produce ionic column densities consistent with the optical and IR absorption line measurements and limits, and that the absorber lies $\sim 100\rm \, pc$ from the central engine. These results suggest that the \ion{He}{1}* lines are produced in an ordinary quasar BAL wind that impacts upon, compresses, and accelerates the nuclear starburst's dusty effluent (feedback in action), and the \ion{Ca}{2} and \ion{Na}{1} lines are produced in this dusty accelerated gas. This unusual circumstance explains the rarity of \ion{Na}{1} absorption lines; without the compression along our line of sight, Mrk~231 would appear as an ordinary FeLoBAL.
• Effects of Kerr Strong Gravity on Quasar X-ray Microlensing(1211.6487)

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.
• Inclination-Dependent AGN Flux Profiles From Strong Lensing of the Kerr Space-Time(1211.2510)

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 ultraviolet (UV) spectroscopy and photometry of four Type Ia supernovae (SNe 2004dt, 2004ef, 2005M, and 2005cf) obtained with the UV prism of the Advanced Camera for Surveys on the Hubble Space Telescope. This dataset provides unique spectral time series down to 2000 Angstrom. Significant diversity is seen in the near maximum-light spectra (~ 2000--3500 Angstrom) for this small sample. The corresponding photometric data, together with archival data from Swift Ultraviolet/Optical Telescope observations, provide further evidence of increased dispersion in the UV emission with respect to the optical. The peak luminosities measured in uvw1/F250W are found to correlate with the B-band light-curve shape parameter dm15(B), but with much larger scatter relative to the correlation in the broad-band B band (e.g., ~0.4 mag versus ~0.2 mag for those with 0.8 < dm15 < 1.7 mag). SN 2004dt is found as an outlier of this correlation (at > 3 sigma), being brighter than normal SNe Ia such as SN 2005cf by ~0.9 mag and ~2.0 mag in the uvw1/F250W and uvm2/F220W filters, respectively. We show that different progenitor metallicity or line-expansion velocities alone cannot explain such a large discrepancy. Viewing-angle effects, such as due to an asymmetric explosion, may have a significant influence on the flux emitted in the UV region. Detailed modeling is needed to disentangle and quantify the above effects.
• Type II Supernovae as Probes of Cosmology(0902.3142)

Feb. 18, 2009 astro-ph.CO
- Constraining the cosmological parameters and understanding Dark Energy have tremendous implications for the nature of the Universe and its physical laws. - The pervasive limit of systematic uncertainties reached by cosmography based on Cepheids and Type Ia supernovae (SNe Ia) warrants a search for complementary approaches. - Type II SNe have been shown to offer such a path. Their distances can be well constrained by luminosity-based or geometric methods. Competing, complementary, and concerted efforts are underway, to explore and exploit those objects that are extremely well matched to next generation facilities. Spectroscopic follow-up will be enabled by space- based and 20-40 meter class telescopes. - Some systematic uncertainties of Type II SNe, such as reddening by dust and metallicity effects, are bound to be different from those of SNe Ia. Their stellar progenitors are known, promising better leverage on cosmic evolution. In addition, their rate - which closely tracks the ongoing star formation rate - is expected to rise significantly with look- back time, ensuring an adequate supply of distant examples. - These data will competitively constrain the dark energy equation of state, allow the determination of the Hubble constant to 5%, and promote our understanding of the processes involved in the last dramatic phases of massive stellar evolution.
• The New Standard Stellar Population Models (NSSPM) -- The Prologue(astro-ph/0701666)

Jan. 23, 2007 astro-ph
We are developing a brand new stellar population models with flexible chemistry (isochrones plus stellar colors and spectra) in order to set a new standard of completeness and excellence. Here we present preliminary results to assess the effects of stellar evolution models and stellar model atmosphere to the well-known Lick indices at constant heavy element mass fraction Z that self-consistently account for varying heavy element mixtures. We have enhanced chemical elements one by one. Our ultimate goal is to demonstrate 10% absolute mean ages for a sample of local galaxies derived from an integrated light spectrum.
• Joint Efficient Dark-energy Investigation (JEDI): a Candidate Implementation of the NASA-DOE Joint Dark Energy Mission (JDEM)(astro-ph/0507043)

July 14, 2005 astro-ph
JEDI will probe dark energy in three independent ways by measuring the expansion history of the universe: (1) using type Ia supernovae as cosmological standard candles over a range of distances, (2) using baryon oscillations as a cosmological standard ruler over a range of cosmic epochs, (3) mapping the weak gravitational lensing distortion by foreground galaxies of the images of background galaxies at different distances. JEDI will unravel the nature of dark energy with accuracy and precision. JEDI is a 2m-class space telescope with the unique ability of simultaneous wide-field imaging (0.8-4.2 micron in five bands) and multi-slit spectroscopy (0.8-3.2 micron) with a field of view of 1 square degree. What makes JEDI efficient is its ability to simultaneously obtain high signal-to-noise ratio, moderate resolution slit spectra for all supernovae and ~ 5000 galaxies in its wide field of view, and to combine imaging and spectroscopy so that the appropriate balance of time is devoted to each. JEDI will measure the cosmic expansion history H(z) as a free function to < 2% accuracy in redshift bins of 0.2-0.3. Assuming a flat universe and \sigma(\Omega_m)=0.01 (0.03), JEDI could measure a departure from a vanilla \Lambda CDM model (w_0=-1, w'=0) with \sigma(w_0)=0.013 (0.031) and \sigma(w')=0.046 (0.063). JEDI will obtain the well-sampled lightcurves in Z, J, H, K, L bands and spectra of ~ 14,000 type Ia supernovae with redshifts ranging from 0 to 1.7; the redshifts of ~ 10-100 million galaxies to H ~ 23 and z ~ 4 over 1000-10,000 square degrees; and measurements of the shapes of galaxies over 1000-10,000 square degrees in Z,J,H,K,L for > 10^9 galaxies to H \~ 25.
• Non-LTE model atmosphere analysis of the early ultraviolet spectra of nova OS Andromedae 1986(astro-ph/9608199)

Aug. 29, 1996 astro-ph
We have analyzed the early optically thick ultraviolet spectra of Nova OS And 1986 using a grid of spherically symmetric, non-LTE, line-blanketed, expanding model atmospheres and synthetic spectra with the following set of parameters: $5,000\le$ T$_{model}$ $\le 60,000$K, solar abundances, $\rho \propto r^{-3}$, $\v_{max} = 2000\kms$, $L=6 \times 10^{4}\Lsun$, and a statistical or microturbulent velocity of 50 $\kms$. We used the synthetic spectra to estimate the model parameters corresponding to the observed {\it IUE} spectra. The fits to the observations were then iteratively improved by changing the parameters of the model atmospheres, in particular T$_{model}$ and the abundances, to arrive at the best fits to the optically thick pseudo-continuum and the features found in the {\it IUE} spectra. The {\it IUE} spectra show two different optically thick subphases. The earliest spectra, taken a few days after maximum optical light, show a pseudo-continuum created by overlapping absorption lines. The later observations, taken approximately 3 weeks after maximum light, show the simultaneous presence of allowed, semi-forbidden, and forbidden lines in the observed spectra. Analysis of these phases indicate that OS And 86 had solar metallicities except for Mg which showed evidence of being underabundant by as much as a factor of 10. We determine a distance of 5.1 kpc to OS And 86 and derive a peak bolometric luminosity of $\sim$ 5 $\times$ 10$^4$ L$_{\odot}$. The computed nova parameters provide insights into the physics of the early outburst and explain the spectra seen by {\it IUE}. Lastly, we find evidence in the later observations for large non-LTE effects of Fe{\sc ii} which, when included, lead to much better agreement with the observations.