We present optical and near-infrared photometry and spectroscopy of the Type IIn supernova (SN) 2014ab, obtained by the Carnegie Supernova Project II (CSP-II) and initiated immediately after its optical discovery. We also present mid-infrared photometry obtained by the Wide-field Infrared Survey Explorer (WISE) satellite extending from 56 days prior to the optical discovery to over 1600 days. The light curve of SN 2014ab evolves slowly, while the spectra exhibit strong emission features produced from the interaction between rapidly expanding ejecta and dense circumstellar matter. The light curve and spectral properties are very similar to those of SN 2010jl. The estimated mass-loss rate of the progenitor of SN 2014ab is of the order of 0.1 Msun/yr under the assumption of spherically symmetric circumstellar matter and steady mass loss. Although the mid-infrared luminosity increases due to emission from dust, which is characterized by a blackbody temperature close to the dust evaporation temperature (~ 2000 K), no clear signatures of in situ dust formation within the cold dense shell located behind the forward shock are observed in SN 2014ab in early phases. Mid-infrared emission of SN 2014ab may originate from pre-existing dust located within dense circumstellar matter that is heated by the SN shock or shock-driven radiation. Finally, for the benefit of the community, we also present in an Appendix five near-infrared spectra of SN 2010jl obtained between 450 to 1300 days post discovery.

We perform two-dimensional and three-dimensional simulations of cold, dense clouds, which are accelerated by radiation pressure on dust relative to a hot, diffuse background gas. We examine the relative effectiveness of acceleration by ultraviolet and infrared radiation fields, both independently and acting simultaneously on the same cloud. We study clouds that are optically thin to infrared emission but with varying ultraviolet optical depths. Consistent with previous work, we find relatively efficient acceleration and long cloud survival times when the infrared band flux dominates over the ultraviolet flux. However, when ultraviolet is dominant or even a modest percentage ($\sim 5-10$\%) of the infrared irradiating flux, it can act to compress the cloud, first crushing it and then disrupting the outer layers. This drives mixing of outer regions of the dusty gas with the hot diffuse background to the point where most dust is not likely to survive or stay coupled to the gas. Hence, the cold cloud is unable to survive for a long enough timescale to experience significant acceleration before disruption even though efficient infrared cooling keeps the majority of the gas close to radiative equilibrium temperature ($T \lesssim 100$K). We discuss implications for observed systems, concluding that radiation pressure driving is most effective when the light from star-forming regions is efficiently reprocessed into the infrared.

We present a new method to photometrically delineate between various sub-types of type Ia supernovae (SNe Ia). Using the color-stretch parameters, $s_{BV}$ or $s_{gr}$, and the time of i-band primary maximum relative to the B-band or g-band maximum it is demonstrated that 2003fg-like, 1991bg-like, and 2002cx-like SNe Ia can readily be identified. In the cases of these extreme SNe Ia, their primary i-band maximum occurs after the time of the B or g band maxima. We suggest that the timing of the i-band maximum can reveal the physical state of the SN Ia explosion as it traces: i) the speed of the recombination front of iron group elements in the ejecta, ii) the temperature evolution and rate of adiabatic cooling in the ejecta and, iii) the presence of interaction with a stellar envelope. This photometric sub-typing can be used in conjunction with other SNe analysis, such as the Branch diagram, to examine the physics and diversity of SNe Ia. The results here can also be used to screen out non-Ia SNe from cosmological samples that do not have complete spectroscopic typing. Finally, as future surveys like LSST create large databases of light curves of many objects this photometric identification can be used to readily identify and study the rates and bulk properties of peculiar SNe Ia.

We describe the structure and implementation of a radiation hydrodynamic solver for MANGA, the moving-mesh hydrodynamics module of the large-scale parallel code, Charm N-body GrAvity solver (ChaNGa). We solve the equations of time dependent radiative transfer using a reduced speed of light approximation following the algorithm of Jiang et al (2014). By writing the radiative transfer equations as a generalized conservation equation, we solve the transport part of these equations on an unstructured Voronoi mesh. We then solve the source part of the radiative transfer equations following Jiang et al (2014) using an implicit solver, and couple this to the hydrodynamic equations. The use of an implicit solver ensure reliable convergence and preserves the conservation properties of these equations even in situations where the source terms are stiff due to the small coupling timescales between radiation and matter. We present the results of a limited number of test cases (energy conservation, momentum conservation, dynamic diffusion, linear waves, crossing beams, and multiple shadows) to show convergence with analytic results and numerical stability. We also show that it produces qualitatively the correct results in the presence of multiple sources in the optically thin case.

It has now become possible to study directly, via numerical simulation, the evolution of relativistic, radiation-dominated flows around compact objects. With this in mind we set out explicitly covariant forms of the radiative transfer equation that are suitable for numerical integration in curved spacetime or flat spacetime in curvilinear coordinates. Our work builds on and summarizes in consistent form earlier work by Lindquist, Thorne, Morita and Kaneko, and others. We give explicitly the basic equations in spherical-polar coordinates for Minkowski space and the Kerr spacetime in Kerr-Schild coordinates.

We present $81$ near-infrared (NIR) spectra of $30$ Type II supernovae (SNe II) from the Carnegie Supernova Project-II (CSP-II), the largest such dataset published to date. We identify a number of NIR features and characterize their evolution over time. The NIR spectroscopic properties of SNe II fall into two distinct groups. This classification is first based on the strength of the He I $\lambda1.083\,\mu$m absorption during the plateau phase; SNe II are either significantly above (spectroscopically strong) or below $50$ angstroms (spectroscopically weak) in pseudo equivalent width. However between the two groups, other properties, such as the timing of CO formation and the presence of Sr II, are also observed. Most surprisingly, the distinct weak and strong NIR spectroscopic classes correspond to SNe II with slow and fast declining light curves, respectively. These two photometric groups match the modern nomenclature of SNe IIP and IIL. Including NIR spectra previously published, 18 out of 19 SNe II follow this slow declining-spectroscopically weak and fast declining-spectroscopically strong correspondence. This is in apparent contradiction to the recent findings in the optical that slow and fast decliners show a continuous distribution of properties. The weak SNe II show a high-velocity component of helium that may be caused by a thermal excitation from a reverse-shock created by the outer ejecta interacting with the red supergiant wind, but the origin of the observed dichotomy is not understood. Further studies are crucial in determining whether the apparent differences in the NIR are due to distinct physical processes or a gap in the current data set.

Type II supernovae (SNe II) show strong hydrogen features in their spectra throughout their whole evolution while type IIb supernovae (SNe IIb) spectra evolve from dominant hydrogen lines at early times to increasingly strong helium features later on. However, it is currently unclear whether the progenitors of these supernova (SN) types form a continuum in pre-SN hydrogen mass or whether they are physically distinct. SN light-curve morphology directly relates to progenitor and explosion properties such as the amount of hydrogen in the envelope, the pre-SN radius, the explosion energy and the synthesized mass of radioactive material. In this work we study the morphology of the optical-wavelength light curves of hydrogen-rich SNe II and hydrogen-poor SNe IIb to test whether an observational continuum exists between the two. Using a sample of 95 SNe (73 SNe II and 22 SNe IIb), we define a range of key observational parameters and present a comparative analysis between both types. We find a lack of events that bridge the observed properties of SNe II and SNe IIb. Light curve parameters such as rise times and post-maximum decline rates and curvatures clearly separate both SN types and we therefore conclude that there is no continuum, with the two SN types forming two observationally distinct families. In the V-band a rise time of 17 days (SNe II lower, SNe IIb higher), and a magnitude difference between 30 and 40 days post explosion of 0.4 mag (SNe II lower, SNe IIb higher) serve as approximate thresholds to differentiate both types.

We present high-cadence ultraviolet (UV), optical, and near-infrared (NIR) data on the luminous Type II-P supernova SN 2017gmr from hours after discovery through the first 180 days. SN 2017gmr does not show signs of narrow, high-ionization emission lines in the early optical spectra, yet the optical lightcurve evolution suggests that an extra energy source from circumstellar medium (CSM) interaction must be present for at least 2 days after explosion. Modeling of the early lightcurve indicates a ~500R$_{\odot}$ progenitor radius, consistent with a rather compact red supergiant, and late-time luminosities indicate up to 0.130 $\pm$ 0.026 M$_{\odot}$ of $^{56}$Ni are present, if the lightcurve is solely powered by radioactive decay, although the $^{56}$Ni mass may be lower if CSM interaction contributes to the post-plateau luminosity. Prominent multi-peaked emission lines of H$\alpha$ and [O I] emerge after day 154, as a result of either an asymmetric explosion or asymmetries in the CSM. The lack of narrow lines within the first two days of explosion in the likely presence of CSM interaction may be an example of close, dense, asymmetric CSM that is quickly enveloped by the spherical supernova ejecta.

Our recent work demonstrates a correlation between the high-velocity blue edge, $v_{edge}$, of the iron-peak Fe/Co/Ni $H$-band emission feature and the optical light curve shape of normal, transitional and sub-luminous type Ia Supernovae (SNe Ia). We explain this correlation in terms of SN Ia physics. $v_{edge}$ corresponds to the sharp transition between the complete and incomplete silicon burning regions in the ejecta. It measures the point in velocity space where the outer $^{56}$Ni mass fraction, $X_{\rm{Ni}}$, falls to the order of 0.03-0.10. For a given $^{56}$Ni mass, $M(^{56}Ni)$, $v_{edge}$ is sensitive to the specific kinetic energy $E_{\rm kin}$($M(^{56}Ni)/M_{WD}$) of the corresponding region. Combining $v_{edge}$ with light curve parameters (i.e., s$_{BV}$, $\Delta m_{15,s}$ in $B$ and $V$) allows us to distinguish between explosion scenarios. The correlation between $v_{edge}$ and light-curve shape is consistent with explosion models near the Chandrasekhar limit. However, the available sub-$M_{Ch}$ WD explosion model based on SN 1999by exhibits velocities which are too large to explain the observations. Finally, the sub-luminous SN 2015bo exhibits signatures of a dynamical merger of two WDs demonstrating diversity among explosion scenarios at the faint end of the SNe Ia population.

We report the discovery of the diatomic molecule SiO in the gas phase in the environment of the recurrent nova T Coronae Borealis. While some of the SiO is photospheric, a substantial portion must arise in the wind from the red giant component of T CrB. A simple fit to the SiO feature, assuming local thermodynamic equilibrium, suggests a SiO column density of 2.8x10^17 /cm2 and temperature ~1000K; the SiO column density is similar to that present in the winds of field red giants. A search for SiO maser emission is encouraged both before and after the next anticipated eruption. We find that the 12C/13C ratio in the red giant is <9, with a best fit value of ~5, a factor ~18 times lower than the solar value of 89. We find no convincing evidence for the presence of dust in the environment of T CrB, which we attribute to the destructive effects on nucleation sites of hard X-ray emission. When the next eruption of T CrB occurs, the ejected material will shock the wind, producing X-ray and coronal line emission, as is the case for the recurrent nova RS Oph. T CrB is also a good candidate for very high energy gamma-ray emission, as first observed during the 2010 outburst of V407 Cyg. We include in the paper a wide variety of infrared spectroscopic and photometric data.

We present the $H$-band wavelength region of thirty post-maximum light near-infrared (NIR) spectra of fourteen transitional and sub-luminous type Ia supernovae (SNe Ia), extending from $+$5d to +20d relative to the epoch of $B$-band maximum. We introduce a new observable, the blue-edge velocity, $v_{edge}$, of the prominent Fe/Co/Ni-peak $H$-band emission feature which is quantitatively measured. The $v_{edge}$ parameter is found to slowly decrease over sub-type ranging from around $-$13,000km/s for transitional SNe~Ia, down to $-$5,000km/s for the sub-luminous SNe Ia. Furthermore, inspection of the +10$\pm$3d spectra indicates that $v_{edge}$ is correlated with the color-stretch parameter, s$_{BV}$, and hence with peak luminosity. These results follow the previous findings that brighter SNe Ia tend to have $^{56}$Ni located at higher velocities as compared to sub-luminous objects. As $v_{edge}$ is a model-independent parameter, we propose it can be used in combination with traditional observational diagnostics to provide a new avenue to robustly distinguish between leading SNe Ia explosion models.

We use global three dimensional radiation magneto-hydrodynamic simulations to study the properties of inner regions of accretion disks around a 5\times 10^8 solar mass black hole with mass accretion rates reaching 7% and 20% of the Eddington value. This region of the disk is supported by magnetic pressure with surface density significantly smaller than the values predicted by the standard thin disk model but with a much larger disk scale height. The disks do not show any sign of thermal instability over many thermal time scales. More than half of the accretion is driven by radiation viscosity in the optically thin corona region for the lower accretion rate case, while accretion in the optically thick part of the disk is driven by the Maxwell and Reynolds stresses from MRI turbulence. Coronae with gas temperatures > 10^8 K are generated only in the inner \approx 10 gravitational radii in both simulations, being more compact in the higher accretion rate case. In contrast to the thin disk model, surface density increases with increasing mass accretion rate, which causes less dissipation in the optically thin region and a relatively weaker corona. The simulation results may explain the formation of X-ray coronae in Active Galactic Nuclei (AGNs), the compact size of such coronae, and the observed trend of optical to X-ray luminosity with Eddington ratio for many AGNs.

Shifting the focus of Type Ia supernova (SN Ia) cosmology to the near-infrared (NIR) is a promising way to significantly reduce the systematic errors, as the strategy minimizes our reliance on the empirical width-luminosity relation and uncertain dust laws. Observations in the NIR are also crucial for our understanding of the origins and evolution of these events, further improving their cosmological utility. Any future experiments in the rest-frame NIR will require knowledge of the SN Ia NIR spectroscopic diversity, which is currently based on a small sample of observed spectra. Along with the accompanying paper, Phillips et al. (2018), we introduce the Carnegie Supernova Project-II (CSP-II), to follow up nearby SNe Ia in both the optical and the NIR. In particular, this paper focuses on the CSP-II NIR spectroscopy program, describing the survey strategy, instrumental setups, data reduction, sample characteristics, and future analyses on the data set. In collaboration with the Harvard-Smithsonian Center for Astrophysics (CfA) Supernova Group, we obtained 661 NIR spectra of 157 SNe Ia. Within this sample, 451 NIR spectra of 90 SNe Ia have corresponding CSP-II follow-up light curves. Such a sample will allow detailed studies of the NIR spectroscopic properties of SNe Ia, providing a different perspective on the properties of the unburned material, radioactive and stable nickel produced, progenitor magnetic fields, and searches for possible signatures of companion stars.

By fitting synthetic spectral models computed via the TLUSTY code, we examine how the spectra from thin accretion disks are expected to vary in accreting black hole systems. We fit color-corrected blackbody models to our synthetic spectra to estimate the spectral hardening factor f, which parameterizes the departure from blackbody and is commonly used to help interpret multitemperature blackbody fitting results. We find we can define a reasonably robust f value to spectra when the effects of Compton scattering dominate radiation transfer. We examine the evolution of f with black hole mass and accretion rate, typically finding a moderate variation (f ~ 1.4-2) for accretion rates between 1% and 100% of the Eddington rate. Consistent with most previous work, we find f tends to increase with accretion rate, but we also infer a weaker correlation of f with black holes mass. We find that f is rarely much larger than 2 unless the disk becomes photon starved, in contention with some previous calculations. Significant spectral hardening (f > 2) is only found when the disk mass surface density is lower than expected for alpha-disk models unless alpha is near unity or larger.

Context: It has been predicted that the object KIC 9832227 - a contact binary star - will undergo a merger in 2022.2+/-0.7. We describe the near infrared spectrum of this object as an impetus to obtain pre-merger data. Aims: We aim to characterise (i)~the nature of the individual components of the binary and (ii)~the likely circumbinary environment, so that the merger - if and when it occurs - can be interpreted in an informed manner. Methods: We use infrared spectroscopy in the wavelength range 0.7micron - 2.5microns, to which we fit model atmospheres to represent the individual stars. We use the binary ephemeris to determine the orbital phase at the time of observation. Results: We find that the infrared spectrum is best fitted by a single component having effective temperature 5920K, log[g]=4.1 and solar metallicity, consistent with the fact that the system was observed at conjunction. Conclusions: The strength of the infrared H lines is consistent with a high value of log[g], and the strength of the CaII triplet indicates the presence of a chromosphere, as might be expected from rapid stellar rotation. The HeI absorption we observe likely arises in He excited by coronal activity in a circumstellar envelope, suggesting that the weakness of the CaII triplet is also likely chromospheric in origin.

We perform two-dimensional and three-dimensional radiation hydrodynamic simulations to study cold clouds accelerated by radiation pressure on dust in the environment of rapidly star-forming galaxies dominated by infrared flux. We utilize the reduced speed of light approximation to solve the frequency-averaged, time-dependent radiative transfer equation. We find that radiation pressure is capable of accelerating the clouds to hundreds of kilometers per second while remaining dense and cold, consistent with observations. We compare these results to simulations where acceleration is provided by entrainment in a hot wind, where the momentum injection of the hot flow is comparable to the momentum in the radiation field. We find that the survival time of the cloud accelerated by the radiation field is significantly longer than that of a cloud entrained in a hot outflow. We show that the dynamics of the irradiated cloud depends on the initial optical depth, temperature of the cloud, and the intensity of the flux. Additionally, gas pressure from the background may limit cloud acceleration if the density ratio between the cloud and background is $\lesssim 10^{2}$. In general, a 10 pc-scale optically thin cloud forms a pancake structure elongated perpendicular to the direction of motion, while optically thick clouds form a filamentary structure elongated parallel to the direction of motion. The details of accelerated cloud morphology and geometry can also be affected by other factors, such as the cloud lengthscale, the reduced speed of light approximation, spatial resolution, initial cloud structure, and the dimensionality of the run, but these have relatively little affect on the cloud velocity or survival time.

Hydrogen-poor superluminous supernovae (SLSNe-I) have been predominantly found in low-metallicity, star-forming dwarf galaxies. Here we identify Gaia17biu/SN 2017egm as an SLSN-I occurring in a "normal" spiral galaxy (NGC 3191) in terms of stellar mass (several times 10^10 M_sun) and metallicity (roughly Solar). At redshift z=0.031, Gaia17biu is also the lowest redshift SLSN-I to date, and the absence of a larger population of SLSNe-I in dwarf galaxies of similar redshift suggests that metallicity is likely less important to the production of SLSNe-I than previously believed. With the smallest distance and highest apparent brightness for an SLSN-I, we are able to study Gaia17biu in unprecedented detail. Its pre-peak near-ultraviolet to optical color is similar to that of Gaia16apd and among the bluest observed for an SLSN-I while its peak luminosity (M_g = -21 mag) is substantially lower than Gaia16apd. Thanks to the high signal-to-noise ratios of our spectra, we identify several new spectroscopic features that may help to probe the properties of these enigmatic explosions. We detect polarization at the ~0.5% level that is not strongly dependent on wavelength, suggesting a modest, global departure from spherical symmetry. In addition, we put the tightest upper limit yet on the radio luminosity of an SLSN-I with <5.4x10^26 erg/s/Hz (at 10 GHz), which is almost a factor of 40 better than previous upper limits and one of the few measured at an early stage in the evolution of an SLSN-I. This limit largely rules out an association of this SLSNe-I with known populations of gamma-ray burst (GRB) like central engines.

We use global three dimensional radiation magneto-hydrodynamical simulations to study accretion disks onto a $5\times 10^8M_{\odot}$ black hole with accretion rates varying from $\sim 250L_{Edd}/c^2$ to $1500 L_{Edd}/c^2$. We form the disks with torus centered at $50-80$ gravitational radii with self-consistent turbulence initially generated by the magneto-rotational instability. We study cases with and without net vertical magnetic flux. The inner regions of all disks have radiation pressure $\sim 10^4-10^6$ times the gas pressure. Non-axisymmetric density waves that steepen into spiral shocks form as gas flows towards the black hole. In simulations without net vertical magnetic flux, Reynolds stress generated by the spiral shocks are the dominant mechanism to transfer angular momentum. Maxwell stress from MRI turbulence can be larger than the Reynolds stress only when net vertical magnetic flux is sufficiently large. Outflows are formed with speed $\sim 0.1-0.4c$. When the accretion rate is smaller than $\sim 500 L_{Edd}/c^2$, outflows start around $10$ gravitational radii and the radiative efficiency is $\sim 5\%-7\%$ with both magnetic field configurations. With accretion rate reaching $1500 L_{Edd}/c^2$, most of the funnel region close to the rotation axis becomes optically thick and the outflow only develops beyond $50$ gravitational radii. The radiative efficiency is reduced to $1\%$. We always find the kinetic energy luminosity associated with the outflow is only $\sim 15\%-30\%$ of the radiative luminosity. The mass flux lost in the outflow is $\sim 15\%-50\%$ of the net mass accretion rates. We discuss implications of our simulation results on the observational properties of these disks.

We present very early, high-cadence photometric observations of the nearby Type Ia SN 2017cbv. The light curve is unique in that it has a blue bump during the first five days of observations in the U, B, and g bands, which is clearly resolved given our photometric cadence of 5.7 hr during that time span. We model the light curve as the combination of early shocking of the supernova ejecta against a nondegenerate companion star plus a standard SN Ia component. Our best-fit model suggests the presence of a subgiant star 56 solar radii from the exploding white dwarf, although this number is highly model-dependent. While this model matches the optical light curve well, it overpredicts the observed flux in the ultraviolet bands. This may indicate that the shock is not a blackbody, perhaps because of line blanketing in the UV. Alternatively, it could point to another physical explanation for the optical blue bump, such as interaction with circumstellar material or an unusual nickel distribution. Early optical spectra of SN 2017cbv show strong carbon absorption up through day -13 with respect to maximum light, suggesting that the progenitor system contains a significant amount of unburned material. These early results on SN 2017cbv illustrate the power of early discovery and intense follow-up of nearby supernovae to resolve standing questions about the progenitor systems and explosion mechanisms of SNe Ia.

We study dusty winds driven by radiation pressure in the atmosphere of a rapidly star-forming environment. We apply the variable Eddington tensor algorithm to re-examine the two-dimensional radiation hydrodynamic problem of a column of gas that is accelerated by a constant infrared radiation flux. In the absence of gravity, the system is primarily characterized by the initial optical depth of the gas. We perform several runs with different initial optical depth and resolution. We find that the gas spreads out along the vertical direction, as its mean velocity and velocity dispersion increase. In contrast to previous work using flux-limited diffusion algorithm, we find little evolution in the trapping factor. The momentum coupling between radiation and gas in the absence of gravity is similar to that with gravity. For Eddington ratio increasing with the height in the system, the momentum transfer from the radiation to the gas is not merely $\sim L/c$, but amplified by a factor of $1+\eta \tau_{\rm IR}$, where $\tau_{\rm IR}$ is the integrated infrared optical depth through the system, and $\eta\sim0.5-0.9$, decreasing with the optical depth. We apply our results to the atmosphere of galaxies and conclude that radiation pressure may be an important mechanism for driving winds in the most rapidly star-forming galaxies and starbursts.

M82 X-1 is one of the brightest ultraluminous X-ray sources (ULXs) known, which, assuming Eddington-limited accretion and other considerations, makes it one of the best intermediate-mass black hole (IMBH) candidates. However, the ULX may still be explained by super-Eddington accretion onto a stellar-remnant black hole. We present simultaneous NuSTAR, Chandra and Swift/XRT observations during the peak of a flaring episode with the aim of modeling the emission of M82 X-1 and yielding insights into its nature. We find that thin-accretion disk models all require accretion rates at or above the Eddington limit in order to reproduce the spectral shape, given a range of black hole masses and spins. Since at these high Eddington ratios the thin-disk model breaks down due to radial advection in the disk, we discard the results of the thin-disk models as unphysical. We find that the temperature profile as a function of disk radius ($T(r)\propto r^{-p}$) is significantly flatter ($p=0.55^{+ 0.07}_{- 0.04}$) than expected for a standard thin disk ($p=0.75$). A flatter profile is instead characteristic of a slim disk which is highly suggestive of super-Eddington accretion. Furthermore, radiation hydrodynamical simulations of super-Eddington accretion have shown that the predicted spectra of these systems are very similar to what we observe for M82 X-1. We therefore conclude that M82 X-1 is a super-Eddington accretor. Our mass estimates inferred from the inner disk radius imply a stellar-remnant black hole ($M_{\rm BH}=26^{+9}_{-6} M_{\odot}$) when assuming zero spin, or an IMBH ($M_{\rm BH}=125^{+45}_{-30} M_{\odot}$) when assuming maximal spin.

Accretion disks around supermassive black holes have regions where the Rosseland mean opacity can be much larger than the electron scattering opacity primarily due to the large number of bound-bound transitions in iron. We study the effects of this iron opacity "bump" on the thermal stability and vertical structure of radiation pressure dominated accretion disks, utilizing three dimensional radiation magneto-hydrodynamic simulations in the local shearing box approximation. The simulations self-consistently calculate the heating due to MHD turbulence caused by magneto-rotational instability and radiative cooling by using the radiative transfer module based on a variable Eddington tensor in Athena. For a $5\times 10^8$ solar mass black hole with $\sim 3\%$ of the Eddington luminosity, a model including the iron opacity bump maintains its structure for more than $10$ thermal times without showing significant signs of thermal runaway. In contrast, if only electron scattering and free-free opacity are included as in the standard thin disk model, the disk collapses on the thermal time scale. The difference is caused by a combination of 1) an anti-correlation between the total optical depth and the temperature, and 2) enhanced vertical advective energy transport. These results suggest that the iron opacity bump may have a strong impact on the stability and structure of AGN accretion disks, and may contribute to a dependence of AGN properties on metallicity. Since this opacity is relevant primarily in UV emitting regions of the flow, it may help to explain discrepancies between observation and theory that are unique to AGNs.

This paper describes the Polarization Spectroscopic Telescope Array (PolSTAR), a mission proposed to NASA's 2014 Small Explorer (SMEX) announcement of opportunity. PolSTAR measures the linear polarization of 3-50 keV (requirement; goal: 2.5-70 keV) X-rays probing the behavior of matter, radiation and the very fabric of spacetime under the extreme conditions close to the event horizons of black holes, as well as in and around magnetars and neutron stars. The PolSTAR design is based on the technology developed for the Nuclear Spectroscopic Telescope Array (NuSTAR) mission launched in June 2012. In particular, it uses the same X-ray optics, extendable telescope boom, optical bench, and CdZnTe detectors as NuSTAR. The mission has the sensitivity to measure ~1% linear polarization fractions for X-ray sources with fluxes down to ~5 mCrab. This paper describes the PolSTAR design as well as the science drivers and the potential science return.

We study super-Eddington accretion flows onto black holes using a global three dimensional radiation magneto-hydrodynamical simulation. We solve the time dependent radiative transfer equation for the specific intensities to accurately calculate the angular distribution of the emitted radiation. Turbulence generated by the magneto-rotational instability provides self-consistent angular momentum transfer. The simulation reaches inflow equilibrium with an accretion rate ~220L_edd/c^2 and forms a radiation driven outflow along the rotation axis. The mechanical energy flux carried by the outflow is ~20% of the radiative energy flux. The total mass flux lost in the outflow is about 29% of the net accretion rate. The radiative luminosity of this flow is ~10L_edd. This yields a radiative efficiency ~4.5%, which is comparable to the value in a standard thin disk model. In our simulation, vertical advection of radiation caused by magnetic buoyancy transports energy faster than photon diffusion, allowing a significant fraction of the photons to escape from the surface of the disk before being advected into the black hole. We contrast our results with the lower radiative efficiencies inferred in most models, such as the slim disk model, which neglect vertical advection. Our inferred radiative efficiencies also exceed published results from previous global numerical simulations, which did not attribute a significant role to vertical advection. We briefly discuss the implications for the growth of supermassive black holes in the early universe and describe how these results provided a basis for explaining the spectrum and population statistics of ultraluminous X-ray sources.

(Abridged) We describe a new algorithm for solving the coupled frequency-integrated transfer equation and the equations of magnetohydrodynamics when the light-crossing time is only marginally shorter than dynamical timescales. The transfer equation is solved in the mixed frame, including velocity dependent source terms accurate to O(v/c). An operator split approach is used to compute the specific intensity along discrete rays, with upwind monotonic interpolation used along each ray to update the transport terms, and implicit methods used to compute the scattering and absorption source terms. Conservative differencing is used for the transport terms, which ensures the specific intensity (as well as energy and momentum) are conserved along each ray to round-off error. The use of implicit methods for the source terms ensures the method is stable even if the source terms are very stiff. To couple the solution of the transfer equation to the MHD algorithms in the Athena code, we perform direct quadrature of the specific intensity over angles to compute the energy and momentum source terms. We present the results of a variety of tests of the method, such as calculating the structure of a non-LTE atmosphere, an advective diffusion test, linear wave convergence tests, and the well-known shadow test. We use new semi-analytic solutions for radiation modified shocks to demonstrate the ability of our algorithm to capture the effects of an anisotropic radiation field accurately. Since the method uses explicit differencing of the spatial operators, it shows excellent weak scaling on parallel computers. The method is ideally suited for problems in which characteristic velocities are non-relativistic, but still within a few percent or more of the speed of light. The method is an intermediate step towards algorithms for fully relativistic flows.

We use our variable Eddington tensor (VET) radiation hydrodynamics code to perform two-dimensional simulations to study the impact of radiation forces on atmospheres composed of dust and gas. Our setup closely follows that of Krumholz & Thompson, assuming that dust and gas are well-coupled and that the radiation field is characterized by blackbodies with temperatures >~ 80 K, as might be found in ultraluminous infrared galaxies. In agreement with previous work, we find that Rayleigh-Taylor instabilities develop in radiation supported atmospheres, leading to inhomogeneities that limit momentum exchange between radiation and dusty gas, and eventually providing a near balance of the radiation and gravitational forces. However, the evolution of the velocity and spatial distributions of the gas differs significantly from previous work, which utilized a less accurate flux-limited diffusion (FLD) method. Our VET simulations show continuous net acceleration of the gas, with no steady-state reached by the end of the simulation. In contrast, FLD results show little net acceleration of the gas and settle in to a quasi-steady, turbulent state with low velocity dispersion. The discrepancies result primarily from the inability of FLD to properly model the variation of the radiation field around structures that are less than a few optical depths across. We conclude that radiation feedback remains a viable mechanism for driving high-Mach number turbulence. We discuss implications for observed systems and global numerical simulations of feedback, but more realistic setups are needed to make robust observational predictions and assess the prospect of launching outflows with radiation.

A new mechanism to form a magnetic pressure supported, high temperature corona above the photosphere of an accretion disk is explored using three dimensional radiation magneto-hydrodynamic (MHD) simulations. The thermal properties of the disk are calculated self-consistently by balancing radiative cooling through the surfaces of the disk with heating due to dissipation of turbulence driven by magneto-rotational instability (MRI). As has been noted in previous work, we find the dissipation rate per unit mass increases dramatically with height above the mid-plane, in stark contrast to the alpha-disk model which assumes this quantity is a constant. Thus, we find that in simulations with a low surface density (and therefore a shallow photosphere), the fraction of energy dissipated above the photosphere is significant (about 3.4% in our lowest surface density model), and this fraction increases as surface density decreases. When a significant fraction of the accretion energy is dissipated in the optically thin photosphere, the gas temperature increases substantially and a high temperature, magnetic pressure supported corona is formed. The volume-averaged temperature in the disk corona is more than 10 times larger than at the disk mid-plane. Moreover, gas temperature in the corona is strongly anti-correlated with gas density, which implies the corona formed by MRI turbulence is patchy. This mechanism to form an accretion disk corona may help explain the observed relation between the spectral index and luminosity from AGNs, and the soft X-ray excess from some AGNs. It may also be relevant to spectral state changes in X-ray binaries.

AGN SEDs generally show a turnover at lambda 1000A, implying a maximal accretion disc (AD) temperature of T_max~50,000K. Massive O stars display a similar T_max, associated with a sharp rise in a line driven mass loss Mdot_wind with increasing surface temperature. AGN AD are also characterized by similar surface gravity to massive O stars. The Mdot_wind of O stars reaches ~10^-5 Msun/year. Since the surface area of AGN AD can be 10^6 larger, the implied Mdot_wind in AGN AD can reach the accretion rate Mdot. A rise to Mdot_wind Mdot towards the AD center may therefore set a similar cap of T_max~50,000K. To explore this idea, we solve the radial structure of an AD with a mass loss term, and calculate the implied AD emission using the mass loss term derived from observations of O stars. We find that Mdot_wind becomes comparable to Mdot typically at a few 10s of GM/c^2. Thus, the standard thin AD solution is effectively truncated well outside the innermost stable orbit. The calculated AD SED shows the observed turnover at lambda~1000A, which is weakly dependent on the AGN luminosity and black hole mass. The AD SED is generally independent of the black hole spin, due to the large truncation radius. However, a cold AD (low Mdot, high black hole mass) is predicted to be windless, and thus its SED should be sensitive to the black hole spin. The accreted gas may form a hot thick disc with a low radiative efficiency inside the truncation radius, or a strong line driven outflow, depending on its ionization state.

As small solid grains grow into larger ones in protoplanetary nebulae, or in the cloudy atmospheres of exoplanets, they generally form porous aggregates rather than solid spheres. A number of previous studies have used highly sophisticated schemes to calculate opacity models for irregular, porous particles with size much smaller than a wavelength. However, mere growth itself can affect the opacity of the medium in far more significant ways than the detailed compositional and/or structural differences between grain constituents once aggregate particle sizes exceed the relevant wavelengths. This physics is not new; our goal here is to provide a model that provides physical insight and is simple to use in the increasing number of protoplanetary nebula evolution, and exoplanet atmosphere, models appearing in recent years, yet quantitatively captures the main radiative properties of mixtures of particles of arbitrary size, porosity, and composition. The model is a simple combination of effective medium theory with small-particle closed-form expressions, combined with suitably chosen transitions to geometric optics behavior. Calculations of wavelength-dependent emission and Rosseland mean opacity are shown and compared with Mie theory. The model's fidelity is very good in all comparisons we have made, except in cases involving pure metal particles or monochromatic opacities for solid particles with size comparable to the wavelength.

We report on the first phase of our study of cloud irradiation. We study irradiation by means of numerical, two-dimensional time-dependent radiation-hydrodynamic simulations of a cloud irradiated by a strong radiation. We adopt a very simple treatment of the opacity, neglect photoionization and gravity, and instead focus on assessing the role of the type and magnitude of the opacity on the cloud evolution. Our main result is that even relatively dense clouds that are radiatively heated (i.e., with significant absorption opacity) do not move as a whole instead they undergo a very rapid and major evolution in its shape, size and physical properties. In particular, the cloud and its remnants become optical thin within less than one sound crossing time and before they can travel over a significant distance (a distance of a few radii of the initial cloud). We also found that a cloud can be accelerated as a whole under quite extreme conditions, e.g., the opacity must be dominated by scattering. However, the acceleration due to the radiation force is relatively small and unless the cloud is optically thin the cloud quickly changes its size and shape. We discuss implications for the modelling and interpetation broad line regions of active galactic nuclei.

We study the long-term thermal stability of radiation dominated disks in which the vertical structure is determined self-consistently by the balance of heating due to dissipation of MHD turbulence driven by the magneto-rotational instability (MRI), and cooling due to radiation emitted at the photosphere. The calculations adopt the local shearing box approximation, and utilize the recently developed radiation transfer module in the Athena MHD code based on a variable Eddington tensor rather than an assumed local closure. After saturation of the MRI, in many cases the disk maintains a steady vertical structure for many thermal times. However, in every case in which the box size in the horizontal directions is at least one pressure scale height, fluctuations associated with MRI turbulence and dynamo action in the disk eventually trigger a thermal runaway which causes the disk to either expand or contract until the calculation must be terminated. During runaway, the dependence of the heating and cooling rates on total pressure satisfy the simplest criterion for classical thermal instability. We identify several physical reasons why the thermal runaway observed in our simulations differ from the standard alpha disk model, for example the advection of radiation contributes a non-negligible fraction to the vertical energy flux at the largest radiation pressure, most of the dissipation does not happen in the disk mid-plane, and the change of dissipation scale height with mid-plane pressure is slower than the change of density scale height. We discuss how and why our results differ from those published previously. Such thermal runaway behavior might have important implications for interpreting temporal variability in observed systems, but fully global simulations are required to study the saturated state before detailed predictions can be made.

Accretion discs are present around both stellar-mass black holes in X-ray binaries and supermassive black holes in active galactic nuclei. A wide variety of circumstantial evidence implies that many of these discs are warped. The standard Bardeen--Petterson model attributes the shape of the warp to the competition between Lense--Thirring torque from the central black hole and viscous angular-momentum transport within the disc. We show that this description is incomplete, and that torques from the companion star (for X-ray binaries) or the self-gravity of the disc (for active galactic nuclei) can play a major role in determining the properties of the warped disc. Including these effects leads to a rich set of new phenomena. For example, (i) when a companion star is present and the warp arises from a misalignment between the companion's orbital axis and the black hole's spin axis, there is no steady-state solution of the Pringle--Ogilvie equations for a thin warped disc when the viscosity falls below a critical value; (ii) in AGN accretion discs, the warp can excite short-wavelength bending waves that propagate inward with growing amplitude until they are damped by the disc viscosity. We show that both phenomena can occur for plausible values of the black hole and disc parameters, and briefly discuss their observational implications.

The saturation level of the magneto-rotational instability (MRI) in a strongly radiation dominated accretion disk is studied using a new Godunov radiation MHD code in the unstratified shearing box approximation. Since vertical gravity is neglected in this work, our focus is on how the MRI saturates in the optically thick mid-plane of the disk. We confirm that turbulence generated by the MRI is very compressible in the radiation dominated regime, as found by previous calculations using the flux-limited diffusion approximation. We also find little difference in the saturation properties in calculations that use a larger horizontal domain (up to four times the vertical scale height in the radial direction). However, in strongly radiation pressure dominated disks (one in which the radiation energy density reaches 1% of the rest mass energy density of the gas), we find Maxwell stress from the MRI turbulence is larger than the value produced when radiation pressure is replaced with the same amount of gas pressure. At the same time, the ratio between Maxwell stress and Reynolds stress is increased by almost a factor of 8 compared with the gas pressure dominated case. We suggest that this effect is caused by radiation drag, which acts like bulk viscosity and changes the effective magnetic Prandtl number of the fluid. Radiation viscosity significantly exceeds both the microscopic plasma viscosity and resistivity, ensuring that radiation dominated systems occupy the high magnetic Prandtl number regime. Nevertheless, we find radiative shear viscosity is negligible compared to the Maxwell and Reynolds stress in the flow. This may have important implications for the structure of radiation dominated accretion disks.

The non-linear regime of Rayleigh-Taylor instability (RTI) in a radiation supported atmosphere, consisting of two uniform fluids with different densities, is studied numerically. We perform simulations using our recently developed numerical algorithm for multi-dimensional radiation hydrodynamics based on a variable Eddington tensor as implemented in Athena, focusing on the regime where scattering opacity greatly exceeds absorption opacity. We find that the radiation field can reduce the growth and mixing rate of RTI, but this reduction is only significant when radiation pressure significantly exceeds gas pressure. Small scale structures are also suppressed in this case. In the non-linear regime, dense fingers sink faster than rarefied bubbles can rise, leading to asymmetric structures about the interface. By comparing the calculations that use a variable Eddington tensor (VET) versus the Eddington approximation, we demonstrate that anisotropy in the radiation field can affect the non-linear development of RTI significantly. We also examine the disruption of a shell of cold gas being accelerated by strong radiation pressure, motivated by models of radiation driven outflows in ultraluminous infrared galaxies. We find that when the growth rate of RTI is smaller than acceleration time scale, the amount of gas that would be pushed away by the radiation field is reduced due to RTI.

It is generally thought that the light coming from the inner plunging region of black hole accretion discs contributes negligibly to the disc's overall spectrum, i.e. the plunging fluid is swallowed by the black hole before it has time to radiate. In the standard disc model used to fit X-ray observations of accretion discs, the plunging region is assumed to be perfectly dark. However, numerical simulations that include the full physics of the magnetized flow predict that a small fraction of the disc's total luminosity emanates from the plunging region. We investigate the observational consequences of this neglected inner light. We compute radiative transfer based disc spectra that correspond to 3D general relativistic magnetohydrodynamic simulated discs (which produce light inside their plunging regions). In the context of black hole spin estimation, we find that the neglected inner light only has a modest effect (this bias is less than typical observational systematic errors). For rapidly spinning black holes, we find that the combined emission from the plunging region produces a weak power-law tail at high energies. This indicates that infalling matter is the origin for some of the `coronal' emission observed in the thermal dominant and steep power-law states of X-ray binaries.

We describe a numerical algorithm to integrate the equations of radiation magnetohydrodynamics in multidimensions using Godunov methods. This algorithm solves the radiation moment equations in the mixed frame, without invoking any diffusion-like approximations. The moment equations are closed using a variable Eddington tensor whose components are calculated from a formal solution of the transfer equation at a large number of angles using the method of short characteristics. We use a comprehensive test suite to verify the algorithm, including convergence tests of radiation-modified linear acoustic and magnetosonic waves, the structure of radiation modified shocks, and two-dimensional tests of photon bubble instability and the ablation of dense clouds by an intense radiation field. These tests cover a very wide range of regimes, including both optically thick and thin flows, and ratios of the radiation to gas pressure of at least 10^{-4} to 10^{4}. Across most of the parameter space, we find the method is accurate. However, the tests also reveal there are regimes where the method needs improvement, for example when both the radiation pressure and absorption opacity are very large. We suggest modifications to the algorithm that will improve accuracy in this case. We discuss the advantages of this method over those based on flux-limited diffusion. In particular, we find the method is not only substantially more accurate, but often no more expensive than the diffusion approximation for our intended applications.

We describe the implementation of a module for the Athena magnetohydrodynamics (MHD) code which solves the time-independent, multi-frequency radiative transfer (RT) equation on multidimensional Cartesian simulation domains, including scattering and non-LTE effects. The module is based on well-known and well-tested algorithms developed for modeling stellar atmospheres, including the method of short characteristics to solve the RT equation, accelerated Lambda iteration to handle scattering and non-LTE effects, and parallelization via domain decomposition. The module serves several purposes: it can be used to generate spectra and images, to compute a variable Eddington tensor (VET) for full radiation MHD simulations, and to calculate the heating and cooling source terms in the MHD equations in flows where radiation pressure is small compared with gas pressure. For the latter case, the module is combined with the standard MHD integrators using operator-splitting and we describe this approach in detail. Implementation of the VET method for radiation pressure dominated flows is described in a companion paper. We present results from a suite of test problems for both the RT solver itself, and for dynamical problems that include radiative heating and cooling. These tests demonstrate that the radiative transfer solution is accurate, and confirm that the operator split method is stable, convergent, and efficient for problems of interest. We demonstrate there is no need to adopt ad-hoc assumptions of questionable accuracy to solve RT problems in concert with MHD: the computational cost for our general-purpose module for simple (e.g. LTE grey) problems can be comparable to or less than a single timestep of Athena's MHD integrators, and only few times more expensive than that for more general problems. (Abridged)

A recent study of the accretion efficiency of the PG sample AGN, based on the thin accretion disk emission model, finds the accretion radiative efficiency is correlated with the black hole mass (eta Mbh^0.5). A followup study suggests the correlation is an artifact induced by selection effects. Here we point out there are two independent effects. The first is a sample selection effect, which leads to a high L/L_Edd AGN sample. The second effect is the observed small spread in the SED shape, which is not induced by selection effects. The second effect is what leads to the eta Mbh^0.5 relation in the PG sample. The physical reason for the small spread in the SED shape is an intriguing open question.

The compact primary in the X-ray binary Cygnus X-1 was the first black hole to be established via dynamical observations. We have recently determined accurate values for its mass and distance, and for the orbital inclination angle of the binary. Building on these results, which are based on our favored (asynchronous) dynamical model, we have measured the radius of the inner edge of the black hole's accretion disk by fitting its thermal continuum spectrum to a fully relativistic model of a thin accretion disk. Assuming that the spin axis of the black hole is aligned with the orbital angular momentum vector, we have determined that Cygnus X-1 contains a near-extreme Kerr black hole with a spin parameter a/M>0.95 (3\sigma). For a less probable (synchronous) dynamical model, we find a/M>0.92 (3\sigma). In our analysis, we include the uncertainties in black hole mass, orbital inclination angle and distance, and we also include the uncertainty in the calibration of the absolute flux via the Crab. These four sources of uncertainty totally dominate the error budget. The uncertainties introduced by the thin-disk model we employ are particularly small in this case given the extreme spin of the black hole and the disk's low luminosity.

(Abridged) Narrow Line Seyfert 1 (NLS1) galaxies have low mass black holes and mass accretion rates close to (or exceeding) Eddington, so a standard blackbody accretion disc should peak in the EUV. However, the lack of true absorption opacity in the disc means that the emission is better approximated by a colour temperature corrected blackbody, and this colour temperature correction is large enough ($\sim 2.4$) that the bare disc emission from a zero spin black hole can extend into the soft X-ray bandpass. Part of the soft X-ray excess seen in these objects must be intrinsic emission from the disc unless the vertical structure is very different to that predicted. However, the soft excess is much broader than predicted by a bare disc spectrum, indicating some Compton upscattering by cool, optically thick material. We associate this with the disc itself, so it must ultimately be powered by mass accretion. We build an energetically self consistent model assuming that the emission thermalises at large radii, but that at smaller radii the gravitational energy is split between powering optically thick Comptonised disc emission (forming the soft X-ray excess) and an optically thin corona above the disc (forming the tail to higher energies). We show examples of this model fit to the extreme NLS1 REJ1034+396, and to the much lower Eddington fraction Broad Line Seyfert 1 PG1048+231. We use these to guide our fits and interpretations of three template spectra made from co-adding multiple sources to track out a sequence of AGN spectra as a function of $L/L_{Edd}$. The new model is publically available within the {\sc xspec} spectral fitting package.

The optical-UV continuum of quasars is broadly consistent with the emission from a geometrically thin optically thick accretion disk (AD). The AD produces the ionizing continuum which powers the broad and narrow emission lines. The maximum AD effective temperature is given by Teff=fmax(Mdot/M^2)^1/4, where M is the black hole mass, Mdot the accretion rate, and fmax is set by the black hole spin a_*. For a low enough value of Mdot/M^2 the AD may become too cold to produce ionizing photons. Such an object will form a lineless quasar. This occurs for a local blackbody (BB) AD with a luminosity Lopt=10^46 erg/s for M>3.6E9 Msun, when a_*=0, and for M>1.4E10 Msun, when a_*=0.998. Using the AD based Mdot, derived from M and Lopt, and the reverberation based M, derived from Lopt and the Hbeta FWHM, v, gives Teff \propto Lopt^-0.13v^-1.45. Thus, Teff is mostly set by v. Quasars with a local BB AD become lineless for v> 8,000 km/s, when a_*=0, and for v> 16,000 km/s, when a_*=0.998. Higher values of v are required if the AD is hotter than a local BB. The AD becoming non-ionizing may explain why line emitting quasars with v>10,000 km/s are rare. Weak low ionization lines may still be present if the X-ray continuum is luminous enough, and such objects may form a population of weak emission line quasars (WLQ). If correct, such WLQ should show a steeply falling SED at lambda<1000A. Such an SED was observed by Hryniewicz et al. in SDSS J094533.99+100950.1, a WLQ observed down to 570A, which is well modeled by a rather cold AD SED. UV spectroscopy of z~1-2 quasars is required to eliminate potential intervening Lyman limit absorption by the intergalactic medium (IGM), and to explore if the SEDs of lineless quasars and some additional WLQ are also well fit by a cold AD SED.

With an inferred bolometric luminosity exceeding 10^42 erg/s, HLX-1 in ESO 243-49 is the most luminous of ultraluminous X-ray sources and provides one of the strongest cases for the existence of intermediate mass black holes. We obtain good fits to disk-dominated observations of the source with BHSPEC, a fully relativistic black hole accretion disk spectral model. Due to degeneracies in the model arising from the lack of independent constraints on inclination and black hole spin, there is a factor of 100 uncertainty in the best-fit black hole mass M. Nevertheless, spectral fitting of XMM-Newton observations provides robust lower and upper limits with 3000 Msun < M < 3 x 10^5 Msun, at 90% confidence, placing HLX-1 firmly in the intermediate-mass regime. The lower bound on M is entirely determined by matching the shape and peak energy of the thermal component in the spectrum. This bound is consistent with (but independent of) arguments based solely on the Eddington limit. Joint spectral modelling of the XMM-Newton data with more luminous Swift and Chandra observations increases the lower bound to 6000 Msun, but this tighter constraint is not independent of the Eddington limit. The upper bound on M is sensitive to the maximum allowed inclination i, and is reduced to M < 10^5 Msun if we limit i < 75 deg.

In the magnetar model, the quiescent non-thermal soft X-ray emission from Anomalous X-ray Pulsars and Soft-Gamma Repeaters is thought to arise from resonant comptonization of thermal photons by charges moving in a twisted magnetosphere. Robust inference of physical quantities from observations is difficult, because the process depends strongly on geometry and current understanding of the magnetosphere is not very deep. The polarization of soft X-ray photons is an independent source of information, and its magnetospheric imprint remains only partially explored. In this paper we calculate how resonant cyclotron scattering would modify the observed polarization signal relative to the surface emission, using a multidimensional Monte Carlo radiative transfer code that accounts for the gradual coupling of polarization eigenmodes as photons leave the magnetosphere. We employ a globally-twisted, self-similar, force-free magnetosphere with a power-law momentum distribution, assume a blackbody spectrum for the seed photons, account for general relativistic light deflection close to the star, and assume that vacuum polarization dominates the dielectric properties of the magnetosphere. The latter is a good approximation if the pair multiplicity is not much larger than unity. Phase-averaged polarimetry is able to provide a clear signature of the magnetospheric reprocessing of thermal photons and to constrain mechanisms generating the thermal emission. Phase-resolved polarimetry, in addition, can characterize the spatial extent and magnitude of the magnetospheric twist angle at ~100 stellar radii, and discern between uni- or bidirectional particle energy distributions, almost independently of every other parameter in the system. We discuss prospects for detectability with GEMS.

The X-ray spectra of accretion discs of eight stellar-mass black holes have been analyzed to date using the thermal continuum fitting method, and the spectral fits have been used to estimate the spin parameters of the black holes. However, the underlying model used in this method of estimating spin is the general relativistic thin-disc model of Novikov & Thorne, which is only valid for razor-thin discs. We therefore expect errors in the measured values of spin due to inadequacies in the theoretical model. We investigate this issue by computing spectra of numerically calculated models of thin accretion discs around black holes, obtained via three-dimensional general relativistic magnetohydrodynamic (GRMHD) simulations. We apply the continuum fitting method to these computed spectra to estimate the black hole spins and check how closely the values match the actual spin used in the GRMHD simulations. We find that the error in the dimensionless spin parameter is up to about 0.2 for a non-spinning black hole, depending on the inclination. For black holes with spins of 0.7, 0.9 and 0.98, the errors are up to about 0.1, 0.03 and 0.01 respectively. These errors are comparable to or smaller than those arising from current levels of observational uncertainty. Furthermore, we estimate that the GRMHD simulated discs from which these error estimates are obtained correspond to effective disc luminosities of about 0.4-0.7 Eddington, and that the errors will be smaller for discs with luminosities of 0.3 Eddington or less, which are used in the continuum-fitting method. We thus conclude that use of the Novikov-Thorne thin-disc model does not presently limit the accuracy of the continuum-fitting method of measuring black hole spin.

A typical galaxy is thought to contain tens of millions of stellar-mass black holes, the collapsed remnants of once massive stars, and a single nuclear supermassive black hole. Both classes of black holes accrete gas from their environments. The accreting gas forms a flattened orbiting structure known as an accretion disk. During the past several years, it has become possible to obtain measurements of the spins of the two classes of black holes by modeling the X-ray emission from their accretion disks. Two methods are employed, both of which depend upon identifying the inner radius of the accretion disk with the innermost stable circular orbit (ISCO), whose radius depends only on the mass and spin of the black hole. In the Fe K method, which applies to both classes of black holes, one models the profile of the relativistically-broadened iron line with a special focus on the gravitationally redshifted red wing of the line. In the continuum-fitting method, which has so far only been applied to stellar-mass black holes, one models the thermal X-ray continuum spectrum of the accretion disk. We discuss both methods, with a strong emphasis on the continuum-fitting method and its application to stellar-mass black holes. Spin results for eight stellar-mass black holes are summarized. These data are used to argue that the high spins of at least some of these black holes are natal, and that the presence or absence of relativistic jets in accreting black holes is not entirely determined by the spin of the black hole.

The radiative efficiency of AGN is commonly estimated based on the total mass accreted and the total AGN light emitted per unit volume in the universe integrated over time (the Soltan argument). In individual AGN, thin accretion disk model spectral fits can be used to deduce the absolute accretion rate Mdot, if the black hole mass M is known. The radiative efficiency {\eta} is then set by the ratio of the bolometric luminosity L_bol to Mdot c^2. We apply this method to determine {\eta} in a sample of 80 PG quasars with well determined L_bol, where Mdot is set by thin accretion disk model fits to the optical luminosity density, and the M determination based on the bulge stellar velocity dispersion (13 objects) or the broad line region (BLR). For the BLR-based masses, we derive a mean log {\eta} = -1.05 +/- 0.52 consistent with the Soltan argument based estimates. We find a strong correlation of {\eta} with M, rising from {\eta} ~ 0.03 at M = 10^7 M{\odot} and L/L_Edd ~ 1 to {\eta} ~ 0.4 at M = 10^9 M{\odot} and L/L_Edd ~ 0.3. This trend is related to the overall uniformity of L_opt/L_bol in our sample, particularly the lack of the expected increase in L_opt/L_bol with increasing M (and decreasing L/L_Edd), which is a generic property of thermal disk emission at fixed {\eta}. The significant uncertainty in the M determination is not large enough to remove the correlation. The rising {\eta} with M may imply a rise in the black hole spin with M, as proposed based on other indirect arguments.

In prior work, {\it Chandra} and Gemini-North observations of the eclipsing X-ray binary M33 X-7 have yielded measurements of the mass of its black hole primary and the system's orbital inclination angle of unprecedented accuracy. Likewise, the distance to the binary is known to a few percent. In an analysis based on these precise results, fifteen {\it Chandra} and {\it XMM-Newton} X-ray spectra, and our fully relativistic accretion disk model, we find that the dimensionless spin parameter of the black hole primary is $a_* = 0.77 \pm 0.05$. The quoted 1-$\sigma$ error includes all sources of observational uncertainty. Four {\it Chandra} spectra of the highest quality, which were obtained over a span of several years, all lead to the same estimate of spin to within statistical errors (2%), and this estimate is confirmed by 11 spectra of lower quality. There are two remaining uncertainties: (1) the validity of the relativistic model used to analyze the observations, which is being addressed in ongoing theoretical work; and (2) our assumption that the black hole spin is approximately aligned with the angular momentum vector of the binary, which can be addressed by a future X-ray polarimetry mission.

The Suzaku observation of LMC X-3 gives the best data to date on the shape of the accretion disk spectrum. This is due to the combination of very low absorbing column density along this line of sight which allows the shape of the disk emisison to be constrained at low energies by the CCD's, while the tail can be simultaneously determined up to 30 keV by the high energy detectors. These data clearly demonstrate that the observed disk spectrum is broader than a simple 'sum of blackbodies', and relativistic smearing of the emission is strongly required. However, the intrinsic emission should be more complex than a (color-corrected) sum of blackbodies as it should also contain photo-electric absorption edges from the partially ionised disk photosphere. These are broadened by the relativistic smearing, but the models predict ~ 3-5 per cent deviations for 1/3- 1 solar abundance around the edge energies, significantly stronger than observed. This indicate that the models need to include more physical processes such as self-irradiation, bound-bound (line) absorption, and/or emission from recombination continuua and/or lines. Alternatively, if none of these match the data, it may instead require that the accretion disk density and/or emissivity profile with height is different to that assumed. Thus these data demonstrate the feasibility of observational tests of our fundamental understanding of the vertical structure of accretion disks.

We examine the effects of density stratification on magnetohydrodynamic turbulence driven by the magnetorotational instability in local simulations that adopt the shearing box approximation. Our primary result is that, even in the absence of explicit dissipation, the addition of vertical gravity leads to convergence in the turbulent energy densities and stresses as the resolution increases, contrary to results for zero net flux, unstratified boxes. The ratio of total stress to midplane pressure has a mean of ~0.01, although there can be significant fluctuations on long (>~50 orbit) timescales. We find that the time averaged stresses are largely insensitive to both the radial or vertical aspect ratio of our simulation domain. For simulations with explicit dissipation, we find that stratification extends the range of Reynolds and magnetic Prandtl numbers for which turbulence is sustained. Confirming the results of previous studies, we find oscillations in the large scale toroidal field with periods of ~10 orbits and describe the dynamo process that underlies these cycles.

We present the results of one and three-dimensional radiative transfer calculations of polarized spectra emerging from snapshots of radiation magnetohydrodynamical simulations of the local vertical structure of black hole accretion disks. The simulations cover a wide range of physical regimes relevant for the high/soft state of black hole X-ray binaries. We constrain the uncertainties in theoretical spectral color correction factors due to the presence of magnetic support of the disk surface layers and strong density inhomogeneities. For the radiation dominated simulation, magnetic support increases the color correction factor by about ten percent, but this is largely compensated by a ten percent softening due to inhomogeneities. We also compute the effects of inhomogeneities and Faraday rotation on the resulting polarization. Magnetic fields in the simulations are just strong enough to produce significant Faraday depolarization near the spectral peak of the radiation field. X-ray polarimetry may therefore be a valuable diagnostic of accretion disk magnetic fields, being able to directly test simulations of magnetorotational turbulence.

The first extragalactic X-ray binary, LMC X-1, was discovered in 1969. In the 1980s, its compact primary was established as the fourth dynamical black-hole candidate. Recently, we published accurate values for the mass of the black hole and the orbital inclination angle of the binary system. Building on these results, we have analyzed 53 X-ray spectra obtained by RXTE and, using a selected sample of 18 of these spectra, we have determined the dimensionless spin parameter of the black hole to be a* = 0.92(-0.07,+0.05). This result takes into account all sources of observational and model-parameter uncertainties. The standard deviation around the mean value of a* for these 18 X-ray spectra, which were obtained over a span of several years, is only 0.02. When we consider our complete sample of 53 RXTE spectra, we find a somewhat higher value of the spin parameter and a larger standard deviation. Finally, we show that our results based on RXTE data are confirmed by our analyses of selected X-ray spectra obtained by the XMM-Newton, BeppoSAX and Ginga missions.

In order for a modified gravity model to be a candidate for cosmological dark energy it has to pass stringent local gravity experiments. We find that a Brans-Dicke (BD) theory with well-defined second order corrections that include the Gauss-Bonnet term possess this feature. We construct the generic second order theory that gives, to linear order, a BD metric solution for a point-like mass source. We find that these theories interpolate between general relativity (GR) and BD gravity. In particular it is found that the relevant Eddington parameter, that is commonly heavily constrained by time delay experiments, can be arbitrarily close to the GR value of 1, with an arbitrary BD parameter. We find the region where the solution is stable to small timelike perturbations.

The accretion flow in the disk dominated state of black hole binaries has peak temperature and luminosity which vary together in such a way as to indicate an approximately constant emitting area. The association of this with the last stable orbit gives one of the few ways to estimate spin when the mass of the black hole is known. However, deriving this radius requires knowledge of how the disk spectrum is modified by radiative transfer through the vertical structure of the disk, as well as special and general relativistic effects on the propagation of this radiation. Here we investigate the extent to which differences in vertical structure change the derived disk spectra by calculating these for a range of different stress prescriptions. We find that at a given mass accretion rate the spectra are almost identical for accretion rates of L/L_Edd <~ 0.1. The spectra are remarkably similar even up to the highest luminosities considered (L/L_Edd ~ 0.6) as long as the stresses do not dissipate more than about 10 per cent of the gravitational energy above the effective photosphere. This is exceeded only by classic alpha disks with alpha >~ 0.1, but these models give spectral variation which is incompatible with existing data. Therefore, we conclude that disk spectral modelling can place interesting constraints on angular momentum transport, but still provide a robust estimate of the spin of the black hole.

Inflation and moduli stabilisation mechanisms work well independently, and many string-motivated supergravity models have been proposed for them. However a complete theory will contain both, and there will be (gravitational) interactions between the two sectors. These give corrections to the inflaton potential, which generically ruin inflation. This holds true even for fine-tuned moduli stabilisation schemes. Following a suggestion by 0712.3460, we show that a viable combined model can be obtained if it is the Kahler functions (G= K+\ln |W|^2) of the two sectors that are added, rather than the superpotentials (as is usually done). Interaction between the two sectors does still impose some restrictions on the moduli stabilisation mechanism, which are derived. Significantly, we find that the (post-inflation) moduli stabilisation scale no longer needs to be above the inflationary energy scale.

Context. In support of NASA's exploration program and the return to the Moon, the polarimetric signature of dispersed individual Lunar regolith dust grains is studied to enable the characterization of the dust exopsheric environment by remote, in-situ, and standoff sensing. Aims. We explore the value of the negative polarization branch (NPB) as a signature for characterizing individual grains to determine if it can be used in the same way as for surfaces of planets and atmosphereless bodies. Results. Calculations show that polarization phase curves for spherical grains exhibit a sharp transition over a narrow range of size parameter between two distinct regimes, one typical of Rayleigh scattering and another dominated by a large NPB. The linear polarimetric signature observed for aggregates is a composite of a) the polarization induced by individual grains composing the aggregate and b) the polarization due to the aggregate as a whole dust grain. The weight of each component varies depending on the porosity of the aggregate. An NPB similar to the one observed for atmosphereless astronomical bodies is present for different ranges of the size parameter depending on the value of the porosity. It appears as a remnant of the negative branch exhibited by the single spherical grains. The sharper, narrow negative branch that is measured for some granular surfaces in the laboratory or seen in astronomical observations is not observed here.

Chaotic inflation predicts a large gravitational wave signal which can be tested by the upcoming Planck satellite. We discuss a SUGRA implementation of chaotic inflation in the presence of moduli fields, and find that inflation does not work with a generic KKLT moduli stabilisation potential. A viable model can be constructed with a fine-tuned moduli sector, but only for a very specific choice of Kahler potential. Our analysis also shows that inflation models satisfying \partial_{i} W_{\rm inf}=0 for all inflation sector fields \phi_i can be combined successfully with a fine-tuned moduli sector.

Corrections to solar system gravity are derived for f(G) gravity theories, in which a function of the Gauss-Bonnet curvature term is added to the gravitational action. Their effects on Newton's law, as felt by the planets, and on the frequency shift of signals from the Cassini spacecraft, are both determined. Despite the fact that the Gauss-Bonnet term is quadratic in curvature, the resulting constraints are substantial. It is shown that they practically rule out f(G) as a natural explanation for the late-time acceleration of the universe. Possible exceptions are when f(G) reduces to something very close to a cosmological constant, or if the form of the function f is exceptionally fine-tuned.

We show that the luminosity of a star forming galaxy is capped by the production and subsequent expulsion of cosmic rays from its interstellar medium. By defining an Eddington luminosity in cosmic rays, we show that the star formation rate of a given galaxy is limited by its mass content and the cosmic ray mean free path. When the cosmic ray luminosity and pressure reaches a critical value as a result of vigorous star formation, hydrostatic balance is lost, a cosmic ray-driven wind develops, and star formation is choked off. Cosmic ray pressure-driven winds are likely to produce wind velocities significantly in excess of the galactic escape velocity. It is possible that cosmic ray feedback results in the Faber-Jackson relation for a plausible set of input parameters that describe cosmic ray production and transport, which are calibrated by observations of the Milky Way's interstellar cosmic rays as well as other galaxies.

A spectral index n_s < 0.95 appears to be a generic prediction of racetrack inflation models. Reducing a general racetrack model to a single-field inflation model with a simple potential, we obtain an analytic expression for the spectral index, which explains this result. By considering the limits of validity of the derivation, possible ways to achieve higher values of the spectral index are described, although these require further fine-tuning of the potential.

Although the Gauss-Bonnet term is a topological invariant for general relativity, it couples naturally to a quintessence scalar field, modifying gravity at solar system scales. We determine the solar system constraints due to this term by evaluating the post-Newtonian metric for a distributional source. We find a mass dependent, 1/r^7 correction to the Newtonian potential, and also deviations from the Einstein gravity prediction for light-bending. We constrain the parameters of the theory using planetary orbits, the Cassini spacecraft data, and a laboratory test of Newton's law, always finding extremely tight bounds on the energy associated to the Gauss-Bonnet term. We discuss the relevance of these constraints to late-time cosmological acceleration.

We present the CMD from deep HST imaging in the globular cluster NGC 6397. The ACS was used for 126 orbits to image a single field in two colors (F814W, F606W) 5 arcmin SE of the cluster center. The field observed overlaps that of archival WFPC2 data from 1994 and 1997 which were used to proper motion (PM) clean the data. Applying the PM corrections produces a remarkably clean CMD which reveals a number of features never seen before in a globular cluster CMD. In our field, the main sequence stars appeared to terminate close to the location in the CMD of the hydrogen-burning limit predicted by two independent sets of stellar evolution models. The faintest observed main sequence stars are about a magnitude fainter than the least luminous metal-poor field halo stars known, suggesting that the lowest luminosity halo stars still await discovery. At the bright end the data extend beyond the main sequence turnoff to well up the giant branch. A populous white dwarf cooling sequence is also seen in the cluster CMD. The most dramatic features of the cooling sequence are its turn to the blue at faint magnitudes as well as an apparent truncation near F814W = 28. The cluster luminosity and mass functions were derived, stretching from the turn off down to the hydrogen-burning limit. It was well modeled with either a very flat power-law or a lognormal function. In order to interpret these fits more fully we compared them with similar functions in the cluster core and with a full N-body model of NGC 6397 finding satisfactory agreement between the model predictions and the data. This exercise demonstrates the important role and the effect that dynamics has played in altering the cluster IMF.

We measure long (2200-4000 ang) and short (1450-2200 ang) wavelength spectral slopes \alpha (F_\nu proportional to \nu^\alpha) for quasar spectra from the Sloan Digital Sky Survey. The long and short wavelength slopes are computed from 3646 and 2706 quasars with redshifts in the z=0.76-1.26 and z=1.67-2.07 ranges, respectively. We calculate mean slopes after binning the data by monochromatic luminosity at 2200 ang and virial mass estimates based on measurements of the MgII line width and 3000 ang continuum luminosity. We find little evidence for mass dependent variations in the mean slopes, but a significant luminosity dependent trend in the near UV spectral slopes is observed with larger (bluer) slopes at higher luminosities. The far UV slopes show no clear variation with luminosity and are generally lower (redder) than the near UV slopes at comparable luminosities, suggesting a slightly concave quasar continuum shape. We compare these results with Monte Carlo distributions of slopes computed from models of thin accretion disks, accounting for uncertainties in the mass estimates. The model slopes produce mass dependent trends which are larger than observed, though this conclusion is sensitive to the assumed uncertainties in the mass estimates. The model slopes are also generally bluer than observed, and we argue that reddening by dust intrinsic to the source or host galaxy may account for much of the discrepancy.

NGC 6397 is the second closest globular star cluster to the Sun. Using 5 days of time on the Hubble Space Telescope, we have constructed the deepest ever color-magnitude diagram for this cluster. We see a clear truncation in each of its two major stellar sequences. Faint red main sequence stars run out well above our observational limit and near to the theoretical prediction for the lowest mass stars capable of stable hydrogen-burning in their cores. We also see a truncation in the number counts of faint blue stars, namely white dwarfs. This reflects the limit to which the bulk of the white dwarfs can cool over the lifetime of the cluster. There is also a turn towards bluer colors in the least luminous of these objects. This was predicted for the very coolest white dwarfs with hydrogen-rich atmospheres as the formation of H2 causes their atmospheres to become largely opaque to infrared radiation due to collision-induced absorption.

As a by-product of high-precision, ultra-deep stellar photometry in the Galactic globular cluster NGC 6397 with the Hubble Space Telescope, we are able to measure a large population of background galaxies whose images are nearly point-like. These provide an extragalactic reference frame of unprecedented accuracy, relative to which we measure the most accurate absolute proper motion ever determined for a globular cluster. We find mu_alpha = 3.56 +/- 0.04 mas/yr and mu_delta = -17.34 +/- 0.04 mas/yr. We note that the formal statistical errors quoted for the proper motion of NGC 6397 do not include possible unavoidable sources of systematic errors, such as cluster rotation. These are very unlikely to exceed a few percent. We use this new proper motion to calculate NGC 6397's UVW space velocity and its orbit around the Milky Way, and find that the cluster has made frequent passages through the Galactic disk.

Models with a scalar field coupled to the Gauss-Bonnet Lagrangian appear naturally from Kaluza-Klein compactifications of pure higher-dimensional gravity. We study linear, cosmological perturbations in the limits of weak coupling and slow-roll, and derive simple expressions for the main observable sub-horizon quantities: the anisotropic stress factor, the time-dependent gravitational constant, and the matter perturbation growth factor. Using present observational data, and assuming slow-roll for the dark energy field, we find that the fraction of energy density associated with the coupled Gauss-Bonnet term cannot exceed 15%. The bound should be treated with caution, as there are significant uncertainies in the data used to obtain it. Even so, it indicates that the future prospects for constraining the coupled Gauss-Bonnet term with cosmological observations are encouraging.

Turbulent Comptonization, a potentially important damping and radiation mechanism in relativistic accretion flows, is discussed. Particular emphasis is placed on the physical basis, relative importance, and thermodynamics of turbulent Comptonization. The effects of metal-absorption opacity on the spectral component resulting from turbulent Comptonization is considered as well.

Microquasars are galactic black hole binary systems with radio jets which can sometimes be spatially resolved to show superluminal motion. The first and best known of this class of objects is GRS 1915+105, the brightest accreting source in our Galaxy. There is persistent speculation that strong jet emission could be linked to black hole spin. If so, the high spin should also be evident in accretion disc spectra. We search the RXTE archive to find disc-dominated X-ray spectra from this object, as these are the only ones which can give reliable spin determinations by this method. Finding these is complicated by the rapid, unique limit cycle variability, but we are able to identify such spectra by going to the shortest possible time resolution (16 s). We fit them with a simple multicolour disc blackbody diskbb, and with the best current model which include full radiative transfer as well as relativistic effects bhspec. Both these models show that the spin is intermediate, neither zero nor maximal. bhspec, the most physical model, gives a value for the dimensionless spin of a* \~ 0.7 for a distance of 12.5 kpc and inclination of 66 deg. This, together with the range of spins 0.1 < a* < 0.8 derived using this method for other black holes, suggests that jet emission is probably fundamentally powered by gravity rather than spin, and implies that high-to-maximal spin is not a pre-requisite for powerful relativistic jets.

Based on a spectral analysis of the X-ray continuum that employs a fully relativistic accretion-disk model, we conclude that the compact primary of the binary X-ray source GRS 1915+105 is a rapidly-rotating Kerr black hole. We find a lower limit on the dimensionless spin parameter of a* greater than 0.98. Our result is robust in the sense that it is independent of the details of the data analysis and insensitive to the uncertainties in the mass and distance of the black hole. Furthermore, our accretion-disk model includes an advanced treatment of spectral hardening. Our data selection relies on a rigorous and quantitative definition of the thermal state of black hole binaries, which we used to screen all of the available RXTE and ASCA data for the thermal state of GRS 1915+105. In addition, we focus on those data for which the accretion disk luminosity is less than 30% of the Eddington luminosity. We argue that these low-luminosity data are most appropriate for the thin alpha-disk model that we employ. We assume that there is zero torque at the inner edge of the disk, as is likely when the disk is thin, although we show that the presence of a significant torque does not affect our results. Our model and the model of the relativistic jets observed for this source constrain the distance and black hole mass and could thus be tested by determining a VLBA parallax distance and improving the measurement of the mass function. Finally, we comment on the significance of our results for relativistic-jet and core-collapse models, and for the detection of gravitational waves.

We present results from spectral modeling of three black hole X-ray binaries: LMC X-3, GRO J1655-40, and XTE J1550-564. Using a sample of disk dominated observations, we fit the data with a range of spectral models that includes a simple multitemperature blackbody (DISKBB), a relativistic accretion disk model based on color-corrected blackbodies (KERRBB), and a relativistic model based on non-LTE atmosphere models within an alpha prescription (BHSPEC). BHSPEC provides the best fit for a BeppoSAX observation of LMC X-3, which has the broadest energy coverage of our sample. It also provides the best fit for multiple epochs of Rossi X-ray Timing Explorer (RXTE) data in this source, except at the very highest luminosity (L/L_{Edd} >~ 0.7), where additional physics must be coming into play. BHSPEC is also the best-fit model for multi-epoch RXTE observations of GRO J1655-40 and XTE J1550-564, although the best-fit inclination of the inner disk differs from the binary inclination. All our fits prefer alpha=0.01 to alpha=0.1, in apparent disagreement with the large stresses inferred from the rapid rise times observed in outbursts of these two sources. In all three sources our fits imply moderate black hole spins (a/M ~ 0.1 - 0.8), but this is sensitive to the reliability of independent measurements of these system parameters and to the physical assumptions which underly our spectral models.

Self-consistent vertical structure models together with non-LTE radiative transfer should produce spectra from accretion disks around black holes which differ from multitemperature blackbodies at levels which may be observed. High resolution, high signal-to-noise observations warrant spectral modeling which both accounts for relativistic effects, and treats the physics of radiative transfer in detail. In Davis et al. (2005) we presented spectral models which accounted for non-LTE effects, Compton scattering, and the opacities due to ions of abundant metals. Using a modification of this method, we have tabulated spectra for black hole masses typical of Galactic binaries. We make them publicly available for spectral fitting as an Xspec model. These models represent the most complete realization of standard accretion disk theory to date. Thus, they are well suited for both testing the theory's applicability to observed systems and for constraining properties of the black holes, including their spins.

Stellar atmosphere models of ionized accretion disks have generally neglected the contribution of magnetic fields to the vertical hydrostatic support, although magnetic fields are widely believed to play a critical role in the transport of angular momentum. Simulations of magnetorotational turbulence in a vertically stratified shearing box geometry show that magnetic pressure support can be dominant in the upper layers of the disk. We present calculations of accretion disk spectra that include this magnetic pressure support, as well as a vertical dissipation profile based on simulation. Magnetic pressure support generically produces a more vertically extended disk atmosphere with a larger density scale height. This acts to harden the spectrum compared to models that neglect magnetic pressure support. We estimate the significance of this effect on disk-integrated spectra by calculating an illustrative disk model for a stellar mass black hole, assuming that similar magnetic pressure support exists at all radii.

We fit X-ray spectral data in the thermal dominant or high soft state of two dynamically confirmed black holes, GRO J1655-40 and 4U 1543-47, and estimate the dimensionless spin parameters a* = a/M of the two holes. For GRO J1655-40, using a spectral hardening factor computed for a non-LTE relativistic accretion disk, we estimate a* ~ 0.75 and ~ 0.65-0.75, respectively, from ASCA and RXTE data. For 4U 1543-47, we estimate a* ~ 0.75-0.85 from RXTE data. Thus, neither black hole has a spin approaching the theoretical maximum a* = 1.

The radiation spectra of many of the brightest ultraluminous X-ray sources (ULXs) are dominated by a hard power law component, likely powered by a hot, optically thin corona that Comptonizes soft seed photons emitted from a cool, optically thick black hole accretion disk. Before its dissipation and subsequent conversion into coronal photon power, the randomized gravitational binding energy responsible for powering ULX phenomena must separate from the mass of its origin by a means other than, and quicker than, electron scattering-mediated radiative diffusion. Therefore, the release of accretion power in ULXs is not necessarily subject to Eddington-limited photon trapping, as long as it occurs in a corona. Motivated by these basic considerations, we present a model of ULXs powered by geometrically thin accretion onto stellar mass black holes. We argue that the radiative efficiency of the flow remains high if the corona is magnetized or optically thin and the majority of the accretion power escapes in the form of radiation rather than an outflow. Within the context of the current black hole X-ray binary paradigm, our ULX model may be viewed as an extension of the very high state observed in Galactic sources. (abridged)

We examine the effects of photon bubble instability in radiation-dominated accretion disks such as those found around black holes in active galactic nuclei and X-ray binary star systems. Two- and 3-D numerical radiation MHD calculations of small patches of disk are used. Modes with wavelengths shorter than the gas pressure scale height grow faster than the orbital frequency in the surface layers. The fastest growth rate observed is five times the orbital frequency and occurs on nearly-vertical magnetic fields. The spectrum of linear modes agrees with a WKB analysis indicating still faster growth at unresolved scales, with a maximum proportional to the gravity and inversely proportional to the gas sound speed. Disturbances reaching non-linear amplitudes steepen into trains of shocks similar to a 1-D periodic non-linear analytic solution. Variations in propagation speed result in merging of adjacent fronts, and over time the shock spacing and amplitude increase. Growth is limited by the strength of the field, and the structure is disrupted when the ram pressure exceeds the magnetic pressure. The largest horizontal density variations are similar to the ratio of magnetic to gas pressure, and in our calculations are more than 100. Under the conditions considered, radiation diffuses through the inhomogeneneous flow 5x faster than through the initial hydrostatic equilibrium, and the net cooling rate is several times greater than in a similar calculation without magnetic fields showing the effects of convection. These results indicate photon bubbles may be important in cooling radiation-dominated disks.

We present calculations of non-LTE, relativistic accretion disk models applicable to the high/soft state of black hole X-ray binaries. We include the effects of thermal Comptonization and bound-free and free-free opacities of all abundant ion species. We present spectra calculated for a variety of accretion rates, black hole spin parameters, disk inclinations, and stress prescriptions. We also consider nonzero inner torques on the disk, and explore different vertical dissipation profiles, including some which are motivated by recent radiation MHD simulations of magnetorotational turbulence. Bound-free metal opacity generally produces significantly less spectral hardening than previous models which only considered Compton scattering and free-free opacity. It also tends to keep the effective photosphere near the surface, resulting in spectra which are remarkably independent of the stress prescription and vertical dissipation profile, provided little dissipation occurs above the effective photosphere. We provide detailed comparisons between our models and the widely used multicolor disk model. Frequency dependent discrepancies exist that may affect the parameters of other spectral components when this simpler disk model is used to fit modern X-ray data. For a given source, our models predict that the luminosity in the high/soft state should approximately scale with the fourth power of the empirically inferred maximum temperature, but with a slight hardening at high luminosities. This is in good agreement with observations. (abridged)

We present XMM-Newton observations of the radio galaxy 3C 120. The hard X-ray spectrum contains a marginally resolved Fe I K-alpha emission line with FWHM=9,000 km/s and an equivalent width of 57 eV. The line arises via fluorescence in a broad-line region with covering fraction of 0.4. There is no evidence of relativistically broad Fe K-alpha, contrary to some previous reports. The normal equivalent widths of the X-ray and optical emission lines exclude a strongly beamed synchrotron component to the hard X-ray and optical continua. There is an excess of 0.3-2 keV soft X-ray continuum over an extrapolation of the hard X-ray power-law, which may arise in a disk corona. Analysis of an archival Chandra image shows that extended emission from the jet and other sources contributes <3% of the total X-ray flux. A break in the X-ray spectrum below 0.6 keV indicates an excess neutral hydrogen column density of N_H=1.57 * 10^21 cm^{-2}. However, the neutral absorber must have an oxygen abundance of <1/50 of the solar value to explain the absence of an intrinsic or intervening O I edge. There is no ionized absorption in the soft X-ray spectrum, but there is a weak, narrow O VIII Ly-alpha emission line. We do not detect previously claimed O VIII absorption from the intervening intergalactic medium. Radio observations at 37 GHz show a fast, high frequency flare starting 8 days after the XMM-Newton observation. However, this has no obvious effect on the X-ray spectrum. The X-ray spectrum, including the soft excess, became harder as the X-ray flux decreased, with an estimated pivot energy of 40 keV. The UV and soft X-ray fluxes are strongly correlated over the 120 ks duration of the XMM-Newton observation. This is qualitatively consistent with Comptonization of UV photons by a hot corona. (Abridged)

This paper presents new results from a photometric study of the main-sequence stars in M4, the closest globular cluster to the sun. Multi-field, multi-epoch observations at approximately 1, 2, and 6 core-radii were obtained with the WFPC2 camera on the HST. The multi-epoch observations allowed us to clean the data on the basis of proper motion. In all the fields the cluster main sequence can be traced to at least V = 27.0 but there remains a trail of stars to the limit of the data near V = 30 in the deepest outer field. There is no evidence that we have reached the end of the hydrogen burning main sequence in any of our fields, however, there is some indication that very few stars remain to be detected in the deepest data. The scatter about the cluster main sequence yields a surprisingly small and variable binary fraction; f_b = 2% in the inner parts of the cluster falling to the 1% range outside. For the currently visible main sequence stars, the cluster mass function is very flat (alpha =0.1) in the outer field and flattens further in the inner fields suggesting well developed mass segregation. The observed variation in the mass function is broadly consistent with isotropic, multi-mass Michie-King models. Because we have a large sample of white dwarfs in the outer field we are able to show that the cluster IMF above 0.8Msun was considerably steeper than the present day mass function for low mass stars. Two appendicies are included in this contribution; a detailed discussion of the techniques used to reduce the data and a direct comparison between the cluster stars and those belonging to the inner spheroid of the Galaxy yielding a Galactocentric distance of 7.5+/-0.6 kpc.

In this note we address the issue of hydrodynamical instabilities in Astrophysical rotating shear flows in the light of recent publications focused on the possibility for differential rotation to trigger and sustain turbulence in the absence of a magnetic field. We wish to present in a synthetic form the major arguments in favor of this thesis along with a simple schematic scenario of the transition to and self-sustenance of such turbulence. We also propose that the turbulent diffusion length scale scales as the local Rossby number of the mean flow. A new prescription for the turbulent viscosity is introduced. This viscosity reduces to the so-called beta-prescription in the case of velocity profiles with a constant Rossby number, which includes Keplerian rotating flows.

Numerical calculations and linear theory of radiation magnetohydrodynamic flows indicate that the photon bubble and magnetorotational instability (MRI) may produce large density inhomogeneities in radiation pressure supported media. We study the effects of the photon bubble instability on accretion disk spectra using 2-D Monte Carlo (MC) and 1-D Feautrier radiative transfer calculations on a snapshot of a 2-D numerical simulation domain. We find an enhancement in the thermalization of the MC spectra over that of the Feautrier calculation. In the inner-most regions of these disks, the turbulent magnetic pressure may greatly exceed that of the gas. It is then possible for bulk turbulent Alfvenic motions driven by the MRI to exceed the thermal velocity making turbulent Comptonization the dominant radiative process. We estimate the spectral distortion due to turbulent Comptonization utilizing a 1-D MC calculation.

In the inner-most regions of radiation pressure supported accretion disks, the turbulent magnetic pressure may greatly exceed that of the gas. If this is the case, it is possible for bulk Alfvenic motions driven by the magnetorotational instability (MRI) to surpass the electron thermal velocity. Bulk rather than thermal Comptonization may then be the dominant radiative process which mediates gravitational energy release. For sufficiently large turbulent stresses, we show that turbulent Comptonization produces a significant contribution to the far-UV and X-ray emission of black hole accretion disks. The existence of this spectral component provides a means of obtaining direct observational constraints on the nature of the turbulence itself. We describe how this component may affect the spectral energy distributions and variability properties of X-ray binaries and active galactic nuclei.

We examine the 3-10 keV EPIC spectra of Mrk 205 and Mrk 509 to investigate their Fe K features. The most significant feature in the spectra of both objects is an emission line at 6.4 keV. The spectra can be adequately modelled with a power law and a relatively narrow (sigma < 0.2 keV) Fe K alpha emission line. Better fits are obtained when an additional Gaussian emission line, relativistic accretion-disk line, or Compton reflection from cold material, is added to the spectral model. We obtain similar goodness of fit for any of these three models, but the model including Compton reflection from cold material offers the simplest, physically self-consistent solution, because it only requires one reprocessing region. Thus the Fe K spectral features in Mrk 205 and Mrk 509 do not present strong evidence for reprocessing in the inner, relativistic parts of accretion disks.

Competing models for broad spectral features in the soft X-ray spectrum of the Seyfert I galaxy Mrk766 are tested against data from a 130 ks XMM-Newton observation. A model including relativistically broadened Ly-alpha emission lines of OVIII, NVII and CVI is a better fit to 0.3-2 keV XMM RGS data than a dusty warm absorber. Moreover, the measured depth of neutral iron absorption lines in the spectrum is inconsistent with the magnitude of the iron edge required to produce the continuum break at 17-18Angstroms in the dusty warm absorber model. The relativistic emission line model can reproduce the broad-band (0.1-12 keV) XMM-EPIC data with the addition of a fourth line to represent emission from ionized iron at 6.7 keV and an excess due to reflection at energies above the iron line. The profile of the 6.7 keV iron line is consistent with that measured for the low energy lines. There is evidence in the RGS data at the 3sigma level for spectral features that vary with source flux. The covering fraction of warm absorber gas is estimated to be ~12%. Iron in the warm absorber is found to be overabundant with respect to CNO compared to solar values.

We present medium-resolution z-, J-, and H-band spectra of four late-type dwarfs with spectral types ranging from M8 to L7.5. In an attempt to determine the origin of numerous weak absorption features throughout their near-infrared spectra, and motivated by the recent tentative identification of the E 4\Pi- A ^4\Pi system of FeH near 1.6 microns in umbral and cool star spectra, we have compared the dwarf spectra to a laboratory FeH emission spectrum. We have identified nearly 100 FeH absorption features in the z-, J-, and H-band spectra of the dwarfs. In particular, we have identified 34 features which dominate the appearance of the H-band spectra of the dwarfs and which appear in the laboratory FeH spectrum. Finally, all of the features are either weaker or absent in the spectrum of the L7.5 dwarf which is consistent with the weakening of the known FeH bandheads in the spectra of the latest L dwarfs.

We report the results of a worldwide campaign to observe WZ Sagittae during its 2001 superoutburst. After a 23-year slumber at V=15.5, the star rose within 2 days to a peak brightness of 8.2, and showed a main eruption lasting 25 days. The return to quiescence was punctuated by 12 small eruptions, of ~1 mag amplitude and 2 day recurrence time; these "echo outbursts" are of uncertain origin, but somewhat resemble the normal outbursts of dwarf novae. After 52 days, the star began a slow decline to quiescence. Periodic waves in the light curve closely followed the pattern seen in the 1978 superoutburst: a strong orbital signal dominated the first 12 days, followed by a powerful /common superhump/ at 0.05721(5) d, 0.92(8)% longer than P_orb. The latter endured for at least 90 days, although probably mutating into a "late" superhump with a slightly longer mean period [0.05736(5) d]. The superhump appeared to follow familiar rules for such phenomena in dwarf novae, with components given by linear combinations of two basic frequencies: the orbital frequency omega_o and an unseen low frequency Omega, believed to represent the accretion disk's apsidal precession. Long time series reveal an intricate fine structure, with ~20 incommensurate frequencies. Essentially all components occurred at a frequency n(omega_o)-m(Omega), with m=1, ..., n. But during its first week, the common superhump showed primary components at n (omega_o)-Omega, for n=1, 2, 3, 4, 5, 6, 7, 8, 9 (i.e., m=1 consistently); a month later, the dominant power shifted to components with m=n-1. This may arise from a shift in the disk's spiral-arm pattern, likely to be the underlying cause of superhumps. The great majority of frequency components ... . (etc., abstract continues)

In brane world scenarios the Friedmann equation is modified, resulting in an increased expansion at early times. This has important effects on cosmological phase transitions which we investigate, elucidating significant differences to the standard case. First order phase transitions require a higher nucleation rate to complete; baryogenesis and particle abundances could be suppressed. Topological defect evolution is also affected, though the current defect densities are largely unchanged. In particular, the increased expansion does not solve the usual monopole and domain wall problems.

The Friedmann equation for a positive tension brane situated between two bulk spacetimes that posses the same 5D cosmological constant, but which does not posses a Z_2 symmetry of the metric itself is derived, and the possible effects of dropping the Z_2 symmetry on the expansion of our Universe are examined; cosmological constraints are discussed. The global solutions for the metric in the infinite extra dimension case are found and comparison with the symmetric case is made. We show that any brane world senario of this type must revert to a Z_2 symmetric form at late times, and hence rule out certain proposed scenarios.

The gravitational properties of cosmic strings in the Randall-Sundrum brane world scenario are investigated. Using a gauge in which the brane remains straight, the leading order corrections to the metric on the brane are determined. In contrast to their non-brane equivalents, these strings have an attractive 1/r^2 potential and a radially dependent deficit angle. These two effects alter many cosmological properties of the string, such as the formation of double images. As a result of the attractive force the string will collect matter as it moves through space.

It is shown that it is possible for bound fermions on a cosmic string to form a superconducting state. Due to the attractive force between them, particles moving in opposite directions along the string form bound pairs. This involves a similar mechanism to superconductivity in metals at low temperatures. The method of Gorkov is used to analyse the system. In contrast to the situation in metals, the unusual properties of the string fermion spectrum allow a massless Abelian gauge field to provide the required attractive force. This results in far stronger superconductivity than usual. A massive gauge field can also be used, in which case the standard results apply.

We investigate the microphysics of supersymmetric cosmic strings. In particular we focus on the vortices admitted by N=1 supersymmetric abelian Higgs models. We find the vortex solutions and demonstrate that the two simplest supersymmetric cosmic string models admit fermionic superconductivity. Further, by using supersymmetry transformations, we show how to solve for the fermion zero modes giving rise to string superconductivity in terms of the background string fields.