The analyses of sunspot observations revealed a fundamental magnetic property of the umbral boundary, the invariance of the vertical component of the magnetic field. We aim to analyse the magnetic properties of the umbra-penumbra boundary in simulated sunspots and thus assess their similarity to observed sunspots. Also, we aim to investigate the role of plasma $\beta$ and the ratio of kinetic to magnetic energy in simulated sunspots on the convective motions. We use a set of non-grey simulation runs of sunspots with the MURaM code. These data are used to synthesise the Stokes profiles that are then degraded to the Hinode spectropolarimeter-like observations. Then, the data are treated like real Hinode observations of a sunspot and magnetic properties at the umbral boundaries are determined. Simulations with potential field extrapolation produce a realistic magnetic field configuration on their umbral boundaries. Two simulations with potential field upper boundary, but different subsurface magnetic field structures, differ significantly in the extent of their penumbrae. Increasing the penumbra width by forcing more horizontal magnetic fields at the upper boundary results in magnetic properties that are not consistent with observations. This implies that the size of the penumbra is given by the subsurface structure of the magnetic field. None of the sunspot simulations is consistent with observed properties of the magnetic field and direction of the Evershed flow at the same time. Strong outward directed Evershed flows are only found in setups with artificially enhanced horizontal component of the magnetic field at the top boundary that are not consistent with the observed magnetic field properties at the UP boundary. We want to stress out that the `photospheric' boundary of simulated sunspots is defined by a magnetic field strength of equipartition field value.

Context. In a recent statistical study of sunspots in 79 active regions, the vertical magnetic field component $B_\text{ver}$ averaged along the umbral boundary is found to be independent of sunspot size. The authors of that study conclude that the absolute value of $B_\text{ver}$ at the umbral boundary is the same for all spots. Aims. We investigate the temporal evolution of $B_\text{ver}$ averaged along the umbral boundary of one long-lived sunspot during its stable phase. Methods. We analysed data from the HMI instrument on-board SDO. Contours of continuum intensity at $I_\text{c}=0.5I_\text{qs}$, whereby $I_\text{qs}$ refers to the average over the quiet sun areas, are used to extract the magnetic field along the umbral boundary. Projection effects due to different formation heights of the Fe I 617.3 nm line and continuum are taken into account. To avoid limb artefacts, the spot is only analysed for heliocentric angles smaller than $60^{\circ}$. Results. During the first disc passage, NOAA AR 11591, $B_\text{ver}$ remains constant at 1693 G with a root-mean-square deviation of 15 G, whereas the magnetic field strength varies substantially (mean 2171 G, rms of 48 G) and shows a long term variation. Compensating for formation height has little influence on the mean value along each contour, but reduces the variations along the contour when away from disc centre, yielding a better match between the contours of $B_\text{ver}=1693$ G and $I_\text{c}=0.5I_\text{qs}$. Conclusions. During the disc passage of a stable sunspot, its umbral boundary can equivalently be defined by using the continuum intensity $I_\text{c}$ or the vertical magnetic field component $B_\text{ver}$. Contours of fixed magnetic field strength fail to outline the umbral boundary.