DIRK Schemes with High Weak Stage Order

Runge-Kutta time-stepping methods in general suffer from order reduction: the observed order of convergence may be less than the formal order when applied to certain stiff problems. Order reduction can be avoided by using methods with high stage order. However, diagonally-implicit Runge-Kutta (DIRK) schemes are limited to low stage order. In this paper we explore a weak stage order criterion, which for initial boundary value problems also serves to avoid order reduction, and which is compatible with a DIRK structure. We provide specific DIRK schemes of weak stage order up to 3, and demonstrate their performance in various examples.Comment: 10 pages, 5 figure

Disruption of a Proto-Planetary Disk by the Black Hole at the Milky Way Centre

Recently, an ionized cloud of gas was discovered plunging toward the supermassive black hole, SgrA*, at the centre of the Milky Way. The cloud is being tidally disrupted along its path to closest approach at ~3100 Schwarzschild radii from the black hole. Here, we show that the observed properties of this cloud of gas can naturally be produced by a proto-planetary disk surrounding a low-mass star, which was scattered from the observed ring of young stars orbiting SgrA*. As the young star approaches the black hole, its disk experiences both photo-evaporation and tidal disruption, producing a cloud. Our model implies that planets form in the Galactic centre, and that tidal debris from proto-planetary disks can flag low mass stars which are otherwise too faint to be detected.Comment: Accepted to Nature Communications; new Figure 4b provides predicted Br-gamma emission as a function of tim

Seroprevalence of HTLV-1 and HTLV-2 amongst mothers and children in Malawi within the context of a systematic review and meta-analysis of HTLV seroprevalence in Africa

OBJECTIVES: Human T‐lymphotropic virus (HTLV)‐1 causes T‐cell leukaemia and myelopathy. Together with HTLV‐2, it is endemic in some African nations. Seroprevalence data from Malawi are scarce, with no reports on associated disease incidence. HTLV seroprevalence and type were tested in 418 healthy mothers from Malawi. In addition, we tested the sera of 534 children to investigate mother‐to‐child transmission. To provide context, we conducted a systematic review and meta‐analysis of HTLV seroprevalence in African women and children. METHODS: Stored samples from a previous childhood cancer and BBV study were analysed. ELISA was used for HTLV screening followed by immunoblot for confirmation and typing. Standard methods were used for the systematic review. RESULTS: HTLV seroprevalence was 2.6% (11/418) in mothers and 2.2% (12/534) in children. Three mothers carried HTLV‐1 alone, seven had HTLV‐2 and one was dually infected. Three children carried HTLV‐1 alone, seven had HTLV‐2 and two were dually infected. Only two corresponding mothers of the 12 HTLV‐positive children were HTLV positive. The systematic review included 66 studies of women and 13 of children conducted in 25 African countries. Seroprevalence of HTLV‐1 varied from 0 to 17% and of HTLV‐2 from 0 to 4%. CONCLUSIONS: In contrast to findings from other studies in Africa, the seroprevalence of HTLV‐2 was higher than that of HTLV‐1 in Malawi and one of the highest for the African region. The lack of mother–child concordance suggests alternative sources of infection among children. Our data and analyses contribute to HTLV prevalence mapping in Africa

The Mass Function of Newly Formed Stars (Review)

The topic of the stellar "original mass function" has a nearly 50 year history,dating to the publication in 1955 of Salpeter's seminal paper. In this review I discuss the many more recent results that have emerged on the initial mass function (IMF), as it is now called, from studies over the last decade of resolved populations in star forming regions and young open clusters.Comment: 9 pages, 1 figure; to appear in "The Dense Instellar Medium in Galaxies -- 4'th Cologne-Bonn-Zermatt-Symposium" editted by S. Pfalzner, C. Kramer, C. Straubmeier and A. Heithausen, Springer-Verlag (2004

Initial fixation placement in face images is driven by top-down guidance

The eyes are often inspected first and for longer period during face exploration. To examine whether this saliency of the eye region at the early stage of face inspection is attributed to its local structure properties or to the knowledge of its essence in facial communication, in this study we investigated the pattern of eye movements produced by rhesus monkeys (Macaca mulatta) as they free viewed images of monkey faces. Eye positions were recorded accurately using implanted eye coils, while images of original faces, faces with scrambled eyes, and scrambled faces except for the eyes were presented on a computer screen. The eye region in the scrambled faces attracted the same proportion of viewing time and fixations as it did in the original faces, even the scrambled eyes attracted substantial proportion of viewing time and fixations. Furthermore, the monkeys often made the first saccade towards to the location of the eyes regardless of image content. Our results suggest that the initial fixation placement in faces is driven predominantly by ‘top-down’ or internal factors, such as the prior knowledge of the location of “eyes” within the context of a face

Degeneracy: a link between evolvability, robustness and complexity in biological systems

A full accounting of biological robustness remains elusive; both in terms of the mechanisms by which robustness is achieved and the forces that have caused robustness to grow over evolutionary time. Although its importance to topics such as ecosystem services and resilience is well recognized, the broader relationship between robustness and evolution is only starting to be fully appreciated. A renewed interest in this relationship has been prompted by evidence that mutational robustness can play a positive role in the discovery of adaptive innovations (evolvability) and evidence of an intimate relationship between robustness and complexity in biology. This paper offers a new perspective on the mechanics of evolution and the origins of complexity, robustness, and evolvability. Here we explore the hypothesis that degeneracy, a partial overlap in the functioning of multi-functional components, plays a central role in the evolution and robustness of complex forms. In support of this hypothesis, we present evidence that degeneracy is a fundamental source of robustness, it is intimately tied to multi-scaled complexity, and it establishes conditions that are necessary for system evolvability

Use of the bootstrap in analysing cost data from cluster randomised trials: some simulation results

BACKGROUND: This work has investigated under what conditions confidence intervals around the differences in mean costs from a cluster RCT are suitable for estimation using a commonly used cluster-adjusted bootstrap in preference to methods that utilise the Huber-White robust estimator of variance. The bootstrap's main advantage is in dealing with skewed data, which often characterise patient costs. However, it is insufficiently well recognised that one method of adjusting the bootstrap to deal with clustered data is only valid in large samples. In particular, the requirement that the number of clusters randomised should be large would not be satisfied in many cluster RCTs performed to date. METHODS: The performances of confidence intervals for simple differences in mean costs utilising a robust (cluster-adjusted) standard error and from two cluster-adjusted bootstrap procedures were compared in terms of confidence interval coverage in a large number of simulations. Parameters varied included the intracluster correlation coefficient, the sample size and the distributions used to generate the data. RESULTS: The bootstrap's advantage in dealing with skewed data was found to be outweighed by its poor confidence interval coverage when the number of clusters was at the level frequently found in cluster RCTs in practice. Simulations showed that confidence intervals based on robust methods of standard error estimation achieved coverage rates between 93.5% and 94.8% for a 95% nominal level whereas those for the bootstrap ranged between 86.4% and 93.8%. CONCLUSION: In general, 24 clusters per treatment arm is probably the minimum number for which one would even begin to consider the bootstrap in preference to traditional robust methods, for the parameter combinations investigated here. At least this number of clusters and extremely skewed data would be necessary for the bootstrap to be considered in favour of the robust method. There is a need for further investigation of more complex bootstrap procedures if economic data from cluster RCTs are to be analysed appropriately

Fast fluorescence microscopy for imaging the dynamics of embryonic development

Live imaging has gained a pivotal role in developmental biology since it increasingly allows real-time observation of cell behavior in intact organisms. Microscopes that can capture the dynamics of ever-faster biological events, fluorescent markers optimal for in vivo imaging, and, finally, adapted reconstruction and analysis programs to complete data flow all contribute to this success. Focusing on temporal resolution, we discuss how fast imaging can be achieved with minimal prejudice to spatial resolution, photon count, or to reliably and automatically analyze images. In particular, we show how integrated approaches to imaging that combine bright fluorescent probes, fast microscopes, and custom post-processing techniques can address the kinetics of biological systems at multiple scales. Finally, we discuss remaining challenges and opportunities for further advances in this field

Circadian Rhythms in Visual Responsiveness in the Behaviorally Arrhythmic Drosophila Clock Mutant ClkJrk

An organism's biological day is characterized by a pattern of anticipatory physiological and behavioral changes that are governed by circadian clocks to align with the 24-h cycling environment. Here, we used flash electroretinograms (ERGs) and steady-state visually evoked potentials (SSVEPs) to examine how visual responsiveness in wild-type Drosophila melanogaster and the circadian clock mutant ClkJrk varies over circadian time. We show that the ERG parameters of wild-type flies vary over the circadian day, with a higher luminance response during the subjective night. The SSVEP response that assesses contrast sensitivity also showed a time-of-day dependence, including 2 prominent peaks within a 24-h period and a maximal response at the end of the subjective day, indicating a tradeoff between luminance and contrast sensitivity. Moreover, the behaviorally arrhythmic ClkJrk mutants maintained a circadian profile in both luminance and contrast sensitivity, but unlike the wild-types, which show bimodal profiles in their visual response, ClkJrk flies show a weakening of the bimodal character, with visual responsiveness tending to peak once a day. We conclude that the ClkJrk mutation mainly affects 1 of 2 functionally coupled oscillators and that the visual system is partially separated from the locomotor circadian circuits that drive bouts of morning and evening activity. As light exposure is a major mechanism for entrainment, our work suggests that a detailed temporal analysis of electrophysiological responses is warranted to better identify the time window at which circadian rhythms are most receptive to light-induced phase shifting