Fourier techniques for an analysis of eclipsing binary light curves: III

The aim of the present paper will be to present a new approach to the light changes of eclipsing binary systems. The light changes have been expanded into the Fourier-Bessel and Dini series. The coefficients of these expansions which are the Hankel transforms of the order of the loss of light (1–l) have been expressed in terms of the eclipse elements. These discrete Hankel transformsH (j _m ) andH ( _m ) valid for only the positive real zerosj _m and _m of the functionsJ and [xJ(x)+J (x)] have been generalized for any positive value of andy. Thus, these general expressions for the Hankel transforms of the light curves which are valid for all types of eclipses, for any arbitrary degree of the adopted limb-darkening law and, moreover, for any positive value of the free parameters andy, may be used for the solution of the elements of eclipsing binary systems.     

Further notes on the associated alpha-functions and related integrals in the theory of the light changes of eclipsing variables

With a view to furthering the theory of the light changes of eclipsing variables, developed before systematically by Z. Kopal, this paper presents a number of new (and computable) expressions for the associated alpha-function _n ⁰ ,(r₁,r₂,) (and also for its partial derivatives), where _n ⁰ ,(r₁,r₂,) represents the fractional loss of light suffered by an eclipse of a circular disc of fractional radiusr ₁ (and darkened at the limb to thenth degree) by an opaque disc of radiusr ₂, with their centres separated by a fractional (projected) distance , provided that the transparency of the occulting disc increases with the angle of foreshortening in the same manner as the limb-darkening of the eclipsed star (that is, when the transparency functiong(, ) of the second aperture is given by Equation (4) below). Many of the explicit expressions derived here are valid for any type of eclipse, occultation or transit, regardless of whetherr ₁>r ₂ orr ₁<r ₂. and for any degreen of the adopted law of limb-darkening. It is also pointed out how some of the results obtained in this paper are related to the various representations given earlier in the literature for the case =0.     

One-armed oscillations of a non-selfgravitating polytropic disk

One-armedglobal oscillations in a non-selfgravitating polytropic disk rotating around a star are investigated. The unperturbed disk is axisymmetric, geometrically thin, and extends infinitely in the radial direction keeping its thickness constant. Perturbations considered are inviscid and adiabatic. It is found that there are one-armed retrograde wave modes which are trapped in an inner region of the disk. The eignefrequency of the lowest order mode is given by _K(r _s)(z ₀/r _s)₂, wherer _s is the radius of the central star,z ₀ is the half-thickness of the disk, and _K(r _s) is the Keplerian angular frequency at the surface of the star.Paper presented at the IAU Third Asian-Pacific Regional Meeting, held in Kyoto, Japan, between 30 September–6 October, 1984.     

The polarization of coronal emission lines

By means of a photographic polarimeter, we attempted to measure both the amount and direction of linear polarization of all emission lines between 3400 and 9000 Å in the inner corona (1.034 r/r ₀ 1.085). Only the green and red coronal lines have been analyzed in detail. Neither of these lines shows polarization exceeding the probable error of 1.0% for 5303 and 1.8% for 6374. None of the other 17 coronal lines observed during the 7 March, 1970 solar eclipse show any obvious (>5%) polarization.Presently at the Sacramento Peak Observatory.     

The geometry of the Roche coordinates

The exact geometry of the Roche curvilinear coordinates (, , ) in which corresponds to the zero-velocity surfaces is investigated numerically in the plane, as well as in the spatial, case for various values of the mass-ratio between the two point-masses (m ₁,m ₂) constituting a binary system.The geometry of zero-velocity surfaces specified by -values at the Lagrangian points are first discussed by taking their intersections with various planes parallel to thexy-, xz- andyz-planes. The intersection of the zero-velocity surface specified by the -value at the Lagrangian equilateral-triangle pointsL _4,5 with the planex=1/2 discloses two invariable curves passing through the pointsL _4,5 and situated symmetrically with respect to thexy-plane whose form is independent of the mass-ratio.The geometry of the remaining two coordinates (, ) orthogonal to the zero-velocity surfaces is investigated in thexy- andxz-planes from extensive numerical integrations of differential equations generated from the orthogonality relations among the coordinates. The curves (x, y)=constant in thexy-plane are found to be separated into three families by definite envelopes acting as boundaries whose forms depend upon the mass-ratio only: the inner -constant curves associated with the masspointm ₁, the inner -constant curves associated with the mass-pointm ₂ and the outer -constant curves. All the -constant curves in thexy-plane coalesce at either of the Lagrangian equilateraltriangle pointsL _4,5, except for a limiting case coincident with thex-axis. The curves (x, z)=constant in thexz-plane are also separated by definite envelopes depending upon the mass-ratio into different families: the inner -constant curves associated with the mass-pointm ₁, the inner -constant curves associated with the mass-pointm ₂ and the outer -constant curves on both sides out of the envelopes. For larger values ofz, the curves =constant tend asymptotically to the line perpendicular to thex-axis and passing through the centre of mass of the system, except for a limiting case coincident with thex-axis. The geometrical aspects of the envelopes for the curves (x, y)=constant in thexy-plane and the curves (x, z)=constant in thexz-plane are also discussed independently.In the three-dimensional space, the Roche coordinates can be conveniently defined in such a way as to correspond to the polar coordinates in the immediate neighbourhood of the origin, and to the cylindrical coordinates at great distances. From numerical integrations of simultaneous differential equations generating spatial curves orthogonal to the zero-velocity surfaces, the surfaces (x, y, z)=constant and the surfaces (x, y, z)=constant are constructed as groups of such spatial curves with common values of some parameters specifying the respective surfaces.On leave of absence from the University of Tokyo as an Honorary Fellow of the Victoria University of Manchester.     

Comparison between some Hα and X-ray flares

In the present paper, H-evolutive curves of chromospheric events are compared with flux evolutive curves of X-ray events observed at the same time in different spectral regions. A correspondence between the emissions E(I _H/I _chr)'s at higher and higher H-intensity levels, and the X-ray fluxes F()'s in harder and harder -ranges is shown. Further, the present observations seem to indicate the existence of a single triggering mechanism during the flash-phase of a flare. It is also shown that these results may be in agreement with Brown's model for chromospheric flares.     

The Energy Equation in the Lower Solar Convection Zone

As a consequence of the Taylor–Proudman balance, a balance between the pressure, Coriolis and buoyancy forces in the radial and latitudinal momentum equations (that is expected to be amply satisfied in the lower solar convection zone), the superadiabatic gradient is determined by the rotation law and by an unspecified function of r, say, S(r), where r is the radial coordinate. If the rotation law and S(r) are known, then the solution of the energy equation, performed in this paper in the framework of the ML formalism, leads to a knowledge of the Reynolds stresses, convective fluxes, and meridional motions. The ML-formalism is an extension of the mixing length theory to rotating convection zones, and the calculations also involve the azimuthal momentum equation, from which an expression for the meridional motions in terms of the Reynolds stresses can be derived. The meridional motions are expanded as U _r(r,)=P ₂(cos)₂(r)/r ²+P ₄(cos)₄(r)/r ² +..., and a corresponding equation for U (r,). Here is the polar angle, is the density, and P ₂(cos), P ₄(cos) are Legendre polynomials. A good approximation to the meridional motion is obtained by setting ₄(r)=–H₂(r) with H–1.6, a constant. The value of ₂(r) is negative, i.e., the P ₂ flow rises at the equator and sinks at the poles. For the value of H obtained in the numerical calculations, the meridional motions have a narrow countercell at the poles, and the convective flux has a relative maximum at the poles, a minimum at mid latitudes and a larger maximum at the equator. Both results are in agreement with the observations.     

Fourier analysis of the light curves of eclipsing variables-XXV: Error analysis in the frequency-domain

The aim of the present paper has been two-fold. In the first part (Sections 1–2), closed algebraic formulae will be set up furnishing the momentsA of the light curves of arbitrary index , and, due to arbitrary type of eclipses, in terms of the coefficientsa of Fourier cosine series obtained by least-squares fit to the given data; and the uncertainty of the momentsA deduced from that of thea 's.In the second part (Sections 3–4) we shall establish the explicit forms of the lincar functions r _1,2, (cosi) and L ₁ for the variation of the respective elements expressible likewise in terms of the Fourier coefficientsa . The probable errors of these elements can then be identified with those of the respective linear functions, and are obtainable from the same matrix of coefficients which furnished the most probable values of the elements.     

Einstein-like field equations with conserved source and decreasing ∧ term; Their cosmological consequences

The consequences of a cosmological term varying asS ^–2 in a spatially isotropic universe with scale factorS and conserved matter tensor are investigated. One finds a perpetually expanding universe with positive and gravitational constantG that increases with time. The hard equation of state 3P>U (U mass-energy density,P scalar pressure) applied to the early universe leads to the expansion lawSt (t cosmic time) which solves the horizon problem with no need of inflation. Also the flatness problem is resolved without inflation. The model does not affect the well known predictions on the cosmic light elements abundance which come from standard big bang cosmology.In the present, matter dominated universe one findsdG/dt=2H/U (H is the Hubble parameter) which is consistent with observations provided <10^–57 cm^–2. Asymptotically (S) the term equalsGU/2, in agreement with other studies.     

Kelvin-Helmholtz instability of composite plasma in an oblique magnetic field with resistivity

Analytic structure of high-density steady isothermal spheres is discussed using the TOV equation of hydrostatic equilibrium which satisfies an equation of state of the kind:P = K _g , = _g c ².Approximate analytical solutions to the Tolman-Oppenheimer-Volkoff (TOV) equations of hydrostatic equilibrium in (, ), (,U) and (u, v) phase planes in concise and simple form useful for short computer programmes or on small calculator, have been given. In Figures 1, 2, and 3, respectively, we display the qualitative behaviours of the ratio of gas density _g to the central density _gc , _g / _gc ; pressureP to the _gc ,P/ _gc ; and the metric componente ^–, for three representative general relativistic (GR) isothermal configurations =0.1, 0.2, and 0.3. Figure 4 shows the solution curve (, ) for =0.1, 0.2, and 0.3 (=0 represents the classical (Newtonian) curve). Numerical values of physical quantitiesv (=4r ² P ^*(r)), in steps ofu (=M(r)/r)=0.03, and the mass functionU, in steps of =0.2 (dimensionless radial distance), are given, respectively, in Tables I and II. Other interesting features of the configurations, such as ratio of gravitational radius 2GM/c ² to the coordinate radiusR, mass distributionM(r)/M, pressure (or density) distributionP/P _c , binding energy (B.E.), etc., have also been incorporated in the text. It has further been shown that velocity of sound inside the configurations is always less than the velocity of light.Part of the work done at Azerbaijan State University, Baku, U.S.S.R., and Mosul University, Mosul, Iraq, 1985-1986     

Photometric and spectroscopic study of R CMa

UBV light curves and spectrograms of R CMa obtained with the 48-inch telescope of Japal-Rangapur Observatory during 1980–82 have been used for deriving the eclipse and orbital elements as well as the absolute dimensions of the components. The primary is found to be a Main-Sequence F2V star of mass 1.52M and the secondary a subgiant star of spectral type G8 and mass 0.20M which fills its Roche lobe, in agreement with Kopal and Shapley (1956) results, Kopal (1959), or Sahade's (1963) results. From a consideration of the possible evolution of this system it is concluded that a large fraction of the original mass of the secondary is lost from the system. A study of the period changes indicates the possible presence of a third component of mass of about 0.5M which is most likely to be anM dwarf.Paper presented at the Lembang-Bamberg IAU Colloquium No. 80 on Double Stars: Physical Properties and Generic Relations, held at Bandung, Indonesia, 3–7 June, 1983.     

Mean free paths and diffusion coefficients for energetic protons at small heliodistances calculated using Helios 1 and 2 data

Pitch angle scattering of energetic particles (100 MeV) in the interplanetary medium are studied using Helios 1 and 2 magnetometer and plasma data during 1976 near the minimum of solar activity. An IMF configuration was used in the computer experiments which allowed the pitch angle diffusion coefficient, D and hence the parallel mean free path, to be determined. The radial mean free path was found to vary as _r r ^-0.9 between 0.4 and 1 AU, but between 0.3 and 0.4 AU it decreases significantly. To reconcile our value of _r at 1 AU, lying between 0.01 and 0.02 AU, with the average prompt solar proton event profile, an increasing value of _r at lower radial distances would be required.     

Coronal magnetic-field model with non-spherical source surface

Previous global models of coronal magnetic fields have used a geometrical construction based on a spherical source surface because of requirements for computational speed. As a result they have had difficulty accounting for (a) the tendency of full magnetohydrodynamic (MHD) models to predict non-radial plasma flow out to r 10r and (b) the appreciable magnitude, 3, of B _r , (the radial component of B) consistently observed at r 1 AU. We present a new modelling technique based on a non-spherical source surface, which is taken to be an isogauss of the underlying potential field generated by currents in or below the photosphere. This modification of the source surface significantly improves the agreement between the geometrical construction and the MHD solution while retaining most of the computational ease provided by a spherical source surface. A detailed comparison between the present source-surface model and the MHD solution is made for the internal dipole case. The resulting B field agrees well in magnitude and direction with the coronal B field derived from the full MHD equations. It shows evidence of the slightly equatorward meridional plasma flow that is characteristic of the MHD solution. Moreover, the B field obtained by using our non-spherical source surface agrees well with that observed by spacecraft in the vicinity of the Earth's orbit. Applied to a solar dipole field with a moment of 1 G-r ³ , the present model predicts that B _r at r 1 AU lies in the range of 1–2 and is remarkably insensitive to heliomagnetic latitude. Our method should be applicable also to more general (i.e., more realistic) configurations of the solar magnetic field. Isogauss surfaces for two representative solar rotations, as calculated from expansions of observed photospheric magnetic-field data, are found to show large and significant deviations from sphericity.     

Observed oddities in the lines H,K, b and Hβ

We compare microphotometer intensity traces perpendicular to dispersion in simultaneous spectrograms of good spatial resolution traced at various 's in each of the lines. Cross correlations between the different traces show the following: (a) For each _K there is a corresponding _{b 1}at which the coefficient of correlation, r, is a maximum, usually > 0.8. (b) No such high correlations are found between H and H. (c) Comparison of traces in the continuum and at all observed 's in K, H, b₁, b₂ show a range of 's in each line over which r is very significantly negative, while H shows no such peculiarity.     

The distinctive feature of relativistic stellar systems

In this paper we adopt the method of relativistic fluid dynamics to examine the number density distribution of stars around a massive black hole in the core of stellar clusters. We obtain extensive results,n(r) r ^–a, 3/2a9/2, which include, respectively, then(r) r ^–7/4 power law obtained by Bahcall and Wolf and then(r) r ^–9/4 power law by Peebles. Sincen(r) is not an observable quantity for star clusters, we also consider general relativity effects, i.e., the consequence of the bending of light, in calculating the projected density of stars in such a system. As an example we employ a massive black hole 10³ M inlaid in the center of a globular cluster and calculate various projected densities of stars. The results show that cusp construction occurs in all cases unless the central black hole massM=0, and the polytropic index does not affect at all the position of the capture radiusr _a. The obvious differences in the surface density is only embodied in the interior of the capture radius. At the outer regions of the core, the surface density of stars declines rapidly with ar ^–5 power law in all cases. These results can be applied to cases of unequal-mass and non-steady state.     

Lower hybrid wave model for aurora

In the framework of non-linear fluid theory we use a lower hybrid (LH) wave of the form as a pump which interacts with the small fluctuations with the low-frequency vibrations _i or =0, where _i , is the hydrogen ion-cyclotron (HIC) gyrofrequency. The ponderomotive force generated by the beating of the high-frequency pump wave ₀ and the sideband LH waves (±₀) produces a non-linear coupling between the high- and low-frequency motions of electrons and ions. Under certain conditions the HIC waves and the zero-frequency waves both become parametrically unstable and start to grow. These excited waves then heat the ions by stochastic acceleration in the transverse direction, thus explaining the formation of ion comics along the auroral field lines. Electrons would be heated in the parallel direction directly by the pump field as well as by low-frequency waves. Thus a single mechanism can explain the existence of ion-cyclotron waves, zero-frequency waves, ion conics, and energetic electrons along the auroral field lines.     

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     

Characteristics of electron and high-energy proton flares

High-energy proton (E _p > 55 MeV) and electron (E _e > 50 keV) events were observed by University of Iowa experiments on the satellites Explorer 33 and 35. The solar X-ray (2–12 Å) flares associated with the energetic proton events were found to have in general higher peak fluxes, considerably longer decay times (t) and smaller rise to decay time ratios (r) than the X-ray flares associated with the electron events. The most common decay times and rise to decay time ratios are: 80 t 100 min, 0.1 r 0.2 for the proton X-ray flares and t 20 min, 0.3 r 0.7 for the electron ones.     

Analysis of red and infrared wide-band photometry of Algol

The basic physical picture of the Algol system is reviewed, and, using collected red and infrared observations, photometric curve fits are investigated by applying numerical quadratures to determine theoretical light curves appropriate to Roche model stars. The contact nature of Algol B appears to be confirmed, and effective temperatures of the three components areT _A=11 400°,T _B=5300°,T _C=7600°. In terms of a Lambert's law approach to the reflection effect, the effective heat-albedo is required to be reduced from unity to one half; and it is also found that the averaged gravity-darkening coefficient is close to a value appropriate for a diffusion type of heat-transfer mechanism operating in sub-photospheric layers.     

Analysis of the light curves of RT Persei by automatized fourier techniques

A new method has been developed by Kopal (1977c, Paper XII) to make use of expressions for the fractional loss of light ₀ ^l of the arbitrarily limb-darkened stars in the form of Hankel transforms of zero-order, in order to evaluate the explicit forms of the ₀ ^l s for different types of eclipse, as well as of the momentsA _2m of the respective light curves in a closed form. The automated method has been tested successfully on the light curves of RT Persei. Also, a photometric curve fit of RT Persei is investigated by application of numerical quadratures to determine the theoretical light curve appropriate for the Roche model. Finally a comparative discussion is given of various methods of light curve analaysis.On leave from Department of Physics, University of Ferdowsi, Mashad, Iran.