There are expected to be physical relationships between the globular clusters (GCs) and stellar substructures in the Milky Way, not all of which have yet been found. We search for such substructures from a combined halo sample of SDSS blue horizontal-branch and SDSS+LAMOST RR Lyrae stars, cross-matched with astrometric information from $Gaia$ DR2. This is a sample of old stars which are also excellent tracers of structures, ideal for searching for ancient relics in the outer stellar halo. By applying the neural-network-based method StarGO to the full 4D dynamical space of our sample, we rediscover the Sagittarius Stream, and find the debris of the $Gaia$-Enceladus-Sausage (GES) and the Sequoia events in the outer halo, as well as their linkages with several GCs. Most importantly, we find a new, low-mass, debris stream associated with a pair of GCs (NGC 5024 and NGC 5053), which we dub LMS-1. This stream has a very polar orbit, and occupies a region between 10 to 20 kpc from the Galactic center. NGC 5024 (M53), the more-massive of the associated GC pair, is very likely the nuclear star cluster of a now-disrupted dwarf galaxy progenitor, based on the results from N-body simulations.

Providing a possible connection between neutrino emission and gravitational-wave (GW) bursts is important to our understanding of the physical processes that occur when black holes or neutron stars merge. In the Daya Bay experiment, using data collected from December 2011 to August 2017, a search has been performed for electron-antineutrino signals coinciding with detected GW events, including GW150914, GW151012, GW151226, GW170104, GW170608, GW170814, and GW170817. We used three time windows of $\mathrm{\pm 10~s}$, $\mathrm{\pm 500~s}$, and $\mathrm{\pm 1000~s}$ relative to the occurrence of the GW events, and a neutrino energy range of 1.8 to 100 MeV to search for correlated neutrino candidates. The detected electron-antineutrino candidates are consistent with the expected background rates for all the three time windows. Assuming monochromatic spectra, we found upper limits (90% confidence level) on electron-antineutrino fluence of $(1.13~-~2.44) \times 10^{11}~\rm{cm^{-2}}$ at 5 MeV to $8.0 \times 10^{7}~\rm{cm^{-2}}$ at 100 MeV for the three time windows. Under the assumption of a Fermi-Dirac spectrum, the upper limits were found to be $(5.4~-~7.0)\times 10^{9}~\rm{cm^{-2}}$ for the three time windows.

We revisit the diagonal ridge feature (diagonal distributions in the $R, v_{\phi}$ plane) found in $Gaia$ and present timing analysis for it between Galactocentric distances of $R=7.5$ and 12 \,kpc, using Main-Sequence-Turn-Off and OB stars selected from the LAMOST Galactic spectroscopic surveys. We recover the ridge pattern in the $R$-$v_{\phi}$ plane color coded by stellar number density and mean radial velocity and find this feature is presented from very young (OB stars, few hundred \,Myr) to very old populations ($\tau > $ 9 \,Gyr), accompanying with the north-south asymmetrical ridges. One of the ridge patterns, with constant angular momentum per unit mass, shows clear variations with different age populations compared. However, the remaining two ones are very stable, clearly implying there might have two kinds of ridge patterns with different dynamical origins and evolution. For the vertical velocity and height, there are no ridge features corresponding to the in-plane Galactoseismogy features, but the metallicity and [$\alpha/Fe$] have some weak signals. This first temporal chemo-kinematical analysis for ridge strongly supports that the ridge pattern might be originated from the phase mixing of disrupting spiral arms implying different ridges have different physical processes and no coupling phenomenon, although other possible scenarios such as Sgr perturbations could not be ruled out.

By combining LAMOST DR4 and Gaia DR2 common red clump stars with age and proper motion, we analyze the amplitude evolution of the stellar warp independently of any assumption with a simple model. The greatest height of the warp disk increases with Galactocentric distance in different populations and it is dependent on the age: the younger stellar populations exhibit stronger warp features than the old ones, accompanied with the warp amplitude $\gamma ({\rm age})$ decreasing with age and its first derivative $\dot{\gamma} ({\rm age})$ is different from zero. The azimuth of line of nodes $\phi _w$ is stable at $-$5 degree without clear time evolution, which perfectly confirms some of previous works. All these self-consistent evidences support that our Galactic warp should most likely be a long lived, but non-steady structure and not a transient one, which is supporting the warp is originated from gas infall onto the disk or other hypotheses that suppose that the warp mainly affects to the gas, and consequently younger populations tracing the gas are stronger than the older ones. In other words, the Galactic warp is induced by the non-gravitational interaction over the disk models.

A large number of new members ($\sim$150) of the Cetus Stream (CS) were identified from their clustering features in dynamical space using 6D kinematic data by combining LAMOST DR5 and Gaia DR2 surveys. They map a diffuse structure that extends over at least 100 degrees in the northern and southern Galactic hemispheres, at heliocentric distances between 20 to 50 kpc. Taking advantage of this expanded dataset, we model the stream with a suite of tailored N-body simulations. Our findings exclude the possibility that the NGC 5824 globular cluster is the core of the progenitor of the stream, as postulated by previous studies. Our best models, which successfully reproduce the features of the CS indicate that the progenitor is likely a dwarf galaxy of $\sim$ 2$\times$10$^9$M$_{\odot}$, with a diffuse disc morphology. The merger occured $\sim$ 5 Gyr ago and since then it has experienced approximately eight apo-center passages. Our results suggest that NGC 5824 was either a globular cluster situated off-centre in the dwarf progenitor or, alternatively, it was the nuclear star cluster of another dwarf galaxy that has very similar orbit as the progenitor of the CS. In both scenarios, the progenitor systems would leave streams around NGC 5824, but with distinct distance distributions. To discriminate between these scenarios, the detection and accurate distance measurements of the predicted stream around the GC are crucial, which will be possible in the upcoming LSST era. Our simulations also predict that part of the Southern Cetus stream is very likely the newly discovered Palca stream, and possibly related to another, more diffuse Southern substructure, the Eridanus-Pheonix overdensity.

Context: In our Paper I, by using statistical deconvolution methods, extended kinematics maps of Gaia-DR2 data have been produced in a range of heliocentric distances that are a factor of two to three larger than those analyzed previously by the Gaia Collaboration with the same data. It added the range of Galactocentric distances between 13 kpc and 20 kpc to the previous maps. Aims: Here, we investigate the dynamical effects produced by different mechanisms that can explain the radial and vertical components of these extended kinematic maps, including a decomposition of bending and breathing of the vertical components. This paper as a whole tries to be a compendium of different dynamical mechanisms whose predictions can be compared to the kinematic maps. Methods: Using analytical methods or simulations, we are able to predict the main dynamical factors and compare them to the predictions of the extended kinematic maps of Gaia-DR2. Results: The gravitational influence of Galactic components that are different from the disk, such as the long bar or bulge, the spiral arms, or a tidal interaction with Sagittarius dwarf galaxy, may explain some features of the velocity maps, especially in the inner parts of the disk. However, they are not sufficient in explaining the most conspicuous gradients in the outer disk. Vertical motions might be dominated by external perturbations or mergers, although a minor component may be due to a warp whose amplitude evolves with time. Here, we show with two different methods, which analyze the dispersion of velocities, that the mass distribution of the disk is flared. Despite these partial explanations, the main observed features can only be explained in terms of out-of-equilibrium models, which are either due to external perturbers or to the fact that the disk has not had time to reach equilibrium since its formation.

Dark matter produced from thermal freeze-out is typically restricted to have masses above roughly 1 MeV. However, if the couplings are small, the freeze-in mechanism allows for production of dark matter down to keV masses. We consider dark matter coupled to a dark photon that mixes with the photon and dark matter coupled to photons through an electric or magnetic dipole moment. We discuss contributions to the freeze-in production of such dark matter particles from standard model fermion-antifermion annihilation and plasmon decay. We also derive constraints on such dark matter from the cooling of red giant stars, horizontal branch stars, and the Sun, carefully evaluating the thermal processes as well as the Compton scattering that dominates for masses above the plasma frequency. For the dark photon portal dark matter, the parameters to obtain the observed relic abundance from freeze-in are excluded below a few tens of keV, depending on the value of the dark gauge coupling constant. For dark matter with an electric or magnetic dipole moment, the freeze-in production parameters are barely constrained through stellar cooling arguments. While laboratory probes are unlikely to probe these freeze-in scenarios in general, we show that for dark matter with an electric or magnetic dipole moment and for dark matter masses above the reheating temperature, the couplings needed for freeze-in to produce the observed relic abundance can be probed partially by upcoming direct-detection experiments.

We perform analysis of the three-dimensional kinematics of Milky Way disk stars in mono-age populations. We focus on stars between Galactocentric distances of $R=6$ and 14 \,kpc, selected from the combined LAMOST DR4 red clump giant stars and Gaia DR2 proper motion catalogue. We confirm the 3D asymmetrical motions of recent works and provide time tagging of the Galactic outer disk asymmetrical motions near the anticenter direction out to Galactocentric distances of 14\,kpc. Radial Galactocentric motions reach values up to 10 km s$^{-1}$, depending on the age of the population, and present a north-south asymmetry in the region corresponding to density and velocity substructures that were sensitive to the perturbations in the early 6 \,Gyr. After that time, the disk stars in this asymmetrical structure have become kinematically hotter, and are thus not sensitive to perturbations, and we find the structure is a relatively younger population. With quantitative analysis, we find stars both above and below the plane at $R\gtrsim 9$ kpc that exhibit bending mode motions of which the sensitive duration is around 8 \,Gyr. We speculate that the in-plane asymmetries might not be mainly caused by a fast rotating bar, intrinsically elliptical outer disk, secular expansion of the disk, or streams. Spiral arm dynamics, out-of-equilibrium models, minor mergers or others are important contributors. Vertical motions might be dominated by bending and breathing modes induced by complicated inner or external perturbers. It is likely that many of these mechanisms are coupled together.

The precise measurement of the spectrum of protons, the most abundant component of the cosmic radiation, is necessary to understand the source and acceleration of cosmic rays in the Milky Way. This work reports the measurement of the cosmic ray proton fluxes with kinetic energies from 40 GeV to 100 TeV, with two and a half years of data recorded by the DArk Matter Particle Explorer (DAMPE). This is the first time an experiment directly measures the cosmic ray protons up to ~100 TeV with a high statistics. The measured spectrum confirms the spectral hardening found by previous experiments and reveals a softening at ~13.6 TeV, with the spectral index changing from ~2.60 to ~2.85. Our result suggests the existence of a new spectral feature of cosmic rays at energies lower than the so-called knee, and sheds new light on the origin of Galactic cosmic rays.

We characterize the kinematic and chemical properties of $\sim$3,000 Sagittarius (Sgr) stream stars, including K-giants, M-giants, and BHBs, select from SEGUE-2, LAMOST, and SDSS separately in Integrals-of-Motion space. The orbit of Sgr stream is quite clear from the velocity vector in $X$-$Z$ plane. Stars traced by K-giants and M-giants present the apogalacticon of trailing steam is $\sim$ 100 kpc. The metallicity distributions of Sgr K-, M-giants, and BHBs present that the M-giants are on average the most metal-rich population, followed by K-giants and BHBs. All of the K-, M-giants, and BHBs indicate that the trailing arm is on average more metal-rich than leading arm, and the K-giants show that the Sgr debris is the most metal-poor part. The $\alpha$-abundance of Sgr stars exhibits a similar trend with the Galactic halo stars at lower metallicity ([Fe/H] $<\sim$ $-$1.0 dex), and then evolve down to lower [$\alpha$/Fe] than disk stars at higher metallicity, which is close to the evolution pattern of $\alpha$-element of Milky Way dwarf galaxies. We find $V_Y$ and metallicity of K-giants have gradients along the direction of line-of-sight from the Galactic center in $X$-$Z$ plane, and the K-giants show that $V_Y$ increases with metallicity at [Fe/H] $>\sim-$1.5 dex. After dividing the Sgr stream into bright and faint stream according to their locations in equatorial coordinate, the K-giants and BHBs show that the bright and faint stream present different $V_Y$ and metallicities, the bright stream is on average higher in $V_Y$ and metallicity than the faint stream.

We investigate the three-dimensional asymmetrical kinematics and present time stamps of the Milky Way disk between Galactocentric distances of $R=12$ and 15 \,kpc, using red clump stars selected from the LAMOST Galactic survey, also with proper motion measurements provided by the Gaia DR2. We discover velocity substructure above the Galactic plane corresponding to a density dip found recently ("south-middle opposite" density structure[R $\sim$ 12-15 \,kpc, Z $\sim$ 1.5 \,kpc] discovered in \citet{Wang2018b, Wang2018c}) in the radial and azimuthal velocity. For the vertical velocity, we detect clear vertical bulk motions or bending mode motions, which has no clear north-south asymmetry corresponding to the in-plane asymmetrical features. In the subsample of stars with different ages, we find that there is little temporal evolution of the in-plane asymmetry from 0$-$14 \,Gyr, which means the structure is sensitive to the perturbations in almost cosmic time possibly. We propose that the possible scenario of this asymmetric velocity structure is caused by the mechanisms generated in-plane, rather than vertical perturbations.

The DArk Matter Particle Explorer (DAMPE), a satellite-based cosmic ray and gamma-ray detector, was launched on December 17, 2015, and began its on-orbit operation on December 24, 2015. In this work we document the on-orbit calibration procedures used by DAMPE and report the calibration results of the Plastic Scintillator strip Detector (PSD), the Silicon-Tungsten tracKer-converter (STK), the BGO imaging calorimeter (BGO), and the Neutron Detector (NUD). The results are obtained using Galactic cosmic rays, bright known GeV gamma-ray sources, and charge injection into the front-end electronics of each sub-detector. The determination of the boundary of the South Atlantic Anomaly (SAA), the measurement of the live time, and the alignments of the detectors are also introduced. The calibration results demonstrate the stability of the detectors in almost two years of the on-orbit operation.

We use a novel cluster identification tool StarGO to explore the metal poor ([Fe/H] $<$ -1.5) outer stellar halo (d $>$ 15 kpc) of the Milky Way using data from Gaia, LAMOST and SDSS. Our method is built using an unsupervised learning algorithm, a self-organizing map, which trains a 2-D neural network to learn the topological structures of a data set from an n-D input space. Using a 4-D space of angular momentum and orbital energy, we identify three distinct groups corresponding to the Sagittarius, Orphan, and Cetus Streams. For the first time we are able to discover a northern counterpart to the Cetus stream. We test the robustness of this new detection using mock data and find that the significance is more than 5-sigma. We also find that the existing southern counterpart bifurcates into two clumps with different radial velocities. By exploiting the visualization power of StarGO, we attach MW globular clusters to the same trained neural network. The Sagittarius stream is found to have five related clusters, confirming recent literature studies, and the Cetus stream has one associated cluster, NGC 5824. This latter association has previously been postulated, but can only now be truly confirmed thanks to the high-precision Gaia proper motions and large numbers of stellar spectra from LAMOST. The large metallicity dispersion of the stream indicates that the progenitor cannot be a globular cluster. Given the mean metallicity of the stream, we propose that the stream is the result of a merger of a low-mass dwarf galaxy that hosted a large nuclear star cluster (NGC 5824).

We construct a large halo K-giant sample by combining the positions, distances, radial velocities, and metallicities of over 13,000 LAMOST DR5 halo K giants with the Gaia DR2 proper motions, which covers a Galactocentric distance range of 5-120 kpc. Using a position-velocity clustering estimator (the 6Distance), we statistically quantify the presence of position-velocity substructure at high significance: K giants have more close pairs in position-velocity space than a smooth stellar halo. We find that the amount of substructure in the halo increases with increasing distance and metallicity. With a percolation algorithm named friends-of-friends (FoF) to identify groups, we identify members belonging to Sagittarius (Sgr) Streams, Monoceros Ring, Virgo overdensity, Hercules-Aquila Cloud, Orphan Streams and other unknown substructures and find that the Sgr streams account for a large part of grouped stars beyond 20 kpc and enhance the increase of substructure with distance and metallicity. For the first time, we identify spectroscopic members of Monoceros Ring in the south and north Galactic hemisphere, which presents a rotation of about 185 km s^{-1} and mean metallicity is -0.66 dex.

In order to search for the line-like signals in the Fermi-LAT data, we have analyzed totally 49152 regions of interest (ROIs) that cover the whole sky. No ROI displays a line signal with test statistic (TS) value above 25, while for 50 ROIs weak line-like excesses with ${\rm TS}>16$ are presented. The intrinsic significances of these potential signals are further reduced by the large trial factor introduced in such kind of analysis. For the largest TS value of 24.3 derived in our analysis, the corresponding global significance is only $0.54\sigma$. We thus do not find any significant line-like signal and set up constraints on the cross section of dark matter annihilating to gamma-ray lines, $\left<\sigma{v}\right>_{\gamma\gamma}$.

We revisit constraints on dark photons with masses below ~ 100 MeV from the observations of Supernova 1987A. If dark photons are produced in sufficient quantity, they reduce the amount of energy emitted in the form of neutrinos, in conflict with observations. For the first time, we include the effects of finite temperature and density on the kinetic-mixing parameter, epsilon, in this environment. This causes the constraints on epsilon to weaken with the dark-photon mass below ~ 15 MeV. For large-enough values of epsilon, it is well known that dark photons can be reabsorbed within the supernova. Since the rates of reabsorption processes decrease as the dark-photon energy increases, we point out that dark photons with energies above the Wien peak can escape without scattering, contributing more to energy loss than is possible assuming a blackbody spectrum. Furthermore, we estimate the systematic uncertainties on the cooling bounds by deriving constraints assuming one analytic and four different simulated temperature and density profiles of the proto-neutron star. Finally, we estimate also the systematic uncertainty on the bound by varying the distance across which dark photons must propagate from their point of production to be able to affect the star. This work clarifies the bounds from SN1987A on the dark-photon parameter space.

The DArk Matter Particle Explorer (DAMPE) developed in China was designed to search for evidence of dark matter particles by observing primary cosmic rays and gamma rays in the energy range from 5 GeV to 10 TeV. Since its launch in December 2015, a large quantity of data has been recorded. With the data set acquired during more than a year of operation in space, a precise time-dependent calibration for the energy measured by the BGO ECAL has been developed. In this report, the instrumentation and development of the BGO Electromagnetic Calorimeter (BGO ECAL) are briefly described. The calibration on orbit, including that of the pedestal, attenuation length, minimum ionizing particle peak, and dynode ratio, is discussed, and additional details about the calibration methods and performance in space are presented.

We present the complete history of structure formation in a simple dissipative dark-sector model. The model has only two particles: a dark electron, which is a subdominant component of dark matter, and a dark photon. Dark-electron perturbations grow from primordial overdensities, become non-linear, and form dense dark galaxies. Bremsstrahlung cooling leads to fragmentation of the dark-electron halos into clumps that vary in size from a few to millions of solar masses, depending on the particle model parameters. In particular, we show that asymmetric dark stars and black holes form within the Milky Way from the collapse of dark electrons. These exotic compact objects may be detected and their properties measured at new high-precision astronomical observatories, giving insight into the particle nature of the dark sector without the requirement of non-gravitational interactions with the visible sector.

We consider the constraints from Supernova 1987A on particles with small couplings to the Standard Model. We discuss a model with a fermion coupled to a dark photon, with various mass relations in the dark sector; millicharged particles; dark-sector fermions with inelastic transitions; the hadronic QCD axion; and an axion-like particle that couples to Standard Model fermions with couplings proportional to their mass. In the fermion cases, we develop a new diagnostic for assessing when such a particle is trapped at large mixing angles. Our bounds for a fermion coupled to a dark photon constrain small couplings and masses <200 MeV, and do not decouple for low fermion masses. They exclude parameter space that is otherwise unconstrained by existing accelerator-based and direct-detection searches. In addition, our bounds are complementary to proposed laboratory searches for sub-GeV dark matter, and do not constrain several "thermal" benchmark-model targets. For a millicharged particle, we exclude charges between 10^(-9) to a few times 10^(-6) in units of the electron charge; this excludes parameter space to higher millicharges and masses than previous bounds. For the QCD axion and an axion-like particle, we apply several updated nuclear physics calculations and include the energy dependence of the optical depth to accurately account for energy loss at large couplings. We rule out a hadronic axion of mass between 0.1 and a few hundred eV, or equivalently bound the PQ scale between a few times 10^4 and 10^8 GeV, closing the hadronic axion window. For an axion-like particle, our bounds disfavor decay constants between a few times 10^5 GeV up to a few times 10^8 GeV. In all cases, our bounds differ from previous work by more than an order of magnitude across the entire parameter space. We also provide estimated systematic errors due to the uncertainties of the progenitor.

In this paper we study sneutrino dark matter in a recently proposed supersymmetric electroweak-scale inverse seesaw model, in which the majority of the sneutrino dark matter particle is a mixture of the right-handed sneutrino $\tilde{N}^c$ and the singlet field $\tilde{S}$. The scalar field $X$ responsible for the generation of neutrino masses can simultaneously play a crucial role for sneutrino annihilation in the early Universe via the pseudoscalar mediator $A_X$ into neutrinos. We focus here on the dominant annihilation channels and provide all the formulas together with analytic estimates in order to identify the relevant parameters. Furthermore, we show that the direct detection scattering cross section is many orders of magnitude below the current limits, and estimate the indirect detection annihilation rate, which is only a few orders of magnitude below the current limits.

We develop a new method StarGO (Stars' Galactic Origin) to identify the galactic origins of halo stars using their kinematics. Our method is based on self-organizing map (SOM), which is one of the most popular unsupervised learning algorithms. StarGO combines SOM with a novel adaptive group identification algorithm with essentially no free parameters. In order to evaluate our model, we build a synthetic stellar halo from mergers of nine satellites in the Milky Way. We construct the mock catalogue by extracting a heliocentric volume of 10 kpc from our simulations and assigning expected observational uncertainties corresponding to bright stars from Gaia DR2 and LAMOST DR5. We compare the results from StarGO against that from a Friends-of-Friends (FoF) based method in the space of orbital energy and angular momentum. We show that StarGO is able to systematically identify more satellites and achieve higher number fraction of identified stars for most of the satellites within the extracted volume. When applied to data from Gaia DR2, StarGO will enable us to reveal the origins of the inner stellar halo in unprecedented detail.

In this work, we search for $\gamma$-ray emission from the directions of some nearby Milky Way dwarf spheroidal galaxies (dSphs) and candidates with the publicly-available Pass 8 data of Fermi-LAT. Our sample includes 12 sources with the distances $<50$ kpc. Very weak $\gamma$-ray excesses ($\sim 2\sigma$) are found in some dSphs/candidates, consistent with that reported in the previous literature. Intriguingly, the peak test statistic (TS) value of the weak emission from the direction of Reticulum II rises continually. If interpreted as dark matter (DM) annihilation, the peak TS value is 13.5 for the annihilation channel of $\chi\chi \rightarrow \tau^{+}\tau^{-}$ and the DM mass of $m_\chi \sim 16$ GeV. The combination of all these nearby sources yields a more significant (with local significance $> 4\sigma$) $\gamma$-ray signal.

High energy cosmic ray electrons plus positrons (CREs), which lose energy quickly during their propagation, provide an ideal probe of Galactic high-energy processes and may enable the observation of phenomena such as dark-matter particle annihilation or decay. The CRE spectrum has been directly measured up to $\sim 2$ TeV in previous balloon- or space-borne experiments, and indirectly up to $\sim 5$ TeV by ground-based Cherenkov $\gamma$-ray telescope arrays. Evidence for a spectral break in the TeV energy range has been provided by indirect measurements of H.E.S.S., although the results were qualified by sizeable systematic uncertainties. Here we report a direct measurement of CREs in the energy range $25~{\rm GeV}-4.6~{\rm TeV}$ by the DArk Matter Particle Explorer (DAMPE) with unprecedentedly high energy resolution and low background. The majority of the spectrum can be properly fitted by a smoothly broken power-law model rather than a single power-law model. The direct detection of a spectral break at $E \sim0.9$ TeV confirms the evidence found by H.E.S.S., clarifies the behavior of the CRE spectrum at energies above 1 TeV and sheds light on the physical origin of the sub-TeV CREs.

The DArk Matter Particle Explorer (DAMPE), one of the four scientific space science missions within the framework of the Strategic Pioneer Program on Space Science of the Chinese Academy of Sciences, is a general purpose high energy cosmic-ray and gamma-ray observatory, which was successfully launched on December 17th, 2015 from the Jiuquan Satellite Launch Center. The DAMPE scientific objectives include the study of galactic cosmic rays up to $\sim 10$ TeV and hundreds of TeV for electrons/gammas and nuclei respectively, and the search for dark matter signatures in their spectra. In this paper we illustrate the layout of the DAMPE instrument, and discuss the results of beam tests and calibrations performed on ground. Finally we present the expected performance in space and give an overview of the mission key scientific goals.

We present new evidence for eight early-type galaxies (ETGs) from the CALIFA Survey that show clear rotation around their major photometric axis ("prolate rotation"). These are LSBCF560-04, NGC 0647, NGC 0810, NGC 2484, NGC 4874, NGC 5216, NGC 6173 and NGC 6338. Including NGC 5485, a known case of an ETG with stellar prolate rotation, as well as UGC 10695, a further possible candidate for prolate rotation, we report ten CALIFA galaxies in total that show evidence for such a feature in their stellar kinematics. Prolate rotators correspond to ~9% of the volume-corrected sample of CALIFA ETGs, a fraction much higher than previously reported. We find that prolate rotation is more common among the most massive ETGs. We investigate the implications of these findings by studying N-body merger simulations, and show that a prolate ETG with rotation around its major axis could be the result of a major polar merger, with the amplitude of prolate rotation depending on the initial bulge-to-total stellar mass ratio of its progenitor galaxies. Additionally, we find that prolate ETGs resulting from this formation scenario show a correlation between their stellar line-of-sight velocity and higher order moment h_3, opposite to typical oblate ETGs, as well as a double peak of their stellar velocity dispersion along their minor axis. Finally, we investigate the origin of prolate rotation in polar galaxy merger remnants. Our findings suggest that prolate rotation in massive ETGs might be more common than previously expected, and can help towards a better understanding of their dynamical structure and formation origin.

This white paper summarizes the workshop "U.S. Cosmic Visions: New Ideas in Dark Matter" held at University of Maryland on March 23-25, 2017.

A software system has been developed for the DArk Matter Particle Explorer (DAMPE) mission, a satellite-based experiment. The DAMPE software is mainly written in C++ and steered using Python script. This article presents an overview of the DAMPE offline software, including the major architecture design and specific implementation for simulation, calibration and reconstruction. The whole system has been successfully applied to DAMPE data analysis, based on which some results from simulation and beam test experiments are obtained and presented.

Milky Way-like galaxies are predicted to host a very large number of dark matter subhalos. Some massive and nearby subhalos could generate detectable gamma-rays, appearing as unidentified, spatially-extended and stable gamma-ray sources. We search for such sources in the third Fermi Large Area Telescope source List (3FGL) and report the identification of a new candidate, 3FGL J1924.8-1034. With the Fermi-LAT Pass 8 data, we find that 3FGL J1924.8-1034 is spatially-extended at a high confidence level of $5.4\sigma$, with a best-fit extension radius of $\sim0.15^{\circ}$. No significant variability has been found and its gamma-ray spectrum is well fitted by the dark matter annihilation into $b\bar{b}$ with a mass of $\sim 43$ GeV. All these facts make 3FGL J1924.8-1034 a possible dark matter subhalo candidate. However, due to the limited angular resolution, the possibility of that the spatial extension of 3FGL J1924.8-1034 is caused by the contamination from the other un-resolved point source can not be ruled out.

N-body simulations predict that galaxies at the Milky Way scale host a large number of dark matter (DM) subhalos. Some of these subhalos, if they are massive enough or close enough to the Earth, might be detectable in $\gamma$ rays due to the DM annihilation. 3FGL J2212.5+0703, an unidentified gamma-ray source, has been suggested to be the counterpart candidate of a DM subhalo by Bertoni et al. (2015, 2016). In this work we analyze the Fermi-LAT Pass 8 data of 3FGL J2212.5+0703 to independently test the DM subhalo hypothesis of this source. In order to suppress the possible contamination from two nearby very-bright blazars, we just take into account the front-converting gamma-rays which have better angular resolutions than that of the back-converting photons. In addition to the spatial distribution analysis, we have extended the spectrum analysis down to the energies of $\sim 100$ MeV, and thoroughly examined the variability of the emission during the past 8 years. We confirm that 3FGL J2212.5+0703 is a steady and spatially-extended gamma-ray emitter at a high confidence level. The spectrum is well consistent with that expected from DM annihilation into $b\bar{b}$. The introduction of a phenomenological LogParabola spectrum just improves the fit slightly. All these results suggest that 3FGL J2212.5+0703 could be indicative of a DM subhalo.

With 91 months of the publicly available Fermi-LAT Pass 8 data, we analyze the gamma-ray emission from the Milky Way satellites to search for potential line signals due to the annihilation of dark matter particles into double photons. The searched targets include a sample of dwarf spheroidal galaxies, the Large Magellanic Cloud (LMC) and Small Magellanic Cloud (SMC). No significant line emission has been found neither in the stacked dwarf galaxy sample nor in the direction of LMC/SMC. The corresponding upper limits on the cross section of DM annihilation into two photons are derived. Compared with results of previous gamma-ray line searches with the Pass 7 data, the current constraints on the line emission from dwarf spheroidal galaxies has been significantly improved in a wide energy range. With the rapid increase of the sample of dwarf spheroidal galaxies (candidates), we expect that the sensitivity of gamma ray line searches will be significantly improved in the near future.

Galaxy clusters are the largest gravitationally bound objects in the universe and may be suitable targets for indirect dark matter searches. With 85 months of Fermi-LAT Pass 8 publicly available data, we analyze the gamma-ray emission in the directions of 16 nearby Galaxy Clusters with an unbinned likelihood analysis. No globally statistically-significant $\gamma-$ray line feature is identified and a tentative line signal may be present at $\sim 43$ GeV. The 95\% confidence level upper limits on the velocity-averaged cross section of dark matter particles annihilating into double $\gamma-$rays (i.e., $\langle \sigma v \rangle_{\chi\chi\rightarrow \gamma\gamma}$) are derived. Unless very optimistic boost factors of dark matter annihilation in these Galaxy Clusters have been assumed, such constraints are much weaker than the bounds set by the Galactic $\gamma-$ray data.

HERD is the High Energy cosmic-Radiation Detection instrument proposed to operate onboard China's space station in the 2020s. It is designed to detect energetic cosmic ray nuclei, leptons and photons with a high energy resolution ($\sim1\%$ for electrons and photons and $20\%$ for nuclei) and a large geometry factor ($>3\,{ m^2\,sr}$ for electrons and diffuse photons and $>2\,{ m^2\,sr}$ for nuclei). In this work we discuss the capability of HERD to detect monochromatic $\gamma$-ray lines, based on simulations of the detector performance. It is shown that HERD will be one of the most sensitive instruments for monochromatic $\gamma$-ray searches at energies between $\sim10$ to a few hundred GeV. Above hundreds of GeV, Cherenkov telescopes will be more sensitive due to their large effective area. As a specific example, we show that a good portion of the parameter space of a supersymmetric dark matter model can be probed with HERD.

Very recently the Dark Energy Survey (DES) Collaboration has released their second group of Dwarf spheroidal (dSph) galaxy candidates. With the publicly-available Pass 8 data of Fermi-LAT we search for $\gamma-$ray emissions from the directions of these eight newly discovered dSph galaxy candidates. No statistically significant $\gamma-$ray signal has been found in the combined analysis of these sources. With the empirically estimated J-factors of these sources, the constraint on the annihilation channel of $\chi\chi \rightarrow \tau^{+}\tau^{-}$ is comparable to the bound set by the joint analysis of fifteen previously known dSphs with kinematically constrained J-factors for the dark matter mass $m_\chi>250$ GeV. In the direction of Tucana III (DES J2356-5935), one of the nearest dSph galaxy candidates that is $\sim 25$ kpc away, there is a weak $\gamma-$ray signal and its peak test statistic (TS) value for the dark matter annihilation channel $\chi\chi\rightarrow \tau^{+}\tau^{-1}$ is $\approx 6.7$ at $m_\chi \sim 15$ GeV. The significance of the possible signal likely increases with time. More data is highly needed to pin down the physical origin of such a GeV excess.

With the accurate cosmic ray (CR) electron and positron spectra (denoted as $\Phi_{\rm e^{-}}$ and $\Phi_{\rm e^{+}}$, respectively) measured by AMS-02 collaboration, the difference between the electron and positron fluxes (i.e., $\Delta \Phi=\Phi_{\rm e^{-}}-\Phi_{\rm e^{+}}$), dominated by the propagated primary electrons, can be reliably inferred. In the standard model, the spectrum of propagated primary CR electrons at energies $\geq 30$ GeV softens with the increase of energy. The absence of any evidence for such a continuous spectral softening in $\Delta \Phi$ strongly suggests a significant `excess' of primary CR electrons and at energies of $100-400$ GeV the identified excess component has a flux comparable to that of the observed positron excess. Middle-age but `nearby' supernova remnants (e.g., Monogem and Geminga) are favored sources for such an excess.

Several groups have analyzed the publicly-available Fermi-LAT data and reported a spatially extended $\gamma-$ray excess of around $1-3$ GeV from the region surrounding the Galactic Center that might originate from annihilation of dark matter particles with a rest mass $m_\chi \sim 30-40$ GeV. In this work we examine the role of the diffuse Galactic gamma ray emission (DGE) templates played in suppressing the GeV excess. For such a purpose, we adopt in total 128 background templates that have been generated by Ackermann et al. \cite{FermiLAT:2012aa} in the study of the {Fermi-LAT} observations of the diffuse gamma ray emission considering the effects of cosmic rays and the interstellar medium. The possible GeV excess, assumed to follow the spatial distribution of the prompt gamma-rays produced in the annihilation of dark matter particles taking a generalized NFW profile with an inner slope $\alpha=1.2$, has been analyzed in some regions of interest. The introduction of such an additional component centered at the Galactic center is found to have improved the goodness of fit to the data significantly in all background template models regardless of whether the excess spectrum is fixed or not. Our results thus suggest that the presence of a statistically significant GeV excess in the inner Galaxy is robust thought its spectrum depends on the DGE model adopted in the analysis. The possible physical origin of the GeV excess component is discussed and in the dark matter model the annihilation cross section of such particles is evaluated.

The China Dark Matter Experiment reports results on light WIMP dark matter searches at the China Jinping Underground Laboratory with a germanium detector array with a total mass of 20 g. The physics threshold achieved is 177 eVee ("ee" represents electron equivalent energy) at 50% signal efficiency. With 0.784 kg-days of data, exclusion region on spin-independent coupling with the nucleon is derived, improving over our earlier bounds at WIMP mass less than 4.6 GeV.

We describe the instrument concept of a high angular resolution telescope dedicated to the sub-GeV (from $\gtrsim$10 MeV to $\gtrsim$1 GeV) gamma-ray photon detection. This mission, named PANGU (PAir-productioN Gamma-ray Unit), has been suggested as a candidate for the joint small mission between the European Space Agency (ESA) and the Chinese Academy of Science (CAS). A wide range of topics of both astronomy and fundamental physics can be attacked with PANGU, covering Galactic and extragalactic cosmic-ray physics, extreme physics of a variety of extended (e.g. supernova remnants, galaxies, galaxy clusters) and compact (e.g. black holes, pulsars, gamma-ray bursts) objects, solar and terrestrial gamma-ray phenomena, and searching for dark matter decay and/or annihilation signature etc. The unprecedented point spread function can be achieved with a pair-production telescope with a large number of thin active tracking layers to precisely reconstruct the pair-produced electron and positron tracks. Scintillating fibers or thin silicon micro-strip detectors are suitable technology for such a tracker. The energy measurement is achieved by measuring the momentum of the electrons and positrons through a magnetic field. The innovated spectrometer approach provides superior photon pointing resolution, and is particular suitable in the sub-GeV range. The level of tracking precision makes it possible to measure the polarization of gamma rays, which would open up a new frontier in gamma-ray astronomy. The frequent full-sky survey at sub-GeV with PANGU's large field of view and significantly improved point spread function would provide crucial information to GeV-TeV astrophysics for current/future missions including Fermi, DAMPE, HERD, and CTA, and other multi-wavelength telescopes.

We use high-resolution Eulerian simulations to study the stability of cold gas flows in a galaxy size dark matter halo (10^12 Msun) at redshift z=2. Our simulations show that a cold stream penetrating a hot gaseous halo is stable against thermal convection and Kelvin-Helmholtz instability. We then investigate the effect of a satellite orbiting the main halo in the plane of the stream. The satellite is able to perturb the stream and to inhibit cold gas accretion towards the center of the halo for 0.5 Gyr. However, if the supply of cold gas at large distances is kept constant, the cold stream is able to re-establish itself after 0.3 Gyr. We conclude that cold streams are very stable against a large variety of internal and external perturbations.

A physical model of Geant4-based simulation of galactic cosmic ray (GCR) particles interaction with the lunar surface matter has been developed to investigate the production rate of cosmogenic nuclei. In this model the GCRs, mainly very high energy protons and $\alpha$ particles, bombard the surface of the Moon and produce many secondary particles such as protons and neutrons. The energies of proton and neutron at different depths are recorded and saved into ROOT files, and the analytical expressions for the differential proton and neutron fluxes are obtained through the best-fit procedure under the ROOT software. To test the validity of this model, we calculate the production rates of long-lived nuclei $^{10}$Be and $^{26}$Al in the Apollo 15 long drill core by combining the above differential fluxes and the newly evaluated spallation reaction cross sections. Numerical results show that the theoretical production rates agree quite well with the measured data. It means that this model works well. Therefore, it can be expected that this model can be used to investigate the cosmogenic nuclei in lunar samples returned by Chinese lunar exploration program and can be extended to study other objects, such as the meteorites and the Earth's atmosphere.

The China Dark matter Experiment collaboration reports the first experimental limit on WIMP dark matter from 14.6 kg-day of data taken with a 994 g p-type point-contact germanium detector at the China Jinping underground Laboratory where the rock overburden is more than 2400 m. The energy threshold achieved was 400 eVee. According to the 14.6 kg-day live data, we placed the limit of N= 1.75 * 10^{-40} cm^{2} at 90% confidence level on the spin-independent cross-section at WIMP mass of 7 GeV before differentiating bulk signals from the surface backgrounds.

The data collected by ATIC, CREAM and PAMELA all display remarkable cosmic-ray-nuclei spectrum hardening above the magnetic rigidity $\sim$ 240 GV. One natural speculation is that the primary electron spectrum also gets hardened (possibly at $\sim 80$ GV) and the hardening partly accounts for the electron/positron total spectrum excess discovered by ATIC, HESS and Fermi-LAT. If it is the case, the increasing behavior of the subsequent positron-to-electron ratio will get flattened and the spectrum hardening should be taken into account in the joint fit of the electron/psoitron data otherwise the inferred parameters will be biased. Our joint fits to the latest AMS-02 positron fraction data together with the PAMELA/Fermi-LAT electron/positron spectrum data suggest that the primary electron spectrum hardening is needed in most though not all modelings. The bounds on dark matter models have also been investigated. In the presence of spectrum hardening of primary electrons, the amount of dark-matter-originated electron/positron pairs needed in the modeling is smaller. Even with such a modification, the annihilation channel $\chi\chi \rightarrow \mu^{+}\mu^{-}$ has been tightly constrained by the Fermi-LAT Galactic diffuse emission data. The decay channel $\chi\rightarrow \mu^{+}\mu^{-}$ is found to be viable.

The CDEX Collaboration has been established for direct detection of light dark matter particles, using ultra-low energy threshold p-type point-contact germanium detectors, in China JinPing underground Laboratory (CJPL). The first 1 kg point-contact germanium detector with a sub-keV energy threshold has been tested in a passive shielding system located in CJPL. The outputs from both the point-contact p+ electrode and the outside n+ electrode make it possible to scan the lower energy range of less than 1 keV and at the same time to detect the higher energy range up to 3 MeV. The outputs from both p+ and n+ electrode may also provide a more powerful method for signal discrimination for dark matter experiment. Some key parameters, including energy resolution, dead time, decay times of internal X-rays, and system stability, have been tested and measured. The results show that the 1 kg point-contact germanium detector, together with its shielding system and electronics, can run smoothly with good performances. This detector system will be deployed for dark matter search experiments.

In this paper we study the tidal stripping process for satellite galaxies orbiting around a massive host galaxy, and focus on its dependence on the morphology of both satellite and host galaxy. For this purpose, we use three different morphologies for the satellites: pure disc, pure bulge and a mixture bulge+disc. Two morphologies are used for the host galaxies: bulge+disc and pure bulge. We find that while the spheroidal stellar component experiences a constant power-law like mass removal, the disc is exposed to an exponential mass loss when the tidal radius of the satellite is of the same order of the disc scale length. This dramatic mass loss is able to completely remove the stellar component on time scale of 100 Myears. As a consequence two satellites with the same stellar and dark matter masses, on the same orbit could either retain considerable fraction of their stellar mass after 10 Gyrs or being completely destroyed, depending on their initial stellar morphology. We find that there are two characteristic time scales describing the beginning and the end of the disc removal, whose values are related to the size of the disc. This result can be easily incorporated in semi-analytical models. We also find that the host morphology and the orbital parameters also have an effect on the determining the mass removal, but they are of secondary importance with respect to satellite morphology. We conclude that satellite morphology has a very strong effect on the efficiency of stellar stripping and should be taken into account in modeling galaxy formation and evolution.

The ratios of fluxes of heavy nuclei from sulfur (Z=16) to chromium (Z=24) to the flux of iron were measured by the ATIC-2 experiment. The ratios are decreasing functions of energy from 5 GeV/n to approximately 80 GeV/n, as expected. However, an unexpected sharp upturn in the ratios are observed for energies above 100 GeV/n for all elements from Z=16 to Z=24. Similar upturn but with lower amplitude was also discovered in the ATIC-2 data for the ratio of fluxes of abundant even nuclei (C, O, Ne, Mg, Si) to the flux of iron. Therefore the spectrum of iron is significantly different from the spectra of other abundant even nuclei.

Weakly Interacting Massive Particles (WIMPs) are the candidates of dark matter in our universe. Up to now any direct interaction of WIMP with nuclei has not been observed yet. The exclusion limits of the spin-independent cross section of WIMP-nucleon which have been experimentally obtained is about 10^{-7}pb at high mass region and only 10^{-5}pb} at low mass region. China Jin-Ping underground laboratory CJPL is the deepest underground lab in the world and provides a very promising environment for direct observation of dark matter. The China Dark Matter Experiment (CDEX) experiment is going to directly detect the WIMP flux with high sensitivity in the low mass region. Both CJPL and CDEX have achieved a remarkable progress in recent two years. The CDEX employs a point-contact germanium semi-conductor detector PCGe whose detection threshold is less than 300 eV. We report the measurement results of Muon flux, monitoring of radioactivity and Radon concentration carried out in CJPL, as well describe the structure and performance of the 1 kg PCGe detector CDEX-1 and 10kg detector array CDEX-10 including the detectors, electronics, shielding and cooling systems. Finally we discuss the physics goals of the CDEX-1, CDEX-10 and the future CDEX-1T detectors.

Recently, several groups identified a tentative $\gamma$-ray line signal with energy $\sim 130$ GeV in the central Galaxy from the Fermi-LAT data. %The morphology study shows that the signal is consistent with dark matter %annihilation, but with an offset $\sim 220$ pc ($1.5^{\circ}$) of the %center from the Galactic center Sgr A$^{\star}$, Such a $\gamma-$ray line can be interpreted as the signal of dark matter annihilation. However, the offset $\sim 220$ pc ($1.5^{\circ}$) of the center of the most prominent signal region from the Galactic center Sgr A$^{\star}$ has been thought to challenge the dark matter annihilation interpretation. Considering the fact that such a 130 GeV $\gamma$-ray line signal consists of only $\sim14$ photons, we suggest that the "imperfect" consistency of these photons with the expected dark matter distribution is due to the limited statistics. The offset will be smaller as more signal photons have been collected in the near future. Our Monte Carlo simulation supports the above speculation.

The Hawking evaporation of small black holes formed by the collapse of dark matter at the center of neutron stars plays a key role in loosing the constraint on the mass of asymmetric bosonic non-interacting dark matter particles. Different from previous works we show that such a kind of dark matter is viable in the mass range from 3.3 GeV to ~ 10 TeV, which covers the most attractive regions, including the preferred asymmetric dark matter mass ~ 5.7 GeV as well as the 5-15 GeV range favored by DAMA and CoGeNT.

Dark matter particles may be captured by a star and then thermalized in the star's core. At the end of its life a massive star collapses suddenly and a compact object is formed. The dark matter particles redistribute accordingly. In the inelastic dark matter model, an extended dense dark matter mini-halo surrounding the neutron star may be formed. Such mini-halos may be common in the Galaxy. The electron/positron flux resulting in the annihilation of dark matter particles, however, is unable to give rise to observable signal unless a nascent mini-halo is within a distance \sim a few 0.1 pc from the Earth.

Very old isolated neutron stars and white dwarfs have been suggested to be probes of dark matter. To play such a role, two requests should be fulfilled, i.e., the annihilation luminosity of the captured dark matter particles is above the thermal emission of the cooling compact objects (request-I) and also dominate over the energy output due to the accretion of normal matter onto the compact objects (request-II). Request-I calls for very dense dark matter medium and the critical density sensitively depends on the residual surface temperature of the very old compact objects. The accretion of interstellar/intracluster medium onto the compact objects is governed by the physical properties of the medium and by the magnetization and rotation of the stars and may outshine the signal of dark matter annihilation. Only in a few specific scenarios both requests are satisfied and the compact objects are dark matter burners. The observational challenges are discussed and a possible way to identify the dark matter burners is outlined.

A strong excess in a form of a wide peak in the energy range of 300-800 GeV was discovered in the first measurements of the electron spectrum in the energy range from 20 GeV to 3 TeV by the balloon-borne experiment ATIC (J. Chang et al. Nature, 2008). The experimental data processing and analysis of the electron spectrum with different criteria for selection of electrons, completely independent of the results reported in (J. Chang et al. Nature, 2008) is employed in the present paper. The new independent analysis generally confirms the results of (J. Chang et al. Nature, 2008), but shows that the spectrum in the region of the excess is represented by a number of narrow peaks. The measured spectrum is compared to the spectrum of (J. Chang et al. Nature, 2008) and to the spectrum of the Fermi/LAT experiment.

The final results of processing the data from the balloon-born experiment ATIC-2 (Antarctica, 2002-2003) for the energy spectra of protons and He, C, O, Ne, Mg, Si, and Fe nuclei, the spectrum of all particles, and the mean logarithm of atomic weight of primary cosmic rays as a function of energy are presented. The final results are based on improvement of the methods used earlier, in particular, considerably increased resolution of the charge spectrum. The preliminary conclusions on the significant difference in the spectra of protons and helium nuclei (the proton spectrum is steeper) and the non-power character of the spectra of protons and heavier nuclei (flattening of carbon spectrum at energies above 10 TeV) are confirmed. A complex structure of the energy dependence of the mean logarithm of atomic weight is found.

The data collected by ATIC, PPB-BETS, FERMI-LAT and HESS all indicate that there is an electron/positron excess in the cosmic ray energy spectrum above $\sim$ 100 GeV, although different instrumental teams do not agree on the detailed spectral shape. PAMELA also reported a clear excess feature of the positron fraction above several GeV, but no excess in anti-protons. Here we review the observational status and theoretical models of this interesting observational feature. We pay special attention to various physical interpretations proposed in the literature, including modified supernova remnant models for the $e^\pm$ background, new astrophysical sources, and new physics (the dark matter models). We suggest that although most models can make a case to interpret the data, with the current observational constraints the dark matter interpretations, especially those invoking annihilation, require much more exotic assumptions than some astrophysical interpretations. Future observations may present some ``smoking-gun'' observational tests to differentiate among different models and to identify the correct interpretation to the phenomenon.

Titanium is a rare, secondary nucleus among Galactic cosmic rays. Using the Silicon matrix in the ATIC experiment, Titanium has been separated. The energy dependence of the Ti to Fe flux ratio in the energy region from 5 GeV per nucleon to about 500 GeV per nucleon is presented.

We have observed cosmic-ray electrons from 10 GeV to 800 GeV by a long duration balloon flight using Polar Patrol Balloon (PPB) in Antarctica. The observation was carried out for 13 days at an average altitude of 35 km in January 2004. The detector is an imaging calorimeter composed of scintillating-fiber belts and plastic scintillators inserted between lead plates with 9 radiation lengths. The performance of the detector has been confirmed by the CERN-SPS beam test and also investigated by Monte-Carlo simulations. New telemetry system using a commercial satellite of Iridium, power supply by solar batteries, and automatic level control using CPU have successfully been developed and operated during the flight. From the long duration balloon observations, we derived the energy spectrum of cosmic-ray electrons in the energy range from 100 GeV to 800 GeV. In addition, for the first time we derived the electron arrival directions above 100 GeV, which is consistent with the isotropic distribution.

The ATIC balloon-borne experiment measures the energy spectra of elements from H to Fe in primary cosmic rays from about 100 GeV to 100 TeV. ATIC is comprised of a fully active bismuth germanate calorimeter, a carbon target with embedded scintillator hodoscopes, and a silicon matrix that is used as the main charge detector. The silicon matrix produces good charge resolution for protons and helium but only partial resolution for heavier nuclei. In the present paper, the charge resolution of ATIC was improved and backgrounds were reduced in the region from Be to Si by using the upper layer of the scintillator hodoscope as an additional charge detector. The flux ratios of nuclei B/C, C/O, N/O in the energy region from about 10 GeV/nucleon to 300 GeV/nucleon obtained from this high-resolution, high-quality charge spectra are presented, and compared with existing theoretical predictions.

The results of ATIC-2 baloon experiment (2002-2003) for energy spectra of protons, He, C, O, Ne, Mg, Si, Fe, some groups of nuclei, and all-particle spectrum in primary cosmic rays are presented in energy region 50GeV--200TeV. The conclusion is that the spectra of protons and helium nuclei are essentially different (the spectrum of protons is steeper) and the spectra of protons and heavy nuclei have no-power form.