New data on the absolute value on the inelastic proton-air cross section at 10 PeV are presented. Conclusion is made about the growth with energy of the inelastic proton-air cross section according to comparisons with experimental data were obtained at the Tien Shan complex array on various components of extensive air showers: hadrons, Cherenkov light and electrons with many different calculated models of cosmic rays interactions at the atmosphere. The analysis showed that the rise conforms to 7-9 % per one order of the energy from 0.2 TeV (accelerators with fixed targets) to 10 PeV (EAS cosmic rays). That corresponds to around 350 mb at 1 PeV and 380 mb at 10 PeV of primary cosmic rays. These data corresponds better to the QGSJET-II-04 model version.

The events of anomalous energy deposit and high number of hadrons in hadron calorimeter, often not accompanied at all, or accompanied by very small EAS of cosmic rays, are observed by means of the complex installation aimed at research of cosmic rays. The data quoted and their analysis point to the real contribution and the essential role of pion Bose condensation in high energy nuclear interactions. The possibility of the existence of such process in the nuclear matter (in the neutron stars and in some kinds of the nuclei) was proposed by A.B.Migdal in 1971.

Experimental data from the Tien Shan complex array on different components of extensive air showers at 0.5-10 PeV primary cosmic rays are compared with results of various calculated models of cosmic rays interactions at the atmosphere. Conclusion is made about the growth with energy of the inelastic proton-air cross section {\sigma}p-air from 0.2 TeV (accelerator experiments with fixed targets) to 10 PeV (cosmic rays). The analysis showed that the rise conforms to (7-9)% per one order of energy. That corresponds to {\sigma}p-air (1 PeV) = 350 mb. These data correspond better to the new QGSJET-II-04 version of the interaction model based on the recent LHC results. This model predicts better the slower rise of the cross section than previous versions of QGSJET-II and some other models.

The primary cosmic ray sources are searched by means of CHRONOTRON - KLARA separate array of the P.N. Lebedev Physical Institute Tien Shan station. It was done on the base of 35 millions observed PCR extensive air showers from 5x10**13 to 5x10**14 eV energies. The data analysis was carried on the method of the direct selection of local areas in equatorial coordinates where the deviation of event numbers exceeded the definite value from normal Gaussian standard. These directions are compared with other arrays observed results and with coordinates of astrophysical sources.

Results of a blind search for localized regions of an excessive flux of cosmic rays in the energy range from 50 TeV to 20 PeV with the data of the FIAN KLARA-Chronotron experiment, the EAS MSU array and the Prototype of the EAS-1000 array are presented. A number of regions with a significant excess of the registered flux over an expected isotropic background are found. Some of the regions are present in at least two of the data sets considered.

With the help of the experimental installation aimed at the investigation of high energy cosmic rays (Tien-Shan high mountain laboratory) the signal of Solar and star optical radiation is registered. The signal is well provided statistically and possesses the strictly expressed maximum in the region of EAS sizes Ne 1.19 106 particles (primary energy Eo 1.33 1015 eV). This signal is the peak from gamma EAS, generated by gamma quanta from decay of pi zero mesons, photo produced by the Primary Cosmic Radiation (PCR) nuclei on the photons of stars and of PCR sources. The assumption is made, that exactly this process provides the main contribution in the formation of so called knee on the primary spectrum. Due to the universality and distinct maximum of this signal, its usage for independent and reliable calibration of the EAS installations, for the mutual calibration of these installations and, possibly, for the merger of experimental data obtained by means of these installations to increase the statistics, is proposed. It is especially vital to day, when the further essential increase of the energy range under investigation is necessary.

Bursts of high energy EAS intensity (the series of EAS), following each other in short intervals of time were observed by means of Tien-Shan high mountain installation. For the lower boundary of EAS, uniting in one series, the size Ne= 10**6 (primary energies of the order of 4x10**15 eV) was taken. The condition of amalgamation into one series was the presence of at least two EAS of Ne >10**7. The number of EAS in a series is from 4 to 9 events, with the mean time interval between them 1-5 minutes. Five such series were found in the material treated (approximately 250 days of pure time of the installation run). For each EAS of each series, all the basic parameters are given: observation date and time, age parameter S, galactic coordinates, coordinates of EAS axes relative to the installation center, energy release in the calorimeter and its distance from EAS axes, fitting the Nishimura- Kamata- Greisen (NKG) function. It is essential, that the frequency of the appearing of such series is much higher, than the probability of the occasional fluctuations in the succession of independent EAS. This probably is in favor of the commonality of their origin.

Analysis of the experimental data of joint Russian - Hungarian research was carried out at Tien Shan station on the search of sources of primary cosmic rays (PCR) at energies 0.1-1.0 PeV. Maps of extensive air shower arrival directions were presented where the standard deviation of mean event numbers in equatorial co-ordinates exceeded a definite value. The excess of EAS number was observed from an Galaxy object - the SN remnant and some other directions. The young pulsar PSR 0656+14 located near the centre of this object.

The dependence E_h (N_e) of Extensive Air Shower (EAS) hadronic component energy flux on the number N_e of particles in EAS is investigated in the primary energy range of the order of 10^{15} - 10^{17} eV. The work was aimed at checking the existence of irregularities of E_h (N_e)/N_e behavior at these energies in several independent experiments. The investigation is carried out using large statistical material obtained at different configurations of experimental apparatus and under different triggering conditions. The existence of irregularities of E_h (N_e)/N_e behavior in the region Ne > 2*10^6 is confirmed. These irregularities have the character of sharp deeps and are located near the same values of N_e regardless of the experimental material and selection conditions used. So, at recent stage of research the existence of these irregularities of E_h (N_e)/N_e behavior in the range of N_e > 2*10^6 may be regarded as reliably established. This fact supports our earlier conclusion on the existence of primary cosmic radiation (PCR) nuclei spectra cutoff effect in the primary energy region 10^{15} - 10^{17} eV.