Radio structure of extragalactic X-ray sources

The radio properties ofUhuru X-ray sources with fairly certain extragalactic identifications are described briefly. Radio to X-ray flux ratios are low for rich clusters of galaxies and high for double radio sources. There is some evidence from the Abell 426 (Perseus) and Abell 1367 clusters that a radio galaxy in a rich cluster may be the centre of extended X-ray emission. Nuclei of galaxies have an enormous range in X-ray luminosity; the known range is from 10³⁰ W for our galaxy to 3×10³⁸ W for 3C 273. Unidentified X-ray sources at high galactic latitudes may include new classes of objects with very low radio to X-ray flux ratios or hard X-ray emission.     

Pulsar dispersion measures and the mean and mean-squared interstellar electron density

Theoretical distributions of pulsar dispersion measures (times sinb ^II) are computed for various assumed pulsar spatial distributions above (and below) the galactic plane, assuming a distribution for the ionized gas. The statistics on the twelve high-latitude pulsars lead to the conclusion that the pulsar distribution inz is at least as broad as the distribution of the ionized gas. The value derived for the local mean interstellar electron density in the galactic plane is 0.12 cm^–3 and is interpreted to be due to a uniform ionized intercloud medium.Interstellar absorption of radio waves at low frequencies and cosmic X-rays at low energies are considered with regard to irregularities in the distribution of ionized gas. It is shown that if the obervations are made with a wide angle receiver the effective absorption optical deph is –²/2 where is the mean value and is the dispersion in . This relation assumes is much larger than . Analysis of recent low-frequency radio measurements from a satellite-borne receiver, however, leads to the conclusion that effects of irregulatirities are large.     

Primeval galaxies: Predicted luminosities

Absorption by gas and dust in circumstellar Hii regions within primeval galaxies could seriously depress the far-ultraviolet continuum radiation emitted by primeval galaxies. This effect might account for the failure of Partridge (1974) and Davis and Wilkinson (1974) to detect the redshifted radiation from primeval galaxies at optical and near-infrared wavelengths. A primeval galaxy becomes very bright only during the final stages of contraction. Provided that dust can form by the time the primeval galaxy reaches peak luminosity, a significant fraction of the stellar far-ultraviolet radiation is converted into far-infrared. Thus an appropriate spectral region to search for the redshifted integrated background from primeval galaxies lies between 350 , where the 2.7 K microwave background radiation becomes important, and 150 , where other extragalactic discrete sources, such as nearby galactic nuclei, may contribute. The expected IR flux is calculated with Kaufman's (1975) model for the star formation rate in the contracting galaxy. Letz _p be the redshift andT _g the grain temperature when the primeval galaxy becomes very bright. Unlessz _p10 orT _g is fairly high, the intensity of the far-infrared radiation from primeval galaxies would be dominated by the high frequency tail of the 2.7 K microwave background. On the other hand, if dust is unimportant, we determine the spectral energy distribution of a primeval galaxy emitted in the range 912 Å to 2050 Å; we find that the luminosities are not very sensitive to the dependence of effective temperatures on metal abundance.     

Probing the Obscuring Medium Around Active Nuclei Using Masers: The Case of 3C 403

We report the first detection of a water megamaser in a radio-loud galaxy, 3C 403, and present a follow-up study using the VLA. 3C 403 has been observed as a part of a small sample of FR II galaxies with evidence of nuclear obscuration. The isotropic luminosity of the maser is 1200 L. With a recessional velocity of cz 17680 km s^–1 it is the most distant water maser so far reported. The line arises from the densest (> 10⁸ cm^–3) interstellar gas component ever observed in a radio-loud galaxy. Two spectral features are identified, likely bracketing the systemic velocity of the galaxy. Our interferometric data clearly indicate that these arise from a location within 0.1 (110 pc) from the active galactic nucleus. We conclude that the maser spots are most likely associated with the tangentially seen parts of a nuclear accretion disk, while an association with dense warm gas interacting with the radio jets cannot yet be ruled out entirely.     

Optically thick accretion onto black holes

Spherically symmetric, steady-state, optically thick accretion onto a nonrotating black hole with the mass of is studied. The gas accreting onto the black hole is assumed to be a fully ionized hydrogen plasma withn ₀=10⁸ cm^–3 andT ₀=10⁴ K far from the black hole, and a new approximate expression for the Eddington factor is introduced. The luminosity is estimated to beL=1.875×10³³ erg s^–1, which primarily arises from the optical surface (1) ofT10⁴ K. The accretion flow is characterized by 1 and (v/c)10. In the optically thin region, the flow remains isothermal, and the increase of temperature occurs at 1. The radiative equilibrium is strictly realized at (v/c)10.     

Quasi-stellar objects as rotating magnetic superstars

In this paper the magnetic superstar model is used to discuss QSO luminosity and density evolution. Our main hypotheses are that (i) mass loss from old stars in massive galaxies cools and then falls into the centre to form a nuclear disc (Bailey, 1980); and (ii) magnetic superstars in galactic nuclei condense out of gaseous material at the centre of a supermassive-magnetised disc (Kundt, 1979). On this generalised model we find that the non-thermal (synchrotron) optical luminosity scales asL _opt L ³ t ^–7/3, whereL is the total blue luminosity of old stars in the galaxy and t is cosmic time. In addition we show that QSO co-moving density follows the lawD(t)exp-(t/t _Evol)^16/15 with an evolution timescalet _Evol = 1.95 × 10⁹ yr. The model as a whole is in good agreement with observations.     

Electromagnetic solutions of higher-dimensional gravity

The variation of spectral index _1.4 ⁵ with luminosity (P _1.4) was investigated for Fanaroff and Riley type II galaxies and also for type I and type I/II galaxies. To reduce the effect of the redshift dependence of luminosity, samples which did not have widely different median redshifts were used and the data was binned into redshift ranges.By a comparison of the median spectral indices in different redshift bins for FRII galaxies, no dependence of spectral index was formed on redshift. However, an increase in spectral index with increasing luminosity was found. The results for FRI and FRI/II galaxies were similar to those for the FRII galaxies.     

The possible origin of the spiral-arm instability

Physical arguments suggest the spiral arms may be manifestations of the galaxy not being in dynamical equilibrium — in the sense that the kinetic energy of tis stars and gas is less relative to its binding energy than that dictated by the virial theorem. Without constant cooling of the galactic disk (i.e., a progressive increase in the binding energy of the galaxy) such a departure from dynamical equilibrium would be corrected and the spiral arms destroyed in about 10⁹ yr due to an increase in the velocity dispersion of the stars in the disk resulting from their interacting with the spiral arms. The rate of cooling required to maintain the spiral arms, about 6×10⁴ L , may be provided by mass loss from stars in the disk population. The cooling arises from the average scale-heights and velocities of these stars being larger than that of the gas in the disk, so that there is a net loss of kinetic energy and an increase in the binding energy of the galaxy due to the ejected gas settling down to a lower terminal velocity and scale-height in the galactic disk.     

On the interpretation of Fraunhofer line Doppler shifts at supergranule boundaries

We have used the SPO tower telescope and echelle spectrograph to study differences in the profiles of three Fei lines, between magnetic network and cells. Ca K slit-jaw pictures were used to identify the network and cell areas, and mean network and cell profiles were computed from digitized spectra for the g = 0 lines 4065, 5434, and the g = 1.5 line 5233. The profile bisectors show that the wings of all three lines are red-shifted in the network by between 75–200 m s^–1 relative to the cell profiles. But the redshift decreases in the line core and becomes less than the standard error of 20 m s^–1 near the line core minimum. This disappearance of the redshift at the cores of all 3 lines formed over the height range 250–500 km above _0.5 = 1, argues against a steady downflow at supergranule boundaries. We show that such red-shifted wings and a relatively unshifted core can result if granular convection is suppressed near the network flux tubes, without implying any downflow in the vicinity of these flux tubes. Our results also indicate that searches for large-scale convective velocity patterns should measure shifts of the line core, rather than the line wings which appear to be very sensitive to inhomogeneities in granule structure.Visiting Astronomers, Sacramento Peak Observatory.     

The distribution of galaxy pair redshifts

High signal to noise neutral hydrogen observations of a complete sample of 132 galaxy pairs show a velocity difference distribution which decreases monotonically from zero. There is no strong indication of redshift periodicity for the entire sample and no indication at all for a subset of 79isolated galaxy pairs. In addition, when redshifts are corrected for the solar motion around the Galactic Center there is no indication of a redshift periodicity of 37.6 km s^–1 in V for galaxy pairs (as suggested by Guthrie and Napier 1996).     

Planetary accretions in jet streams

The problem of planetary accretion in a jet stream is studied using the model developed by Alfvén and Arrhenius. We find that there are basically three types of planetary accretion: namely, fast process _c < _i , slow process _c ~ _i and delayed process _c > _i where _c is the characteristic time of the occurrence of catastrophic accretion and _i the time-scale of mass injection to the planetary system (3×10⁸ yr). These different time scales of accretion are found to be closely related to the primordial thermal profiles and equatorial inclinations of the planets. Finally, Venus' retrograde rotational spin is shown to be a possible result of accretion process in a jet stream.     

Linear hydrodynamical equations coupled with radiative transfer in a non-isothermal atmosphere

A method coupling the hydrodynamical equations and radiative transfer in a realistic solar model atmosphere is described. The influence of the temperature gradient of the model and the radiative dissipation is pointed out.The effect of the large temperature gradient is important in the layers where the optical depth ₅₀₀₀ is greater than 0.5; the ratio between the amplitude of the temperature and the velocity fluctuations decreases with the altitude by a factor 2 between = 1 and = 0.5 and in the case of the acoustic waves, the phase shift between these fluctuations is small.The radiative energy loss in the thick layers ( ₅₀₀₀ = 1) leads to a decrease of the vertical phase velocity of the waves and to a damping of their amplitudes in the layers of intermediate optical depth (10^-2 < ₅₀₀₀ < 0.5). The effect of the dissipation is negligible in the thin layers (₅₀₀₀ < 10^-2).     

A novel methodology for radiative transfer in a planetary atmosphere

A novel methodology for evaluating the field of anisotropically scattered radiation within a homogeneous slab atmosphere of arbitrary optical thickness is provided. It departs from the traditional radiative transfer approach in first considering that the atmosphere is illuminated by an isotropic light source. From the solution of this problem, it subsequently proceeds to that for the more conventional case of monodirectional illumination. The azimuthal dependence of the field is separated in the usual manner by an harmonic expansion, leaving a problem in four dimensions (=optical depth, ₀=thickness, , =directions of incidence and scattering) which, as is well known, is numerically extremely inconvenient. Two auxiliary radiative transfer formulations of increasing dimensionality are considered: (i) a transfer equation for the newly introduced functionb ^m(,,₀) with Sobolev's function ^m(,₀) playing the role of a source-function. Because the incident direction does not intervene, ^m is simply expressed as a single integral term involvingb ^m. For bottom illumination, an analogous equation holds for the other new functionh ^m(,,₀). However, simple reciprocity relations link the two functions so that it is only necessary to considerb ^m; (ii) a transfer equation for the other new functiona ^m(,,,₀) with a source-function provided by Sobolev's functionD ^m(,,₀). For bottom illumination, another functionf ^m(,,,₀) is introduced; by a similar argument using reciprocity relations,f ^m is reduced toa ^m rendering necessary only the consideration ofa ^m. However, a fundamental decomposition formula is obtained which shows thata ^m is expressible algebraically in terms of functions of a single angular variable. The functions ^m andD ^m are shown to be the values in the horizontal plane ofb ^m anda ^m, respectively. The other auxiliary functionsX ^m andY ^m are also expressed algebraically in terms ofb ^m. These results enable one to proceed to the final step of evaluating the radiation field for monodirectional illumination. The above reductions toalgebraic relations involving only the functionb ^m appear to be more advantageous than Sobolev's (1972) recent approach; they also circumvent some basic numerical difficulties in it. We believe the present approach may likewise prove to be superior to most (if not all) other methods of solution known heretofore.This paper presents the results of one phase of research carried out at the Jet Propulsion Laboratory under Contract No. NAS-7-100 sponsored by the National Aeronautics and Space Administration.     

Massive accretion disks

Recent high resolution near infrared (HST-NICMOS) and mm-interferometric imaging have revealed dense gas and dust accretion disks in nearby ultra-luminous galactic nuclei. In the best studied ultraluminousIR galaxy, Arp 220, the 2m imaging shows dust disks in both of the merging galactic nuclei and mm-CO line imaging indicates molecular gasmasses 10⁹M for each disk. The two gas disks in Arp 220 are counterrotating and their dynamical masses are 2×10⁹ M , that is, only slightly largerthan the gas masses. These disks have radii 100 pc and thickness 10-50 pc. The high brightness temperatures of the CO lines indicatethat the gas in the disks has area filling factors 25-50% and mean densitiesof 10⁴ cm^-3. Within these nuclear disks, the rate of massive star formation is undoubtedly prodigious and, given the high viscosity of the gas, there will also be high radial accretion rates, perhaps 10 M yr ^-1. If this inflow persists to very small radii, it is enough to feed even the highest luminosity AGNs.     

Quasars in Regions of Rich Clusters of Galaxies

Certain characteristics of quasars observed in the regions of three rich clusters of galaxies are investigated. It is shown that the luminosity functions constructed for samples of quasars from the regions of the Virgo and Fornax clusters are similar to each other and show a very sharp maximum at a stellar magnitude of 18^m.0-18^m.3. The maximum of the luminosity function for quasars from the region of the Coma cluster (A 1656) is shifted toward fainter magnitudes and lies at 21^m. This effect is explained by the affiliation of the majority of the quasars to the respective clusters. It is shown that the absolute stellar magnitude of these quasars is –13^m, which corresponds to the luminosity of dwarf galaxies of low luminosity. It is suggested that the local quasars are formed in two ways: by ejection from galactic nuclei and in quasar associations without a parent galaxy, directly from protostellar matter.     

Physical Properties of Dust in the Martian Atmosphere: Analysis of Contradictions and Possible Ways of Their Resolution

Atmospheric aerosols play an important role in forming the Martian climate. However, the basic physical properties of the Martian aerosols are still poorly known; there are many contradictions in their estimates. We present an analytical overview of the published results and potentialities of various methods. We consider mineral dust. Zonally averaged data obtained from mapping IR instruments (TES and IRTM) give the optical thickness of mineral aerosols ₉ = 0.05–0.1 in the 9-m band for quite atmospheric conditions. There is a problem of comparing these estimates with those obtained in the visible spectral range. We suggest that the commonly used ratio _vis/₉ >2 depends on the interpretation and it may actually be smaller. The ratio _vis/₉ 1 is in better agreement with the IRIS data (materials like montmorillonite). If we assume that _vis/₉ = 1 and take into account the nonspherical particle shape, then the interpretation of ground-based integrated polarimetric observations ( < 0.04) can be reconciled with IR measurements from the orbit. However, for thin layers, the sensitivity of both methods to the optical thickness is poorly understood: on the one hand, polarimetry depends on the cloud cover and, on the other hand, the interpretation of IR measurements requires that the atmospheric temperature profile and the surface temperature and emissivity be precisely known. For quite atmospheric conditions, the local optical-thickness estimates obtained by the Bouguer–Lambert–Beer method and from the sky brightness measured from Viking 1 and 2 and Mars Pathfinder landers are much larger: = 0.3–0.6. Estimates of the contrasts in images from theVikingorbiters yield the same values. Thus, there is still a factor of 3 to 10 difference between different groups of optical-thickness estimates for the quiet atmosphere. This difference is probably explained by the contribution of condensation clouds and/or by local/time variations.     

30 Doradus as the active centre of the Large Magellanic Cloud

Evidence is presented for the hypothesis that the supergiantHii-complex 30 Doradus (NGC 2070) is the mildly active galactic nucleus of the Large Magellanic Cloud. For this purpose the general properties of galactic nuclei and the characteristics of active nuclei are reviewed (Section 2). Examination of 30 Doradus shows that it plasy the same exceptional role among allHii-regions of the LMC as Sgr A among those of our Galaxy, and has all the properties of a galactic nucleus (luminosity, emission spectrum, IR source, semistellar central object R 136, symmetry centre of an starting point of spiral structure). Evidence for the activity is given by the peculiar filamentary structure (Figure 1), the young spiral filaments superposed on old, broad and smooth near-circular arms (Figure 2), the splitting of the [Oiii] 5007 profile in two components corresponding to an expansion velocity of 50 km s^–1, and the strong non-thermal component (Section 3). The mass loss of 30 Dor is estimated at 0.05M /a. It is speculated that the nucleus of a galaxy may be wandering due to explosive events.     

Evolutionary models of nucleosynthesis in the galaxy

A model of the galaxy is constructed and evolved in which the integrated influence of stellar and supernova nucleosynthesis on the composition of the interstellar gas is traced numerically. Our detailed assumptions concerning the character of the matter released from evolving stars and supernovae are guided by the results of recent stellar evolutionary calculations and hydrodynamic studies of supernova events. Stars of main sequence mass in the range 4M8M are assumed to give rise to supernova events, leaving remnants we identify with neutron stars and pulsars and forming both the carbon-to-iron nuclei and ther-process heavy elements in the explosive ejection of the core material. For more massive stars, we assume the core implosion will result in the formation of a Schwarzschild singularity, that is, a black hole or collapsar. The straightforward assumptions (1) that the gas content of the galaxy decreases exponentially with time to its present level of 5% and (2) that the luminosity function characteristic of young clusters and the solar neighborhood is appropriate throughout galactic history, lead to the prediction that 20% of the unevolved stars of approximately one solar mass (M ) in the galaxy today should have metal compositionsZ0.1Z . As Schmidt has argued from similar reasoning, this is quite inconsistent with current observations; an early generation dominated by more massive stars—which would by now have evolved—is suggested by this difficulty. Many of these massive stars, according to our assumptions, will end their lives as collapsed black hole remnants. It is difficult to visualize an epoch of massive star formation in the collapsing gas cloud which formed our galaxy which would enrich the gas rapidly enough to account for the level of heavy element abundances in halo population stars; we have therefore proposed a stage of star formation which is entirely pregalactic in character. We suggest that the Jeans' length-sized initial condensations in the expanding universe discussed by Peebles and Dicke may provide the appropriate setting for this first generation of stars. Guided by these considerations, and by the need for a substantial quantity of unseen mass to bind our local group of galaxies, we have constructed a model of the galaxy in which this violent early phase of massive star formation produces both (1) approximately 25% of the level of heavy elements observed in the solar system and (2) an enormous unseen mass in the form of black holes. The implications of our model for other features of the galaxy, including supernova nucleosynthesis, the cosmic ray production of the light elements, and cosmochronology, are discussed in detail.     

High order symplectic integrators for perturbed Hamiltonian systems

A family of symplectic integrators adapted for the integration of perturbed Hamiltonian systems of the form H=A+B was given in (McLachlan, 1995). We give here a constructive proof that for all integer p, such integrator exists, with only positive steps, and with a remainder of order O(^p + ²²), where is the stepsize of the integrator. Moreover, we compute the analytical expressions of the leading terms of the remainders at all orders. We show also that for a large class of systems, a corrector step can be performed such that the remainder becomes O(^p +⁴²). The performances of these integrators are compared for the simple pendulum and the planetary three-body problem of Sun–Jupiter–Saturn.