A multiphysics and multiscale software environment for modeling astrophysical systems

We present MUSE, a software framework for combining existing computational tools for different astrophysical domains into a single multiphysics, multiscale application. MUSE facilitates the coupling of existing codes written in different languages by providing inter-language tools and by specifying an interface between each module and the framework that represents a balance between generality and computational efficiency. This approach allows scientists to use combinations of codes to solve highly coupled problems without the need to write new codes for other domains or significantly alter their existing codes. MUSE currently incorporates the domains of stellar dynamics, stellar evolution and stellar hydrodynamics for studying generalized stellar systems. We have now reached a “Noah’s Ark” milestone, with (at least) two available numerical solvers for each domain. MUSE can treat multiscale and multiphysics systems in which the time- and size-scales are well separated, like simulating the evolution of planetary systems, small stellar associations, dense stellar clusters, galaxies and galactic nuclei. In this paper we describe three examples calculated using MUSE: the merger of two galaxies, the merger of two evolving stars, and a hybrid N-body simulation. In addition, we demonstrate an implementation of MUSE on a distributed computer which may also include special-purpose hardware, such as GRAPEs or GPUs, to accelerate computations. The current MUSE code base is publicly available as open source at http://muse.li.     

A Full-Sun Magnetic Index from Helioseismology Inferences

Solar magnetic indices are used to model the solar irradiance and ultimately to forecast it. However, the observation of such indices is generally limited to the Earth-facing hemisphere of the Sun. Seismic maps of the far side of the Sun have proven their capability to locate and track medium–large active regions at the non-visible hemisphere. We present here the possibility of using the average signal from these seismic far-side maps, combined with similarly calculated near-side maps, as a proxy to the full-Sun magnetic activity.     

Photometric Study of the W Uma-Type Eclipsing Binary System GSC 4832.400

We present a photometric analysis of the eclipsing binary systemGSC4832.400 based on 732 photometric observations obtained on sixteen clear nights during 1999–2000 in the V-Johnson and yb-Strömgren filters. The observed times of the minima yielded a linear ephemeris showing that the system has a constant period of 0.313150 ± 0.000002days. The light curve corresponds to a circular orbit for the system and clearly shows that GSC4832.400 has the defining characteristics of a W UMa-type eclipsing binary. It also shows a small O'Connell effect, which could be due to the presence of a hot spot in the primary star. A preliminary photometric solution based on the assumption of an over contact systemsatis factorily reproduced the observed light curve.     

Low velocity shock-cloud encounters I.

Shocks propagating in the interstellar medium (ISM) play an important role in the life of molecular clouds. Through a theoretical study of interaction between clouds and shocks we can understand, for example, the density distribution of observed molecular clouds and the first steps of star formation. The only way to study of interaction in detail is via a numerical hydrodynamical simulation. In this paper we present the first results of a hydrocode which is able to follow the processes after the collision between the cloud and shock front.Our main theoretical result is that the chemical processes (e.g. H₂ dissociation) can affect the dynamical processes significantly. Global parameters of the cloud are calculated for the comparision of the simulation and the observations.     

The search for hot spots associated with the Cygnus X-3 relativistic jet

We report and shortly discuss here the observational work carried out in order to test the possibility that two previously detected radio sources, in the vicinity of the well known microquasar Cygnus X-3, could be hot spot tracers of interaction between its relativistic jet and the interstellar medium (ISM). The motivation behind this search is in part justified considering recent theoretical models of high energy γ-ray emission which strongly rely on the interaction sites of galactic relativistic jets with nearby ISM clouds. The results presented in this paper include an improved radio exploration of the several arc-minute field around Cygnus X-3 using the Very Large Array (VLA), as well as deep near infrared (NIR) imaging with the Calar Alto 3.5 m telescope. We anticipate here that our observations do not appear to support the initial hot spot hypothesis. Instead, the resulting images suggest that the two radio sources, originally believed to be hot spot candidates, are most likely background or foreground objects.     

Dynamics of a satellite orbiting a planet with an inhomogeneous gravitational field

We study the dynamics of a satellite (artificial or natural) orbiting an Earth-like planet at low altitude from an analytical point of view. The perturbation considered takes into account the gravity attraction of the planet and in particular it is caused by its inhomogeneous potential. We begin by truncating the equations of motion at second order, that is, incorporating the zonal and the tesseral harmonics up to order two. The system is formulated as an autonomous Hamiltonian and has three degrees of freedom. After three successive Lie transformations, the system is normalised with respect to two angular co-ordinates up to order five in a suitable small parameter given by the quotient between the angular velocity of the planet and the mean motion of the satellite. Our treatment is free of power expansions of the eccentricity and of truncated Fourier series in the anomalies. Once these transformations are performed, the truncated Hamiltonian defines a system of one degree of freedom which is rewritten as a function of two variables which generate a phase space which takes into account all of the symmetries of the problem. Next an analysis of the system is achieved obtaining up to six relative equilibria and three types of bifurcations. The connection with the original system is established concluding the existence of various families of invariant 3-tori of it, as well as quasiperiodic and periodic trajectories. This is achieved by using KAM theory techniques.     

Observational Evidence for Rising Penumbral Flux Tubes?

On 13 May 2000 parts of a penumbra were observed in an active region NOAA 8990 with the La Palma Stokes Polarimeter attached to the Swedish Vacuum Solar Telescope. The stratification over the solar atmosphere of different physical parameters is retrieved from these data by using the Stokes inversion based on response functions. The results confirm the previous findings of the penumbral structure. In general, the magnetic field becomes weaker and more horizontal with increasing distance from the umbra and the line-of-sight velocities are increasing towards the outer boundary of the penumbra. The results also suggest the existence of the unresolved fine structure of the penumbra. The stratifications of the temperature and of the magnetic field strength indicate the presence of rising flux tubes, which were predicted theoretically by Schlichenmaier, Jahn and Schmidt (1998, Astron. Astrophys. 337, 897).     

Fourier Power Spectra of Solar Noise Storms

We analyzed three noise storms recorded on 200?–?400 MHz Trieste Callisto radio spectra on 2 July 2012, 8 July 2012, and 16 July 2012 by the Fourier method. We divided intervals of the noise storms into five-minute intervals, and in these intervals we computed the mean Fourier spectra as a function of the wave numbers in the frequency and height-scale spaces. We found that these Fourier spectra, where the spectrum from the quiet-activity interval was subtracted, are power-law spectra. The mean power-law index of these spectra in the range \(\ln(k_{z}) = [1.8, 2.9]\) (where \(k_{z}\) is the wave number in the height-scale space) is \(-1.7\pm0.14\), \(-1.6\pm0.14\), and \(-1.5 \pm0.12\) for the 2 July 2012, the 8 July 2012, and the 16 July 2012 noise storms, respectively. It appears that as the number of Type-I bursts in the studied interval increases, the power-law index becomes closer to \(-5/3\); this is known as the Kolmogorov spectral index. The power-law index of the noise storms is very similar to that of the narrowband dm-spikes found in our previous studies. Furthermore, we found a break in the power spectra at \(\ln(k_{z}) \approx2.9\), and the mean power-law index values above this break are \(-2.9\pm0.46\), \(-3.1\pm0.65\), and \(-3.4\pm0.98\), respectively.     

Fitting the integrated spectral energy distributions of galaxies

Fitting the spectral energy distributions (SEDs) of galaxies is an almost universally used technique that has matured significantly in the last decade. Model predictions and fitting procedures have improved significantly over this time, attempting to keep up with the vastly increased volume and quality of available data. We review here the field of SED fitting, describing the modelling of ultraviolet to infrared galaxy SEDs, the creation of multiwavelength data sets, and the methods used to fit model SEDs to observed galaxy data sets. We touch upon the achievements and challenges in the major ingredients of SED fitting, with a special emphasis on describing the interplay between the quality of the available data, the quality of the available models, and the best fitting technique to use in order to obtain a realistic measurement as well as realistic uncertainties. We conclude that SED fitting can be used effectively to derive a range of physical properties of galaxies, such as redshift, stellar masses, star formation rates, dust masses, and metallicities, with care taken not to over-interpret the available data. Yet there still exist many issues such as estimating the age of the oldest stars in a galaxy, finer details of dust properties and dust-star geometry, and the influences of poorly understood, luminous stellar types and phases. The challenge for the coming years will be to improve both the models and the observational data sets to resolve these uncertainties. The present review will be made available on an interactive, moderated web page (), where the community can access and change the text. The intention is to expand the text and keep it up to date over the coming years.     

A minimal dynamical model for tidal synchronization and orbit circularization

We study tidal synchronization and orbit circularization in a minimal model that takes into account only the essential ingredients of tidal deformation and dissipation in the secondary body. In previous work we introduced the model (Escribano et al. in Phys. Rev. E, 78:036216, 2008); here we investigate in depth the complex dynamics that can arise from this simplest model of tidal synchronization and orbit circularization. We model an extended secondary body of mass m by two point masses of mass m/2 connected with a damped spring. This composite body moves in the gravitational field of a primary of mass M ≫ m located at the origin. In this simplest case oscillation and rotation of the secondary are assumed to take place in the plane of the Keplerian orbit. The gravitational interactions of both point masses with the primary are taken into account, but that between the point masses is neglected. We perform a Taylor expansion on the exact equations of motion to isolate and identify the different effects of tidal interactions. We compare both sets of equations and study the applicability of the approximations, in the presence of chaos. We introduce the resonance function as a resource to identify resonant states. The approximate equations of motion can account for both synchronization into the 1:1 spin-orbit resonance and the circularization of the orbit as the only true asymptotic attractors, together with the existence of relatively long-lived metastable orbits with the secondary in p:q (p and q being co-prime integers) synchronous rotation.