We simulate the stresses induced by temperature changes in a putative hard layer near the surface of comet 67P/Churyumov--Gerasimenko with a thermo-viscoelastic model. Such a layer could be formed by the recondensation or sintering of water ice (and dust grains), as suggested by laboratory experiments and computer simulations, and would explain the high compressive strength encountered by experiments on board the Philae lander. Changes in temperature from seasonal insolation variation penetrate into the comet's surface to depths controlled by the thermal inertia, causing the material to expand and contract. Modelling this with a Maxwellian viscoelastic response on a spherical nucleus, we show that a hard, icy layer with similar properties to Martian permafrost will experience high stresses: up to tens of MPa, which exceed its material strength (a few MPa), down to depths of centimetres to a metre. The stress distribution with latitude is confirmed qualitatively when taking into account the comet's complex shape but neglecting thermal inertia. Stress is found to be comparable to the material strength everywhere for sufficient thermal inertia ($\gtrsim50$ J m$^{-2}$ K$^{-1}$ s$^{-1/2}$) and ice content ($\gtrsim 45\%$ at the equator). In this case, stresses penetrate to a typical depth of $\sim0.25$ m, consistent with the detection of metre-scale thermal contraction crack polygons all over the comet. Thermal fracturing may be an important erosion process on cometary surfaces which breaks down material and weakens cliffs.

We study the effect of angular momentum on the surface density profiles of disc galaxies, using high resolution simulations of major mergers whose remnants have downbending radial density profiles (type II). As described in the previous papers of this series, in this scenario, most of the disc mass is acquired after the collision via accretion from a hot gaseous halo. We find that the inner and outer disc scalelengths, as well as the break radius, correlate with the total angular momentum of the initial merging system, and are larger for high angular momentum systems. We follow the angular momentum redistribution in our simulated galaxies, and find that, like the mass, the disc angular momentum is acquired via accretion, i.e. to the detriment of the gaseous halo. Furthermore, high angular momentum systems give more angular momentum to their discs, which affects directly their radial density profile. Adding simulations of isolated galaxies to our sample, we find that the correlations are valid also for disc galaxies evolved in isolation. We show that the outer part of the disc at the end of the simulation is populated mainly by inside-out stellar migration, and that in galaxies with higher angular momentum, stars travel radially further out. This, however, does not mean that outer disc stars (in type II discs) were mostly born in the inner disc. Indeed, generally the break radius increases over time, and not taking this into account leads to overestimating the number of stars born in the inner disc.

Context: In a series of papers, we study the major merger of two disk galaxies in order to establish whether or not such a merger can produce a disc galaxy. Aims: Our aim here is to describe in detail the technical aspects of our numerical experiments. Methods: We discuss the initial conditions of our major merger, which consist of two protogalaxies on a collision orbit. We show that such merger simulations can produce a non-realistic central mass concentration, and we propose simple, parametric, AGN-like feedback as a solution to this problem. Our AGN-like feedback algorithm is very simple: at each time-step we take all particles whose local volume density is above a given threshold value and increase their temperature to a preset value. We also compare the GADGET3 and GIZMO codes, by applying both of them to the same initial conditions. Results: We show that the evolution of isolated protogalaxies resembles the evolution of disk galaxies, thus arguing that our protogalaxies are well suited for our merger simulations. We demonstrate that the problem with the unphysical central mass concentration in our merger simulations is further aggravated when we increase the resolution. We show that our AGN-like feedback removes this non-physical central mass concentration, and thus allows the formation of realistic bars. Note that our AGN-like feedback mainly affects the central region of a model, without significantly modifying the rest of the galaxy. We demonstrate that, in the context of our kind of simulation, GADGET3 gives results which are very similar to those obtained with the PSPH (density independent SPH) flavor of GIZMO. Moreover, in the examples we tried, the differences between the results of the two flavors of GIZMO, namely PSPH, and MFM (mesh-less algorithm) are similar to and, in some comparisons, larger than the differences between the results of GADGET3 and PSPH.

We use N-body chemo-dynamic simulations to study the coupling between morphology, kinematics and metallicity of the bar/bulge region of our Galaxy. We make qualitative comparisons of our results with available observations and find very good agreement. We conclude that this region is complex, since it comprises several stellar components with different properties -- i.e. a boxy/peanut bulge, thin and thick disc components, and, to lesser extents, a disky pseudobulge, a stellar halo and a small classical bulge -- all cohabiting in dynamical equilibrium. Our models show strong links between kinematics and metallicity, or morphology and metallicity, as already suggested by a number of recent observations. We discuss and explain these links.

Using three fiducial Nbody+SPH simulations, we follow the merging of two disk galaxies with a hot gaseous halo component each, and examine whether the merger remnant can be a spiral galaxy. The stellar progenitor disks are destroyed by violent relaxation during the merging and most of their stars form a classical bulge, while the remaining form a thick disk and its bar. A new stellar disk forms subsequently and gradually in the remnant from the gas accreted mainly from the halo. It is vertically thin and well extended in its equatorial plane. A bar starts forming before the disk is fully in place, contrary to what is assumed in idealised simulations of isolated bar-forming galaxies. It has morphological features such as ansae and boxy/peanut bulges. Stars of different ages populate different parts of the box/peanut. A disky pseudobulge forms also, so that by the end of the simulation, all three types of bulges coexist. The oldest stars are found in the classical bulge, followed by those of the thick disk, then by those in the thin disk. The youngest stars are in the spiral arms and the disky pseudobulge. The disk surface density profiles are of type II (exponential with downbending), and the circular velocity curves are flat and show that the disks are submaximum in these examples: two clearly so and one near-borderline between maximum and submaximum. On average, only roughly between 10 and 20% of the stellar mass is in the classical bulge of the final models, i.e. much less than in previous simulations.

The statistical validation of transiting exoplanets proved to be an efficient technique to secure the nature of small exoplanet signals which cannot be established by purely spectroscopic means. However, the spectroscopic diagnoses are providing us with useful constraints on the presence of blended stellar contaminants. In this paper, we present how a contaminating star affects the measurements of the various spectroscopic diagnoses as function of the parameters of the target and contaminating stars using the model implemented into the PASTIS planet-validation software. We find particular cases for which a blend might produce a large radial velocity signal but no bisector variation. It might also produce a bisector variation anti-correlated with the radial velocity one, as in the case of stellar spots. In those cases, the full width half maximum variation provides complementary constraints. These results can be used to constrain blend scenarios for transiting planet candidates or radial velocity planets. We review all the spectroscopic diagnoses reported in the literature so far, especially the ones to monitor the line asymmetry. We estimate their uncertainty and compare their sensitivity to blends. Based on that, we recommend the use of BiGauss which is the most sensitive diagnosis to monitor line-profile asymmetry. In this paper, we also investigate the sensitivity of the radial velocities to constrain blend scenarios and develop a formalism to estimate the level of dilution of a blended signal. Finally, we apply our blend model to re-analyse the spectroscopic diagnoses of HD16702, an unresolved face-on binary which exhibits bisector variations.

We propose a new algorithm for an adaptive optics system control law, based on the Linear Quadratic Gaussian approach and a Kalman Filter adaptation with localizations. It allows to handle non-stationary behaviors, to obtain performance close to the optimality defined with the residual phase variance minimization criterion, and to reduce the computational burden with an intrinsically parallel implementation on the Extremely Large Telescopes (ELTs).

We propose a new algorithm for an adaptive optics system control law which allows to reduce the computational burden in the case of an Extremely Large Telescope (ELT) and to deal with non-stationary behaviors of the turbulence. This approach, using Ensemble Transform Kalman Filter and localizations by domain decomposition is called the local ETKF: the pupil of the telescope is split up into various local domains and calculations for the update estimate of the turbulent phase on each domain are performed independently. This data assimilation scheme enables parallel computation of markedly less data during this update step. This adapts the Kalman Filter to large scale systems with a non-stationary turbulence model when the explicit storage and manipulation of extremely large covariance matrices are impossible. First simulation results are given in order to assess the theoretical analysis and to demonstrate the potentiality of this new control law for complex adaptive optics systems on ELTs.

We demonstrate that in N-body simulations of isolated disc galaxies there is numerical vertical heating which slowly increases the vertical velocity dispersion and the disc thickness. Even for models with over a million particles in a disc, this heating can be significant. Such an effect is just the same as in numerical experiments by Sellwood (2013). We also show that in a stellar disc, outside a boxy/peanut bulge, if it presents, the saturation level of the bending instability is rather close to the value predicted by the linear theory. We pay attention to the fact that the bending instability develops and decays very fast, so it couldn't play any role in secular vertical heating. However the bending instability defines the minimal value of the ratio between the vertical and radial velocity dispersions $\sigma_z / \sigma_R \approx 0.3$ (so indirectly the minimal thickness) which could have stellar discs in real galaxies. We demonstrate that observations confirm last statement.

We follow the formation and evolution of bars in N-body simulations of disc galaxies with gas and/or a triaxial halo. We find that both the relative gas fraction and the halo shape play a major role in the formation and evolution of the bar. In gas-rich simulations, the disc stays near-axisymmetric much longer than in gas-poor ones, and, when the bar starts growing, it does so at a much slower rate. Due to these two effects combined, large-scale bars form much later in gas-rich than in gas-poor discs. This can explain the observation that bars are in place earlier in massive red disc galaxies than in blue spirals. We also find that the morphological characteristics in the bar region are strongly influenced by the gas fraction. In particular, the bar at the end of the simulation is much weaker in gas-rich cases. In no case did we witness bar destruction. Halo triaxiality has a dual influence on bar strength. In the very early stages of the simulation it induces bar formation to start earlier. On the other hand, during the later, secular evolution phase, triaxial haloes lead to considerably less increase of the bar strength than spherical ones. The shape of the halo evolves considerably with time. The inner halo parts may become more elongated, or more spherical, depending on the bar strength. The main body of initially triaxial haloes evolves towards sphericity, but in initially strongly triaxial cases it stops well short of becoming spherical. Part of the angular momentum absorbed by the halo generates considerable rotation of the halo particles that stay located relatively near the disc for long periods of time. Another part generates halo bulk rotation, which, contrary to that of the bar, increases with time but stays small.

The baryonic discs of galaxies are believed to alter the shapes of the dark matter haloes in which they reside. We perform a set of hydrodynamical N-body simulations of disc galaxies with triaxial dark matter haloes, using elliptical discs with a gaseous component as initial conditions. We explore models of different halo triaxiality and also of different initial gas fractions, which allows us to evaluate how each affects the formations of the bar. Due to star formation, models of all halo shapes and of all initial gas fractions reach approximately the same gas content at the end of the simulation. Nevertheless, we find that the presence of gas in the early phases has important effects on the subsequent evolution. Bars are generally weaker for larger initial gas content and for larger halo triaxiality. The presence of gas, however, is a more efficient factor in inhibiting the formation of a strong bar than halo triaxiality is.

We present results of a strong-lensing analysis of MACS J0717.5+3745 (hereafter MACS J0717), an extremely X-ray luminous galaxy cluster at z=0.55. Observations at different wavelengths reveal a complex and dynamically very active cluster, whose core is connected to a large scale filament extended over several Mpc. Using multi-passband imaging data obtained with the Hubble Space Telescope Advanced Camera for Surveys (ACS), we identify 15 multiply imaged systems across the full field of view of ACS, five of which we confirmed spectroscopically in groundbased follow-up observations with the Keck telescope. We use these multiply imaged systems to constrain a parametric model of the mass distribution in the cluster core, employing a new parallelized version of the Lenstool software. The main result is that the most probable description of the mass distribution comprises four cluster-scale dark matter haloes. The total mass distribution follows the light distribution but strongly deviates from the distribution of the intra-cluster gas as traced by the X-ray surface brightness. This confirms the complex morphology proposed by previous studies. We interpret this segregation of collisional and collisionless matter as strong evidence of multiple mergers and ongoing dynamical activity. MACS J0717 thus constitutes one of the most disturbed clusters presently known and, featuring a projected mass within the ACS field of view (R=150"=960 kpc) of 2.11 +- 0.23 x 10^{15} M_{sun}, the system is also one of the most massive known.

Aims: We aim to develop an algorithm for constructing equilibrium initial conditions for simulations of disk galaxies with a triaxial halo and/or a gaseous component. This will pave the way for N-body simulations of realistic disk galaxies. Methods: We use the iterative method, which we presented in a previous article. The idea of this method is very simple. It relies on constrained evolution. Results: We develop an algorithm for constructing equilibrium models of disk galaxies including a gaseous disk and a triaxial or axisymmetric halo. We discuss two test models. The first model consists of a spherical halo, a stellar disk, and an isothermal gaseous disk. The second model consists of a triaxial halo, a stellar disk, and a star-forming gaseous disk. We demonstrate that both test models are very close to equilibrium, as we had intended.

We present dynamical models of NGC 4494, which we built using our iterative method presented in a previous paper. These models are live N-body models consisting of equal mass particles, and they are steady state as confirmed by a fully self-consistent evolution. Our goals were twofold. The first one -- namely to test whether our iterative method could indeed be used to construct galactic models following given observational constraints, both photometric and kinematic -- was fully achieved. Our method allowed us to go beyond a simple spherical model and to make full sets of rotating, axisymmetric models without any limitations to the velocity distribution. Our second goal was to understand better the structure of NGC 4494, and more specifically to set constraints on its halo mass. For this we tried three families of models: without halo, with a light halo and with a heavy halo, respectively. Our models reproduce well the photometry and the kinematics, the latter except specific regions where some non-equilibrium or non-axisymmetric structure could be present in the galaxy (e.g. the kinematically decoupled core). However, the lower order moments of the velocity distribution (up to and including the second order) do not allow us to discriminate between the three halos. On the other hand, when we extend the comparison to the higher order moments of the velocity distribution obtained from the long-slit data, we find that our light halo model fits the data better than the no halo, or the heavy halo models. They also reproduce the shape of the angular dependence of the PNe velocity dispersion in the outermost parts of the galaxy, but not the amplitude of its azimuthal variation. This may imply that a yet more general class of models, such as triaxial, may be necessary for a yet better fit.

We present a new method for constructing equilibrium phase models for stellar systems, which we call the iterative method. It relies on constrained, or guided evolution, so that the equilibrium solution has a number of desired parameters and/or constraints. This method is very powerful, to a large extent due to its simplicity. It can be used for mass distributions with an arbitrary geometry and a large variety of kinematical constraints. We present several examples illustrating it. Applications of this method include the creation of initial conditions for N-body simulations and the modelling of galaxies from their photometric and kinematic observations.

An iterative approach is used to construct spherically symmetric equilibrium models with an anisotropic velocity distribution. The potentialities of the method have been tested on models with known distribution functions, the Osipkov-Merritt models. It is shown that models that differ significantly from the Osipkov-Merritt models can be constructed. An N-body model of a dark halo with a density distribution that approximates the results of cosmological simulations (the Navarro--Frenk--White model) has been constructed. The anisotropy profile has been taken to be similar to that yielded by cosmological simulations. The constructed models can serve as direct input data for investigating the dynamics and stability of such systems in N-body simulations.

We present a new iterative method for constructing equilibrium phase models of stellar systems. Importantly, this method can provide phase models with arbitrary mass distributions. The method is based on the following principle. Our task is to generate an equilibrium N-body system with a given mass distribution. For this purpose, we let the system reach equilibrium through its dynamical evolution. During this evolution we hold mass distribution in this system. This principle is realized in our method by means of an iterative procedure. We have used our method to construct a phase model of the disc of our Galaxy. In our method, we use the mass distribution in the Galaxy as input data. Here we used two Galactic density models (suggested by Flynn, Sommer-Larsen and Christensen and by Dehnen and Binney). For a fixed-mass model of the Galaxy we can construct a one-parameter family of equilibrium models of the Galactic disc. We can, however, choose a unique model using local kinematic parameters that are known from Hipparcos data. We show that the phase models constructed using our method are close to equilibrium. The problem of uniqueness for our models is discussed, and we discuss some further applications of our method.

One widely used technique for the construction of equilibrium models of stellar disks is based on the Jeans equations and the moments of velocity distribution functions derived using these equations. Stellar disks constructed using this technique are shown to be "not entirely" in equilibrium. Our attempt to abandon the epicyclic approximation and the approximation of infinite isothermal layers, which are commonly adopted in this technique, failed to improve the situation substantially. We conclude that the main drawback of techniques based on the Jeans equations is that the system of equations employed is not closed, and therefore requires adopting an essentially ad hoc additional closure condition. A new iterative approach to constructing equilibrium N-body models with a given density distribution is proposed. The main idea behind this approach is that a model is first constructed using some approximation method, and is then allowed to adjust to an equilibrium state with the specified density and the required parameters of the velocity distribution remaining fixed in the process. This iterative approach was used to construct isotropic, spherically symmetric models and models of stellar disks embedded in an external potential. The numerical models constructed prove to be close to equilibrium. It is shown that the commonly adopted assumption that the profile of the radial velocity dispersion is exponential may be wrong. The technique proposed can be applied to a wide range of problems involving the construction of models of stellar systems with various geometries.

We analyze the relationship between the mass of a spherical component and the minimum possible thickness of stable stellar disks. This relationship for real galaxies allows the lower limit on the dark halo mass to be estimated (the thinner the stable stellar disk is, the more massive the dark halo must be). In our analysis, we use both theoretical relations and numerical N-body simulations of the dynamical evolution of thin disks in the presence of spherical components with different density profiles and different masses. We conclude that the theoretical relationship between the thickness of disk galaxies and the mass of their spherical components is a lower envelope for the model data points. We recommend using this theoretical relationship to estimate the lower limit for the dark halo mass in galaxies. The estimate obtained turns out to be weak. Even for the thinnest galaxies, the dark halo mass within four exponential disk scale lengths must be more than one stellar disk mass.

A criterion for the choice of optimal softening length $\epsilon$ and time-step $dt$ for $N$-body simulations of a collisionless stellar system is analyzed. Plummer and Hernquist spheres are used as models to follow how changes in various parameters of an initially equilibrium stable model depend on $\epsilon$ and $dt$. These dependences are used to derive a criterion for choosing $\epsilon$ and $dt$. The resulting criterion is compared to Merritt's criterion for choosing the softening length, which is based on minimizing the mean irregular force acting on a particle with unit mass. Our criterion for choosing $\epsilon$ and $dt$ indicate that $\epsilon$ must be a factor of 1.5-2 smaller than the mean distance between particles in the densest regions to be resolved. The time-step must always be adjusted to the chosen $\epsilon$ (the particle must, on average, travel a distance smaller than $0.5\epsilon$ during one time-step). An algorithm for solving N-body problems with adaptive variations of the softening length is discussed in connection with the task of choosing $\epsilon$, but is found not to be promising.

The saturation conditions for bending modes in inhomogeneous thin stellar disks that follow from an analysis of the dispersion relation are compared with those derived from $N$-body simulations. In the central regions of inhomogeneous disks, the reserve of disk strength against the growth of bending instability is smaller than that for a homogeneous layer. The spheroidal component (a dark halo, a bulge) is shown to have a stabilizing effect. The latter turns out to depend not only on the total mass of the spherical component, but also on the degree of mass concentration toward the center. We conclude that the presence of a compact (not necessarily massive) bulge in spiral galaxies may prove to be enough to suppress the bending perturbations that increase the disk thickness. This conclusion is corroborated by our $N$-body simulations in which we simulated the evolution of almost equilibrium, but unstable finite-thickness disks in the presence of spheroidal components. The final disk thickness at the same total mass of the spherical component (dark halo + bulge) has been found to be much smaller than that in the simulations where a concentrated bulge is present.

We investigate the nonlinear growth stages of bending instability in stellar disks with exponential radial density profiles.We found that the unstable modes are global (the wavelengths are larger than the disk scale lengths) and that the instability saturation level is much higher than that following from a linear criterion. The instability saturation time scales are of the order of one billion years or more. For this reason, the bending instability can play an important role in the secular heating of a stellar disk in the $z$ direction. In an extensive series of numerical $N$-body simulations with a high spatial resolution, we were able to scan in detail the space of key parameters (the initial disk thickness $z_0$, the Toomre parameter $Q$, and the ratio of dark halo mass to disk mass $M_{\rm h} / M_{\rm d}$). We revealed three distinct mechanisms of disk heating in the $z$ direction: bending instability of the entire disk, bending instability of the bar, and heating on vertical inhomogeneities in the distribution of stellar matter.