The evidence of an early stellar encounter in Edgeworth-Kuiper belt

We investigate the influence of a stellar fly-by encounter on the Edgeworth-Kuiper belt objects through numerical orbital calculations, in order to explain both mass depletion and high orbital inclinations of the classical Edgeworth-Kuiper belt (CEKB) objects, which have semimajor axis of 42-48 AU and perihelia beyond 35 AU. The observationally inferred total mass of the CEKB is ∼1/10 Earth masses, which is only ∼0.02 of that extrapolated from the minimum-mass solar nebula model. The CEKB consists of bimodal population: “hot population” with inclinations i?0.2-0.6 radians and “cold population” with i?0.1. The observationally suggested difference in size and color of objects between the two populations may imply different origins of the two populations. We find that both the depletion of solid materials in the CEKB and the formation of the hot population are accounted for by a single close stellar encounter with pericenter distance of 80-100 AU and inclination relative to the initial protoplanetary disk ?50°-70°. Such a stellar encounter highly pumps up eccentricities of most objects in the CEKB and then their perihelia migrate within 35 AU. These objects would be removed by Neptune's perturbations after Neptune is formed at or migrates to the current position (30 AU). Less than 10% of the original objects remain in stable orbits with small eccentricities and perihelion distances larger than 35 AU, in the CEKB, which is consistent with the observation. We find that i of the remaining objects are as large as that of the observed hot population. The only problem is how to stop Neptune's migration at ∼30 AU, which is addressed in a separate paper. The depletion by the stellar encounter extends deeply into ∼30-35 AU, which provides the basis of the formation model for the cold population through Neptune's outward migration by Levison and Morbidelli (2003, Nature, 426, 419-421). The combination of our model with Levison and Morbidelli's model could consistently explain the mass depletion, truncation at 50 AU, bimodal distribution in i, and differences in size and color between the hot and the cold populations in the CEKB.     

The layered structure of lunar maria: Identification of the HF-radar reflector in Mare Serenitatis using multiband optical images

Comparison of the Lunar Radar Sounder (LRS) data to the Multiband Imager (MI) data is performed to identify the subsurface reflectors in Mare Serenitatis. The LRS is FM-CW radar (4–6 MHz) and the 2 MHz bandwidth leads to the range resolution of 75 m in a vacuum, whereas the sampling interval in the flight direction is about 75 m when an altitude of the spacecraft with polar orbit is nominal (100 km). Horizontally continuous reflectors were clearly detected by LRS in limited areas that consist of about 9% of the whole maria. The typical depth of the reflectors is estimated to be a few hundred meters. Layered structures of mare basalts are also discernible on some crater walls in the MI data of the visible bands (VIS). The VIS range has nine wavelengths of 415, 750, 900, 950, and 1000 nm, and their spatial resolution is 20 m/pixel at a nominal altitude. The stratigraphies around Bessel and Bessel-H craters in Mare Serenitatis are examined in this paper. It was revealed that the subsurface reflectors lie on the boundaries between basalt units with different chemical compositions. In addition, model calculations using the simplified radar equation indicate that the subsurface reflectors are not compositional interfaces but layer boundaries with a high-porosity contrast. These results suggest that the detected reflectors in Mare Serenitatis are regolith accumulated during so long hiatus of mare volcanisms enough for chemical composition of magma to change, not instantaneously. Therefore combination of the LRS and MI data has a potential to reveal characteristics of a series of magmatism forming each lithostratigraphic unit in Mare Serenitatis and other maria.     

On the encounter desorption of hydrogen atoms on an ice mantle

At low temperatures(10 K),hydrogen atoms can diffuse quickly on grain ice mantles and frequently encounter hydrogen molecules,which cover a notable fraction of grain surface.The desorption energy of H atoms on H2 substrates is much less than that on water ice.The H atom encounter desorption mechanism is adopted to study the enhanced desorption of H atoms on H_(2 )substrates.Using a small reaction network,we show that the steady-state surface H abundances predicted by the rate equation model that includes H atom encounter desorption agree reasonably well with the results from the more rigorous microscopic Monte Carlo method.For a full gas-grain model,H atom encounter desorption can reduce surface H abundances.Therefore,if a model adopts the encounter desorption of H atoms,it becomes more difficult for hydrogenation products such as methanol to form,but it is easier for C,O and N atoms to bond with each other on grain surfaces.     

On the depths and shapes of the freshest kilometer-scale simple craters on the lunar maria: A new crater shape model

Recent work on the shapes of small, simple impact craters on the Moon has shown that the parabolic ideal does not well represent the vast majority of these craters. They are hyperbolic in shape and usually resemble a cone more than a parabola. A parabolic shape also does not fit the most commonly held archetype for simple craters in general (Linné), which is also hyperbolic. In addition, Linné itself may not be the best model for fresh simple craters, in terms of cross-sectional shape, although shape data to compare it to have heretofore been lacking. Here, the “free shadowfront method” for determining the shapes of simple craters is used to measure 64 fresh simple craters on five lunar maria to test both assumptions. Laser altimetry cross sections, available for many of the craters measured herein, are used to complement and spot-check the shadow measurement results, and thereby demonstrate the efficacy of the free shadowfront method. A new shape model is established, and two craters that better fit this model than Linné are identified. These are located at 24.45° N/328.12° E and 31.35° N/296.46° E and have diameters of 1.40 and 2.73 km, respectively. An apparent dichotomy between fresh simple craters smaller than 2.5 km and those larger than this is observed. Flat floors are found to be ubiquitous among the larger craters, but rare and small in extent in smaller ones. A slide in one crater which appears to be an incipient flat floor suggests a major mode of formation for these flat floors.     

The circular polarization of sources of synchrotron radiation

The degree of circular polarizationp _c is calculated for two models of a source of synchrotron radiation:

1. A source with an inhomogeneous magnetic field and isotropic angular distribution of the electrons with respect to the magnetic field;
2. A source with a homogeneous magnetic field and anisotropic angular distribution of the electrons in which the anisotropy of angular distribution substantially increases with the electron energy.

The first model can be used to describe extended radio-sources; and the second, to describe compact radio-sources. For those sources, whose observed polarization properties correspond to the first model, we obtain an integral equation which connects the observed distribution of the sources with the extent of their linear and circular polarization (p _l andp _c ) and the unknown distribution of the sources over the strengthB and the degree of homogeneity ?=(B ₀/B)² of the magnetic field;B ₀ is a homogenous field,B ₀?B. A solution of the integral equation obtained is found for a particular case. This solution makes it possible to determine the distribution of different types of sources over ? if the distribution of these sources in the extent of linear polarization is known. The formulae obtained make it possible to indicate which sources with a known degree of linear polarization should be expected to exhibit highest circular polarization. In the discussion of the first model the question is raised as to the information one can get about the magnetic field by using observations of both linear and circular polarization for a separate source, and for a number of sources. It is shown that the determination of the most probable values ofB and ? in a separate source based on the known values ofp _l andp _c for the source, is possible only if one knows the distribution overB and ? of the sources of the type to which the source in question belongs. The observational data now available make it possible to find the distribution of the sources only over ?. Since the distribution overB and ? is at present unknown, even a very strong upper limit forp _c in the case of a separate source does not enable us to give an exact upper limit for the strength of the magnetic field in this source. In the first model the upper limit for the magnetic field can be obtained only if the upper limit ofp _c is known for a certain number of sourcesN, withN?1. This limit allows for much stronger fields than are usually admitted. This last fact should be taken into consideration when one deals with the results of observations of circular polarization in sources with strong magnetic fields. The first model presents some difficulties when we compare it with observations of some compact sources. The second model can explain why one observes in these sources a violation of the lawp _c ～v ^?1/2 and a change of sign inp _c when the frequency of the observationsv changes.     

Dynamical friction on a satellite of a disk galaxy: The circular orbit

In a previous paper, we have studied dynamical friction during a parabolic passage of a companion galaxy past a disk galaxy. This paper continues with the study of satellites in circular orbits around the disk galaxy. Simulations of orbit decay in a self gravitating disk are compared with estimates based on two-body scattering theories; the theories are found to give a satisfactory explanation of the orbital changes. The disk friction is strongly dependent on the sense of rotation of the companion relative to the rotation of the disk galaxy as well as on the amount of mass in a spherical halo. The greatest amount of dynamical friction occurs in direct motion if no spherical halo is present. Then the infall time from the edge of the disk is about one half of the orbital period of the disk edge. A halo twice as massive as the disk increases the infall time four fold. The results of Quinn and Goodman, obtained with a non-self-gravitating method, agree well with our experiments with massive halos (Q ₀ 1.5), but are not usable in a more general case. We give analytic expressions for calculating the disk friction in galaxies of different disk/halo mass ratios.     

The dust distribution within the inner coma of comet P/Halley 1982i: encounter by Giotto's impact detectors

Analysis of the data from Giotto's Dust Impact Detection System experiment (DIDSY) is presented. These data represent measurement of the size of dust grains incident on the Giotto dust shield along its trajectory through the coma of comet P/Halley on 1986 March 13/14. First detection occurred at some 287000 km distance from the nucleus on the inbound leg; the majority of the DIDSY subsystems remained operational after closest approach (604 km) yielding the last detection at about 202000 km from the nucleus. In order to improve the data coverage (and especially for the smallest grains, to approximately 10(-19) kg particle mass), data from the PIA instrument has been combined with DIDSY data. Flux profiles are presented for the various mass channels showing, to a first approximation, a 1/R2 flux dependence, where R is the distance of the detection point from the cometary nucleus, although significant differences are noted. Deviations from this dependence are observed, particularly close to the nucleus. From the flux profiles, mass and geometrical area distributions for the dust grains are derived for the trajectory through the coma. Groundbased CCD imaging of the dust continuum in the inner coma at the time of encounter is also used to derive the area of grains intercepted by Giotto. The results are consistent with the area functions derived by Giotto data and the low albedo of the grains deduced from infrared emission. For the close encounter period (-5 min to +5 min), the cumulative mass distribution function has been investigated, initially in 20 second periods; there is strong evidence from the data for a steepening of the index of the mass distribution for masses greater than 10(-13) kg during passage through dust jets which is not within the error limits of statistical uncertainty. The fluences for dust grains along the entire trajectory is calculated; it is found that extrapolation of the spectrum determined at intermediate masses (cumulative mass index alpha = 0.85) is not able to account for the spacecraft deceleration as observed by the Giotto Radio Science Experiment and by ESOC tracking operations. Data at large masses (>10(-8) kg) recently analysed from the DIDSY data set show clear evidence of a decrease in the mass distribution index at these masses within the coma, and it is shown that such a value of the mass index can provide sufficient mass for consistency with the observed deceleration. The total particulate mass output from the nucleus of comet P/Halley at the time of encounter would be dependent on the maximum mass emitted if this change in slope observed in the coma were also applicable to the emission from the nucleus; this matter is discussed in the text. The flux time profiles have been converted through a simple approach to modeling of the particle trajectories to yield an indication of nucleus surface activity. There is indication of an enhancement in flux at t approximately -29 s corresponding to crossing of the dawn terminator, but the flux detected prior to crossing of the dawn terminator is shown to be higher than predicted by simple modelling. Further enhancements corresponding to jet activity are detected around +190 s and +270 s.     

The first large format Si:Sb BIB arrays

The first 128×128 Si:Sb blocked impurity band (BIB) detectors, manufactured by Rockwell International, are sensitive detectors from 10 to at least 40 m. While further work is required to make these arrays suitable for the low backgrounds of space infrared telescopes, they can be used now for observations from the ground and aircraft.     

The product of the globular cluster collapse

In this paper we calculate the number of close binaries formed during the evolution process of a globular cluster core. The globular cluster core is assumed to contain a massive black hole at its center. We show that the central black hole can drive binaries formation in the core and the rate of binaries formation depends on the mass of the black hole at its center. When the massM of the black hole is between 10² M and 3×10³ M , there will be a few binaries formed. When the mass of the black hole is 4×10³ M M6×10³ M , the number of binary star formation will suddenly increase with a jump to the maximum value 58. When the mass of the black hole is 7×10³ M M9×10³ M , the number of binary star will immediately decrease. Whether cluster X-ray is produced mainly by the central black hole or by binaries in the core depends on the mass of the central black hole. Therefore, two cases arise: namely, black hole accretion domination and binaries radiation domination. We do think that we cannot exclude the possibility of the existence of a central black hole even when binary radiation characteristics have been observed in globular cluster X-ray sources.     

The alteration of the magnetic field of transient variations by an axisymmetric island in a large ocean

It is shown how the density of electric currents flowing in a horizontal conducting slab (the sea) may be calculated by B.G. Galerkin's method of moments, when a uniform electric field acts in the slab. The slab is uniform except for a hole, symmetrical about a vertical axis (an island). The cross-section of the hole may be arbitrary. Electric current streamlines are drawn for an island whose shape is a frustum of a cone. The magnetic field of the disturbing currents is numerically computed using successive over-relaxation. These results are compared with a two dimensional model which is a modification of the model due to Ashour and Chapman (1965). It is shown that introducing the vertical current component in the sea affects the magnetic field at the surface of the island by only a few per cent of the total magnetic field. The computer model is applied to the study of the island effect on Oshima and Jarvis Islands. The computed results for Oshima Island explain the observed enhancement of Zat some points and the reversal, in fact, of Z at other points on the island.     

Interactive image retrieval in large astronomical archives: The Aspect system

We present a method to access efficiently and effectively images in large astronomical archives. This method allows users to interactively participate in the retrieval process. They may virtually 'navigate' through the archive for choosing images that are relevant for their purposes. The navigation process is performed by switching back and forth between three retrieval-system modes: a query mode, a browse mode and an inspection mode. The method has been applied to an image retrieval system to access solar radio spectrograms archived at ETH in Zurich. The system, called ASPECT, allows access to over 50000 data sets on-line.     

The stability of vertical motion in the N-body circular Sitnikov problem

We present results about the stability of vertical motion and its bifurcations into families of 3-dimensional (3D) periodic orbits in the Sitnikov restricted N-body problem. In particular, we consider ν = N ? 1 equal mass primary bodies which rotate on a circle, while the Nth body (of negligible mass) moves perpendicularly to the plane of the primaries. Thus, we extend previous work on the 4-body Sitnikov problem to the N-body case, with N = 5, 9, 15, 25 and beyond. We find, for all cases we have considered with N ≥ 4, that the Sitnikov family has only one stability interval (on the z-axis), unlike the N = 3 case where there is an infinity of such intervals. We also show that for N = 5, 9, 15, 25 there are, respectively, 14, 16, 18, 20 critical Sitnikov periodic orbits from which 3D families (no longer rectilinear) bifurcate. We have also studied the physically interesting question of the extent of bounded dynamics away from the z-axis, taking initial conditions on x, y planes, at constant z(0) = z ₀ values, where z ₀ lies within the interval of stable rectilinear motions. We performed a similar study of the dynamics near some members of 3D families of periodic solutions and found, on suitably chosen Poincaré surfaces of section, “islands” of ordered motion, while away from them most orbits become chaotic and eventually escape to infinity. Finally, we solve the equations of motion of a small mass in the presence of a uniform rotating ring. Studying the stability of the vertical orbits in that case, we again discover a single stability interval, which, as N grows, tends to coincide with the stability interval of the N-body problem, when the values of the density and radius of the ring equal those of the corresponding system of N ? 1 primary masses.     

The general equilibria of a spinning satellite in a circular orbit

Let a rigid satellite move in a circular orbit about a spherically symmetric central body, taking into account only the main term of the gravitational torque. We shall investigate and find all solutions of the following problem: Let the satellite be permitted to spin about an axis that is fixed in the orbit frame; the satellite need not be symmetric, the spin not uniform, and the spin axis not a principal axis of inertia. The complete discussion of this type of spin reveals that the cases found by Lagrange and by Pringle - and the well-known spin about a principal axis of inertia orthogonal to the orbit plane — are essentially the only ones possible; the only further (degenerate) case is uniform spin of a two-dimensional, not necessarily symmetric satellite about certain axes that are orthogonal to the plane containing the body and to the orbit of the satellite around the central body.     

Noble gas study of a philippinite with an unusually large bubble

Abstract— A philippinite with a large bubble (volume ?5 cm³) was crushed in a specifically constructed device under vacuum. The total pressure, the major gas components, and the concentration and the isotopic composition of noble gases in the trapped gas were measured. The total gas pressure was found to be low, ～10^?4 atm. Compared to terrestrial atmosphere, the abundances of He and Ne are considerably more enriched than the concentrations of the heavy noble gases. Unusually high ²⁰Ne/²²Ne and ²¹Ne/²²Ne ratios were measured in the bubble gas. These ratios are higher than values expected for any steady state process and are interpreted to reflect very fast nonsteady state diffusion in the early stages of tektite formation.