We report evidence for a quasi-periodic oscillation (QPO) in the optical light curve of KIC 9650712, a narrow-line Seyfert 1 galaxy in the original Kepler field. After the development and application of a pipeline for Kepler data specific to active galactic nuclei (AGN), one of our sample of 21 AGN selected by infrared photometry and X-ray flux demonstrates a peak in the power spectrum at log $\nu$ = -6.58 Hz, corresponding to a temporal period of t=44 days. We note that although the power spectrum is well-fit by a model consisting of a Lorentzian and a single power law, alternative continuum models cannot be ruled out. From optical spectroscopy, we measure the black hole mass of this AGN as log $M_{\mathrm{BH}} / M_\odot$ = 8.17. We find that this frequency lies along a correlation between low-frequency QPOs and black hole mass from stellar and intermediate mass black holes to AGN, similar to the known correlation in high-frequency QPOs.

Self-organization is a property of dissipative nonlinear processes that are governed by an internal driver and a positive feedback mechanism, which creates regular geometric and/or temporal patterns and decreases the entropy, in contrast to random processes. Here we investigate for the first time a comprehensive number of 16 self-organization processes that operate in planetary physics, solar physics, stellar physics, galactic physics, and cosmology. Self-organizing systems create spontaneous {\sl order out of chaos}, during the evolution from an initially disordered system to an ordered stationary system, via quasi-periodic limit-cycle dynamics, harmonic mechanical resonances, or gyromagnetic resonances. The internal driver can be gravity, rotation, thermal pressure, or acceleration of nonthermal particles, while the positive feedback mechanism is often an instability, such as the magneto-rotational instability, the Rayleigh-B\'enard convection instability, turbulence, vortex attraction, magnetic reconnection, plasma condensation, or loss-cone instability. Physical models of astrophysical self-organization processes involve hydrodynamic, MHD, and N-body formulations of Lotka-Volterra equation systems.

The Large Synoptic Survey Telescope is designed to provide an unprecedented optical imaging dataset that will support investigations of our Solar System, Galaxy and Universe, across half the sky and over ten years of repeated observation. However, exactly how the LSST observations will be taken (the observing strategy or "cadence") is not yet finalized. In this dynamically-evolving community white paper, we explore how the detailed performance of the anticipated science investigations is expected to depend on small changes to the LSST observing strategy. Using realistic simulations of the LSST schedule and observation properties, we design and compute diagnostic metrics and Figures of Merit that provide quantitative evaluations of different observing strategies, analyzing their impact on a wide range of proposed science projects. This is work in progress: we are using this white paper to communicate to each other the relative merits of the observing strategy choices that could be made, in an effort to maximize the scientific value of the survey. The investigation of some science cases leads to suggestions for new strategies that could be simulated and potentially adopted. Notably, we find motivation for exploring departures from a spatially uniform annual tiling of the sky: focusing instead on different parts of the survey area in different years in a "rolling cadence" is likely to have significant benefits for a number of time domain and moving object astronomy projects. The communal assembly of a suite of quantified and homogeneously coded metrics is the vital first step towards an automated, systematic, science-based assessment of any given cadence simulation, that will enable the scheduling of the LSST to be as well-informed as possible.

Diskoseismology, the theoretical study of normal mode oscillations in geometrically thin, optically thick accretion disks, is a strong candidate to explain some QPOs in the power spectra of many black hole X-ray binary systems. The existence of g-modes, presumably the most robust and visible of the modes, depends on general relativistic gravitational trapping in the hottest part of the disk. As the existence of the required cavity in the presence of magnetic fields has been put into doubt by theoretical calculations, we will explore in greater generality what the inclusion of magnetic fields has to say on the existence of g-modes. We use an analytical perturbative approach on the equations of MHD to assess the impact of such effects. Our main conclusion is that there appears to be no compelling reason to discard g-modes. In particular, the inclusion of a non-zero {\it radial} component of the magnetic field enables a broader scenario for cavity non-destruction, especially taking into account recent simulations' saturation values for the magnetic field.

We compare the determinations of the angular momentum of stellar mass black holes via the continuum and line methods with those from diskoseismology. The assumption being tested is that one of the QPOs (quasi-periodic oscillations) in each binary X-ray source is produced by the fundamental g-mode. This should be the most robust and visible normal mode of oscillation of the accretion disk, and therefore its absence should rule out diskoseismology as the origin of QPOs. The comparisons are consistent with the second highest frequency QPO being produced by this g-mode, but are not consistent with models in which one QPO frequency is that of the innermost stable circular orbit.

The relationship between the metric and nonrelativistic matter distribution depends on the theory of gravity and additional fields, providing a possible way of distinguishing competing theories. With the assumption that the geometry and kinematics of the homogeneous universe have been measured to sufficient accuracy, we present a procedure for understanding and testing the relationship between the cosmological matter distribution and metric perturbations (along with their respective evolution) using the ratio of the physical size of the perturbation to the size of the horizon as our small expansion parameter. We expand around Newtonian gravity on linear, subhorizon scales with coefficient functions in front of the expansion parameter. Our framework relies on an ansatz which ensures that (i) the Poisson equation is recovered on small scales (ii) the metric variables (and any additional fields) are generated and supported by the nonrelativistic matter overdensity. The scales for which our framework is intended are small enough so that cosmic variance does not significantly limit the accuracy of the measurements and large enough to avoid complications from nonlinear effects and baryon cooling. The coefficient functions provide a general framework for contrasting the consequences of Lambda CDM and its alternatives. We calculate the coefficient functions for general relativity with a cosmological constant and dark matter, GR with dark matter and quintessence, scalar-tensor theories, f(R) gravity and braneworld models. We identify a possibly unique signature of braneworld models. Constraining the coefficient functions provides a streamlined approach for testing gravity in a scale dependent manner. We briefly discuss the observations best suited for an application of our framework.

We demonstrate that the corotation resonance affects only some non-axisymmetric g-mode oscillations of thin accretion disks, since it is located within their capture zones. Using a more general (weaker radial WKB approximation) formulation of the governing equations, such g-modes, treated as perfect fluid perturbations, are shown to formally diverge at the position of the corotation resonance. A small amount of viscosity adds a small imaginary part to the eigenfrequency which has been shown to induce a secular instability (mode growth) if it acts hydrodynamically. The g-mode corotation resonance divergence disappears, but the mode magnitude can remain largest at the place of the corotation resonance. For the known g-modes with moderate values of the radial mode number and axial mode number (and any vertical mode number), the corotation resonance lies well outside their trapping region (and inside the innermost stable circular orbit), so the observationally relevant modes are unaffected by the resonance. The axisymmetric g-mode has been seen by Reynolds & Miller in a recent inviscid hydrodynamic accretion disk global numerical simulation. We also point out that the g-mode eigenfrequencies are approximately proportional to m for axial mode numbers |m|>0.

This paper studies the hydrodynamical problem of normal modes of small adiabatic oscillations of relativistic barotropic thin accretion disks around black holes (and compact weakly magnetic neutron stars). Employing WKB techniques, we obtain the eigenfrequencies and eigenfunctions of the modes for different values of the mass and angular momentum of the central black hole. We discuss the properties of the various types of modes and show for the first time that modes covering the whole extension of the disk (full--disk p--modes) can exist within the studied thin disk model. However, these modes have relatively short wavelengths.

Diskoseismology, the theoretical study of small adiabatic hydrodynamical global perturbations of geometrically thin, optically thick accretion disks around black holes (and other compact objects), is a potentially powerful probe of the gravitational field. For instance, the frequencies of the normal mode oscillations can be used to determine the elusive angular momentum parameter of the black hole. The general formalism developed by diskoseismologists for relativistic systems can be readily applied to the Newtonian case of cataclysmic variables (CVs). Some of these systems (e.g., the dwarf nova SS Cygni) show rapid oscillations in the UV with periods of tens of seconds and high coherence. In this paper, we assess the possibility that these dwarf nova oscillations (DNOs) are diskoseismic modes. Besides its importance in investigating the physical origin of DNOs, the present work could help us to answer the following question. To what extent are the similarities in the oscillation phenomenology of CVs and X-ray binaries (XRBs) indicative of a common physical mechanism?

We extend our investigation of the normal modes of small adiabatic oscillations of relativistic barotropic thin accretion disks to the inertial-pressure (p) modes. We focus here on the lowest frequency fundamental p-modes, those with no axial or vertical nodes in their distribution. Through a variety of analyses, we obtain closed-form expressions for the eigenfrequencies and eigenfunctions. These depend on the luminosity and viscosity parameter of the disk as well as the mass and angular momentum of the black hole via detailed formulae for the speed of sound. The effect of a torque on the inner edge of the disk is also included. We compare the p-mode properties to those of the g- and c-modes.

The coupled evolution of the spin frequency, core temperature, and r-mode amplitude of an accreting neutron star is calculated. We focus on those conditions that can produce persistent gravitational radiation from the r-mode. During X-ray quiescent phases of transient LMXBs, one may be able to identify the constant contribution of the gravitational wave emission to the spindown rate. Another signature is the r-mode contribution to the heating.

The gravitational-wave and accretion driven evolution of the angular velocity, core temperature, and (small) amplitude of an r-mode of neutron stars in low mass X-ray binaries and similar systems is investigated. The conditions required for evolution to a stable equilibrium state (with gravitational wave flux proportional to average X-ray flux) are determined. In keeping with conclusions derived from observations of neutron star cooling, the core neutrons are taken to be normal while the core protons and hyperons and the crust neutrons are taken to be singlet superfluids. The dominant sources of damping are then hyperon bulk viscosity (if much of the core is at least 2--3 times nuclear density) and (e-e and n-n) shear (and possibly magnetic) viscosity within the core--crust boundary layer. It is found that a stable equilibrium state can be reached if the superfluid transition temperature of the hyperons is sufficiently small ($\lesssim 2\times 10^9$ K), allowing the gravitational radiation from Sco X-1 and several other neutron stars in low-mass X-ray binaries to be potentially detectable by the second generation LIGO (and VIRGO) arrays.

We compare our calculations of the frequencies of the fundamental g, c, and p--modes of relativistic thin accretion disks with recent observations of high frequency QPOs in X-ray binaries with black hole candidates. These classes of modes encompass all adiabatic perturbations of such disks. The frequencies of these modes depend mainly on only the mass and angular momentum of the black hole; their weak dependence on disk luminosity is also explicitly indicated. Identifying the recently discovered relatively stable QPO pairs with the fundamental g and c modes provides a determination of the mass and angular momentum of the black hole. For GRO J1655-40, $M=5.9\pm 1.0 M_\sun$, $J=(0.917\pm 0.024)GM^2/c$, in agreement with spectroscopic mass determinations. For GRS 1915+105, $M=42.4\pm 7.0 M_\sun$, $J=(0.926\pm 0.020)GM^2/c$ or (less favored) $M=18.2\pm 3.1 M_\sun$, $J=(0.701\pm 0.043)GM^2/c$. We briefly address the issues of the amplitude, frequency width, and energy dependence of these QPOs.

The gravitational-wave and accretion driven evolution of neutron stars in low mass X-ray binaries and similar systems is analyzed, while the amplitude of the radiating perturbation (here assumed to be an r-mode) remains small. If most of the star is superfluid, with (temperature independent) mutual friction dominating the ordinary (temperature dependent) shear viscosity, the amplitude of the mode and the angular velocity of the star oscillate about their equilibrium values with a period of at least a few hundred years. The resulting oscillation of the neutron star temperature is also computed. For temperature dependent viscosity, the general conditions for the equilibrium to be stable are found.

We report the development of a new tool to determine the origin of quasi-periodic oscillations (QPOs) in accretion disk systems. The technique uses the source energy spectrum and the energy dependence of the QPO fractional amplitude to restrict the location of the emission region of the modulated photons, which are assumed to originate in the inner accretion disk. Both Monte Carlo and semi-analytical methods are presented. We assume the accretion disk is enshrouded by a slab atmosphere of hot electrons in which unsaturated Compton scattering produces the high-energy spectrum. Properties of the atmosphere, in particular the electron temperature, are assumed functions of radius from the central compact object. We show that our model reproduces the observed energy dependence of the fractional amplitude of the 67 Hz QPO in GRS 1915+105 if the QPO is assumed to originate within five gravitational radii from the innermost stable circular orbit.

We first briefly review how we investigate the modes of oscillation trapped within the inner region of accretion disks by the strong-field gravitational properties of a black hole (or a compact, weakly-magnetized neutron star). Then we focus on the `corrugation'(c)-modes, nearly incompressible perturbations of the inner disk. The fundamental c-modes have eigenfrequencies (ordered by radial mode number) which correspond to the Lense-Thirring frequency, evaluated at the outer trapping radius of the mode, in the slow rotation limit. This trapping radius is a decreasing function of the black hole angular momentum, so a significant portion of the disk is modulated only for slowly rotating black holes. The eigenfrequencies are thus strongly increasing functions of black hole angular momentum. The dependence of the eigenfrequencies on the speed of sound within (or the luminosity) within the disk is very weak, except for slowly rotating black holes.

The perturbations of weakly-viscous, barotropic, non-self-gravitating, Newtonian rotating fluids are analyzed via a single partial differential equation. The results are then used to find an expression for the viscosity-induced normal-mode complex eigenfrequency shift, with respect to the case of adiabatic perturbations. However, the effects of viscosity are assumed to have been incorporated in the unperturbed (equilibrium) model. This paper is an extension of the normal-mode formalism developed by Ipser & Lindblom for adiabatic pulsations of purely-rotating perfect fluids. The formulas derived are readily applicable to the perturbations of thin and thick accretion disks. We provide explicit expressions for thin disks, employing results from previous relativistic analyses of adiabatic normal modes of oscillation. In this case, we find that viscosity causes the fundamental p- and g- modes to grow while the fundamental c-mode could have either sign of the damping rate.

Compton scattering by a locally plane parallel atmosphere is analyzed within the Kompaneets approximation. In contrast to almost all previous analytic calculations of inverse Compton scattering, we assume neither homogeneity nor isotropy of the photon distribution. However, we do adopt an assumed angular distribution of the photons. For any specified energy dependence of the flux entering the base of the atmosphere, an analytic expression is obtained for the flux at any optical depth. Only integrals over the input flux are involved. If the thickness of the atmosphere is much less than its radial position on an accretion disk, the contribution from all radii may be summed. The initial motivation (addressed in a subsequent paper) is the computation of the energy spectrum of (possibly modulated) photons originating from a specified radial region of an accretion disk, relative to the total energy spectrum.

Synthetic spectra generated with the parameterized supernova synthetic-spectrum code SYNOW are compared to observed photospheric-phase spectra of the Type Ic supernova 1994I. The observed optical spectra can be well matched by synthetic spectra that are based on the assumption of spherical symmetry. We consider the identification of the infrared absorption feature observed near 10,250 \AA, which previously has been attributed to He I $\lambda10830$ and regarded as strong evidence that SN 1994I ejected some helium. We have difficulty accounting for the infrared absorption with He I alone. It could be a blend of He I and C I lines. Alternatively, we find that it can be fit by Si I lines without compromising the fit in the optical region. In synthetic spectra that match the observed spectra, from 4 days before to 26 days after the time of maximum brightness, the adopted velocity at the photosphere decreases from 17,500 to 7000 \kms. Simple estimates of the kinetic energy carried by the ejected mass give values that are near the canonical supernova energy of $10^{51}$ ergs. The velocities and kinetic energies that we find for SN 1994I in this way are much lower than those that we find elsewhere for the peculiar Type Ic SNe 1997ef and 1998bw, which therefore appear to have been hyper-energetic.

We study the asymptotic behavior at late times of Friedmann-Robertson-Walker (uniform density) cosmological models within scalar-tensor theories of gravity. Particularly, we analyze the late time behavior in the present (matter dominated) epoch of the universe. The result of Damour and Nordtvedt that for a massless scalar in a flat cosmology the Universe evolves towards a state indistinguishable from general relativity is generalized. We first study a massless scalar field in an open universe. It is found that, while the universe tends to approach a state with less scalar contribution to gravity, the attractor mechanism is not effective enough to drive the theory towards a final state indistinguishable from general relativity. For the self-interacting case it is found that the scalar field potential dominates the late time behavior. In most cases this makes the attractor mechanism effective, thus resulting in a theory of gravity with vanishingly small scalar contribution even for the open Universe.

We will summarize results of calculations of the modes of oscillation trapped within the inner region of accretion disks by the strong-field gravitational properties of a black hole (or a compact, weakly-magnetized neutron star). Their driving and damping will also be addressed. The focus will be on the most observable class: the analogue of internal gravity modes in stars. Their frequencies which corrrespond to the lowest mode numbers depend almost entirely upon only the mass and angular momentum of the black hole. Such a feature may have been detected in the X-ray power spectra of two galactic `microquasars', allowing the angular momentum of the black hole to be determined in one case.

We study the cosmological evolution of massless single-field scalar-tensor theories of gravitation from the time before the onset of $e^+e^-$ annihilation and nucleosynthesis up to the present. The cosmological evolution together with the observational bounds on the abundances of the lightest elements (those mostly produced in the early universe) place constraints on the coefficients of the Taylor series expansion of $a(\phi)$, which specifies the coupling of the scalar field to matter and is the only free function in the theory. In the case when $a(\phi)$ has a minimum (i.e., when the theory evolves towards general relativity) these constraints translate into a stronger limit on the Post-Newtonian parameters $\gamma$ and $\beta$ than any other observational test. Moreover, our bounds imply that, even at the epoch of annihilation and nucleosynthesis, the evolution of the universe must be very close to that predicted by general relativity if we do not want to over- or underproduce $^{4}$He. Thus the amount of scalar field contribution to gravity is very small even at such an early epoch.

The Rossi X-ray Timing Explorer (RXTE) has made feasible for the first time the search for high-frequency (~ 100 Hz) periodic features in black hole candidate (BHC) systems. Such a feature, with a 67 Hz frequency, recently has been discovered in the BHC GRS 1915+105 (Morgan, Remillard, & Greiner). This feature is weak (rms variability ~0.3%-1.6%), stable in frequency (to within ~2 Hz) despite appreciable luminosity fluctuations, and narrow (quality factor Q ~ 20). Several of these properties are what one expects for a ``diskoseismic'' g-mode in an accretion disk about a 10.6 M_sun (nonrotating) - 36.3 M_sun (maximally rotating) black hole (if we are observing the fundamental mode frequency). We explore this possibility by considering the expected luminosity modulation, as well as possible excitation and growth mechanisms---including turbulent excitation, damping, and ``negative'' radiation damping. We conclude that a diskoseismic interpretation of the observations is viable.

We generalize previous calculations to a fully relativistic treatment of adiabatic oscillations which are trapped in the inner regions of accretion disks by non-Newtonian gravitational effects of a black hole. We employ the Kerr geometry within the scalar potential formalism of Ipser and Lindblom, neglecting the gravitational field of the disk. This approach treats perturbations of arbitrary stationary, axisymmetric, perfect fluid models. It is applied here to thin accretion disks. Approximate analytic eigenfunctions and eigenfrequencies are obtained for the most robust and observable class of modes, which corresponds roughly to the gravity (internal) oscillations of stars. The dependence of the oscillation frequencies on the mass and angular momentum of the black hole is exhibited. These trapped modes do not exist in Newtonian gravity, and thus provide a signature and probe of the strong-field structure of black holes. Our predictions are relevant to observations which could detect modulation of the X-ray luminosity from stellar mass black holes in our galaxy and the UV and optical luminosity from supermassive black holes in active galactic nuclei.

The high state of black hole candidates is characterized by a quasi- thermal emission component at $kT \sim 1$ keV. In addition, this state tends to have very low variability which indicates that it is relatively stable, at least on {\it short} time scales. Most models of the high state imply that the bulk of the emission comes from an optically thick accretion disc; therefore, this state may be an excellent laboratory for testing our ideas about the physics of accretion discs. In this work we consider the implications of assuming that accretion disc viscosity arises from some form of turbulence. Specifically, we consider the simple case of three dimensional hydrodynamic turbulence. It is found that the coupling of such turbulence to acoustic modes in the disc can alter the disc emission. We calculate the amplitude and frequencies of this modulation, and we express our results in terms of the X-ray power spectral density. We compare our calculations with observations of the black hole candidate GS 1124-683, and show that for certain parameters we can reproduce some of the high frequency power. We then briefly explore mechanisms for producing the low frequency power, and note the difficulty that a single variability mechanism has in reproducing the full range of observed variability. In addition, we outline ways in which future spacecraft missions -- such as USA and XTE -- can further constrain our model, especially at frequencies above $\sim 10^2$ Hz.

We extend the approximate radiative transfer analysis of Hershkowitz, Linder, and Wagoner (1986) to a more general class of supernova model atmospheres, using a simple fit to the effective continuum opacity produced by lines (Wagoner, Perez, and Vasu 1991). At the low densities considered, the populations of the excited states of hydrogen are governed mainly by photoionization and recombination, and scattering dominates absorptive opacity. We match the asymptotic expressions for the spectral energy density $\Jnu$ at the photosphere, whose location at each frequency is determined by a first-order calculation of the deviation of $\Jnu$ from the Planck function $\Bnu$. The emergent spectral luminosity then assumes the form $ L_\nu = 4\pi^2 \r*^2 \zeta^2 \Bnu(T_p) $, where $T_p(\nu)$ is the photospheric temperature, $\zeta$ is the dilution factor, and $\r*$ is a fiducial radius [ultimately taken to be the photospheric radius $r_p (\nu)$]. The atmosphere is characterized by an effective temperature $T_e$ ($\propto L^{1/4} r_*^{-1/2}$) and hydrogen density $n_H=n_* (r_*/r)^\alpha$ ; and less strongly by the heavy element abundance and velocity gradient. Our major result is the dependence of $\zeta$ on frequency $\nu$ and the parameters $T_p$, $r_p$, and $\alpha$.

We study scalar perturbations to a Robertson-Walker cosmological metric in terms of a pseudo-Newtonian potential, which emerges naturally from the solution of the field equations. This potential is given in terms of a Green function for matter density fluctuations of arbitrary amplitude whose time and spatial dependence are assumed known. The results obtained span both the linearized and Newtonian limits, and do not explicitly depend on any kind of averaging procedure, but make the valid assumption that the global expansion rate is that of a Friedmann-Robertson-Walker model. In addition, we discuss the similarity to diffusive processes in the evolution of the potential, and possible applications.