Anelastic dynamo models with variable electrical conductivity: an application to gas giants

The observed surface dynamics of Jupiter and Saturn is dominated by a banded system of zonal winds. Their depth remains unclear but they are thought to be confined to the very outer envelopes where hydrogen remains molecular and the electrical conductivity is small. The dynamo maintaining the dipole-dominated magnetic fields of both gas giants likely operates in the deeper interior where hydrogen assumes a metallic state. Here, we present numerical simulations that attempt to model both the zonal winds and the interior dynamo action in an integrated approach. Using the anelastic version of the MHD code MagIC, we explore the effects of density stratification and radial electrical conductivity variation. The electrical conductivity is mostly assumed to remain constant in the thicker inner metallic region and it decays exponentially towards the outer boundary throughout the molecular envelope. Our results show that the combination of stronger density stratification and weaker conducting outer layer is essential for reconciling dipole dominated dynamo action and a fierce equatorial zonal jet. Previous simulations with homogeneous electrical conductivity show that both are merely exclusive, with solutions either having strong zonal winds and multipolar magnetic fields or weak zonal winds and dipole-dominated magnetic fields. All jets tend to be geostrophic and therefore reach right through the convective shell in our simulations. The particular setup explored here allows a strong equatorial jet to remain confined to the weaker conducting outer region where it does not interfere with the deeper seated dynamo action. The flanking mid to high latitude jets, on the other hand, have to remain faint to yield a strongly dipolar magnetic field. The fiercer jets on Jupiter and Saturn only seem compatible with the observed dipolar fields when they remain confined to a weaker conducting outer layer.Comment: 16 pages, 11 figures, 2 tables, submitted to PEP

Stellar populations in superclusters of galaxies

A catalogue of superclusters of galaxies is used to investigate the influence of the supercluster environment on galaxy populations, considering galaxies brighter than M$_r<$-21+5$\log$ h. Empirical spectral synthesis techniques are applied to obtain the stellar population properties of galaxies which belong to superclusters and representative values of stellar population parameters are attributed to each supercluster. We show that richer superclusters present denser environments and older stellar populations. The galaxy populations of superclusters classified as filaments and pancakes are statistically similar, indicating that the morphology of superclusters does not have a significative influence on the stellar populations. Clusters of galaxies within superclusters are also examined in order to evaluate the influence of the supercluster environment on their galaxy properties. Our results suggest that the environment affects galaxy properties but its influence should operate on scales of groups and clusters, more than on the scale of superclusters.Comment: 7 pages, 4 figures; accepted to MNRA

Bridging planets and stars using scaling laws in anelastic spherical shell dynamos

Dynamos operating in the interiors of rapidly rotating planets and low-mass stars might belong to a similar category where rotation plays a vital role. We quantify this similarity using scaling laws. We analyse direct numerical simulations of Boussinesq and anelastic spherical shell dynamos. These dynamos represent simplified models which span from Earth-like planets to rapidly rotating low-mass stars. We find that magnetic field and velocity in these dynamos are related to the available buoyancy power via a simple power law which holds over wide variety of control parameters.Comment: 2 pages; Proceedings of IAUS 302: Magnetic fields throughout stellar evolution (August 2013, Biarritz, France

A calcareous nannofossil and organic geochemical study of marine palaeoenvironmental changes across the Sinemurian/Pliensbachian (early Jurassic, ~191Ma) in Portugal

The Sinemurian/Pliensbachian boundary (~ 191 Ma) is acknowledged as one of the most important steps in the radiation of planktonic organisms, especially primary producers such as dinoflagellates and coccolithophores. To date, there is no detailed study documenting changes in planktonic assemblages related to palaeoceanographic changes across this boundary. The aim of this study is to characterize the palaeoenvironmental changes occurring across the Sinemurian/Pliensbachian boundary at the São Pedro de Moel section (Lusitanian Basin, Portugal) using micropalaeontology and organic geochemistry approaches. Combined calcareous nannofossil assemblage and lipid biomarker data document for a decrease in primary productivity in relation to a major sea-level rise occurring above the boundary. The Lusitanian Basin was particularly restricted during the late Sinemurian with a relatively low sea level, a configuration that led to the recurrent development of black shales. After a sharp sea-level fall, the basin became progressively deeper and more open during the earliest Pliensbachian, subsequently to a major transgression. This sea-level increase seems to have been a global feature and could have been related to the opening of the Hispanic Corridor that connected the Tethys and palaeo-Pacific oceans. The palaeoceanographic and palaeoclimatic changes induced by this opening may have played a role in the diversification of coccolithophores with the first occurrence or colonization of Tethyan waters by placolith-type coccoliths

Retired galaxies: not to be forgotten in the quest of the star formation -- AGN connection

We propose a fresh look at the Main Galaxy Sample of the Sloan Digital Sky Survey by packing the galaxies in stellar mass and redshift bins. We show how important it is to consider the emission-line equivalent widths, in addition to the commonly used emission-line ratios, to properly identify retired galaxies (i.e. galaxies that have stopped forming stars and are ionized by their old stellar populations) and not mistake them for galaxies with low-level nuclear activity. We find that the proportion of star-forming galaxies decreases with decreasing redshift in each mass bin, while that of retired galaxies increases. Galaxies with $M_\star > 10^{11.5} M_\odot$ have formed all their stars at redshift larger than 0.4. The population of AGN hosts is never dominant for galaxy masses larger than $10^{10} M_\odot$. We warn about the effects of stacking galaxy spectra to discuss galaxy properties. We estimate the lifetimes of active galactic nuclei (AGN) relying entirely on demographic arguments --- i.e. without any assumption on the AGN radiative properties. We find upper-limit lifetimes of about 1--5 Gyr for detectable AGN in galaxies with masses between $10^{10}$--$10^{12} M_\odot$. The lifetimes of the AGN-dominated phases are a few $10^8$ yr. Finally, we compare the star-formation histories of star-forming, AGN and retired galaxies as obtained by the spectral synthesis code STARLIGHT. Once the AGN is turned on it inhibits star formation for the next $\sim$ 0.1 Gyr in galaxies with masses around $10^{10} M_\odot$, $\sim$ 1 Gyr in galaxies with masses around $10^{11} M_\odot$.Comment: accepted for MNRAS figure resolution has been degraded with respect to what will be published in MNRA

Multiscale Partition of Unity

We introduce a new Partition of Unity Method for the numerical homogenization of elliptic partial differential equations with arbitrarily rough coefficients. We do not restrict to a particular ansatz space or the existence of a finite element mesh. The method modifies a given partition of unity such that optimal convergence is achieved independent of oscillation or discontinuities of the diffusion coefficient. The modification is based on an orthogonal decomposition of the solution space while preserving the partition of unity property. This precomputation involves the solution of independent problems on local subdomains of selectable size. We deduce quantitative error estimates for the method that account for the chosen amount of localization. Numerical experiments illustrate the high approximation properties even for 'cheap' parameter choices.Comment: Proceedings for Seventh International Workshop on Meshfree Methods for Partial Differential Equations, 18 pages, 3 figure

Clone size distributions in networks of genetic similarity

We build networks of genetic similarity in which the nodes are organisms sampled from biological populations. The procedure is illustrated by constructing networks from genetic data of a marine clonal plant. An important feature in the networks is the presence of clone subgraphs, i.e. sets of organisms with identical genotype forming clones. As a first step to understand the dynamics that has shaped these networks, we point up a relationship between a particular degree distribution and the clone size distribution in the populations. We construct a dynamical model for the population dynamics, focussing on the dynamics of the clones, and solve it for the required distributions. Scale free and exponentially decaying forms are obtained depending on parameter values, the first type being obtained when clonal growth is the dominant process. Average distributions are dominated by the power law behavior presented by the fastest replicating populations.Comment: 17 pages, 4 figures. One figure improved and other minor changes. To appear in Physica