We present the first spectroscopic measurements of the ATLAS and Aliqa Uma streams from the Southern Stellar Stream Spectroscopic Survey ($S^5$), in combination with the photometric data from the Dark Energy Survey and astrometric data from $Gaia$. From the coherence of spectroscopic members in radial velocity and proper motion, we find out that these two systems are physically one stream with discontinuity in morphology and density on the sky (the "kink" feature). We refer to this entire stream as ATLAS-Aliqa Uma stream, or AAU stream. We perform a comprehensive exploration of the effect of baryonic substructures and find that only an encounter with the Sagittarius dwarf $\sim 0.5$ Gyr ago can create a feature similar to the observed "kink". In addition, we also identify two gaps in the ATLAS component associated with the broadening in the stream width (the "broadening" feature). These gaps have likely been created by small mass perturbers, such as dark matter halos, as AAU stream is the most distant cold stream known with severe variations in both the stream surface density and the stream track on the sky. With the stream track, stream distance and kinematic information, we determine the orbit of the AAU stream and find that it has been affected by the Large Magellanic Cloud, resulting in a misalignment between the proper motion and stream track. Together with the Orphan-Chenab Stream, AAU is the second stream pair that has been found to be a single stream separated into two segments by external perturbation.

We report the discovery of J1521-3538, a bright (V=12.2), very metal-poor ([Fe/H]=-2.8) strongly r-process enhanced field horizontal branch star, based on a high-resolution, high signal-to-noise Magellan/MIKE spectrum. J1521-3538 shows the largest r-process element over-abundance in any known r-process-enhanced star, with [Eu/Fe]=+2.2, and its chemical abundances of 22 neutron-capture elements closely match the scaled solar r-process pattern. J1521-3538 is also one of few known carbon-enhanced metal-poor stars with r-process enhancement (CEMP-r stars), as found after correcting the measured C abundance for the star's evolutionary status. We propose to extend the existing classification of moderately enhanced (+0.3<=[Eu/Fe]<=+1.0) r-I and strongly r-process enhanced ([Eu/Fe]>+1.0) r-II stars to include an r-III class, for r-process stars such as J1521-3538, with [Eu/Fe]>+2.0 and [Ba/Eu]<-0.5, or >100 times the solar ratio of europium to iron. Using cosmochronometry, we estimate J1521-3538 to be 12.5+-5 Gyr and 8.9+-5 Gyr old, using two different sets of initial production ratios. These ages are based on measurements of the Th line at 4019 A and other r-process element abundances. This is broadly consistent with the old age of a low-mass metal-poor field red horizontal branch star. J1521-3538 likely originated in a low-mass dwarf galaxy that was later accreted by the Milky Way, as evidenced by its highly eccentric orbit.

We provide a statistical mechanical derivation of relativistic magnetohydrodynamics on the basis of the $(3+1)$-dimensional quantum electrodynamics; the system endowed with the magnetic one-form symmetry. The conservation laws and the constitutive relations are presented in a manifestly covariant way with respect to the general coordinate transformation. The method of the local Gibbs ensemble (or nonequilibrium statistical operator) combined with the path-integral formula for the thermodynamic functional enables us to obtain an exact form of the constitutive relations. Applying the derivative expansion to the exact formula, we derive the first-order constitutive relations for the relativistic magnetohydrodynamics. The result for the QED plasma preserving the parity and charge-conjugation symmetries is equipped with two electrical resistivities and five (three bulk and two shear) viscosities. We also show that those transport coefficients satisfy the Onsager's reciprocal relation and a set of inequalities, indicating the semi-positivity of the entropy production rate consistent with the local second law of thermodynamics.

We search for dynamical substructures in the LAMOST DR3 very metal-poor (VMP) star catalog. After cross-matching with Gaia DR2, there are 3300 VMP stars with available high-quality astrometric information that have halo-like kinematics. We apply a method based on self-organizing maps to find groups clustered in the 4D space of orbital energy and angular momentum. We identify 57 dynamically tagged groups, which we label DTG-1 to DTG-57. Most of them belong to existing substructures in the nearby halo, such as the $Gaia$ Sausage or Sequoia. The stream identified by Helmi et al. is recovered, but the two disjoint portions of the substructure have distinct dynamical properties. The very retrograde substructure Rg5 found previously by Myeong et al. is also retrieved. We report 6 new DTGs with highly retrograde orbits, 2 with very prograde orbits, and 12 with polar orbits. By mapping other datasets (APOGEE halo stars, and catalogs of r-process-enhanced and CEMP stars) onto the trained neuron map, we can associate stars with detailed chemical abundances to the DTGs, and look for associations with chemically peculiar stars. The highly eccentric $Gaia$ Sausage groups contain representatives both of debris from the satellite itself (which is $\alpha$-poor) and the Splashed Disk, sent up into eccentric halo orbits from the encounter (and is $\alpha$-rich). The new prograde substructures also appear to be associated with the Splashed Disk. The DTGs belonging to the $Gaia$ Sausage host two relatively metal-rich $r$-II stars and six CEMP stars in different sub-classes, consistent with the idea that the $Gaia$ Sausage progenitor is a massive dwarf galaxy. Rg5 is dynamically associated with two highly $r$-process-enhanced stars with [Fe/H] $\sim -$3. This finding indicates that its progenitor might be an ultra-faint dwarf galaxy that has experienced $r$-process enrichment from neutron star mergers.

The rotation curve of the Milky Way is commonly used to estimate the local dark matter density $\rho_{{\rm DM},\odot}$. However, the estimates are subject to the choice of the distribution of baryons needed in this type of studies. In this work we explore several Galactic mass models that differ in the distribution of baryons and dark matter, in order to determine $\rho_{{\rm DM},\odot}$. For this purpose we analyze the precise circular velocity curve measurement of the Milky Way up to $\sim 25$ kpc from the Galactic centre obtained from Gaia DR2 [1]. We find that the estimated value of $\rho_{{\rm DM},\odot}$ stays robust to reasonable changes in the spherical dark matter halo. However, we show that $\rho_{{\rm DM},\odot}$ is affected by the choice of the model for the underlying baryonic components. In particular, we find that $\rho_{{\rm DM},\odot}$ is mostly sensitive to uncertainties in the disk components of the Galaxy. We also show that, when choosing one particular baryonic model, the estimate of $\rho_{{\rm DM},\odot}$ has an uncertainty of only about $10\%$ of its best-fit value, but this uncertainty gets much bigger when we also consider the variation of the baryonic model. In particular, the rotation curve method does not allow to exclude the presence of an additional very thin component, that can increase $\rho_{{\rm DM},\odot}$ by more than a factor of 8 (the thin disk could even be made of dark matter). Therefore, we conclude that exclusively using the rotation curve of the Galaxy is not enough to provide a robust estimate of $\rho_{{\rm DM},\odot}$. For all the models that we study without the presence of an additional thin component, our resulting estimates of the local dark matter density take values in the range $\rho_{{\rm DM},\odot} \simeq \text{0.3--0.4}\,\mathrm{GeV/cm^3}$, consistent with many of the estimates in the literature.

A recently discovered young, high-velocity giant star J01020100-7122208 is a good candidate of hypervelocity star ejected from the Galactic center, although it has a bound orbit. If we assume that this star was ejected from the Galactic center, it can be used to constrain the Galactic potential, because the deviation of its orbit from a purely radial orbit informs us of the torque that this star has received after its ejection. Based on this assumption, we estimate the flattening of the dark matter halo of the Milky Way by using the Gaia DR2 data and the circular velocity data from Eilers et al. (2019). Our Bayesian analysis shows that the orbit of J01020100-7122208 favors a prolate dark matter halo within $\sim$ 10 kpc from the Galactic center. The posterior distribution of the density flattening $q$ shows a broad distribution at $q \gtrsim1$ and peaks at $q \simeq 1.5$. Also, 98.5\% of the posterior distribution is located at $q>1$, highly disfavoring an oblate halo.

The Carrington event is considered to be one of the most extreme space weather events in observational history within a series of magnetic storms caused by extreme interplanetary coronal mass ejections (ICMEs) from a large and complex active region (AR) emerged on the solar disk. In this article, we study the temporal and spatial evolutions of the source sunspot active region and visual aurorae, and compare this storm with other extreme space weather events on the basis of their spatial evolution. Sunspot drawings by Schwabe, Secchi, and Carrington describe the position and morphology of the source AR at that time. Visual auroral reports from the Russian Empire, Iberia, Ireland, Oceania, and Japan fill the spatial gap of auroral visibility and revise the time series of auroral visibility in mid to low magnetic latitudes (MLATs). The reconstructed time series is compared with magnetic measurements and shows the correspondence between low to mid latitude aurorae and the phase of magnetic storms. The spatial evolution of the auroral oval is compared with those of other extreme space weather events in 1872, 1909, 1921, and 1989 as well as their storm intensity, and contextualizes the Carrington event, as one of the most extreme space weather events, but likely not unique.

Large sunspots can be observed with the unaided eye under suitable atmospheric seeing conditions. Such observations are of particular value because the frequency of their appearance provides an approximate indication of the prevailing level of solar activity. Unaided-eye sunspot (UES) observations can be traced back well before the start of telescopic observations of the Sun, especially in the East Asian historical records. It is therefore important to compare more modern, UES observations with the results of telescopic sunspot observations, to gain a better understanding of the nature of the UES records. A previous comparison of Chinese UES records and Greenwich photo-heliographic results between 1874 and 1918 indicated that a few of the UES were apparently not supported by direct photographic evidence of at least one sunspot with a large area. This article reveals that one of the Chinese unaided-eye observations had possibly captured the transit of Venus on 9 December 1874. The Chinese sunspot records on this date are compared with Western sunspot observations on the same day. It is concluded that sunspots on the solar disk were quite small and the transit of Venus was probably misinterpreted as a sunspot by the Chinese local scholars. This case indicates that sunspots or comparable "obscuring" objects with an area as large as 1000 millionths of the solar disk could easily have been seen with the unaided eye under suitable seeing conditions. It also confirms the visibility of sunspots near the solar limb with the unaided eye. This study provides an explanation of the apparent discrepancy between the Chinese UES observation on 9 December 1874 and the Western sunspot observations using telescopes, as well as a basis for further discussion on the negative pairs in 1900 and 1911, apparently without sufficiently large area.

Although knowing the occurrence frequency of severe space weather events is important for a modern society, it is insufficiently known due to the lack of magnetic or sunspot observations, before the Carrington event in 1859 known as one of the largest events during the last two centuries. Here, we show that a severe magnetic storm occurred on 8 March 1582 based on auroral records in East Asia. The equatorward boundary of auroral visibility reached 28.8{\deg} magnetic latitude. The equatorward boundary of the auroral oval is estimated to be 33.0{\deg} invariant latitude (ILAT), which is comparable to the storms on 25/26 September 1909 (~31.6{\deg} ILAT, minimum Dst of -595 nT), 28/29 August 1859 (~36.5{\deg} ILAT), and 13/14 March 1989 (~40{\deg} ILAT, minimum Dst of -589 nT). Assuming that the equatorward boundary is a proxy for the scale of magnetic storms, we presume that the storm on March 1582 was severe. We also found that the storm on March 1582 lasted, at least, for three days by combining European records. The auroral oval stayed at mid-latitude for the first two days and moved to low-latitude (in East Asia) for the last day. It is plausible that the storm was caused by a series of ICMEs (interplanetary coronal mass ejections). We can reasonably speculate that a first ICME could have cleaned up interplanetary space to make the following ICMEs more geo-effective, as probably occurred in the Carrington and Halloween storms.

We report that LAMOST-HVS1 is a massive hyper-runaway subgiant star with mass of 8.3 Msun and super-Solar metallicity, ejected from the inner stellar disk of the Milky Way $\sim$ 33 Myr ago with the intrinsic ejection velocity of $568^{+19}_{-17}$ km/s (corrected for the streaming motion of the disk), based on the proper motion data from Gaia Data Release 2 (DR2) and high-resolution spectroscopy. The extremely large ejection velocity indicates that this star was not ejected by the supernova explosion of the binary companion. Rather, it was probably ejected by a 3- or 4-body dynamical interaction with more massive objects in a high-density environment. Such a high-density environment may be attained at the core region of a young massive cluster with mass of $\gtrsim 10^4$ Msun. The ejection agent that took part in the ejection of LAMOST-HVS1 may well be an intermediate mass black hole ($\gtrsim$ 100 Msun), a very massive star ($\gtrsim$ 100 Msun), or multiple ordinary massive stars ($\gtrsim$ 30 Msun). Based on the flight time and the ejection location of LAMOST-HVS1, we argue that its ejection agent or its natal star cluster is currently located near the Norma spiral arm. The natal star cluster of LAMOST-HVS1 may be an undiscovered young massive cluster near the Norma spiral arm.

The origin of the Hercules stream, the most prominent velocity substructure in the Solar neighbour disc stars, is still under debate. Recent accurate measurements of position, velocity, and metallicity provided by Tycho Gaia Astrometric Solution (TGAS) and RAdial Velocity Experiments (RAVE) have revealed that the Hercules stream is most clearly seen in the metal-rich region ([Fe/H] > 0), while it is not clearly seen in lower metallicity region ([Fe/H] < -0.25). By using a large number of chemo-dynamical 2D test-particle simulations with a rotating bar and/or spiral arms, we find that the observed [Fe/H] dependence of the Hercules stream is a natural consequence of the inside-out formation of the stellar disc and the existence of highly non-closed orbits in the rotating frame of the bar or spiral arms. Our successful models that reproduce the observed properties of the Hercules stream include not only fast-bar-only and fast-bar+spiral models, but also slow-bar+spiral models. This indicates that it is very difficult to estimate the pattern speed of the bar or spiral arms based only on the observations of the Hercules stream in the Solar neighbourhood. As a by-product of our simulations, we make some predictions about the locations across the Galactic plane where we can observe velocity bimodality that is not associated with the Hercules stream. These predictions can be tested by the Gaia Data Release 2, and such a test will improve our understanding of the evolution of the Milky Way stellar disc.

Gravitational 3-body interaction among binary stars and the supermassive black hole (SMBH) at the center of the Milky Way occasionally ejects a hypervelocity star (HVS) with a velocity of ~1000 km/s. Due to the ejection location, such a HVS initially has negligible azimuthal angular momentum Lz ~ 0 kpc km/s. Even if the halo is mildly triaxial, Lz of a recently ejected nearby HVS remains negligible, since its flight time from the Galactic Center is too short to accumulate noticeable torque. However, if we make a wrong assumption about the Solar position and velocity, such a HVS would apparently have noticeable non-zero azimuthal angular momentum, due to the wrong reflex motion of the Sun. Conversely, with precise astrometric data for a nearby HVS, we can measure the Solar position and velocity by assuming that the HVS has zero azimuthal angular momentum. Based on this idea, here we propose a method to estimate the Galactocentric distance of the Sun R0 and the Galactocentric Solar azimuthal velocity Vsun by using a HVS. We demonstrate with mock data for a nearby HVS candidate that the Gaia astrometric data, along with the currently available constraint on Vsun/R0 from the proper motion measurement of Sgr A*, can constrain R0 and Vsun with uncertainties of ~0.27 kpc and ~7.8 km/s (or fractional uncertainties of 3 percent), respectively. Our method will be a promising tool to constrain (R0, Vsun), given that Gaia is expected to discover many nearby HVSs in the near future.

We report the discovery of 30 stars with extreme space velocities ($>$ 480 km/s) in the Gaia-DR2 archive. These stars are a subset of 1743 stars with high-precision parallax, large tangential velocity ($v_{tan}>$ 300 km/s), and measured line-of-sight velocity in DR2. By tracing the orbits of the stars back in time, we find at least one of them is consistent with having been ejected by the supermassive black hole at the Galactic Center. Another star has an orbit that passed near the Large Magellanic Cloud (LMC) about 200 Myr ago. Unlike previously discovered blue hypervelocity stars, our sample is metal-poor (-1.5 $<$ [Fe/H] $<$ -1.0) and quite old ($>$ 1 Gyr). We discuss possible mechanisms for accelerating old stars to such extreme velocities. The high observed space density of this population, relative to potential acceleration mechanisms, implies that these stars are probably bound to the Milky Way (MW). If they are bound, the discovery of this population would require a local escape speed of around $\sim$ 600 km/s and consequently imply a virial mass of $M_{200} \sim 1.4 \times 10^{12} M_\odot$ for the MW.

We present kinematics of 35 highly r-process-enhanced ([Eu/Fe] >= +0.7) metal-poor (-3.8 < [Fe/H] < -1.4) field stars. We calculate six-dimensional positions and velocities, evaluate energies and integrals of motion, and compute orbits for each of these stars using parallaxes and proper motions from the second Gaia data release and published radial velocities. All of these stars have halo kinematics. Most stars (66%) remain in the inner regions of the halo (< 13 kpc), and many (51%) have orbits that pass within 2.6 kpc of the Galactic center. Several stars (20%) have orbits that extend beyond 20 kpc, including one with an orbital apocenter larger than the Milky Way virial radius. We apply three clustering methods to search for structure in phase space, and we identify eight groups. No abundances are considered in the clustering process, but the [Fe/H] dispersions of the groups are smaller than would be expected by random chance. The orbital properties, clustering in phase space and metallicity, and lack of highly r-process-enhanced stars on disk-like orbits indicate that such stars likely were accreted from disrupted satellites. Comparison with the galaxy luminosity-metallicity relation suggests M_V > -9 for most of the progenitor satellites, characteristic of ultra-faint or low-luminosity classical dwarf spheroidal galaxies. Environments with low rates of star formation and Fe production, rather than the nature of the r-process site, may be key to obtaining the [Eu/Fe] ratios found in highly r-process-enhanced stars.

Aims. Historical records provide evidence of extreme magnetic storms with equatorward auroral extensions before the epoch of systematic magnetic observations. One significant magnetic storm occurred on February 15, 1730. We scale this magnetic storm with auroral extension and contextualise it based on contemporary solar activity. Methods. We examined historical records in East Asia and computed the magnetic latitude (MLAT) of observational sites to scale magnetic storms. We also compared them with auroral records in Southern Europe. We examined contemporary sunspot observations to reconstruct detailed solar activity between 1729 and 1731. Results. We show 29 auroral records in East Asian historical documents and 37 sunspot observations. Conclusions. These records show that the auroral displays were visible at least down to 25.8{\deg} MLAT throughout East Asia. In comparison with contemporary European records, we show that the boundary of the auroral display closest to the equator surpassed 45.1{\deg} MLAT and possibly came down to 31.5{\deg} MLAT in its maximum phase, with considerable brightness. Contemporary sunspot records show an active phase in the first half of 1730 during the declining phase of the solar cycle. This magnetic storm was at least as intense as the magnetic storm in 1989, but less intense than the Carrington event.

We report a measurement of the B-mode polarization power spectrum in the cosmic microwave background (CMB) using the POLARBEAR experiment in Chile. The faint B-mode polarization signature carries information about the Universe's entire history of gravitational structure formation, and the cosmic inflation that may have occurred in the very early Universe. Our measurement covers the angular multipole range 500 < l < 2100 and is based on observations of an effective sky area of 25 square degrees with 3.5 arcmin resolution at 150 GHz. On these angular scales, gravitational lensing of the CMB by intervening structure in the Universe is expected to be the dominant source of B-mode polarization. Including both systematic and statistical uncertainties, the hypothesis of no B-mode polarization power from gravitational lensing is rejected at 97.1% confidence. The band powers are consistent with the standard cosmological model. Fitting a single lensing amplitude parameter A_BB to the measured band powers, A_BB = 1.12 +/- 0.61 (stat) +0.04/-0.12 (sys) +/- 0.07 (multi), where A_BB = 1 is the fiducial WMAP-9 LCDM value. In this expression, "stat" refers to the statistical uncertainty, "sys" to the systematic uncertainty associated with possible biases from the instrument and astrophysical foregrounds, and "multi" to the calibration uncertainties that have a multiplicative effect on the measured amplitude A_BB.

The study of historical great geomagnetic storms is crucial for assessing the possible risks to the technological infrastructure of a modern society, caused by extreme space-weather events. The normal benchmark has been the great geomagnetic storm of September 1859, the so-called "Carrington Event". However, there are numerous records of another great geomagnetic storm in February 1872. This storm, about 12 years after the Carrington Event, resulted in comparable magnetic disturbances and auroral displays over large areas of the Earth. We have revisited this great geomagnetic storm in terms of the auroral and sunspot records in the historical documents from East Asia. In particular, we have surveyed the auroral records from East Asia and estimated the equatorward boundary of the auroral oval to be near 24.3 deg invariant latitude (ILAT), on the basis that the aurora was seen near the zenith at Shanghai (20 deg magnetic latitude, MLAT). These results confirm that this geomagnetic storm of February 1872 was as extreme as the Carrington Event, at least in terms of the equatorward motion of the auroral oval. Indeed, our results support the interpretation of the simultaneous auroral observations made at Bombay (10 deg MLAT). The East Asian auroral records have indicated extreme brightness, suggesting unusual precipitation of high-intensity, low-energy electrons during this geomagnetic storm. We have compared the duration of the East Asian auroral displays with magnetic observations in Bombay and found that the auroral displays occurred in the initial phase, main phase, and early recovery phase of the magnetic storm.

Inflation is the leading theory of the first instant of the universe. Inflation, which postulates that the universe underwent a period of rapid expansion an instant after its birth, provides convincing explanation for cosmological observations. Recent advancements in detector technology have opened opportunities to explore primordial gravitational waves generated by the inflation through B-mode (divergent-free) polarization pattern embedded in the Cosmic Microwave Background anisotropies. If detected, these signals would provide strong evidence for inflation, point to the correct model for inflation, and open a window to physics at ultra-high energies. LiteBIRD is a satellite mission with a goal of detecting degree-and-larger-angular-scale B-mode polarization. LiteBIRD will observe at the second Lagrange point with a 400 mm diameter telescope and 2,622 detectors. It will survey the entire sky with 15 frequency bands from 40 to 400 GHz to measure and subtract foregrounds. The U.S. LiteBIRD team is proposing to deliver sub-Kelvin instruments that include detectors and readout electronics. A lenslet-coupled sinuous antenna array will cover low-frequency bands (40 GHz to 235 GHz) with four frequency arrangements of trichroic pixels. An orthomode-transducer-coupled corrugated horn array will cover high-frequency bands (280 GHz to 402 GHz) with three types of single frequency detectors. The detectors will be made with Transition Edge Sensor (TES) bolometers cooled to a 100 milli-Kelvin base temperature by an adiabatic demagnetization refrigerator.The TES bolometers will be read out using digital frequency multiplexing with Superconducting QUantum Interference Device (SQUID) amplifiers. Up to 78 bolometers will be multiplexed with a single SQUID amplidier. We report on the sub-Kelvin instrument design and ongoing developments for the LiteBIRD mission.

The velocity anisotropy parameter, beta, is a measure of the kinematic state of orbits in the stellar halo which holds promise for constraining the merger history of the Milky Way (MW). We determine global trends for beta as a function of radius from three suites of simulations, including accretion only and cosmological hydrodynamic simulations. We find that both types of simulations are consistent and predict strong radial anisotropy (<beta>~0.7) for Galactocentric radii greater than 10 kpc. Previous observations of beta for the MW's stellar halo claim a detection of an isotropic or tangential "dip" at r~20 kpc. Using the N-body+SPH simulations, we investigate the temporal persistence, population origin, and severity of "dips" in beta. We find dips in the in situ stellar halo are long-lived, while dips in the accreted stellar halo are short-lived and tied to the recent accretion of satellite material. We also find that a major merger as early as z~1 can result in a present day low (isotropic to tangential) value of beta over a wide range of radii and angular expanse. While all of these mechanisms are plausible drivers for the beta dip observed in the MW, in the simulations, each mechanism has a unique metallicity signature associated with it, implying that future spectroscopic surveys could distinguish between them. Since an accurate knowledge of beta(r) is required for measuring the mass of the MW halo, we note significant transient dips in beta could cause an overestimate of the halo's mass when using spherical Jeans equation modeling.

Dim red aurora at low magnetic latitudes is a visual and recognized manifestation of geomagnetic storms. The great low-latitude auroral displays seen throughout East Asia on 16-18 September 1770 are considered to manifest one of the greatest storms. Recently found 111 historical documents in East Asia attest that these low-latitude auroral displays were succeeding for almost 9 nights during 10-19 September 1770 in the lowest magnetic latitude areas (< 30{\deg}). This suggests that the duration of the great magnetic storm is much longer than usual. Sunspot drawings from 1770 reveals the fact that sunspots area was twice as large as those observed in another great storm of 1859, which substantiates this unusual storm activities in 1770. These spots likely ejected several huge, sequential magnetic structures in short duration into interplanetary space, resulting in spectacular world-wide aurorae in mid-September 1770. These findings provide new insights about the history, duration, and effects of extreme magnetic storms that may be valuable for those who need to mitigate against extreme events.

A continuously rotating half-wave plate (CRHWP) is a promising tool to improve the sensitivity to large angular scales in cosmic microwave background (CMB) polarization measurements. With a CRHWP, single detectors can measure three of the Stokes parameters, $I$, $Q$ and $U$, thereby avoiding the set of systematic errors that can be introduced by mismatches in the properties of orthogonal detector pairs. We focus on the implementation of CRHWPs in large aperture telescopes (i.e. the primary mirror is larger than the current maximum half-wave plate diameter of $\sim$0.5 m), where the CRHWP can be placed between the primary mirror and focal plane. In this configuration, one needs to address the intensity to polarization ($I{\rightarrow}P$) leakage of the optics, which becomes a source of 1/f noise and also causes differential gain systematics that arise from CMB temperature fluctuations. In this paper, we present the performance of a CRHWP installed in the POLARBEAR experiment, which employs a Gregorian telescope with a 2.5 m primary illumination pattern. The CRHWP is placed near the prime focus between the primary and secondary mirrors. We find that the $I{\rightarrow}P$ leakage is larger than the expectation from the physical properties of our primary mirror, resulting in a 1/f knee of 100 mHz. The excess leakage could be due to imperfections in the detector system, i.e. detector non-linearity in the responsivity and time-constant. We demonstrate, however, that by subtracting the leakage correlated with the intensity signal, the 1/f noise knee frequency is reduced to 32 mHz ($\ell \sim$39 for our scan strategy), which is very promising to probe the primordial B-mode signal. We also discuss methods for further noise subtraction in future projects where the precise temperature control of instrumental components and the leakage reduction will play a key role.

Determining the velocity distribution of halo stars is essential for estimating the mass of the Milky Way and for inferring its formation history. Since the stellar halo is a dynamically hot system, the velocity distribution of halo stars is well described by the 3-dimensional velocity dispersions $(\sigma_r, \sigma_\theta, \sigma_\phi)$, or by the velocity anisotropy parameter $\beta=1-(\sigma_\theta^2+\sigma_\phi^2)/(2\sigma_r^2)$. Direct measurements of $(\sigma_r, \sigma_\theta, \sigma_\phi)$ consistently suggest $\beta =0.5$-$0.7$ for nearby halo stars. In contrast, the value of $\beta$ at large Galactocentric radius $r$ is still controversial, since reliable proper motion data are available for only a handful of stars. In the last decade, several authors have tried to estimate $\beta$ for distant halo stars by fitting the observed line-of-sight velocities at each radius with simple velocity distribution models (local fitting methods). Some results of local fitting methods imply $\beta<0$ at $r \gtrsim 20 \;\rm{kpc}$, which is inconsistent with recent predictions from cosmological simulations. Here we perform mock-catalogue analyses to show that the estimates of $\beta$ based on local fitting methods are reliable only at $r \leq 15 \;\rm{kpc}$ with the current sample size ($\sim10^3$ stars at a given radius). As $r$ increases, the line-of-sight velocity (corrected for the Solar reflex motion) becomes increasingly closer to the Galactocentric radial velocity, so that it becomes increasingly more difficult to estimate tangential velocity dispersion $(\sigma_\theta, \sigma_\phi)$ from line-of-sight velocity distribution. Our results suggest that the forthcoming Gaia data will be crucial for understanding the velocity distribution of halo stars at $r \geq 20\;\rm{kpc}$.

POLARBEAR-2 (PB-2) is a cosmic microwave background (CMB) polarization experiment that will be located in the Atacama highland in Chile at an altitude of 5200 m. Its science goals are to measure the CMB polarization signals originating from both primordial gravitational waves and weak lensing. PB-2 is designed to measure the tensor to scalar ratio, r, with precision {\sigma}(r) < 0.01, and the sum of neutrino masses, {\Sigma}m{\nu}, with {\sigma}({\Sigma}m{\nu}) < 90 meV. To achieve these goals, PB-2 will employ 7588 transition-edge sensor bolometers at 95 GHz and 150 GHz, which will be operated at the base temperature of 250 mK. Science observations will begin in 2017.

We constrain anisotropic cosmic birefringence using four-point correlations of even-parity $E$-mode and odd-parity $B$-mode polarization in the cosmic microwave background measurements made by the POLARization of the Background Radiation (POLARBEAR) experiment in its first season of observations. We find that the anisotropic cosmic birefringence signal from any parity-violating processes is consistent with zero. The Faraday rotation from anisotropic cosmic birefringence can be compared with the equivalent quantity generated by primordial magnetic fields if they existed. The POLARBEAR nondetection translates into a 95% confidence level (C.L.) upper limit of 93 nanogauss (nG) on the amplitude of an equivalent primordial magnetic field inclusive of systematic uncertainties. This four-point correlation constraint on Faraday rotation is about 15 times tighter than the upper limit of 1380 nG inferred from constraining the contribution of Faraday rotation to two-point correlations of $B$-modes measured by Planck in 2015. Metric perturbations sourced by primordial magnetic fields would also contribute to the $B$-mode power spectrum. Using the POLARBEAR measurements of the $B$-mode power spectrum (two-point correlation), we set a 95% C.L. upper limit of 3.9 nG on primordial magnetic fields assuming a flat prior on the field amplitude. This limit is comparable to what was found in the Planck 2015 two-point correlation analysis with both temperature and polarization. We perform a set of systematic error tests and find no evidence for contamination. This work marks the first time that anisotropic cosmic birefringence or primordial magnetic fields have been constrained from the ground at subdegree scales.

The readout of transition-edge sensor (TES) bolometers with a large multiplexing factor is key for the next generation Cosmic Microwave Background (CMB) experiment, Polarbear-2, having 7,588 TES bolometers. To enable the large arrays, we have been developing a readout system with a multiplexing factor of 40 in the frequency domain. Extending that architecture to 40 bolometers requires an increase in the bandwidth of the SQUID electronics above 4 MHz. This paper focuses on cryogenic readout and shows how it affects cross talk and the responsivity of the TES bolometers. A series resistance, such as equivalent series resistance (ESR) of capacitors for LC filters, leads to non-linear response of the bolometers. A wiring inductance modulates a voltage across the bolometers and causes cross talk. They should be controlled well to reduce systematic errors in CMB observations. We have been developing a cryogenic readout with a low series impedance and have tuned bolometers in the middle of their transition at a high frequency (> 3 MHz).

We present an overview of the design and status of the \Pb-2 and the Simons Array experiments. \Pb-2 is a Cosmic Microwave Background polarimetry experiment which aims to characterize the arc-minute angular scale B-mode signal from weak gravitational lensing and search for the degree angular scale B-mode signal from inflationary gravitational waves. The receiver has a 365~mm diameter focal plane cooled to 270~milli-Kelvin. The focal plane is filled with 7,588 dichroic lenslet-antenna coupled polarization sensitive Transition Edge Sensor (TES) bolometric pixels that are sensitive to 95~GHz and 150~GHz bands simultaneously. The TES bolometers are read-out by SQUIDs with 40 channel frequency domain multiplexing. Refractive optical elements are made with high purity alumina to achieve high optical throughput. The receiver is designed to achieve noise equivalent temperature of 5.8~$\mu$K$_{CMB}\sqrt{s}$ in each frequency band. \Pb-2 will deploy in 2016 in the Atacama desert in Chile. The Simons Array is a project to further increase sensitivity by deploying three \Pb-2 type receivers. The Simons Array will cover 95~GHz, 150~GHz and 220~GHz frequency bands for foreground control. The Simons Array will be able to constrain tensor-to-scalar ratio and sum of neutrino masses to $\sigma(r) = 6\times 10^{-3}$ at $r = 0.1$ and $\sum m_\nu (\sigma =1)$ to 40 meV.

The dynamics of stellar streams in rotating barred potentials is explained for the first time. Naturally, neighbouring stream stars reach pericentre at slightly different times. In the presence of a rotating bar, these neighbouring stream stars experience different bar orientations during pericentric passage and hence each star receives a different torque from the bar. These differing torques reshape the angular momentum and energy distribution of stars in the stream, which in turn changes the growth rate of the stream. For a progenitor orbiting in the same sense as the bar's rotation and satisfying a resonance condition, the resultant stream can be substantially shorter or longer than expected, depending on whether the pericentric passages of the progenitor occur along the bar's minor or major axis respectively. We present a full discussion of this phenomenon focusing mainly on streams confined to the Galactic plane. In stark contrast with the evolution in static potentials, which give rise to streams that grow steadily in time, rotating barred potentials can produce dynamically old, short streams. This challenges the traditional viewpoint that the inner halo consists of well phase-mixed material whilst the tidally-disrupted structures in the outer halo are more spatially coherent. We argue that this mechanism plays an important role in explaining the mysteriously short Ophiuchus stream that was recently discovered near the bulge region of the Milky Way.

Atmosphere is one of the most important noise sources for ground-based cosmic microwave background (CMB) experiments. By increasing optical loading on the detectors, it amplifies their effective noise, while its fluctuations introduce spatial and temporal correlations between detected signals. We present a physically motivated 3d-model of the atmosphere total intensity emission in the millimeter and sub-millimeter wavelengths. We derive a new analytical estimate for the correlation between detectors time-ordered data as a function of the instrument and survey design, as well as several atmospheric parameters such as wind, relative humidity, temperature and turbulence characteristics. Using an original numerical computation, we examine the effect of each physical parameter on the correlations in the time series of a given experiment. We then use a parametric-likelihood approach to validate the modeling and estimate atmosphere parameters from the POLARBEAR-I project first season data set. We derive a new 1.0% upper limit on the linear polarization fraction of atmospheric emission. We also compare our results to previous studies and weather station measurements. The proposed model can be used for realistic simulations of future ground-based CMB observations.

We discuss properties of photons in extremely strong magnetic fields induced by the relativistic heavy-ion collisions. We investigate the vacuum birefringence, the real-photon decay, and the photon splitting which are all forbidden in the ordinary vacuum, but become possible in strong magnetic fields. These effects potentially give rise to anisotropies in photon and dilepton spectra.

We explored a method to reconstruct the distribution function of the Galactic thick disc within the action space where nearby thick-disc stars are distributed. By applying this method to 127 chemically-selected thick-disc stars in the Solar neighbourhood, we found that the vertical velocity dispersion that corresponds to the reconstructed distribution function declines approximately as $\exp (-R/R_s)$ at 4.5 kpc < R < 9.5 kpc, with $R_s$ = 8.3 $\pm$ 1.1 (rand.) $\pm$ 1.6 (sys.) kpc. Also, we found that the vertical velocity dispersion $\sigma_z$ of our local thick-disc stars shows only weak dependency on radial and azimuthal velocities $(v_R, v_\phi)$. We discuss possible implications of these results on the global structure of the Milky Way thick disc.

POLARBEAR-2 is a next-generation receiver for precision measurements of the polarization of the cosmic microwave background (Cosmic Microwave Background (CMB)). Scheduled to deploy in early 2015, it will observe alongside the existing POLARBEAR-1 receiver, on a new telescope in the Simons Array on Cerro Toco in the Atacama desert of Chile. For increased sensitivity, it will feature a larger area focal plane, with a total of 7,588 polarization sensitive antenna-coupled Transition Edge Sensor (TES) bolometers, with a design sensitivity of 4.1 uKrt(s). The focal plane will be cooled to 250 milliKelvin, and the bolometers will be read-out with 40x frequency domain multiplexing, with 36 optical bolometers on a single SQUID amplifier, along with 2 dark bolometers and 2 calibration resistors. To increase the multiplexing factor from 8x for POLARBEAR-1 to 40x for POLARBEAR-2 requires additional bandwidth for SQUID readout and well-defined frequency channel spacing. Extending to these higher frequencies requires new components and design for the LC filters which define channel spacing. The LC filters are cold resonant circuits with an inductor and capacitor in series with each bolometer, and stray inductance in the wiring and equivalent series resistance from the capacitors can affect bolometer operation. We present results from characterizing these new readout components. Integration of the readout system is being done first on a small scale, to ensure that the readout system does not affect bolometer sensitivity or stability, and to validate the overall system before expansion into the full receiver. We present the status of readout integration, and the initial results and status of components for the full array.

We describe the design of a new polarization sensitive receiver, SPT-3G, for the 10-meter South Pole Telescope (SPT). The SPT-3G receiver will deliver a factor of ~20 improvement in mapping speed over the current receiver, SPTpol. The sensitivity of the SPT-3G receiver will enable the advance from statistical detection of B-mode polarization anisotropy power to high signal-to-noise measurements of the individual modes, i.e., maps. This will lead to precise (~0.06 eV) constraints on the sum of neutrino masses with the potential to directly address the neutrino mass hierarchy. It will allow a separation of the lensing and inflationary B-mode power spectra, improving constraints on the amplitude and shape of the primordial signal, either through SPT-3G data alone or in combination with BICEP-2/KECK, which is observing the same area of sky. The measurement of small-scale temperature anisotropy will provide new constraints on the epoch of reionization. Additional science from the SPT-3G survey will be significantly enhanced by the synergy with the ongoing optical Dark Energy Survey (DES), including: a 1% constraint on the bias of optical tracers of large-scale structure, a measurement of the differential Doppler signal from pairs of galaxy clusters that will test General Relativity on ~200 Mpc scales, and improved cosmological constraints from the abundance of clusters of galaxies.

Gravitational lensing due to the large-scale distribution of matter in the cosmos distorts the primordial Cosmic Microwave Background (CMB) and thereby induces new, small-scale $B$-mode polarization. This signal carries detailed information about the distribution of all the gravitating matter between the observer and CMB last scattering surface. We report the first direct evidence for polarization lensing based on purely CMB information, from using the four-point correlations of even- and odd-parity $E$- and $B$-mode polarization mapped over $\sim30$ square degrees of the sky measured by the POLARBEAR experiment. These data were analyzed using a blind analysis framework and checked for spurious systematic contamination using null tests and simulations. Evidence for the signal of polarization lensing and lensing $B$-modes is found at 4.2$\sigma$ (stat.+sys.) significance. The amplitude of matter fluctuations is measured with a precision of $27\%$, and is found to be consistent with the Lambda Cold Dark Matter ($\Lambda$CDM) cosmological model. This measurement demonstrates a new technique, capable of mapping all gravitating matter in the Universe, sensitive to the sum of neutrino masses, and essential for cleaning the lensing $B$-mode signal in searches for primordial gravitational waves.

We reconstruct the gravitational lensing convergence signal from Cosmic Microwave Background (CMB) polarization data taken by the POLARBEAR experiment and cross-correlate it with Cosmic Infrared Background (CIB) maps from the Herschel satellite. From the cross-spectra, we obtain evidence for gravitational lensing of the CMB polarization at a statistical significance of 4.0$\sigma$ and evidence for the presence of a lensing $B$-mode signal at a significance of 2.3$\sigma$. We demonstrate that our results are not biased by instrumental and astrophysical systematic errors by performing null-tests, checks with simulated and real data, and analytical calculations. This measurement of polarization lensing, made via the robust cross-correlation channel, not only reinforces POLARBEAR auto-correlation measurements, but also represents one of the early steps towards establishing CMB polarization lensing as a powerful new probe of cosmology and astrophysics.

The entire evolution of the Milky Way, including its mass-assembly and star-formation history, is imprinted onto the chemo-dynamical distribution function of its member stars, f(x, v, [X/H]), in the multi-dimensional phase space spanned by position, velocity, and elemental abundance ratios. In particular, the chemo-dynamical distribution functions for low-mass stars (e.g., G- or K-type dwarfs) are precious tracers of the earliest stages of the Milky Way's formation, since their main-sequence lifetimes approach or exceed the age of the universe. A basic tenet of essentially all previous analyses is that the stellar metallicity, usually parametrized as [Fe/H], is conserved over time for main-sequence stars (at least those that have not been polluted due to mass transfer from binary companions). If this holds true, any correlations between metallicity and kinematics for long-lived main-sequence stars of different masses, effective temperatures, or spectral types must strictly be the same, since they reflect the same mass-assembly and star-formation histories. By analyzing a sample of nearby metal-poor halo and thick-disk stars on the main sequence, taken from Data Release 8 of the Sloan Digital Sky Survey, we find that the median metallicity of G-type dwarfs is systematically higher (by about 0.2 dex) than that of K-type dwarfs having the same median rotational velocity about the Galactic center. If it can be confirmed, this finding may invalidate the long-accepted assumption that the atmospheric metallicities of long-lived stars are conserved over time.

LiteBIRD is a next-generation satellite mission to measure the polarization of the cosmic microwave background (CMB) radiation. On large angular scales the B-mode polarization of the CMB carries the imprint of primordial gravitational waves, and its precise measurement would provide a powerful probe of the epoch of inflation. The goal of LiteBIRD is to achieve a measurement of the characterizing tensor to scalar ratio $r$ to an uncertainty of $\delta r=0.001$. In order to achieve this goal we will employ a kilo-pixel superconducting detector array on a cryogenically cooled sub-Kelvin focal plane with an optical system at a temperature of 4~K. We are currently considering two detector array options; transition edge sensor (TES) bolometers and microwave kinetic inductance detectors (MKID). In this paper we give an overview of LiteBIRD and describe a TES-based polarimeter designed to achieve the target sensitivity of 2~$\mu$K$\cdot$arcmin over the frequency range 50 to 320~GHz.

The POLARBEAR-2 CosmicMicrowave Background (CMB) experiment aims to observe B-mode polarization with high sensitivity to explore gravitational lensing of CMB and inflationary gravitational waves. POLARBEAR-2 is an upgraded experiment based on POLARBEAR-1, which had first light in January 2012. For POLARBEAR-2, we will build a receiver that has 7,588 Transition Edge Sensor (TES) bolometers coupled to two-band (95 and 150 GHz) polarization-sensitive antennas. For the large array's readout, we employ digital frequency-domain multiplexing and multiplex 32 bolometers through a single superconducting quantum interference device (SQUID). An 8-bolometer frequency-domain multiplexing readout has been deployed on POLARBEAR-1 experiment. Extending that architecture to 32 bolometers requires an increase in the bandwidth of the SQUID electronics to 3 MHz. To achieve this increase in bandwidth, we use Digital Active Nulling (DAN) on the digital frequency multiplexing platform. In this paper, we present requirements and improvements on parasitic inductance and resistance of cryogenic wiring and capacitors used for modulating bolometers. These components are problematic above 1 MHz. We also show that our system is able to bias a bolometer in its superconducting transition at 3 MHz.

We investigate decay and production of neutral pions in strong magnetic fields. In the presence of strong magnetic fields, transition between a neutral pion and a virtual photon becomes possible through the triangle diagram relevant for the chiral anomaly. We find that the decay mode of a neutral pion into two photons cannot persist in the dominant mode in strong magnetic fields, and that decay into a dilepton instead dominates over the other modes. We also investigate the effects of magnetic fields on prompt virtual photons created in ultrarelativistic heavy-ion collisions. There is no anisotropy in the spectrum at the stage of creation of prompt virtual photons, but after traversing the strong magnetic field that is induced perpendicularly to the reaction plane, virtual photons turn into neutral pions, leading to an anisotropic spectrum of dileptons as a feasible signature in the measurement.

We compute the refractive indices of a photon propagating in strong magnetic fields on the basis of the analytic representation of the vacuum polarization tensor obtained in our previous paper. When the external magnetic field is strong enough for the fermion one-loop diagram of the polarization tensor to be approximated by the lowest Landau level, the propagating mode in parallel to the magnetic field is subject to modification: The refractive index deviates from unity and can be very large, and when the photon energy is large enough, the refractive index acquires an imaginary part indicating decay of a photon into a fermion-antifermion pair. We study dependences of the refractive index on the propagating angle and the magnetic-field strength. It is also emphasized that a self-consistent treatment of the equation which defines the refractive index is indispensable for accurate description of the refractive index. This self-consistent treatment physically corresponds to consistently including the effects of back reactions of the distorted Dirac sea in response to the incident photon.

The orbital motions of halo stars in the Milky Way reflect the orbital motions of the progenitor systems in which they formed, making it possible to trace the mass-assembly history of the Galaxy. Direct measurement of three-dimensional velocities, based on accurate proper motions and line-of-sight velocities, has revealed that the majority of halo stars in the inner-halo region move on eccentric orbits. However, our understanding of the motions of distant, in-situ halo-star samples is still limited, due to the lack of accurate proper motions for these stars. Here we explore a model-independent analysis of the line-of-sight velocities and spatial distribution of a recent sample of 1865 carefully selected halo blue horizontal-branch (BHB) stars within 30 kpc of the Galactic center. We find that the mean rotational velocity of the very metal-poor ([Fe/H] < -2.0) BHB stars significantly lags behind that of the relatively more metal-rich ([Fe/H] > -2.0) BHB stars. We also find that the relatively more metal-rich BHB stars are dominated by stars with eccentric orbits, as previously observed for other stellar samples in the inner-halo region. By contrast, the very metal-poor BHB stars are dominated by stars on rounder, lower-eccentricity orbits. Our results indicate that the motion of the progenitor systems of the Milky Way that contributed to the stellar populations found within 30 kpc correlates directly with their metal abundance, which may be related to their physical properties such as gas fractions. These results are consistent with the existence of an inner/outer halo structure for the halo system, as advocated by Carollo et al. (2010).

Photons propagating in strong magnetic fields are subject to a phenomenon called the "vacuum birefringence" where refractive indices of two physical modes both deviate from unity and are different from each other. We compute the vacuum polarization tensor of a photon in a static and homogeneous magnetic field by utilizing Schwinger's proper-time method, and obtain a series representation as a result of double integrals analytically performed with respect to proper-time variables. The outcome is expressed in terms of an infinite sum of known functions which is plausibly interpreted as summation over all the Landau levels of fermions. Each contribution from infinitely many Landau levels yields a kinematical condition above which the contribution has an imaginary part. This indicates decay of a sufficiently energetic photon into a fermion-antifermion pair with corresponding Landau level indices. Since we do not resort to any approximation, our result is applicable to the calculation of refractive indices in the whole kinematical region of a photon momentum and in any magnitude of the external magnetic field.

The Polarbear Cosmic Microwave Background (CMB) polarization experiment is currently observing from the Atacama Desert in Northern Chile. It will characterize the expected B-mode polarization due to gravitational lensing of the CMB, and search for the possible B-mode signature of inflationary gravitational waves. Its 250 mK focal plane detector array consists of 1,274 polarization-sensitive antenna-coupled bolometers, each with an associated lithographed band-defining filter. Each detector's planar antenna structure is coupled to the telescope's optical system through a contacting dielectric lenslet, an architecture unique in current CMB experiments. We present the initial characterization of this focal plane.

We present the design and characterization of the POLARBEAR experiment. POLARBEAR will measure the polarization of the cosmic microwave background (CMB) on angular scales ranging from the experiment's 3.5 arcminute beam size to several degrees. The experiment utilizes a unique focal plane of 1,274 antenna-coupled, polarization sensitive TES bolometers cooled to 250 milliKelvin. Employing this focal plane along with stringent control over systematic errors, POLARBEAR has the sensitivity to detect the expected small scale B-mode signal due to gravitational lensing and search for the large scale B-mode signal from inflationary gravitational waves. POLARBEAR was assembled for an engineering run in the Inyo Mountains of California in 2010 and was deployed in late 2011 to the Atacama Desert in Chile. An overview of the instrument is presented along with characterization results from observations in Chile.

We present theoretical calculations for the differential distribution of stellar orbital eccentricity for a sample of solar-neighbour halo stars. Two types of static, spherical gravitational potentials are adopted to define the eccentricity e for given energy E and angular momentum L, such as an isochrone potential and a Navarro-Frenk-White potential that can serve as two extreme ends covering in-between any realistic potential of the Milky Way halo. The solar-neighbour eccentricity distribution \Delta N(e) is then formulated, based on a static distribution function of the form f(E,L) in which the velocity anisotropy parameter \beta monotonically increases in the radial direction away from the galaxy center, such that beta is below unity (near isotropic velocity dispersion) in the central region and asymptotically approaches \sim 1 (radially anisotropic velocity dispersion) in the far distant region of the halo. We find that \Delta N(e) sensitively depends upon the radial profile of \beta, and this sensitivity is used to constrain such profile in comparison with some observational properties of \Delta N_{obs}(e) recently reported by Carollo et al. (2010). Especially, the linear e-distribution and the fraction of higher-e stars for their sample of solar-neighbour inner-halo stars rule out a constant profile of \beta, contrary to the opposite claim by Bond et al. (2010). Our constraint of \beta \lesssim 0.5 at the galaxy center indicates that the violent relaxation that has acted on the inner halo is effective within a scale radius of \sim 10 kpc from the galaxy center. We discuss that our result would help understand the formation and evolution of the Milky Way halo.

We observed diffuse cosmic and atmospheric gamma rays at balloon altitudes with the Sub-MeV gamma-ray Imaging Loaded-on-balloon Experiment I (SMILE-I) as the first step toward a future all-sky survey with a high sensitivity. SMILE-I employed an electron-tracking Compton camera comprised of a gaseous electron tracker as a Compton-scattering target and a scintillation camera as an absorber. The balloon carrying the SMILE-I detector was launched from the Sanriku Balloon Center of the Institute of Space and Astronomical Science/Japan Space Exploration Agency on September 1, 2006, and the flight lasted for 6.8 hr, including level flight for 4.1 hr at an altitude of 32-35 km. During the level flight, we successfully detected 420 downward gamma rays between 100 keV and 1 MeV at zenith angles below 60 degrees. To obtain the flux of diffuse cosmic gamma rays, we first simulated their scattering in the atmosphere using Geant4, and for gamma rays detected at an atmospheric depth of 7.0 g cm-2, we found that 50% and 21% of the gamma rays at energies of 150 keV and 1 MeV, respectively, were scattered in the atmosphere prior to reaching the detector. Moreover, by using Geant4 simulations and the QinetiQ atmospheric radiation model, we estimated that the detected events consisted of diffuse cosmic and atmospheric gamma rays (79%), secondary photons produced in the instrument through the interaction between cosmic rays and materials surrounding the detector (19%), and other particles (2%). The obtained growth curve was comparable to Ling's model, and the fluxes of diffuse cosmic and atmospheric gamma rays were consistent with the results of previous experiments. The expected detection sensitivity of a future SMILE experiment measuring gamma rays between 150 keV and 20 MeV was estimated from our SMILE-I results and was found to be ten times better than that of other experiments at around 1 MeV.

We present theoretical calculations for the differential distribution of stellar orbital eccentricity in a galaxy halo, assuming that the stars constitute a spherical, collisionless system in dynamical equilibrium with a dark matter halo. In order to define the eccentricity e of a halo star for given energy E and angular momentum L, we adopt two types of gravitational potential, such as an isochrone potential and a Navarro-Frenk-White potential, that could form two ends covering in-between any realistic potential of dark matter halo. Based on a distribution function of the form f(E,L) that allows constant anisotropy in velocity dispersions characterized by a parameter \beta, we find that the eccentricity distribution is a monotonically increasing function of e for the case of highly radially anisotropic velocity dispersions (\beta > 0.6), while showing a hump-like shape for the cases from radial through tangential velocity anisotropy (\beta < 0.6). We also find that when the velocity anisotropy agrees with that observed for the Milky Way halo stars (\beta = 0.5-0.7), a nearly linear eccentricity distribution of N(e) \alpha e results at e < 0.7, largely independent of the potential adopted. Our theoretical eccentricity distribution would be a vital tool of examining how far out in the halo the dynamical equilibrium has been achieved, through comparison with kinematics of halo stars sampled at greater distances. Given that large surveys of the SEGUE and Gaia projects would be in progress, we discuss how our results would serve as a new guide in exploring the formation and evolution of the Milky Way halo.

A direction-sensitive dark matter search experiment at Kamioka underground laboratory with the NEWAGE-0.3a detector was performed. The NEWAGE- 0.3a detector is a gaseous micro-time-projection chamber filled with CF4 gas at 152 Torr. The fiducial volume and target mass are 20*25*31 cm3 and 0.0115 kg, respectively. With an exposure of 0.524 kgdays, improved spin-dependent weakly interacting massive particle (WIMP)-proton cross section limits by a direction-sensitive method were achieved including a new record of 5400 pb for 150 GeV/c2 WIMPs. We studied the remaining background and found that ambient gamma-rays contributed about one-fifth of the remaining background and radioactive contaminants inside the gas chamber contributed the rest.

We present the case for a dark matter detector with directional sensitivity. This document was developed at the 2009 CYGNUS workshop on directional dark matter detection, and contains contributions from theorists and experimental groups in the field. We describe the need for a dark matter detector with directional sensitivity; each directional dark matter experiment presents their project's status; and we close with a feasibility study for scaling up to a one ton directional detector, which would cost around $150M.

Different numbers of self-gravitating particles (in different types of periodic motion) are most likely to generate very different shapes of gravitational waves, some of which, however, can be accidentally almost the same. One such example is a binary and a three-body system for Lagrange's solution. To track the evolution of these similar waveforms, we define a chirp mass to the triple system. Thereby, we show that the quadrupole waveforms cannot distinguish the sources. It is suggested that waveforms with higher $\ell$-th multipoles will be important for classification of them (with a conjecture of $\ell \leq N$ for N particles).

We developed a three-dimensional gaseous tracking device and performed a direction-sensitive dark matter search in a surface laboratory. By using 150 Torr carbon-tetrafluoride (CF_4 gas), we obtained a sky map drawn with the recoil directions of the carbon and fluorine nuclei, and set the first limit on the spin-dependent WIMP (Weakly Interacting Massive Particles)-proton cross section by a direction-sensitive method. Thus, we showed that a WIMP-search experiment with a gaseous tracking device can actually set limits. Furthermore, we demonstrated that this method will potentially play a certain role in revealing the nature of dark matter when a low-background large-volume detector is developed.

We have developed a position-sensitive scintillation camera with a large area absorber for use as an advanced Compton gamma-ray camera. At first we tested GSO(Ce) crystals. We compared light output from the GSO(Ce) crystals under various conditions: the method of surface polishing, the concentration of Ce, and co-doping Zr. As a result, we chose the GSO(Ce) crystals doped with only 0.5 mol% Ce, and its surface polished by chemical etching as the scintillator of our camera. We also made a 16$\times$16 cm$^2$ scintillation camera which consisted of 9 position-sensitive PMTs (PSPMTs Hamamatsu flat-panel H8500), the each of which had 8$\times$8 anodes with a pitch of 6 mm and coupled to 8$\times$8 arrays of pixelated 6$\times6\times$13 mm$^3$ GSO(Ce) scintillators. For the readout system of the 576 anodes of the PMTs, we used chained resistors to reduce the number of readout channels down to 48 to reduce power consumption. The camera has a position resolution of less than 6mm and a typical energy resolution of 10.5% (FWHM) at 662 keV at each pixel in a large area of 16$\times$16 cm$^2$. %to choose the best scintillator for our project. Furthermore we constructed a 16$\times$16 array of 3$\times3\times$13 mm$^3$ pixelated GSO(Ce) scintillators, and glued it to a PMT H8500. This camera had the position resolution of less than 3mm, over an area of 5$\times$5 cm$^2$, except for some of the edge pixels; the energy resolution was typically 13% (FWHM) at 662 keV.

The Time Projection Chamber (TPC) has been recognized as a potentially powerful detector for the search of WIMPs by measuring the directions of nuclear recoils, in which the most convincing signature of WIMPs, caused by the Earth's motion around the Galaxy, appears. We report on the first results of a performance study of the neutron exposure of our prototype micro-TPC with Ar-C$_2$H$_6$ (90:10) and CF$_4$ gas of 150 Torr.

A prototype of the MeV gamma-ray imaging camera based on the full reconstruction of the Compton process has been developed. This camera consists of a micro-TPC that is a gaseous Time Projection Chamber (TPC) and scintillation cameras. With the information of the recoil electrons and the scattered gamma-rays, this camera detects the energy and incident direction of each incident gamma-ray. We developed a prototype of the MeV gamma-ray camera with a micro-TPC and a NaI(Tl) scintillator, and succeeded in reconstructing the gamma-rays from 0.3 MeV to 1.3 MeV. Measured angular resolutions of ARM (Angular Resolution Measure) and SPD (Scatter Plane Deviation) for 356 keV gamma-rays were $18^\circ$ and $35^\circ$, respectively.