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.     

Tianjin photographic zenith tube observations

In this paper we analyse the observational data obtained by the Chinese-made PZT in the two periods 1979 Feb – 1980 May and 1981 Dec – 1983 March. The internal accuracy of single star is found to be m_u = ±13.0 ms, m_φ = ±0. “144 for the first period, and m_u = ±14.6 ms, m_φ = ±0.” 152 for the second. Correction of star position is found by the chain method. The systematic error caused by the sealed window of the evacuated chamber and the temperature effect of the plate scale are investigated. Monthly means of time and latitude are given.     

Cyclotron amplification of geomagnetic micropulsations Pc1 in the magnetosphere

A numerical analysis of cyclotron instabilities is carried out by computing the dispersion relation for a three component cold plasma-beam system. Rates of growth and damping for various values of the stream density are calculated from the dispersion relation. The rates of growth and damping increase monotonically as the number density of the proton stream increases. It is found that the frequencies at the rates of maximum growth and the damping decrease slightly to lower frequencies and a sharp peak at these frequencies becomes blunt. The minimum e-folding times of an ion cyclotron wave for (a) σ_s = 10⁻⁴, σ_i = 10⁻² and (b) σ_s = 10⁻¹, σ_i = 10⁻² are about 3·84 and 0·16 sec respectively in the vicinity of the equatorial plane at 6 Re, where σ_s and σ_i are the ratios of the beam density N_s and the helium ion (H₆⁺) density N_i to the total positive ions in the plasma-beam system.     

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°.     

Small-scale auroral arc distortions

Small-scale spatially periodic distortions of auroral forms have been studied utilizing low-light level television observations made at various locations in the Northern and Southern Hemispheres. The most commonly observed features were folds and vortex-like curl formations. The curls, identified here with the Kelvin-Helmholtz instability due to fluid shear, invariably had a counterclockwise rotational shape and motion when viewed in a direction anti-parallel to the Earth's magnetic field. The typical measured wavelength (5 km) and measured growth rate (4.2 sec⁻¹) were used to evaluate the Kelvin-Helmholtz dispersion relation for the apparent shear ω_s = ∂ ν_x/ ∂y (28 sec⁻¹). The apparent horizontal velocities of both folds (0–5 km/sec) and curls (0–22 km/sec) were invariably observed to be counterclockwise with respect to the multiple arc centre when viewed antiparallel to B. Consistent agreement between rotational shape and rotational motion suggests that the apparent growth rate and the apparent horizontal velocities closely approximate the actual values. If the shear results from E×B drifts in a space charge field, the calculated value for ω_s, implies an unneutralized electron density 0–1 cm⁻³ and a ΔE across the arc element 500mV/m. The velocity measurements indicate that the ΔE values for individual elements can combine to produce transient electric fields at the edges of multiple arcs as high as 1000 mV/m.     

Cataclysmic variables III. AC Cancri

AC Cnc is a nova-like, eclipsing binary of period 7^h13^m. I chose it for observation because its eclipses are rather symmetrical. A photometric solution gives inclination i = 74.5° ± 0.8°, mass of white dwarf M₁ = 0.74 ± 0.07 M, mass of the late-type companion, M₂ = 0.97 ± 0.8 M. Temperature of the accretion disk varies approximately as inverse half-power of the radial distance, the temperature at the edge of the disk is 7600 K. Rate of mass transfer from the late-type star to the white dwarf is 7(−9) M/yr. The distance of AC Cnc is 500 ± 100 pc.     

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.     

Changes in thermospheric molecular oxygen abundance inferred from twilight 6300 A airglow

When the local solar zenith angle, χ_L, is < 105° the 6300 A line is much stronger than expected on the basis of F region ionic recombination alone. Between 95 and 105° the additional intensity is quantitatively explained by production of O(¹D) from photolysis of O₂ in the Schumann-Runge continuum, (λλ 1300–1750 A) using current values for solar flux, atmospheric composition and quenching of O(¹D) by N₂. The Schumann-Runge (SR) component exhibits a large seasonal variation with a maximum in summer. We interpret this variation as implying a seasonal change in thermospheric O₂ abundance; the change seems largely to reflect a variation in O₂ density at the base of the diffusive regime although some contribution may come from changes in thermospheric temperature structure. Large changes in the SR component exist from day to day and with a 27 day period following a major magnetic storm. The photodissociation source becomes inadequate when x_l < 95°; at 90° more than half of the intensity comes from still another source which we identify as local photoelectron excitation of O atoms.     

The embedded cluster DBSB 48 in the nebula Hoffleit 18: Comparison with Trumpler 14

We derive fundamental parameters of the embedded cluster DBSB 48 in the southern nebula Hoffleit 18 and the very young open cluster Trumpler 14, by means of deep JHK_s infrared photometry. We build colour–magnitude and colour–colour diagrams to derive reddening and age, based on main sequence and pre-main sequence distributions. Radial stellar density profiles are used to study cluster structure and guide photometric diagram extractions. Field-star decontamination is applied to uncover the intrinsic cluster sequences in the diagrams. Ages are inferred from K-excess fractions. A prominent pre-main sequence population is present in DBSB 48, and the K-excess fraction f_K = 55 ± 6% gives an age of 1.1 ± 0.5 Myr. A mean reddening of A_Ks=0.9±0.03 was found, corresponding to A_V = 8.2 ± 0.3. The cluster CMD is consistent with the far kinematic distance of 5 kpc for Hoffleit 18. For Trumpler 14 we derived similar parameters as in previous studies in the optical, in particular an age of 1.7 ± 0.7 Myr. The fraction of stars with infrared excess in Trumpler 14 is f_K = 28 ± 4%. Despite the young ages, both clusters are described by a King profile with core radii R_core = 0.46 ± 0.05 pc and R_core = 0.35 ± 0.04 pc, respectively, for DBSB 48 and Trumpler 14. Such cores are smaller than those of typical open clusters. Small cores are probably related to the cluster formation and/or parent molecular cloud fragmentation. In DBSB 48, the magnitude extent of the upper main sequence is ΔK_s ≈ 2 mag, while in Trumpler 14 it is ΔK_s ≈ 5 mag, consistent with the estimated ages.     

The relation between low-latitude neutral density variations near 400 km and magnetic activity indices

Neutral density data were obtained near 400km (1600 LT) from a microphone density gauge on OGO-6 from 0°G to 40°N magnetic latitude for 25 September–3 October 1969. Several geomagnetic storms occurred during this period (a_p varied from 0 to 207). Least-squares fits were made to data points on density-a_p and density-D_st scatter diagrams, where the density values selected were delayed in time behind a_p and D_st. An equation representing the least-squares fit was computed for each delay time. The equation of best fit (and the corresponding time delay between the density and the magnetic index which resulted in this best fit) was found by choosing the equation that gave the minimum standard error. For example, the best fit at 10°N geomagnetic latitude occurred for a_p at t — 3 hr, where t is the time of the density values. The implications of the time differences associated with the best fits at various latitudes and longitudes are discussed with regard to the time delays involved in geomagnetic heating of the neutral upper atmosphere.

A low-latitude density bulge has been found between 0°N and 40°N whose magnitude varies with a_p. DeVries (1972b) has independently discovered this daytime phenomenon. If the bulge is a semi-permanent feature near the equinoxes because of the enhanced geomagnetic activity, this may help explain the semi-annual effect in density, which was uncovered first in the drag data from low inclination satellites.     

Analysis of the variation in inclination of Intercosmos 13 Rocket, 1975-22B

The orbit of Intercosmos 13 rocket (1975-22B) has been determined at 103 epochs between 30 April 1975 and 10 April 1980 from almost 7000 observations. One hundred and three values of inclination have been determined and corrections incoporated for the effects due to zonal harmonic, lunisolar and tesseral harmonic perturbations, precession, and solid Earth tides. The modified data have been analysed to yield values of the atmospheric rotation rate, Λ rev day⁻¹, viz. Λ = 0.94 ± 0.10 at an average height of 322 ± 6 km and Λ = 1.27 ± 0.02 at 288 km. Analysis of the inclination near 14th-order resonance has indicated lumped harmonic values 10^{9 1.0}_1.4 = − 76.13 ± 12.47, 10^{9 1,0}₁₄ = − 29.89 ± 32.64, 10^{9 −1.2}₁₄ = − 63.11 ± 15.44 10^{9 −1.2}₁₄ = − 32.52 ± 26.96, for inclination 82.952°.     

Constraining cosmological models with cluster power spectra

Using extensive N-body simulations we estimate redshift space power spectra of clusters of galaxies for different cosmological models (SCDM, TCDM, CHDM, ΛCDM, OCDM, BSI, τCDM) and compare the results with observational data for Abell–ACO clusters. Our mock samples of galaxy clusters have the same geometry and selection functions as the observational sample which contains 417 clusters of galaxies in a double cone of galactic latitude |b|>30° up to a depth of 240 h⁻¹ Mpc. The power spectrum has been estimated for wave numbers k in the range 0.03k0.2 h Mpc⁻¹. For k>k_max0.05 h Mpc⁻¹ the power spectrum of the Abell–ACO clusters has a power-law shape, P(k)∝kⁿ, with n≈−1.9, while it changes sharply to a positive slope at k<k_max. By comparison with the mock catalogues SCDM, TCDM (n=0.9), and also OCDM with Ω₀=0.35 are rejected. Better agreement with observation can be found for the ΛCDM model with Ω₀=0.35 and h=0.7 and the CHDM model with two degenerate neutrinos and Ω_HDM=0.2 as well as for a CDM model with broken scale invariance (BSI) and the τCDM model. As for the peak in the Abell–ACO cluster power spectrum, we find that it does not represent a very unusual finding within the set of mock samples extracted from our simulations.     

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.     

XMM-Newton observations of the nearby brown dwarf LP 944-20

The nearby (d=5.0 pc) brown dwarf LP 944-20 was observed with the XMM-Newton satellite on 07 January 2001. The target was detected with the Optical Monitor (V=16.736±0.081), but it was not detected during the ≈48 ks observation with the X-ray telescopes. We determine a 3σ upper limit for the X-ray emission from this object of L_X<3.1×10²³ ergs·s⁻¹, equivalent to a luminosity ratio upper limit of log(L_X/L_bol)≤−6.28. This measurement improves by a factor of three the previous Chandra limit on the quiescent X-ray flux. This is the most sensitive limit ever obtained on the quiescent X-ray emission of a brown dwarf. Combining the XMM-Newton data with previous ROSAT and Chandra data, we derive flare duty cycles as a function of their luminosities. We find that very strong flares [Log(L_X/L_bol)>−2.5] are very rare (less than 0.7% of the time). Flares like the one detected by Chandra [Log(L_X/L_bol)=−4.1] have a duty cycle of about 6%, which is lower than the radio flare duty cycle (13%). When compared with other M dwarfs, LP 944-20 appears to be rather inactive in X-rays despite of its relative youth, fast rotation and its moderately strong activity at radio wavelengths.     

Discussion of trignometric parallaxes determined along right ascension and declination

We discuss the data of trignometric parallaxes determined along and δ published by the US Naval Observatory and Van Vleck Observatory. We found that the absolute difference δπ = π_x−π_y shows a very similar periodic variation with for both observations. This variation is shown to be due to errors in π_y and to be related to the current practice of maximising the parallax factor for parallax determination in right ascension. Therefore, if possible, we should increase observations where the parallax factor in declination is maximised.     

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.     

Whistler wave propagation in magnetospheric ducts (in the equatorial region)

An analytical theory of whistler wave propagation in axially symmetric field-aligned density ducts is developed. Both enhancements and rarefactions of the density (crests and troughs) are considered. Simple equations giving the dependence of the number of modes n upon the angular frequency ω are derived. From these results it follows that in density crests n decreases when ω approaches ω_c/2 for ω < ω_c/2. The limiting frequency of the wave trapping is calculated. An analytical investigation of wave attenuation in a density crest due to wave leakage is presented. An analysis of the whistler modes in density troughs for ω < ω_c/2, ω > ω_c/2, and ω → ω_c/2 shows that the number of modes is of the same order of magnitude in these three cases.     

The cosmic ray primary composition in the “knee” region through the EAS electromagnetic and muon measurements at EAS-TOP

The evolution of the cosmic ray primary composition in the energy range 10⁶–10⁷ GeV (i.e. the “knee” region) is studied by means of the e.m. and muon data of the Extensive Air Shower EAS-TOP array (Campo Imperatore, National Gran Sasso Laboratories). The measurement is performed through: (a) the correlated muon number (N_μ) and shower size (N_e) spectra, and (b) the evolution of the average muon numbers and their distributions as a function of the shower size. From analysis (a) the dominance of helium primaries at the knee, and therefore the possibility that the knee itself is due to a break in their energy spectrum (at E_k^He=(3.5±0.3)×10⁶ GeV) are deduced. Concerning analysis (b), the measurement accuracies allow the classification in terms of three mass groups: light (p,He), intermediate (CNO), and heavy (Fe). At primary energies E₀≈10⁶ GeV the results are consistent with the extrapolations of the data from direct experiments. In the knee region the obtained evolution of the energy spectra leads to: (i) an average steep spectrum of the light mass group (γ_p,He>3.1), (ii) a spectrum of the intermediate mass group harder than the one of the light component (γ_CNO2.75, possibly bending at E_k^CNO≈(6–7)×10⁶ GeV), (iii) a constant slope for the spectrum of the heavy primaries (γ_Fe2.3–2.7) consistent with the direct measurements. In the investigated energy range, the average primary mass increases from lnA=1.6–1.9 at E₀1.5×10⁶ GeV to lnA=2.8–3.1 at E₀1.5×10⁷ GeV. The result supports the standard acceleration and propagation models of galactic cosmic rays that predict rigidity dependent cut-offs for the primary spectra of the different nuclei. The uncertainties connected to the hadronic interaction model (QGSJET in CORSIKA) used for the interpretation are discussed.     

Origin of the ten degree Solar System dust bands

The Solar System dust bands discovered by IRAS are toroidal distributions of dust particles with common proper inclinations. It is impossible for particles with high eccentricity (approximately 0.2 or greater) to maintain a near constant proper inclination as they precess, and therefore the dust bands must be composed of material having a low eccentricity, pointing to an asteroidal origin. The mechanism of dust band production could involve either a continual comminution of material associated with the major Hirayama asteroid families, the equilibrium model (Dermott et al. (1984) Nature 312, 505–509) or random disruptions in the asteroid belt of small, single asteroids (Sykes and Greenberg (1986) Icarus 65, 51–69). The IRAS observations of the zodiacal cloud from which the dust band profiles are isolated have excellent resolution, and the manner in which these profiles change around the sky should allow the origin of the bands, their radial extent, the size-frequency distribution of the material and the optical properties of the dust itself to be determined. The equilibrium model of the dust bands suggests Eos as the parent of the 10° band pair. Results from detailed numerical modeling of the 10° band pair are presented. It is demonstrated that a model composed of dust particles having mean semimajor axis, proper eccentricity and proper inclination equal to those of the Eos family member asteroids, but with a dispersion in proper inclination of 2.5°, produces a convincing match with observations. Indeed, it is impossible to reproduce the observed profiles of the 10° band pair without imposing such a dispersion on the dust band material. Since the dust band profiles are matched very well with Eos, Themis and Koronis type material alone, the result is taken as strong evidence in favor of the equilibrium model. The effects of planetary perturbations are included by imposing the appropriate forced elements on the dust particle orbits (these forced elements vary with heliocentric distance). A subsequent model in which material is allowed to populate the inner solar system by a Poynting-Robertson drag distribution is also constructed. A dispersion in proper inclination of 3.5° provides the best match with observations, but close examination of the model profiles reveals that they are slightly broader than the observed profiles. If the variation of the number density of asteroidal material with heliocentric distance r is given by an expression of the form 1/r^τ then these results indicate that γ < 1 compared with γ = 1 expected for a simple Poynting-Robertson drag distribution. This implies that asteroidal material is lost from the system as it spirals in towards the Sun, owing to interparticle collisions.     

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.