We review recent advancements in modeling the stellar to substellar transition. The revised molecular opacities, solar oxygen abundances and cloud models allow to reproduce the photometric and spectroscopic properties of this transition to a degree never achieved before, but problems remain in the important M-L transition characteristic of the effective temperature range of characterizable exoplanets. We discuss of the validity of these classical models. We also present new preliminary global Radiation HydroDynamical M dwarfs simulations.

High-resolution images of the solar surface show a granulation pattern of hot rising and cooler downward-sinking material -- the top of the deep-reaching solar convection zone. Convection plays a role for the thermal structure of the solar interior and the dynamo acting there, for the stratification of the photosphere, where most of the visible light is emitted, as well as for the energy budget of the spectacular processes in the chromosphere and corona. Convective stellar atmospheres can be modeled by numerically solving the coupled equations of (magneto)hydrodynamics and non-local radiation transport in the presence of a gravity field. The CO5BOLD code described in this article is designed for so-called "realistic" simulations that take into account the detailed microphysics under the conditions in solar or stellar surface layers (equation-of-state and optical properties of the matter). These simulations indeed deserve the label "realistic" because they reproduce the various observables very well -- with only minor differences between different implementations. The agreement with observations has improved over time and the simulations are now well-established and have been performed for a number of stars. Still, severe challenges are encountered when it comes to extending these simulations to include ideally the entire star or substellar object: the strong stratification leads to completely different conditions in the interior, the photosphere, and the corona. Simulations have to account for a large range of spatial and time scales and also non-equilibrium processes. Last but not least, realistic simulations are based on detailed microphysics and depend on the quality of the input data. This article provides an overview of the physical problem and the numerical solution and the capabilities of CO5BOLD, illustrated with a number of applications.

Three-dimensional magnetohydrodynamic simulations of the surface layers of the Sun intrinsically produce a predominantly horizontal magnetic field in the photosphere. This is a robust result in the sense that it arises from simulations with largely different initial and boundary conditions for the magnetic field. While the disk-center synthetic circular and linear polarization signals agree with measurements from Hinode, their center-to-limb variation sensitively depends on the height variation of the horizontal and the vertical field component and they seem to be at variance with the observed behavior.

Results of realistic simulations of solar surface convection on the scale of supergranules (96 Mm wide by 20 Mm deep) are presented. The simulations cover only 10% of the geometric depth of the solar convection zone, but half its pressure scale heights. They include the hydrogen, first and most of the second helium ionization zones. The horizontal velocity spectrum is a power law and the horizontal size of the dominant convective cells increases with increasing depth. Convection is driven by buoyancy work which is largest close to the surface, but significant over the entire domain. Close to the surface buoyancy driving is balanced by the divergence of the kinetic energy flux, but deeper down it is balanced by dissipation. The damping length of the turbulent kinetic energy is 4 pressure scale heights. The mass mixing length is 1.8 scale heights. Two thirds of the area is upflowing fluid except very close to the surface. The internal (ionization) energy flux is the largest contributor to the convective flux for temperatures less than 40,000 K and the thermal energy flux is the largest contributor at higher temperatures. This data set is useful for validating local helioseismic inversion methods. Sixteen hours of data are available as four hour averages, with two hour cadence, at steinr.msu.edu/~bob/96averages, as idl save files. The variables stored are the density, temperature, sound speed, and three velocity components. In addition, the three velocity components at 200 km above mean continuum optical depth unity are available at 30 sec. cadence.

Observations with the Hinode space observatory led to the discovery of predominantly horizontal magnetic fields in the photosphere of the quiet internetwork region. Here we investigate realistic numerical simulations of the surface layers of the Sun with respect to horizontal magnetic fields and compute the corresponding polarimetric response in the Fe I 630 nm line pair. We find a local maximum in the mean strength of the horizontal field component at a height of around 500 km in the photosphere, where it surpasses the vertical component by a factor of 2.0 or 5.6, depending on the initial and boundary conditions. From the synthesized Stokes profiles we derive a mean horizontal field component that is, respectively, 1.6 and 4.3 times stronger than the vertical component. This is a consequence of both the intrinsically stronger flux density of, and the larger area occupied by the horizontal fields. We find that convective overshooting expels horizontal fields to the upper photosphere, making the Poynting flux positive in the photosphere, while this quantity is negative in the convectively unstable layer below it.

With numerical experiments we explore the feasibility of using high frequency waves for probing the magnetic fields in the photosphere and the chromosphere of the Sun. We track a plane-parallel, monochromatic wave that propagates through a non-stationary, realistic atmosphere, from the convection-zone through the photosphere into the magnetically dominated chromosphere, where it gets refracted and reflected. We compare the wave travel time between two fixed geometrical height levels in the atmosphere (representing the formation height of two spectral lines) with the topography of the surface of equal magnetic and thermal energy density (the magnetic canopy or beta=1 contour) and find good correspondence between the two. We conclude that high frequency waves indeed bear information on the topography of the `magnetic canopy'.