Physics of the gravitational effect of the galactic bar and spiral structure is presented. Physical equations differ from the conventionally used equations. Application to the motion of the Sun is treated. The speed of the Sun is taken to be consistent with the Oort constants. Galactic radial migration of the Sun is less than +- 0.4 kpc for the four-armed spiral structure. The Sun remains about 75 % of its existence within galactocentric distances (7.8 - 8.2) kpc and the results are practically independent on the spiral structure strength. Thus, the radial distance changes only within 5 % from the value of 8 kpc. Galactic radial migration of the Sun is less than +- (0.3 - 1.2) kpc, for the two-armed spiral structure. The Sun remains (29 - 95) % of its existence within galactocentric distances (7.8 - 8.2) kpc and the results strongly depend on the spiral structure strength and the angular speed of the spiral arms. The radial distance changes within (3.8 - 15.0) % from the value of 8 kpc. If observational arguments prefer relevant radial migration of the Sun, then the Milky Way is characterized by the two-arm spiral structure.

Laboratory experiments show that dusty bodies in a gaseous environment eject dust particles if they are illuminated. We find that even more intense dust eruptions occur when the light source is turned off. We attribute this to a compression of gas by thermal creep in response to the changing temperature gradients in the top dust layers. The effect is studied at a light flux of 13 kW/(m*m) and 1 mbar ambient pressure. The effect is applicable to protoplanetary disks and Mars. In the inner part of protoplanetary disks, planetesimals can be eroded especially at the terminator of a rotating body. This leads to the production of dust which can then be transported towards the disk edges or the outer disk regions. The generated dust might constitute a significant fraction of the warm dust observed in extrasolar protoplanetary disks. We estimate erosion rates of about 1 kg/s for 100 m parent bodies. The dust might also contribute to subsequent planetary growth in different locations or on existing protoplanets which are large enough not to be susceptible to particle loss by light induced ejection. Due to the ejections, planetesimals and smaller bodies will be accelerated or decelerated and drift outward or inward, respectively. The effect might also explain the entrainment of dust in dust devils on Mars, especially at high altitudes where gas drag alone might not be sufficient.

Orbital evolution of spherical interplanetary dust particles in the Edgeworth-Kuiper belt zone is treated for semimajor axes 30-50 AU. Besides solar gravity, also solar electromagnetic and corpuscular radiation, and, fast interstellar gas flow are important forces influencing motion of the particles. The solar electromagnetic radiation is represented by the Poynting-Robertson effect and the solar corpuscular radiation corresponds to the solar wind. Time-variability of the non-radial solar wind can significantly increase dust lifetime in the zone. The average time for the particle stay in the zone is more than 30-times greater than the conventional case of constant (time independent) radial solar wind offers, for the particles of tens micrometers in size. This holds for the most realistic material properties of the particles: $\bar{Q}'_{pr} =$ 1, where $\bar{Q}'_{pr}$ is the dimensionless efficiency factor for electromagnetic radiation pressure. If $\bar{Q}'_{pr} =$ 1/2, then the average time of the particle stay in the zone is only 4-times the conventional value. The approach used in the paper illustrates the relevance of the solar wind action in comparison with the Poynting-Robertson effect. The results have an important consequence for our understanding of the structure of dust distribution in the Edgeworth-Kuiper belt zone, and, also, of dust belts in other stellar systems.

Orbital evolution of an interplanetary dust particle under action of an interstellar gas flow is investigated. Secular time derivatives of the particle orbital elements, for arbitrary orbit orientation, are presented. An important result concerns secular evolution of semi-major axis. Secular semi-major axis of the particle on a bound orbit decreases under the action of fast interstellar gas flow. Possible types of evolution of other Keplerian orbital elements are discussed. The paper compares influences of the Poynting-Robertson effect, the radial solar wind and the interstellar gas flow on dynamics of the dust particle in outer planetary region of the Solar System and beyond it, up to 100 AU. Evolution of putative dust ring in the zone of the Edgeworth-Kuiper belt is studied. Also non-radial solar wind and gravitational effect of major planets may play an important role. Low inclination orbits of micron-sized dust particles in the belt are not stable due to fast increase of eccentricity caused by the interstellar gas flow and subsequent planetary perturbations - the increase of eccentricity leads to planet crossing orbits of the particles. Gravitational and non-gravitational effects are treated in a way which fully respects physics. As a consequence, some of the published results turned out to be incorrect. Moreover, the paper treats the problem in a more general way than it has been presented up to now. The influence of the fast interstellar neutral gas flow might not be ignored in modeling of evolution of dust particles beyond planets.

The paper deals with several properties of the Oort cloud of comets. Sun, Galaxy (and Jupiter) gravitationally act on the comets. New physical model of galactic tide is considered. The main results can be summarized as follows: 1. Mass of the Oort cloud of comets is less than 1 mass of the Earth ($M_{E}$), probably not greater than 1/2 $M_{E}$. 2. Theoretical number of long-period comets with perihelion distance $q$ $<$ 5 AU is about 50-times greater than the conventional approach yields. Gravity of Jupiter was taken into account in finding this result. 3. Semi-major axis $a$ and period of oscillations $P$ of eccentricity (and other orbital elements) are related as $a^{3}$ $P$ $=$ 1 in natural units for a moving Solar System in the Galaxy. The natural unit for time is the orbital period of the Solar System revolution around the galactic center and the natural unit for measuring the semi-major axis is its maximum value for the half-radius of the Solar System corresponding to the half-radius of the Oort cloud. The relation holds for the cases when comets approach the inner part of the Solar System, e.g., perihelion distances are less than $\approx$ 100 AU. 4. The minimum value of semi-major axis for the Oort cloud is $a_{min}$ $\ll$ 1 $\times$ 10$^{4}$ AU. This condition was obtained both from the numerical results on cometary evolution under the action of the galactic tides and from the observational distribution of long-period comets. If the density function of semi-major axis is approximated by proportionality $a^{\alpha}$, then $\alpha$ is - 1/2, approximately. 5. The magnitude of the change in perihelion distance per orbit, $\Delta q$, of a comet due to galactic tides is a strong function of semi-major axis $a$, proportional to $a^{8.25}$.

Equation of motion for the galactic tide is treated for the case of a comet situated in the Oort cloud of comets. We take into account that galactic potential and mass density depend on a distance from the galactic equator and on a distance from the rotational axis of the Galaxy. Secular evolution of orbital elements is presented. New terms generated by the Sun's oscillation about the galactic plane are considered. The inclusion of the new terms into the equation of motion of the comet leads to orbital evolution which may be significantly different from the conventional approach. The usage of the secular time derivatives is limited to the cases when orbital period of the comet is much less than i) the period of oscillations of the Sun around the galactic equator, and, ii) the orbital period of the motion of the Sun around the galactic center.

Equation of motion and the vector of perturbing acceleration (force) for the galactic tide in a noninertial frame of reference is derived. The noninertial reference frame is rotating with a fixed angular velocity $\vec{\omega}$ $=$ $-$ $\omega_{0}$ $\hat{\vec{z}}$ with respect to the inertial frame of reference of the Galaxy. $\vec{\omega}$ is the angular velocity of the solar rotation (rotation of the Local Standard of Rest) around the galactic center, the unit vector $\hat{\vec{z}}$ is oriented toward the north pole of the Galaxy: the Sun is always situated in the plane $y'$ $=$ 0 ($x' - z'$-plane). The equation of motion can be applied to the dynamics of the Oort cloud of comets. Relations for calculation of the osculating orbital elements are presented and a new integral of motion is derived for the conventional approach in modelling of the effect of the galactic tidal field.

Equation of motion for the galactic tide is derived under the assumption of cylindrically symmetric gravitational potential of the Galaxy. The paper considers galactic tide both for the galactic plane $x-$ and $y-$ components and also for the normal $z-$ component. The $x-$ and $y-$ components of the acceleration come not only from the $x-$ and $y-$ components of the position of a body, but also from its $z-$component of the position vector. Values of the Oort constants are $A$ $=$ (14.2 $\pm$ 0.5) $km s^{-1} kpc^{-1}$ and $B$ $=$ ($-$ 12.4 $\pm$ 0.5) $km s^{-1} kpc^{-1}$. %(the values hold for the galactocentric distance of the Sun). Mass density in the solar neighborhood, 30 $pc$ above the galactic equatorial plane, equals to (0.117 $\pm$ 0.005) $M_{\odot} pc^{-3}$. The result for the acceleration is written in the form easily applicable to Solar System studies, to the evolution of comets in the Oort cloud.

Equation of motion for a comet in the Oort cloud is numerically solved. Orbital evolution of the comet under the action of the gravity of the Sun and the Galaxy is presented for various initial conditions. Oscillations of the Sun with respect to the galactic equatorial plane are taken into account. Real values of physical quantities concerning the gravitational action of the galactic neighbourhood of the Sun are important. The results are compared with currently used more simple models of the galactic tide. It turns out that physically improved models yield results which significantly differ from the results obtained on the basis of the conventional models. E.g., the number of returns of the comets into the inner part of the Solar System are about two times greater than it is in the conventional models. It seems that a comet from the Oort cloud can be a source of the dinosaurs extinction at about 65 Myr ago. A close encounter of a star or an interstellar cloud disturbed a comet of the Oort cloud in the way that its semi-major axis increased/decreased above the value 5 $\times$ 10$^{4}$ AU and the comet hit the Earth.

The effect of solar wind on dust particle, meteoroid, is investigated. Rotation of the particle is also considered. Detail derivations of both equation of motion and Euler's dynamical equations, are presented. The most simple form of the action of the solar wind is considered and the derived equation of motion reduces to "pseudo-Poynting-Robertson effect" for non-rotating spherical particle. The final equation of motion contains also angular velocity term and it may cause change of orbital plane of the meteoroid. Analysis of the Euler's dynamical equations shows that rotational axis of the meteoroid changes in time and equilibrium state for the orientation of the axis of rotation does not exist for spherical particle.

Action of solar wind on arbitrarily shaped interplanetary dust particle is investigated. The final relativistically covariant equation of motion of the particle contains both orbital evolution and change of particle's mass. Non-radial solar wind velocity vector is also included. The covariant equation of motion reduces to the Poynting-Robertson effect in the limiting case when spherical particle is treated, the speed of the incident solar wind corpuscles tends to the speed of light and the corpuscles spread radially from the Sun. The results of quantum mechanics have to be incorporated into the physical considerations, in order to obtain the limiting case. The condition for the solar wind effect on motion of spherical interplanetary dust particle is $\vec{p}'_{out}$ $=$ (1 $-$ $\sigma'_{pr} / \sigma'_{tot}$) $\vec{p}'_{in}$, where $\vec{p}'_{in}$ and $\vec{p}'_{out}$ are incoming and outgoing radiation momenta (per unit time) measured in the proper frame of reference of the particle; $\sigma'_{pr}$ and $\sigma'_{tot}$ are solar wind pressure and total scattering cross sections. Real flux density of solar wind energy produces shift of perihelion of interplanetary dust particles. This result significantly differs from the standard treatment of the action of the solar wind on dust particles, when analogy with the Poynting-Robertson effect is stressed. Moreover, the evolution of the shift of perihelion depends on orbital position of the parent body at the time of ejection of the particle.

Physics of the Poynting-Robertson (P-R) effect is discussed and compared with the statements published in the past thirty years. Relativistically covariant formulation reveals the essence of the P-R effect and points out to nonphysical explanations in scientific papers and monographs. Although the final equation of motion $m$ $d\vec{v} / dt$ $=$ ($S A'\bar{Q'}_{pr}$ $ / $ $c$) ${(1 - \vec{v} \cdot \vec{e} / c) \vec{e} - \vec{v} / c}$ has been usually correctly presented and used, its derivation and explanation of its essence is frequently incorrect. The relativistically covariant form of the equation of motion yields the P-R effect as an action of the radiation pressure force on a moving spherical body. No "P-R drag", as a particular relativistically covariant equation of motion, exists. Omission of the nonphysical term "P-R drag" excludes any confusion in the published definitions. The difference between the effects of solar electromagnetic and corpuscular (solar wind) radiation is stressed. The force acting on the particle due to the solar wind (the simple case of radial solar wind velocity is considered) is $\vec{F}_{sw}$ $=$ $F_{sw}$ [ (1 $-$ $\vec{v} \cdot \vec{e} / v_{sw}$) $\vec{e}$ $-$ $x'$ $\vec{v} / v_{sw}$ ], where $F_{sw}$ is the force on the stationary particle, $v_{sw}$ is the heliocentric solar-wind speed, and, the value of $x'$ depends on material properties of the particle (1 $<$ $x'$ $<$ 3). Secular evolution of orbital elements is presented. Initial conditions are included.

Spherical meteoroids orbiting the Sun experience orbital decay due to the Poynting-Robertson effect and can become trapped in commensurability resonances with a planet. Detail conditions for trapping in the planar circular restricted three-body problem with action of solar electromagnetic radiation on spherical grain are presented. The results are given in terms of the radiation presure factor $\beta$. In general, the greater value of $\beta$, the shorter capture time in a resonance. Captures of particles may be even several times longer than it is presented in literature. Analytical formula for maximum capture times is given. Analytical theory for minimal capture eccentricities is presented. The obtained analytical results differ from those presented in the literature. However, the known analytical results are not very consistent with our detail numerical computational results for $\beta-$values 0.01 and 0.05 and for the planet Earth. Minimum eccentricities for captures into resonances are also smaller, even in order of magnitude, than the published values. The presented results enable comparison with calculations for nonspherical particles.

Equation of motion of an uncharged arbitrarily shaped dust particle under the effects of (stellar) electromagnetic radiation and thermal emission is derived. The resulting relativistically covariant equation of motion is expressed in terms of standard optical parameters. Relations between energy and mass of the incoming and outgoing radiation are obtained, together with relations between radiation energy and mass of the particle. The role of the diffraction nicely fits the relativistic formulation of the momentum of the outgoing radiation. The inequality 0 < $\bar{Q}'_{pr, 1} / \bar{Q}'_{ext}$ < 2 is a simple relativistic consequence for the Poynting-Robertson (P-R) effect ($\bar{Q}'_{ext}$ and $\bar{Q}'_{pr, 1}$ are dimensionless efficiency factors for the extinction and radial direction of the radiation pressure, integrated over stellar spectrum). The condition for the P-R effect is $\vec{p}'_{o}$ = (1 - $\bar{Q}'_{pr, 1} / \bar{Q}'_{ext}$) $\vec{p}'_{i}$, where $\vec{p}'_{i}$ and $\vec{p}'_{o}$ are incoming and outgoing radiation momenta (per unit time) measured in the proper frame of reference of the particle. The case of "perfectly absorbing spherical dust particle", within geometrical optics approximation, corresponds to the condition $\vec{p}'_{o}$ = 0.5 $\vec{p}'_{i}$. As for arbitrarily shaped dust particle, the condition 0 < $\bar{C}'_{pr, 1}$ / $\bar{C}'_{ext}$ < 2 / ($1 ~+~ \sum_{j=2}^{3} \bar{C}'_{pr, j} / \bar{C}'_{pr, 1}$) holds for cross sections of extinction and radiation pressure components. The condition can add a new information to the results obtained from observations, measurements and numerical calculations of the optical properties of the particle.

A paradox associated with the astrophysical Poynting-Robertson effect is presented. The paradox arises when relativity theory and Mie's solution of Maxwell's equations are confronted with the statements on the Poynting-Robertson effect. Although the relevant physics has been known already for a century (Poynting 1903, Einstein 1905, Mie 1908), nobody has been aware of the inconsistency between the theories.

Effect of stellar electromagnetic radiation on motion of spherical dust particle in mean-motion orbital resonances with a planet is investigated. Planar circular restricted three-body problem with the Poynting-Robertson (P-R) effect yields monotonous secular evolution of eccentricity when the particle is trapped in the resonance. Elliptically restricted three-body problem with the P-R effect enables nonmonotonous secular evolution of eccentricity and the evolution of eccentricity is qualitatively consistent with the published results for the complicated case of interaction of electromagnetic radiation with nonspherical dust grain. Thus, it is sufficient to allow either nonzero eccentricity of the planet or nonsphericity of the grain and the orbital evolutions in the resonances are qualitatively equal for the two cases. This holds both for exterior and interior mean-motion orbital resonances. Evolutions of longitude of pericenter in the planar circular and elliptical restricted three-body problems are shown. Our numerical integrations suggest that any analytic expression for secular time derivative of the particle's longitude of pericenter does not exist, if a dependence on semi-major axis, eccentricity and longitude of pericenter is considered (the P-R effect and mean-motion resonance with the planet in circular orbit is taken into account). Change of optical properties of the spherical grain with the heliocentric distance is also considered. The change of the optical properties: i) does not have any significant influence on secular evolution of eccentricity, ii) causes that the shift of pericenter is mainly in the same direction/orientation as the particle motion around the Sun. The statements hold both for circular and noncircular planetary orbits.

We consider the motion of uncharged dust grains of arbitrary shape including the effects of electromagnetic radiation and thermal emission. The resulting relativistically covariant equation of motion is expressed in terms of standard optical parameters. Explicit expressions for secular changes of osculating orbital elements are derived in detail for the special case of the Poynting-Robertson effect. Two subcases are considered: (i) central acceleration due to gravity and the radial component of radiation pressure independent of the particle velocity, (ii) central acceleration given by gravity and the radiation force as the disturbing force. The latter case yields results which may be compared with secular orbital evolution in terms of orbital elements for an arbitrarily shaped dust particle. The effects of solar wind are also presented.

We consider the motion of uncharged dust grains of arbitrary shape including the effects of electromagnetic radiation and thermal emission. The resulting relativistically covariant equation of motion is expressed in terms of standard optical parameters. Explicit expressions for secular changes of osculating orbital elements are derived in detail for the special case of the Poynting-Robertson effect. Two subcases are considered: (i) central acceleration due to gravity and the radial component of radiation pressure independent of the particle velocity, (ii) central acceleration given by gravity and the radiation force as the disturbing force. The latter case yields results which may be compared with secular orbital evolution in terms of orbital elements for an arbitrarily shaped dust particle. The effects of solar wind are also presented.

The Poynting-Robertson (P-R) effect has been applied to meteoroids for several decades. It is well known that the P-R effect produces only changes in the orbital plane of the particle. The differential equations governing the secular changes in both eccentricity and semi-major axis are known since the time of Robertson, provided the initial orbits are not near circular. A disadvantage of this type of osculating orbital elements is that it cannot be applied to motion of arbitrarily shaped particles. Relevant type of osculating orbital elements is discussed. It's application to the P-R effect is presented and simple analytical formulae for secular changes in both eccentricity and semi-major axis are derived.

Relativistically covariant form of equation of motion for arbitrarily shaped dust particle (neutral in charge) under the action of electromagnetic radiation is derived -- emission, scattering and absorption of radiation is considered. The result is presented in the form of optical quantities used in optics of dust particles. The obtained equation of motion represents a generalization of the Poynting-Robertson (P-R) effect, which is standardly used in orbital evolution of dust particles in astrophysics. Simultaneous action of electromagnetic radiation and gravitational fields of the central body -- star -- on the motion of the particle is discussed.

The force acting on a test particle moving around a body that is gravitating and emitting electromagnetic radiation at the same time is computed up to the second order in $1/c$, and the shift of perihelion is found. The total shift is positive (the perihelion moves in the direction of particle revolution) and it increases from zero to infinity as the radiation pressure to gravity ratio increases from zero to 1. The additional shift caused by the radiation pressure contains a term inversely proportional to eccentricity, therefore it may dominate over the general relativistic shift if the orbit is round enough.

Equation of motion for real dust particle under the action of electromagnetic radiation is more general than equation of motion corresponding to standardly used Poynting-Robertson effect (P-R effect). As a consequence, orbital evolution of particles may significantly differ from that corresponding to the P-R effect. The paper discusses recently published (Icarus, June 2002) derivation of equation of motion, which is in contradiction with known relativistically covariant formulation. The ``new'' derivation does not respect fundamental physical laws (law of conservation of energy, law of conservation of momentum) which must hold in any frame of reference. Application of the derived ``general'' equation of motion to the special case treated by Einstein in 1905 yields result which is not consistent with Einstein's result. Correct solution is presented.

Arguments of astronomers against equation of motion for arbitrarily shaped particle under action of electromagnetic radiation are discussed. Each of the arguments is commented in detail from the point of view of the required physics. It is shown that the arguments of astronomers, including referees in several astronomical and astrophysical journals, are unacceptable from the physical point of view. Detail explanations should help astronomers in better physical understanding of the equation of motion. Relativistically covariant equation of motion for real dust particle under the action of electromagnetic radiation is derived. The particle is neutral in charge. Equation of motion is expressed in terms of particle's optical properties, standardly used in optics for stationary particles.

Discussion of different types of osculating orbital elements for motion of real dust particle under the action of electromagnetic radiation in the central gravitational field is presented. It is shown that physically correct access is based on gravitational acceleration as the only radial acceleration -- ``radiation pressure'' is not included in the radial acceleration.

Relativistically covariant form of equation of motion for real particle (body) under the action of electromagnetic radiation is derived. Equation of motion in the proper frame of the particle uses the radiation pressure cross section 3 $\times$ 3 matrix. Obtained covariant equation of motion is compared with another covariant equation of motion which was presented more than one year ago.

Relativistically covariant form of equation of motion for real particle (neutral in charge) under the action of electromagnetic radiation is derived. Various formulations of the equation of motion in the proper frame of reference of the particle are used. Main attention is devoted to the reformulation of the equation of motion in the general frame of reference, e. g., in the frame of reference of the source of electromagnetic radiation. This is the crucial form of equation of motion in applying it to motion of particles (cosmic dust, asteroids, ...) in the Universe if electromagnetic radiation acts on the particles. General relativistic equation of motion is presented.

Detailed discussion of university textbook statements (Harwit 1988) concerning the Poynting-Robertson effect (P-R effect) is presented. Discussion is concentrated on better physical understanding of the P-R effect. References to complete correct equation of motion for real dust particle are also presented.

Covariant form of equation of motion for arbitrarily shaped particle in the electromagnetic radiation field is presented. Equation of motion in the proper frame of the particle uses the radiation pressure cross section 3 $\times$ 3 matrix. The obtained equation of motion is compared with known result.

Relativistically covariant equation of motion for real dust particle under the action of electromagnetic radiation is derived. The particle is neutral in charge. Equation of motion is expressed in terms of particle's optical properties, standardly used in optics for stationary particles.

Correct and complete (to terms of $\vec{v} / c$ -- $\vec{v}$ is particle's velocity, $c$ is the speed of light) derivation of equation of motion for real dust particle under the action of electromagnetic radiation is derived. The effect of aberration of light is used. Equation of motion is expressed in terms of particle's optical properties, standardly used in optics for stationary particles.

Effects of solar wind and solar electromagnetic radiation on motion of asteroids are discussed. The results complete the statements presented in Vokrouhlick\'{y} and Milani (2000). As for the effect of electromagnetic radiation, the complete equation of motion is presented to the first order in $v/c$ -- the shape of asteroid (spherical body is explicitly presented) and surface distribution of albedo should be taken into account. Optical quantities must be calculated in proper frame of reference.

Comments to the statements of Srikanth (1999) are presented. As for the standard definition of the Poynting-Robertson drag, the results of Srikanth (1999) are incorrect. Srikanth's statements about the isothermality condition and ``red-shift'' are also discussed. Srikanth's results are generalized for the most general case of the equation of motion when momentum loss per unit time is proportional to -v/c (higher orders are neglected), where v is (heliocentric) velocity of the particle, c is the speed of light.

Criteria for the membership of individual meteors in meteor streams are discussed from the point of view of their mathematical and also physical properties. Discussion is also devoted to the motivation. It is shown that standardly used criteria (mainly D-criterion of Southworth and Hawkins, 1963) have unusual mathematical properties in the sense of a term ``distance'', between points in a phase space, and, physical motivation and realization for the purpose of obtaining their final form is not natural and correct, and, moreover, they lead also to at least surprising astrophysical results. General properties of possible criteria are discussed. A new criterion for the membership in meteor streams is suggested. It is based on probability theory. Finally, a problem of meteor orbit determination for known parent body is discussed.

A simple model for motion of dust particle (meteoroid) under the action of (solar) electromagnetic radiation is presented. The particle of the form of plane mirror is taken into account and exact analytical results are presented. As for long-term orbital evolution, particle may spiral outwards the central body (Sun); initial conditions are important. As a consequence, motion of real dust particles may differ from that generally considered.

Solutions of the two-body problem with the simultaneous action of the solar electromagnetic radiation in the form of the Poynting-Robertson effect are discussed. Special attention is devoted to pseudo-circular orbits and terminal values of osculating elements. The obtained results complete those of Klacka and Kaufmannova (1992) and Breiter and Jackson (1998). Terminal values of osculating elements presented in Breiter and Jackson (1998) are of no physical sense due to the fact that relativistic equation of motion containing only first order of $\vec{v}/c$ was used in the paper.

The problem of (non)random distribution of points on the sphere is investigated. Published procedures for obtaining preferred direction and preferred plane for points on the sphere (in the sky) are discussed. It is shown that the published methods are incorrect, and, as a consequence, the results obtained by these methods cannot be considered to be significant. The correct methods and their applications on real data will be presented in other papers of this set of papers.

The problem of (non)random distribution of points on sphere and in space is investigated. The procedure for obtaining preferred direction (and plane) for points on the sphere (in the sky) and in the space is discussed. At present, directions of perihelia of the observed long-period comets cannot be considered to be distributed uniformly in the sky (contrary of the statement made by Neslu\v{s}an 1996). The action of galactic disk should increase the concentration in the zone of galactic latitude $b \approx 0^{\circ}$, unless the observational selection effects play an important role.

The problem of concentration of cometary perihelia on the sphere is discussed from the point of view of physics and astronomy. The important facts which may be crucial in understanding the observed concentration of cometary perihelia in the direction of solar apex are presented.

Criterion for the membership of individual meteors to meteoroid streams presented by Valsecchi {\it et. al.} (1999) and Jopek {\it et. al.} (1999) is discussed. The authors characterize and use their criterion as a distance function. However, it is not a distance function. Some practical aspects are also discussed. Correct criterion is presented.

Application of the theoretical results of the author are presented for the case of meteor streams of comet Encke. Theoretical determination of meteor orbit for comet Encke, as a parent body, is presented. Four significant theoretical meteor streams, corresponding to Tauds N, Tauds S, $\beta$ Tauds and $\xi$ Perds are found. The meteor stream membership criterion is applied to the photographic orbits of The IAU Meteor Data Center in Lund: Taurid meteor stream is found for several possible areas in phase-space of orbital elements.

Problem of meteor orbit determination for a given parent body is discussed. Some of the published methods for obtaining meteoroid's orbital elements at the moment of intersecting Earth's orbit on the basis of geometrical variation of parent body's orbital elements are discussed. The main result concerns the following two facts: i) in real situations physical quantities for the change of orbital elements of the parent body must be used, and, ii) the usage of Southworth and Hawkins (1963) D-criterion yields results not corresponding to observations.

The importance of the Poynting-Robertson effect on the motion of interplanetary dust particles is discussed. Precise numerical calculations for real dust particle show that condition for the validity of the Poynting-Robertson effect is not fulfilled. The interaction of the (solar) electromagnetic radiation with really shaped dust particle is different from that which yields the Poynting-Robertson effect. The magnitude of the Poynting-Robertson effect's deceleration term is in one to two orders in magnitude (it depends on particle's size) less important than terms corresponding to nonforward (or, nonbackward) scattering.

Motion of a point mass in gravitational fields of the Sun and of the galactic disk is studied. Fundamental features of the motion are found by investigating the time-averaged differential equations for orbital evolution. Several types of possible orbits are mathematically exactly derived in a strictly analytical way. The relation $a^{3} ~ P^{2} = f (e_{0}, i_{0}, \omega_{0})$ between semimajor axis a and period P of the change of osculating orbital elements is found (the index 0 denotes initial values of the quantities). Due to conservation of energy in potential fields a is a constant. Moreover, the component of angular momentum perpendicular to the galactic plane is conserved. Due to these facts the system of equations reduces to two equations for either (e, $\omega$), or (i, $\omega$) (the length of the ascending node does not enter the equations for a, e, i, $\omega$ and is not solved here).

The effect of nonradial component of solar wind is discussed from the qualitative point of view. It is shown that the direction of nonradial component is opposite in comparison with the direction used in papers dealing with orbital evolution of meteoroids.

The method of constructing of rotation curve for our galactic disc is discussed. Simple physical model shows that the velocity of the Sun around the center of the Galaxy must be less than currently accepted IAU value 220 km/s. Independently obtained result for neutral hydrogen (\lambda = 21.1 cm) confirms the previous statement. The values of the Sun's velocity is about 110 km/s. As a consequence, the rotation curve for the disc of our Galaxy is not flat, but it is a decreasing function of the distance from the center of the Galaxy. This conclusion seems to be consistent also with other galaxies.