The missing mass problem has not been solved decisively yet. Observations show that if gravity is to be modified, then the MOND theory is its excellent approximation on galactic scales. MOND suggests an adjustments of the laws of physics in the limit of low accelerations. Comparative simulations of interacting galaxies in MOND and Newtonian gravity with dark matter revealed two principal differences: 1) galaxies can have close flybys without ending in mergers in MOND because of weaker dynamical friction, and 2) tidal dwarf galaxies form very easily in MOND. When this is combined with the fact that many interacting galaxies are observed at high redshift, we obtain a new perspective on tidal features: they are often formed by non-merging encounters and tidal disruptions of tidal dwarf galaxies. Here we present the results from our self-consistent MOND $N$-body simulation of a close flyby of two galaxies similar to the Milky Way. It turns out that most types of the structures that are traditionally assigned to galaxy mergers can be formed by non-merging encounters, including tidal arms, bridges, streams, shells, disk warps, thick disks, and most probably also disks of satellites. The success of MOND in explaining the dynamics of galaxies hints us that this way of formation of tidal structures should be considered seriously.

(abridged) The Milky Way (MW) and Andromeda (M31) galaxies possess rotating planes of satellites. Their formation has not been explained satisfactorily yet. It was suggested that the MW and M31 satellites are ancient tidal dwarf galaxies, which could explain their configuration. This suggestion gained support by an analytic backward calculation of the relative MW-M31 orbit in the MOND modified dynamics paradigm by Zhao et al. (2013) implying their close flyby 7-11 Gyr ago. Here we explore the Local Group history in MOND in more detail using a simplified first-ever self-consistent simulation. We note the features induced by the encounter in the simulation and identify their possible real counterparts. The simulation was set to approximately reproduce the observed MW and M31 masses, effective radii, separation, relative velocity and disk inclinations. We used the publicly available adaptive-mesh-refinement code Phantom of RAMSES. In the simulation, matter was transferred from the MW to M31 along a tidal tail. The encounter induced formation of several structures resembling the peculiarities of the Local Group. Most notably: 1) A rotating planar structure formed around M31 from the transferred material. It had a size similar to the observed satellite plane and was oriented edge-on to the simulated MW, just as the real one. 2) The same structure also resembled the tidal features observed around M31 by its size and morphology. 3) A warp in the MW developed with an amplitude and orientation similar to that observed. 4) A cloud of particles formed around the simulated MW, with the extent of the actual MW satellite system. The encounter did not end by merging in a Hubble time. The simulation thus demonstrated that MOND can possibly explain many peculiarities of the Local Group and, moreover, that non-merging galaxy encounters in MOND can produce tidal features in galaxies.

Whether BDs form as stars through gravitational collapse ("star-like") or BDs and some very low-mass stars constitute a separate population which form alongside stars comparable to the population of planets, e.g. through circumstellar disk ("peripheral") fragmentation, is one of the key questions of the star-formation problem. For young stars in Taurus-Auriga the binary fraction is large with little dependence on primary mass above ~0.2Msun, while for BDs it is <10%. We investigate a case in which BDs in Taurus formed dominantly through peripheral fragmentation. The decline of the binary frequency in the transition region between star-like and peripheral formation is modelled. A dynamical population synthesis model is employed in which stellar binary formation is universal. Peripheral objects form separately in circumstellar disks with a distinctive initial mass function (IMF), own orbital parameter distributions for binaries and a low binary fraction. A small amount of dynamical processing of the stellar component is accounted for as appropriate for the low-density Taurus-Auriga embedded clusters. The binary fraction declines strongly between the mass-limits for star-like and peripheral formation. The location of characteristic features and the steepness depend on these mass-limits. Such a trend might be unique to low density regions hosting dynamically unprocessed binary populations. The existence of a strong decline in the binary fraction -- primary mass diagram will become verifiable in future surveys on BD and VLMS binarity in the Taurus-Auriga star forming region. It is a test of the (non-)continuity of star formation along the mass-scale, the separateness of the stellar and BD populations and the dominant formation channel for BDs and BD binaries in regions of low stellar density hosting dynamically unprocessed populations.

According to the standard model of cosmology, galaxies are embedded in dark matter halos which are made of particles beyond the standard model of particle physics, thus extending the mass and the size of the visible baryonic matter by typically two orders of magnitude. The observed gas distribution throughout the nearby M81 group of galaxies shows evidence for past significant galaxy--galaxy interactions but without a merger having occurred. This group is here studied for possible dynamical solutions within the dark-matter standard model. In order to cover a comprehensive set of initial conditions, the inner three core members M81, M82 and NGC3077 are treated as a three-body model based on Navarro-Frenk-White profiles. The possible orbits of these galaxies are examined statistically taking into account dynamical friction. Long living, non-merging initial constellations which allow multiple galaxy-galaxy encounters comprise unbound galaxies only, which are arriving from a far distance and happen to simultaneously encounter each other within the recent 500 Myr. Our results are derived by the employment of two separate and independent statistical methods, namely a Markov chain Monte Carlo method and the genetic algorithm using the SAP system environment. The conclusions reached are confirmed by high-resolution simulations of live self-consistent systems (N-body calculations). Given the observed positions of the three galaxies the solutions found comprise predictions for their proper motions.

Since its publication 1983, Milgromian dynamics (aka MOND) has been very successful in modeling the gravitational potential of galaxies from baryonic matter alone. However, the dynamical modeling has long been an unsolved issue. In particular, the setup of a stable galaxy for Milgromian N-body calculations has been a major challenge. Here, we show a way to set up disc galaxies in MOND for calculations in the PHANTOM OF RAMSES (PoR) code by L\"ughausen (2015) and Teyssier (2002). The method is done by solving the QUMOND Poisson equations based on a baryonic and a phantom dark matter component. The resulting galaxy models are stable after a brief settling period for a large mass and size range. Simulations of single galaxies as well as colliding galaxies are shown.

The separation distribution for M-dwarf binaries in the ASTRALUX survey is narrower and peaking at smaller separations than the distribution for solar-type binaries. This is often interpreted to mean that M-dwarfs constitute a continuous transition from brown dwarfs (BDs) to stars. Here a prediction for the M-dwarf separation distribution is presented, using a dynamical population synthesis (DPS) model in which "star-like" binaries with late-type primaries ($\lesssim1.5 M_{\rm sun}$) follow universal initial distribution functions and are dynamically processed in their birth embedded clusters. A separate "BD-like" population has both its own distribution functions for binaries and initial mass function (IMF), which overlaps in mass with the IMF for stars. Combining these two formation modes results in a peak on top of a wider separation distribution for late M-dwarfs consistent with the late ASTRALUX sample. The DPS separation distribution for early M-dwarfs shows no such peak and is in agreement with the M-dwarfs in Multiples (MinMS) data. We note that the latter survey is potentially in tension with the early ASTRALUX data. Concluding, the ASTRALUX and MinMS data are unable to unambiguously distinguish whether or not BDs are a continuous extension of the stellar IMF. Future observational efforts are needed to fully answer this interesting question. The DPS model predicts that binaries outside the sensitivity range of the ASTRALUX survey remain to be detected. For application to future data, we present a means to observationally measure the overlap of the putative BD-like branch and the stellar branch. We discuss the meaning of universal star formation and distribution functions.

The stellar initial mass function (IMF) is a key property of stellar populations. There is growing evidence that the classical star-formation mechanism by the direct cloud fragmentation process has difficulties reproducing the observed abundance and binary properties of brown dwarfs and very-low-mass stars. In particular, recent analytical derivations of the stellar IMF exhibit a deficit of brown dwarfs compared to observational data. Here we derive the residual mass function of brown dwarfs as an empirical measure of the brown dwarf deficiency in recent star-formation models with respect to observations and show that it is compatible with the substellar part of the Thies-Kroupa IMF and the mass function obtained by numerical simulations. We conclude that the existing models may be further improved by including a substellar correction term that accounts for additional formation channels like disk or filament fragmentation. The term "peripheral fragmentation" is introduced here for such additional formation channels. In addition, we present an updated analytical model of stellar and substellar binarity. The resulting binary fraction and the dynamically evolved companion mass-ratio distribution are in good agreement with observational data on stellar and very-low-mass binaries in the Galactic field, in clusters, and in dynamically unprocessed groups of stars if all stars form as binaries with stellar companions. Cautionary notes are given on the proper analysis of mass functions and the companion mass-ratio distribution and the interpretation of the results. The existence of accretion disks around young brown dwarfs does not imply that these form just like stars in direct fragmentation.

We present the adaptable MUESLI code for investigating dynamics and erosion processes of globular clusters (GCs) in galaxies. MUESLI follows the orbits of individual clusters and applies internal and external dissolution processes to them. Orbit integration is based on the self-consistent field method in combination with a time-transformed leapfrog scheme, allowing us to handle velocity-dependent forces like triaxial dynamical friction. In a first application, the erosion of globular cluster systems (GCSs) in elliptical galaxies is investigated. Observations show that massive ellipticals have rich, radially extended GCSs, while some compact dwarf ellipticals contain no GCs at all. For several representative examples, spanning the full mass scale of observed elliptical galaxies, we quantify the influence of radial anisotropy, galactic density profiles, SMBHs, and dynamical friction on the GC erosion rate. We find that GC number density profiles are centrally flattened in less than a Hubble time, naturally explaining observed cored GC distributions. The erosion rate depends primarily on a galaxy's mass, half-mass radius and radial anisotropy. The fraction of eroded GCs is nearly 100% in compact, M 32 like galaxies and lowest in extended and massive galaxies. Finally, we uncover the existence of a violent tidal disruption dominated phase which is important for the rapid build-up of halo stars.

The current knowledge on the stellar IMF is documented. It appears to become top-heavy when the star-formation rate density surpasses about 0.1Msun/(yr pc^3) on a pc scale and it may become increasingly bottom-heavy with increasing metallicity and in increasingly massive early-type galaxies. It declines quite steeply below about 0.07Msun with brown dwarfs (BDs) and very low mass stars having their own IMF. The most massive star of mass mmax formed in an embedded cluster with stellar mass Mecl correlates strongly with Mecl being a result of gravitation-driven but resource-limited growth and fragmentation induced starvation. There is no convincing evidence whatsoever that massive stars do form in isolation. Various methods of discretising a stellar population are introduced: optimal sampling leads to a mass distribution that perfectly represents the exact form of the desired IMF and the mmax-to-Mecl relation, while random sampling results in statistical variations of the shape of the IMF. The observed mmax-to-Mecl correlation and the small spread of IMF power-law indices together suggest that optimally sampling the IMF may be the more realistic description of star formation than random sampling from a universal IMF with a constant upper mass limit. Composite populations on galaxy scales, which are formed from many pc scale star formation events, need to be described by the integrated galactic IMF. This IGIMF varies systematically from top-light to top-heavy in dependence of galaxy type and star formation rate, with dramatic implications for theories of galaxy formation and evolution.

Recent discoveries of strongly misaligned transiting exoplanets pose a challenge to the established planet formation theory which assumes planetary systems to form and evolve in isolation. However, the fact that the majority of stars actually do form in star clusters raises the question how isolated forming planetary systems really are. Besides radiative and tidal forces the presence of dense gas aggregates in star-forming regions are potential sources for perturbations to protoplanetary discs or systems. Here we show that subsequent capture of gas from large extended accretion envelopes onto a passing star with a typical circumstellar disc can tilt the disc plane to retrograde orientation, naturally explaining the formation of strongly inclined planetary systems. Furthermore, the inner disc regions may become denser, and thus more prone to speedy coagulation and planet formation. Pre-existing planetary systems are compressed by gas inflows leading to a natural occurrence of close-in misaligned hot Jupiters and short-period eccentric planets. The likelihood of such events mainly depends on the gas content of the cluster and is thus expected to be highest in the youngest star clusters.

Most stars are born in clusters and the resulting gravitational interactions between cluster members may significantly affect the evolution of circumstellar discs and therefore the formation of planets and brown dwarfs. Recent findings suggest that tidal perturbations of typical circumstellar discs due to close encounters may inhibit rather than trigger disc fragmentation and so would seem to rule out planet formation by external tidal stimuli. However, the disc models in these calculations were restricted to disc radii of 40 AU and disc masses below 0.1 M_sun. Here we show that even modest encounters can trigger fragmentation around 100 AU in the sorts of massive (~0.5 M_sun), extended (>=100 AU) discs that are observed around young stars. Tidal perturbation alone can do this, no disc-disc collision is required. We also show that very-low-mass binary systems can form through the interaction of objects in the disc. In our computations, otherwise non-fragmenting massive discs, once perturbed, fragment into several objects between about 0.01 and 0.1 M_sun, i.e., over the whole brown dwarf mass range. Typically these orbit on highly eccentric orbits or are even ejected. While probably not suitable for the formation of Jupiter- or Neptune-type planets, our scenario provides a possible formation mechanism for brown dwarfs and very massive planets which, interestingly, leads to a mass distribution consistent with the canonical substellar IMF. As a minor outcome, a possible explanation for the origin of misaligned extrasolar planetary systems is discussed.

The empirical binary properties of brown dwarfs (BDs) differ from those of normal stars suggesting BDs form a separate population. Recent work by Thies & Kroupa revealed a discontinuity of the initial mass function (IMF) in the very-low-mass star regime under the assumption of a low multiplicity of BDs of about 15 per cent. However, previous observations had suggested that the multiplicity of BDs may be significantly higher, up to 45 per cent. This contribution investigates the implication of a high BD multiplicity on the appearance of the IMF for the Orion Nebula Cluster, Taurus-Auriga, IC 348 and the Pleiades. We show that the discontinuity remains pronounced even if the observed MF appears to be continuous, even for a BD binary fraction as high as 60%. We find no evidence for a variation of the BD IMF with star-forming conditions. The BD IMF has a power-law index alpha = +0.3 and about two BDs form per 10 low-mass stars assuming equal-mass pairing of BDs.

The origin of brown dwarfs (BDs) is still an unsolved mystery. While the standard model describes the formation of BDs and stars in a similar way recent data on the multiplicity properties of stars and BDs show them to have different binary distribution functions. Here we show that proper treatment of these uncovers a discontinuity of the multiplicity-corrected mass distribution in the very-low-mass star (VLMS) and BD mass regime. A continuous IMF can be discarded with extremely high confidence. This suggests that VLMSs and BDs on the one hand, and stars on the other, are two correlated but disjoint populations with different dynamical histories. The analysis presented here suggests that about one BD forms per five stars and that the BD-star binary fraction is about 2%-3% among stellar systems.

We present a parameter study of the possibility of tidally triggered disk instability. Using a restricted N-body model which allows for a survey of an extended parameter space, we show that a passing dwarf star with a mass between 0.1 and 1 M_sun can probably induce gravitational instabilities in the pre-planetary solar disk for prograde passages with minimum separations below 80-170 AU for isothermal or adiabatic disks. Inclined and retrograde encounters lead to similar results but require slightly closer passages. Such encounter distances are quite likely in young moderately massive star clusters (Scally & Clarke 2001; Bonnell et al. 2001). The induced gravitational instabilities may lead to enhanced planetesimal formation in the outer regions of the protoplanetary disk, and could therefore be relevant for the existence of Uranus and Neptune, whose formation timescale of about 100 Myr (Wuchterl, Guillot & Lissauer 2000) is inconsistent with the disk lifetimes of about a few Myr according to observational data by Haisch, Lada & Lada (2001). The relatively small gas/solid ratio in Uranus and Neptune can be matched if the perturbing fly-by occurred after early gas depletion of the solar system, i.e. when the solar system was older than about 5 Myr. We also confirm earlier results by Heller (1993) that the observed 7 degree tilt of the solar equatorial plane relative to the ecliptic plane could be the consequence of such a close encounter.