Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “galactic”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Crystallography of rare galactic honeycomb structure near supernova 1987a

Near supernova 1987a, the rare honeycomb structure of 20-30 galactic bubbles measures 30 x 90 light years. Its remarkable regularity in bubble size suggests a single-event origin which may correlate with the nearby supernova. To test the honeycomb's regularity in shape and size, the formalism of statistical crystallography is developed here for bubble sideness. The standard size-shape relations (Lewis's law, Desch's law, and Aboav-Weaire's law) govern area, perimeter and nearest neighbor shapes. Taken together, they predict a highly non-equilibrium structure for the galactic honeycomb which evolves as a bimodal shape distribution without dominant bubble perimeter energy.

Noever, David A.↗
NASA NTRS1994

The ultraviolet emission properties of five low-redshift active galactic nuclei at high signal-to-noise ratio and spectral resolution

We analyze the ultraviolet (UV) emission line and continuum properties of five low-redshift active galactic nuclei (four luminous quasars: PKS 0405-123, H1821 + 643, PG 0953 + 414, and 3C 273, and one bright Seyfert 1 galaxy: Mrk 205). The HST spectra have higher signal-to-noise ratios (typically approximately 60 per resolution element) and spectral resolution (R = 1300) than all previously published UV spectra used to study the emission characteristics of active galactic nuclei. We include in the analysis ground-based optical spectra covering H beta and the narrow (O III) lambda lambda 4959, 5007 doublet. New results are obtained and presented.

Laor, Ari↗
NASA NTRS1994

X-ray color analysis of the spectra of active galactic nuclei

The identification and detection of X-ray absorption and emission features depends on the resolution and the signal-to-noise ratio (S/N) of the observation, the understanding of the instrument response, and the Galactic line-of-sight absorption. Since many of the active galactic nucleus (AGN) data sets are limited in their S/N and full modeling of the physical conditions is rather complicated, we suggest a new analysis method based on 'X-ray colors.' The two sets of X-ray colors, defined for low (ROSAT Position Sensitive Proportional Counter (PSPC)) and medium (Broad Band X-Ray Telescope (BBXRT)) and ASCA Solid-State Imaging Spectrometers (SIS) resolution experiments, are used to separate regions of different physical conditions in a two-dimensional color-color plane. They are similar but superior to previous methods using the X-ray 'hardness ratio' in being able to reveal more of the physical properties of the source. We illustrate the use of such diagrams by studying a number of AGNs suspected of showing absorption features. A sample of 14 AGNs observed by the ROSAT PSPC is presented which includes several objects with suspected 'warm absorbers' along the line-of-sight to the nucleus, several others exhibiting intrinsic continuum variations, and a number of control objects thought to be featureless. Our new observations show, for the first time, the color variation as a function of time for three of the Seyfert 1 sources: NGC 4051, Mrk 335, and Mrk 766. The variations suggest that in two sources we are witnessing real changes in continuum shape, while one (NGC 4051) is consistent with having a warm absorber. Four of the objects observed by BBXRT are reanalyzed using our X-ray colors. Out of these, we discuss in detail the case of NGC 4151 and show that the color-color analysis agrees very well with previous, detailed spectral fitting methods. In particular, we confirm that the observed BBXRT observation of this source is not consistent with the warm absorber explanation and must be interpreted as a combination of a variable hard component and another component. We explain the limitations in measuring low-energy emission features in ASCA SIS observations and present color-color diagrams which can be used as a powerful diagnostic tool in the analysis of data collected by this experiment.

Netzer, Hagai↗
NASA NTRS1995

Elliptical accretion disks in active galactic nuclei

We present a calculation of the profiles of emission lines originating in a relativistic, eccentric disk, and show examples of the resulting model profiles. Our calculations are motivated by the fact that in about one-quarter of the double-peaked emission lines observed in radio-loud active galactic nuclei (and in the mildly active nucleus of NGC 1097), the red peak is stronger than the blue peak, which is contrary to the prediction of relativistic, circular disk models. Using the eccentric disk model we fit some of the observed profiles that cannot be fitted with a circular disk model. We propose two possible scenarios for the formation of an eccentric disk in an active galactic nucleus: (a) tidal perturbation of the disk around a supermassive black hole by a smaller binary companion, and (b) formation of an elliptical disk from the debris resulting from the tidal disruption of a star by the central black hole. In the former case we show that the eccentricity can be long-lived because of the presence of the binary companion. In the latter case, although the inner parts of the disk may circularize quickly, we estimate that the outer parts will maintain their eccentricity for times much longer than the local viscous time. We suggest that it may be possible to detect profile variability on much shorter timescales than those ranging from a decade to several centuries by comparing the evolution of the line profile with detailed model predictions. We argue that line-profile variability may also be the most promising discriminant among competing models for the origin of asymmetric, double-peaked emission lines.

Eracleous, Michael↗
NASA NTRS1994

Dynamics of massive black holes as a possible candidate of Galactic dark matter

If the dark halo of the Galaxy is comprised of massive black holes (MBHs), then those within approximately 1 kpc will spiral to the center, where they will interact with one another, forming binaries which contract, owing to further dynamical friction, and then possibly merge to become more massive objects by emission of gravitational radiation. If successive mergers would invariably lead, as has been proposed by various authors, to the formation of a very massive nucleus of 10(exp 8) solar mass, then the idea of MBHs as a dark matter candidate could be excluded on observational grounds, since the observed limit (or value) for a Galactic central black hole is approximately 10(exp 6.5) solar mass. But, if successive mergers are delayed or prevented by other processes, such as the gravitational slingshot or rocket effect of gravitational radiation, then a large mass accumulation will not occur. In order to resolve this issue, we perform detailed N-body simulations using a modfied Aarseth code to explore the dynamical behavior of the MBHs, and we find that for a 'best estimate' model of the Galaxy a runaway does not occur. The code treates the MBHs as subject to the primary gravitational forces of one another and to the smooth stellar distribution, as well as the secondary perturbations in their orbits due to another and to the smooth stellar distribution, as well as the secondary perturbations in their orbits due to dynamical friction and gravitational radiation. Instead of a runaway, three-body interactions between hard binaries and single MBHs eject massive objects before accumulation of more than a few units, so that typically the center will contain zero, one, or two MBHs. We study how the situation depends in detail on the mass per MBH, the rotation of the halo, the mass distribution within the Galaxy, and other parameters. A runaway will most sensitively depend on the ratio of initial (spheroid/halo) central mass densities and secondarily on the typical values for the mass per MBH, with the rough dividing line, using Galactic parameters, being M(sub BH) less than or = 10(exp 6.5) solar mass. Using parameters from Lacey & Ostriker (1985) and our most accurate model for Galaxy, no runaway occurs.

Xu, Guohong↗
NASA NTRS1994

Highly ionized gas in the Galactic halo

We reexamine the values of electron density n(sub e) and gas pressure P/k in the interstellar medium (ISM) of the Galactic halo, as inferred from C IV emission and absorption lines and using current C IV atomic data. In a homogeneous model with 4.7 less than or equal to log T less than or equal to 5.3, the data are consistent with 0.01 less than or equal to n(sub e) less than or equal to 0.02/cu cm and 2200 less than or equal to P/k less than or equal to 3700/cu cm K, a factor of 2-3 higher than advocated by Martin & Bowyer (1990) and comparable to the thermal pressure in the disk. If some of the C IV absorption arises from nonemitting, photoionized gas, then the inferred density and pressure will increase accordingly. The volume filling factor for homogeneous models ranges from 0.5% to 5%. Because of the constraints arising from filling factor and radiated power, most of the C IV must arise from gas near the peak of the cooling curve, at log t less than or equal to 5.6. We relate both emission-line and absorption-line observations to recent models in which turbulent mixing layers and isobarically cooling supernova remnants (SNRs) provide significant amounts of halo gas at approximately 10(exp 5.3) K and process 20-40 solar mass/yr with a power of approximately 10(exp 41) ergs/sec. Since the observed C IV and N V absorption scale heights have been reported to differ, at 4.9 kpc and 1.6 kpc, respectively, we examine inhomogeneous models with different exponential scale heights of T, P, and SN energy input. The ISM may change its character with distance above the Galactic plane, as superbubbles and mixing layers dominate over isolated SNRs as the source of the C IV. For appropiate scale heights, the midplane pressure is twice the homogeneous values quoted above. The O IV lambda 1034 diffuse emission line, which can be used as a temperature diagnostic of the hot gas, is predicted to be comparable in strength to that of C IV lambda 1549 (approximately 6000 photons/cu cm/s/sr. The ions C IV, N V, and O VI are predicted to show a decreasing trend in vertical extent, with the O IV scale height perhaps as low as approximately 1 kpc.

Shull, J. Michael↗
NASA NTRS1995

Microlensing by the galactic bar

We compute the optical depth and duration distribution of microlensing events towrd Baade's window in a model composed of a Galactic disk and a bar. The bar model is a self-consistent dynamical model built out of individual orbits that has been populated to be consistent with the COBE maps of the Galaxy and kinematic observations of the Galactic bulge. We find that most of the lenses are in the bulge with a line-of-sight distance 6.25 kpc (adopting R(sub 0) = 8 kpc). The microlensing optical depth of a 2 x 10(exp 10) solar mass bar plus a truncated disk is (2.2 +/- 0.45) x 10(exp -6), consistent with the large optical depth (3.2 +/- 1.2) x 10(exp -6) found by Udalski et al. (1994). This model optical depth is enhanced over the predictions of axisymmetric models by Kiraga & Paczynski (1994) by slightly more than a factor of 2, since the bar is elongated along the line of sight. The large Einstein radius and small transverse velocity dispersion also predict a longer event duration in the self-consistent bar model than in the Kiraga-Paczynski model. The event rate and duration distribution also depend on the lower mass cutoff of the lens mass function. With a 0.1 solar mass cutoff, five to seven events (depending on the contribution of disk lenses) with a logarithmic mean duration of 20 days are expected for the Optical Gravitational Lensing Experiment (OGLE) according to our model, while Udalski et al. (1994) observed nine events with durations from 8 to 62 days. On the other hand, if most of the lenses are brown dwarfs, our model predicts too many short-duration events. A Kolmogorov-Smirnov test finds only 7% probability for the model with 0.01 solar mass cutoff to be consistent with current data.

Zhao, Hongsheng↗
NASA NTRS1995

Galactic cosmic-ray mediation of a spherical solar wind flow. 1: The steady state cold gas hydrodynamical approximation

Realistic models of the outer heliosphere should consider that the interstellar cosmic-ray pressure becomes comparable to pressures in the solar wind at distances more than 100 AU from the Sun. The cosmic-ray pressure dynamically affects solar wind flow through deceleration. This effect, which occurs over a scale length of the order of the effective diffusion length at large radial distances, has important implications for cosmic-ray modulation and acceleration. As a first step toward solution of this nonlinear problem, a steady state numerical model was developed for a relatively cold spherical solar wind flow which encounters the confining isotropic pressure of the surrounding Galactic medium. This pressure is assumed to be dominated by energetic particles (Galactic cosmic rays). The system of equations, which are solved self-consistently, includes the relevant hydrodynamical equations for the solar wind flow and the spherical cosmic-ray transport equation. To avoid the closure parameter problem of the two-fluid model, the latter equation is solved for the energy-dependent cosmic-ray distribution function.

Le Roux, J. A.↗
NASA NTRS1995

Discovery and imaging of a Galactic cirrus cloud with the far ultraviolet space telescope

We present new far-ultraviolet (1400-1800 A) data concerning a Galactic cirrus cloud G251.2+73.3 near the north Galactic pole obtained with the space-borne imaging telescope FAUST (Far Ultraviolet Space Telescope). We obtain a good correlation between the far-ultraviolet (FUV) and IRAS 100 micrometers surface brightnesses, their relation being I(sub FUV) = (128 +/- 3) I(sub 100 micrometers) - (264 +/- 9), where the I(sub FUV) flux is given in units of photon/s/sq cm/A/sr and I(sub 100 micrometers) in MJy/sr. Using uvbyH-beta photometry, we get a distance of 120 pc and a visual extinction in the center of the cloud of 0.39 mag corresponding to an extinction of 1.0 mag at 1565 A. We have performed a multiple scattering calculation for the scattered light using the Monte Carlo method. These calculations provide restrictions on the FUV scattering properties of the interstellar dust.

Haikala, Lauri K.↗
NASA NTRS1995

Search for Compton-backscattered annihilation radiation from the galactic center with the OSSE

An emission feature near 170 keV, interpreted as Compton-backscattered 511 keV positron-annihilation radiation, has been reported twice by balloon-borne germanium spectrometers from within approximately 15 deg of the Galactic center (Leventhal, MacCallum, & Stang 1978; Smith et al. 1993). Upper limits on this feature set by HEAO 3 (Mahoney, Ling, & Wheaton 1993) and other instruments indicate that it must be transient. We have searched data from the Oriented Scintillation Spectrometer Experiment (OSSE) on the Compton Gamma Ray Observatory (CGRO) for this feature, using daily spectral accumulations from all pointings near the Galactic center up to 1993 August, and covering most of the region viewed by the balloon instruments. We find no evidence for backscatter emission. Under the hypothesis that the source is 1E 1740.7-2942, the OSSE data set (186 days) disagrees with the balloon measurements with 99.3% confidence. The average daily 3 sigma OSSE upper limit on bakscatter flux from 1E 1740.7-2942 is 6.8 x 10(exp -4) photons/sq cm/s, compared to the 1.3 x 10(exp -3) photons/sq cm reported by the balloon observations. We also saw no evidence in 186 days for linelike emission from the point source EXS 1737.9-2952 recently discovered by Grindlay, Covault, & Manandhar (1993). This source exhibited bright emission from 83-111 keV, which has been interpreted as doubly backscattered 511 keV radiation. The average daily 3 sigma upper limit from OSSE for this line is 9.8 x 10(exp -4) photons/sq cm/s, or approximately 8% of the reported flux.

Smith, D. M.↗
NASA NTRS1995

An HST study of galactic inerstellar zinc and chromium

We present a survey of interstellar Zn II and Cr II absorption extracted from the Hubble Space Telescope Goddard High Resolution Spectrograph (HST GHRS) data archive. We find clear evidence for an enhanced depletion of Zn from the gas phase with increasing fractional abundance of molecular hydrogen f(H2). Our lower limit to the Galactic interstellar metallicity is approximately 65% of the solar value as determined by the measured Zn abundances in the lowest f(H2) sightlines, (N(Zn)/N(H(sup 0)(sub tot)) = -0.19 +/- 0.04. The correspondingly high depletion of Cr with respect to solar (N(Cr/N(H(sup 0)(sub tot)) = -1.44 +/- 0.26 indicates that there are significant amounts of dust present in these lines of sight. The Galactic abundances of Zn and Cr in the ISM provide a fundamental reference point which is used to understand the metal enrichment and dust formation history of damped Lyman alpha QSO absorption-line systems, normally believed to arise from intervening precursors to modern disk galaxies. Although the spread in Zn abundances is large for both the local ISM and in damped Lyman alpha systems, we still find a substantial difference (factor of 4-10) in metallicity between the two sets. This survey and future observations of more distant objects which probe the full extent of the Milky Way halo provide a more complete picture of the enrichment and depletion characteristics of present-day galaxies.

Roth, Katherine C.↗
NASA NTRS1996

Shielding against galactic cosmic rays

Ions of galactic origin are modified but not attenuated by the presence of shielding materials. Indeed, the number of particles and the absorbed energy behind most shield materials increases as a function of shield thickness. The modification of the galactic cosmic ray composition upon interaction with shielding is the only effective means of providing astronaut protection. This modification is intimately conntected with the shield transport porperties and is a strong function of shield composition. The systematic behavior of the shield properites in terms of microscopic energy absorption events will be discussed. The shield effectiveness is examined with respect to convectional protection practice and in terms of a biological endpoint: the efficiency for reduction of the probability of transformation of shielded C3H1OT1/2 mouse cells. The relative advantage of developing new shielding technologies is discussed in terms of a shield performance as related to biological effect and the resulting uncertainty in estimating astronaut risk.

Schimmerling, W.↗
NASA NTRS1995

KAO and AAT observations of the Galactic Center filaments

We have used the Kuiper Airborne Observatory (KAO) and the Anglo-Australian Telescope (AAT) to investigate the nature of the filamentary radio emission from the Galactic center region. KAO observations of the FIR line and continuum emission from the radio peak G0.095+0.012 and the E2 thermal radio filament northeast of the Galactic center can be produced by numerous nearby stars with T(sub eff) approx. 35,000 K; these can account for both the FIR luminosity and the excitation of the gas. Much of the FIR continuum and most of the strong (Si II) (34.8 micron) line emission are probably produced in the ionized gas of the filament. The FIR (O III) 52 and 88 micron lines imply an electron density of a few hundred; when compared with the radio emission measure, this implies the filament is roughly tubular or somewhat flattened in the plane of the sky. The (O III) and (S III) lines show higher excitation associated with the filament, and suggest that exciting stars may be located within the filaments and/or southeast of the E2 filament. AAT observations in the near infrared (NIR) in fact reveal a nearby cluster of hot stars southeast of the E2 filament. Additional hot stars, not identifiable from their NIR spectra, are likely to be present. These stars and those in the cluster can plausibly produce the observed radio and FIR emission in the region. The morphology of the filament is not explained by existing information however.

Erickson, Edwin F.↗
NASA NTRS1997

Spectral and Temporal Properties of Galactic Black Hole Systems

Kusunose, Mineshige & Yamada (1996; hereafter KMY) extended the model of Kusunose & Mineshige (1995) to the Galactic black hole candidates by considering nonthermal electron injection with gamma(EQ\0(,\s\up2(less than),\s\do-l(_))) 10. The effects of pair escape and advection on the disk structure and general relativistic effects on the emission spectrum were also examined. They found that the energy spectral index (alpha)(sub x) of the power law X-rays is about-0.8 and-2.0 when 1(sub soft)/1 = 0.2 and 2, respectively, where 1(sub soft)/1 is the ratio of the compactness of the injected soft photons to that of the gravitational energy. The power law index was found to be nearly independent of the mass accretion which is consistent with the observed luminosity independence. The model with small 1(sub soft)/1 (less than 1) shows promise for explaining the low state observed in Galactic black hole candidates. Model fits were provided for GX339-4 and Cyg X- 1 data from COMPTEL and OSSE on the Compton Gamma Ray Observatory. The difference in emission spectra between thermal disks and the model of KMY appears only in the energy range greater than 100 keV. Li, Kusunose and Liang (1996) studied stochastic particle acceleration to produce nonthermal particle distributions which then were used in the model of Kusunose & Mineshige (1995) to model the spectrum above 1 Mev from GBHC's. Under certain conditions, stochastic electron acceleration overcame Coulomb and Compton losses resulting in a suprathermal electron population. Good fits were obtained by COMPTEL and OSSE observations of Cyg X-1 and GRO J0422+23. Kusunose & Mineshige (1996a) examined the role of electron-positron pairs in advection-dominated disks. They found that the results for advection-dominated disks without pairs are not qualitatively changed by including pairs. Summaries of work sponsored by this grant are given in Wheeler, Kim, Moscoso, Kusunose & Mineshige (1996) and Kusunose (1996) Work was also done on developing a model for an e(+-) pair wind from the inner disk region of a black hole. The model consists of three zones: a pair production/annihilation zone at the base of the wind, a pair annihilation zone slightly further out from the disk and a wind zone extending to infinity where no annihilation occurs. The model assumes an input X-ray / gamma-ray power-law spectrum as a function of photon energy and angular distribution. Pairs are created in the pair production/annihilation zone via photon-photon, photonparticle and particle-particle collisions. The bulk velocity of the pairs is obtained from the momentum component perpendicular to the disk taking into account the radiation pressure acceleration of the pairs. Energy balance in the pair production/annihilation zone is calculated by assuming that the momentum component parallel to the disk is thermalized and by taking into account heating/cooling via anisotropic thermal Compton scattering along with the cooling processes of bremsstrahlung and synchrotron radiation. The pair density is also calculated. The outer boundary of the pair production/annihilation zone is where the optical depth for photon-photon pair produciton is unity. Beyond this point only pair annihilation will occur in the pair annihilation zone. When the pair density becomes small, pair annihilation will be negligible and the pairs will flow freely to infinity. In this model we have found that the X-ray power-law spectral index, alpha(sub x) is the primary parameter which determines the density, temperature and velocity of the pair production/annihilation zone (Moscoso, Kusunose & Wheeler 1996).

Wheeler, J. Craig↗
NASA NTRS1998

Variability Analysis and the Structure of Active Galactic Nuclei

This five-year Long-Term Space Astrophysics grant provided the support for several major steps in advancing our knowledge of the internal structure of active galactic nuclei. The single largest portion of this program had to do with the development and application of techniques for "reverberation mapping", the use of spectral monitoring of several different bands related by radiation reprocessing to infer the internal geometry of sources. Major steps were taken in this regard, particularly in establishing the distribution in radius of emission line material, and in relating the apparent reprocessing of continuum bands to the underlying structure of the accretion disk. Another major effort built directly upon these results. Once the case for continuum reprocessing was made by the monitoring, it next behooved us to understand the spectral output of AGN as a result of this reprocessing. As a result, our view of continuum production in AGN is now much better focussed on the key problems. A third focus of effort had to do with the nature of X-ray variability in AGN, and what it can tell us about the dynamics of extremely hot material in the immediate outskirts of the supermassive black holes that form the central engines of active galactic nuclei. In addition to these primary efforts, this grant also supported many other, smaller projects. Several of these were demonstrations of how the material spewed out of AGN in relativistic.ets generate the radiation by which we observe them. J Finally, the portion of this study that had to do with continuum production by accretion disks in AGN led naturally to several papers in which new developments were presented having to do with "advection-dominated accretion disks", those disks in which accretion appears to proceed at a substantial rate, but in which radiation processes are weak.

Krolik, Julian H.↗
NASA NTRS1998

A Three-Dimensional Analysis of the Galactic Gamma-Ray Emission Resulting from Cosmic-Ray Interaction with the Interstellar Gas and Radiation Fields

This final report outlines the progress of the contractor's support for the analysis of data under ADP (NRA 96-ADP-08; Proposal No. 167-96 adp). The primary task object was to construct a 3-D model for the distribution of high-energy (20 MeV-30 GeV) gamma-ray emission in the Galactic disk. Under this task the contractor was to utilize data from the EGRET instrument on the Compton Gamma-Ray Observatory, HI and CO surveys, radio-continuum surveys at 408 MHz, 1420 MHz, 5 GHz, and 19 GHz, the COBE Diffuse Infrared Background Experiment (DIRBE) all-sky maps from I to 240 um, and ground-based B,V, J, H and K photometry. The respective contributions to the high-latitude gamma-ray emission from cosmic ray-matter interactions, inverse Compton scattering, and extra-galactic emission were to be determined.

Sodraski, Thomas J.↗
NASA NTRS1999

Interpreting Broad Double-Peaked Emission Lines in Active Galactic Nuclei

The principal objectives of this project were to probe the inner regions of active galactic nuclei and to test general relativity in the strong-field limit. The approach takes advantage of broad atomic line emission observed from material deep in the potential well of an active galactic nucleus which contains key information as to the physics of the system. Line profiles in a wide range of wavebands from optical to X-ray have provided compelling evidence of the existence of a relativistic accretion disk around a supermassive black hole in a number of galaxies. The simplest model posits a geometrically thin disk in Keplerian orbit, with general relativistic effects in evidence. This model is the point of departure for the proposed work. We developed a high-performance numerical code to calculate photon trajectories in a Schwarzschild or Kerr metric and implemented it on parallel supercomputers. This code includes a general purpose ray tracer that calculates line profiles, light curves, and other observable quantities for a wide variety of emitter configurations. The versatility comes from the fact that the ray tracing algorithm does not depend on any symmetries regarding emitter locations. The speed comes from parallel implementation which enables us to sample hitherto unattainable volumes of disk model parameter space. During the period 1 March 1997 through 28 February 1998, two papers, supported in whole or in part by this grant, were published in refereed journals. They are reproduced in their entirety in the next two sections of this report.

Halpern, Jules↗
NASA NTRS1999

Research on the Nature of Accretion in Active Galactic Nuclei

he purpose of this grant was to support theoretical research on the nature of accretion in active galactic nuclei. In the brief time of the award (one year), four papers that appeared in refereed journals were written, as well as two invited reviews in conference proceedings These papers significantly advanced our understanding of the structure of the most important parts of bright accretion disks around accreting black holes, such as active galactic nuclei.

Krolik, J. H.↗