The variational method with constraints recently developed by Verkley and Gerkema to describe maximum-entropy atmospheric profiles is generalized to ideal gases but with temperature-dependent specific heats. In so doing, an extended and non standard potential temperature is introduced that is well suited for tackling the problem under consideration. This new formalism is successfully applied to the atmosphere of Venus. Three well defined regions emerge in this atmosphere up to a height of $100 km$ from the surface: the lowest one up to about $35 km$ is adiabatic, a transition layer located at the height of the cloud deck and finally a third region which is practically isothermal.

We reconsider the scenario in which the knee in the cosmic ray spectrum is explained as due to a change in the escape mechanism of cosmic rays from the Galaxy from one dominated by transverse diffusion to one dominated by drifts. We solve the diffusion equations adopting realistic galactic field models and using diffusion coefficients appropriate for strong turbulence (with a Kolmogorov spectrum of fluctuations) and consistent with the assumed magnetic fields. We show that properly taking into account these effects leads to a natural explanation of the knee in the spectrum, and a transition towards a heavier composition above the knee is predicted.

The increasing cosmic ray statistics collected by present experiments and the future prospects with new large arrays demand accurate calculations of the extensive air shower parameters. The energy of the primary particle is estimated by ground arrays fitting a lateral distribution function (LDF) to the particle densities at a given observing level. However, the lack of appropriate parameterization for these distributions, able to reproduce the data collected from all arrival directions, makes it difficult the experimental analysis. Starting from analytical LDF for vertical showers we study the atmospheric depth dependence of the shower parameters, which is necessary to obtain the corresponding particle density distributions for inclined showers. In our approach both geometrical and atmospheric attenuation effects producing a lateral asymmetry in non-vertical showers are included. The resulting electron and muon LDF fit very well the Monte Carlo simulated data at all zenith angles.

Despite many efforts to find a reasonable explanation, the origin of the "knee" in the cosmic ray spectrum at energy around 10**15.5 eV remains mysterious. In this letter we suggest that the "knee" may be due to a GZK-like effect of cosmic rays interacting with massive neutrinos in the galactic halo. Simple kinematics connects the location of the "knee" with the mass of the neutrinos, and, while the required interaction cross section is larger than that predicted by the Standard Model, it can be accommodated by a small neutrino magnetic dipole moment. The values for the neutrino parameters obtained from the analysis of existing experimental data are compatible with present laboratory bounds.

We explore in some detail the scenario proposed to explain the observed knee of the cosmic ray (CR) spectrum as due to the effects of photodisintegration of the CR nuclei by interactions with optical and soft UV photons in the source region. We show that the photon column densities needed to explain the experimental data are significantly lower than those obtained in previous estimations which neglected multinucleon emission in the photodisintegration process. We also treat more accurately the photodisintegration thresholds, we discuss the effects of photopion production processes and the neutron escape mechanism, identifying the physical processes responsible for the qualitative features of the results. This scenario would require the CR nuclei to traverse column densities of $\sim 5 \times 10^{27}- 2 \times 10^{28}$ eV/cm$^2$ after being accelerated in order to reproduce the observed knee, and predicts that the CR composition should become lighter above $\sim 10^{16}$ eV.

We discuss in detail the possibility of observing pairs of simultaneous parallel air showers produced by the fragments of cosmic ray nuclei which disintegrated in collisions with solar photons. We consider scenarios with different cosmic ray compositions, exploring the predicted rates for existing and planned detectors and looking for methods to extract information on the initial composition from the characteristics of the signal. In particular, we find that fluorescence detectors, such as HiRes or the Telescope Array, due to their low threshold ($\sim 10^{17}$ eV) and large area ($\sim 10^4 km^2$) may observe several events per year if cosmic rays at those energies are indeed heavy nuclei. The possibility of exploiting the angular orientation of the plane containing the two showers to further constrain the cosmic ray composition is also discussed.

We study the propagation of ultra-high energy cosmic ray nuclei through the background of cosmic microwave and intergalactic infrared photons, using recent re-estimates for the density of the last ones. We perform a detailed Monte Carlo simulation to follow the disintegration histories of nuclei starting as Fe and reaching the Earth from extragalactic sources. We obtain the maximum energies of the arriving nuclear fragments as well as the mass composition as a function of the distance traveled. Cosmic rays with energies in excess of $2\times 10^{20}$ eV cannot originate from Fe nuclei produced in sources beyond 10 Mpc.

We show that the photodisintegration of heavy cosmic ray nuclei with energies above 10^20 eV is dominated by interactions with photons from the cosmic microwave background radiation, rather than from infrared ones. This implies that the observed air shower events with energies 2-3 10^20 eV cannot originate from Fe nuclei coming from distances beyond 10 Mpc

We re-examine the interaction of ultra high energy nuclei with the microwave background radiation. We find that the giant dipole resonance leaves a new signature in the differential energy spectrum of iron sources located around 3 Mpc: A depression before the bump which is followed by the expected cutoff.