The Large High Altitude Air Shower Observatory (LHAASO) project is a new generation multi-component instrument, to be built at 4410 meters of altitude in the Sichuan province of China, with the aim to study with unprecedented sensitivity the spec trum, the composition and the anisotropy of cosmic rays in the energy range between 10$^{12}$ and 10$^{18}$ eV, as well as to act simultaneously as a wide aperture (one stereoradiant), continuously-operated gamma ray telescope in the energy range between 10$^{11}$ and $10^{15}$ eV. The experiment will be able of continuously surveying the TeV sky for steady and transient sources from 100 GeV to 1 PeV, t hus opening for the first time the 100-1000 TeV range to the direct observations of the high energy cosmic ray sources. In addition, the different observables (electronic, muonic and Cherenkov/fluorescence components) that will be measured in LHAASO will allow to investigate origin, acceleration and propagation of the radiation through a measurement of energy spec trum, elemental composition and anisotropy with unprecedented resolution. The remarkable sensitivity of LHAASO in cosmic rays physics and gamma astronomy would play a key-role in the comprehensive general program to explore the High Energy Universe. LHAASO will allow important studies of fundamental physics (such as indirect dark matter search, Lorentz invariance violation, quantum gravity) and solar and heliospheric physics. In this document we introduce the concept of LHAASO and the main science goals, providing an overview of the project.

Solar energetic particles acceleration by a shock wave accompanying a coronal mass ejection (CME) is studied. The description of the accelerated particle spectrum evolution is based on the numerical calculation of the diffusive transport equation with a set of realistic parameters. The relation between the CME and the shock speeds, which depend on the initial CME radius, is determined. Depending on the initial CME radius, its speed, and the magnetic energy of the scattering Alfven waves, the accelerated particle spectrum is established during 10-60 minutes from the beginning of CME motion. The maximum energies of particles reach 0.1-10 GeV. The CME radii of 3-5 $R_\odot$ and the shock radii of 5-10 $R_\odot$ agree with observations. The calculated particle spectra agree with the observed ones in events registered by ground-based detectors if the turbulence spectrum in the solar corona significantly differs from the Kolmogorov one.

A nonlinear kinetic theory, combining cosmic-ray (CR) acceleration in supernova remnants (SNRs) with their gas dynamics, is used to re-examine the nonthermal properties of the remnant of SN 1987A for an extended evolutionary period of 5-50 yr. This spherically symmetric model is approximately applied to the different features of the SNR which consist of (i) a blue supergiant wind and bubble, and (ii) of the swept-up red supergiant (RSG) wind structures in the form of an H II region, an equatorial ring (ER), and an hourglass region. The RSG wind involves a mass loss rate that decreases significantly with elevation above and below the equatorial plane. The model adapts recent three-dimensional hydrodynamical simulations by Potter et al. in 2014 that use a significantly smaller ionized mass of the ER than assumed in the earlier studies by the present authors. The SNR shock recently swept up the ER, which is the densest region in the immediate circumstellar environment. Therefore, the expected gamma-ray energy flux density at TeV energies in the current epoch has already reached its maximal value of $\sim 10^{-13}$ erg cm$^{-2}$ s$^{-1}$. This flux should decrease by a factor of about two over the next 10 years.

We present the energy spectrum of antiproton cosmic ray (CR) component calculated on the basis of the nonlinear kinetic model of CR production in supernova remnants (SNR). The model includes reacceleration of already existing in interstellar medium antiprotons as well as creation of antiprotons in nuclear collisions of accelerated protons with gas nuclei and their subsequent acceleration by SNR shock. It is shown that antiprotons production in SNRs produces considerable effect in their resultant energy spectrum making it essentially flatter above 10 GeV so that the spectrum at TeV-energies increases by a factor of five. Calculated antiproton spectrum is well consistent with the PAMELA data, which correspond to energies below 100 GeV. As a consistency check we have also calculated within the same model the energy spectra of secondary nuclei and show that the measured boron-to-carbon ratio is consistent with the significant SNR contribution.

The nature of the recently detected HE and VHE gamma-ray emission of Tycho's supernova remnant (SNR) is studied. A nonlinear kinetic theory of cosmic ray (CR) acceleration in supernova remnants (SNRs) is employed to investigate the properties of Tycho's SNR and their correspondence to the existing experimental data, taking into account that the ambient interstellar medium (ISM) is expected to be clumpy. It is demonstrated that the overall steep gamma-ray spectrum observed can be interpreted as the superposition of two spectra produced by the CR proton component in two different ISM phases: The first gamma-ray component, extending up to about $10^{14}$ eV, originates in the diluted warm ISM, whereas the second component, extending up to 100 GeV, comes from numerous dense, small-scale clouds embedded in this warm ISM. Given the consistency between acceleration theory and the observed properties of the nonthermal emission of Tycho's SNR, a very efficient production of nuclear CRs in Tycho's SNR is established. The excess of the GeV gamma-ray emission due to the clouds' contribution above the level expected in the case of a purely homogeneous ISM, is inevitably expected in the case of type Ia SNe.

The properties of the Galactic supernova remnant (SNR) SN 1006 are theoretically re-analyzed in the light of the recent H.E.S.S. results. Nonlinear kinetic theory is used to determine the momentum spectrum of cosmic rays (CRs) in space and time in the supernova remnant SN 1006. The physical parameters of the model - proton injection rate, electron-to-proton ratio and downstream magnetic field strength - are determined through a fit of the result to the observed spatially-integrated synchrotron emission properties. The only remaining unknown astronomical parameter, the circumstellar gas number density, is determined by a normalization of the amplitude of the gamma-ray flux to the observed amplitude. The bipolar morphology of both nonthermal X-ray and gamma-ray emissions is explained by the preferential injection of suprathermal nuclei and subsequent magnetic field amplification in the quasi-parallel regions of the outer supernova shock. The above parameters provide an improved fit to all existing nonthermal emission data, including the TeV emission spectrum recently detected by H.E.S.S., with the circumstellar hydrogen gas number density $N_H\approx 0.06 cm^{-3}$ close to values derived from observations of thermal X-rays. The hadronic and leptonic gamma-ray emissions are of comparable strength. The overall energy of accelerated CRs at the present epoch is of the order of 5 % of the total hydrodynamic explosion energy, and is predicted to rise with time by a factor of $\approx 2$. The relevance of CR escape from the SNR for the spectrum of the gamma-ray emission is demonstrated. The sum of the results suggests that SN 1006 is a CR source with a high efficiency of nuclear CR production, as required for the Galactic CR sources, both in flux as well as in cutoff energy.

A nonlinear kinetic theory of cosmic ray (CR) acceleration in supernova remnants is employed to re-examine the nonthermal properties of the remnant of SN 1987A for an extended evolutionary period of 5--100 yr. It is shown that an efficient production of nuclear CRs leads to a strong modification of the outer supernova remnant shock and to a large downstream magnetic field $B_\mathrm{d}\approx 20$ mG. The shock modification and the strong field are required to yield the steep radio emission spectrum observed, as well as to considerable synchrotron cooling of high energy electrons which diminishes their X-ray synchrotron flux. These features are also consistent with the existing X-ray observations. The expected \gr energy flux at TeV-energies at the current epoch is nearly $\epsilon_{\gamma}F_{\gamma}\approx 4\times 10^{-13}$ erg cm$^2$s$^{-1}$ under reasonable assumptions about the overall magnetic field topology and the turbulent perturbations of this field. The general nonthermal strength of the source is expected to increase roughly by a factor of two over the next 15 to 20 yrs; thereafter it should decrease with time in a secular form.

The properties of the - presumably - youngest Galactic supernova remnant (SNR) G1.9+0.3 are investigated within the framework of nonlinear kinetic theory of cosmic ray acceleration in SNRs. The observed angular size and expansion speed as well as the radio and X-ray emission measurements are used to determine relevant physical parameters of this SNR. Under the assumption that SNR G1.9+0.3 is the result of a Type Ia supernova near the Galactic center (at the distance d=8.5 kpc) the nonthermal properties are calculated. In particular, the expected TeV gamma-ray spectral energy density is predicted to be as low as $\epsilon_{\gamma}F_{\gamma} \approx 5\times 10^{-15}$ erg cm$^{-2}$ s$^{-1}$, strongly dependent ($F_{\gamma}\propto d^{-11}$) upon the source distance d.

SS433, located at the center of the supernova remnant W50, is a close proximity binary system consisting of a compact star and a normal star. Jets of material are directed outwards from the vicinity of the compact star symmetrically to the east and west. Non-thermal hard X-ray emission is detected from lobes lying on both sides. Shock accelerated electrons are expected to generate sub-TeV gamma rays through the inverse-Compton process in the lobes. Observations of the western X-ray lobe region of SS433/W50 system have been performed to detect sub-TeV gamma-rays using the 10m CANGAROO-II telescope in August and September, 2001, and July and September, 2002. The total observation times are 85.2 hours for ON source, and 80.8 hours for OFF source data. No significant excess of sub-TeV gamma rays has been found at 3 regions of the western X-ray lobe of SS433/W50 system. We have derived 99% confidence level upper limits to the fluxes of gamma rays and have set constraints on the strengths of the magnetic fields assuming the synchrotron/inverse-Compton model for the wide energy range of photon spectrum from radio to TeV. The derived lower limits are 4.3 microgauss for the center of the brightest X-ray emission region and 6.3 microgauss for the far end from SS433 in the western X-ray lobe. In addition, we suggest that the spot-like X-ray emission may provide a major contribution to the hardest X-ray spectrum in the lobe.

The properties of the Galactic supernova remnant SN 1006 are theoretically reanalysed. Nonlinear kinetic theory is used to determine the acceleration efficiency of cosmic rays (CRs) in the supernova remnant SN 1006. The known range of astronomical parameters and the existing measurements of nonthermal emission are examined in order to define the values of the relevant physical parameters which determine the CR acceleration efficiency. It is shown that the parameter values -- proton injection rate, electron to proton ratio and downstream magnetic field strength -- are determined with the appropriate accuracy. In particular also the observed azimuthal variations in the gamma-ray morphology agree with the theoretical expectation. These parameter values, together with the reduction of the gamma-ray flux relative to a spherically symmetric acceleration geometry, allow a good fit to the existing data, including the recently detected TeV emission by H.E.S.S. SN 1006 represents the first example where a high efficiency of nuclear CR production, required for the Galactic CR sources, is consistently established.

We theoretically and phenomenologically investigate the question whether the gamma-ray emission from the remnants of the type Ia supernovae SN 1006, Tycho's SN and Kepler's SN can be the result of electron acceleration alone. The observed synchrotron spectra of the three remnants are used to determine the average momentum distribution of nonthermal electrons as a function of the assumed magnetic field strength. Then the inverse Compton emission spectrum in the Cosmic Microwave Background photon field is calculated and compared with the existing upper limits for the very high energy gamma-ray flux from these sources. It is shown that the expected interstellar magnetic fields substantially overpredict even these gamma-ray upper limits. Only rather strongly amplified magnetic fields could be compatible with such low gamma-ray fluxes. However this would require a strong component of accelerated nuclear particles whose energy density substantially exceeds that of the synchrotron electrons, compatible with existing theoretical acceleration models for nuclear particles and electrons. Even though the quantitative arguments are simplistic, they appear to eliminate simplistic phenomenological claims in favor of a inverse Compton gamma-ray scenario for these sources.

The consequences of a newly suggested value for the SN explosion energy 1.2x10^{51} erg are explored for the case of Tycho's supernova remnant (SNR). A nonlinear kinetic theory of cosmic ray (CR) acceleration in SNRs is employed to investigate the properties of Tycho's SNR and their correspondence to the existing experimental data. It is demonstrated that the large mean ratio between the radii of the contact discontinuity and the forward shock is consistent with the very effective acceleration of nuclear energetic particles at the forward shock. It is also argued that consistency of the value E_{sn} = 1.2x10^{51} erg with the gas dynamics, acceleration theory, and the existing gamma-ray measurements requires the source distance to be greater than 3.3 kpc. The corresponding ambient gas number density is lower than 0.4 cm^{-3}. Since the expected gamma-ray flux strongly depends on the source distance, F_{\gamma}\propto d^{-7}, a future experimental determination of the actual gamma-ray flux from Tycho's SNR will make it possible to determine the values of the source distance and of the mean ambient gas density. A simple inverse Compton model without a dominant population of nuclear CRs is not compatible with the present upper limit for the gamma-ray emission for any reasonable ambient interstellar B-field. Given the consistency between acceleration theory and overall, as well as internal, gas dynamics, a future gamma-ray detection would make the case for nuclear particle acceleration in Tycho's SNR incontrovertible in our view.

For the case of Tycho's supernova remnant (SNR) we present the relation between the blast wave and contact discontinuity radii calculated within the nonlinear kinetic theory of cosmic ray (CR) acceleration in SNRs. It is demonstrated that these radii are confirmed by recently published Chandra measurements which show that the observed contact discontinuity radius is so close to the shock radius that it can only be explained by efficient CR acceleration which in turn makes the medium more compressible. Together with the recently determined new value $E_{sn}=1.2\times 10^{51}$ erg of the SN explosion energy this also confirms our previous conclusion that a TeV gamma-ray flux of $(2-5)\times 10^{-13}$ erg/(cm$^2$s) is to be expected from Tycho's SNR. Chandra measurements and the HEGRA upper limit of the TeV gamma-ray flux together limit the source distance $d$ to $3.3\leq d\leq 4$ kpc.

A nonlinear kinetic theory of cosmic ray (CR) acceleration in supernova remnants is employed to investigate the properties of the remnant SN 1987A. It is shown that a large downstream magnetic field ~10 mG is required to fit the existing observational data. Such a strong field together with the strong shock modification due to CR backreaction provides the steep and concave radioemission spectrum and considerable synchrotron cooling of high energy electrons which diminish their X-ray synchrotron flux below the observed Chandra flux which has to be considered as an upper limit for nonthermal X-ray emission. The expected gamma-ray energy flux at TeV-energies at the current epoch is 2x10^{-13} erg/(cm^2 s).

Nonlinear kinetic theory of cosmic ray (CR) acceleration in supernova remnants (SNRs) is used to investigate the properties of Kepler's SNR and, in particular, to predict the gamma-ray spectrum expected from this SNR. Observations of the nonthermal radio and X-ray emission spectra as well as theoretical constraints for the total supernova (SN) explosion energy E_sn are used to constrain the astronomical and particle acceleration parameters of the system. Under the assumption that Kepler's SN is a type Ia SN we determine for any given explosion energy E_sn and source distance d the mass density of the ambient interstellar medium (ISM) from a fit to the observed SNR size and expansion speed. This makes it possible to make predictions for the expected gamma-ray flux. Exploring the expected distance range we find that for a typical explosion energy E_sn=10^51 erg the expected energy flux of TeV gamma-rays varies from 2x10^{-11} to 10^{-13} erg/(cm^2 s) when the distance changes from d=3.4 kpc to 7 kpc. In all cases the gamma-ray emission is dominated by \pi^0-decay gamma-rays due to nuclear CRs. Therefore Kepler's SNR represents a very promising target for instruments like H.E.S.S., CANGAROO and GLAST. A non-detection of gamma-rays would mean that the actual source distance is larger than 7 kpc.

We use nonlinear kinetic theory to study the remnant dynamics and the particle acceleration as well as the properties of the nonthermal emission from the supernova remnant SN 1006. The known range of astronomical parameters is examined to determine whether it encompasses the existing synchrotron emission data. Given the present-day spatial extent and expansion rate of the object, it is shown that the hadronic gamma-ray flux is very sensitive to the ambient gas density N_H and that the existing H.E.S.S. upper limit requires N_H < 0.1 cm^{-3}. The strength of the amplified magnetic field downstream of the shock is about 150 \mu G.

We have detected gamma-ray emission at the 6sigma level at energies greater than 500GeV from the supernova remnant RX J0852.0-4622 (G266.2-1.2) using the CANGAROO-II Imaging Atmospheric Cherenkov Telescope (IACT). The flux was 0.12 times of that of Crab at 1TeV. The signal centroid is consistent with the peak of the X-ray emission in the north-west rim of the remnant.

It is shown that amplification of the magnetic field in supernova remnants (SNRs) occurs in all six objects where morphological measurements are presently available in the hard X-ray continuum at several keV. For the three archetypical objects (SN 1006, Cas A and Tycho's SNR) to which nonlinear time-dependent acceleration theory has been successfully applied up to now, the global theoretical and the local observational field strengths agree very well, suggesting in addition that all young SNRs exhibit the amplification effect as a result of very efficient acceleration of nuclear cosmic rays (CRs) at the outer shock. Since this appears to be empirically the case, we may reverse the argument and consider field amplification as a measure of nuclear CR acceleration and it has indeed been argued that acceleration in the amplified fields allows the CR spectrum from SNRs to reach the knee in the spectrum or, in special objects, even beyond. The above results are furthermore used to investigate the time evolution of field amplification in young SNRs. Although the uncertainties in the data do not allow precise conclusions regarding this point, they rather clearly show that the ratio of the magnetic field energy density and the kinetic energy density of gas flow into the shock is of the order of a few percent if the shock speed is high enough V_s > 10^3 km/s, and this ratio remains nearly constant during the SNR evolution. The escape of the highest energy nuclear particles from their sources becomes progressively important with age, reducing also the cutoff in the \pi^0 -decay gamma-ray emission spectrum with time after the end of the sweep-up phase. Simultaneously the leptonic gamma-ray channels will gain in relative importance with increasing age of the sources.

We have detected sub-TeV gamma-ray emission from the direction of the Galactic Center (GC) using the CANGAROO-II Imaging Atmospheric Cherenkov Telescope (IACT). We detected a statistically signicant excess at energies greater than 250GeV. The flux was one order of magnitude lower than that of Crab at 1 TeV with a soft spectrum E^{-4.6+-0.5}. The signal centroid is consistent with the GC direction and the observed prole is consistent with a point-like source. Our data suggests that the GeV source 3EG J1746-2851 is identical with this TeV source and we study the combined spectra to determine the possible origin of the gamma-ray emission. We also obtain an upper limit on the cold dark-matter density in the Galactic halo.

It is shown that the nonlinear kinetic model of cosmic ray (CR) acceleration in supernova remnants (SNRs) fits the shell-type nonthermal X-ray morphology, obtained in Chandra observations, in a satisfactory way. The set of empirical parameters is the same which reproduces the dynamical properties of the SNR and the spectral characteristics of the emission produced by CRs. The extremely small spatial scales of the observed X-ray distribution are due to the large effective magnetic field $B_{\mathrm d}\sim 100$ $\mu$G in the interior, which is also required to give a good fit for the spatially integrated radio and X-ray synchrotron spectra. The only reasonably thinkable condition for the production of such a large effective field strength is an efficiently accelerated nuclear CR component. Therefore the Chandra data confirm the inference that SN 1006 indeed accelerates nuclear CRs with the high efficiency required for SNRs to be considered as the main Galactic CR sources.

The injection rate of suprathermal protons into the diffusive shock acceleration process should vary strongly over the surface of supernova remnant shocks. These variations and the absolute value of the injection rate are investigated. In the simplest case, like for SN 1006, the shock can be approximated as being spherical in a uniform large-scale magnetic field. The injection rate depends strongly on the shock obliquity and diminishes as the angle between the ambient field and the shock normal increases. Therefore efficient particle injection, which leads to conversion of a significant fraction of the kinetic energy at a shock surface element, arises only in relatively small regions near the "poles", reducing the overall CR production. The sizes of these regions depend strongly on the random background field and the Alfven wave turbulence generated due to CR streaming instability. For the cases of SN 1006 and Tycho's SNR they correspond to about 20, and for Cas A to between 10 and 20 percent of the entire shock surface. In first approximation, the CR production rate, calculated under the assumption of spherical symmetry, has therefore to be renormalized by this factor, while the shock as such remains roughly spherical.

We study the production of cosmic rays (CRs) in supernova remnants (SNRs), including the reacceleration of background galactic cosmic rays (GCRs) - thus refining the early considerations by Blandford & Ostriker (1980) and Wandel et al. (1987) - and the effects of the nuclear spallation inside the sources (the SNRs). This combines for the first time nuclear spallation inside CR sources and in the diffuse interstellar medium, as well as reacceleration, with the injection and subsequent acceleration of suprathermal particles from the postshock thermal pool. Selfconsistent CR spectra are calculated on the basis of the nonlinear kinetic model. It is shown that GCR reacceleration and CR spallation produce a measurable effect at high energies, especially in the secondary to primary (s/p) ratio, making its energy-dependence substantially flatter than predicted by the standard model. Quantitatively, the effect depends strongly upon the density of the surrounding circumstellar matter. GCR reacceleration dominates secondary CR production for a low circumstellar density. It increases the expected s/p ratio substantially and flattens its spectrum to an almost energy-independent form for energies larger than 100 GeV/n if the supernovae explode on average into a hot dilute medium with hydrogen number density $N_H=0.003$cm$^{-3}$. The contribution of CR spallation inside SNRs to the s/p ratio increases with increasing circumstellar density and becomes dominant for $N_H\gsim 1$ cm$^{-3}$, leading at high energies to a flat s/p ratio which is only by a factor of three lower than in the case of the hot medium. Measurements of the boron to carbon ratio at energies above 100 GeV/n could be used in comparison with the values predicted here as a consistency test for the supernova origin of the GCRs.

We searched for TeV gamma-rays from the remnant of SN 1987A around 5400 days after the supernova. The observations were carried out in 2001, from November 16 to December 11, using the CANGAROO-II Imaging Atmospheric Cherenkov Telescope. In total, 708 minutes of ON- and 1019 minutes of OFF-source data were obtained under good conditions. The detection threshold was estimated to be 1 TeV, due to the mean zenith angle of 39$^\circ$. The upper limits for the gamma-ray flux were obtained and compared with the previous observations and theoretical models. The observations indicate that the gamma-ray luminosity is lower than $1\times 10^{37}$ erg s$^{-1}$ at $\sim 10$ TeV.

A nonlinear kinetic model of cosmic ray (CR) acceleration in supernova remnants (SNRs) is used to describe the properties of Tycho's SNR. Observations of the expansion characteristics and of the nonthermal radio and X-ray emission spectra, assumed to be of synchrotron origin, are used to constrain the overall dynamical evolution and the particle acceleration parameters of the system, in addition to what is known from independent estimates of the distance and thermal X-ray observations. It is shown that a very efficient production of nuclear cosmic rays, leading to strong shock modification, and a large downstream magnetic field strength B_d approx 240muG are required to reproduce the observed synchrotron emission from radio to X-ray frequencies. This field strength is still well within the upper bound for the effective magnetic field, consistent with the acceleration process. The pi^0-decay gamma-ray flux turns out to be somewhat greater than the inverse Compton (IC) flux off the Cosmic Microwave Background at energies below 1 TeV, dominating it strongly at 10 TeV. The predicted TeV gamma-ray flux is consistent with but close to the very low upper limit recently obtained by HEGRA. A future detection at E_gamma ~ 10 TeV would clearly indicate hadronic emission.

The nonlinear kinetic model of cosmic ray (CR) acceleration in supernova remnants (SNRs) is used to describe the properties of the remnant of SN 1006. It is shown, that the theory fits the existing data in a satisfactory way within a set of parameters which is consistent with the idea that SN 1006 is a typical source of Galactic CR nucleons, although not necessarily of CR electrons. The adjusted parameters are those that are not very well determined by present theory or not directly amenable to astronomical observations. The calculated expansion law and the radio-, X-ray and gamma-ray emissions produced by the accelerated CRs in SN 1006 agree quite well with the observations. A rather large interior magnetic field B_d\approx 100 \muG is required to give a good fit for the radio and X-ray synchrotron emission. In the predicted TeV gamma-ray flux from SN 1006, the pi^0-decay gamma-rays, generated by the nuclear CR component, dominate over the inverse Compton (IC) gamma-rays, generated by the CR electrons in the cosmic microwave background. The predicted source morphology in high energy gamma-rays roughly corresponds to that of the synchrotron emission. The predicted integral gamma-ray flux F \propto e^-1 extends up to energies 100 TeV if CR diffusion is as strong as the Bohm limit. Only if the interior magnetic field is much lower in the SNR, B_d\approx 10 \muG, then the observed gamma-ray emission can be due to the accelerated electron component alone. In this case, not plausible physically in our view, the lowest permissible value of the electron to proton ratio is rather high, and the maximum individual energy and total energy content of accelerated nucleons so small, that SN 1006 can not be considered as a typical source of the nuclear Galactic CRs.

Spectrum and chemical composition of cosmic rays accelerated in supernova remnants are studied on a basis of kinetic approach. The cosmic ray transport equation with the Bohm diffusion coefficient has been numerically solved self-consistently with gas dynamic equations for the underlying flow. Comparison with observational results gives some indications that galactic cosmic rays are produced by supernovae shocks which expanded in the relatively low temperature partially ionized interstellar medium.