We report the detection of a double planetary system around the evolved intermediate-mass star HD 47366 from precise radial-velocity measurements at Okayama Astrophysical Observatory, Xinglong Station, and Australian Astronomical Observatory. The star is a K1 giant with a mass of 1.81+-0.13M_sun, a radius of 7.30+-0.33R_sun, and solar metallicity. The planetary system is composed of two giant planets with minimum mass of 1.75^{+0.20}_{-0.17}Mjup and 1.86^{+0.16}_{-0.15}Mjup, orbital period of 363.3^{+2.5}_{-2.4} d and 684.7^{+5.0}_{-4.9} d, and eccentricity of 0.089^{+0.079}_{-0.060} and 0.278^{+0.067}_{-0.094}, respectively, which are derived by a double Keplerian orbital fit to the radial-velocity data. The system adds to the population of multi-giant-planet systems with relatively small orbital separations, which are preferentially found around evolved intermediate-mass stars. Dynamical stability analysis for the system revealed, however, that the best-fit orbits are unstable in the case of a prograde configuration. The system could be stable if the planets were in 2:1 mean-motion resonance, but this is less likely considering the observed period ratio and eccentricity. A present possible scenario for the system is that both of the planets have nearly circular orbits, namely the eccentricity of the outer planet is less than ~0.15, which is just within 1.4sigma of the best-fit value, or the planets are in a mutually retrograde configuration with a mutual orbital inclination larger than 160 degree.

We report the discovery of a substellar companion around the intermediatemass giant HD 175679. Precise radial velocity data of the star from Xinglong Station and Okayama Astrophysical Observatory (OAO) revealed a Keplerian velocity variation with an orbital period of 1366.8 \pm 5.7 days, a semiamplitude of 380.2 \pm 3.2m s.1, and an eccentricity of 0.378 \pm 0.008. Adopting a stellar mass of 2.7 \pm 0.3 M\odot, we obtain the minimum mass of the HD 175679 b is 37.3 \pm 2.8 MJ, and the semimajor axis is 3.36 \pm 0.12 AU. This discovery is the second brown dwarf companion candidate from a joint planet-search program between China and Japan.

We construct a new sample of 32 obscured active galactic nuclei (AGNs) selected from the Second XMM-Newton Serendipitous Source Catalogue to investigate their multiwavelength properties in relation to the "scattering fraction", the ratio of the soft X-ray flux to the absorption-corrected direct emission. The sample covers a broad range of the scattering fraction (0.1%-10%). A quarter of the 32 AGNs have a very low scattering fraction (smaller than 0.5%), which suggests that they are buried in a geometrically thick torus with a very small opening angle. We investigate correlations between the scattering fraction and multiwavelength properties. We find that AGNs with a small scattering fraction tend to have low [O III]lambda5007/X-ray luminosity ratios. This result agrees with the expectation that the extent of the narrow-line region is small because of the small opening angle of the torus. There is no significant correlation between scattering fraction and far-infrared luminosity. This implies that a scale height of the torus is not primarily determined by starburst activity. We also compare scattering fraction with black hole mass or Eddington ratio and find a weak anti-correlation between the Eddington ratio and scattering fraction. This implies that more rapidly growing supermassive black holes tend to have thicker tori.

We present the results of X-ray spectral analysis of 22 active galactic nuclei (AGNs) with a small scattering fraction selected from the Second XMM-Newton Serendipitous Source Catalogue using hardness ratios. They are candidates of buried AGNs, since a scattering fraction, which is a fraction of scattered emission by the circumnuclear photoionized gas with respect to direct emission, can be used to estimate the size of the opening part of an obscuring torus. Their X-ray spectra are modeled by a combination of a power law with a photon index of 1.5-2 absorbed by a column density of 10^23-24 cm^-2, an unabsorbed power law, narrow Gaussian lines, and some additional soft components. We find that scattering fractions of 20 among 22 objects are less than a typical value (3%) for Seyfert2s observed so far. In particular, those of eight objects are smaller than 0.5%, which are in the range for buried AGNs found in recent hard X-ray surveys. Moreover, [O III] lambda5007 luminosities at given X-ray luminosities for some objects are smaller than those for Seyfert2s previously known. This fact could be interpreted as a smaller size of optical narrow emission line regions produced in the opening direction of the obscuring torus. These results indicate that they are strong candidates for the AGN buried in a very geometrically thick torus.

We derive analytic formulas for the radiation power output when electrons are accelerated by a relativistic comoving kinetic Poynting flux, and validate these analytic results with Particle-In-Cell simulations. We also derive analytically the critical frequency of the radiation spectrum. Potential astrophysical applications of these results are discussed. A quantitative model of gamma-ray bursts based on the breakout of kinetic Poynting flux is presented.

Balmer and Paschen continuum emission as well as Balmer series lines of P Cygni-type profile from H_gamma through H_23 are revealed in the violet spectra of BM Gem, a carbon star associated with an oxygen-rich circumstellar shell (`silicate carbon star') observed with the high dispersion spectrograph (HDS) on the Subaru telescope. The blue-shifted absorption in the Balmer lines indicates the presence of an outflow, the line of sight velocity of which is at least 400 km s^-1, which is the highest outflow velocity observed to date in a carbon star. We argue that the observed unusual features in BM Gem are strong evidence for the presence of a companion, which should form an accretion disk that gives rise to both an ionized gas region and a high velocity, variable outflow. The estimated luminosity of ~0.2 (0.03-0.6) L_sun for the ionized gas can be maintained by a mass accretion rate to a dwarf companion of ~10^-8 M_sun yr^-1, while ~10^-10 M_sun yr^-1 is sufficient for accretion to a white dwarf companion. These accretion rates are feasible for some detached binary configurations on the basis of the Bond-Hoyle type accretion process. We concluded that the carbon star BM Gem is in a detached binary system with a companion of low mass and low luminosity. However, we are unable to determine whether this companion object is a dwarf or a white dwarf. The upper limits for binary separation are 210 AU and 930 AU for a dwarf and a white dwarf, respectively. We also note that the observed features of BM Gem mimic those of Mira (omi Cet), which may suggest actual similarities in their binary configurations and circumstellar structures.

Stellar parameters and abundances of Na, Mg, Al, K, Ca, Sr, Ba, and Eu are determined for four very metal-poor stars (-2.66 < [Fe/H] < -2.15) based on non-LTE line formation and analysis of high-resolution (R ~60000 and 90000) high signal-to-noise (S/N > 200) observed spectra. A model atom for H I is presented. An effective temperature was obtained from the Balmer Halpha and Hbeta line wing fits, the surface gravity from the Hipparcos parallax if available and the non-LTE ionization balance between Ca I and Ca II. Based on the hyperfine structure affecting the Ba II resonance line, the fractional abundance of the odd isotopes of Ba was derived for HD 84937 and HD 122563 from a requirement that Ba abundances from the resonance line and subordinate lines of Ba II must be equal. For each star, non-LTE leads to a consistency of Teff from two Balmer lines and to a higher temperature compared to the LTE case, by up to 60 K. Non-LTE effects are important in spectroscopic determination of surface gravity from Ca I/Ca II. For each star with a known trigonometric gravity, non-LTE abundances from the lines of two ionization stages agree within the error bars, while a difference in the LTE abundances consists of 0.23 dex to 0.40 dex for different stars. Departures from LTE are found to be significant for the investigated atoms, and they strongly depend on stellar parameters. For HD 84937, the Eu/Ba ratio is consistent with the relative solar system r-process abundances, and the fraction of the odd isotopes of Ba, f_odd, equals 0.43+-0.14. The latter can serve as a constraint on r-process models. The lower Eu/Ba ratio and f_odd = 0.22+-0.15 found for HD 122563 suggest that the s-process or the unknown process has contributed significantly to the Ba abundance in this star.

We derive analytic formulas for the power output and critical frequency of radiation by electrons accelerated by relativistic kinetic Poynting flux, and validate these results with Particle-In-Cell plasma simulations. We find that the in-situ radiation power output and critical frequency are much below those predicted by the classical synchrotron formulae. We discuss potential astrophysical applications of these results.

We report the detection of a substellar companion orbiting the intermediate-mass giant star 11 Com (G8 III). Precise Doppler measurements of the star from Xinglong station and Okayama Astrophysical Observatory (OAO) revealed Keplerian velocity variations with an orbital period of 326.03 +/- 0.32 days, a semiamplitude of 302.8 +/- 2.6 m/s, and an eccentricity of 0.231 +/- 0.005. Adopting a stellar mass of 2.7 +/- 0.3 M_solar, the minimum mass of the companion is 19.4 +/- 1.5 M_Jup, well above the deuterium burning limit, and the semimajor axis is 1.29 +/- 0.05 AU. This is the first result from the joint planet search program between China and Japan aiming at revealing statistics of substellar companions around intermediate-mass giants. 11 Com b emerged from 300 targets of the planet search program at OAO. The current detection rate of a brown dwarf candidate seems to be comparable to that around solar-type stars within orbital separations of $\sim$3 AU.

We summarize latest PIC simulation results on the radiation from Poynting jets and strongly magnetized collisionless shocks when a Poynting jet runs into cold ambient medium. We find that in all cases the radiative power output is much below that predicted by synchrotron radiation and the critical frequency is also much lower than the synchrotron critical frequency. We discuss the implications for the interpretation of GRB prompt emission data.

We report the Compton scattering emission from the Poynting flux acceleration of electron- positron plasma simulated by the 2-1/2 dimensional particle-in-cell(PIC) code. We show these and other remarkable properties of Poynting flux acceleration and Compton spectral output, and discuss the agreement with the observed spectra of GRBs and XRFs.

We review recent developments in particle acceleration by Poynting flux using plasma kinetic simulations, and discuss their potential applications to gamma-ray burst phenomenology

We review PIC simulation results of GRB emissions for both the Poynting flux and shock scenarios

Plasma outflows from gamma-ray bursts (GRB), pulsar winds, relativistic jets, and ultra-intense laser targets radiate high energy photons. However, radiation damping is ignored in conventional PIC simulations. In this letter, we study the radiation damping effect on particle acceleration via Poynting fluxes in two-and-half-dimensional particle-in-cell (PIC) plasma simulation of electron-positron plasmas. Radiation damping force is self-consistently calculated for each particle and reduces the net acceleration force. The emitted radiation is peaked within a few degrees from the direction of Poynting flux and strongly linear-polarized.

We study the radiation damping effect on the relativistic acceleration of electron-positron plasmas with two-and-half-dimensional particle-in-cell (PIC) simulation. Particles are accelerated by Poynting flux via the diamagnetic relativistic pulse accelerator (DRPA), and decelerated by the self-consistently solved radiation damping force. With $\Omega_{ce}/\omega_{pe}\geq 10$, the Lorentz factor of the highest energy particles reaches gamma>100, and the acceleration still continues. The emitted radiation is peaked within few degrees from the direction of Poynting flux and strongly linearly polarized, which may be detectable in gamma-ray burst(GRB) observations. We also show that the DRPA is insensitive to the initial supporting currents.

The radiation damping effect on the diamagnetic relativistic pulse accelerator (DRPA) is studied in two-and-half dimensional Particle-in-Cell (PIC) simulation with magnetized electron-positron plasmas. Self-consistently solved radiation damping force converts particle energy to radiation energy. The DRPA is still robust with radiation, and the Lorentz factor of the most high energy particles reach more than two thousand before they decouple from the electromagnetic pulse. Resulted emitted power from the pulse front is lower in the radiative case than the estimation from the non-radiative case due to the radiation damping. The emitted radiation is strongly linearly polarized and peaked within few degrees from the direction of Poynting flux.

We have obtained high-resolution (R $\simeq$ 50,000 or 90,000), high-quality (S/N $\ga$ 100) spectra of 22 very metal-poor stars ([Fe/H] $\la$ -2.5) with the High Dispersion Spectrograph fabricated for the 8.2m Subaru Telescope. The spectra cover the wavelength range from 3500 to 5100 {\AA}; equivalent widths are measured for isolated lines of numerous elemental species, including the $\alpha$ elements, the iron-peak elements, and the light and heavy neutron-capture elements. Errors in the measurements and comparisons with previous studies are discussed. These data will be used to perform detailed abundance analyses in the following papers of this series. Radial velocities are also reported, and are compared with previous studies. At least one moderately r-process-enhanced metal-poor star, HD 186478, exhibits evidence of a small-amplitude radial velocity variation, confirming the binary status noted previously. During the course of this initial program, we have discovered a new moderately r-process-enhanced, very metal-poor star, CS 30306-132 ([Fe/H] $= -2.4$; [Eu/Fe] $= +0.85$), which is discussed in detail in the companion paper.

All experiments, which have been proposed so far to model the magnetorotational instability (MRI) in the laboratory, involve a Couette flow of liquid metals in a rotating annulus. All liquid metals have small magnetic Prandtl numbers, Pm, of about 10^{-6} (the ratio of kinematic viscosity to magnetic diffusivity). With plasmas both large and small Pm are achievable by varying the temperature and the density of plasma. Compressibility and fast rotation of the plasma result in radial stratification of the equilibrium plasma density. Evolution of perturbations in radially stratified viscous and resistive plasma permeated by an axial uniform magnetic field is considered. The differential rotation of the plasma is induced by the ExB drift in applied radial electric field. Global unstable eigenmodes are calculated by our newly developed matrix code. The plasma is shown to be MRI unstable for parameters easily achievable in experimental setup.

Despite the importance of the magnetorotational instability (MRI) as a fundamental mechanism for angular momentum transport in magnetized accretion disks, it has yet to be demonstrated in the laboratory. A liquid sodium alpha-omega dynamo experiment at the New Mexico Institute of Mining and Technology provides an ideal environment to study the MRI in a rotating metal annulus (Couette flow). A local stability analysis is performed as a function of shear, magnetic field strength, magnetic Reynolds number, and turbulent Prandtl number. The later takes into account the minimum turbulence induced by the formation of an Ekman layer against the rigidly rotating end walls of a cylindrical vessel. Stability conditions are presented and unstable conditions for the sodium experiment are compared with another proposed MRI experiment with liquid gallium. Due to the relatively large magnetic Reynolds number achievable in the sodium experiment, it should be possible to observe the excitation of the MRI for a wide range of wavenumbers and further to observe the transition to the turbulent state.

We have obtained a high resolution (R~90,000) spectrum of the carbon- and s-process-element-rich, very metal-poor ([Fe/H]=-2.7) subgiant LP625-44, as well as that of HD140283 (a metal-poor subgiant with normal abundance ratio) for comparison, with the High Dispersion Spectrograph (HDS) for the Subaru Telescope for detailed abundance study. The oxygen abundance derived from the O I triplet around 7770A is uncertain, but the excess of oxygen in LP625-44 seems remarkable (perhaps by nearly a factor 10), in comparison with that of HD140283 derived from the same lines. The Na enhancement in LP625-44 is by about a factor 50, suggesting hydrogen burning in the 22Ne-rich layer in an asymptotic giant branch star which produces the abundance pattern of this object. In our new spectrum of LP625-44, the Pb I lambda 3683A line has been detected, as well as the Pb I lambda 4057A line which has already been studied, confirming the Pb abundance (log epsilon(Pb)~1.9) derived by the previous work. The abundance ratio of s-process elements at the second peak (e.g., La, Ce and Nd) to that at the third peak (Pb) in LP 625-44 is significantly higher (by a factor 5) than that in other three s-process element-rich objects recently studied by van Eck et al.. Recent theoretical works have modeled the s-process nucleosynthesis in the radiative layer of asymptotic giant branch stars in the inter-pulse phase, and the above results means that these processes produced a large scatter in the abundance ratios. Another possiblity is that different processes (e.g., s-process nucleosynthesis during thermal pulses) have contributed to heavy elements in the early Galaxy.

The stability of nonaxisymmetric perturbations in differentially rotating astrophysical accretion disks is analyzed by fully incorporating the properties of shear flows. We verify the presence of discrete unstable eigenmodes with complex and pure imaginary eigenvalues, without any artificial disk edge boundaries, unlike Ogilvie & Pringle(1996)'s claim. By developing the mathematical theory of a non-self-adjoint system, we investigate the nonlocal behavior of eigenmodes in the vicinity of Alfven singularities at omega_D=omega_A, where omega_D is the Doppler-shifted wave frequency and omega_A=k_// v_A is the Alfven frequency. The structure of the spectrum of discrete eigenmodes is discussed and the magnetic field and wavenumber dependence of the growth rate are obtained. Exponentially growing modes are present even in a region where the local dispersion relation theory claims to have stable eigenvalues. The velocity field created by an eigenmode is obtained, which explains the anomalous angular momentum transport in the nonlinear stage of this stability.