The light curved in the CM field

In this paper we introduce the CM field in Sections 2 and 3 based on the paper by Wang and Peng (1985), and calculate the light curved in the CM field in Section 4. The result shows thatP makes _CM larger than _C at , and smaller at . Under a special circumstance which source, CM lens, and observer are in the same line, if we get | ₀₌₀ , and | _=/2 , we can determine theP(M) andQ(M) of the CM lens,M is the mass of the CM lens.     

Static spherical configurations of cold matter in the Einstein-Cartan theory of gravitation

We investigate static, spherical configurations of cold catalized matter in the Einstein-Cartan theory of gravitation. Assuming that density of spin is proportional to the number density of baryonsn and using an equation of state of a degenerate, relativistic Fermi gas, we numerically integrated the relativistic equation of equilibrium. We have also studied the stability of those configurations. Configurations with central number densityn _c such that where is the effective pressure, are very similar to general relativistic configurations with the same central density. In the Einstein-Cartan theory there exists another disjoint family of equilibrium configurations for which but . Those configurations have very small masses 10^–6 g and raddi 10^–34 cm and are unstable.Supported in part by Research Grant MR-I-7.     

A new technique for the determination of coronal magnetic fields: A fixed mesh solution to Laplace's equation using line-of-sight boundary conditions

A new method for computing potential magnetic field configurations in the solar atmosphere is described. A discrete approximation to Laplace's equation is solved in the domain R r R ₁, 0 , 0 2 (R ₁being an arbitrary radial distance from the solar center). The method utilizes the measured line-of-sight magnetic fields directly as the boundary condition at the solar surface and constrains the field to become radial at the outer boundary, R ₁. First the differential equation and boundary conditions are reduced to a set of two-dimensional equations in r, by Fourier transforming out the periodic dependence. Next each transformed boundary condition is converted to a Dirichlet surface condition. Then each two-dimensional equation with standard Dirichlet-Dirichlet boundary conditions is solved for the Fourier coefficient it determines. Finally, the solution of the original three dimensional equation is obtained through inverse Fourier transformation. The primary numerical tools in this technique are the use of a finite fast Fourier transform technique and also a generalized cyclic reduction algorithm developed at NCAR. Any extraneous monopole component present in the data can be removed if so desired.The code was developed for the HAO solar-interplanetary modeling effort in response to the following specific requirements:

Charged particle diffusion in the presence of Alfvén waves: The perpendicular particle flux

The diffusion of charged particles through a weak stochastic electro-magnetic field which is superimposed on a constant background magnetic field is considered. The stochastic electromagnetic fields are assumed to consist of unpolarized Alfvén waves propagating at arbitrary angles to the direction . When the Alfvén waves are propagating in directions other than and the particle gyro-radius,r _g, is sufficiently large (but may be smaller than the correlation length of the stochastic fields) it is shown that the particle flux perpendicular to the direction is , wherev is the particle speed andf the particle density. The expression forK differs from those calculated by previous authors. For small particle gyro-radii the flux S has a different functional form and is identical to that found by Urch (1977) to describe particle diffusion when the Alfvén waves only propagate in the direction .     

Invariant properties of families of moon-to-earth trajectories

Analytical techniques are employed to demonstrate certain invariant properties of families of moon-to-earth trajectories. The analytical expressions which demonstrate these properties have been derived from an earlier analytical solution of the restricted three-body problem which was developed by the method of matched asymptotic expansions. These expressions are given explicitly to orderµ ^1/2 where is the dimensionless mass of the moon. It is also shown that the inclusion of higher order corrections does not affect the nature of the invariant properties but only increases the accuracy of the analytic expressions.The results are compared with the work of Hoelker, Braud, and Herring who first discovered invariant properties of earth-to-moon trajectories by exact numerical integration of the equations of motion. (Similar properties for moon-to-earth trajectories follow from the principle of reflection). In each instance the analytical expressions result in properties which are equivalent, to orderµ ^1/2, with those found by numerical integration. Some quantitative comparisons are presented which show the analytical expressions to be quite accurate for calculating particular geometrical characteristics.

Nomenclature

Orbital Elements near the Moon energy - angular momentum - semi-major axis - eccentricity - inclination - argument of node - argument of pericynthion Orbital Elements near the Earth h _e energy - l _e angular momentum - i inclination - argument of node - argument of perigee - t _f time of flight Other symbols parameters used in matehing - U a function of the energy near the earth - a function of the angular momentum near the earth - r _p perigee radius - perincynthion radius - radius at node near moon - true anomaly of node near moon - initial angle between node near moon and earth-moon line - a function ofU, , andi - earth phase angle - dimensionless mass of the moon - U ₀, U₁ U=U ₀+U ₁ - i ₀, i_1/2, i₁ i=i ₀+µ ^1/2 i _1/2+µ i ₁ - ₀, _1/2, ₁ = ₀+µ ^1/2 i _1/2+µ i ₁ - _p longitude of vertex line - _n latitude of vertex line - R _o ,S _o ,N _o functions ofU ₀ and - a function ofU ₀, and      

The critical and the saturation content of magnetic monopoles in rotating relativistic objects

Both the critical content _c ( N _m /N _B , whereN _m ,N _B are the total numbers of monopoles and nucleons, respectively, contained in the object), and the saturation content _s of monopoles in a rotating relativistic object are found in this paper. The results are:

Two Sets of Improved Approximate Expressions of the Gyrosynchrotron Radiation

In this paper two sets of improved approximate expressions of emissivity , absorptivity , effective temperature T_eff, and frequency of peak brightness _p of gyrosynchrotron radiation are presented respectively for the ranges from 5 to 10 and 10 to 100 of harmonic numbers s(= /_B). The expressions are designed for the range from 20° to 80° of viewing angle , and the range 2 to 7 of electron energy spectral index . They are expressed by a power-law function in which the indexes are fitted by polynomial expressions of . Their statistical errors are, respectively, 24% and 32% for and for and 28% for . Their accuracies are much better than those of linear fitting of the power-law index.     

The First DIRECT Distance to a Detached Eclipsing Binary in M33

We present the first direct distance determination to a detached eclipsing binary in M33, which was found by the DIRECT Project. Located in the OB 66 association, it was one of the most suitable detached eclipsing binaries found by DIRECT for distance determination, given its 4.8938 day period. We obtained follow-up BV photometry and spectroscopy from which we determined the parameters of the system. It contains two O7 main sequence stars with masses of and and radii of and , respectively. We derive temperatures of K and K and determine the reddening . Using HST photometry for flux calibration in the V band, we obtain a preliminary distance modulus of mag ( kpc). The photometry and thus distance is subject to revision in the final paper.     

Inertia coefficient considerations and the structure of Jovian planets

For Jupiter, an overall density model of the form= ₀(1–x ⁿ ), withn1/3 and , is consistent with information presently at hand; for Saturn, however, such a density law would lead to unacceptably high densities in the vicinity of the centre. The limiting cases of the previous law are shown to ben=+, corresponding to a homogeneous sphere, andn=–3, corresponding to a particular central particle model, investigated by a number of astronomers over the last hundred years. Forn0, the central density becomes +. Another possible representation, valid both for Jupiter and Saturn, is the density law= ₀(1–x) ^m ), with in the case of Jupiter, and in the case of Saturn. Graber's density law based on a maximum entropy principle leads to unacceptably high surface densities, both for Jupiter and Saturn. Finally, the paper investigates the problems involved in fitting two-layered parametrically simple density laws to theoretically derived much more elaborate models of the Jovian planets.     

3He enrichments by magnetosonic waves in current-driven instabilities

A theoretical model for³He enrichments in solar energetic particles is developed. First, current-driven, electrostatic instabilities that have frequencies ( is the cyclotron frequency of³He) are investigated for a plasma consisting of H,⁴He,³He, and electrons with the density of³He much lower than those of H and⁴He. It is found that in many cases the oblique ion-acoustic waves can have positive growth rates at frequencies and, at the same time, negative growth rates at and at _H. This can occur near the marginal state of the instability. The wave damping at these frequencies is caused by the cyclotron resonances of⁴He and H. The cyclotron damping at is negligible, however, because the abundance of³He is very small. The H cyclotron waves can be unstable at for a wide region of plasma parameters; the electron-to-ion temperature ratio must beT _e /T _H 1.5. To destabilize the₄He cyclotron waves with , high⁴He density and high electron temperature are both required. Then,³He enrichments are studied on the basis of the theory of nonlinear magnetosonic waves, which can promptly accelerate ions. The current-driven electrostatic waves with can enhance fluctuation velocities of³He. Thus, in the presence of these waves, magnetosonic waves can selectively accelerate³He particles to high energies. Finally, cyclotron resonances of heavy ions with the waves or are briefly discussed.     

Epstein weight functions for non-radial oscillations of polytropes

Weight functions for the determination of the periods of linear adiabatic non-radial oscillations have been calculated in the same manner as Epstein's classic treatment of purely radial oscillations. Quadrupole (l=2) oscillations for thef and lower orderp andg-modes were considered. One group of static models were polytropes in the range 1.0n4.0 with ; thus included were configurations that were convectively stable, unstable and neutrally stable throughout. Another group consisted ofn=3.0 polytropes with convective shells or convective cores; ₁ was set at different values in each region in order to produce stability ( ) or instability ( ). The weight function provides a pictorial means for assessing the relative importance of each region of a given static model with respect to generating a given non-radial mode.     

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.     

On the uniqueness of the determination of the coronal potential magnetic field from line-of-sight boundary conditions

We consider a simple model in which the coronal magnetic field B is assumed to be potential in the region between the solar surface _o and an exterior source-surface ₁ of arbitrary shape. We prove that the boundary value problem that determines B from the value B _lof its component on ₀ along either (orthoradial direction) or (fixed direction) has at most one solution. On the other hand, we show that a solution can exist only if B _lsatisfies some solubility conditions.     

Graphisches Konstruktionsschema zur Bestimmung der Bewegungsparameter eines drallstabilisierten Flugkörpers mit Nutationsdämpfung unter dem Einflusz von Lagekorrekturimpulsen

Zusammenfassung Es wird gezeigt, daß die unter der Einwirkung einer Momentenimpulsserie entstehende Bewegung eines rotierenden Flugkörpers mit Nutationsdämpfung sich vollständig einem regelmäßigen Polygon entnehmen läßt, das durch das Trägheitsmomentenverhältnis, den Integralwert eines Einzelimpulses, den Drall und eine die Dämpfung charakterisierende KonstanteK ₀ bestimmt ist.Die Bewegung setzt sich aus logarithmischen Spiralen zusammen, derenn-ten Anfangsradius man erhält, indem man den Teilungspunkt des im VerhältnisK ₀:1 geteilten (n–1)-ten Radius mit der (n+1)-ten Polygonecke verbindet.Es wird bewiesen, daß das Konstruktionsnetz zu einem im äußeren Polygon liegenden ähnlichen inneren Polygon konvergiert, das gegenüber ersterem gedreht ist.Einfache Beziehungen zur Bewegungsbestimmung mit dem Polygonschema werden für Pulsfrequenzen angegeben, die ganzzahlige Vielfache oder Bruchteile der Spinfrequenz sind.

---

It is shown that the motion of a spinning body with nutation damping due to a series of torque pulses can be completely derived from a regular polygon determined by the ratio of inertias, the integral of one pulse, the momentum and a constantK ₀ characterizing damping.The motion is composed of spirals thenth initial radius of which is obtained by connecting the dividing point of the (n–1)th radius with the (n+1)th polygon corner. Each dividing point divides the respective radius in the ratioK ₀:1. The net of construction lines converges into an inner polygon turned against the outer one and having the same shape.Simple rules are shown for the application of the scheme on pulse frequencies which are multiples or fractions of spin frequency.

---

Symbole 1-2-3 Achsen des flugkörperfesten Koordinatensystems - a,b,c Hilfsgrößen zur Bestimmung der Iterationsgrößen - E _i i-te Polygonecke - H Drall des Flugkörpers - K _i Verhältnis deri-ten Drehzeigerlängen zu Beginn und am Ende eines Impulses - M Iterationsmatrix - Integralwert des Momentenimpulses - P ₀ Äußeres Polygon - P ₁ Spitze des Drehzeigersr _00e - P Drehpunkt des Drehzeigersr ₀₀ - P Konvergierendes Polygon - P _i Teilungspunkt des [i–1]-ten Zeigers - r _0i Drehzeiger aufgrund desi-ten Impulses allein - r _0ia Zeigerr _0i in Anfangslage - r _0ie Zeigerr _0i in Endlage - r _i i-ter Summenzeiger - r _ia Zeigerr _i in Anfangslage - r _ie Zeigerr _i in Endlage - T Dauer einer Flugkörperumdrehung - t,t, Zeitargumente - x-y-z Achsen eines raumfesten Koordinatensystems - x _i ,y _i Iterationskoordinaten - _n Phase desn-ten Radius gegenüber der anliegenden Polygonseite - Drehung des inneren Polygons gegenüber dem äußeren - Abklingkonstante - Phasenänderung des Drehzeigers innerhalb einer Flugkörperumdrehung - ₀ Anteil der über 2 hinausgehenden Phasenänderung des Drehzeigers - ₃ Trägheitsmoment um die Spinachse - ₁₂ Trägheitsmoment um die Querachsen - Zahl der Ecken des Konstruktionspolygons - _1,2 Eigenwerte der Iterationsmatrix - Zahl der vollen Umläufe des Konstruktionspolygons - Fortbewegungsachse des Drallvektors - ₀ Ausgangsphasenwinkel - _i Phasenlage desi-ten Summenzeigers - _x, _y Drehwinkel nach Einzelimpuls fürt - , Funktionen der Iterationsgrößen - , Drehwinkel umx-bzw.y-Achse - Drehgeschwindigkeit der Spinachse um den Drallvektor - Fiktive Größen bei Pulsfrequenzen kleiner als Spinfrequenz - Fiktive Größen bei Pulsfrequenzen größer als Spinfrequenz     

Der offene sternhaufen NGC 2362

The very young open star cluster NGC 2362 was investigated by the strip method on charts of two photographs taken with the 1-m Schmidt telescope of the European Southern Observatory. Up to the limiting magnitudeM _v ^* =5 _. ^m 8 the cluster contains 100 stars and can be described by the Gaussian density law (6). Further results are: Mass = 246 , central mass density ₀ = 43.1 = 246 pc^-3 , radiusR2.6 pc, mean velocity of the stars = 0.64 km s^–1.

---

---

Mitteilungen Serie A.     

A Reconnection Model for the Distribution of Flare Energies

A physically based explanation is given for the distribution of flare energies N(E)E ^– where 1.5. In contrast to previous approaches, the present treatment is based on a physical theory of the flare reconnection site. The central assumption is that topological flare energy, although released explosively, is slowly accumulated over several hundred Alfvén timescales. When coupled to the geometric properties of the reconnective flare source, this assumption is shown to lead naturally to a deduction of the flare energy distribution. Current sheet models yield the exponent whereas more compact current structures imply steeper spectra .     

On the peculiarities in the rotational frequency evolution of isolated neutron stars

The measurements of pulsar frequency second derivatives have shown that they are 10²−10⁶ times larger than expected for standard pulsar spin-down law, and are even negative for about half of pulsars. We explain these paradoxical results on the basis of the statistical analysis of the rotational parameters ν, and of the subset of 295 pulsars taken mostly from the ATNF database. We have found a strong correlation between and for both and , as well as between ν and . We interpret these dependencies as evolutionary ones due to being nearly proportional to the pulsars’ age. The derived statistical relations as well as “anomalous” values of are well described by assuming the long-time variations of the spin-down rate. The pulsar frequency evolution, therefore, consists of secular change of ν _ev(t), and according to the power law with n≈5, the irregularities, observed within a timespan as a timing noise, and the variations on the timescale larger than that—several decades. This work has been supported by the Russian Foundation for Basic Research (grant No 04-02-17555), Russian Academy of Sciences (program “Evolution of Stars and Galaxies”), and by the Russian Science Support Foundation. The authors would also like to thank the anonymous referee for valuable comments.     

Compatibility conditions for a non-quadratic integral of motion

Conditions are found which are satisfied by the coefficients of the expression being a second integral of the motion of an autonomous dynamical system with two degrees of freedom. The coefficientsA, B. , ,E are differentiable functions of the cartesian position coordinatesx, y. The velocity components are denoted by . It is shown that must be constant andB must be of the formB =f(x+y) +g(x-y) wheref, g are arbitrary.Given andB one can always find the remaining coefficientsA, E and also the corresponding potential and second integral. Depending on the specifica case at hand a certain number of arbitrary constants (or arbitrary functions) enter into the potential and the second integral. To each potential (which may be of the separable or nonseparable type in the coordinatesx andy)there corresponds one integral of the above form.     

Distribution of Flare Energies Based on Independent Reconnecting Structures

A model is presented to explain the observed frequency distribution of flare energies, based on independent flaring at a number of distinct topological structures (separators) within active-region magnetic fields. The model is a modification and generalization of a recent model due to Craig (2001), and reconciles that model with the observed flare waiting-time distribution, and the observed absence of a flare waiting-time versus energy relationship. The basic assumptions of the model are that flares of energy E ² occur at separators of length , and that the frequency of flaring at a separator is defined by the Alfvén transit time of the structure. To reproduce the observed distribution of flare energies the model requires a probability distribution P( ) ^–1 of separator lengths within active regions. This prediction of the model is in principle testable. A theoretical origin for this distribution is also discussed.