The paper presents the summary data of the authors' research within the framework of the OLIMPIYA project (the Russian acronym of {\bf OLI}viny iz {\bf M}eteoritov --- {\bf P}oisk tyazholykh {\bf I} sverkhtyazholykh {\bf YA}der / Olivines from meteorites: Search for heavy and superheavy nuclei) and results of track analysis for heavy cosmic ray nuclei (\emph{Z} = 26--129) in olivine crystals from meteorites using an original processing technique. A total of 21,743 tracks of nuclei heavier than iron have been identified in meteoritic matter to date to form the largest database within this charge range. The database includes three tracks of superheavy nuclei with the lifetimes of about a few decades, which can be considered as direct experimental evidence for the existence of natural superheavy nuclei from the "island of stability". Comprehensive comparative analysis of data from two meteorites with different cosmic ray exposure ages, Marjalahti (from 178 to 205 Myr) and Eagle Station (from 35 to 71 Myr), is presented for the first time. The results are discussed within the existing concepts of nuclei formation in astrophysical processes.

Escape of an early atmosphere from Titan, during which time NH3 could be converted by photolysis into the present N2 dominated atmosphere, is an important problem in planetary science. Recently Gilliam and Lerman (2014) estimated escape driven by the surface temperature and pressure, which we show gave loss rates that are orders of magnitude too large. Their model, related to Jeans escape from an isothermal atmosphere, was used to show that escape driven only by surface heating would deplete the atmospheric inventory of N for a suggested Titan accretion temperature of ~355 K. Therefore, they concluded that the accretion temperature must be lower in order to retain the present Titan atmosphere. Here we show that the near surface atmospheric temperature is essentially irrelevant for determining the atmospheric loss rate from Titan and that escape is predominantly driven by solar heating of the upper atmosphere. We also give a rough estimate of the escape rate in the early solar system (~10^4 kg/s) consistent with an inventory of nitrogen being available over the time period suggested by Atreya et al. (1978) for conversion of NH3 into N2.

Observations indicate that some of the largest Kuiper Belt Objects (KBOs) have retained volatiles in the gas phase, which implies the presence of an atmosphere that can affect their reflectance spectra and thermal balance. Volatile escape rates driven by solar heating of the surface were estimated by Schaller and Brown (2007) (SB) and Levi and Podolak (2009)(LP) using Jeans escape from the surface and a hydrodynamic model respectively. Based on recent molecular kinetic simulations these rates can be hugely in error (e.g., a factor of $\sim 10^{16}$ for the SB estimate for Pluto). In this paper we estimate the loss of primordial N$_2$ for several large KBOs guided by recent molecular kinetic simulations of escape due to solar heating of the surface and due to UV/EUV heating of the upper atmosphere. For the latter we extrapolate simulations of escape from Pluto (Erwin et al. 2013) using the energy limited escape model recently validated for the KBOs of interest by molecular kinetic simulations (Johnson et al. 2013). Unless the N$_2$ atmosphere is thin ($\lesssim 10^{18}$ N$_2$/cm$^2$) and/or the radius small ($\lesssim 200-300$ km), we find that escape is primarily driven by the UV/EUV radiation absorbed in the upper atmosphere rather than the solar heating of the surface. This affects the previous interpretations of the relationship between atmospheric loss and the observed surface properties. The long-term goal is to connect detailed atmospheric loss simulations with a model for volatile transport (e.g., Young, 2014) for individual KBOs.

The equations of gas dynamics are extensively used to describe atmospheric loss from solar system bodies and exoplanets even though the boundary conditions at infinity are not uniquely defined. Using molecular-kinetic simulations that correctly treat the transition from the continuum to the rarefied region, we confirm that the energy-limited escape approximation is valid when adiabatic expansion is the dominant cooling process. However, this does not imply that the outflow goes sonic. In fact in the sonic regime, the energy limited approximation can significantly under estimate the escape rate. Rather large escape rates and concomitant adiabatic cooling can produce atmospheres with subsonic flow that are highly extended. Since this affects the heating rate of the upper atmosphere and the interaction with external fields and plasmas, we give a criterion for estimating when the outflow goes transonic in the continuum region. This is applied to early terrestrial atmospheres, exoplanet atmospheres, and the atmosphere of the ex-planet, Pluto, all of which have large escape rates. The paper and its erratum, combined here, are published: ApJL 768, L4 (2013); ApJ, 779, L30 (2013).

A combined fluid/kinetic model is developed to calculate thermally driven escape of N2 from Pluto's atmosphere for two solar heating conditions: no heating above 1450 km and solar minimum heating conditions. In the combined model, one-dimensional fluid equations are applied for the dense part of the atmosphere, while the exobase region is described by a kinetic model and calculated by the direct simulation Monte Carlo method. Fluid and kinetic parts of the model are iteratively solved in order to maintain constant total mass and energy fluxes through the simulation region. Although the atmosphere was found to be highly extended, with an exobase altitude at ~6000 km at solar minimum, the outflow remained subsonic and the escape rate was within a factor of two of the Jeans rate for the exobase temperatures determined. This picture is drastically different from recent predictions obtained solely using a fluid model which, in itself, requires assumptions about atmospheric density, flow velocity and energy flux carried away by escaping molecules at infinity. Gas temperature, density, velocity and heat flux versus radial distance are consistent between the hydrodynamic and kinetic model up to the exobase, only when the energy flux across the lower boundary and escape rate used to solve the hydrodynamic equations is obtained from the kinetic model. This limits the applicability of fluid models to atmospheric escape problems. Finally, the recent discovery of CO at high altitudes, the effect of Charon and the conditions at the New Horizon encounter are briefly considered.

Thermally-driven atmospheric escape evolves from an organized outflow (hydrodynamic escape) to escape on a molecule by molecules basis (Jeans escape) with increasing Jeans parameter, the ratio of the gravitational to thermal energy of molecules in a planet's atmosphere. This transition is described here using the direct simulation Monte Carlo method for a single component spherically symmetric atmosphere. When the heating is predominantly below the lower boundary of the simulation region, R0, and well below the exobase, this transition is shown to occur over a surprisingly narrow range of Jeans parameters evaluated at R0: {\lambda}0 ~ 2-3. The Jeans parameter {\lambda}0 ~ 2.1 roughly corresponds to the upper limit for isentropic, supersonic outflow and for {\lambda}0 >3 escape occurs on a molecule by molecule basis. For {\lambda}0 > ~6, it is shown that the escape rate does not deviate significantly from the familiar Jeans rate evaluated at the nominal exobase, contrary to what has been suggested. Scaling by the Jeans parameter and the Knudsen number, escape calculations for Pluto and an early Earth's atmosphere are evaluated, and the results presented here can be applied to thermally-induced escape from a number of solar and extrasolar planetary bodies.