Exact solution to radiation-filled Brans-Dicke closed space cosmology

Under the assumption of a power law (k·R ⁿ=C,C=const.) between the gravitational constantk and the radius of curvatureR of the Universe and forP=1/3 the exact solution is sought for the cosmological equations of Brans and Dicke. The solution turns out to be valid for closed space and the parameter of the scalar-tensor theory is necessarily negative. The radius of curvature increases linearly with respect to the age of the Universe while the gravitational constant grows with the square of the radius of curvature. It has been shown (Lessner, 1974) that in this case (KR ²) the spatial component of the field equations is independent of the remaining equations. However, our solution satisfies this independent equation. This solution for the radiation-dominated era corresponds to the solution for the matter-dominated era found by Dehnen and one of the authors (Dehnen and Obregón, 1971). Our solution, as is the solution previously obtained for the matter-dominated era, is in contradiction to Dirac's hypothesis in which the gravitational constant should decrease with time in an expanding Universe.     

Comments on the BDT-FRW perfect fluid solutions

It is shown that the BDT-FRW perfect fluid solutions given some years ago by Chauvet and Obregón (1979) are identical to the solutions first given by Jordan (1955) and later by Brill (1962).     

Exact cosmological solutions in Brans and Dicke's scalar-tensor theory,II

In addition to our previous paper (Dehnen and Obregón, 1971) the exact cosmological solutions of Brans and Dicke's scalar-tensor theory allowing a power law between the gravitational constant and the radius of curvatureR of the universe are sought in case that — in contrast to our previous paper — the initial conditionR(t=0)=0 is avoided. There exist two different types of solutions especially for the closed space of positive curvature and for positive values of the freely available parameter of the scalar-tensor theory. The radius of curvature and also the gravitational constant increase at first with respect to time and decrease after reaching a maximum value, in contradiction to Dirac's hypothesis whereafter the gravitational constant should decrease with time in an expanding universe. The age of the universe following from these solutions is in accordance with the observations.     

(Super) cosmological solution for the Kasner model

In the theory of supergravity (N=1), the supersymmetric version of general relativity, and for the Kasner cosmological model (Bianchi type I) we find a non-trivial solution (for the metric and spinor-vector) under the most simple assumption =₁₁ + ₂₂; ₁₂+₂₁=0 and for a special choosed gaugeN=1,N ^j=0, ₀=0. This method could be also applied to other cosmological metrics and extended to enlarged Grassmann basis.O. Obregón was partially supported by the Alexander von Humboldt Stiftung.     

Comment on the LRS-BDT-Bianchi type-V vacuum solution

We wish to point out that the anisotropic LRS-BDT-Bianchi type-V vacuum solution recently given by Ram and Singh is wrong. Moreover, the correct solution is nothing but an isotropic FRW (k=–1) solution given by us previously.     

The rate of growth of density perturbations in simple big-bang model universes

In this paper we present exact solutions to the density perturbation equation derived by Bonnor for the cases where = _c ,k=1 and = – _c ,k=–1. The solutions when =0,k=1 and =0,k=–1 have been previously published. Using these solutions a quantitative analysis has been carried out that has enabled us to estimate the size of the fluctuations that must be postulated at decoupling in order to explain the formation of the galaxies in these model universes.     

Some special non-flat BDT-FRW solutions with a cosmological constant

We consider the non-flatk ²=1 Friedmann-Robertson-Walker models with a non-vanishing cosmological constant in the Brans-Dicke theory of gravitation. Some special solutions are given in case of =0, which may be considered as exact solutions of the projective theory of Jordan and Thiry.     

The application of a non-static spherically-symmetric metric to clusters of galaxies

An approximate metric is found which represents a sphere of matter embedded in a background of dust. The use of this metric in conjunction with the Friedmann equations gives values of for the three possible values ofk as +6×10^–36 (k=+1), +3×10^–35 (k=0), +10^–36 (k=–1). These values depend on data regarding clusters of galaxies, and are probably accurate to within an order of magnitude given the correctness of the assumptions on which their derivation rests.     

Analysis of Neckel and Chini standard stars in theUBVRI photometric system

This paper studies the list of 60 faint standard stars of Neckel and Chini (1980) in theUBVRI Johnson photometric system. These standards have been used due to the collaboration in the HIPPARCOS project of our department, but this work was made in parallel of this and the results were not used.We have re-reduced the standard stars using our own method of reduction from a large number of selected observations made at the Observatorio Astronómico Nacional and at the Centro Astronómico Hispano-Alemán, both at Calar Alto, Spain, and at the Observatorio del Roque de los Muchachos, La Palma Island, Spain.The reduction method used is explained and the results are given for the four colours used: (B-V),(U-B), (V-R), and (V-I), together with theV magnitude, that is, we have corrected many standard stars for each colour and magnitude. These results are supported with graphs where residual colour (o magnitude), is plotted against corresponding colour (or magnitude).Finally, our (V-R) and (V-I) results are compared to those of Tayloret al. (1989), with a good correspondence.Observations were made at Observatorio Astronómico Nacional and Centro Astronómico Hispano-Alemán, at Calar Alto, and in the Observatorio del Roque de los Muchachos at La Palma. This work was supported, in part, by the Comisión Interministerial de Ciencia y Tecnologia of Spain, under contract PB85-0017.     

Third and fourth order resonances in Hamiltonian systems

The stability of the origin of an autonomous Hamiltonian system is investigated when the system possesses a third or fourth-order resonance.H ₂, the quadratic part ofH isH ₂=ⁿ _{i=1 i} J _i and the resonance condition is ⁿ _i=1 k _{i i} where thek 0,i = 1, 2, ...,n are the natural or fundamental frequencies. It is shown that the only case in which the origin can be unstable is ifk _i0,i=1,2,...,n. The condition for instability is then given in terms of the coefficients of the higher order terms in the Hamiltonian. The transfer of energy between modes is also investigated when a near-resonant condition exists.     

A self-consistent world model

The field equations of the generalized field theory constructed by Mikhail and Wanas have been applied to a well-established geometrical structure given earlier by H. P. Robertson in connection with the cosmological problem. A unique solution, representing a specified expanding Universe (withq ₀=0, ₀=0.75,k=–1) has been obtained. The model obtained has been compared with cosmological observations and with FRW-models of relativistic cosmology. It has been shown that the suggested model is free of particle horizons. The existence of singularities has been discussed.The two cases, when the associated Riemannian-space has a definite or indefinite metric are considered. The case of indefinite metric with signature (+ – – –) is found to be characterized byk=–1, while the case of +ve definite metric is characterized byk=+1. Apart from that difference, the two cases give rise to the same cosmological parameters. It has been shown that energy conditions are satisfied by the material contents in both cases.     

Macro-and micro-turbulent filter functions for weak lines in stellar atmospheres

Following a similar discussion given earlier for the solar case (De Jager, 1972) we compute in this paper spectral line profiles for the spatial wavelengths in which a stellar motion field can be decomposed, and thereafter the macro-and micro-turbulent filter functionsf _M(k) andf (k), where is the optical scale height andf ²(k) dk the fraction of the energy of the turbulent motions between wavenumbersk andk+dk of the spectrum of turbulence that contributes to either kind of turbulence. If micro-and macro-turbulent velocity components are known for a certain star, and if the spectrum of turbulence is sharp enough, the ratiof _M/f would enable one to derive the average size of the turbulent elements in the star's atmosphere. The computations apply to weak lines in idealized stellar atmospheres, and refer to two cases: isotropic turbulence, and radial pulsations. These filters can be suitably used in a diagnostic method for the analysis of the motion field in the solar and stellar atmospheres. Some examples of applications to stars of very different kinds are given.     

Apsidal motion in the eclipsing binary system V 451 Ophiuchi

MulticolourWBVR photoelectric observations of the eclipsing binary V 451 Oph were carried out, and a highly accurate light curve was obtained. The angular velocity of the orbital rotation, =2.1 deg yr^–1, and the apsidal motion constantk ₂=0.0045 are given.     

Limits of stability for an area-preserving polynomial mapping

The Henon-Heiles mappingx=x+a(y–y ³), y=y–a(x–x³) has been studied, with the aim of finding where the unstable regions of the (x, y) plane are. When this mapping is put into the normal form, it is found to be a typical twist mapping. The criteria of Moser (1971) are used to obtain an upper limit to the size of a stable region around the origin, and this limit decreases to zero as the value of the parameter a increases toward 2.0. However, direct calculation fora=1.99 shows that there is a fairly large region insidex=0.412,y=0, from which escape from near the outer boundary requires at least 160 mappings. The region of high stability thus appears to be much larger than any region of absolute stability predicted by the KAM theorem.A general survey has been made of instability regions for the parameter valuea=1.0, this survey having been carried out to the extent which is allowed by a computer with 18-decimal-place accuracy. First, for all thex-axis fixed points (of the above mapping) deemed to be representative and significant, both the locations and variational matrix traces have been calculated. (The latter show whether the fixed point is elliptic or hyperbolic.) Ifn is the number of mappings andk is the number of circuits around the origin, then the listing (Table IV) is for fractionsk/n between 1/6 and 1/22, inclusive. (This covers the range 0x<0.96, withx=0 the fixed point forn=6,k=1).Escape toward infinity can be rapid, with less than 200 mappings necessary to reach the vicinity of then=1 fixed points (atx=±1,y=0 andx=0,y=±1) from outer regions of the (x, y) plane, such as for |x|>0.93,y=0. In this case, the unstable regions may be tongues encircling the origin. However, as the distance from the origin is decreased, the tongues can be replaced by exceedingly fine threads rapidly becoming less than say 10^–16 in thickness. Such a thread issues fromx=0.905468199,y=0 and requires of the order of 40 000 mappings to escape. It does so by spiralling about the origin and penetrating through several series of loops associated with various fixed points at successively greater (absolute) values ofx(y=0). The region between this thread and the origin is therefore highly stable. Practical stability of a region may be regarded as attained when the region is interior to a series of loops for which the trace of the variational matrix is close to 2.0. This occurs forn=53,k=4, with fixed point atx=0.819786,y=0 and Trace=2.0000 0004.If an invariant curve does in fact exist, then one must be able to show that the outward spiralling from a given series of loops is brought to a halt at some stage. This does not occur in the region where direct computation is possible, as we show in this article, and it remains to be seen under what conditions it can take place.     

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 .     

Tidal parameters from the variation of inclination of GEOS-1 and GEOS-2

Analysis of the luni-solar tidal perturbations of the inclination of GEOS-1 (1965-89A,i=59°) and GEOS-2 (1968-002A,i=106°.) has yielded the valuesk ₂=0.22 (=0.02) and 0.31 (=0.01) respectively for the apparent second degree Love number. For GEOS-1 a new purely numerical method involvingosculating elements was employed. For GEOS-2 it was necessary to analyze the variations of themean elements because of the very long period (450 days) of the dominant solar tidal perturbation. An additional analysis of the variation of the mean elements of GEOS-1 confirmed the value ofk ₂ obtained from the osculating elements. The disparate values indicate that the simple 2nd degree zonal harmonic model of the tidal potential employed by ourselves and all other investigators is accommodating other effects in addition to those caused by the solid Earth tides. A recent paper by Lambecket al. (1973) indicates that ocean tide effects have significant perturbations on satellite orbits and cannot be neglected.     

General integral transform of the exponential integrals

General integral transform of the exponential integralsE _n is considered and will be denoted asB _(k) ⁿ (). Different expressions and the equations satisfied byB _(k) ⁿ are developed. Two-term recurrence formula forB _(k) ⁿ (0) and three-term recurrence formula forB _(k) ⁿ (); 0 will be established for a givenk1 andn=2,3, ...,N. The computational algorithms based on these formulae are also constructed for the casesk=1,2,3, andn2. Finally the numerical results fork=2,3 andn=2(1)25 are presented to 15-digit accuracy     

Viscous fluid cosmological models in the presence of the zero-mass scalar field

The curvature-free (k=0) FRW expanding cosmological model is developed corresponding to interacting viscous fluids and zero-mass scalar fields. In the absence of non-static scalar fields the model exhibits the existence of the initial singularity (Q=0). However, with non-negative coefficient of shear viscosity, in the presence of non-static scalar fields we find thatQ has a minimum value (0). If this epoch is treated as the initial one, it may be said that the presence of scalar fields avoids the initial singularity. Other physical behaviour that the model exhibit has been discussed.     

Comment on the BDT-Bianchi type-VI0 stiff matter solution

We wish to point out that the Brans-Dicke-Bianchi type-VI₀ stiff matter solution recently given by Ram and Singh is not the most general solution of the corresponding field equations. Moreover, this solution has no vacuum limit and the general relativistic limit is obtained only after making an asymptotic expansion. In this paper we rediscuss the entire problem in a different way. In the limit =const., , we obtain both the stiff matter and the vacuum relativistic limit first given by Ellis and MacCallum (1969) in analytic form.