Hierarchical code coupling strategies make it possible to combine the results of individual numerical solvers into a self-consistent symplectic solution. We explore the possibility of allowing such a coupling strategy to be non-intrusive. In that case, the underlying numerical implementation is not affected by the coupling itself, but its functionality is carried over in the interface. This method is efficient for solving the equations of motion for a self-gravitating system over a wide range of scales. We adopt a dedicated integrator for solving each particular part of the problem and combine the results to a self-consistent solution. In particular, we explore the possibilities of combining the evolution of one or more microscopic systems that are embedded in a macroscopic system. The here presented generalizations of Bridge include higher-order coupling strategies (from the classic 2nd order up to 10th-order), but we also demonstrate how multiple bridges can be nested and how additional processes can be introduced at the bridge time-step to enrich the physics, for example by incorporating dissipative processes. Such augmentation allows for including additional processes in a classic Newtonian N-body integrator without alterations to the underlying code. These additional processes include for example the Yarkovsky effect, dynamical friction or relativistic dynamics. Some of these processes operate on all particles whereas others apply only to a subset. The presented method is non-intrusive in the sense that the underlying methods remain operational without changes to the code (apart from adding the get- and set-functions to enable the bridge operator). As a result, the fundamental integrators continue to operate with their internal time step and preserve their local optimizations and parallelism. ... abridged ...

We perform a simulation using the Astrophysical Multipurpose Software Environment of the Orion Trapezium star cluster in which the evolution of the stars and the dynamics of planetary systems are taken into account. The initial conditions from earlier simulations were selected in which the size and mass distributions of the observed circumstellar disks in this cluster are satisfactorily reproduced. Four, five, or size planets per star were introduced in orbit around the 500 solar-like stars with a maximum orbital separation of 400au. Our study focuses on the production of free-floating planets. A total of 357 become unbound from a total of 2522 planets in the initial conditions of the simulation. Of these, 281 leave the cluster within the crossing timescale of the star cluster; the others remain bound to the cluster as free-floating intra-cluster planets. Five of these free-floating intra-cluster planets are captured at a later time by another star. The two main mechanisms by which planets are lost from their host star, ejection upon a strong encounter with another star or internal planetary scattering, drive the evaporation independent of planet mass of orbital separation at birth. The effect of small perturbations due to slow changes in the cluster potential are important for the evolution of planetary systems. In addition, the probability of a star to lose a planet is independent of the planet mass and independent of its initial orbital separation. As a consequence, the mass distribution of free-floating planets indistinguishable from the mass distribution of planets bound to their host star.

Gaia Data Release 2 (Gaia DR2) contains results for 1693 million sources in the magnitude range 3 to 21 based on observations collected by the European Space Agency Gaia satellite during the first 22 months of its operational phase. We describe the input data, models, and processing used for the astrometric content of Gaia DR2, and the validation of these results performed within the astrometry task. Some 320 billion centroid positions from the pre-processed astrometric CCD observations were used to estimate the five astrometric parameters (positions, parallaxes, and proper motions) for 1332 million sources, and approximate positions at the reference epoch J2015.5 for an additional 361 million mostly faint sources. Special validation solutions were used to characterise the random and systematic errors in parallax and proper motion. For the sources with five-parameter astrometric solutions, the median uncertainty in parallax and position at the reference epoch J2015.5 is about 0.04 mas for bright (G<14 mag) sources, 0.1 mas at G=17 mag, and 0.7 mas at G=20 mag. In the proper motion components the corresponding uncertainties are 0.05, 0.2, and 1.2 mas/yr, respectively. The optical reference frame defined by Gaia DR2 is aligned with ICRS and is non-rotating with respect to the quasars to within 0.15 mas/yr. From the quasars and validation solutions we estimate that systematics in the parallaxes depending on position, magnitude, and colour are generally below 0.1 mas, but the parallaxes are on the whole too small by about 0.03 mas. Significant spatial correlations of up to 0.04 mas in parallax and 0.07 mas/yr in proper motion are seen on small (<1 deg) and intermediate (20 deg) angular scales. Important statistics and information for the users of the Gaia DR2 astrometry are given in the appendices.

If the Sun was born in a relatively compact open cluster, it is quite likely that a massive (10MSun) star was nearby when it exploded in a supernova. The repercussions of a supernova can be rather profound, and the current Solar System may still bear the memory of this traumatic event. The truncation of the Kuiper belt and the tilt of the ecliptic plane with respect to the Sun's rotation axis could be such signatures. We simulated the effect of a nearby supernova on the young Solar System using the Astronomical Multipurpose Software Environment. Our calculations are realized in two subsequent steps in which we study the effect of the supernova irradiation on the circumstellar disk and the effect of the impact of the nuclear blast-wave which arrives a few decades later. We find that the blastwave of our adopted supernova exploding at a distance of $0.15$--$0.40$\,pc and at an angle of $35^\circ$--$65^\circ$ with respect to the angular-momentum axis of the circumsolar disk would induce a misalignment between the Sun's equator and its disk to $5^\circ.6\pm1^\circ.2$, consistent with the current value. The blast of a supernova truncates the disk at a radius between $42$ and $55$\,au, which is consistent with the current edge of the Kuiper belt. For the most favored parameters, the irradiation by the supernova as well as the blast wave heat the majority of the disk to $\sim 1200$\,K, which is sufficiently hot to melt chondrules in the circumstellar disk. The majority of planetary system may have been affected by a nearby supernova, some of its repercussions, such as truncation and tilting of the disk, may still be visible in their current planetary system's topology. The amount of material from the supernova blast wave that is accreted by the circumstellar disk is too small by several orders of magnitude to explain the current abundance of the short live radionuclide $^{26}$Al.

Coming high-cadence wide-field optical telescopes will image hundreds of thousands of sources per minute. Besides inspecting the near real-time data streams for transient and variability events, the accumulated data archive is a wealthy laboratory for making complementary scientific discoveries. The goal of this work is to optimise column-oriented database techniques to enable the construction of a full-source and light-curve database for large-scale surveys, that is accessible by the astronomical community. We adopted LOFAR's Transients Pipeline as the baseline and modified it to enable the processing of optical images that have much higher source densities. The pipeline adds new source lists to the archive database, while cross-matching them with the known cataloged sources in order to build a full light-curve archive. We investigated several techniques of indexing and partitioning the largest tables, allowing for faster positional source look-ups in the cross matching algorithms. We monitored all query run times in long-term pipeline runs where we processed a subset of IPHAS data that have image source density peaks over $170,000$ per field of view ($500,000$ deg$^{-2}$). Our analysis demonstrates that horizontal table partitions of declination widths of one-degree control the query run times. Usage of an index strategy where the partitions are densily sorted according to source declination yields another improvement. Most queries run in sublinear time and a few (<20%) run in linear time, because of dependencies on input source-list and result-set size. We observed that for this logical database partitioning schema the limiting cadence the pipeline achieved with processing IPHAS data is 25 seconds.

Due to the high stellar densities in young clusters, planetary systems formed in these environments are likely to have experienced perturbations from encounters with other stars. We carry out direct $N$-body simulations of multi-planet systems in star clusters to study the combined effects of stellar encounters and internal planetary dynamics. These planetary systems eventually become part of the Galactic field population the parental cluster dissolves, which is where most presently-known exoplanets are observed. We show that perturbations induced by stellar encounters lead to distinct signatures in the field planetary systems, most prominently, the excited orbital inclinations and eccentricities. Planetary systems that form within the cluster's half-mass radius are more prone to such perturbations. The orbital elements are most strongly excited in the outermost orbit, but the effect propagates to the entire planetary system through secular evolution. Planet ejections may occur long after a stellar encounter. The surviving planets in these reduced systems tend to have, on average, higher inclinations and larger eccentricities compared to systems that were perturbed less strongly. As soon as the parental star cluster dissolves, external perturbations stop affecting the escaped planetary systems, and further evolution proceeds on a relaxation time scale. The outer regions of these ejected planetary systems tend to relax so slowly that their state carries the memory of their last strong encounter in the star cluster. Regardless of the stellar density, we observe a robust anticorrelation between multiplicity and mean inclination/eccentricity. We speculate that the "Kepler dichotomy" observed in field planetary systems is a natural consequence of their early evolution in the parental cluster.

Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.

Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.

Context. This paper presents an overview of the photometric data that are part of the first Gaia data release. Aims. The principles of the processing and the main characteristics of the Gaia photometric data are presented. Methods. The calibration strategy is outlined briefly and the main properties of the resulting photometry are presented. Results. Relations with other broadband photometric systems are provided. The overall precision for the Gaia photometry is shown to be at the milli-magnitude level and has a clear potential to improve further in future releases.

The Galactic center hosts several hundred early-type stars, about 20% of which lie in the so-called clockwise disk, while the remaining 80% do not belong to any disks. The circumnuclear ring (CNR), a ring of molecular gas that orbits the supermassive black hole (SMBH) with a radius of 1.5 pc, has been claimed to induce precession and Kozai-Lidov oscillations onto the orbits of stars in the innermost parsec. We investigate the perturbations exerted by a gas ring on a nearly-Keplerian stellar disk orbiting a SMBH by means of combined direct N-body and smoothed particle hydrodynamics simulations. We simulate the formation of gas rings through the infall and disruption of a molecular gas cloud, adopting different inclinations between the infalling gas cloud and the stellar disk. We find that a CNR-like ring is not efficient in affecting the stellar disk on a timescale of 3 Myr. In contrast, a gas ring in the innermost 0.5 pc induces precession of the longitude of the ascending node Omega, significantly affecting the stellar disk inclination. Furthermore, the combined effect of two-body relaxation and Omega-precession drives the stellar disk dismembering, displacing the stars from the disk. The impact of precession on the star orbits is stronger when the stellar disk and the inner gas ring are nearly coplanar. We speculate that the warm gas in the inner cavity might have played a major role in the evolution of the clockwise disk.

Astronomical phenomena are governed by processes on all spatial and temporal scales, ranging from days to the age of the Universe (13.8,Gyr) as well as from km size up to the size of the Universe. This enormous range in scales is contrived, but as long as there is a physical connection between the smallest and largest scales it is important to be able to resolve them all, and for the study of many astronomical phenomena this governance is present. Although covering all these scales is a challenge for numerical modelers, the most challenging aspect is the equally broad and complex range in physics, and the way in which these processes propagate through all scales. In our recent effort to cover all scales and all relevant physical processes on these scales we have designed the Astrophysics Multipurpose Software Environment (AMUSE). AMUSE is a Python-based framework with production quality community codes and provides a specialized environment to connect this plethora of solvers to a homogeneous problem solving environment.

We present the open source Astrophysical Multi-purpose Software Environment (AMUSE, www.amusecode.org), 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.

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.

We describe AMUSE, the Astrophysical Multipurpose Software Environment, a programming framework designed to manage multi-scale, multi-physics simulations in a hierarchical, extensible, and internally consistent way. Constructed as a collection of individual modules, AMUSE allows computational tools for different physical domains to be easily combined into a single task. It facilitates the coupling of modules written in different languages by providing inter-language tools and a standard programming interface that represents a balance between generality and computational efficiency. The framework currently incorporates the domains of stellar dynamics, stellar evolution, gas dynamics, and radiative transfer. We present some applications of the framework and outline plans for future development of the package.

We introduce a general-purpose framework for interconnecting scientific simulation programs using a homogeneous, unified software interface. Our framework is intrinsically parallel, and conveniently separates all components in memory. It performs unit conversion between different modules automatically and defines common data structures to communicate across different codes. We use the framework to simulate embedded star clusters. For this purpose we couple solvers for gravitational dynamics, stellar evolution and hydrodynamics to self consistently resolve the dynamical evolution simultaneousy with the internal nuclear evolution of the stars and the hydrodynamic response of the gas. We find, in contrast to earlier studies, that the survival of a young star cluster depends only weakly on the efficiency of star formation. The main reason for this weak dependency is the asymmetric expulsion of the embedding gas from the cluster.