Automated Stellar Classification for Large Surveys with EKF and RBF Neural Networks

An automated classification technique for large size stellar surveys is proposed. It uses the extended Kalman filter as a feature selector and pre-classifier of the data, and the radial basis function neural networks for the classification. Experiments with real data have shown that the correct classification rate can reach as high as 93%, which is quite satisfactory. When different system models are selected for the extended Kalman filter, the classification results are relatively stable. It is shown that for this particular case the result using extended Kalman filter is better than using principal component analysis.     

Free-convection flow in the stokes problem for a vertical plate with variable surface temperature

The free-convection flow of an incompressible viscous fluid past an infinite, vertical plate for impulsive as well as uniformly accelerated motion of the plate is discussed when the plate temperatuve varies as the square root of time. The Laplace transform technique is used to obtain the expressions for velocity and skin-friction. The influence of the various parameters, entering into problem, on the velocity field and skin-friction is extensively discussed.     

On the concept of the russian information and analytical center for asteroid and comet hazards

The concept is presented of the Russian information and analytical center for asteroid and comet hazards (IAC ACH) as a part of the national IAC for space threats. An overview is given of the center’s goals and objectives, as well as the requirements for the center, its structure, and the possible cooperation between the institutes of the Russian Academy of Sciences and other organizations in the creation of the center.     

Solar rotation and the giant cells

Departures from the mean solar differential rotation rate as a function of latitude, longitude, and epoch of the solar cycle, together with variations in the rotation rate as determined by spectroscopic and tracer measurements are reviewed. It is shown that, if giant convection cells do exist as predicted, real variations in the subsurface rotation rate should occur and that this may be responsible for the observed surface anomalies.In terms of this hypothesis, a simple account is given for the anomalous rotation rates of sunspots. Furthermore, the torsional oscillations are identified as a modulation of the differential rotation produced by a system of toroidal convective rolls generated near the poles and propagating towards the equator. It is suggested that, as these rolls progress through lower latitudes, they break up into a system of cells which are the long sought for giant cells of the convection zone. Thus the torsional oscillations are identified as direct surface evidence for the existence of these cells.Solar Cycle Workshop Paper.     

Radiative Transfer in Inhomogeneous Atmospheres. III

The approach proposed in the previous parts of this series of papers is used to solve the radiative transfer problem in scattering and absorbing multicomponent atmospheres. Linear recurrence relations are obtained for both the reflectance and transmittance of these kinds of atmospheres, as well as for the emerging intensities when the atmosphere contains energy sources. Spectral line formation in a one-dimensional inhomogeneous atmosphere is examined as an illustration of the possibility of generalizing our approach to the matrix case. It is shown that, in this case as well, the question reduces to solving an initial value problem for linear differential equations. Some numerical calculations are presented.     

Solar Si XI Line Ratios Observed by the Normal Incidence Spectrometer on SOHO CDS

New measurements of line intensity ratios in the Be-like ion Si xi are presented for observations of the quiet Sun, active regions, coronal holes and above-limb regions obtained using the Coronal Diagnostic Spectrometer on the Solar and Heliospheric Observatory. A model ion, constructed using the best available atomic data, is used to predict the line intensity ratios for a wide range of electron temperatures and densities. Comparisons of the theoretical ratios with the new intensity ratios as well as with those from previous solar observations and laboratory measurements are given. The usefulness of the ratios for electron temperature and density diagnostics, as well as for spectrometer calibration, is discussed.     

Non-thermal ionization and recombination processes during solar flares

Ionization and recombination processes are studied for a plasma of which the electrons follow a power-law energy distribution.The rates for collisional ionization, radiative and dielectronic recombination and for autoionization are evaluated.Numerical computations are performed for H-like, He-like and Li-like ions from neon to nickel as a function of the spectral index of the electron distribution. The ionization equilibrium is evaluated as well as the ratios of fluxes emitted in two lines pertaining to two successive ionization stages of the same element. A comparison with a few experimental data is made and the possibility of a non-thermal interpretation of X-ray line emission during solar flares is discussed.     

Periodicities of the number of solar flares for the period 1966–1988

In the present paper a study is made of the mean monthly number of grouped solar flaresf for the time period 1966–1988. Corresponding dataF were taken from the catalogue published by Coffey (1989). Periodicities of 140, 120, 48, 18, 12, and 11 months as well as shorter periodicities of the order of 6 and 3 months for the solar flares have been found. The emphasis is given as far as it concerns the period of 48 months, i.e., 4 years which is for the first time revealed by the present investigation.     

Kernel regression for determining photometric redshifts from Sloan broad-band photometry

We present a new approach, namely kernel regression, to determine photometric redshifts for 399 929 galaxies in the Fifth Data Release of the Sloan Digital Sky Survey (SDSS). Kernel regression is a weighted average of spectral redshifts of the neighbours for a query point, and higher weights are associated with points that are closer to the query point. One important design decision when using kernel regression is the choice of bandwidth. We apply 10-fold cross-validation to choose the optimal bandwidth, which is obtained as the cross-validation error approaches its minimum. The results show that the optimal bandwidth is different for different input patterns: the lowest rms error of photometric redshift estimation arrives at 0.019 using colour+eClass as the inputs, the lowest rms errors comes to 0.020 using ugriz +eClass as the inputs. Where eClass is a galaxy spectral type, and 0.021 using colour+ r as the inputs. Thus, in addition to parameters such as magnitude and colour, eClass is a valid parameter with which to predict photometric redshifts. Moreover, the results suggest that the accuracy of estimating photometric redshifts is improved when the sample is divided into early-type and late-type galaxies; in particular, for early-type galaxies, the rms scatter is 0.016 with colour+eClass as the inputs. In addition, kernel regression achieves high accuracy when predicting the photometric eClass  (σ_rms= 0.034)  using colour+ r as the input pattern. For kernel regression, the accuracy of the photometric redshifts does not always increase with the number of parameters considered, but is satisfactory only when appropriate parameters are chosen. Kernel regression is a comprehensible and accurate regression method. Experiments reveal the superiority of kernel regression over other empirical training approaches.     

A new way to conserve total energy for Eulerian hydrodynamic simulations with self-gravity

We propose a new method to conserve the total energy to round-off error in grid-based codes for hydrodynamic simulations with self-gravity. A formula for the energy flux due to the work done by the self-gravitational force is given, so the change in total energy can be written in conservative form. Numerical experiments with the code Athena show that the total energy is indeed conserved with our new algorithm and the new algorithm is second order accurate. We have performed a set of tests that show the numerical errors in the traditional, non-conservative algorithm can affect the dynamics of the system. The new algorithm only requires one extra solution of the Poisson equation, as compared to the traditional algorithm which includes self-gravity as a source term. If the Poisson solver takes a negligible fraction of the total simulation time, such as when FFTs are used, the new algorithm is almost as efficient as the original method. This new algorithm is useful in Eulerian hydrodynamic simulations with self-gravity, especially when results are sensitive to small energy errors, as for radiation pressure dominated flow.     

A re-assessment of the nearest neighbour alignment of the X-ray isophotes of galaxy clusters

Alignment is defined as the tendency of the distribution of pointing angles between the major axes of clusters and their nearest neighbours to be more concentrated towards small values for small nearest neighbour distances, whereas the distribution is expected to be uniform over all angles at larger distances. Conflicting pronouncements on the reality of this effect have been published in the astronomy literature. A re-assessment of the evidence for alignment is presented, based on three recently published X-ray data sets. We find that whereas there is evidence for alignment, it is not as convincing as previously claimed. In particular, the scale to which the effect has been claimed to extend seems to have been severely overstated.     

A generalised principle of invariance and radiation field in multilayer atmospheres

A generalized principle of invariance is derived for a plane-parallel atmosphere. On the basis of this principle a method for determining the radiation field in a multilayer atmosphere is proposed. This method, the first part of which is the well-known adding method, permits the application to problems involving optically finite as well as semi-infinite atmospheres. The reflecting boundaries may be incorporated, though in that case it is not possible to use the adding method.Some numerical results are given for the standard and Milne problems and for the problem with internal sources.     

Homochirality as the signature of life: the SETH Cigar

A characteristic hallmark of life is its homochirality: all biomolecules are usually of one hand, e.g. on Earth life uses only L-amino acids for protein synthesis and not their D mirror images. It is therefore suggested that a search for extra-terrestrial life can be approached as a Search for Extra-Terrestrial Homochirality (SETH). A novel miniaturized space polarimeter, called the SETH Cigar, is described which could he used to detect optical rotation as the homochiral signature of life on other planets. Moving parts are avoided by replacing the normal rotating polarizer by multiple fixed polarizers at different angles as in the eye of the bee. It is believed that homochirality will be found in the subsurface layers on Mars as a relic of extinct life.     

Inferences of plasma parameters from coronal hole observations

The temperature in the acceleration region of the solar wind remains one of the most elusive parameters to measure. Knowledge of the temperature as well as its gradient in the inner corona is fundamental for placing constraints on physical mechanisms thought to be responsible for the coronal heating process, as well as for understanding the flow properties of the solar wind. Estimates of the helium abundance is essential for understanding the puzzling behavior of heavier ions in the solar wind. As an illustration of the difficulties and uncertainties involved in the inferences of plasma parameters in the wolar wind acceleration region, The inference of electron temperature and helium abundance will be described. Prospects for future observations will be briefly discussed.     

Tilt-angle variations of active regions

An examination of the tilt angles of multi-spot sunspot groups and plages shows that on average they tend to rotate toward the average tilt angle in each hemisphere. This average tilt angle is about twice as large for plages as it is for sunspot groups. The larger the deviation from the average tilt angle, the larger, on average, is the rotation of the magnetic axis in the direction of the average tilt angle. The rate of rotation of the magnetic axis is about twice as fast for sunspot groups as it is for plages. Growing plages and spot groups rotate their axes significantly faster than do decaying plages and spot groups. There is a latitude dependence of this effect that follows Joy's law. The fact that these tilt angles move toward the average tilt angle and not toward 0 deg (the east-west orientation), combined with other results presented here, suggest that a commonly accepted view of the origin of active region magnetic flux at the solar surface may have to be re-examined.Operated by the Association of Universities for Research in Astronomy, Inc., under Cooperative Agreement with the National Science Foundation.     

The age distribution of galactic cosmic rays

The distribution of ages for secondary and interacting cosmic-ray nuclei is derived for a class of steady-state, bounded models of diffusion in the galaxy. Results are presented in detail for the model in which diffusion is in one dimension, sources are uniformly distributed throughout the scattering region, and particles are observed near the central plane of the galaxy. The leakage-lifetime approximation is shown to be accurate as long as the ratio of the cosmic-ray diffusive lifetime to the mean interaction or decay lifetime is less than about ten.     

Removing cosmic-ray hits from CCD images in real-time mode by means of an artificial neural network

A feed-forward artificial neural network has been implemented to the problem of removing cosmic-ray hits (CRH) from CCD images. The results of a number of tests demonstrate the effectiveness of this method especially for undersampled stellar profiles. The problem of optimal and low price preparing of training data, which could enable real-time or at least fast post-processing filtering out of CRH is discussed. The training and test ensembles were composed of a number of synthetic stellar profiles involving different S/N ratios and CRH images taken from real data. Certain aspects of the network’s architecture and its training efficiency for different modes of the back-propagation procedure as well as for the pre-process normalization of data have been examined. It is shown that for training set composed of stellar images and CRH at a ratio of 1:2 recognition can reach 99% in the case of stars and 96% for CRH. To determine the extent to which the cognition power of a network trained using an ensemble of circular symmetric stellar profiles of a given radius can be generalised the test data included stellar profiles of different radii, as well as elongated profiles. The goal was to mimic temporal changes in seeing as well as such problems as image defocusing, the lack of isoplanatism and improper sideral tracking of a telescope. The experiments provided us with the conclusion that for S/N > 10 excellent classification property is maintained in cases where the change in the radius of a circular profile is up to 30%, as well as for elongated profiles where the longest dimension is almost double that of the shortest one. Moreover, the generalization capability has been investigated for test images of synthetic pairs of overlapping stars with different distances between components. Almost 99% recognition efficiency was achieved even if the separation was nearly three times the radius of the stellar profile, a case when two stars could be analyzed by appropriate software as separate objects. The example of removal of CRH from real CCD images is presented to give an idea of how an algorithm based on a neural network can work in practice. The result of such an experiment appears fully consistent with the conclusions drawn from the tests made on synthetic data.     

Investigation of some of the principal geometric effects on planetary polarization

Three major geometric factors which are likely to influence theoretical interpretation of planetary polarization measurements, viz., observer—planet distance, horizontal inhomogeneity of planetary disk, and deviation from a spherical body, are investigated.The distance effect is examined for regional as well as global polarizations. For convenience of analysis, the expressions for zenith and azimuth angles of incident and emergent light appropriate for a snap-shot observation are derived as explicit functions of distance between observer and planet. Sample computations for Venus indicate that regional polarization near the planetary limb is significantly affected by the observer's distance. This effect should be particularly noticeable when an observation is made at a phase angle around which the single scattering polarization of atmospheric scattering agents exhibits a steep variation. The global polarization at large phase angles (measured at disk-center) is gradually moved toward smaller phase angles, as the observer approaches the planet. Any narrow polarization features such as rainbow and glory at small phase angles are heavily smoothed out.The effects of horizontal inhomogeneity are investigated with a planetary disk having highly polarizing regions at high latitudes. Comparison of theoretical global polarization computed for such a disk with the Pioneer Venus OCPP measurements shows a possible change in cloud-haze stratification approximately at 50° latitude, consistent with other imaging observations. An approximate analytical representation of residual polarization at zero phase angle is then derived to compare to the numerical results for Venus. An attempt is also made to explain the relatively large magnitude of residual polarization observed on Jupiter.Finally, to study the effects of nonsphericity of planetary body, the global polarizations are computed for a spheroidal planet. The global polarization tends to increase as the planet's oblateness increases. However, for Jupiter and Saturn, such effect may be of secondary importance.     

Analysis of Extra-Solar Planetary Transits Using Tools from Eclipsing Binaries

Analytical formulae for the computation of the photometric elements of extra-solar planetary transits are presented. They were initially derived for the study of well-detached eclipsing binaries and are valid for any degree of limb darkening and type of transit as well as for eccentric orbits. The only assumption is that the projections of the star and the planet on the plane of the sky are well represented by circular discs. The procedure to get valid ranges for the involved parameters, as well as to make precise estimations of initial parameters, using the analytical nature of the equations is given together with some discussion on their practical application. Examples are shown for OGLE-TR-113, representative of a light curve obtained with ground-based telescopes, and HD 209458 with a precise light curve obtained with the Hubble Space Telescope.     

太阳系天体地面VLBI观测的相对论时延模型

以太阳系质心参考系为基础，根据太阳系的质心参考系和非旋转地球质心参考系的坐标转换关系，推导了太阳系天体地面VLBI观测的相对论时间延迟模型，给出了一个通用的解析表达式．根据这一公式可以得到平劲松博士所采用的公式，以及当地心与源的距离无限大时，可得河外射电源VLBI观测的Zhu—Groten模型、Shapiro模型和IERS(92，96)推荐模型．所推导的公式严格解析且无误差，在实际应用中建议采用这一公式．同时详细地讨论了所推导公式的实用范围和各种舍掉项的量级估计，并详细给出了时间延迟理论模型的计算步骤．