Recent Developments in Parallelization of the Multidimensional Integration Package DICE

DICE is a general purpose multidimensional numerical integration package. There can be two ways in the parallelization of DICE, "distributing random numbers into workers" and "distributing hypercubes into workers". Furthermore, there can be the combination of both ways. So far, we had developed the parallelization code using the former way and reported it in ACAT2002 in Moscow. Here, we will present the recent developments of parallelized DICE in the latter way as the 2nd stage of our parallelization activities.Comment: 5 pages, 2 figures, Talk given at the X International Workshop on Advanced Computing and Analysis Techniques in Physics Research, ACAT 2005, DESY-Zeuthen, Germany, 22-27 May 200

Monte Carlo integration on GPU

We use a graphics processing unit (GPU) for fast computations of Monte Carlo integrations. Two widely used Monte Carlo integration programs, VEGAS and BASES, are parallelized on GPU. By using $W^{+}$ plus multi-gluon production processes at LHC, we test integrated cross sections and execution time for programs in FORTRAN and C on CPU and those on GPU. Integrated results agree with each other within statistical errors. Execution time of programs on GPU run about 50 times faster than those in C, and more than 60 times faster than the original FORTRAN programs.Comment: 6 pages, 2 figure

A massive star origin for an unusual helium-rich supernova in an elliptical galaxy

The unusual helium-rich (type Ib) supernova SN 2005E is distinguished from any supernova hitherto observed by its faint and rapidly fading light curve, prominent calcium lines in late-phase spectra and lack of any mark of recent star formation near the supernova location. These properties are claimed to be explained by a helium detonation in a thin surface layer of an accreting white dwarf (Perets et al. 2010). Here we report on observations of SN 2005cz appeared in an elliptical galaxy, whose observed properties resemble those of SN 2005E in that it is helium-rich and unusually faint, fades rapidly, shows much weaker oxygen emission lines than those of calcium in the well-evolved spectrum. We argue that these properties are best explained by a core-collapse supernova at the low-mass end ($8-12 M_{\odot}$) of the range of massive stars that explode (Smartt 2009). Such a low mass progenitor had lost its hydrogen-rich envelope through binary interaction, having very thin oxygen-rich and silicon-rich layers above the collapsing core, thus ejecting a very small amount of radioactive $^{56}$Ni and oxygen. Although the host galaxy NGC 4589 is an elliptical, some studies have revealed evidence of recent star-formation activity (Zhang et al. 2008), consistent with the core-collapse scenario.Comment: Accepted by Nature (24 March 2010), 32 pages including Supplementary Informatio

Optical I-band Linear Polarimetry of the Magnetar 4U 0142+61 with Subaru

The magnetar 4U~0142+61 has been well studied at optical and infrared wavelengths and is known to have a complicated broad-band spectrum over the wavelength range. Here we report the result from our linear imaging polarimetry of the magnetar at optical $I$-band. From the polarimetric observation carried out with the 8.2-m Subaru telescope, we determine the degree of linear polarization $P=1.0\pm$3.4\%, or $P\leq$5.6\% (90\% confidence level). Considering models suggested for optical emission from magnetars, we discuss the implications of our result. The upper limit measurement indicates that different from radio pulsars, magnetars probably would not have strongly polarized optical emission if the emission arises from their magnetosphere as suggested.Comment: 5 pages, 1 figure, accepted for publication on Ap

Large Polarization Degree of Comet 2P/Encke Continuum Based on Spectropolarimetric Signals During Its 2017 Apparition

Spectropolarimetry is a powerful technique for investigating the physical properties of gas and solid materials in cometary comae without mutual contamination, but there have been few spectropolarimetric studies to extract each component. We attempt to derive the continuum polarization degree of comet 2P/Encke, free from influence of molecular emissions. The target is unique in that it has an orbit dynamically decoupled from Jupiter like main-belt asteroids, while ejecting gas and dust like ordinary comets. We observed the comet using the Higashi-Hiroshima Optical and Near-Infrared Camera attached to the Cassegrain focus of the 150-cm Kanata telescope on UT 2017 February 21 when the comet was at the solar phase angle of 75.7 deg. We find that the continuum polarization degree with respect to the scattering plane is 33.8+/-2.7 % at the effective wavelength of 0.815 um, which is significantly higher than those of cometary dust in a high-Pmax group at similar phase angles. Assuming that an ensemble polarimetric response of 2P/Encke's dust as a function of phase angle is morphologically similar with those of other comets, its maximum polarization degree is estimated to > 40 % at the phase angle of ~100 deg. In addition, we obtain the polarization degrees of the C2 swan bands (0.51-0.56 um), the NH2 alpha bands (0.62-0.69 um) and the CN-red system (0.78-0.94 um) in a range of 3-19 %, which depend on the molecular species and rotational quantum numbers of each branch. The polarization vector aligns nearly perpendicularly to the scattering plane with the average of 0.4 deg over a wavelength range of 0.50-0.97 um. From the observational evidence, we conjecture that the large polarization degree of 2P/Encke would be attributable to a dominance of large dust particles around the nucleus, which have remained after frequent perihelion passages near the Sun.Comment: 9 pages, 4 figures, accepted for publication in Astronomy & Astrophysic

Computed Tomography Appearance of Hem-O-Lok Clips in Patients Who Have Undergone Laparoscopic Nephrectomy or Nephroureterectomy

Hem-O-Lok clips are radiopaque on computed tomography in patients who have undergone laparoscopic nephrectomy and nephroureterectomy

Calculation of HELAS amplitudes for QCD processes using graphics processing unit (GPU)

We use a graphics processing unit (GPU) for fast calculations of helicity amplitudes of quark and gluon scattering processes in massless QCD. New HEGET ({\bf H}ELAS {\bf E}valuation with {\bf G}PU {\bf E}nhanced {\bf T}echnology) codes for gluon self-interactions are introduced, and a C++ program to convert the MadGraph generated FORTRAN codes into HEGET codes in CUDA (a C-platform for general purpose computing on GPU) is created. Because of the proliferation of the number of Feynman diagrams and the number of independent color amplitudes, the maximum number of final state jets we can evaluate on a GPU is limited to 4 for pure gluon processes ($gg\to 4g$), or 5 for processes with one or more quark lines such as $q\bar{q}\to 5g$ and $qq\to qq+3g$. Compared with the usual CPU-based programs, we obtain 60-100 times better performance on the GPU, except for 5-jet production processes and the $gg\to 4g$ processes for which the GPU gain over the CPU is about 20