Characterizing multiplicity in the very low mass (VLM) domain is a topic of much current interest and fundamental importance. Here we report on a near-IR AO imaging survey of 31 young brown dwarfs and VLM stars, 28 of which are in Chamaeleon I, using the ESO VLT. Our survey is sensitive enough to detect equal mass binaries down to separations of 0.04-0.07" (6-10 AU at 160 pc) and, typically, companions with mass ratios as low as 0.2 outside of 0.2" (~30 AU). We resolve the suspected 0.16" (~26 AU) binary Cha_Halpha 2 and present two new binaries, Hn 13 and CHXR 15, with separations of 0.13" (~20 AU) and 0.30" (~50 AU); the latter is one of the widest VLM systems discovered to date. We do not find companions around the majority of our targets, giving a binary frequency of 11% [+9,-6], thus confirming the trend for a lower binary frequency with decreasing mass. By combining our work with previous surveys, we arrive at the largest sample of young VLMOs (72) with high angular resolution imaging to date. Its multiplicity fraction is in statistical agreement with that for VLMOs in the field. In addition we note that many field stellar binaries with lower binding energies and/or wider cross sections have survived dynamical evolution and that statistical models suggest tidal disruption by passing stars is unlikely to affect the binary properties of our systems. Thus, we argue that there is no significant evolution of multiplicity with age among brown dwarfs and VLM stars in OB and T associations between a few Myr to several Gyr. Instead, the observations to date suggest that VLM objects are either less likely to be born in fragile multiple systems than solar mass stars or such systems are disrupted very early.

We perform numerical simulations of a disc-planet system using various grid-based and smoothed particle hydrodynamics (SPH) codes. The tests are run for a simple setup where Jupiter and Neptune mass planets on a circular orbit open a gap in a protoplanetary disc during a few hundred orbital periods. We compare the surface density contours, potential vorticity and smoothed radial profiles at several times. The disc mass and gravitational torque time evolution are analyzed with high temporal resolution. There is overall consistency between the codes. The density profiles agree within about 5% for the Eulerian simulations while the SPH results predict the correct shape of the gap although have less resolution in the low density regions and weaker planetary wakes. The disc masses after 200 orbital periods agree within 10%. The spread is larger in the tidal torques acting on the planet which agree within a factor 2 at the end of the simulation. In the Neptune case the dispersion in the torques is greater than for Jupiter, possibly owing to the contribution from the not completely cleared region close to the planet.

It has now been known for over a decade that low-mass stars located in star-forming regions are very frequently members of binary and multiple systems, even more so than main sequence stars in the solar neighborhood. This high multiplicity rate has been interpreted as the consequence of the fragmentation of small molecular cores into a few seed objects that accrete to their final mass from the remaining material and dynamically evolve into stable multiple systems, possibly producing a few ejecta in the process. Analyzing the statistical properties of young multiple systems in a variety of environments therefore represents a powerful approach to place stringent constraints on star formation theories. In this contribution, we first review a number of recent results related to the multiplicity of T Tauri stars. We then present a series of studies focusing on the multiplicity and properties of optically-undetected, heavily embedded protostars. These objects are much younger than the previously studied pre-main sequence stars, and they therefore offer a closer look at the primordial population of multiple systems. In addition to these observational avenues, we present new results of a series of numerical simulations that attempt to reproduce the fragmentation of small molecular cores into multiple systems, and compare these results to the observations.

In this paper we discuss the nature and the possible formation scenarios of the companion of the brown dwarf 2MASS 1207334-393254. We initially discuss the basic physical properties of this object and conclude that, although from its absolute mass ($5M_{\rm Jup}$), it is a planetary object, in terms of its mass ratio $q$ and of its separation $a$ with respect to the primary brown dwarf, it is consistent with the statistical properties of binaries with higher primary mass. We then explore the possible formation mechanism for this object. We show that the standard planet formation mechanism of core accretion is far too slow to form this object within 10 Myr, the observed age of the system. On the other hand, the alternative mechanism of gravitational instability (proposed both in the context of planet and of binary formation) may, in principle, work and form a system with the observed properties.

We have undertaken a series of hydrodynamic + N-body simulations in order to explore the properties of young stars. Our results suggest that the IMF may be sensitive to environment in its substellar region, with more brown dwarfs being formed where clusters are denser or more compact. We find that multiple stars are a natural outcome of collapsing turbulent flows, with a high incidence of N > 2 multiples. We find a positive correlation of multiplicity with primary mass but a companion frequency that decreases with age. Binary brown dwarfs are rarely formed, in conflict with observations. Brown dwarfs as companions are predominantly found orbiting binaries or triples at large separations.

Abridged/ We present numerical results on the properties of young binary and multiple stellar systems. Our analysis is based on a series of SPH + Nbody simulations of the fragmentation of small molecular clouds, that fully resolve the opacity limit for fragmentation. We have produced a statistically significant number of stable multiple systems, with components separations in the range 1-1000 AU. At the end of the hydrodynamical evolution (0.5 Myr) we find that ~60% of stars and brown dwarfs are members of multiples systems, with about a third of these being low mass, weakly bound outliers in wide eccentric orbits. Our results imply that in the stellar regime most stars are in multiples (~80%) and that this fraction is an increasing function of primary mass. After Nbody integration to 10.5 Myr, the percentage of bound objects has dropped to ~40%, as most very low mass stars and brown dwarfs have been released to the field. Brown dwarfs are never found to be very close companions to stars (brown dwarf desert at very small separations), but one case exists of a brown dwarf companion at intermediate separations (10 AU). Our simulations can accommodate the existence of brown dwarf companions at large separations, but only if the primaries of these systems are themselves multiples. We have compared the outcome of our simulations with the properties of real stellar systems as deduced from the IR CM diagram of the Praesepe cluster and from spectroscopic and high-resolution imaging surveys of young clusters and the field.

Abridged: We have undertaken a series of hydrodynamical simulations of multiple star formation in small turbulent molecular clouds. Our goal is to determine the sensitivity of the properties of the resulting stars and brown dwarfs to variations in the initial conditions imposed. In this paper we report on the results obtained by applying two different initial turbulent velocity fields. The slope of the turbulent power-law spectrum alpha is set to -3 in half of the calculations and to -5 in the other half. We find that, whereas the stellar mass function seems to only be weakly dependent on the value of alpha, the sub-stellar mass function turns out to be more sensitive to the initial slope of the velocity field. We argue that, since the role of turbulence is to create substructure from which gravitational instabilities may grow, variations in other initial conditions that also determine the fragmentation process are likely to affect the shape of the sub-stellar mass function as well. The absence of many planetary mass `free-floaters' in our simulations, especially in the mass range 1-10 MJ, suggests that, if these objects are abundant, they are likely to form by similar mechanisms to those thought to operate in quiescent accretion discs, instead of via instabilities in gravitationally unstable discs. Finally, we find that multiple and single stars share comparable kinematical properties, both populations being able to attain velocities in the range 1-10 km/s. From these values we draw the conclusion that only low-mass star-forming regions such as Taurus-Auriga or Ophiuchus, where the escape speed is low, might have suffered some depletion of its single and binary stellar population.

We present results from high-resolution hydrodynamical simulations which explore the effects of small scale clustering in star-forming regions. A large ensemble of small-N clusters with 5 stellar seeds have been modelled and the resulting properties of stars and brown dwarfs statistically derived and compared with observational data. Close dynamical interactions between the protostars and competitive accretion driven by the cloud collapse are shown to produce a distribution of final masses which is bimodal, with most of the mass residing in the binary components. When convolved with a suitable core mass function, the final distribution of masses resembles the observed IMF, both in the stellar and sub-stellar regime. Binaries and single stars are found to constitute two kinematically distinct populations, with about half of the singles attaining velocities > 2 km/s, which might deprive low mass star-forming regions of their lightest members in a few crossing times. The eccentricity distribution of binaries and multiples is found to follow a distribution similar to that of observed long period (uncircularized) binaries. We predict that many close binary stars should have wide brown dwarfs companions. Brown dwarfs and, in general, very low mass stars, would be rare as pure binary companions. The binary fraction should be a decreasing function of primary mass, with low-mass or sub-stellar primaries being scarce. Where such binaries exist, they are either expected to be close enough (semi-major axis ~ 10 AU) to survive strong interactions with more massive binaries or else born in very small molecular cloud cores.

We present results from a survey of dwarf galaxy candidates in the vicinity of strongly interacting galaxies. The goal of the survey was a test of the hypothesis that massive condensations of stars and HI in tidal tails of large interacting galaxies may be a significant source of independent, self-gravitating dwarf galaxies. These so called tidal dwarf galaxies (TDG) can be expected to resemble the blue luminous knots found in tidal tails, but also might appear as redder, evolved systems if formed much before they are being observed. For the present study, a homogeneous subsample of 6 fields was selected from the catalog of extended objects in fields around 15 strongly interacting galaxies by Deeg et al. (1998). Criteria for the subsample were: similar redshifts of the central interacting galaxies, and photometric completeness of the extended objects in V and R. The number density of these TDG candidates was compared with expected background galaxy densities. Within the statistical errors, background galaxies account for most, if not all of the extended objects. There is no evidence for a substantial locally formed dwarf galaxy population. Thus, we conclude that field galaxy-galaxy interactions are likely to result in the formation of only a few long-lived TDG.

Brown dwarfs may have such low masses because they are prematurely ejected from small unstable multiple systems, while the members are still actively building up their masses. We demonstrate that this scenario is consistent with all currently existing observations of brown dwarfs, and propose further observational tests. We review the status of the latest realistic numerical simulations of disintegrating small N clusters, which show that many of the ejected members end up with masses that are substellar, drifting away from their birth region with velocities rarely exceeding 2 km/s.