We revisit the sunspot drawings made by the Japanese astronomer Kunitomo Toubei during 1835-1836 and recount the sunspot group number for each image. There are two series of drawings, preliminary (P, containing 17 days with observations) and summary (S, covering 156 days with observations), all made using brush and ink. S is a compilation of drawings for the period from February 1835, to March 1836. Presently, the P drawings are available only for one month, September 1835; those of other periods have presumably been lost. Another drawing (I) lets us recover the raw group count (RGC) for 25 September 1836, on which the RGC has not been registered in the existing catalogs. We also revise the RGCs from P and S using the Zurich classification and determine that Kunitomo's results tend to yield smaller RGCs than those of other contemporary observers. In addition, we find that Kunitomo's RGCs and spot areas have a correlation (0.71) that is not very different from the contemporary observer Schwabe (0.82). Although Kunitomo's spot areas are much larger than those determined by Schwabe due to skill and instrument limitations, Kunitomo at least captured the growing trend of the spot activity in the early phase of the Solar Cycle 8. We also determine the solar rotation axis to estimate the accurate position (latitude and longitude) of the sunspot groups in Kunitomo's drawings.

A flare kernel associated with a C4 class flare was observed in a spectral window including the He I triplet 10830 A and Si I 10827 A with a spectropolarimeter on the Domeless Solar Telescope at Hida Observatory on August 9th, 2015. Observed Stokes profiles of the He I triplet in the flare kernel in its post-maximum phase are well reproduced through inversions considering the Zeeman and the Paschen-Back effects with a three slab model of the flare kernel, in which two slabs having up and downward velocities produce emissions and one slab produces an absorption. The magnetic field strength inferred from the emission components of the He I line is 1400 G, which is significantly stronger than 690 G that is observed at the same location in the same line 6.5 hours before the flare. In addition, photospheric magnetic field vector derived from the Si I10827 A is similar to that of the flare kernel. To explain this results, we suggest that the emission in the He I triplet during the flare is produced in the deep layer, around which bombardment of non-thermal electrons leads to the formation of a coronal temperature plasma. Assuming a hydrogen column density at the location where the He I emissions are formed, and a power-law index of non-thermal electron energy distribution, we derived the low-energy cutoff of the non-thermal electron as 20 - 30 keV, which is consistent with that inferred from hard X-ray data obtained by RHESSI.

The total solar fluxes at 1, 2, 3.75, and 9.4 GHz were observed continuously from 1957 to 1994 at Toyokawa, and from 1994 until now at Nobeyama, Japan with the current Nobeyama Radio Polarimeters. We examined the multi-frequency and long-term datasets, and found that not only the microwave solar flux but also its monthly standard deviation well indicates the long-term variation of solar activity. Furthermore, we found that the microwave spectra at the solar minima of Cycle 20~24 agree with each other. These results show that the average atmospheric structure above the upper chromosphere in the quiet Sun has not varied for half a century, and suggest that the energy input for atmospheric heating from the sub-photosphere to the corona has not changed in the quiet Sun despite significantly differing strengths of magnetic activity in the last five solar cycles.

We investigate the upper chromosphere and the transition region of the sunspot umbra using the radio brightness temperature at 34 GHz (corresponding to 8.8-mm observations) as observed by the Nobeyama Radioheliograph (NoRH). Radio free-free emission in the longer millimeter range is generated around the transition region, and its brightness temperature yields the region's temperature and density distribution. We use the NoRH data at 34 GHz by applying the Steer-CLEAN image synthesis. These data and the analysis method enable us to investigate the chromospheric structures in the longer millimeter range with high spatial resolution and sufficient visibilities. We also perform simultaneous observations of one sunspot using the NoRH and the Nobeyama 45-m telescope operating at 115 GHz. We determine that 115-GHz emission mainly originates from the lower chromosphere while 34-GHz emission mainly originates from the upper chromosphere and transition region. These observational results are consistent with the radio emission characteristics estimated from the current atmospheric models of the chromosphere. On the other hand, the observed brightness temperature of the umbral region is almost the same as that of the quiet region. This result is inconsistent with the current sunspot models, which predict a considerably higher brightness temperature of the sunspot umbra at 34 GHz. This inconsistency suggests that the temperature of the region at which the 34 GHz radio emission becomes optically thick should be lower than that predicted by the models.

We estimated the accuracy of coronal magnetic fields derived from radio observations by comparing them to potential field calculations and the DEM measurements using EUV observations. We derived line of sight component of the coronal magnetic field from polarization observations of the thermal bremsstrahlung in the NOAA active region 11150, observed around 3:00 UT on February 3, 2011 using the Nobeyama Radioheliograph at 17 GHz. Because the thermal bremsstrahlung intensity at 17 GHz includes both chromospheric and coronal components, we extracted only the coronal component by measuring the coronal emission measure in EUV observations. In addition, we derived only the radio polarization component of the corona by selecting the region of coronal loops and weak magnetic field strength in the chromosphere along the line of sight. The upper limit of the coronal longitudinal magnetic fields were determined as 100 - 210 G. We also calculated the coronal longitudinal magnetic fields from the potential field extrapolation using the photospheric magnetic field obtained from the Helioseismic and Magnetic Imager (HMI). However, the calculated potential fields were certainly smaller than the observed coronal longitudinal magnetic field. This discrepancy between the potential and the observed magnetic field strengths can be explained consistently by two reasons; (a) the underestimation of the coronal emission measure resulting from the limitation of the temperature range of the EUV observations, (b) the underestimation of the coronal magnetic field resulting from the potential field assumption.

We derived the coronal magnetic field, plasma density, and temperature from the observation of polarization and intensity of radio thermal free-free emission using the Nobeyama Radioheliograph (NoRH) and extreme ultraviolet (EUV) observations. We observed a post-flare loop on the west limb 11 April 2013. The line-of-sight magnetic field was derived from the circularly polarized free-free emission observed by NoRH. The emission measure and temperature were derived from the Atmospheric Imaging Assembly (AIA) onboard Solar Dynamics Observatory (SDO). The derived temperature was used to estimate the emission measure from the NoRH radio free-free emission observations. The derived density from NoRH was larger than that determined using AIA, which can be explained by the fact that the low temperature plasma is not within the temperature coverage of the AIA filters used in this study. We also discuss the other observation of the post-flare loops by the EUV Imager onboard the Solar Terrestrial Relations Observatory (STEREO), which can be used in future studies to reconstruct the coronal magnetic field strength. The derived plasma parameters and magnetic field were used to derive the plasma beta, which is a ratio between the magnetic pressure and the plasma pressure. The derived plasma beta is about 5.7*10^(-4) to 7.6*10^(-4) at the loop top region.

We study the Sunyaev-Zeldovich effect for clusters of galaxies. We explore the relativistic corrections to the Kompaneets equation in terms of two different expansion approximation schemes, namely, the Fokker-Planck expansion approximation and delta function expansion approximation. We show that two expansion approximation formalisms are equivalent under the Thomson approximation, which is extremely good approximation for the CMB photon energies. This will clarify the situation for existing theoretical methods to analyse observation data.

We study the Sunyaev-Zeldovich effect for clusters of galaxies. The Boltzmann equations for the cosmic microwave background photon distribution function are studied in three Lorentz frames. We extend the previous work and derive analytic expressions for the integrated photon redistribution functions over the photon frequency. We also derive analytic expressions in the power series expansion approximation. By combining two formulas, we offer a simple and accurate tool to analyse observation data. These formulas are applicable to the non-thermal electron distributions as well as the standard thermal distribution. The Boltzmann equation is reduced to a single integral form of the electron velocity.

We study the Sunyaev-Zeldovich effect for clusters of galaxies. The Boltzmann equations for the CMB photon distribution function are studied in three Lorentz frames. We clarify the relations of the SZ effects among the different Lorentz frames. We derive analytic expressions for the photon redistribution functions. These formulas are applicable to the nonthermal electron distributions as well as the standard thermal distribution. We show that the Fokker-Planck expansion of the Boltzmann equation can be expanded by the power series of the diffusion operator of the original Kompaneets equation.

We study the inverse Compton scattering of the CMB photons off nonthermal high-energy electrons. In the previous study, assuming the power-law distribution for electrons, we derived the analytic expression for the spectral intensity function $I(\omega)$ in the Thomson approximation, which was applicable up to the photon energies of $\omega <$ O(GeV). In the present paper, we extend the previous work to higher photon energies of $\omega >$ O(GeV) by taking into account the terms dropped in the Thomson approximation, i.e., the Klein-Nishina formula. The analytic expression for $I(\omega)$ is derived with the Klein-Nishina formula. It is shown that $I(\omega)$ has a "knee" structure at $\omega =$ O(PeV). The knee, if exists, should be accessible with gamma-ray observatories such as Fermi-LAT. We propose simple analytical formulae for $I(\omega)$ which are applicable to wide photon energies from Thomson region to extreme Klein-Nishina region.

We have carried out 12CO(J =2-1) and 12CO(J =3-2) observations at spatial resolutions of 1.0-3.8 pc toward the entirety of loops 1 and 2 and part of loop 3 in the Galactic center with NANTEN2 and ASTE. These new results revealed detailed distributions of the molecular gas and the line intensity ratio of the two transitions, R3-2/2-1. In the three loops, R3-2/2-1 is in a range from 0.1 to 2.5 with a peak at ~ 0.7 while that in the disk molecular gas is in a range from 0.1 to 1.2 with a peak at 0.4. This supports that the loops are more highly excited than the disk molecular gas. An LVG analysis of three transitions, 12CO J =3-2 and 2-1 and 13CO J =2-1, toward six positions in loops 1 and 2 shows density and temperature are in a range 102.2 - 104.7 cm-3 and 15-100 K or higher, respectively. Three regions extended by 50-100 pc in the loops tend to have higher excitation conditions as characterized by R3-2/2-1 greater than 1.2. The highest ratio of 2.5 is found in the most developed foot points between loops 1 and 2. This is interpreted that the foot points indicate strongly shocked conditions as inferred from their large linewidths of 50-100 km s-1, confirming the suggestion by Torii et al. (2010b). The other two regions outside the foot points suggest that the molecular gas is heated up by some additional heating mechanisms possibly including magnetic reconnection. A detailed analysis of four foot points have shown a U shape, an L shape or a mirrored-L shape in the b-v distribution. It is shown that a simple kinematical model which incorporates global rotation and expansion of the loops is able to explain these characteristic shapes.

We study a covariant formalism for the Sunyaev-Zeldovich effects developed in the previous papers by the present authors, and derive analytic expressions for the redistribution functions in the Thomson approximation. We also explore another covariant formalism recently developed by Poutanen and Vurm. We show that the two formalisms are mathematically equivalent in the Thomson approximation which is fully valid for the cosmic microwave background photon energies. The present finding will establish a theoretical foundation for the analysis of the Sunyaev-Zeldovich effects for the clusters of galaxies.

Fukui et al. (2006) discovered two huge molecular loops in the Galactic center located in (l, b) ~ (355 deg-359 deg, 0 deg-2 deg) in a large velocity range of -180-40 km s^-1. Following the discovery, we present detailed observational properties of the two loops based on NANTEN 12CO(J=1-0) and 13CO(J=1-0) datasets at 10 pc resolution including a complete set of velocity channel distributions and comparisons with HI and dust emissions as well as with the other broad molecular features. We find new features on smaller scales in the loops including helical distributions in the loop tops and vertical spurs. The loops have counterparts of the HI gas indicating that the loops include atomic gas. The IRAS far infrared emission is also associated with the loops and was used to derive an X-factor of 0.7(+/-0.1){\times}10^20 cm^-2 (K km s^-1)^-1 to convert the 12CO intensity into the total molecular hydrogen column density. From the 12CO, 13CO, H I and dust datasets we estimated the total mass of loops 1 and 2 to be ~1.4 {\times} 106 Msun and ~1.9 {\times} 10^6 Msun, respectively, where the H I mass corresponds to ~10-20% of the total mass and the total kinetic energy of the two loops to be ~10^52 ergs. An analysis of the kinematics of the loops yields that the loops are rotating at ~47 km s-1 and expanding at ~141 km s^-1 at a radius of 670 pc from the center. Fukui et al. (2006) presented a model that the loops are created by the magnetic flotation due to the Parker instability with an estimated magnetic field strength of ~150 {\mu}G. We present comparisons with the recent numerical simulations of the magnetized nuclear disk by Machida et al. (2009) and Takahashi et al. (2009) and show that the theoretical results are in good agreements with the observations. The helical distributions also suggest that some magnetic instability plays a role similarly to the solar helical features.

Fukui et al. (2006) discovered two molecular loops in the Galactic center and argued that the foot points of the molecular loops, two bright spots at both loops ends, represent the gas accumulated by the falling motion along the loops, subsequent to magnetic flotation by the Parker instability. We have carried out sensitive CO observations of the foot points toward l=356 deg at a few pc resolution in the six rotational transitions of CO; 12CO(J=1-0, 3-2, 4-3, 7-6), 13CO(J=1-0) and C18O(J=1-0). The high resolution image of 12CO (J=3-2) has revealed the detailed distribution of the high excitation gas including U shapes, the outer boundary of which shows sharp intensity jumps accompanying strong velocity gradients. An analysis of the multi-J CO transitions shows that the temperature is in a range from 30-100 K and density is around 10^3-10^4 cm^-3, confirming that the foot points have high temperature and density although there is no prominent radiative heating source such as high mass stars in or around the loops. We argue that the high temperature is likely due to the shock heating under C-shock condition caused by the magnetic flotation. We made a comparison of the gas distribution with theoretical numerical simulations and note that the U shape is consistent with numerical simulations. We also find that the region of highest temperature of ~100 K or higher inside the U shape corresponds to the spur having an upward flow, additionally heated up either by magnetic reconnection or bouncing in the interaction with the narrow neck at the bottom of the U shape. We note these new findings further reinforce the magnetic floatation interpretation.

In this paper, results of 2.5-dimensional magnetohydrodynamical simulations are reported for the magnetic reconnection of non-perfectly antiparallel magnetic fields. The magnetic field has a component perpendicular to the computational plane, that is, guide field. The angle theta between magnetic field lines in two half regions is a key parameter in our simulations whereas the initial distribution of the plasma is assumed to be simple; density and pressure are uniform except for the current sheet region. Alfven waves are generated at the reconnection point and propagate along the reconnected field line. The energy fluxes of the Alfven waves and magneto-acoustic waves (slow mode and fast mode) generated by the magnetic reconnection are measured. Each flux shows the similar time evolution independent of theta. The percentage of the energies (time integral of energy fluxes) carried by the Alfven waves and magneto-acoustic waves to the released magnetic energy are calculated. The Alfven waves carry 38.9%, 36.0%, and 29.5% of the released magnetic energy at the maximum (theta=80^\circ) in the case of beta=0.1, 1, and 20 respectively, where beta is the plasma beta (the ratio of gas pressure to magnetic pressure). The magneto-acoustic waves carry 16.2% (theta=70^\circ), 25.9% (theta=60^\circ), and 75.0% (theta=180^\circ) of the energy at the maximum. Implications of these results for solar coronal heating and acceleration of high-speed solar wind are discussed.

We study the inverse Compton scattering of the CMB photons off high-energy nonthermal electrons. We extend the formalism obtained by the previous paper to the case where the electrons have non-zero bulk motions with respect to the CMB frame. Assuming the power-law electron distribution, we find the same scaling law for the probability distribution function P_{1,K}(s) as P_{1}(s) which corresponds to the zero bulk motions, where the peak height and peak position depend only on the power-index parameter. We solved the rate equation analytically. It is found that the spectral intensity function also has the same scaling law. The effect of the bulk motions to the spectral intensity function is found to be small. The present study will be applicable to the analysis of the X-ray and gamma-ray emission models from various astrophysical objects with non-zero bulk motions such as radio galaxies and astrophysical jets.

Based upon the rate equations for the photon distribution function obtained in the previous paper, we study the inverse Compton scattering process for high-energy nonthermal electrons. Assuming the power-law electron distribution, we find a scaling law in the probability distribution function P_1(s), where the peak height and peak position depend only on the power index parameter. We solved the rate equation analytically. It is found that the spectral intensity function also has the scaling law, where the peak height and peak position depend only on the power index parameter. The present study will be particularly important to the analysis of the X-ray and gamma-ray emission models from various astrophysical objects such as radio galaxies and supernova remnants.

We have discovered a molecular dome-like feature towards $355^{\circ} \leq l \leq 359^{\circ}$ and $0^{\circ} \leq b \leq 2^{\circ}$. The large velocity dispersions of 50--100 km s$^{-1}$ of this feature are much larger than those in the Galactic disk and indicate that the feature is located in the Galactic center, probably within $\sim1$ kpc of Sgr A$^{*}$. The distribution has a projected length of $\sim600$ pc and height of $\sim300$ pc from the Galactic disk and shows a large-scale monotonic velocity gradient of $\sim130$ km s $^{-1}$ per $\sim600$ pc. The feature is also associated with HI gas having a more continuous spatial and velocity distribution than that of $^{12}$CO. We interpret the feature as a magnetically floated loop similar to loops 1 and 2 and name it "loop 3". Loop 3 is similar to loops 1 and 2 in its height and length but is different from loops 1 and 2 in that the inner part of loop 3 is filled with molecular emission. We have identified two foot points at the both ends of loop 3. HI, $^{12}$CO and $^{13}$CO datasets were used to estimate the total mass and kinetic energy of loop 3 to be $\sim3.0 \times 10^{6} \Mo$ and $\sim1.7 \times 10^{52}$ ergs. The huge size, velocity dispersions and energy are consistent with the magnetic origin the Parker instability as in case of loops 1 and 2 but is difficult to be explained by multiple stellar explosions. We argue that loop 3 is in an earlier evolutionary phase than loops 1 and 2 based on the inner-filled morphology and the relative weakness of the foot points. This discovery indicates that the western part of the nuclear gas disk of $\sim1$ kpc radius is dominated by the three well-developed magnetically floated loops and suggests that the dynamics of the nuclear gas disk is strongly affected by the magnetic instabilities.

Starting from a covariant formalism of the Sunyaev-Zeldovich effect for the thermal and non-thermal distributions, we derive the frequency redistribution function identical to Wright's method assuming the smallness of the photon energy (in the Thomson limit). We also derive the redistribution function in the covariant formalism in the Thomson limit. We show that two redistribution functions are mathematically equivalent in the Thomson limit which is fully valid for the cosmic microwave background photon energies. We will also extend the formalism to the kinematical Sunyaev-Zeldovich effect. With the present formalism we will clarify the situation for the discrepancy existed in the higher order terms of the kinematical Sunyaev-Zeldovich effect.

Based upon the rate equations for the photon distribution function obtained in the previous paper, we study the formal solutions in three different representation forms for the Sunyaev-Zeldovich effect. By expanding the formal solution in the operator representation in powers of both the derivative operator and electron velocity, we derive a formal solution that is equivalent to the Fokker-Planck expansion approximation. We extend the present formalism to the kinematical Sunyaev-Zeldovich effect. The properties of the frequency redistribution functions are studied. We find that the kinematical Sunyaev-Zeldovich effect is described by the redistribution function related to the electron pressure. We also solve the rate equations numerically. We obtain the exact numerical solutions, which include the full-order terms in powers of the optical depth.

We carried out two-dimensional magnetohydrodynamic simulations of the Galactic gas disk to show that the dense loop-like structures discovered by the Galactic center molecular cloud survey by NANTEN 4 m telescope can be formed by the buoyant rise of magnetic loops due to the Parker instability. At the initial state, we assumed a gravitationally stratified disk consisting of the cool layer ($T \sim 10^3$ K), warm layer ($T \sim 10^4$ K), and hot layer ($T \sim 10^5$ K). Simulation box is a local part of the disk containing the equatorial plane. The gravitational field is approximated by that of a point mass at the galactic center. The self-gravity, and the effects of the galactic rotation are ignored. Numerical results indicate that the length of the magnetic loops emerging from the disk is determined by the scale height of the hot layer ($\sim$ 100 pc at 1 kpc from the Galactic center). The loop length, velocity gradient along the loops and large velocity dispersions at their foot points are consistent with the NANTEN observations. We also show that the loops become top-heavy when the curvature of the loop is sufficiently small, so that the rising loop accumulates the overlying gas faster than sliding it down along the loop. This mechanism is similar to that of the formation of solar chromospheric arch filaments. The molecular loops emerge from the low temperature layer just like the dark filaments observed in the H$\alpha$ image of the emerging flux region of the sun.

A survey for the molecular clouds in the Galaxy with NANTEN mm telescope has discovered molecular loops in the Galactic center region. The loops show monotonic gradients of the line of sight velocity along the loops and the large velocity dispersions towards their foot points. It is suggested that these loops are explained in terms of the buoyant rise of magnetic loops due to the Parker instability. We have carried out global three-dimensional magneto-hydrodynamic simulations of the gas disk in the Galactic center. The gravitational potential is approximated by the axisymmetric potential proposed by Miyamoto & Nagai (1975). At the initial state, we assume a warm (~ 10^4 K) gas torus threaded by azimuthal magnetic fields. Self-gravity and radiative cooling of the gas are ignored. We found that buoyantly rising magnetic loops are formed above the differentially rotating, magnetically turbulent disk. By analyzing the results of global MHD simulations, we have identified individual loops, about 180 in the upper half of the disk, and studied their statistical properties such as their length, width, height, and velocity distributions along the loops. Typical length and height of a loop are 1kpc and 200pc, respectively. The line of sight velocity changes linearly along a loop and shows large dispersions around the foot-points. Numerical results indicate that loops emerge preferentially from the region where magnetic pressure is large. We argue that these properties are consistent with those of the molecular loops discovered by NANTEN.

The heating of the solar chromosphere and corona is a long-standing puzzle in solar physics. Hinode observations show the ubiquitous presence of chromospheric anemone jets outside sunspots in active regions. They are typically 3 to 7 arc seconds = 2000 to 5000 kilometers long and 0.2 to 0.4 arc second = 150 to 300 kilometers wide, and their velocity is 10 to 20 kilometers per second. These small jets have an inverted Y-shape, similar to the shape of x-ray anemone jets in the corona. These features imply that magnetic reconnection similar to that in the corona is occurring at a much smaller spatial scale throughout the chromosphere and suggest that the heating of the solar chromosphere and corona may be related to small-scale ubiquitous reconnection.

The second Born corrections to the electrical and thermal conductivities are calculated for the dense matter in the liquid metal phase for various elemental compositions of astrophysical importance. Inclusion up to the second Born corrections is sufficiently accurate for the Coulomb scattering of the electrons by the atomic nuclei with Z < 26. Our approach is semi-analytical, and is in contrast to that of the previous authors who have used fully numerical values of the cross section for the Coulomb scattering of the electron by the atomic nucleus. The merit of the present semi-analytical approach is that this approach affords us to obtain the results with reliable Z-dependence and \rho-dependence. The previous fully numerical approach has made use of the numerical values of the cross section for the scattering of the electron off the atomic nucleus for a limited number of Z-values, Z=6, 13, 29, 50, 82, and 92, and for a limited number of electron energies, 0.05MeV, 0.1MeV, 0.2MeV, 0.4MeV, 0.7MeV, 1MeV, 2MeV, 4MeV, and 10MeV. Our study, however, has confirmed that the previous results are sufficiently accurate. They are recovered, if the terms higher than the second Born terms are taken into account. We make a detailed comparison of the present results with those of the previous authors. The numerical results are parameterized in a form of analytic formulae that would facilitate practical uses of the results. We also extend our calculations to the case of mixtures of nuclear species. The corresponding subroutine can be retrieved from http://www.ph.sophia.ac.jp/~itoh-ken/subroutine/subroutine.htm

We calculate the screening corrections to the electron capture rates in dense stars by the relativistically degenerate electron liquid. In order to calculate the screening corrections we adopt the linear response theory which is widely used in the field of solid state physics and liquid metal physics. In particular, we use the longitudinal dielectric function for the relativistically degenerate electron liquid derived by Jancovici. We calculate the screening potential at the position of the nucleus. By using this screening potential one can calculate the screening corrections to the electron capture rates. We will present accurate analytic fitting formulae which summarize our numerical results. These fitting formulae will facilitate the application of the present results. The screening corrections to the electron capture rates are typically a few percent.

We calculate the ion-ion correlation effect on the neutrino-nucleus scattering in supernova cores, which is an important opacity source for the neutrinos and plays a vital role in the supernova explosion. In order to calculate the ion-ion correlation effect we use the results of the improved hypernetted-chain method calculations of the classical one-component plasma. As in the preceding studies on this effect, we find a dramatic decrease of the effective neutrino-nucleus scattering cross section for relatively low energy neutrinos with E < 20MeV. As a matter of fact, our calculation shows a much more dramatic reduction of the effective neutrino-nucleus scattering cross section for the low energy neutrinos with E < 10MeV than the results of Horowitz. Therefore, the ion-ion correlation effect will be more important than has hitherto been recognized. We present an accurate analytic fitting formula that summarizes our numerical results. This fitting formula will facilitate the application of the present results to the supernova explosion simulations.

We present accurate analytic fitting formulae which summarize the results of the recent calculation by Nozawa, Itoh and Kohyama on the relativistic thermal bremsstrahlung Gaunt factor for the intracluster plasma. The fitting is carried out for Z = 1 to 28, 6.0 < log T < 8.5, -4.0 < log (omega/kT) < 1.0, where Z is the charge number of the ion, T is the temperature, and omega is the angular frequency of the emitted photon. The present analytic fitting formulae will be useful for the analysis of the X-ray emission which comes from the intracluster plasma as well as the other X-ray sources.

We present an accurate analytic fitting formula for the thermal bremsstrahlung Gaunt factor in the nonrelativistic limit. The fitting formula excellently reproduces the numerical results of the calculation carried out by the present authors using the method of Karzas and Latter. The present analytic fitting formula will be useful for the analysis of the X-ray emission which comes from the intracluster plasmas with relatively low temperatures as well as the other X-ray sources.

We present accurate analytic fitting formulae for the frequency integrated Gaunt factors for the relativistic as well as nonrelativistic thermal bremsstrahlung which is relevant to the intracluster plasma. The fitting is carried out for 1 =< Zj =< 28, 6.0 < log T < 8.5, where Zj is the charge number of the ion, and T is the electron temperature in kelvins. The present analytic fitting formulae will be useful for the analysis of the X-ray emission which comes from the intracluster plasma as well as the other X-ray sources.

With the launch of the Chandra X-Ray Observatory and XMM-Newton, X-ray astronomy has entered an era of precision science. The precision observational data taken by these satellites require precision basic physics data for their analysis. The present authors have recently published a series of papers in which they have reported the results of the accurate calculations of the electron-ion thermal bremsstrahlung which are particularly suited for the analysis of the high-temperature galaxy clusters. In this paper we will present accurate analytic fitting formulae for the nonrelativistic electron-electron thermal bremsstrahlung Gaunt factor which will be useful for the analysis of the X-ray observational data taken by the Chandra X-Ray Observatory and XMM-Newton. With the publication of the present paper X-ray astronomers will be provided with an accurate thermal bremsstrahlung Gaunt factor which has an overall accuracy better than 1% for the temperature range 6.0 <log T[K] < 8.5.

We improve the calculation of Nozawa, Itoh, and Kohyama (1998) to provide a formula for relativistic corrections to the thermal and kinematical Sunyaev-Zeldovich effects that is accurate to fourth order in theta_e = kT_e/m_ec^2, T_e and m_e being the electron temperature and electron mass, respectively. We also carry out a direct numerical integration of the Boltzmann collision term and confirm the excellent accuracy of the analytic formula. This formula will be useful for the analysis of the observational data of the forthcoming experiments of the kinematical Sunyaev-Zeldovich effect for clusters of galaxies.

We extend the formalism of the relativistic thermal and kinematical Sunyaev-Zeldovich effects to the observer's system (the Solar System) moving with a velocity $\beta_{S}$ with respect to the cosmic microwave background radiation. We confirm the results recently obtained by Chluba, Huetsi, and Sunyaev in the lowest order of the observer's velocity $\beta_{S}. We give a more general analytic expression for the thermal and kinematical Sunyaev-Zeldovich effects corresponding to the observer's system with the power series expansion approximation in terms of $\theta_{e}$. It is found that the effect of the motion of the observer on the Sunyaev-Zeldovich effect will become important in the future high precision observation projects of the Sunyaev-Zeldovich effect.

We present an accurate numerical table for the relativistic corrections to the thermal Sunyaev-Zeldovich effect for clusters of galaxies. The numerical results for the relativistic corrections have been obtained by numerical integration of the collision term of the Boltzmann equation. The numerical table is provided for the ranges 0.002 <= theta_e <= 0.100 and 0 <= X <= 20, where theta_e = kT_e/mc^2, X = hbar omega/kT_0, T_e is the electron temperature, omega is the angular frequency of the photon, and T_0 is the temperature of the cosmic microwave background radiation. We also present an accurate analytic fitting formula that reproduces the numerical results with high precision.

We examine r-process nucleosynthesis in a "prompt supernova explosion" from an 8-10 Msun progenitor star, as an alternative scenario to the "neutrino wind" mechanism. In the present model, the progenitor star has formed an oxygen-neon-magnesium core at its center. The core-collapse simulations are performed with a one-dimension, Newtonian hydrodynamic code. We obtain a very weak prompt explosion, in which no r-processing occurs. We further simulate energetic prompt explosions by enhancement of the shock-heating energy, in order to investigate conditions necessary for the production of r-process nuclei in such events. The highly neutronized ejecta (Ye = 0.14-0.20) leads to robust production of r-process nuclei; their relative abundances are in excellent agreement with the solar r-process pattern. Our results suggest that prompt explosions of 8-10 Msun stars with oxygen-neon-magnesium cores can be a promising site of r-process nuclei.

The enhancement factor of the resonant thermonuclear reaction rates is calculated for the extremely dense stellar plasmas in the liquid phase. In order to calculate the enhancement factor we use the screening potential which is deduced from the numerical experiment of the classical one-component plasma. It is found that the enhancement is tremendous for white dwarf densities if the ^{12}C + ^{12}C fusion cross sections show resonant behavior in the astrophysical energy range. We summarize our numerical results by accurate analytic fitting formulae.

The discovery of the second highly $r$-process-enhanced, extremely metal-poor star, CS 31082-001 ([Fe/H] $= -2.9$) has provided a powerful new tool for age determination, by virtue of the detection and measurement of the radioactive species uranium and thorium. One of the serious limitations of this approach, however, is that predictions of the production ratio of U and Th have not been made in the context of a realistic astrophysical model of the $r$-process. We have endeavored to produce such a model, based on the ``neutrino winds'' that are expected to arise from the nascent neutron star of a core-collapse supernova. The mass-integrated $r$-process yields, obtained by assuming a simple time evolution of the neutrino luminosity, are compared to the available spectroscopic elemental abundance data of CS 31082-001. As a result, the ``age'' of this star is determined to be $14.1 \pm 2.5$ Gyr, in excellent agreement with lower limits on the age of the universe estimated by other dating techniques, as well as with other stellar radioactive age estimates. Future measurements of Pt and Pb in this star, as well as expansion of searches for additional $r$-process-enhanced, metal-poor stars (especially those in which both U and Th are measurable), are of special importance to constrain the current astrophysical models for the $r$-process.

We extend the formalism for the calculation of the relativistic corrections to the Sunyaev-Zel'dovich effect for clusters of galaxies and include the multiple scattering effects in the isotropic approximation. We present the results of the calculations by the Fokker-Planck expansion method as well as by the direct numerical integration of the collision term of the Boltzmann equation. The multiple scattering contribution is found to be very small compared with the single scattering contribution. For high-temperature galaxy clusters of kT=15keV, the ratio of the both contributions is -0.2% in the Wien region. In the Rayleigh--Jeans region the ratio is -0.03%. Therefore the multiple scattering contribution is safely neglected for the observed galaxy clusters.

We have succeeded in obtaining a precision analytic fitting formula for the exact numerical results of the relativistic corrections to the thermal Sunyaev-Zeldovich effect for clusters of galaxies which has a 1% accuracy for the crossover frequency region where the thermal thermal Sunyaev-Zeldovich effect signal changes from negative to positive sign. The fitting has been carried out for the ranges 0.020 < theta < 0.035 and 0 < X < 15, where theta= kTe/mc^{2}, X = omega/kT0, Te is the electron temperature, omega is the angular frequency of the photon, and T0 is the temperature of the cosmic microwave background radiation. The overall accuracy of the fitting is better than 0.1%. The present analytic fitting formula will be useful for accurate analyses of the thermal Sunyaev-Zeldovich effect for clusters of galaxies.

We extend the formalism for the calculation of the relativistic corrections to the Sunyaev-Zel'dovich effect for clusters of galaxies and include the multiple scattering effects. We present a systematic method for the inclusion of the multiple scattering effects. The multiple scattering contribution is found to be very small compared with the single scattering contribution. For high-temperature galaxy clusters of k_{B}T_{e} = 15keV, the ratio of the both contributions is -0.3% in the Wien region. In the Rayleigh--Jeans region the ratio is -0.03%. Therefore the multiple scattering contribution is safely neglected for the observed galaxy clusters.

We present an accurate analytic fitting formula for the numerical results for the relativistic corrections to the thermal Sunyaev-Zel'dovich effect for clusters of galaxies. The numerical results for the relativistic corrections have been obtained by numerical integration of the collision term of the Boltzmann equation. The fitting is carried out for the ranges 0.02 < theta_{e} < 0.05 and 0 < X < 20, where theta_{e} = k_{B}T_{e}/m_{e}c^{2}, X = omega/k_{B}T_{0}, T_{e} is the electron temperature, omega is the angular frequency of the photon, and T_{0} is the temperature of the cosmic microwave background radiation. The accuracy of the fitting is generally better than 0.1%. The present analytic fitting formula will be useful for the analyses of the thermal Sunyaev-Zel'dovich effect for high-temperature galaxy clusters.

We extend the formalism of the relativistic thermal and kinematic Sunyaev-Zel'dovich effects and include the polarization of the cosmic microwave background photons. We consider the situation that a cluster of galaxies is moving with a velocity \vec{v} with respect to the cosmic microwave background radiation. In the present formalism, polarization of the scattered cosmic microwave background radiation caused by the proper motion of a cluster of galaxies is naturally derived as a special case of the kinematic Sunyaev-Zel'dovich effect. The relativistic corrections are included also in a natural way. Our results are in complete agreement with the recent results of the relativistic corrections obtained by Challinor, Ford and Lasenby with an entirely different method as well as the nonrelativistic limit obtained by Sunyaev and Zel'dovich. The relativistic correction becomes to be significant in the Wien region.

We extend the formalism of the relativistic thermal Sunyaev-Zel'dovich effect to the system moving with a velocity beta (= v/c) with respect to the cosmic microwave background radiation. In the present formalism, the kinematic Sunyaev-Zel'dovich effect for the cluster of galaxies with a peculiar velocity is derived in a straightforward manner by the Lorentz boost of the generalized Kompaneets equation. We give an analytic expression for the kinematic Sunyaev-Zel'dovich effect which is valid up to O(beta^{2}) with the power series expansion approximation in terms of theta_{e} = k_{B} T_{e}/ mc^{2}, where T_{e} and m are the electron temperature and the electron mass, respectively. It is found that the relativistic corrections to the kinematic Sunyaev-Zel'dovich effect are significant. For a typical electron temperature 10keV, one obtains -8.2% and +1.3% corrections from the O(beta theta_{e}) and O(beta theta_{e}^{2}) contributions, respectively. The O(beta^{2}) correction is extremely small, +0.2% for beta=1/300 at 10keV. Therefore it can be safely neglected. These relativistic corrections are directly reflected on the determination of the peculiar velocity of the cluster of galaxies with the observation of the kinematic Sunyaev-Zel'dovich effect.

We calculate the relativistic thermal bremsstrahlung Gaunt factor for the high-temperature plasma which exists in clusters of galaxies. We calculate the Gaunt factor by employing the Bethe-Heitler cross section corrected by the Elwert factor. The calculations in this paper are made for the fully ionized plasma for the following cases: Z = 10 (Ne), 12 (Mg), 14 (Si), 16 (S), 26 (Fe). We also calculate the Gaunt factor by using the Coulomb-distorted wave functions for nonrelativistic electrons following the method of Karzas and Latter. By comparing the Gaunt factors calculated by these two different methods, we carefully assess the accuracy of the calculation. We present the numerical results in the form of tables.

We calculate the relativistic thermal bremsstrahlung Gaunt factor for the high-temperature plasma which exists in clusters of galaxies. We calculate the Gaunt factor by employing the Bethe-Heitler cross section corrected by the Elwert factor. We also calculate the Gaunt factor by using the Coulomb-distorted wave functions for nonrelativistic electrons following the method of Karzas and Latter. By comparing the Gaunt factors calculated by these two different methods, we carefully assess the accuracy of the calculation. We present the numerical results in the form of tables.

We study the generalized Kompaneets equation (kinetic equation for the photon distribution function taking into account the Compton scattering by electrons) using a relativistically covariant formalism. Using the generalized Kompaneets equation, we derive an analytic expression for the Sunyaev-Zel'dovich effect, which takes into account up to O(theta**5) terms, where theta=kT/mc**2 is the relativistic expansion parameter and T is the electron temperature. We carefully study the applicable region of the obtained analytic expression by comparing with the result of the direct numerical integration. We conclude that the present analytic expression can be reliably applied to the calculation of the Sunyaev-Zel'dovich effect for kT < 15keV, which is the highest electron temperature in the presently known clusters of galaxies. Therefore, the present analytic expression can be applied to all known clusters of galaxies.