Heating and Cooling of Hot Accretion Flows by Non Local Radiation

We consider non-local effects which arise when radiation emitted at one radius of an accretion disk either heats or cools gas at other radii through Compton scattering. We discuss three situations: 1. Radiation from the inner regions of an advection-dominated flow Compton cooling gas at intermediate radii and Compton heating gas at large radii. 2. Soft radiation from an outer thin accretion disk Compton cooling a hot one- or two-temperature flow on the inside. 3. Soft radiation from an inner thin accretion disk Compton cooling hot gas in a surrounding one-temperature flow. We describe how previous results are modified by these non-local interactions. We find that Compton heating or cooling of the gas by the radiation emitted in the inner regions of a hot flow is not important. Likewise, Compton cooling by the soft photons from an outer thin disk is negligible when the transition from a cold to a hot flow occurs at a radius greater than some minimum $R_{tr,min}$. However, if the hot flow terminates at $R < R_{tr,min}$, non-local cooling is so strong that the hot gas is cooled to a thin disk configuration in a runaway process. In the case of a thin disk surrounded by a hot one-temperature flow, we find that Compton cooling by soft radiation dominates over local cooling in the hot gas for \dot{M} \gsim 10^{-3} \alpha \dot{M}_{Edd}, and R \lsim 10^4 R_{Schw}. As a result, the maximum accretion rate for which an advection-dominated one-temperature solution exists, decreases by a factor of $\sim 10$, compared to the value computed under an assumption of local energy balance.Comment: LaTeX aaspp.sty, 25 pages, and 6 figures; to appear in Ap

Formal Specification and Testing of a Management Architecture

The importance of network and distributed systems management to supply and maintain services required by users has led to a demand for management facilities. Open network management is assisted by representing the system resources to be managed as objects, and providing standard services and protocols for interrogating and manipulating these objects. This paper examines the application of formal description techniques to the specification of managed objects by presenting a case study in the specification and testing of a management architecture. We describe a formal specification of a management architecture suitable for scheduling and distributing services across nodes in a distributed system. In addition, we show how formal specifications can be used to generate conformance tests for the management architecture

X-ray Images of Hot Accretion Flows

We consider the X-ray emission due to bremsstrahlung processes from hot, low radiative-efficiency accretion flows around supermassive and galactic black holes. We calculate surface brightness profiles and Michelson visibility functions for a range of density profiles, rho ~ r^(-3/2+p), with 0 < p < 1, to allow for the presence of outflows. We find that although the 1 keV emitting region in these flows can always extend up to 10^6 Schwarzschild radii (R_S), their surface brightness profiles and visibility functions are strongly affected by the specific density profile. The advection-dominated solutions with no outflows (p=0) lead to centrally peaked profiles with characteristic sizes of only a few tens of R_S. Solutions with strong outflows (p~1) lead to flat intensity profiles with significantly larger characteristic sizes of up to 10^6 R_S. This implies that low luminosity galactic nuclei, such as M87, may appear as extended X-ray sources when observed with current X-ray imaging instruments. We show that X-ray brightness profiles and their associated visibility functions may be powerful probes for determining the relevant mode of accretion and, in turn, the properties of hot accretion flows. We discuss the implications of our results for observations with the Chandra X-ray Observatory and the planned X-ray interferometer MAXIM.Comment: 14 pages, 4 figures, accepted by The Astrophysical Journal, minor change

Reaction–diffusion chemistry implementation of associative memory neural network

Unconventional computing paradigms are typically very difficult to program. By implementing efficient parallel control architectures such as artificial neural networks, we show that it is possible to program unconventional paradigms with relative ease. The work presented implements correlation matrix memories (a form of artificial neural network based on associative memory) in reaction–diffusion chemistry, and shows that implementations of such artificial neural networks can be trained and act in a similar way to conventional implementations

Distinguishing Solar Flare Types by Differences in Reconnection Regions

Observations show that magnetic reconnection and its slow shocks occur in solar flares. The basic magnetic structures are similar for long duration event (LDE) flares and faster compact impulsive (CI) flares, but the former require less non-thermal electrons than the latter. Slow shocks can produce the required non-thermal electron spectrum for CI flares by Fermi acceleration if electrons are injected with large enough energies to resonate with scattering waves. The dissipation region may provide the injection electrons, so the overall number of non-thermal electrons reaching the footpoints would depend on the size of the dissipation region and its distance from the chromosphere. In this picture, the LDE flares have converging inflows toward a dissipation region that spans a smaller overall length fraction than for CI flares. Bright loop-top X-ray spots in some CI flares can be attributed to particle trapping at fast shocks in the downstream flow, the presence of which is determined by the angle of the inflow field and velocity to the slow shocks.Comment: 15 pages TeX and 2 .eps figures, accepted to Ap.J.Let

Are Particles in Advection-Dominated Accretion Flows Thermal?

We investigate the form of the momentum distribution function for protons and electrons in an advection-dominated accretion flow (ADAF). We show that for all accretion rates, Coulomb collisions are too inefficient to thermalize the protons. The proton distribution function is therefore determined by the viscous heating mechanism, which is unknown. The electrons, however, can exchange energy quite efficiently through Coulomb collisions and the emission and absorption of synchrotron photons. We find that for accretion rates greater than \sim 10^{-3} of the Eddington accretion rate, the electrons have a thermal distribution throughout the accretion flow. For lower accretion rates, the electron distribution function is determined by the electron's source of heating, which is primarily adiabatic compression. Using the principle of adiabatic invariance, we show that an adiabatically compressed collisionless gas maintains a thermal distribution until the particle energies become relativistic. We derive a new, non-thermal, distribution function which arises for relativistic energies and provide analytic formulae for the synchrotron radiation from this distribution. Finally, we discuss its implications for the emission spectra from ADAFs.Comment: 29 pages (Latex), 3 Figures. Submitted to Ap

Spectrum of Optically Thin Advection Dominated Accretion Flow around a Black Hole: Application to Sgr A*

The global structure of optically thin advection dominated accretion flows which are composed of two-temperature plasma around black holes is calculated. We adopt the full set of basic equations including the advective energy transport in the energy equation for the electrons. The spectra emitted by the optically thin accretion flows are also investigated. The radiation mechanisms which are taken into accout are bremsstrahlung, synchrotron emission, and Comptonization. The calculation of the spectra and that of the structure of the accretion flows are made to be completely consistent by calculating the radiative cooling rate at each radius. As a result of the advection domination for the ions, the heat transport from the ions to the electrons becomes practically zero and the radiative cooling balances with the advective heating in the energy equation of the electrons. Following up on the successful work of Narayan et al. (1995), we applied our model to the spectrum of Sgr A*. We find that the spectrum of Sgr A* is explained by the optically thin advection dominated accretion flow around a black hole of the mass M_bh=10^6 M_sun. The parameter dependence of the spectrum and the structure of the accretion flows is also discussed.Comment: AAS LaTeX file; 26 pages; 12 ps figures; to be published in ApJ. PDF files are obtainable via following anonymous ftp. ftp://ftp.kusastro.kyoto-u.ac.jp/pub/manmoto/preprint/spec_sgrA.tar.g

Pair plasma relaxation time scales

By numerically solving the relativistic Boltzmann equations, we compute the time scale for relaxation to thermal equilibrium for an optically thick electron-positron plasma with baryon loading. We focus on the time scales of electromagnetic interactions. The collisional integrals are obtained directly from the corresponding QED matrix elements. Thermalization time scales are computed for a wide range of values of both the total energy density (over 10 orders of magnitude) and of the baryonic loading parameter (over 6 orders of magnitude). This also allows us to study such interesting limiting cases as the almost purely electron-positron plasma or electron-proton plasma as well as intermediate cases. These results appear to be important both for laboratory experiments aimed at generating optically thick pair plasmas as well as for astrophysical models in which electron-positron pair plasmas play a relevant role.Comment: Phys. Rev. E, in pres

Scaling Laws for Advection Dominated Flows: Applications to Low Luminosity Galactic Nuclei

We present analytical scaling laws for self-similar advection dominated flows. The spectra from these systems range from 10$^{8}$ - 10$^{20}$ Hz, and are determined by considering cooling of electrons through synchrotron, bremsstrahlung, and Compton processes. We show that the spectra can be quite accurately reproduced without detailed numerical calculations, and that there is a strong testable correlation between the radio and X-ray fluxes from these systems. We describe how different regions of the spectrum scale with the mass of the accreting black hole, $M$, the accretion rate of the gas, $\dot{M}$, and the equilibrium temperature of the electrons, $T_e$. We show that the universal radio spectral index of 1/3 observed in most elliptical galaxies (Slee et al. 1994) is a natural consequence of self-absorbed synchrotron radiation from these flows. We also give expressions for the total luminosity of these flows, and the critical accretion rate, $\dot{M}_{crit}$, above which the advection solutions cease to exist. We find that for most cases of interest the equilibrium electron temperature is fairly insensitive to $M$, $\dot{M}$, and parameters in the model. We apply these results to low luminosity black holes in galactic nuclei. We show that the problem posed by Fabian & Canizares (1988) of whether bright elliptical galaxies host dead quasars is resolved, as pointed out recently by Fabian & Rees (1995), by considering advection-dominated flows.Comment: 30 pages, 5 postscript files. Accepted to ApJ. Also available http://cfa-www.harvard.edu/~rohan/publications.htm

Engineering simulations for cancer systems biology

Computer simulation can be used to inform in vivo and in vitro experimentation, enabling rapid, low-cost hypothesis generation and directing experimental design in order to test those hypotheses. In this way, in silico models become a scientific instrument for investigation, and so should be developed to high standards, be carefully calibrated and their findings presented in such that they may be reproduced. Here, we outline a framework that supports developing simulations as scientific instruments, and we select cancer systems biology as an exemplar domain, with a particular focus on cellular signalling models. We consider the challenges of lack of data, incomplete knowledge and modelling in the context of a rapidly changing knowledge base. Our framework comprises a process to clearly separate scientific and engineering concerns in model and simulation development, and an argumentation approach to documenting models for rigorous way of recording assumptions and knowledge gaps. We propose interactive, dynamic visualisation tools to enable the biological community to interact with cellular signalling models directly for experimental design. There is a mismatch in scale between these cellular models and tissue structures that are affected by tumours, and bridging this gap requires substantial computational resource. We present concurrent programming as a technology to link scales without losing important details through model simplification. We discuss the value of combining this technology, interactive visualisation, argumentation and model separation to support development of multi-scale models that represent biologically plausible cells arranged in biologically plausible structures that model cell behaviour, interactions and response to therapeutic interventions