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.     

Magnetic zones of Mars: Deformation‐controlled origin of magnetic anomalies

Abstract— Intense magnetic anomalies over Martian surface suggest preservation of large volumes of very old crust (>3 Gyr) that formed in the presence of a global magnetic field. The global distribution of the magnetic intensities observed above the Martian crust suggests a division into three zones. Zone 1 is where the magnetic signature is negligible or of relatively low intensity at Mars Global Surveyor (MGS) satellite mapping altitude (400 km). Zone 2 is the region of intermediate crustal magnetic amplitudes and zone 3 is where the highest magnetic intensities are measured. Crater demagnetization near zone 3 reveals the presence of rocks with both high magnetic intensity and coercivity. Magnetic analyses of terrestrial rocks show that compositional banding in orogenic zones significantly enhances both magnetic coercivity and thermal remanent magnetization (TRM) efficiency. Such enhancement offers a novel explanation for the anomalously large intensities inferred of magnetic sources on Mars. We propose that both large magnetic coercivity and intensity near the South Pole is indicative of the presence of a large degree of deformation. Associated compositional zoning creates conditions for large scale magnetic anisotropy allowing magnetic minerals to acquire magnetization more efficiently, thereby causing the distinct magnetic signatures in zone 3, expressed by intense magnetic anomalies. We use a simple model to verify the magnetic enhancement. We hypothesize that magnetically enhanced zone would reside over the down welling plume at the time of magnetization acquisition.     

A model for binary-binary close encounters and collisions from a dynamical point of view

The goal of this paper is to provide a model for binary-binary interactions in star clusters, which is based on simultaneous binary collision of a special case of the one-dimensional 4-body problem where four masses move symmetrically about the center of mass. From the theoretical point of view, the singularity due to binary collisions between point masses can be handled by means of regularization theory. Our main tool is a change of coordinates due to McGehee by which we blow-up the singular set associated to total collision and replace it with an invariant manifold which includes binary and simultaneous binary collisions, and then gain a complete picture of the local behavior of the solutions near to total collision via the homothetic orbit.     

Comparison of four meteorite penetration funnels in the Campo del Cielo crater field,Argentina

Abstract– More craters may be discovered in the future, but as it is currently known, the Campo del Cielo crater field is 18 km long by 4 km at its widest point. Such a distribution of craters suggests that the parent meteoroid entered and traversed the atmosphere at a very low angle relative to horizontal. The crater field contains at least 20 small craters produced by the larger fragments of the parent meteoroid. Four of these are explosion analog craters and the rest are penetration funnels. During four field seasons, we have constructed topographic and magnetic maps of four of the penetration funnels as found, and then dug trenches across them to learn their original structures and recover meteorites preserved within them. Structures of these penetration funnels indicate very low angles of impact, i.e., 9–16° relative to horizontal. This supports the idea that the parent meteoroid traversed the atmosphere at a low angle. Data given here for the four penetration funnels include projectile masses, lengths, widths, depths, and estimates of impact angles and azimuths. One of the penetration funnels described here (No. 6) can almost be classified as an explosion analog crater.     

A Snapshot of the Sun Near Solar Minimum: The Whole Heliosphere Interval

We present an overview of the data and models collected for the Whole Heliosphere Interval, an international campaign to study the three-dimensional solar?Cheliospheric?Cplanetary connected system near solar minimum. The data and models correspond to solar Carrington Rotation 2068 (20 March??C?16 April 2008) extending from below the solar photosphere, through interplanetary space, and down to Earth??s mesosphere. Nearly 200 people participated in aspects of WHI studies, analyzing and interpreting data from nearly 100 instruments and models in order to elucidate the physics of fundamental heliophysical processes. The solar and inner heliospheric data showed structure consistent with the declining phase of the solar cycle. A closely spaced cluster of low-latitude active regions was responsible for an increased level of magnetic activity, while a highly warped current sheet dominated heliospheric structure. The geospace data revealed an unusually high level of activity, driven primarily by the periodic impingement of high-speed streams. The WHI studies traced the solar activity and structure into the heliosphere and geospace, and provided new insight into the nature of the interconnected heliophysical system near solar minimum.     

An algorithm to resolve γ-rays from charged cosmic rays with DAMPE

The DArk Matter Particle Explorer(DAMPE),also known as Wukong in China,which was launched on 2015 December 17,is a new high energy cosmic ray and γ-ray satellite-borne observatory.One of the main scientific goals of DAMPE is to observe Ge V-Te V high energy γ-rays with accurate energy,angular and time resolution,to indirectly search for dark matter particles and for the study of high energy astrophysics. Due to the comparatively higher fluxes of charged cosmic rays with respect to γ-rays,it is challenging to identify γ-rays with sufficiently high efficiency,minimizing the amount of charged cosmic ray contamination. In this work we present a method to identify γ-rays in DAMPE data based on Monte Carlo simulations,using the powerful electromagnetic/hadronic shower discrimination provided by the calorimeter and the veto detection of charged particles provided by the plastic scintillation detector. Monte Carlo simulations show that after this selection the number of electrons and protons that contaminate the selected γ-ray events at～10 Ge V amounts to less than 1% of the selected sample.Finally,we use flight data to verify the effectiveness of the method by highlighting known γ-ray sources in the sky and by reconstructing preliminary light curves of the Geminga pulsar.     

The ages of pre-main-sequence stars

The position of pre-main-sequence or protostars in the Hertzsprung–Russell diagram is often used to determine their mass and age by comparison with pre-main-sequence evolution tracks. On the assumption that the stellar models are accurate, we demonstrate that, if the metallicity is known, the mass obtained is a good estimate. However, the age determination can be very misleading, because it is significantly (generally different by a factor of 2 to 5) dependent on the accretion rate and, for ages less than about 10⁶ yr, the initial state of the star. We present a number of accreting protostellar tracks that can be used to determine age if the initial conditions can be determined and the underlying accretion rate has been constant in the past. Because of the balance established between the Kelvin–Helmholtz, contraction time-scale and the accretion time-scale, a pre-main-sequence star remembers its accretion history. Knowledge of the current accretion rate, together with an HR-diagram position, gives information about the rate of accretion in the past, but does not necessarily improve any age estimate. We do not claim that ages obtained by comparison with these particular accreting tracks are likely to be any more reliable than those from comparisons with non-accreting tracks. Instead, we stress the unreliability of any such comparisons, and use the disparities between various tracks to estimate the likely errors in age and mass estimates. We also show how a set of coeval accreting objects do not appear coeval when compared with non-accreting tracks. Instead, accreting pre-main-sequence stars of around a solar mass are likely to appear older than those of either smaller or larger mass.     

Mechanism of generation of the emission bands in the dynamic spectrum of the Crab pulsar

We show that the proportionately spaced emission bands in the dynamic spectrum of the Crab pulsar fit the oscillations of the square of a Bessel function whose argument exceeds its order. This function has already been encountered in the analysis of the emission from a polarization current with a superluminal distribution pattern: a current whose distribution pattern rotates (with an angular frequency ω) and oscillates (with a frequency  Ω > ω  differing from an integral multiple of ω) at the same time. Using the results of our earlier analysis, we find that the dependence on frequency of the spacing and width of the observed emission bands can be quantitatively accounted for by an appropriate choice of the value of the single free parameter  Ω/ω  . In addition, the value of this parameter, thus implied by Hankins & Eilek's data, places the last peak in the amplitude of the oscillating Bessel function in question at a frequency  (∼Ω³/ω²)  that agrees with the position of the observed ultraviolet peak in the spectrum of the Crab pulsar. We also show how the suppression of the emission bands by the interference of the contributions from differing polarizations can account for the differences in the time and frequency signatures of the interpulse and the main pulse in the Crab pulsar. Finally, we put the emission bands in the context of the observed continuum spectrum of the Crab pulsar by fitting this broad-band spectrum (over 16 orders of magnitude of frequency) with that generated by an electric current with a superluminally rotating distribution pattern.     

The luminosity oscillation of X-ray bursters and recurrent novae

Thermo-mechanical oscillations of a radiating spherically-symmetric shell containing a fermion gas (i.e., oscillation where only mechanical, thermodynamical, and eventually radiation phenomena are taken into account) are studied. With a reasonable choice of the three relevant parameters all other observational data are similar to the ones of X-ray bursters and recurrent novae.Therefore, thermo-mechanical oscillations could play an important role in the oscillation spectrum of these astronomical objects.