A closed-form solution to the minimum $${\Delta V_{\rm tot}^2}$$ Lambert’s problem

A closed form solution to the minimum DV_tot²{\Delta V_{\rm tot}^2} Lambert problem between two assigned positions in two distinct orbits is presented. Motivation comes from the need of computing optimal orbit transfer matrices to solve re-configuration problems of satellite constellations and the complexity associated in facing these problems with the minimization of DV_tot{\Delta V_{\rm tot}}. Extensive numerical tests show that the difference in fuel consumption between the solutions obtained by minimizing DV_tot²{\Delta V_{\rm tot}^2} and DV_tot{\Delta V_{\rm tot}} does not exceed 17%. The DV_tot²{\Delta V_{\rm tot}^2} solution can be adopted as starting point to find the minimum DV_tot{\Delta V_{\rm tot}}. The solving equation for minimum DV_tot²{\Delta V_{\rm tot}^2} Lambert problem is a quartic polynomial in term of the angular momentum modulus of the optimal transfer orbit. The root selection is discussed and the singular case, occurring when the initial and final radii are parallel, is analytically solved. A numerical example for the general case (orbit transfer “pork-chop” between two non-coplanar elliptical orbits) and two examples for the singular case (Hohmann and GTO transfers) are provided.     

Investigation of the Neupert Effect in the Various Intervals of Solar Flares

The Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) gives us a chance to investigate the theoretical Neupert effect using the correlation between the thermal-energy derivative and the nonthermal energy, or the thermal energy and the integral nonthermal energy. Based on this concept, we analyze four M-class RHESSI flares on 13 November 2003, 4 November 2004, 3 and 25 August 2005. According to the evolution of the temperature [T], emission measure [EM], and thermal energy [E _th], each event is divided into three phases during the nonthermal-energy input [ \frac dE_nthdt\frac {\mathrm{d}E_{\mathrm{nth}}}{\mathrm{d}t} in the units of erg s⁻¹]. Phase 1 is identified as the interval before the temperature maximum, while after the thermal-energy maximum is phase 3, between them is phase 2. We find that these four flares show the Neupert effect in phase 1, but not in phase 3. The Neupert effect still works well in the second phase, although the cooling becomes slightly important. We define the parameter μ in the relation of \fracdE_thdt=m\fracdE_nth(t)dt\frac{\mathrm {d}E_{\mathrm{th}}}{\mathrm{d}t}=\mu\frac{\mathrm{d}E_{\mathrm {nth}}(t)}{\mathrm{d}t} or E_th(t₀)=mò₀^t₀\fracdE_nth(t)dt dtE_{\mathrm{th}}(t_{0})=\mu\int_{0}^{t_{0}}\frac{\mathrm{d}E_{\mathrm{nth}}(t)}{\mathrm{d}t}\,\mathrm{d}t when the cooling is ignored in phase 1. Considering the uncertainties in estimating the energy from the observations, it is not possible to precisely determine the fraction of the known energy in the nonthermal electrons transformed into the thermal energy of the hottest plasma observed by RHESSI. After a rough estimate of the flare volume and the assumption of the filling factor, we investigate the parameter μ in these four events. Its value ranges from 0.02 to 0.20, indicating that a small fraction (2% – 20%) of the nonthermal energy can be efficiently transformed into thermal energy, which is traced by the soft X-ray emission, and the bulk of the energy is lost possibly due to cooling.     

Small-Scale Flux Emergence Observed Using <Emphasis Type="Italic">Hinode</Emphasis>/SOT

The aim of this paper is to determine the flux emergence rate due to small-scale magnetic features in the quiet Sun using high-resolution Hinode SOT NFI data. Small-scale magnetic features are identified in the data using two different feature identification methods (clumping and downhill); then three methods are applied to detect flux emergence events. The distribution of the intranetwork peak emerged fluxes is determined. When combined with previous emergence results, from ephemeral regions to sunspots, the distribution of all fluxes are found to follow a power-law distribution which spans nearly seven orders of magnitude in flux (10¹⁶ – 10²³ Mx) and 18 orders of magnitude in frequency. The power-law fit to all these data is of the form

The Variation of Solar Wind Correlation Lengths Over Three Solar Cycles

We present the results of a study of solar wind velocity and magnetic field correlation lengths over the last 35 years. The correlation length of the magnetic field magnitude λ _|B| increases on average by a factor of two at solar maxima compared to solar minima. The correlation lengths of the components of the magnetic field l_BXYZ\lambda_{B_{XYZ}} and of the velocity l_VYZ\lambda_{V_{YZ}} do not show this change and have similar values, indicating a continual turbulent correlation length of around 1.4×10⁶ km. We conclude that a linear relation between λ _|B|, VB ², and Kp suggests that the former is related to the total magnetic energy in the solar wind and an estimate of the average size of geoeffective structures, which is, in turn, proportional to VB ². By looking at the distribution of daily correlation lengths we show that the solar minimum values of λ _|B| correspond to the turbulent outer scale. A tail of larger λ _|B| values is present at solar maximum causing the increase in mean value.     

The Correlation between the Magnetic and Velocity Fields on the Full Solar Disk

A possible correlation between the magnetic and velocity fields has been analyzed based on the SOHO/MDI magnetograms and Dopplergrams. It is found that the observed large-scale weak magnetic field (weaker than 50 G (gauss)) is correlated with the velocity statistically. The curves of u⋅b with latitude, where u and b are the velocity and magnetic fields in a rectangular region (±15^○ in longitude, ±45^○ in latitude) on the Sun, show the same patterns in the years 2000, 2004, and 2007. The patterns indicate that u and b are positively correlated near the equator but are anti-correlated at the middle latitudes. For a strong magnetic field between 50 G and 3000 G, the curves of u⋅b with latitude show the same tendencies at the middle latitudes. Near the equator, however, the slope of the curve is positive in 2000 and is negative in 2004 and 2007. In addition, we give an estimation for the amplitude of the cross helicity h _χ (h_c=[`(u·b)]h_{\chi}=\overline{\mathbf{u}\cdot\mathbf{b}}) inferred from the MDI data, which is of the order of 10³ G m s⁻¹ near the center of the solar disk.     

Periodic solutions for resonance 4/3: Application to the construction of an intermediary orbit for Hyperion’s motion

The motion of Hyperion is an almost perfect application of second kind and second genius orbit, according to Poincaré’s classification. In order to construct such an orbit, we suppose that Titan’s motion is an elliptical one and that the observed frequencies are such that 4n _H−3n _T+3n _ω=0, where n _H, n _T are the mean motions of Hyperion and Titan, n _ω is the rate of rotation of Hyperion’s pericenter. We admit that the observed motion of Hyperion is a periodic motion such as . Then, .N _H, N _T, k∈ N ⁺. With that hypothesis we show that Hyperion’s orbit tends to a particular periodic solution among the periodic solutions of the Keplerian problem, when Titan’s mass tends to zero. The condition of periodicity allows us to construct this orbit which represents the real motion with a very good approximation. This revised version was published online in July 2006 with corrections to the Cover Date.     

Charged analogues of Schwarzschild interior solution in terms of pressure

Recently, Bijalwan (Astrophys. Space Sci., doi:, 2011a) discussed charged fluid spheres with pressure while Bijalwan and Gupta (Astrophys. Space Sci. 317, 251–260, 2008) suggested using a monotonically decreasing function f to generate all possible physically viable charged analogues of Schwarzschild interior solutions analytically. They discussed some previously known and new solutions for Schwarzschild parameter u( = \fracGMc²a ) ￡ 0.142u( =\frac{GM}{c^{2}a} ) \le 0.142, a being radius of star. In this paper we investigate wide range of u by generating a class of solutions that are well behaved and suitable for modeling Neutron star charge matter. We have exploited the range u≤0.142 by considering pressure p=p(ω) and f = ( f₀(1 - \fracR²(1 - w)a²) +f_a\fracR²(1 - w)a² )f = ( f_{0}(1 - \frac{R^{2}(1 - \omega )}{a^{2}}) +f_{a}\frac{R^{2}(1 - \omega )}{a^{2}} ), where w = 1 -\fracr²R²\omega = 1 -\frac{r^{2}}{R^{2}} to explore new class of solutions. Hence, class of charged analogues of Schwarzschild interior is found for barotropic equation of state relating the radial pressure to the energy density. The analytical models thus found are well behaved with surface red shift z _s ≤0.181, central red shift z _c ≤0.282, mass to radius ratio M/a≤0.149, total charge to total mass ratio e/M≤0.807 and satisfy Andreasson’s (Commun. Math. Phys. 288, 715–730, 2009) stability condition. Red-shift, velocity of sound and p/c ² ρ are monotonically decreasing towards the surface while adiabatic index is monotonically increasing. The maximum mass found to be 1.512 M _Θ with linear dimension 14.964 km. Class of charged analogues of Schwarzschild interior discussed in this paper doesn’t have neutral counter part. These solutions completely describe interior of a stable Neutron star charge matter since at centre the charge distribution is zero, e/M≤0.807 and a typical neutral Neutron star has mass between 1.35 and about 2.1 solar mass, with a corresponding radius of about 12 km (Kiziltan et al., [astro-ph.GA], 2010).     

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.     

Interacting Ricci dark energy with logarithmic correction

Motivated by the holographic principle, it has been suggested that the dark energy density may be inversely proportional to the area A of the event horizon of the universe. However, such a model would have a causality problem. In this work, we consider the entropy-corrected version of the holographic dark energy model in the non-flat FRW universe and we propose to replace the future event horizon area with the inverse of the Ricci scalar curvature. We obtain the equation of state (EoS) parameter ω _Λ, the deceleration parameter q and W_D￠\Omega_{D}' in the presence of interaction between Dark Energy (DE) and Dark Matter (DM). Moreover, we reconstruct the potential and the dynamics of the tachyon, K-essence, dilaton and quintessence scalar field models according to the evolutionary behavior of the interacting entropy-corrected holographic dark energy model.     

Recent glitches detected in the Crab pulsar

From 2000 to 2010, monitoring of radio emission from the Crab pulsar at Xinjiang Observatory detected a total of nine glitches. The occurrence of glitches appears to be a random process as described by previous researches. A persistent change in pulse frequency and pulse frequency derivative after each glitch was found. There is no obvious correlation between glitch sizes and the time since last glitch. For these glitches Δν _p and D[(n)\dot]_p\Delta\dot{\nu}_{p} span two orders of magnitude. The pulsar suffered the largest frequency jump ever seen on MJD 53067.1. The size of the glitch is ∼6.8×10⁻⁶ Hz, ∼3.5 times that of the glitch occurred in 1989 glitch, with a very large permanent changes in frequency and pulse frequency derivative and followed by a decay with time constant ∼21 days. The braking index presents significant changes. We attribute this variation to a varying particle wind strength which may be caused by glitch activities. We discuss the properties of detected glitches in Crab pulsar and compare them with glitches in the Vela pulsar.     

Line profile analysis of a coronal formation observed near a quiescent prominence: Intensities,temperatures and velocity fields

A technique developed for analysing line profiles with both speed and high accuracy was used to study the physical conditions of a coronal formation near a quiescent prominence. Detailed analyses of five coronal lines (Fe xiv λ 5303, Fe x λ 6374, Ni xv λ 6702, Fe xv λ 7059, and Fe xi λ 7892) provided total intensities, Doppler width temperatures, ionization temperatures, and velocities. Dissimilar spatial fluctuations in intensity are obvious for ions grouped according to (low vs high) ionization potentials. The intensity of the green line shows a local minimum around the observed quiescent prominence; a corresponding but much more diffuse pattern is visible in the red line intensity. Large differences are observed in temperatures derived by different means. In particular, , while , and . The differences between and are taken as direct evidence of temperature inhomogeneity. One can thus put little significance in T _e (xi/x). T _D(λ5303) and T _e (xv/xiv) fluctuate nearly in parallel at each slit height, with a weak local minimum evident around the prominence. The discrepancy between these two can be removed if a non-thermal turbulent motion of 6–16 km s⁻¹ is assumed. Variations with height of both T _D(λ5303) and T _e (xv/xiv) suggest that the coronal temperature maximum is located no more than 15000 km above the top of spicules. A negative gradient of about 6 deg km⁻¹ is found in the height variation of T _D(λ5303). The height variation of the green line wavelength shows that the majority of coronal material in this region is flowing from west to east on the Sun, with the highest velocity of 12 km s⁻¹ found at the lowest heights. This motion is in the same sense as that of the nearby coronal rain, as determined both from the spectra and wavelength-shifted Hα filtergrams. Superposed on the above flow is a systematic velocity field of up to ±5 km s⁻¹. This field similarly reaches maximum amplitudes at lowest heights showing a local maximum around the prominence. On leave from Institute of Earth Science and Astrophysics, Shiga University, Ohtsu 520, Japan, as 1973–75 National Academy of Science/National Research Council Senior Post-Doctoral Research Associate at Sacramento Peak Observatory.     

Accelerated expansion from?a?modified-quadratic gravity

We investigate the late-time dynamics of a four-dimensional universe based on modified scalar field gravity in which the standard Einstein-Hilbert action R is replaced by f(φ)R+f(R) where f(φ)=φ ² and f(R)=AR ²+BR _μν R ^μν,(A,B)∈ℝ. We discussed two independent cases: in the first model, the scalar field potential is quartic and for this special form it was shown that the universe is dominated by dark energy with equation of state parameter w≈−0.2 and is accelerated in time with a scale factor evolving like a(t)∝t ^5/3 and B+3A≈0.036. When, B+3A→∞ which corresponds for the purely quadratic theory, the scale factor evolves like a(t)∝t ^1/2 whereas when B+3A→0 which corresponds for the purely scalar tensor theory we found when a(t)∝t ^1.98. In the second model, we choose an exponential potential and we conjecture that the scalar curvature and the Hubble parameter vary respectively like R=hH[(f)\dot]/f,h ? \mathbbRR=\eta H\dot{\phi}/\phi,\eta\in\mathbb{R} and H=g[(f)\dot]^c,(g,c) ? \mathbbRH=\gamma\dot{\phi}^{\chi},(\gamma,\chi)\in\mathbb{R}. It was shown that for some special values of  χ, the universe is free from the initial singularity, accelerated in time, dominated by dark or phantom energy whereas the model is independent of the quadratic gravity corrections. Additional consequences are discussed.     

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.     

Study of Non-spinning black holes with reference to the change in energy and entropy

In this research paper, we have derived the formula for both the changes in energy (δE) and entropy (δS) and thereafter calculated the change in entropy (δS) with corresponding change in energy (δE) taking account the first law of the black hole mechanics relating the change in mass M, angular momentum J, horizon area A and charge Q, of a stationary black hole, when it is perturbed, given by formula satisfying in the vacuum as dM = \frack8p dA + WdJ - udQ\delta M = \frac{k}{8\pi} \delta A + \Omega\delta J - \upsilon\delta Q, specially for Non-spinning black holes.     

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.     

Completely Integrable Systems Connected with Lie Algebras

In a recent paper Ballersteros and Ragnisco (1998) have proposed a new method of constructing integrable Hamiltonian systems. A new class of integrable systems may be devised using the following sequence: , where A is a Lie algebra is a Lie–Poisson structure on R ₃, C is a Casimir for is a reduced Poisson bracket and (A, ▵) is a bialgebra. We study the relation between a Lie-Poisson stucture Λ and a reduced Poisson bracket , which is a key element in using the Lie algebra A to constructing this sequence. New examples of Lie algebras and their related integrable Hamiltonian systems are given. This revised version was published online in July 2006 with corrections to the Cover Date.     

The application of the quickest descent method to the analysis of the Hα line in loops

The quickest descent method and the multi-layer model are proposed to analyse symmetrical and asymmetrical line profiles of loops in this paper. By fitting observed lines, the physical parameters such as Doppler width, Δλ _D , the optical thickness at line centre,τ _λ0, and the line displacement, Δλ ₀, can be directly obtained, Using this method, some symmetrical and asymmetrical line profiles of both the loop prominence on February 18, 1984 and the loop on August 17, 1989 are fitted. It is found that the method mentioned in this paper is characterized by rapid convergence, high precision, and less scatter.     

Constraining the gravitational redshift of a neutron star from the Σ-meson well depth

The effect of the Σ-meson well depth on the gravitational redshift is examined within the framework of relativistic mean field theory for the baryon octet system. It is found that, for a stable neutron star, the gravitational redshift increases with the central energy density increase or with the mass increase but decreases as the radius increases. Considering a change of U_S^(N)U_{\Sigma}^{(N)} from −30 MeV to 30 MeV, for a stable neutron star the gravitational redshift near to the maximum mass increases. In addition, it is also found that the growth of the U_S^(N)U_{\Sigma}^{(N)} makes the gravitational redshift as a function of M _max/R increase, the higher the U_S^(N)U_{\Sigma}^{(N)} the less the change in the gravitational redshift.     

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.