On equivalent widths in a penumbral spectrum (a test of the Moe-Maltby model)

From new observational material we made a curve of growth analysis of the penumbra of a large, stable sunspot. The analysis was done relative to the undisturbed photosphere and gave the following results (⊙ denotes photosphere, * denotes penumbra): $$\begin{gathered} (\theta ^ * - \theta ^ \odot )_{exe} = 0.051 \pm 0.007 \hfill \\ {{\xi _t ^ * } \mathord{\left/ {\vphantom {{\xi _t ^ * } {\xi _t }}} \right. \kern-\nulldelimiterspace} {\xi _t }}^ \odot = 1.3 \pm 0.1 \hfill \\ {{P_e ^ * } \mathord{\left/ {\vphantom {{P_e ^ * } {P_e ^ \odot = 0.6 \pm 0.1}}} \right. \kern-\nulldelimiterspace} {P_e ^ \odot = 0.6 \pm 0.1}} \hfill \\ {{P_g ^ * } \mathord{\left/ {\vphantom {{P_g ^ * } {P_g }}} \right. \kern-\nulldelimiterspace} {P_g }}^ \odot = 1.0 \pm 0.2 \hfill \\ \end{gathered} $$ The results of the analysis are in satisfactory agreement with the penumbral model as published by Kjeldseth Moe and Maltby (1969). Additionally we tested this model by computing the equivalent widths of 28 well selected lines and comparing them with our observations.     

Changes in binding energy of interacting galaxies

It is shown that the fractional increase in binding energy of a galaxy in a fast collision with another galaxy of the same size can be well represented by the formula $$\xi _2 = 3({G \mathord{\left/ {\vphantom {G {M_2 \bar R}}} \right. \kern-\nulldelimiterspace} {M_2 \bar R}}) ({{M_1 } \mathord{\left/ {\vphantom {{M_1 } {V_p }}} \right. \kern-\nulldelimiterspace} {V_p }})^2 e^{ - p/\bar R} = \xi _1 ({{M_1 } \mathord{\left/ {\vphantom {{M_1 } {M_2 }}} \right. \kern-\nulldelimiterspace} {M_2 }})^3 ,$$ whereM ₁,M ₂ are the masses of the perturber and the perturbed galaxy, respectively,V _p is the relative velocity of the perturber at minimum separationp, and \(\bar R\) is the dynamical radius of either galaxy.     

Astrophysical thermonuclear functions

Stars are gravitationally stabilized fusion reactors changing their chemical composition while transforming light atomic nuclei into heavy ones. The atomic nuclei are supposed to be in thermal equilibrium with the ambient plasma. The majority of reactions among nuclei leading to a nuclear transformation are inhibited by the necessity for the charged participants to tunnel through their mutual Coulomb barrier. As theoretical knowledge and experimental verification of nuclear cross sections increases it becomes possible to refine analytic representations for nuclear reaction rates. Over the years various approaches have been made to derive closed-form representations of thermonuclear reaction rates (Critchfield, 1972; Haubold and John, 1978; Haubold, Mathai and Anderson, 1987). They show that the reaction rate contains the astrophysical cross section factor and its derivatives which has to be determined experimentally, and an integral part of the thermonuclear reaction rate independent from experimental results which can be treated by closed-form representation techniques in terms of generalized hypergeometric functions. In this paper mathematical/statistical techniques for deriving closed-form representations of thermonuclear functions, particularly the four integrals $$\begin{gathered} I_1 (z,v)\mathop = \limits^{def} \int\limits_0^\infty {y^v e^{ - y} e^{ - zy^{ - {1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-\nulldelimiterspace} 2}} } dy,} \hfill \\ I_2 (z,d,v)\mathop = \limits^{def} \int\limits_0^\infty {y^v e^{ - y} e^{ - zy^{ - {1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-\nulldelimiterspace} 2}} } dy,} \hfill \\ I_3 (z,t,v)\mathop = \limits^{def} \int\limits_0^\infty {y^v e^{ - y} e^{ - z(y + 1)^{ - {1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-\nulldelimiterspace} 2}} } dy,} \hfill \\ I_4 (z,\delta ,b,v)\mathop = \limits^{def} \int\limits_0^\infty {y^v e^{ - y} e^{ - by^\delta } e^{ - zy^{ - {1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-\nulldelimiterspace} 2}} } dy,} \hfill \\ \end{gathered} $$ will be summarized and numerical results for them will be given. The separation of thermonuclear functions from thermonuclear reaction rates is our preferred result. The purpose of the paper is also to compare numerical results for approximate and closed-form representations of thermonuclear functions. This paper completes the work of Haubold, Mathai, and Anderson (1987).     

A comparison of theoretical Fe xii emission line strengths with EUV observations of a solar active region

New theoretical electron-density-sensitive Fe xii emission line ratios $$R_1 = I(3s^2 3p^3 {}^4S_{3/2} - 3s3p^4 {}^4P_{5/2} )/I(3s^2 3p^3 {}^2P_{3/2} - 3s3p^4 D_{5/2} )$$ and $$R_2 = I(3s^2 3p^3 {}^2P_{3/2} - 3s3p^4 {}^2D_{5/2} )/I(3s^2 3p^3 {}^4S_{3/2} - 3s3p^2 P_{3/2} )$$ are derived using R-matrix electron impact excitation rate calculations. We have identified the Fexii \(3s^2 3p^3 {}^4S_{3/2} - 3s3p^4 {}^4P_{5/2} ,{\text{ }}3s^2 3p^3 {}^2P_{3/2} - 3s^3 3p^4 {}^2D_{5/2} ,{\text{ }}3s^2 3p^3 S_{3/2} - 3s^2 3p^3 P_{3/2} \) and \(3s^2 3p^3 {}^4S_{3/2} - 3s^2 3p^3 {}^2P_{1/2}\) transitions in an active region spectrum obtained with the Harvard S-055 spectrometer on board Skylab at wavelengths of 364.0, 382.8, 1241.7, and 1349.4 Å, respectively. Electron densities determined from the observed values of R ₁ (log N _e ? 11.0) and R ₂(log N _e ? 11.4) are significantly larger than the typical active region measurements, but are similar to those derived from some active region spectra observed with the Skylab 2082A instrument, which provides observational support for the atomic data adopted in the line ratio calculations, and also for the identification of the Fe xii transitions in the S-055 spectrum. However the observed value of R ₃ = I(1349.4 Å)/I(1241.7 Å) is approximately a factor of two larger than one would expect from theory which, considering that the 1349.4 Å line lies at the edge of the S-055 wavelength coverage, may reflect errors in the instrument efficiency curve. Another possibility is that the 1349.4 Å transition is blended, probably with Si ii 1350.1 Å.     

Apsidal motion in the close binary IT cassiopeiae

In 1982 and 1993, we carried out highly accurate photoelectric WBVR measurements for the close binary IT Cas. Based on these measurements and on the observations of other authors, we determined the apsidal motion $\left[ {\dot \omega _{obs} = {{(11\mathop .\limits^ \circ 0 \pm 2\mathop .\limits^ \circ 5)} \mathord{\left/ {\vphantom {{(11\mathop .\limits^ \circ 0 \pm 2\mathop .\limits^ \circ 5)} {100 years}}} \right. \kern-0em} {100 years}}} \right]$ . This value is in agreement with the theoretically calculated apsidal motion for these stars $\left[ {\dot \omega _{th} = {{(14^\circ \pm 3^\circ )} \mathord{\left/ {\vphantom {{(14^\circ \pm 3^\circ )} {100 years}}} \right. \kern-0em} {100 years}}} \right]$ .     

Ion cyclotron instability in the solar wind

An ion cyclotron instability, arising because of the relative drift between the beam and the main components of the proton distribution function in the solar wind at 1 AU, is studied. The instability is excited in a bounded range of wave numbers provided the relative drift exceeds a certain minimum value called instability threshold. For 1, the instability threshold is smaller than or equal to the threshold of magnetosonic and Alfvén instabilities. The growth rates are enhanced by increasing relative drift and ratio of beam to main proton number density and by decreasing the wave numbers.     

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:

Deceleration of the Earth's rotation from old solar tables

The solar tables of ibn Yunis and of King Alfonso, those of Kepler, of G. D. and J.-J. Cassini, and of Lalande are compared with Newcomb's theory of the Sun to determine the deceleration of the Earth's rotation. Comparisons of mean motion and of longitude lead to separate determinations. A value for the deceleration, , assumed to be independent of time is obtained for the period — 146 through 1892. This result is based on the assumption that Newcomb's theory converts the Earth's orbital motion into a perfect clock, and the result is independent of lunar and tidal theory. The deceleration seems greater than that obtained from the literature for the tidal deceleration caused by the Moon.     

Characteristic actions ħ(s) in the structure of the universe

It is suggested that gravitationally bound systems in the Universe can be characterized by a set of actions ?^(s). The actions $$\hbar ^{\left( s \right)} = \left( {{\hbar \mathord{\left/ {\vphantom {\hbar {\frac{1}{{2\pi }}\frac{{C^5 }}{{GH_0^2 }}}}} \right. \kern-\nulldelimiterspace} {\frac{1}{{2\pi }}\frac{{C^5 }}{{GH_0^2 }}}}} \right)^{s/6} \left( {\frac{1}{{2\pi }}\frac{{C^5 }}{{GH_0^2 }}} \right)$$ ,derived from general theoretical consideration, are only determined by the fundamental physical constants (Planck's action ?, the velocity of lightC, gravitational constantG, and Hubble's constantH ₀) and a scale parameters. It is shown thats=1, 2, and 3 correspond, respectively, to the scales of galaxies, stars, and larger asteroids. The spectra of the characteristic angular momenta and masses for gravitationally bound systems in the Universe are estimated byJ ^(s) andM ^(s) =(? ^(s) Cα/G)^1/2. Taken together, an angular momentum-mass relation is obtained,J ^(s)=A(M^(s))², where \(A = G/C\alpha ,{\text{ }}\alpha \simeq \tfrac{{\text{1}}}{{{\text{137}}}}\) , for the astronomical systems observed on every scale. ThisJ-M relation is consistent with Brosche's empirical relation (Brosche, 1974).     

Solar-wind properties at the Earth as predicted by the two-fluid model

The two-fluid equations for the solar wind are written down in a simplified form, similar to that suggested by Roberts (1971) for the one-fluid model. The equations are shown to depend only on one parameter, $$K = GM\kappa _e m_p (\varepsilon _\infty T_0 )^{{3 \mathord{\left/ {\vphantom {3 2}} \right. \kern-\nulldelimiterspace} 2}} /4k^2 Fe,$$ , where G is the gravitational constant, M the mass of the star, κ _e the thermal electron conductivity, m _p the proton mass, k the Boltzman constant, k? _∞T₀ the residual energy per particle at infinity and F _e the electron-particle flux. For a variety of values of the density and temperature at the base of the corona we compute the solutions of the two-fluid solar wind model and compare the predicted and observed solar wind parameters at the Earth.     

The comparison of the magnetographic magnetic field measured in different spectral lines

The discrepancies in the values of longitudinal magnetic field obtained from magnetographic records in different spectral lines are considered. On the basis of extensive data including 60 pairs of magnetographic maps for 11 spectral lines, obtained simultaneously for one of these lines and 6103 with the aid of the Crimean double channel magnetograph, the following conclusions have been reached. The relative field strength (6103) depends partly on the distance from the center of the disk (Figure 4) and mainly on the magnetic sensitivity of the line g² (Figure 3), pointing to the primary role of saturation effect. The possible influence of line asymmetry on these discrepancies is also suggested.     

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,

Persistence of solar activity on small scales: Hurst analysis of time series coming from Hα flares

We calculated the Hurst exponent H for the daily averaged intensity Q of optical flares, an index which describes the solar activity. We found that H0.74±0.02 in the range of scales from about 20 days up to 450 days. This value is well beyond H= , expected for a stochastic Brownian process, thus indicating that the solar cycle could show persistence on small scales, in agreement with what has been found using other indices of the solar cycle.     

Simulations of cosmic-ray particle diffusion

The diffusion of charged particles in a stochastic magnetic field (strengthB) which is superimposed on a uniform magnetic fieldB ₀ k is studied. A slab model of the stochastic magnetic field is used. Many particles were released into different realizations of the magnetic field and their subsequent displacements z in the direction of the uniform magnetic field numerically computed. The particle trajectories were calculated over periods of many particle scattering times. The ensemble average was then used to find the parallel diffusion coefficient . The simulations were performed for several types of stochastic magnetic fields and for a wide range of particle gyro-radius and the parameterB/B ₀. The calculations have shown that the theory of charged particle diffusion is a good approximation even when the stochastic magnetic field is of the same strength as the uniform magnetic field.     

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

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

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} $$      

A photometric and polarimetric study of the moon's surface

The results of the photometric and polarimetric observation of the Moon's surface with high resolution from October 1969 to March 1971 are discussed. It is found that, for a wide part of the Moon's illuminated surface, the brightnessm, expressed by visual magnitude per square second of arc, can be expressed as $$m = 5.0 - 2.5\log H + 2.5\log \{ {{P_v } \mathord{\left/ {\vphantom {{P_v } P}} \right. \kern-\nulldelimiterspace} P}(\alpha )\} ,$$ whereH is the variable part of Hapke's formula,P _v is the observed polarization degree, andP(α) is a function of the phase angleα. The color shows a tendency to reddening at the enhancement of the lunar brightness, and this is expressed by the relation $$B - V = - 0.2\Delta m + 1.7,$$ where Δm =m + 2.5 logH which is a sort of residual from the brightness calculated by Hapke's formula. Remarkable enhancements of the lunar brightness correspond to solar flares which appeared just before the observations. Because of this the above formula for the observed brightness can be interpreted by assuming the luminescence on the lunar surface stimulated by the solar activity.     

The cosmogonic shadow effect

We consider the Alfvén-Arrhenius fall-down mechanism and describe an approximate model for the infall, capture and distribution of dust particles on a given magnetic field line and their possible neutralization at the ‘2’/3 points, the points at which the field aligned compnents of the gravitational and centrifugal forces are equal and opposite. We find that a small fraction (<10%) of an incoming particle distribution will actually contribute to the above ‘2’/3 fall-down process. We also show that if at the 2/3 points, the ratio of dust to plasma density is $$\frac{{n_D \left( {\tfrac{2}{3}} \right)}}{{n_p \left( {\tfrac{2}{3}} \right)}} > \frac{{10^{ - 3} }}{{r_{g_\mu } T_{eV} }}$$ . (r _gμ=radius of a grain in microns,T=plasma temperature in eV), then the dust particles will lose their charge, decouple from the field line and follow Keplerian orbits in accordance with the Alfvén-Arrhenius mechanism. We then determine the limits on the plasma parameters in order that rotation of a quasi-neutral plasma in thermal equilibrium be possible in the gravitational and dipole field of a rotating central body. The constraints imposed by the above conditions are rather weak, and the plasma parameters can have a wide range of values. For a plasma corotating with an angular velocity Ω～10^?4s^?1, we show that the plasma temperature and density must satisfy $$10^{ - 1}<< T_{(eV)}<< 10^2 ,10T_{eV}^2<< n^p \left( {cm^3 } \right)<< 10^6 $$ .