The averaged resonant equations of motion for the planar circular restricted three-body problem are solved on the linearization basis taking into account also non-gravitational effects. The averaged resonant equations are derived from Lagrange's planetary equations with additional Gauss's terms caused by the non-gravitational effects. The time depending solution has the standard form with exponential, quadratic, linear and constant terms. The existence of a rotational symmetry in the action of the non-gravitational effects around the star determines the order of a characteristic equation of the linearized system. In the symmetrical case (order 3) the considered non-gravitational effects are the stellar electromagnetic radiation and the radial stellar wind (stellar radiation). In the asymmetrical case (order 4) the stellar radiation and interstellar gas flow are considered. It is investigated how well the linearization solution describes real solution obtained from an equation of motion by a comparison of the resonant libration frequency found analytically and numerically. It is found that from initial values of the evolving orbital parameters (semimajor axis, eccentricity, longitude of pericenter, and resonant angular variable) in the averaged phase space the linearization frequency depends most sensitively on the initial value of the resonant angular variable. For small libration amplitudes of the resonant angular variable the best match of the real libration frequency and the linearization frequency is located approximately at the solution of the resonant condition ($da / dt$ $=$ 0). If the initial averaged conditions are chosen close to the solution of resonant condition, then the linearization frequency for practically all simple oscillatory evolutions matches the real libration frequency and the linearization solution very well approximates the real evolution.

A structure formed by dust particles ejected from the debris ring around HD 61005 is observed in the scattered light. The main aim here is to constrain interstellar wind parameters that lead to shape morphology in the vicinity of HD 61005 using currently available observational data for the debris ring. Equation of motion of 2 $\times$ 10$^5$ dust particles ejected from the debris ring under the action of the electromagnetic radiation, stellar wind, and interstellar wind is solved. A two-dimensional (2D) grid is placed in a given direction for accumulation of the light scattered on the dust particles in order to determine the shape morphology. The interaction of the interstellar wind and the stellar wind is considered. Groups of unknown properties of the interstellar wind that create the observed morphology are determined. A relation between number densities of gas components in the interstellar wind and its relative velocity is found. Variations of the shape morphology caused by the interaction with the interstellar clouds of various temperatures are studied. When the interstellar wind velocity is tilted from debris ring axis a simple relation between the properties of the interstellar wind and an angle between the line of sight and the interstellar wind velocity exists. Dust particles that are most significantly influenced by stellar radiation move on the boundary of observed structure. Observed structure at HD 61005 can be explained as a result of dust particles moving under the action of the interstellar wind. Required number densities or velocities of the interstellar wind are much higher than that of the interstellar wind entering the Solar system.

The equations of secular evolution for dust grains in mean motion resonances with a planet are solved for stationary points. This is done including both Poynting-Robertson effect and stellar wind. The solutions are stationary in semimajor axis, eccentricity, and resonant angle, but allow the pericentre to advance. The semimajor axis of stationary solutions can be slightly shifted from the exact resonant value. The periodicity of the stationary solutions in a reference frame orbiting with the planet is analytically proved. The existence of periodic solutions in mean motion resonances means that analytical theory enables for dust particles also infinitely long capture times. The stationary solutions are periodic motions to which the eccentricity asymptotically approaches and around which the libration occurs. Using numerical integration of equation of motion are successfully found initial conditions corresponding to the stationary solutions. Numerically and analytically determined shifts of the semimajor axis form the exact resonance for the stationary solutions are in excellent agreement. The stationary solutions can be plotted by locations of pericenters in the reference frame orbiting with the planet. The pericenters are distributed in the space according to properties of dust particles.

Context. The eclipsing binary GU Mon is located in the star-forming cluster Dolidze 25, which has the lowest metallicity measured in a Milky Way young cluster. Aims. GU Mon has been identified as a short-period eclipsing binary with two early B-type components. We set out to derive its orbital and stellar parameters. Methods. We present a comprehensive analysis, including B and V light curves and 11 high-resolution spectra, to verify the orbital period and determine parameters. We use the stellar atmosphere code FASTWIND to obtain stellar parameters and create templates for cross-correlation. We obtain a model to fit the light and radial-velocity curves using the Wilson-Devinney code iteratively and simultaneously. Results. The two components of GU Mon are identical stars of spectral type B1 V, with the same mass and temperature. The lightcurves are typical of an EW-type binary. The spectroscopic and photometric analyses agree on a period of 0.896640 +- 0.000007 d. We determine a mass of 9.0 +- 0.6 Msun for each component and temperatures of 28 000 +- 2 000 K. Both values are consistent with the spectral type. The two stars are overfilling their respective Roche lobes, sharing a common envelope, and therefore the orbit is synchronised and circularised. Conclusions. The GU Mon system has a fill-out factor above 0.8, containing two dwarf B-type stars on the main sequence. The two stars are in a very advanced stage of interaction, with their extreme physical similarity likely due to the common envelope. The expected evolution of such a system will very probably lead to a merger while still on the main sequence.

Context. The early-type binary MY Cam belongs to the young open cluster Alicante 1, embedded in Cam OB3. Aims. MY Cam consists of two early-O type main-sequence stars and shows a photometric modulation suggesting an orbital period slightly above one day. We intend to confirm this orbital period and derive orbital and stellar parameters. Methods. Timing analysis of a very exhaustive (4607 points) light curve indicates a period of 1.1754514 +- 0.0000015 d. High- resolution spectra and the cross-correlation technique implemented in the TODCOR program were used to derive radial velocities and obtain the corresponding radial velocity curves for MY Cam. Modelling with the stellar atmosphere code FASTWIND was used to obtain stellar parameters and create templates for cross-correlation. Stellar and orbital parameters were derived using the Wilson-Devinney code, such that a complete solution to the binary system could be described. Results. The determined masses of the primary and secondary stars in MY Cam are 37.7 +- 1.6 and 31.6 +- 1.4 Msol, respectively. The corresponding temperatures, derived from the model atmosphere fit, are 42 000 and 39 000 K, with the more massive component being hotter. Both stars are overfilling their Roche lobes, sharing a common envelope. Conclusions. MY Cam contains the most massive dwarf O-type stars found so far in an eclipsing binary. Both components are still on the main sequence, and probably not far from the zero-age main sequence. The system is a likely merger progenitor, owing to its very short period.

For a body with negligible mass moving in the gravitational field of a star and one planet in a circular orbit (the circular restricted three-body problem) five equilibrium points exist and are known as the Lagrangian points. The positions of the Lagrangian points are not valid for dust particles because in the derivation of the Lagrangian points is assumed that no other forces beside the gravitation act on the body with negligible mass. Here we determined positions of the equilibrium points for the dust particles in the circular restricted three-body problem with radiation. The equilibrium points are located on curves connecting the Lagrangian points in the circular restricted three-body problem. The equilibrium points for Jupiter are distributed in large interval of heliocentric distances due to its large mass. The equilibrium points for the Earth explain a cloud of dust particles trailing the Earth observed with the Spitzer Space Telescope. The dust particles moving in the equilibrium points are distributed in interplanetary space according to their properties.

Circumstellar dust particles can be captured in a mean motion resonance with a planet and simultaneously be affected by non-gravitational effects. It is possible to describe the secular variations of a particle orbit in the mean motion resonance analytically using averaged resonant equations. We derive the averaged resonant equations from the equations of motion in near-canonical form. The secular variations of the particle orbit depending on the orientation of the orbit in space are taken into account. The averaged resonant equations can be derived/confirmed also from Lagrange's planetary equations. We apply the derived theory to the case when the non-gravitational effects are the Poynting--Robertson effect, the radial stellar wind, and an interstellar wind. The analytical and numerical results obtained are in excellent agreement. We found that the types of orbits correspond to libration centers of the conservative problem. The averaged resonant equations can lead to a system of equations which hold for stationary points in a subset of resonant variables. Using this system we show analytically that for the considered non-gravitational effects, all stationary points should correspond to orbits which are stationary in interplanetary space after an averaging over a synodic period. In an exact resonance, the stationary orbits are stable. The stability is achieved by a periodic repetition of the evolution during the synodic period. Numerical solutions of this system show that there are no stationary orbits for either the exact or non-exact resonances.

The orbital evolution of a dust particle captured in a mean motion resonance with a planet in circular orbit under the action of the Poynting-Robertson effect, radial stellar wind and an interstellar gas flow of is investigated. The secular time derivative of Tisserand parameter is analytically derived for arbitrary orbit orientation. From the secular time derivative of Tisserand parameter a general relation between the secular time derivatives of eccentricity and inclination is obtained. In the planar case (the case when the initial dust particle position vector, initial dust particle velocity vector and interstellar gas velocity vector lie in the planet orbital plane) is possible to calculate directly the secular time derivative of eccentricity. Using numerical integration of equation of motion we confirmed our analytical results in the three-dimensional case and also in the planar case. Evolutions of eccentricity of the dust particle captured in an exterior mean motion resonance under the action of the Poynting-Robertson effect, radial stellar wind for the cases with and without the interstellar gas flow are compared. Qualitative properties of the orbital evolution in the planar case are determined. Two main groups of the secular orbital evolutions exist. In the first group the eccentricity and argument of perihelion approach to some values. In the second group the eccentricity oscillates and argument of perihelion rapidly shifts.

The acceleration of a spherical dust particle caused by an interstellar gas flow depends on the drag coefficient which is, for the given particle and flow of interstellar gas, a specific function of the relative speed of the dust particle with respect to the interstellar gas. We investigate the motion of a dust particle in the case when the acceleration caused by the interstellar gas flow represent a small perturbation to the gravity of a central star. We present the secular time derivatives of the Keplerian orbital elements of the dust particle under the action of the acceleration from the interstellar gas flow for arbitrary orbit orientation. The semimajor axis of the dust particle is a decreasing function of time for an interstellar gas flow acceleration with constant drag coefficient and also for such an acceleration with the linearly variable drag coefficient. The decrease of the semimajor axis is slower for the interstellar gas flow acceleration with the variable drag coefficient. The minimal and maximal values of the decrease of the semimajor axis are determined. In the planar case, when the interstellar gas flow velocity lies in the orbital plane of the particle, the orbit always approaches the position with the maximal value of the transversal component of the interstellar gas flow velocity vector measured at perihelion. The properties of the orbital evolution derived from the secular time derivatives are consistent with numerical integrations of the equation of motion. If the interstellar gas flow speed is much larger than the speed of the dust particle, then the linear approximation of dependence of the drag coefficient on the relative speed of the dust particle with respect to the interstellar gas is usable for practically arbitrary (no close to zero) values of the molecular speed ratios (Mach numbers).

The orbital evolution of a dust particle under the action of a fast interstellar gas flow is investigated. The secular time derivatives of Keplerian orbital elements and the radial, transversal, and normal components of the gas flow velocity vector at the pericentre of the particle's orbit are derived. The secular time derivatives of the semi-major axis, eccentricity, and of the radial, transversal, and normal components of the gas flow velocity vector at the pericentre of the particle's orbit constitute a system of equations that determines the evolution of the particle's orbit in space with respect to the gas flow velocity vector. This system of differential equations can be easily solved analytically. From the solution of the system we found the evolution of the Keplerian orbital elements in the special case when the orbital elements are determined with respect to a plane perpendicular to the gas flow velocity vector. Transformation of the Keplerian orbital elements determined for this special case into orbital elements determined with respect to an arbitrary oriented plane is presented. The orbital elements of the dust particle change periodically with a constant oscillation period or remain constant. Planar, perpendicular and stationary solutions are discussed. The applicability of this solution in the Solar system is also investigated. We consider icy particles with radii from 1 to 10 micrometers. The presented solution is valid for these particles in orbits with semi-major axes from 200 to 3000 AU and eccentricities smaller than 0.8, approximately. The oscillation periods for these orbits range from 10^5 to 2 x 10^6 years, approximately.

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

Orbital motion of a body can be found from Newtonian equation of motion. However, it is useful to express the motion through time derivatives of Keplerian orbital elements, mainly if the motion is perturbed by small perturbing force. The first set of equations for the time derivatives of the orbital elements can be derived from the equation of motion using Lagrange brackets. The second one by using equation of motion and perturbation acceleration decomposed to radial, transversal and normal components. This paper shows that the second type of the perturbation equations can be derived from the first type using simple mathematical operations.

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.

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.