Context. Pluto's tenuous nitrogen (N2) atmosphere undergoes strong seasonal effects due to high obliquity and orbital eccentricity, and has been recently (July 2015) observed by the New Horizons spacecraft. Goals are (i) construct a well calibrated record of the seasonal evolution of surface pressure on Pluto and (ii) constrain the structure of the lower atmosphere using a central flash observed in 2015. Method: eleven stellar occultations by Pluto observed between 2002 and 2016 are used to retrieve atmospheric profiles (density, pressure, temperature) between $\sim$5 km and $\sim$380 km altitude levels (i.e. pressures from about 10 microbar to 10 nanobar). Results: (i) Pressure has suffered a monotonic increase from 1988 to 2016, that is compared to a seasonal volatile transport model, from which tight constraints on a combination of albedo and emissivity of N2 ice are derived; (ii) A central flash observed on 2015 June 29 is consistent with New Horizons REX profiles, provided that (a) large diurnal temperature variations (not expected by current models) occur over Sputnik Planitia and/or (b) hazes with tangential optical depth of about 0.3 are present at 4-7 km altitude levels and/or (c) the nominal REX density values are overestimated by an implausibly large factor of about 20% and/or (d) higher terrains block part of the flash in the Charon facing hemisphere.

Flarix is a radiation-hydrodynamical (RHD) code for modeling of the response of the chromosphere to a beam bombardment during solar flares. It solves the set of hydrodynamic conservation equations coupled with non-LTE equations of radiative transfer. The simulations are driven by high energy electron beams. We present results of the Flarix simulations of a flaring loop relevant to the problem of continuum radiation during flares. In particular we focus on properties of the hydrogen Balmer continuum which was recently detected by IRIS.

Recent measurements of vertical extents and positions of the chromospheric hard X-ray (HXR) flare sources based on Ramaty High-Energy Spectroscopic Imager (RHESSI) observations show a significant inconsistency with the theoretical predictions based on the standard collisional thick target model (CTTM). Using a hybrid flare code Flarix, we model simultaneously and self-consistently the propagation, scattering and energy losses of electron beams with power-law energy spectra and various initial pitch-angle distributions in a purely collisional approximation and concurrently the dynamic response of the heated chromosphere on timescales typical for RHESSI image reconstruction. The results of the simulations are used to model the time evolution of the vertical distribution of chromospheric HXR source within a singular (compact) loop. Adopting the typical RHESSI imaging times scales, energy dependent vertical sizes and positions as could be observed by RHESSI are presented.

We study two modifications of the collisional thick-target model (CTTM) based on the global and local re-acceleration of non-thermal electrons by static and stochastic electric fields during their transport from the coronal acceleration site to the thick-target region in the chromosphere. We concentrate on a comparison of the non-thermal electron distribution functions, chromospheric energy deposits, and HXR spectra obtained for both considered modifications with the CTTM itself. The results were obtained using a relativistic test-particle approach. We simulated the transport of non-thermal electrons with a power-law spectrum including the influence of scattering, energy losses, magnetic mirroring, and also the effects of the electric fields corresponding to both modifications of the CTTM. We show that both modifications of the CTTM change the outcome of the chromospheric bombardment in several aspects. The modifications lead to an increase in chromospheric energy deposit, change of its spatial distribution, and a substantial increase in the corresponding HXR spectrum intensity.

We investigate the role of non-thermal electrons in the formation regions of Halpha, Hbeta, and Hgamma lines in order to unfold their influence on the formation of these lines. We concentrate on pulse-beam heating varying on a subsecond timescale. Furthermore, we theoretically explore possibility that a new diagnostic tool exists indicating the presence of non-thermal electrons in the flaring chromosphere based on observations of optical hydrogen lines. To model the evolution of the flaring atmosphere and the time-dependent hydrogen excitation and ionisation, we used a 1-D radiative hydrodynamic code combined with a test-particle code that simulates the propagation, scattering, and thermalisation of a power-law electron beam in order to obtain the flare heating and the non-thermal collisional rates due to the interaction of the beam with the hydrogen atoms. All calculated models have shown a time-correlated response of the modelled Balmer line intensities on a subsecond timescale, with a subsecond timelag behind the beam flux. Depending on the beam parameters, both line centres and wings can show pronounced intensity variations. The non-thermal collisional rates generally result in an increased emission from a secondary region formed in the chromosphere.

With the invention of new far-infrared (FIR) and radio mm and sub-mm instruments (DESIR on SMESE satellite, ESO-ALMA), there is a growing interest in observations and analysis of solar flares in this so far unexplored wavelength region. Two principal radiation mechanisms play a role: the synchrotron emission due to accelerated particle beams moving in the magnetic field and the thermal emission due to the energy deposit in the lower atmospheric layers. In this contribution we explore the time-dependent effects of beams on thermal FIR and radio continua. We show how and where these continua are formed in the presence of time dependent beam heating and non-thermal excitation/ionisation of the chromospheric hydrogen plasma.

Using 1D NLTE radiative hydrodynamics we model the influence of the particle beams on the Halpha line profile treating the beam propagation and the atmosphere evolution self-consistently. We focus on the influence of the non-thermal collisional rates and the return current. Based on our results, we propose a diagnostic method for determination of the particle beam presence in the formation regions of the Halpha line.