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.     

Radio imaging of the Subaru/XMM–Newton Deep Field – II. The 37 brightest radio sources

We study the 37 brightest radio sources in the Subaru/ XMM–Newton Deep Field. We have spectroscopic redshifts for 24 of 37 objects and photometric redshifts for the remainder, yielding a median redshift z _med for the whole sample of   z _med≃ 1.1  and a median radio luminosity close to the 'Fanaroff–Riley type I/type II (FR I/FR II)' luminosity divide. Using mid-infrared (mid-IR) ( Spitzer MIPS 24 μm) data we expect to trace nuclear accretion activity, even if it is obscured at optical wavelengths, unless the obscuring column is extreme. Our results suggest that above the FR I/FR II radio luminosity break most of the radio sources are associated with objects that have excess mid-IR emission, only some of which are broad-line objects, although there is one clear low-accretion-rate object with an FR I radio structure. For extended steep-spectrum radio sources, the fraction of objects with mid-IR excess drops dramatically below the FR I/FR II luminosity break, although there exists at least one high-accretion-rate 'radio-quiet' QSO. We have therefore shown that the strong link between radio luminosity (or radio structure) and accretion properties, well known at z ∼ 0.1, persists to z ∼ 1. Investigation of mid-IR and blue excesses shows that they are correlated as predicted by a model in which, when significant accretion exists, a torus of dust absorbs ∼30 per cent of the light, and the dust above and below the torus scatters ≳1 per cent of the light.     

XUV Photometer System (XPS): Improved Solar Irradiance Algorithm Using CHIANTI Spectral Models

Solar soft X-ray (XUV) radiation is highly variable on all time scales and strongly affects Earth’s ionosphere and upper atmosphere; consequently, the solar XUV irradiance is important for atmospheric studies and for space weather applications. Although there have been several recent measurements of the solar XUV irradiance, detailed understanding of the solar XUV irradiance, especially its variability during flares, has been hampered by the broad bands measured in the XUV range. In particular, the simple conversion of the XUV photometer signal into irradiance, in which a static solar spectrum is assumed, overestimates the flare variations by more than a factor of two as compared to the atmospheric response to the flares. To address this deficiency in the simple conversion, an improved algorithm using CHIANTI spectral models has been developed to process the XUV Photometer System (XPS) measurements with its broadband photometers. Model spectra representative of quiet Sun, active region, and flares are combined to match the signals from the XPS and produce spectra from 0.1 to 40 nm in 0.1-nm intervals for the XPS Level 4 data product. The two XPS instruments are aboard NASA’s Solar Radiation and Climate Experiment (SORCE) and Thermosphere, Ionosphere, Mesosphere, Energetics, and Dynamics (TIMED) satellites. In addition, the XPS responsivities have been updated for the latest XPS data processing version. The new XPS results are consistent with daily variations from the previous simple conversion technique used for XPS and are also consistent with spectral measurements made at wavelengths longer than 27 nm. Most importantly, the XPS flare variations are reduced by factors of 2 – 4 at wavelengths shorter than 14 nm and are more consistent, for the first time, with atmospheric response to solar flares. Along with the details of the new XPS algorithm, several comparisons to dayglow and photoelectron measurements and model results are also presented to help verify the accuracy of the new XUV irradiance spectra.     

Long-term hydrology and aquatic biogeochemistry data from H. J. Andrews Experimental Forest,Cascade Mountains,Oregon

The H. J. Andrews Experimental Forest (HJA) encompasses the 6400 ha Lookout Creek watershed in western Oregon, USA. Hydrologic, chemistry and precipitation data have been collected, curated, and archived for up to 70 years. The HJA was established in 1948 to study the effects of harvest of old-growth conifer forest and logging-road construction on water quality, quantity and vegetation succession. Over time, research questions have expanded to include terrestrial and aquatic species, communities and ecosystem dynamics. There are nine small experimental watersheds and 10 gaging stations in the HJA, including both reference and experimentally treated watersheds. Gaged watershed areas range from 8.5 to 6242 ha. All gaging stations record stage height, water conductivity, water temperature and above-stream air temperature. At nine of the gage sites, flow-proportional water samples are collected and composited over 3-week intervals for chemical analysis. Analysis of stream and precipitation chemistry began in 1968. Analytes include dissolved and particulate species of nitrogen and phosphorus, dissolved organic carbon, pH, specific conductance, suspended sediment, alkalinity, and major cations and anions. Supporting climate measurements began in the 1950s in association with the first small watershed experiments. Over time, and following the initiation of the Long Term Ecological Research (LTER) grant in 1980, infrastructure expanded to include a set of benchmark and secondary meteorological stations located in clearings spanning the elevation range within the Lookout Creek watershed, as well as a large number of forest understory temperature stations. Extensive metadata on sensor configurations, changes in methods over time, sensor accuracy and precision, and data quality control flags are associated with the HJA data.     

A reconnaissance study of oxygen, hydrogen and strontium isotopes in geochemically diverse lakes, Western Nebraska, USA

Reconnaissance ¹⁸O,, D, and ⁸⁷Sr data for fifteen lakes in the Western Lakes Region of the Sand Hills of Nebraska indicate dynamic hydrologic systems. The rather narrow range of ⁸⁷Sr from lake water (1.1 to 2.1) and groundwater (0.9 to 1.7) indicates that the groundwater is generally unradiogenic. Groundwater residence times and relatively unradiogenic volcanic ash within the dune sediments control the ⁸⁷Sr values. Based on the mutual variations of ¹⁸O and D, the lakes can be divided into three groups. In Group 1, both ¹⁸O and D values increase from spring to fall. The ¹⁸O and D values in Group 2 decreased from spring to fall. Group 3 are ephemeral lakes that went dry some time during 1992. The data and isotopic modeling show that variations in the ratio of evaporation relative to groundwater inflow, local humidity conditions, and the _a has substantial influence on the isotopic composition. In addition, isotopic behavior in ephemeral lakes can be rather unusual because of the changing activities of water and mineral precipitation and redissolution. The annual and interannual isotopic variability of these lakes which is reflected in the paleonvironmental indicators may be the rule rather than the exception in these types of systems.     

Mercury: Internal structure and thermal evolution

Astronomical observations and cosmochemical calculations suggest that the planet Mercury may be composed of materials which condensed at relatively high temperatures in the primitive solar nebula and may have a basaltic crust similar to parts of the moon. These findings, plus the long standing inference that Mercury is much richer in metallic iron than the other terrestrial planets, provide important constraints which we apply to models of the thermal evolution and density structure of the planet. The thermal history calculations include explicitly the differing thermal properties of iron and silicates and account for core segregation, melting and differentiation of heat sources, and simulated convection during melting. If the U and Th abundances of Mercury are taken from the cosmochemical model of Lewis, then the planet would have fully differentiated a metal core from the silicate mantle for all likely initial temperature distributions and heat transfer properties. Density distributions for the planet are calculated from the mean density and estimates of the present-day temperature. For the fully differentiated model, the moment of inertia C/MR² is 0.325 (J₂=0.302×10^?6). For models with lower heat source abundances, the planet may not yet have differentiated. The density profiles for such models give C/MR²=0.394 (J₂=0.487×10^?6). These results should be useful for preliminary interpretation of the Mariner 10 measurements of Mercury's gravitational field.     

Sounding rocket measurements of the solar soft X-ray irradiance

Measurements of the solar soft X-ray (XUV: 2 nm to 30 nm) irradiance were performed from a sounding rocket payload flown from White Sands Missile Range, New Mexico on 4 October 1993 and again on 3 November 1994. The soft X-ray instrumentation comprised of silicon photodiodes with thin films deposited directly onto their active areas. The deposited material and its thickness in conjunction with the sensitivity of an uncoated diode determine the passband and the sensitivity of these photometric devices. The measurements are interpreted in terms of appropriate SERF 1 (Hinteregger, Fukui, and Gilson, 1981) model solar spectra. It is found that the data are consistent with a solar spectrum that is on average approximately a factor of two times the model solar spectra. It is shown that the measured irradiances are in reasonable agreement with other experiments.     

A complete sample of radio sources in the North Ecliptic Cap, selected at 38 MHz – III. Further imaging observations and the photometric properties of the sample

Further imaging observations of a sample of radio sources in the North Ecliptic Cap are presented and a number of new identifications are made. Using redshifts from spectroscopic data presented in a companion paper by Lacy et al., the photometric properties of the galaxies in the sample are discussed. It is shown that: (1) out to at least z ≈0.6 radio galaxies are good standard candles irrespective of radio luminosity; (2) for 0.6≲ z ≲1 a large fraction of the sample has magnitudes and colours consistent with a non-evolving giant elliptical, and (3) at higher redshifts, where the R -band samples the rest-frame UV flux, most objects have less UV luminosity than expected if they form their stellar populations at a constant rate from a high redshift to z ∼1 in unobscured star-forming regions (assuming an Einstein–de Sitter cosmology). The consequences of these observations are briefly discussed.