Contribution of Vacuum Field to Angular Deviation of Light Path and Radar Echo Delay

The discovery of ‘twin quasistellar objects’ arose interests among astronomers and astrophysicists to study gravitational lensing problems. The deviation of light from its straight line path is caused by two sources according to the general theory of relativity: (i) the presence of massive objects, i.e. the presence of gravitational field and (ii) the presence of a ‘vacuum field’ which arises because there is a non-zero cosmological vacuum energy. Recently, the research on the relationship between cosmological constant and gravitational lensing process is rather active (see reference [1, 2, 3]. According to the Kottler space time metric, we have deduced an explicit representation of the angular deviation of light path. The deviation term is found to be simply , where M is the mass of the ‘astronomical lens’, r_min is the distance between the point of nearest approach and the centre of M, other symbols have their usual meaning. The presence of this term may be meaningful to the study of cosmological constant using the concept of gravitational lensing; however more sophisticated analysis awaits. Consider a signal radar to be sent from one planet to another. We have found that the radar echo delay contributed by the existence of the cosmological constant Λ is expressible as This revised version was published online in July 2006 with corrections to the Cover Date.     

Possible semi-annual variation of the Newtonian constant of gravitation

A possible semi-annual variation of the Newtonian constant of gravitationG is established. For the aphelion and perihelion points of the Earth's orbit we find, respectively,

Addendum to: Strength of the Solar Coronal Magnetic Field – A Comparison of Independent Estimates Using Contemporaneous Radio and White-Light Observations

This addendum uses an alternate fit for the electron density distribution \(N(r)\) (see Figure 1) and estimates the coronal magnetic field using the new model. We find that the estimates of the magnetic field are in close agreement using both the models.

We have fit the \(N(r)\) distribution obtained from STEREO-A/COR1 and SOHO/LASCO-C2 using a fifth-order polynomial (see Figure 1). The expression can be written as

$$\begin{aligned} N_{\text{cor}}(r) &= 1.43 \times 10^{9} r^{-5} - 1.91 \times 10^{9} r^{-4} + 1.07 \times 10^{9} r^{-3} - 2.87 \times 10^{8} r^{-2} \\ &\quad {} + 3.76 \times 10^{7} r^{-1} - 1.91 \times 10^{6} , \end{aligned}$$

(1)

where \(N_{\text{cor}}(r)\) is in units of cm^?3 and \(r\) is in units of \(\mathrm{R}_{\odot}\). The background coronal electron density is enhanced by a factor of 5.5 at 2.63 \(\mathrm{R}_{\odot}\) during the coronal mass ejection (CME). The estimated coronal magnetic field strength (\(B\)) using radio data indicates that \(B(r) \approx(0.51\text{\,--\,}0.48) \pm 0.02\ \mathrm{G}\) in the range \(r \approx2.65\text{\, --\,}2.82\ \mathrm{R}_{\odot}\). The field strengths for STEREO-A/COR1 and SOHO/LASCO-C2 are ≈?0.32 G at \(r \approx 3.11\ \mathrm{R}_{\odot}\) and ≈?0.12 G at \(r \approx 4.40\ \mathrm{R}_{\odot}\), respectively.

     

Identification of Emissions in the Spectrum of Comet C/1989 VI (Skoritchenko–George) Obtained with the 6-m SAO RAS Reflector

The emission spectrum of comet Skoritchenko–George (C/1989 VI), unusual in its information content, was obtained on February 26.7 UT, 1990, with the use of a TV scanner installed on the 6-m BTA reflector of the Special Astronomical Observatory of the Russian Academy of Sciences (SAO RAS) in Nizhni Arkhyz. Detailed identification of the emission lines of this comet was made. The observed spectrum contains 311 emission lines, including those of the molecules. Among others, the lines of the negative carbon C ₂ ^- ion and the lines corresponding to the electron transition in the neutral CO molecule are discovered. The presence of a large number of lines of the neutral CO molecule (the Asundi bands and the triplet bands) in the visible region is one of the uncommon features of the emission spectrum of this comet. The triplet lines : 15–3, 13–2, 11–2, 9–1, 8–1, 7–1, 7–0, 5–0, 4–0; : 7–0, 6–0, 5–0; and a" : 11–1 (K = 3, 4); 16–4 (K= 0, 1, 2, 4); 9-0 (K= 0, 1, 2); 8–0 (K= 0) were identified for the first time. Prior to this work, the lines of CO in the visible range were observed only in the spectrum of comet C/1979 VI (Bradfield) in 1989.     

Gamma-ray burst measurements with A 6.3 m2 array

Two flights from Alice Springs, Australia, were achieved in November 1977 and November 1978 with a plastic scintillator -burst detector, effective area 6.3 m², thickness 5 cm, energy response in the range 50 keV to 2 MeV. In 33 hr of good, high altitude data, two bursts were detected, yielding a rate corrected to an isotropic flux of at a size of 8.5×10^–9 erg cm^–2. One event, seen at 22.14 on 15 Nov 1978, was confirmed by spacecraft measurements. The second, too small to be detected by spacecraft, arrived from 0 hr RA, –13.2° Decl. ±12° and possibly comes from a confirmed -burst source location. A galactic origin with a source distribution originating from a relatively thick disk, is favoured by these results.     

On the periodic solutions of a particular Lame equation

We consider the Hill's equation: % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGGipm0dc9vqaqpepu0xbbG8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaSaaaeaaca% WGKbWaaWbaaSqabeaacaaIYaaaaOGaeqOVdGhabaGaamizaiaadsha% daahaaWcbeqaaiaaikdaaaaaaOGaey4kaSYaaSaaaeaacaWGTbGaai% ikaiaad2gacqGHRaWkcaaIXaGaaiykaaqaaiaaikdaaaGaam4qamaa% CaaaleqabaGaaGOmaaaakiaacIcacaWG0bGaaiykaiabe67a4jabg2% da9iaaicdaaaa!4973!\[\frac{{d^2 \xi }}{{dt^2 }} + \frac{{m(m + 1)}}{2}C^2 (t)\xi = 0\]Where C(t) = Cn (t, {frbuilt|1/2}) is the elliptic function of Jacobi and m a given real number. It is a particular case of theame equation. By the change of variable from t to defined by: % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGGipm0dc9vqaqpepu0xbbG8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaqcaawaaOWaaiqaaq% aabeqaamaalaaajaaybaGaamizaGGaaiab-z6agbqaaiaadsgacaWG% 0baaaiabg2da9OWaaOaaaKaaGfaacaGGOaqcKbaG-laaigdajaaycq% GHsislkmaaleaajeaybaGaaGymaaqaaiaaikdaaaqcaaMaaeiiaiaa% bohacaqGPbGaaeOBaOWaaWbaaKqaGfqabaGaaeOmaaaajaaycqWFMo% GrcqWFPaqkaKqaGfqaaaqcaawaaiab-z6agjab-HcaOiab-bdaWiab% -LcaPiab-1da9iab-bdaWaaakiaawUhaaaaa!51F5!\[\left\{ \begin{array}{l}\frac{{d\Phi }}{{dt}} = \sqrt {(1 - {\textstyle{1 \over 2}}{\rm{ sin}}^{\rm{2}} \Phi )} \\\Phi (0) = 0 \\\end{array} \right.\]it is transformed to the Ince equation: (1 + · cos(2)) y + b · sin(2) · y + (c + d · cos(2)) y = 0 where % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGGipm0dc9vqaqpepu0xbbG8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaqcaawaaiaadggacq% GH9aqpcqGHsislcaWGIbGaeyypa0JcdaWcgaqaaiaaigdaaeaacaaI% ZaGaaiilaiaabccacaWGJbGaeyypa0Jaamizaiabg2da9aaacaqGGa% WaaSaaaKaaGfaacaWGTbGaaiikaiaad2gacqGHRaWkcaaIXaGaaiyk% aaqaaiaaiodaaaaaaa!4777!\[a = - b = {1 \mathord{\left/{\vphantom {1 {3,{\rm{ }}c = d = }}} \right.\kern-\nulldelimiterspace} {3,{\rm{ }}c = d = }}{\rm{ }}\frac{{m(m + 1)}}{3}\]In the neighbourhood of the poles, we give the expression of the solutions.The periodic solutions of the Equation (1) correspond to the periodic solutions of the Equation (3). Magnus and Winkler give us a theory of their existence. By comparing these results to those of our study in the case of the Hill's equation, we can find the development in Fourier series of periodic solutions in function of the variable and deduce the development of solutions of (1) in function of C(t).     

Velocity-height relation for antimatter meteors

A general velocity-height relation for both antimatter and ordinary matter meteor is derived. This relation can be expressed as % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaSaaaeaacq% aHfpqDdaWgaaWcbaGaamOEaaqabaaakeaacqaHfpqDdaWgaaWcbaGa% eyOhIukabeaaaaGccqGH9aqpcaqGLbGaaeiEaiaabchacaqGGaWaam% WaaeaacqGHsisldaWcaaqaaiaadkeaaeaacaWGHbaaaiaabwgacaqG% 4bGaaeiCaiaabIcacaqGTaGaamyyaiaadQhacaGGPaaacaGLBbGaay% zxaaGaeyOeI0YaaSaaaeaacaWGdbaabaGaamOqaiabew8a1naaBaaa% leaacqGHEisPaeqaaaaakmaacmaabaGaaGymaiabgkHiTiaabwgaca% qG4bGaaeiCamaadmaabaGaeyOeI0YaaSaaaeaacaWGcbaabaGaamyy% aaaacaqGLbGaaeiEaiaabchacaqGOaGaaeylaiaadggacaWG6bGaai% ykaaGaay5waiaaw2faaaGaay5Eaiaaw2haaiaacYcaaaa!64FD!\[\frac{{\upsilon _z }}{{\upsilon _\infty }} = {\text{exp }}\left[ { - \frac{B}{a}{\text{exp( - }}az)} \right] - \frac{C}{{B\upsilon _\infty }}\left\{ {1 - {\text{exp}}\left[ { - \frac{B}{a}{\text{exp( - }}az)} \right]} \right\},\]where _z is the velocity of the meteoroid at height z, its velocity before entrance into the Earth's atmosphere, is the scale-height, and C parameter proportional to the atom-antiatom annihilation cross- section, which is experimentally unknown. The parameter B (B = DA₀/m) is the well known parameter for koinomatter (ordinary matter) meteors, D is the drag factor, ₀ is the air density at sea level, A is the cross sectional area of the meteoroid and m its mass.When the annihilation cross-section is zero — in the case of ordinary meteors — the parameter C is also zero and the above derived equation becomes % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaSaaaeaacq% aHfpqDdaWgaaWcbaGaamOEaaqabaaakeaacqaHfpqDdaWgaaWcbaGa% eyOhIukabeaaaaGccqGH9aqpcaqGLbGaaeiEaiaabchacaqGGaWaam% WaaeaacqGHsisldaWcaaqaaiaadkeaaeaacaWGHbaaaiaabwgacaqG% 4bGaaeiCaiaabIcacaqGTaGaamyyaiaadQhacaGGPaaacaGLBbGaay% zxaaGaaiilaaaa!4CF5!\[\frac{{\upsilon _z }}{{\upsilon _\infty }} = {\text{exp }}\left[ { - \frac{B}{a}{\text{exp( - }}az)} \right],\]which is the well known velocity-height relation for koinomatter meteors.In the case in which the Universe contains antimatter in compact solid structure, the velocity-height relation can be found useful.Work performed mainly at the Nuclear Physics Laboratory of the National University of Athens, Greece.     

On the number of isolating integrals in systems with three degrees of freedom

Dynamical systems with three degrees of freedom can be reduced to the study of a fourdimensional mapping. We consider here, as a model problem, the mapping given by the following equations: $$\left\{ \begin{gathered} x_1 = x_0 + a_1 {\text{ sin (}}x_0 {\text{ + }}y_0 {\text{)}} + b{\text{ sin (}}x_0 {\text{ + }}y_0 {\text{ + }}z_{\text{0}} {\text{ + }}t_{\text{0}} {\text{)}} \hfill \\ y_1 = x_0 {\text{ + }}y_0 \hfill \\ z_1 = z_0 + a_2 {\text{ sin (}}z_0 {\text{ + }}t_0 {\text{)}} + b{\text{ sin (}}x_0 {\text{ + }}y_0 {\text{ + }}z_{\text{0}} {\text{ + }}t_{\text{0}} {\text{) (mod 2}}\pi {\text{)}} \hfill \\ t_1 = z_0 {\text{ + }}t_0 \hfill \\ \end{gathered} \right.$$ We have found that as soon asb≠0, i.e. even for a very weak coupling, a dynamical system with three degrees of freedom has in general either two or zero isolating integrals (besides the usual energy integral).     

On the influence of gradients in the angular velocity on the solar meridional motions

If fluctuations in the density are neglected, the large-scale, axisymmetric azimuthal momentum equation for the solar convection zone (SCZ) contains only the velocity correlations and where u are the turbulent convective velocities and the brackets denote a large-scale average. The angular velocity, , and meridional motions are expanded in Legendre polynomials and in these expansions only the two leading terms are retained (for example, where is the polar angle). Per hemisphere, the meridional circulation is, in consequence, the superposition of two flows, characterized by one, and two cells in latitude respectively. Two equations can be derived from the azimuthal momentum equation. The first one expresses the conservation of angular momentum and essentially determines the stream function of the one-cell flow in terms of : the convective motions feed angular momentum to the inner regions of the SCZ and in the steady state a meridional flow must be present to remove this angular momentum. The second equation contains also the integral indicative of a transport of angular momentum towards the equator.With the help of a formalism developed earlier we evaluate, for solid body rotation, the velocity correlations and for several values of an arbitrary parameter, D, left unspecified by the theory. The most striking result of these calculations is the increase of with D. Next we calculate the turbulent viscosity coefficients defined by whereC _ro ⁰ and C _o ⁰ are the velocity correlations for solid body rotation. In these calculations it was assumed that ₂ was a linear function of r. The arbitrary parameter D was chosen so that the meridional flow vanishes at the surface for the rotation laws specified below. The coefficients v _ro ⁱ and v _0o ⁱ that allow for the calculation of C _ro and C _0o for any specified rotation law (with the proviso that ₂ be linear) are the turbulent viscosity coefficients. These coefficients comply well with intuitive expectations: v _ro ¹ and –v _0o ³ are the largest in each group, and v _0o ³ is negative.The equations for the meridional flow were first solved with ₀ and ₂ two linear functions of r ( ₀ ¹ = – 2 × 10 ^–12 cm ^–1) and ( ₂ ¹ = – 6 × 10 ¹² cm ^–1). The corresponding angular velocity increases slightly inwards at the poles and decreases at the equator in broad agreement with heliosismic observations. The computed meridional motions are far too large ( 150m s^–1). Reasonable values for the meridional motions can only be obtained if _o (and in consequence ), increase sharply with depth below the surface. The calculated meridional motion at the surface consists of a weak equatorward flow for gq < 29° and of a stronger poleward flow for > 29°.In the Sun, the Taylor-Proudman balance (the Coriolis force is balanced by the pressure gradient), must be altered to include the buoyancy force. The consequences of this modification are far reaching: is not required, now, to be constant along cylinders. Instead, the latitudinal dependence of the superadiabatic gradient is determined by the rotation law. For the above rotation laws, the corresponding latitudinal variations of the convective flux are of the order of 7% in the lower SCZ.     

A few observations of and remarks on V1500 Cygni

The development of the post-nova light curve of V1500 Cyg inUBV andHβ, for 15 nights in September and October 1975 are presented. We confirm previous reports that superimposed on the steady decline of the light curve are small amplitude cyclic variations. The times of maxima and minima are determined. These together with other published values yield the following ephemerides from JD 2 442 661 to JD 2 442 674: $$\begin{gathered} {\text{From}} 17 {\text{points:}} {\text{JD}}_{ \odot \min } = 2 442 661.4881 + 0_{^. }^{\text{d}} 140 91{\text{n}} \hfill \\ \pm 0.0027 \pm 0.000 05 \hfill \\ {\text{From}} 15 {\text{points:}} {\text{JD}}_{ \odot \max } = 2 442 661.5480 + 0_{^. }^{\text{d}} 140 89{\text{n}} \hfill \\ \pm 0.0046 \pm 0.0001 \hfill \\ \end{gathered} $$ with standard errors of the fits of ±0 _. ^d 0052 for the minima and ±0 _. ^d 0091 for the maxima. Assuming V1500 Cyg is similar to novae in M31, we foundr=750 pc and a pre-nova absolute photographic magnitude greater than 9.68.     

Intrinsic orbital coordinates on a Hamiltonian surface

A reversible dynamical system with two degrees-of-freedom is reduced to a second-order, Hamiltonian system under a change of independent variable. In certain circumstances, the reduced order system may be integrated following an orthogonal curvilinear transformation from Cartesian x,y to intrinsic orbital coordinates , . Solutions for the orbit position and true time variables are expressed by: % MathType!MTEF!2!1!+- % feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn % hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr % 4rNCHbGeaGGipm0dc9vqaqpepu0xbbG8F4rqqrFfpeea0xe9Lq-Jc9 % vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x % fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiEaiabg2 % da9iaadAgacaGGOaGaeqOVdGNaaiilaiabeE7aOjaacMcacaGGSaGa % aeiiaiaadMhacqGH9aqpcaWGNbGaaiikaiabe67a4jaacYcacqaH3o % aAcaGGPaGaaiilaiaabccacaWGKbGaamiDaiabg2da9iabgglaXoaa % dmaabaWaaSaaaeaacaWGibWaa0baaKqaahaacqaH+oaEaeaacaqGYa % aaaOGaam4raiabgUcaRiaadIeadaqhaaqcbaCaaiabeE7aObqaaiaa % ikdaaaGccaWGfbaabaGaaGOmaiaacIcacaWGibGaey4kaSIaamyvai % aacMcaaaaacaGLBbGaayzxaaWaaWbaaSqabKqaGhaacaaIXaGaai4l % aiaaikdaaaGccaWGKbGaeqiXdqhaaa!6498! \[ x = f(\xi ,\eta ),{\rm{ }}y = g(\xi ,\eta ),{\rm{ }}dt = \pm \left[ {\frac{{_\xi ^{\rm{2}} {\ie} + _\eta ^2 }}{{2( + U)}}} \righ \]1446 1040 where U is the potential function, and z is the new independent variable. The functions f, g may be expressed by quadratures when the metric coefficients {\er},{\ie} are specified. Two second-order, partial differential equations specify {\er}, {\ie} and Hamiltonian {\tH}. Auxiliary conditions are needed because the solutions are underdetermined. For example, both sets of curvilinear coordinate lines are orbits when certain dynamical compatibility conditions between U and {\ie} (or {\er}) are satisfied. Alternatively, when orbits cross the parametric curves, the auxiliary condition {\er} = {\ie} specifies a conformal transformation, and the partial differential equation for {\tH} may be reduced to an ordinary differential equation for the orbit curve. In either case, integrability is guaranteed for Lionville dynamical systems. Specific applications are presented to illustrate direct solution for the orbit (e.g., two fixed centers) and inverse solution for the potential.     

On pulsar-driven isothermal blast-wave models of supernova remanants

We investigate the ‘equilibrium’ and stability of spherically-symmetric self-similar isothermal blast waves with a continuous post-shock flow velocity expanding into medium whose density varies asr ^?ω ahead of the blast wave, and which are powered by a central source (a pulsar) whose power output varies with time ast ^ω?3. We show that:

1. for ω<0, no physically acceptable self-similar solution exists;
2. for ω>3, no solution exists since the mass swept up by the blast wave is infinite;
3. ? must exceed zero in order that the blast wave expand with time, but ?<2 in order that the central source injects a finite total energy into the blast wave;
4. for 3>ω_min(?)>ω>ω_max(?)>0, where $$\begin{gathered} \omega _{\min } (\varphi ){\text{ }} = {\text{ }}2[5{\text{ }} - {\text{ }}\varphi {\text{ }} + {\text{ }}(10{\text{ }} + {\text{ 4}}\varphi {\text{ }} - {\text{ 2}}\varphi ^2 )^{1/2} ]^2 [2{\text{ }} + {\text{ (10 }} + {\text{ 4}}\varphi {\text{ }} - {\text{ 2}}\varphi ^2 {\text{)}}^{{\text{1/2}}} ]^{ - 2} , \hfill \\ \omega _{\max } (\varphi ){\text{ }} = {\text{ }}2[5{\text{ }} - {\text{ }}\varphi {\text{ }} - {\text{ }}(10{\text{ }} + {\text{ 4}}\varphi {\text{ }} - {\text{ 2}}\varphi ^2 )^{1/2} ]^2 [2{\text{ }} - {\text{ (10 }} + {\text{ 4}}\varphi {\text{ }} - {\text{ 2}}\varphi ^2 {\text{)}}^{{\text{1/2}}} ]^{ - 2} , \hfill \\ \end{gathered} $$ two critical points exist in the flow velocity versus position plane. The physically acceptable solution must pass through the origin with zero flow speed and through the blast wave. It must also pass throughboth critical points if \(\varphi > \tfrac{5}{3}\) , while if \(\varphi< \tfrac{5}{3}\) it must by-pass both critical points. It is shown that such a solution exists but a proper connection at the lower critical point (for ?>5/3) (through whichall solutions pass with thesame slope) has not been established;
5. for 3>ω>ω_min(?) it is shown that the two critical points of (iv) disappear. However a new pair of critical points form. The physically acceptable solution passing with zero flow velocity through the origin and also passing through the blast wave mustby-pass both of the new critical points. It is shown that the solution does indeed do so;
6. for 3>ω_min(?)>ω_max(?)>ω it is shown that the dependence of the self-similar solution on either ω or ? is non-analytic and therefore, inferences drawn from any solutions obtained in ω>ω_max(?) (where the dependence of the solutionis analytic on ω and ?) are not valid when carried over into the domain 3>ω_min(?)>ω_max(?)>ω;
7. all of the physically acceptable self-similar solutions obtained in 3>ω>0 are unstable to short wavelength, small amplitude but nonself-similar radial velocity perturbations near the origin, with a growth which is a power law in time;
8. the physical self-similar solutions are globally unstable in a fully nonlinear sense to radial time-dependent flow patterns. In the limit of long times, the nonlinear growth is a power law in time for 5<ω+2?, logarithmic in time for 5>ω+2?, and the square of the logarithm in time for 5=ω+2?.

The results of (vii) and (viii) imply that the memory of the system to initial and boundary values does not decay as time progresses and so the system does not tend to a self-similar form. These results strongly suggest that the evolution of supernova remnants is not according to the self-similar form.     

Tidal effects in the earth–moon system and the earth's rotation

Analytical expressions for tidal torques induced by a tide‐arising planet which perturbs rotation of a nonrigid body are derived. Corresponding expressions both for secular and periodic perturbations of the Euler's angles are given for the case of the earth's rotation. Centennial secular rates of the nutation angle θ and of the earth's angular velocity ω, as well as the centennial logarithmic decrement ν of the Chandler wobble are evaluated:  mas, . In the Universal Time (UT) a large out‐of‐phase (sine) dissipative term with the period 18.6 years and the amplitude 2.3 ms is found. Corrections to nutation coefficients, which presumably have not been taken into account in IAU theory, are given. This revised version was published online in July 2006 with corrections to the Cover Date.     

Meridional flow in thick accretion disks

We have studied the effect of the flow in the accretion disk. The specific angular momentum of the disk is assumed to be constant and the polytropic relation is used. We have solved the structure of the disk and the flow patterns of the irrotational perfect fluid.As far as the obtained results are concerned, the flow does not affect the shape of the configuration in the bulk of the disk, although the flow velocity reaches even a half of the sound velocity at the inner edge of the disk. Therefore, in order to study accretion disk models with the moderate mass accretion rate—i.e.,

Non-singular de Sitter universe inR+R 2−2Λ gravity

For the theory described by the action and taking the FRW flat space metric we find an exact non-singular de Sitter model universe exp(t ²), with . It is also proved that the standard general relativity de Sitter cosmology , >0 is also a model of this higher derivative theory of gravity. If the metric is conformally flatS could describe a consistent quantum theory and its classical solutions would correspond to cosmological models in this theory.This work was supported in part by CONACYT grand P228CCOX891723, and DGICSA SEP grant C90-03-0347.     

A note on first-order normalizations of perturbed Keplerian systems

Techniques are developed to facilitate the transformation of a perturbed Keplerian system into Deláunay normal form at first order. The implicit dependence of the Hamiltonian on 1, the mean anomaly, through the explicit variable f, the true anomaly, or E, the eccentric anomaly, is removed through first order for terms of the form:

A note on theH-functions of transfer problems in multiplying media

Some useful results and remodelled representations ofH-functions corresponding to the dispersion function $$T\left( z \right) = 1 - 2z^2 \sum\limits_1^n {\int_0^{\lambda r} {Y_r } \left( x \right){\text{d}}x/\left( {z^2 - x^2 } \right)} $$ are derived, suitable to the case of a multiplying medium characterized by $$\gamma _0 = \sum\limits_1^n {\int_0^{\lambda r} {Y_r } \left( x \right){\text{d}}x > \tfrac{1}{2} \Rightarrow \xi = 1 - 2\gamma _0< 0} $$      

Exact solution of the equation of transfer in a finite exponential atmosphere by the method of laplace transform and linear singular operators

An approximate solution of the transfer equation for coherent scattering in stellar atmospheres with Planck's function as a nonlinear function of optical depth, viz.,

Apsidal motion in the eclipsing binary system of BW AQR

The first photoelectric light curve of the eclipsing binary system BW Aqr (F71V+F81V;P=6^d.7;V=10 ^m .31), discovered by Miss Leavitt at the beginning of the century, was obtained. The photometric elements were detemined. The components of this system are considerably evolved stars: the age of the system is about 2×10⁹ yr. It follows from the photometric data that the secondary component should have greater mass than the primary one The zero-age spectral classes of components were F2V and F1V. The system has an elliptical orbit with the eccentricitye=0.18. The angular rate of the apsidal motion (_obs = 0.070 deg yr^–1) and the corresponding value of the apsidal parameterk ₂=0.0090 (the relativistic term included) were found. The derived valuek ₂ exceeds by more than a factor of 2 the theoretical coefficient obtained from the modern theory of internal structure of stars with moderate masses .     

On numerical evaluation of theH-functions of transport problems by kernel approximation for the albedo 0<ω≤1

Das Gupta represented theH-functions of transport problems for the albedo [0, 1] in the formH(z)=R(z)–S(z) (see Das Gupta, 1977) whereR(z) is a rational function ofz andS(z) is regular on [–1, 0] ^c . In this paper we have representedS(z) through a Fredholm integral equation of the second kind with a symmetric real kernelL(y, z) as . The problem is then solved as an eigenvalue problem. The kernel is converted into a degenerate kernel through finite Taylor's expansion and the integral equation forS(z) takes the form: (which is solved by the usual procedure) where _r 's are the discrete eigenvalues andF _r 's the corresponding eigenfunctions of the real symmetric kernelL(y, z).