-
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.
-
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 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.
-
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 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.
-
- 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.
-
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.
-
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.
-
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.