A kinematic study of the Galaxy

Using the proper motion and parallax data for 1011 O-B stars in the Hipparcos Catalogue we have derived the Oort constants, A = 17.60 ± 0.21 (km/s)/kpc, B = −14.62 ± 0.20 (km/s)/kpc, and a solar velocity V = 16.7 ± 0.10 km/s in the direction l = 45.3° ± 2.8°, b = 21.0° ± 2.3°. For a galactocentric distance of the sun of R₀ = 8.5 kpc, we then get a galactic rotational velocity of the solar neighbourhood of V_lsr = 273.9 km/s, obviously much higher than the IAU published value of 220 km/s. We have investigated the cause for this difference.     

Galactic kinematics derived from classical cepheids

On the basis of radial velocity and Hipparcos proper motion data, we have analyzed the galactic kinematics of classical Cepheids. Using the 3-D Ogorodnikov-Milne model we have determined the rotational velocity of the Galaxy to be V₀ = 240.5 ± 10.2 km/s, on assuming a glactocentric distance of the Sun of R₀ = 8.5 kpc. The results clearly indicate a contracting motion in the solar neighbourhood of (∂V_θ∂θ)/R = −2.60 ± 1.07 km s⁻¹ kpc⁻¹, along the direction of galactic rotation. Possible reason for this motion is discussed. The solar motion found here is S = 18.78 ± 0.86 km/s in the direction l = 54.4° ± 2.9° and b = +26.6° ± 2.6°.     

Properties of the infrared dark cloud G79.2+0.38

The MSX infrared dark cloud G79.2+0.38 has been observed over a 11′×′ region simultaneously in the J=1-0 rotational transition lines of the ¹²CO and its isotopic molecules ¹³CO and ¹⁸CO. The dense molecular cores defined by the C¹⁸O line are found to be associated with the two high-extinction patches shown in the MSX A-band image. The two dense cores have the column density N (H₂) (5 – 12) × 10²² cm⁻² and the mean number density n (3 ± 1) × 10⁴ cm⁻³. Their sizes are 1.7 and 1.2 pc in ¹³CO(1-0) line, 1.2 and 0.6 pc in C¹⁸O(1-0) line, respectively. The masses of these cloud cores are estimated to be in the range from 2 × 10² to 2 × 10³ M. The profile of radial mean density of the cloud core can be described by the exponential function ¯n(p) p^{−0.34±0.02}. Compared with the cases of typical optical dark clouds, the abundances of the CO isotopic molecules ¹³CO and C¹⁸O in this MSX infrared dark cloud appear to be depleted by a factor of 4–11, but at present there is no evidence for any obvious variation of the relative abundance ratio X_13/18 between ¹³CO and C¹⁸O with the column density.     

Mangeto-optical effects in sunspots and their vector magnetic field

The system of transfer equations of the four Stokes parameters I, Q, U, V under the action of the magneto-optical effect (i.e. the Unno-Beckers equations) are numerically solved in this paper for the magneto-sensitive lines FeI λλ 6302.499 and 5324.191 using an appropriate sunspot model. The errors in the expressions for the coefficients _r and _W in Beckers' paper [2] have been corrected for. From the results of calculations, features of the profiles of the Stokes parameters dependent on the magnetic vector have been isolated. Our computations also show that the magneto-optical effect should be taken into consideration in the measurement of the vector magnetic fields.

In the fourth section of this paper we have established a simple and convenient method for obtaining-information on the magnetic vector (including the field strength B, its inclination to the line of sight γ and its azimuth χ) from the profiles of the Stokes parameters. It consists of three steps: (1) The value of B is determined from the distance of the highest point in the V-profile from the central line. (2) γ is then found from V_max, i.e maximum value of V. (3) Lastly, the angle χ is found from Q₀, i.e. the value of Q at line centre.     

Electron excitation and auroral emission parameters

Auroral luminosities of the main emission lines in the aurora have been calculated for excitation by an isotopic primary electron flux with spectra of the form J(E) = AE exp (−E/E₁) + B(E₂)E exp (−E/E₁). The variation of emissions from O and N₂⁺ with height are shown, as are the variations of column integrated intensities and pertinent intensity ratios with the characteristic energy E₂, this leading to a method of estimating the electron spectrum from ground observation.     

A prediction for three neutrino masses and mixings and similarity between quark and lepton mixings

If neutrinos have mass, we give reasons for a possible pattern of three (squaed) mass eigenvalues: m₁² (2.8−5.8) (eV)², m₂² 0.01 (eV)², m₃² (1.5−1) × 10⁻⁴ (eV)². The flavor states ν_μ and ν_e are mixtures of the eigenstates with m₂ and m₃ with a significant mixing, corresponding to an effective mixing angle of about 0.45. The ν_τ is nearly the state with m₁; the other two effective mixing angles are about an order of magnitude smaller than 0.45. There is a marked similarity to mixing in the quark sector.     

A hybrid simulation study of magnetic reconnection in anisotropic plasmas

The process of magnetic reconnection in anisotropic plasmas is studied numerically using a 2-dimensional, 3-component hybrid simulation. The results of the calculation show that, when the plasma pressure in the direction perpendicular to magnetic field is larger than that in the parallel direction (e.g. P/P = 1.5), instability may greatly increase, speeding up the rate of reconnection. When P is smaller than P, (e.g., when P/P = 0.6), fire hose instability appears, which will restrain the tearing mode instability and the process of magnetic reconnection.     

The odd zonal harmonics in the Earth's gravitational potential

The odd zonal harmonics in the Earth's gravitational potential are determined by analysing the changes in the eccentricities of six satellites having orbital inclinations spaced as uniformly as possible between 28° and 96°. The most satisfactory representation of the potential is found to be in terms of four coefficients, and their values are, in the usual notation: 10⁶J₃ = −2.56, 10⁶J₅ = −0.15, 10⁶J₇ = −0.44, 10⁶J₉ = 0.12. The resulting potential is compared with that obtained by other authors. Three and five-coefficient solutions are also presented.     

Auroral N2(c′4Σu → aΠg emission

Inspection of recent spectra presented by Sivjee (1983) show evidence of the 0–4 and 0–5 bands of the N₂(c′₄¹Σ_u⁺ → a¹Π_g) Gaydon-Herman system. In conjunction with earlier spectra, it is now possible that this band system is a significant auroral component, with an intensity approx. 7% that of the N2 2P system. The absence in aurorae of the potentially far stronger N₂(c′₄¹Σ_u⁺ → X¹Π_g) system is discussed. It is that the O₂(A³Σ_u⁺ → X³Σ_g⁻) band system is indiscernible in Sivjee's auroral spectra, under conditio the foreground nightglow is expected to be clearly visible. On the other hand, at least one relatively strong O₂(A′³Δ_u → a¹Δ_g) band appears to be present in these spectra.     

Ionospheric electric field and current distribution associated with high altitude electric field inhomogeneities

A quantitative theoretical analysis of electric field and current distributions in the ionosphere is given assuming certain time variable convection field profiles at an altitude of 1250 km. A number of idealized assumptions regarding the ionospheric characteristics are defined and discussed. A qualitative discussion of a quasi-stationary configuration with an approximately curl free electric field is also given. Geomagnetically field aligned current densities i of the order 10⁻⁵−10⁻⁴A/m² are consistent with quite reasonable assumptions about the convection field E. Oscillations in E with periods of the order of 10 sec should readily be generated when σ is large. In the quasi-stationary case there may be a mechanism that strengthens and concentrates i locally under certain conditions. It is found that a number of recent high altitude observations of convection field reversals may be consistent with large potential drops along the magnetic field lines. The solutions obtained as well as some of the basic assumptions are compared with observations.     

Nature of solar-cycle and heliomagnetic-polarity dependence of cosmic rays, inferred from their correlation with heliomagnetic spherical surface harmonics in the period 1976–1985

Correlation of cosmic-ray intensity (I) with the solar magnetic field expanded into the spherical surface harmonics, B_ns(n 9), by Hoeksema and Scherrer has been studied using the following regression equation:

, where are subgroups of B_ns classified in ascending order of n, and τ_i is the time lag of I behind correlation coefficient between the observed and simulated intensities (I_obs, I_sml) in the period 1976–1985 is 0.87 and considerably better than that derived from any single index of solar activity. The lag time τ₃ is greater than others, indicating that the higher order magnetic disturbances effective to the cosmic-ray modulation have a longer lifetime in space than the lower order disturbances. The rigidity spectrum of the cosmic-ray intensity variation responsible for A_I due to the dipole moment is harder than those for others (A₂,A₃), indicating that the lowest order (i.e. largest scale) magnetic disturbances can modulate cosmic rays more effectively than the higher order disturbances. As another result of the present analysis, it has been found that the intensity depends also on the polarity of the polar magnetic field of the Sun; the residual (I_obs−I_sml) of the simulation changes its sign from positive to negative with a time lag (0–5 Carrington rotation periods) behind the directional change of the solar magnetic dipole moment from northward to southward, and has a softer rigidity spectrum than A_iS. The dependence is consistent with the result having been obtained in the previous period, 1936–1976, by one (K.N.) of the present authors. The polarity dependence can be found also in the 22-year variation of the time lags obtained every solar cycle in the period 1936–1985. The theoretical interpretation of these polarity dependences is discussed on the basis of the diffusion-convection-drift model.     

Density profiles of dark matter halos with anisotropic velocity tensors

We present density profiles, that are solutions of the spherical Jeans equation, derived under the following two assumptions: (i) the coarse grained phase-density follows a power-law of radius, ρ/σ³ ∝ r⁻, and (ii) the velocity anisotropy parameter is given by the relation β_a(r)=β₁+2β₂ (r/r_*)/[1+(r/r_*)²] where β₁, β₂ are parameters and r_* equals twice the virial radius, r_vir, of the system. These assumptions are well motivated by the results of N-body simulations. Density profiles have increasing logarithmic slopes γ, defined by γ=−d ln ρ/d ln r. The values of γ at r=10^−2.5r_vir, a distance where the systems could be resolved by large N-body simulations, lie in the range 1.0–1.6. These inner values of γ increase for increasing β₁ and for increasing concentration of the system. On the other hand, slopes at r=r_vir lie in the range 2.42–3.82. A model density profile that fits well the results at radial distances between 10⁻³r_vir and r_vir and connects kinematic and structural characteristics of spherical systems is described.     

Application of the complete linearization method to the analysis of solar prominence emission spectrum

We analysed the emission spectra of solar prominences using the complete linearization method [5] and found simultaneously the optical depth at the line centre τ₀, the doppler width of the line Δλ_D and the damping width a. The results show 1) that the complete linearization method has a larger radius of convergence, 2) that we must consider the variation of the source function with depth, when determining τ₀, and 3) that the calculated values of the damping constant for the H, Hβ of hydrogen and H and K lines of Calcium are all much greater than the theoretical values from doppler broadening and radiation damping, showing that other mechanisms besides these two also contribute to the broadening of prominence lines.     

The H-alpha variations of Epsilon Aurigae from the second to the fourth contact

Epsilon Aurigae was observed in Mar., Aug., 1983 and Jan., Mar., 1984 with Reticon at coudé focus of the McDonald 2.7m and 2.1m telescopes. Fifty-six observations were taken over eight nights during 2nd to 4th contact.

The H shows noticeable variations in profile, radial velocity, and equivalent width of both absorption and emission components. A similar phenomenon was detected during the corresponding phase of the 1955 – 1957 eclipse.

A model is proposed, in which the primary is a FO supergiant surrounded by a thin ring (or disk) of radius R = 450R and rotating at velocity V_sini = 60 – 70 kms⁻¹. This is the major source of the H emission component. The secondary is a type B star surrounded by a very extended envelope (R = 1000 R) and its rotational speed is similar to that of the ring of the primary. The B star heats a portion of the envelope comparable in size with the primary with enough hydrogen atoms in the lowest excited states, which cause the absorption of the emission from the primary.     

Spectrum of the relic neutrino background from past supernovae and cosmological models

It is greatly expected that the relic neutrino background from past supernovae will be detected by Superkamiokande (SK) which is now under construction. We calculate the spectrum and the event rate at SK systematically by using the results of simulations of a supernova explosion and reasonable supernova rates. We also investigate the effect of a cosmological constant, Λ, on the spectrum, since some recent cosmological observations strongly suggest the existence of Λ. We find following results. (1) The spectrum has a peak at about 3 MeV, which is much lower than that of previous estimates (6–10 MeV). (2) The event rate at SK in the range from 10 MeV to 50 MeV, where the relic neutrinos from past supernovae are dominant, is about 25h₅₀²(R_SN/0.1 yr⁻¹)(n_Gh₅₀⁻³/0.02 Mpc⁻³) events per year, where R_SN is the supernova rate in a galaxy, n_G is the number density of galaxies, and h₅₀ = H₀/(50 km/s Mpc), where H₀ is the Hubble constant. (3) The event rate is almost insensitive to Λ. The flux increases in the low energy side (< 10 MeV) with increasing Λ, but decreases in the high energy side (> 10 MeV) in models in which the integrated number of supernovae in one galaxy is fixed.     

A further study of the method of identifying absorption line systems in QSO spectra

The method of identifying absorption line systems in QSO spectra (Cui et al. 1983; Chen et al. 1983) is further developed here. Certain limitations of the method and their improvements are discussed. Certain other problems requiring further study are pointed out. The improved method is applied to PKS 0528-250, and gives two new absorption line systems Z_a = 0.065 and 0.0345 in addition to the four systems Z_a = 2.8110, 2.8130, 2.5275, 2.1410, consistent with the systems A₁A₂, B, C of Norton et al. (1980). However, the systems D₁, D₂, E, F and G of Chen and Norton (1984) are not recovered. The reason for this discrepancy is discussed.     

Lithium depletion in the solar atmosphere

Using a complete non-local convection theory, we carried out the theoretical calculations of ⁷Li depletion of the solar convective envelope models with different convective parameters c₁ and c₂, and got a model of the solar convection zone consistent with the observed ⁷Li abundance and the depth of the solar convection zone determined by helioseismic techniques. The overshooting distance of effective non-local convective mixing of ⁷Li is very extensive, which is about 1.07H_P or 0.09R. However, the super-radiative temperature zone is much narrower, and it is only 0.20H_P or 0.016R.     

Constraining the power spectrum using clusters

We analyze an extended redshift sample of Abell/ACO clusters and compare the results with those coming from numerical simulations of the cluster distribution, based on the truncated Zel'dovich approximation (TZA), for a list of eleven dark matter (DM) models. For each model we run several realizations, so that we generate a set of 48 independent mock Abell/ACO cluster samples per model, on which we estimate cosmic variance effects. Other than the standard CDM model, we consider (a) Ω₀ = 1 CDM models based on lowering the Hubble parameter and/or on tilting the primordial spectrum; (b) Ω₀ = 1 Cold + Hot DM models with 0.1 ≤Ω_ν ≤0.5; (c) low-density flat ΛCDM models with 0.3 ≤Ω₀ ≤0.5. We compare real and simulated cluster distributions by analysing correlation statistics, the probability density function, and supercluster properties from percolation analysis. We introduce a generalized definition of the spectrum shape parameter Γ in terms of σ₂₅/σ₈, where σ_ris the rms fluctuation amplitude within a sphere of radius r. As a general result, we find that the distribution of galaxy clusters provides a constraint only on the shape of the power spectrum, but not on its amplitude: a shape parameter 0.18 Γ 0.25 and an effective spectral index at the 20 h⁻¹ Mpc scale −1.1 n_eff −0.9 are required by the Abell/ACO data. In order to obtain complementary constraints on the spectrum amplitude, we consider the cluster abundance as estimated using the Press-Schechter approach, whose reliability is explicitly tested against N-body simulations. By combining results from the analysis of the distribution and the abundance of clusters we conclude that, of the cosmological models considered here, the only viable models are either Cold + Hot DM ones with 0.2 Ω_ν 0.3, better if shared between two massive ν species, and ΛCDM ones with 0.3 Ω₀0.5.     

A spectroscopic orbit for the two—Spectrum eclipsing binary RT Persei

Intensified Reticon spectra have been obtained at a high dispersion for the Algol system, RT Persei. They were measured by the cross-correlation technique. The spectroscopic elements, revised for the primary component and determined for the secondary for the first time, are: T₀ = HJD 2,446,038.9332, K₁ = 55.0, K₂ = 194.7, V₀ = −8.3 km/s. A mass ratio q = m₂/m₁ = 0.282 is deduced. A circular orbit is adopted. The spectrum of the primary is F5V, and the secondary is a subgiant. With the elements determined here and the published photometric parameters, the absolute dimensions of the binary are: A = 4.20, R₁ = 1.20, R₂ = 1.08 R; M₁ = 1.08, M₂ = 0.30 M.     

Galactic cosmic-ray anisotropy and its heliospheric modulation, inferred from the sidereal semidiurnal variations observed in the rigidity range 300–600 GV with multidirectional muon telescope at Sakashita underground station

The existence of sidereal semidiurnal variation of cosmic-ray intensity in a rigidity region 10²-10³ GV has been reported by many researchers, but there is no consensus of opinion on its origin. In this paper, using the observed semidiurnal variations in a rigidity range (300–600 GV) with 10 directional muon telescopes at Sakashita underground station (geog. lat. = 36°, long. = 138°E, DEPTH = 80 m.w.e.), the authors determine the magnitudes (η₁, η₂) and directions (a₁, a₂) of the first- and second-order anisotropies in the following galactic cosmic-ray intensity distribution (j)

jdp = j₀{1 + η₁P₁(cos χ₁) + η₂P₂(cos χ₂)}dp

, where P_nis the nth order spherical function and χ_n is the pitch angle of cosmic rays with respect to a_n. For the determination, the influence of cosmic-ray's heliomagnetospheric modulation, geomagnetic deflection and nuclear interaction with the terrestrial material and also of the geometric configuration of the telescopes are taken into account. Usually, the semidiurnal variation is produced by the second-order anisotropy. The present observation, however, requires also the first-order anisotropy which usually produces only the diurnal variation, but can produce also the semidiurnal variation as a result of the heliospheric modulation. The first- and second-order anisotropies are characterized with η₁) > 0 and η₂ < 0 have almost the same direction (a₁ a₂) specified by the right ascension ( 0.75 h) and declination (δ 50°S) and, therefore, they can be expressed, as a whole, by an axis-symmetric anisotropy of loss-cone type (i.e. deficit intensities in a cone). It is noteworthy that this anisotropy approximately coincides with that inferred from the air shower observation at Mt Norikura in the rigidity region 10⁴ GV.