We present the open source Astrophysical Multi-purpose Software Environment (AMUSE), a component library for performing astrophysical simulations involving different physical domains and scales. It couples existing codes within a Python framework based on a communication layer using MPI. The interfaces are standardized for each domain and their implementation based on MPI guarantees that the whole framework is well-suited for distributed computation. It includes facilities for unit handling and data storage. Currently it includes codes for gravitational dynamics, stellar evolution, hydrodynamics and radiative transfer. Within each domain the interfaces to the codes are as similar as possible. We describe the design and implementation of AMUSE, as well as the main components and community codes currently supported and we discuss the code interactions facilitated by the framework. Additionally, we demonstrate how AMUSE can be used to resolve complex astrophysical problems by presenting example applications.

Context. Young radio-loud active galactic nuclei form an important tool to investigate the evolution of extragalactic radio sources. To study the early phases of expanding radio sources, we have constructed CORALZ, a sample of 25 compact (θ < 2 ) radio sources associated with nearby (z < 0.16) galaxies. Aims. In this paper we determine the morphologies, linear sizes, and place first constraints on the lobe expansion speeds of the sources in the sample. Methods. We observed the radio sources from the CORALZ sample with MERLIN at 1.4 GHz or 1.6 GHz, the EVN at 1.6 GHz, and global VLBI at 1.6 GHz and/or 5.0 GHz. Results. Radio maps, morphological classifications, and linear sizes are presented for all sources in the CORALZ sample. We have determined a first upper limit to the expansion velocity of one of the sources, which is remarkably low compared to the brighter GPS sources at higher redshifts, indicating a relation between radio luminosity and expansion speed, in agreement with analytical models. In addition we present further strong evidence that the spectral turnovers in GPS and CSS sources are caused by synchrotron self-absorption (SSA): the CORALZ sources are significantly offset from the well-known correlation between spectral peak frequency and angular size, but this correlation is recovered after correcting for the flux-density dependence, as predicted by SSA theory.

The evolution of triples has not attracted much attention in the literature, although their evolution can be dramatically different from binaries and single stars. Triples are quite common, and we find that for about 1% of the triples in the Tokovinin catalogue of multiple stellar systems in the solar neighbourhood, the tertiary star will overflow its Roche lobe at some time in its evolution, before any of the inner stars leave the main sequence. For two of these systems, ξ Tau and HD 97131, we simulate in detail this phase of mass transfer, during which stellar evolution, gravitational dynamics and hydrodynamics all play an important role. We have used the Astrophysical Multi-purpose Software Environment (AMUSE) to solve these physical processes in a self-consistent way. The resulting evolution, mass transfer and the effects on the inner as well as on the outer orbit are profound, although it is not trivial to predict the eventual consequence of the phase of mass transfer and the appearance of the resulting system.

We introduce a general-purpose framework for interconnecting scientific simulation programs using a homogeneous, unified interface. Our framework is intrinsically parallel, and conveniently separates all component numerical modules in memory. This strict separation allows automatic unit conversion, distributed execution of modules on different cores within a cluster or grid, and orderly recovery from errors. The framework can be efficiently implemented and incurs an acceptable overhead. In practice, we measure the time spent in the framework to be less than 1% of the wall-clock time. Due to the unified structure of the interface, incorporating multiple modules addressing the same physics in different ways is relatively straightforward. Different modules may be advanced serially or in parallel. Despite initial concerns, we have encountered relatively few problems with this strict separation between modules, and the results of our simulations are consistent with earlier results using more traditional monolithic approaches. This framework provides a platform to combine existing simulation codes or develop new physical solver codes within a rich "ecosystem" of interchangeable modules.

Context. Gigahertz Peaked Spectrum (GPS) radio galaxies are generally thought to be the young counterparts of classical extended radio sources. Statistically complete samples of GPS sources are vital for studying the early evolution of radio-loud AGN and the trigger of their nuclear activity. The "Parkes half-Jansky" sample of GPS radio galaxies is such a sample, representing the southern counterpart of the 1998 Stanghellini sample of bright GPS sources. Aims. As a first step of the investigation of the sample, the host galaxies need to be identified and their redshifts determined. Methods. Deep R-band VLT-FORS1 and ESO 3.6 m EFOSC II images and long slit spectra have been taken for the unidentified sources in the sample. Results. We have identified all twelve previously unknown host galaxies of the radio sources in the sample. Eleven have host galaxies in the range 21.0 < R C < 23.0, while one object, PKS J0210+0419, is identified in the near infrared with a galaxy with K s = 18.3. The redshifts of 21 host galaxies have been determined in the range 0.474 < z < 1.539, bringing the total number of redshifts to 39 (80%). Analysis of the absolute magnitudes of the GPS host galaxies show that at z > 1 they are on average a magnitude fainter than classical 3C radio galaxies, as found in earlier studies. However their restframe UV luminosities indicate that there is an extra light contribution from the AGN, or from a population of young stars.

No abstract

Young radio-loud active galactic nuclei form an important tool to investigate the evolution of extragalactic radio sources. In this paper we present a summary of our recent work on the CORALZ sample of young radio sources at low redshift. We have found strong evidence that the radio spectral turnovers in GPS and CSS sources are caused by synchrotron selfabsorption (SSA): the CORALZ sources follow the well established relation between radio spectral peak frequency and largest angular size, but with significantly smaller sizes at any particular peak frequency, compared to more powerful and more distant GPS/CSS sources, as expected from SSA theory combined with simple self-similar evolution models. Current models that invoke free-free absorption (FFA) to explain the spectral turnovers in GPS and CSS sources can not easily explain the relatively small angular sizes found for the CORALZ sources. We therefore conclude that, although FFA may play a role in some sources, the radio spectral turnovers in GPS and CSS radio sources are generally caused by SSA.

We investigate the effect of different initial virial temperatures, Q, on the dynamics of star clusters. We find that the virial temperature has a strong effect on many aspects of the resulting system, including among others: the fraction of bodies escaping from the system, the depth of the collapse of the system, and the strength of the mass segregation. These differences deem the practice of using "cold" initial conditions no longer a simple choice of convenience. The choice of initial virial temperature must be carefully considered as its impact on the remainder of the simulation can be profound. We discuss the pitfalls and aim to describe the general behavior of the collapse and the resultant system as a function of the virial temperature so that a well reasoned choice of initial virial temperature can be made. We make a correction to the previous theoretical estimate for the minimum radius, R min , of the cluster at the deepest moment of collapse to include a Q dependency, R min ≈ Q+N (−1/3) , where N is the number of particles.We use our numerical results to infer more about the initial conditions of the young cluster R136. Based on our analysis, we find that R136 was likely formed with a rather cool, but not cold, initial virial temperature (Q ≈ 0.13). Using the same analysis method, we examined 15 other young clusters and found the most common initial virial temperature to be between 0.18 and 0.25.