幂律电子谱轫致辐射积分的快速计算

SSW(Solar SoftWare)的能量电子产生X光子的轫致辐射积分计算发展到版本2时,其性能相比初始的版本1提高很多.在版本2的基础上,对这个积分进一步改进.通过对比几种轫致辐射积分方案,结果显示,最终的方案性能上比版本2可以快约2～5倍.在积分的精确性上比版本1及版本2均改进了很多,在缺省的积分控制精度下也不再产生光子谱的尖刺现象.而且,积分耗时不再敏感于积分上限取值.由于积分性能的提高,使得利用精确的轫致辐射截面计算轫致积分成为可能.结果显示,用精确轫致辐射截面比先前的近似截面积分的结果光子流量略小(≤4%),积分时间大约比先前使用近似截面多30%.     

A Statistical Analysis of X1- and X2- or Higher-Class Flares during Solar Cycles 21∼23

The X1- and X2- or higher class ?ares in solar cycles 21, 22, and 23 from 1986 to 2008 have been analyzed statistically in this paper. It is found in the statistical study that the number of the X1-class ?ares accounted for 52.71% of total X- and higher class ?ares, while, the number of the X2- and higher class ?ares accounted for 47.29% of total X- and higher class ?ares. No matter whether the X1- and X2- or higher class ?ares, most of them occured in the descending phases of the solar cycles. Moreover, the weaker the intensity of the solar cycle, the higher the ratio of the ?ares occurred in the descending phase of the solar cycle, and the stronger the intensity of solar ?ares, the higher the ratio of the ?ares occurred in the descending phases of the solar cycles. In addition, the phase difference between the peak of the smoothed monthly mean number of sunspots and that of the X-class ?ares has been calculated, which shows that the smoothed monthly mean number of the X1-class ?ares had a very noticeable time advance of 1 month with respect to that of sunspots in the cycles 21 and 22, but there was a time lag of 13 months in the cycle 23, while, for the X2- and higher class ?ares, there was a time lag of 9 months in the cycle 21, but a one-month time advance existed in the cycle 22, and again a time lag of 32 months appeared in the cycle 23.     

Hα Surges Initiated by Newly-emerging Satellite Magnetic Fields

On July 22, 2011 and in the active region NOAA 11259 there ap- peared the event of the ejection of solar atmospheric Hα surges. According to the full-disc Hα observations of the Big Bear Solar Observatory in United States, three consecutive surges at one and the same place in the north of the main spot of the active region were discovered. The trajectories of these three surges exhib- ited the ?gure of straight lines, and their integral con?guration is like an inverted Eiffel Tower. The ?rst two surges are quite similar, and in each of them there appeared two bright points in the northern part of the main spot. After several minutes, the surges appeared in the midst of bright points. When the bright- ness of the bright points attained the maximum value, the surges spouted out from the midst of bright points. And after reaching the maximum altitude, they quickly vanished. Before the ejection of the third surge took place, no bright points appeared. Besides, its maximal altitude is merely one half of that of the ?rst two surges. Via a comparison with the SDO/HMI (Solar Dynamics Obser- vatory/Helioseismic and Magnetic Imager) data of radial magnetic ?elds, it is found that in more than one hour before the appearance of the ?rst surge there emerged bipolar magnetic ?elds in the region of ejection. Besides, in several min- utes before the ejection of each Hα surge the magnetic ?uxes of positive polarity diminished. Via our analysis it is found that there appeared reconnections be- tween the newly emerging satellite magnetic ?elds and the preexisting magnetic ?elds in the spot, and this caused the continuous ejections of Hα surges.     

A Research on the Characteristics of Light Variation of Blazar 3C 273 in 2-10 keV

From the original observed data of RXTE (Rossi X-ray timing explorer), the data of 3C 273 in the X-ray band of 2-10 keV from February 2, 1996 (MJD 50115) to August 27, 2007 (MJD 54339) are analyzed. The photon spectral indexes and corresponding fluxes of the observational data from 1010 observations are obtained in total. The average time of each observation is 1666.76 seconds. By analyzing the spectral variation and behavior of the light variation of 3C 273 in the hard X-ray band of 2-10 keV during observational periods, it is found that there was a significant anti-correlation between the photon spectral index Γ and flux lgF_{2∼10 keV} in February 2000, March and April 2003, February, March and July 2004, as well as 2006 and 2007. During the observational periods in 1999 and from 2000 to 2007, there was also a significant anti-correlation between the monthly average behaviors of Γ and lgF_{2∼10 keV}. From the quantitative analysis of the behavior of light variation in the entire observational period, several relatively large light variations are found. The evidence of the existence of the Fe emission line with average width of 93.85±21.49 eV is also found by fitting a part of the spectra. Through analyzing the intensity of the light variation and features of Fe emission line, it is found that 3C 273 has some characteristics of blazars and Seyfert galaxies simultaneously.     

A Theoretical Probe for Excitation Mechanisms of Solar-like and Mira-like Oscillations of Stars

The linear stability analysis of the radial and non-radial oscillations for the evolutionary model of a star with the mass of 0.6∼3 M8 has been per- formed by using the nonlocal and time-dependent convection theory. The results show that the unstable low-temperature stars on the right side of the instabil- ity strip in the HR diagram can be divided into two groups. One is of the stars of solar-like oscillations, composed of the main-sequence dwarfs, subgiants, and the red giants with low- and intermediate-luminosity, which are unstable in the intermediate- and high-order (nr ≥ 12) p-modes, but stable in the low- order (nr ≤ 5) p-modes. Another is of the Mira-like stars, composed of the luminous red giants and AGB stars, which are just contrary to the solar-like stars, unstable in the low-order (nr ≤ 5) p-modes, but stable in the intermediate- and high-order (nr ≥ 12) p-modes. On the red edge of Cepheid (δ Scuti) insta- bility strip, the oscillations of solar-like and Mira-like stars can be explained uniformly by the coupling between convection and oscillations (CCO). For the low-temperature stars on the right side of the instability strip, the CCO is the dominant excitation and damping mechanism for the low- and intermediate-order p-modes, and the stochastic excitation of turbulence becomes important only for the high-order p-modes of solar-like oscillations.     

A Study on the CdZnTe Array Detector

The CdZnTe array detector is a new type of semiconductor detector being rapidly developed in recent years. It possesses a high spatial resolution and a high energy resolution, and it can work at room temperatures. This paper describes the physical properties and working principle of the CdZnTe array detector, as well as the manufacturing technology, including the chip pretreatment, passivation, ohmic electrode preparation, array template selection, and array packaging technology (micro-interconnection). For evaluating the perfor-mance of the detector, the authors have developed successfully a 4 pixel×4 pixel CdZnTe array and an 8 pixel×8 pixel CdZnTe array (with the thicknesses of 5 mm and 2 mm, the pixel size of 2 mm×2 mm, and the gaps of 0.15 mm and 0.2 mm, respectively) in cooperation with the partner. A multi-channel electronic readout system based on the ASIC (Application Specific Integrated Circuit) chip is devel-oped independently for the charge measurement of the 4 pixel×4 pixel CdZnTe array. The energy spectra and corresponding energy resolutions of the 16 pixels are obtained with the ¹³⁷Cs radiative source, among them the best resolution is 4.8%@662 kev.     

The SHASTA Code Modified by Self-adaptive Mesh and Numerical Experiment of Magnetic Reconnections

SHASTA (Sharp and Smooth Transport Algorithm) is a code with single mesh to solve the 2-dimensional magnetohydrodynamic (MHD) equations. When SHASTA is used to the numerical simulation of magnetic reconnection problem, it is modified to be the code which adopts the method of the selfadaptive mesh. The modified code can carry out refined calculations in diffusion regions. In the process of the self-adaptive calculations with SHASTA, a “plugand-play” strategy is adopted and the original algorithm to solve 2-dimensional MHD partial differential equations is treated as an independent cell. In addition, the hierarchical data structure is used in this modification and parameters in each refined level are described by a 2-dimensional variable array. The regions where the distributions of magnetic field and pressure exhibit steep variations are marked as the refined regions. Then, the distributions of physical quantities and the boundary conditions in the grid points of refined levels are deduced via interpolation method. Finally, the refined calculated results of refined regions are assigned to the previous level of mesh and the existing results are updated. The numerical experiment of magnetic reconnections which adopts refined calculations indicates that compared with the code with single mesh, the resolution of details is improved and the corresponding increment of computing time is related to the selection of parameters in the simulation. The calculation accuracy and effect on instability, which are caused by a part of the self-adaptive code, depend on the boundary settings, push strategy over each single step as well as the interpolation algorithm.     

K-shell photoionization cross-sections of N v and O vi

An 11-state R -matrix calculation is performed to obtain cross-sections for the photoionization of the 1s²2s ²S ground states of N  v and O  vi . For both ions at low photon energies, corresponding solely to the 2s photoionization cross-section, excellent agreement is found with previous calculations. For higher photon energies, permitting inner-shell photoionization of the 1s electron, excellent accord is found with the recent calculation of Nahar for O  vi , whilst for N  v resonance structure is found which previous calculations omitted.     

Radiation mechanisms of the plasma near point L1 affected by CMEs and associated microwave bursts

When coronal mass ejections (CMEs) interact with the solar corona and the interplanetary medium, emissions at different wavelengths occur. On the basis of study of the various radiation mechanisms of space plasma in the case of absence of CMEs, the radio radiation mechanisms of the plasma close to the Lagrange point L1 and affected by large CMEs from February to August 1999 are statistically analyzed. As shown by the results, the main radiation mechanisms are the Bremsstrahlung, a small amount of cyclotron radiation and a still weaker recombination radiation. Also, solar microwave bursts which are associated with CMEs in the same period are investigated. The results show that the microwave bursts are of the gradual type as well as spike bursts, and that the chief radiation mechanisms are the Bremsstrahlung, cyclotron resonance radiation, plasma radiation and electronic cyclotron maser radiation.     

The Orbit Computation of Binary Black Hole OJ287—A 3.5-th Order Post-Newtonian Approximation

The massive binary black hole OJ287 is a source of intense gravita- tional radiation. To detect the signal of its gravitational waves, a knowledge of the signal waveform will be of great help, and this is mainly determined by the orbital motion of the binary. For this, we carry out a detailed calculation on the orbital motion of OJ287, using the post-Newtonian (PN) approximation up to the 3.5th order within the framework of general relativity. Our result is one order higher than the previous work made by others. As in the process of radiation, there is a time delay from the instance when the secondary black hole impacts on the accretion disk of the primary to the moment of the optical outburst. This time delay has to be taken into consideration when we try to fit the calculating orbit with the observed times of outbursts. Adopting a linear relation between the time delay and the impact distance as an empirical model, we fit the cal- culating orbit with the recent 7 outbursts of OJ287, and obtain the solution of its orbital motion, as well as its averaged orbital parameters. By analyzing the result of 3.5 PN order calculation of the binary system, we find some interesting features. In the early period, the rate of precession of the secondary black hole increases, while in the late period approximate to merging, the rate of precession attains its maximum. Afterwards it diminishes, and finally becomes negative. At present we cannot determine whether this behavior is due to the insuffcient accuracy of the 3.5-th order approximation. For the term of dissipative radiation in the equation of motion, the coeffcients of the 2.5 and 3.5 PN orders possess opposite signs. This implies that the 3.5-th order term represents the absorption of energy from outside. However, the sum of the 2.5-th order and 3.5-th order terms still behaves as radiating gravitational waves outward, the rate of energy variation of the system is negative. The calculated result of this work may be useful for more accurate calculations of the gravitational radiation of OJ287.     

Maximum-likelihood TOA Estimation of X-ray Pulsar Signals on the Basis of Poison Model

The cycle-stationary Poison model for the photon arrival of X-ray pulsars and the estimation of the pulse's time of arrival (TOA) are discussed. Based on this model, the maximum likelihood estimation (MLE) of the TOA, as well as the Cramer-Rao boundary (CRB), are presented, especially the approximate formulation of the likelihood function and the Cramer-Rao boundary are derived for the low-R_SN (signal-to-noise ratio) cases. By using the analytical pulse profile, we have made the Monte-Carlo simulation on the TOA estimation of PSR B1821-24, discussed the estimate error for different observation times and signal-to-noise ratios, and presented the corresponding R_SN-thresholds. The results show that this method of data analysis can estimate effectively the timing accuracy of X-ray pulsar pulses and help to evaluate the corresponding performances in other applications.     

A Modified Gooding Method to Calculate Inclination Functions

In this paper, the expression of V k_lm(I) in the Gooding method is rewritten to be the form convenient for calculation, and a standard recursive lm procedure is used to calculate Ak_lm(I). We have rewritten the Gooding's program under the assumption that l and k have the same odd-even parity, this makes the program be shorten for one half, the computational effciency and readability of the program be raised, the computing time be shortened for 41%, and the computational accuracy and stability are also slightly improved.     

Analyses of the Light Variations of OJ 287, 0716+714 and BL Lacertae

OJ 287 is a BL Lac object which exhibits intense activities of low peak-frequencies. Its energy spectrum in low frequency band is quite similar with those of two other TeV BL Lac objects (i.e., 0716+714 and BL Lacertae). However, the Cerenkov telescope did not detect its TeV rays. By using the observational data of these three heavenly bodies and comparing the discrepan- cies of their minimal periods of light variations and delays at 22 GHz, 37 GHz and B-waveband, we have further investigated the possible reason why the TeV gamma-rays of OJ 287 have not been observed. The results of analyses are as fol- lows. (1) For the minimal periods of light variations, the periods of OJ 287 at 37 GHZ and B-waveband are short. At 22 GHz the results of OJ 287 and 0716+714 are comparable, but the period of OJ 287 is much shorter in comparison with that of BL Lacertae, and this shows that its activity is more intense. However, the TeV gamma-rays of OJ 287 have not been detected, which implies that the radiation of OJ 287 in TeV waveband may have no connection with the minimal periods of light variations in these three low-energy wavebands. (2) In respect of delays, the delay of OJ 287 in the B waveband with respect to 37 GHz is longer than that of 0716+714, but shorter than that of BL Lacertae. Its delay at 37 GHz in respect to 22 GHz is shorter than that of 0716+714. Meanwhile, the delay of BL Lacertae at 37 GHz in respect to 22 GHz is negative, which implies that 22 GHz precedes 37 GHz. Via the comparison and analysis of delays, no obvious differences between OJ 287 and 0716+714 as well as BL Lacertae have been found. On the side of energy spectra, it is quite possible that due to the steep energy spectrum of OJ 287 in TeV waveband, the Cerenkov telescope has not detected the gamma radiation of OJ 287. However, nowadays it is still not clear whether the steep energy spectrum in TeV energy range has some influence on the light variations in low energy realm.     

Improvement on Atmospheric Hydrostatic Delay Correction at Low Elevation Angles

Observation made at low evelation angles is the trend of development of GPS (Global Positioning System) meteorology, wherein the development of highly accurate atmospheric hydrostatic delay corrections at low elevating angles is the main key technique. The comparison among three methods for calculating the atmospheric hydrostatic delay correction of the radio waves from space to the ground-based receiver is made: (1) the atmospheric hydrostatic delay obtained from the path integration of the sounding balloon data under the assumption of atmospheric spherical symmetry, (2) the atmospheric hydrostatic delay acquired from the reanalyzed data of the NCEP (National Centers for Environmental Prediction) under the assumption of atmospheric spherical symmetry and (3) the atmospheric hydrostatic delay got from Niell's atmospheric hydrostatic mapping function. The results of the comparison of them with the atmospheric hydrostatic mapping function obtained from the calculation carried out by taking advantage of the data acquired at 89 sounding balloon stations in China in 2001 show that the accuracy of the method of the path integration of the reanalyzed data of NCEP at low elevating angles (bellow 5°) is about 5 times better than that of the Niell mapping function model.     

The use of integrals in numerical integrations of theN-body problem

The numerical integration of systems of differential equations that possess integrals is often approached by using the integrals to reduce the number of degrees of freedom or by using the integrals as a partial check on the resulting solution, retaining the original number of degrees of freedom.Another use of the integrals is presented here. If the integrals have not been used to reduce the system, the solution of a numerical integration may be constrained to remain on the integral surfaces by a method that applies corrections to the solution at each integration step. The corrections are determined by using linearized forms of the integrals in a least-squares procedure.The results of an application of the method to numerical integrations of a gravitational system of 25-bodies are given. It is shown that by using the method to satisfy exactly the integrals of energy, angular momentum, and center of mass, a solution is obtained that is more accurate while using less time of calculation than if the integrals are not satisfied exactly. The relative accuracy is ascertained by forward and backward integrations of both the corrected and uncorrected solutions and by comparison with more accurate integrations using reduced step-sizes.     

A Definite Integration Method for Calculating Inclination Function and Its Derivative

This paper gives a de?nite integration method for calculating the inclination function and its derivative, which has a very simple expression, and the accuracies as high as 10^-15 for the inclination function, and 10^-13 for its derivative, comparable with the accuracy of Gooding's method. By through a lot of numerical simulations, it is proved that this method has a good stability and an wide applicable range of inclinations, hence it can be used to calculate the inclination function to the maximum order of L_max ≤ 50.     

Studies in the application of recurrence relations to special perturbation methods

Using the rectangular equations of motion for the restricted three-body problem a comparison is made of the Encke and Cowell methods of integration. Each set of differential equations is integrated using Taylor series expansions where the coefficients of the powers of time are determined by recurrence relationships. It is shown that for fairly highly eccentric orbits in which the perturbing force is less than one thousandth of the two-body force the Encke method achieves a considerable saving in machine time. This is also true for almost circular orbits when low or moderate accuracy is required. When very high accuracy is required, however, the Cowell method is faster unless the perturbing force is less than 10^–6 of the two-body force. There is little difference in the accuracy of the two methods, the Cowell method being slightly more accurate when a low or moderate accuracy criterion is imposed.     

The OSCROX stellar oscillaton code

This paper describes the OSCROX stellar oscillation code for the calculation of the adiabatic oscillations of low degree ? of a spherical star. There are two principal versions: one in Lagrangian variables (oscroxL), the second in Eulerian variables (oscroxE). The Lagrangian code does not require values of the Brunt Väisälä frequency or equivalently the density gradient. For ?=1 the oscillation equations have both an exact integral and an exact partial wave solution, and codes oscroxL1 and oscroxE1 incorporate these exact solutions. The difference in the frequencies obtained with the various codes gives some estimate of the uncertainty in the results due both to limited accuracy of hydrostatic support of the stellar model, and the limited accuracy of the integration of the oscillation equations. We compare the results of the different methods by calculating the frequencies in the range 20–2500 μHz of a model of a 1.5 M_⊙ main-sequence star (ModelJC) kindly provided by J. Christensen-Dalsgaard for the purposes of cross comparison of codes, a modified version of this model (ModelJCA) with improved hydrostatic support, and of a highly accurate n=3 polytropic model of a star with the same mass and radius. For the polytropic model the frequencies as calculated by all codes agree to within 0.001 μHz, whereas for the 1.5 M_⊙ main sequence model the frequency differences reach a maximum of 0.04 μHz, due primarily to the limited accuracy of hydrostatic support in the model; this is reduced to 0.01 μHz for ModelJCA.     

Evolutionary Characteristics of Total Solar Irradiance

The significant periods of total solar irradiance are 35 d and 26 d in the 23rd and 24th solar activity cycles, respectively. It is inferred that the solar quasi-rotation periods are also 35 d and 26 d in the 23rd and 24th solar activity cycles, respectively. The value of total solar irradiance around the 24th solar activity minimum may be close to the value of Maunder minimum. On the timescales from one solar rotation period to several months, sunspots are the main reason to cause the variation of total solar irradiance, but not the unique one, and the variation of total solar irradiance are not correlated with the Mg II index on the timescales from a few days to one solar rotation period.