Young massive clusters (YMCs) are the most intense regions of star formation in galaxies. Formulating a model for YMC formation whilst at the same time meeting the constraints from observations is highly challenging however. We show that forming YMCs requires clouds with densities $\gtrsim$ 100 cm$^{-3}$ to collide with high velocities ($\gtrsim$ 20 km s$^{-1}$). We present the first simulations which, starting from moderate cloud densities of $\sim100$ cm$^{-3}$, are able to convert a large amount of mass into stars over a time period of around 1 Myr, to produce dense massive clusters similar to those observed. Such conditions are commonplace in more extreme environments, where YMCs are common, but atypical for our Galaxy, where YMCs are rare.

We investigate the alignment of haloes with the filaments of the cosmic web using an unprecedently large sample of dark matter haloes taken from the P-Millennium $\Lambda$CDM cosmological N-body simulation. We use the state-of-the-art NEXUS morphological formalism which, due to its multiscale nature, simultaneously identifies structures at all scales. We find strong and highly significant alignments, with both the major axis of haloes and their peculiar velocity tending to orient along the filament. However, the spin - filament alignment displays a more complex trend changing from preferentially parallel at low masses to preferentially perpendicular at high masses. This "spin flip" occurs at an average mass of $5\times10^{11}~h^{-1}M_\odot$. This mass increases with increasing filament diameter, varying by more than an order of magnitude between the thinnest and thickest filament samples. We also find that the inner parts of haloes have a spin flip mass that is several times smaller than that of the halo as a whole. These results confirm that recent accretion is responsible for the complex behaviour of the halo spin - filament alignment. Low-mass haloes mainly accrete mass along directions perpendicular to their host filament and thus their spins tend to be oriented along the filaments. In contrast, high-mass haloes mainly accrete along their host filaments and have their spins preferentially perpendicular to them. Furthermore, haloes located in thinner filaments are more likely to accrete along their host filaments than haloes of the same mass located in thicker filaments.

We present simulations of the dynamical evolution of young embedded star clusters. Our initial conditions are directly derived from X-ray, infrared, and radio observations of local systems, and our models evolve both gas and stars simultaneously. Our regions begin with both clustered and extended distributions of stars, and a gas distribution which can include a filamentary structure in addition to gas surrounding the stellar subclusters. We find that the regions become spherical, monolithic, and smooth quite quickly, and that the dynamical evolution is dominated by the gravitational interactions between the stars. In the absence of stellar feedback, the gas moves gently out of the centre of our regions but does not have a significant impact on the motions of the stars at the earliest stages of cluster formation. Our models at later times are consistent with observations of similar regions in the local neighbourhood. We conclude that the evolution of young proto-star clusters is relatively insensitive to reasonable choices of initial conditions. Models with more realism, such as an initial population of binary and multiple stars and ongoing star formation, are the next step needed to confirm these findings.

Context. J1407 (1SWASP J140747.93-394542.6 in full) is a young star in the Scorpius-Centaurus OB association that underwent a series of complex eclipses over 56 days in 2007. To explain these, it was hypothesised that a secondary substellar companion, J1407b, has a giant ring system filling a large fraction of the Hill sphere, causing the eclipses. Observations have not successfully detected J1407b, but do rule out circular orbits for the companion around the primary star. Aims. We test to what degree the ring model of J1407b could survive in an eccentric orbit required to fit the observations. Methods. We run N-body simulations under the AMUSE framework to test the stability of Hill radius-filling systems where the companion is on an eccentric orbit. Results. We strongly rule out prograde ring systems and find that a secondary of 60 to 100MJup with an 11 year orbital period and retrograde orbiting material can survive for at least $10^4$ orbits and produce eclipses with similar durations as the observed one.

The large-scale structure of the universe is comprised of virialized blob-like clusters, linear filaments, sheet-like walls and huge near empty three-dimensional voids. Characterizing the large scale universe is essential to our understanding of the formation and evolution of galaxies. The density range of clusters, walls and voids are relatively well separated, when compared to filaments, which span a relatively larger range. The large scale filamentary network thus forms an intricate part of the cosmic web. In this paper, we describe Felix, a topology based framework for visual exploration of filaments in the cosmic web. The filamentary structure is represented by the ascending manifold geometry of the 2-saddles in the Morse-Smale complex of the density field. We generate a hierarchy of Morse-Smale complexes and query for filaments based on the density ranges at the end points of the filaments. The query is processed efficiently over the entire hierarchical Morse-Smale complex, allowing for interactive visualization. We apply Felix to computer simulations based on the heuristic Voronoi kinematic model and the standard $\Lambda$CDM cosmology, and demonstrate its usefulness through two case studies. First, we extract cosmic filaments within and across cluster like regions in Voronoi kinematic simulation datasets. We demonstrate that we produce similar results to existing structure finders. Filaments that form the spine of the cosmic web, which exist in high density regions in the current epoch, are isolated using Felix. Also, filaments present in void-like regions are isolated and visualized. These filamentary structures are often over shadowed by higher density range filaments and are not easily characterizable and extractable using other filament extraction methodologies.

We study the formation and evolution of the cosmic web, using the high-resolution CosmoGrid \$\Lambda\$CDM simulation. In particular, we investigate the evolution of the large-scale structure around void halo groups, and compare this to observations of the VGS-31 galaxy group, which consists of three interacting galaxies inside a large void. The structure around such haloes shows a great deal of tenuous structure, with most of such systems being embedded in intra-void filaments and walls. We use the Nexus+ algorithm to detect walls and filaments in CosmoGrid, and find them to be present and detectable at every scale. The void regions embed tenuous walls, which in turn embed tenuous filaments. We hypothesize that the void galaxy group of VGS-31 formed in such an environment.

We present a method to couple N-body star cluster simulations to a cosmological tidal field, using the Astrophysical Multipurpose Software Environment. We apply this method to star clusters embedded in the CosmoGrid dark matter-only LambdaCDM simulation. Our star clusters are born at z = 10 (corresponding to an age of the Universe of about 500 Myr) by selecting a dark matter particle and initializing a star cluster with 32,000 stars on its location. We then follow the dynamical evolution of the star cluster within the cosmological environment. We compare the evolution of star clusters in two Milky-Way size haloes with a different accretion history. The mass loss of the star clusters is continuous irrespective of the tidal history of the host halo, but major merger events tend to increase the rate of mass loss. From the selected two dark matter haloes, the halo that experienced the larger number of mergers tends to drive a smaller mass loss rate from the embedded star clusters, even though the final masses of both haloes are similar. We identify two families of star clusters: native clusters, which become part of the main halo before its final major merger event, and the immigrant clusters, which are accreted upon or after this event; native clusters tend to evaporate more quickly than immigrant clusters. Accounting for the evolution of the dark matter halo causes immigrant star clusters to retain more mass than when the z=0 tidal field is taken as a static potential. The reason for this is the weaker tidal field experienced by immigrant star clusters before merging with the larger dark matter halo.

We study the formation and evolution of filamentary configurations of dark matter haloes in voids. Our investigation uses the high-resolution LambdaCDM simulation CosmoGrid to look for void systems resembling the VGS_31 elongated system of three interacting galaxies that was recently discovered by the Void Galaxy Survey (VGS) inside a large void in the SDSS galaxy redshift survey. HI data revealed these galaxies to be embedded in a common elongated envelope, possibly embedded in intravoid filament. In the CosmoGrid simulation we look for systems similar to VGS_31 in mass, size and environment. We find a total of eight such systems. For these systems, we study the distribution of neighbour haloes, the assembly and evolution of the main haloes and the dynamical evolution of the haloes, as well as the evolution of the large-scale structure in which the systems are embedded. The spatial distribution of the haloes follows that of the dark matter environment. We find that VGS_31-like systems have a large variation in formation time, having formed between 10 Gyr ago and the present epoch. However, the environments in which the systems are embedded evolved resemble each other substantially. Each of the VGS_31-like systems is embedded in an intra-void wall, that no later than z = 0.5 became the only prominent feature in its environment. While part of the void walls retain a rather featureless character, we find that around half of them are marked by a pronounced and rapidly evolving substructure. Five haloes find themselves in a tenuous filament of a few Mpc/h long inside the intra-void wall. Finally, we compare the results to observed data from VGS_31. Our study implies that the VGS_31 galaxies formed in the same (proto)filament, and did not meet just recently. The diversity amongst the simulated halo systems indicates that VGS_31 may not be typical for groups of galaxies in voids.

We present the results of the "Cosmogrid" cosmological N-body simulation suites based on the concordance LCDM model. The Cosmogrid simulation was performed in a 30Mpc box with 2048^3 particles. The mass of each particle is 1.28x10^5 Msun, which is sufficient to resolve ultra-faint dwarfs. We found that the halo mass function shows good agreement with the Sheth & Tormen fitting function down to ~10^7 Msun. We have analyzed the spherically averaged density profiles of the three most massive halos which are of galaxy group size and contain at least 170 million particles. The slopes of these density profiles become shallower than -1 at the inner most radius. We also find a clear correlation of halo concentration with mass. The mass dependence of the concentration parameter cannot be expressed by a single power law, however a simple model based on the Press-Schechter theory proposed by Navarro et al. gives reasonable agreement with this dependence. The spin parameter does not show a correlation with the halo mass. The probability distribution functions for both concentration and spin are well fitted by the log-normal distribution for halos with the masses larger than ~10^8 Msun. The subhalo abundance depends on the halo mass. Galaxy-sized halos have 50% more subhalos than ~10^{11} Msun halos have.

We present results from our cosmological N-body simulation which consisted of 2048x2048x2048 particles and ran distributed across three supercomputers throughout Europe. The run, which was performed as the concluding phase of the Gravitational Billion Body Problem DEISA project, integrated a 30 Mpc box of dark matter using an optimized Tree/Particle Mesh N-body integrator. We ran the simulation up to the present day (z=0), and obtained an efficiency of about 0.93 over 2048 cores compared to a single supercomputer run. In addition, we share our experiences on using multiple supercomputers for high performance computing and provide several recommendations for future projects.