We present initial results from a new 440-ks Chandra HETG GTO observation of the canonical Seyfert 2 galaxy NGC 1068. The proximity of NGC 1068, together with Chandra's superb spatial and spectral resolution, allow an unprecedented view of its nucleus and circumnuclear NLR. We perform the first spatially resolved high-resolution X-ray spectroscopy of the `ionization cone' in any AGN, and use the sensitive line diagnostics offered by the HETG to measure the ionization state, density, and temperature at discrete points along the ionized NLR. We argue that the NLR takes the form of outflowing photoionized gas, rather than gas that has been collisionally ionized by the small-scale radio jet in NGC 1068. We investigate evidence for any velocity gradients in the outflow, and describe our next steps in modeling the spatially resolved spectra as a function of distance from the nucleus.

We present initial results from a new 440-ks Chandra HETG GTO observation of the canonical Seyfert 2 galaxy NGC 1068. The proximity of NGC 1068, together with Chandra's superb spatial and spectral resolution, allow an unprecedented view of its nucleus and circumnuclear NLR. We perform the first spatially resolved high-resolution X-ray spectroscopy of the `ionization cone' in any AGN, and use the sensitive line diagnostics offered by the HETG to measure the ionization state, density, and temperature at discrete points along the ionized NLR. We argue that the NLR takes the form of outflowing photoionized gas, rather than gas that has been collisionally ionized by the small-scale radio jet in NGC 1068. We investigate evidence for any velocity gradients in the outflow, and describe our next steps in modeling the spatially resolved spectra as a function of distance from the nucleus.

We present results from modeling of quasi-simultaneous broad band (radio through X-ray) observations of the galactic stellar black hole (BH) transient X-ray binary (XRB) systems XTE J1118+480 and GX 339-4 using an irradiated disc + compact jet model. In addition to quantifying the physical properties of the jet, we have developed a new irradiated disc model which also constrains the geometry and temperature of the outer accretion disc by assuming a disc heated by viscous energy release and X-ray irradiation from the inner regions. For the source XTE J1118+480, which has better spectral coverage of the two in optical and near-IR (OIR) wavelengths, we show that the entire broad band continuum can be well described by an outflow-dominated model + an irradiated disc. The best-fit radius of the outer edge of the disc is consistent with the Roche lobe geometry of the system, and the temperature of the outer edge of the accretion disc is similar to those found for other XRBs. Irradiation of the disc by the jet is found to be negligible for this source. For GX 339-4, the entire continuum is well described by the jet-dominated model only, with no disc component required. For the two XRBs, which have very different physical and orbital parameters and were in different accretion states during the observations, the sizes of the jet base are similar and both seem to prefer a high fraction of non-thermal electrons in the acceleration/shock region and a magnetically dominated plasma in the jet. These results, along with recent similar results from modeling other galactic XRBs and AGNs, may suggest an inherent unity in diversity in the geometric and radiative properties of compact jets from accreting black holes.

Astronomical data generally consists of 2 or more high-resolution axes, e.g., X,Y position on the sky or wavelength and position-along-one-axis (long-slit spectrometer). Analyzing these multi-dimension observations requires combining 3D source models (including velocity effects), instrument models, and multi-dimensional data comparison and fitting. A prototype of such a "Beyond-XSPEC" (Noble & Nowak, 2008) system is presented here using Chandra imag- ing and dispersed HETG grating data. Techniques used include: Monte Carlo event generation, chi-squared comparison, conjugate gradient fitting adapted to the Monte Carlo characteristics, and informative visualizations at each step. These simple baby steps of progress only scratch the surface of the computational potential that is available these days for astronomical analysis.

We describe two means by which XSTAR, a code which computes physical conditions and emission spectra of photoionized gases, has been parallelized. The first is pvm_xstar, a wrapper which can be used in place of the serial xstar2xspec script to foster concurrent execution of the XSTAR command line application on independent sets of parameters. The second is PModel, a plugin for the Interactive Spectral Interpretation System (ISIS) which allows arbitrary components of a broad range of astrophysical models to be distributed across processors during fitting and confidence limits calculations, by scientists with little training in parallel programming. Plugging the XSTAR family of analytic models into PModel enables multiple ionization states (e.g., of a complex absorber/emitter) to be computed simultaneously, alleviating the often prohibitive expense of the traditional serial approach. Initial performance results indicate that these methods substantially enlarge the problem space to which XSTAR may be applied within practical timeframes.

We present a quantitative comparison between software features of the defacto standard X-ray spectral analysis tool, XSPEC, and ISIS, the Interactive Spectral Interpretation System. Our emphasis is on customized analysis, with ISIS offered as a strong example of configurable software. While noting that XSPEC has been of immense value to astronomers, and that its scientific core is moderately extensible--most commonly via the inclusion of user contributed "local models"--we identify a series of limitations with its use beyond conventional spectral modeling. We argue that from the viewpoint of the astronomical user, the XSPEC internal structure presents a Black Box Problem, with many of its important features hidden from the top-level interface, thus discouraging user customization. Drawing from examples in custom modeling, numerical analysis, parallel computation, visualization, data management, and automated code generation, we show how a numerically scriptable, modular, and extensible analysis platform such as ISIS facilitates many forms of advanced astrophysical inquiry.

Motivated by the emergence of multicore architectures, and the reality that parallelism is rarely used for analysis in observational astronomy, we demonstrate how general users may employ tightly-coupled multiprocessors in scriptable research calculations while requiring no special knowledge of parallel programming. Our method rests on the observation that much of the appeal of high-level vectorized languages like IDL or MatLab stems from relatively simple internal loops over regular array structures, and that these loops are highly amenable to automatic parallelization with OpenMP. We discuss how ISIS, an open-source astrophysical analysis system embedding the S-Lang numerical language, was easily adapted to exploit this pattern. Drawing from a common astrophysical problem, model fitting, we present beneficial speedups for several machine and compiler configurations. These results complement our previous efforts with PVM, and together lead us to believe that ISIS is the only general purpose spectroscopy system in which such a range of parallelism -- from single processors on multiple machines to multiple processors on single machines -- has been demonstrated.

We introduce SLIRP, a module generator for the S-Lang numerical scripting language, with a focus on its vectorization capabilities. We demonstrate how both SLIRP and S-Lang were easily adapted to exploit the inherent parallelism of high-level mathematical languages with OpenMP, allowing general users to employ tightly-coupled multiprocessors in scriptable research calculations while requiring no special knowledge of parallel programming. Motivated by examples in the ISIS astrophysical modeling & analysis tool, performance figures are presented for several machine and compiler configurations, demonstrating beneficial speedups for real-world operations.

A review of the literature reveals that while parallel computing is sometimes employed by astronomers for custom, large-scale calculations, no package fosters the routine application of parallel methods to standard problems in astronomical data analysis. This paper describes our attempt to close that gap by wrapping the Parallel Virtual Machine (PVM) as a scriptable S-Lang module. Using PVM within ISIS, the Interactive Spectral Interpretation System, we've distributed a number of representive calculations over a network of 25+ CPUs to achieve dramatic reductions in execution times. We discuss how the approach applies to a wide class of modeling problems, outline our efforts to make it more transparent for common use, and note its growing importance in the context of the large, multi-wavelength datasets used in modern analysis.

In this paper we describe "VWhere," a S-Lang guilet which augments the computational power of "where" with the point-click ease of a Gtk-based visual interface. The result is a new mechanism for interacting with datasets, which unifies the constraint evaluation cycle implicit within observational analysis and provides a number of compelling advantages over traditional tool-based methods.