SpECTRE: A Task-based Discontinuous Galerkin Code for Relativistic Astrophysics

We introduce a new relativistic astrophysics code, SpECTRE, that combines a discontinuous Galerkin method with a task-based parallelism model. SpECTRE's goal is to achieve more accurate solutions for challenging relativistic astrophysics problems such as core-collapse supernovae and binary neutron star mergers. The robustness of the discontinuous Galerkin method allows for the use of high-resolution shock capturing methods in regions where (relativistic) shocks are found, while exploiting high-order accuracy in smooth regions. A task-based parallelism model allows efficient use of the largest supercomputers for problems with a heterogeneous workload over disparate spatial and temporal scales. We argue that the locality and algorithmic structure of discontinuous Galerkin methods will exhibit good scalability within a task-based parallelism framework. We demonstrate the code on a wide variety of challenging benchmark problems in (non)-relativistic (magneto)-hydrodynamics. We demonstrate the code's scalability including its strong scaling on the NCSA Blue Waters supercomputer up to the machine's full capacity of 22,380 nodes using 671,400 threads.Comment: 41 pages, 13 figures, and 7 tables. Ancillary data contains simulation input file

Computational Relativistic Astrophysics With Adaptive Mesh Refinement: Testbeds

We have carried out numerical simulations of strongly gravitating systems based on the Einstein equations coupled to the relativistic hydrodynamic equations using adaptive mesh refinement (AMR) techniques. We show AMR simulations of NS binary inspiral and coalescence carried out on a workstation having an accuracy equivalent to that of a $1025^3$ regular unigrid simulation, which is, to the best of our knowledge, larger than all previous simulations of similar NS systems on supercomputers. We believe the capability opens new possibilities in general relativistic simulations.Comment: 7 pages, 16 figure

Pure O-sequences and matroid h-vectors

We study Stanley's long-standing conjecture that the h-vectors of matroid simplicial complexes are pure O-sequences. Our method consists of a new and more abstract approach, which shifts the focus from working on constructing suitable artinian level monomial ideals, as often done in the past, to the study of properties of pure O-sequences. We propose a conjecture on pure O-sequences and settle it in small socle degrees. This allows us to prove Stanley's conjecture for all matroids of rank 3. At the end of the paper, using our method, we discuss a first possible approach to Stanley's conjecture in full generality. Our technical work on pure O-sequences also uses very recent results of the third author and collaborators.Comment: Contains several changes/updates with respect to the previous version. In particular, a discussion of a possible approach to the general case is included at the end. 13 pages. To appear in the Annals of Combinatoric

Edge-Based Compartmental Modeling for Infectious Disease Spread Part III: Disease and Population Structure

We consider the edge-based compartmental models for infectious disease spread introduced in Part I. These models allow us to consider standard SIR diseases spreading in random populations. In this paper we show how to handle deviations of the disease or population from the simplistic assumptions of Part I. We allow the population to have structure due to effects such as demographic detail or multiple types of risk behavior the disease to have more complicated natural history. We introduce these modifications in the static network context, though it is straightforward to incorporate them into dynamic networks. We also consider serosorting, which requires using the dynamic network models. The basic methods we use to derive these generalizations are widely applicable, and so it is straightforward to introduce many other generalizations not considered here

Vasoactive intestinal polypeptide mediates circadian rhythms in Mammalian olfactory bulb and olfaction

Accumulating evidence suggests that the olfactory bulbs (OBs) function as an independent circadian system regulating daily rhythms in olfactory performance. However, the cells and signals in the olfactory system that generate and coordinate these circadian rhythms are unknown. Using real-time imaging of gene expression, we found that the isolated olfactory epithelium and OB, but not the piriform cortex, express similar, sustained circadian rhythms in PERIOD2 (PER2). In vivo, PER2 expression in the OB of mice is circadian, approximately doubling with a peak around subjective dusk. Furthermore, mice exhibit circadian rhythms in odor detection performance with a peak at approximately subjective dusk. We also found that circadian rhythms in gene expression and odor detection performance require vasoactive intestinal polypeptide (VIP) or its receptor VPAC2R. VIP is expressed, in a circadian manner, in interneurons in the external plexiform and periglomerular layers, whereas VPAC2R is expressed in mitral and external tufted cells in the OB. Together, these results indicate that VIP signaling modulates the output from the OB to maintain circadian rhythms in the mammalian olfactory system.Fil: Kang Miller, Jae Eun. Washington University in St. Louis; Estados UnidosFil: Granados Fuentes, Daniel. Washington University in St. Louis; Estados UnidosFil: Wang, Thomas. Washington University in St. Louis; Estados UnidosFil: Marpegan, Luciano. Washington University in St. Louis; Estados Unidos. Universidad Nacional de Quilmes. Departamento de Ciencia y Tecnología; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas; ArgentinaFil: Holy, Timothy E.. Washington University in St. Louis; Estados UnidosFil: Herzog, Erik D.. Washington University in St. Louis; Estados Unido

Caenorhabditis elegans ETR-1/CELF has broad effects on the muscle cell transcriptome, including genes that regulate translation and neuroblast migration

Migration of neuroblasts and neurons from their birthplace is central to the formation of neural circuits and networks. ETR-1 is the Caenorhabditis elegans homolog of the CELF1 (CUGBP, ELAV-like family 1) RNA-processing factor involved in neuromuscular disorders. etr-1 regulates body wall muscle differentiation. Our previous work showed that etr-1 in muscle has a non-autonomous role in neuronal migration, suggesting that ETR-1 is involved in the production of a signal emanating from body wall muscle that controls neuroblast migration and that interacts with Wnt signaling. etr-1 is extensively alternatively-spliced, and we identified the viable etr-1(lq61) mutant, caused by a stop codon in alternatively-spliced exon 8 and only affecting etr-1 isoforms containing exon 8. We took advantage of viable etr-1(lq61) to identify potential RNA targets of ETR-1 in body wall muscle using a combination of fluorescence activated cell sorting (FACS) of body wall muscles from wild-type and etr-1(lq61) and subsequent RNA-seq. This analysis revealed genes whose splicing and transcript levels were controlled by ETR-1 exon 8 isoforms, and represented a broad spectrum of genes involved in muscle differentiation, myofilament lattice structure, and physiology. Genes with transcripts underrepresented in etr-1(lq61) included those involved in ribosome function and translation, similar to potential CELF1 targets identified in chick cardiomyocytes. This suggests that at least some targets of ETR-1 might be conserved in vertebrates, and that ETR-1 might generally stimulate translation in muscles. As proof-of-principle, a functional analysis of a subset of ETR-1 targets revealed genes involved in AQR and PQR neuronal migration. One such gene, lev-11/tropomyosin, requires ETR-1 for alternative splicing, and another, unc-52/perlecan, requires ETR-1 for the production of long isoforms containing 3′ exons. In sum, these studies identified gene targets of ETR-1/CELF1 in muscles, which included genes involved in muscle development and physiology, and genes with novel roles in neuronal migration

Radio Sources from a 31 GHz Sky Survey with the Sunyaev-Zel'dovich Array

We present the first sample of 31-GHz selected sources to flux levels of 1 mJy. From late 2005 to mid 2007, the Sunyaev-Zel'dovich Array (SZA) observed 7.7 square degrees of the sky at 31 GHz to a median rms of 0.18 mJy/beam. We identify 209 sources at greater than 5 sigma significance in the 31 GHz maps, ranging in flux from 0.7 mJy to ~200 mJy. Archival NVSS data at 1.4 GHz and observations at 5 GHz with the Very Large Array are used to characterize the sources. We determine the maximum-likelihood integrated source count to be N(>S) = (27.2 +- 2.5) deg^-2 x (S_mJy)^(-1.18 +- 0.12) over the flux range 0.7 - 15 mJy. This result is significantly higher than predictions based on 1.4-GHz selected samples, a discrepancy which can be explained by a small shift in the spectral index distribution for faint 1.4-GHz sources. From comparison with previous measurements of sources within the central arcminute of massive clusters, we derive an overdensity of 6.8 +- 4.4, relative to field sources.Comment: 13 pages, 5 figure

Non-Detection of Gravitationally Redshifted Absorption Lines in the X-ray Burst Spectra of GS 1826-24

During a 200 ks observation with the XMM-Newton Reflection Grating Spectrometer, we detected 16 type-I X-ray bursts from GS 1826-24. We combined the burst spectra in an attempt to measure the gravitational redshifts from the surface of the neutron star. We divided the composite GS 1826-24 burst spectrum into three groups based on the blackbody temperature during the bursts. The spectra do not show any obvious discrete absorption lines. We compare our observations with those of EXO 0748-676.Comment: 4 pages, 4 figures; accepted for publication in ApJ

Impaired β-glucocerebrosidase activity and processing in frontotemporal dementia due to progranulin mutations.

Loss-of-function mutations in progranulin (GRN) are a major autosomal dominant cause of frontotemporal dementia. Most pathogenic GRN mutations result in progranulin haploinsufficiency, which is thought to cause frontotemporal dementia in GRN mutation carriers. Progranulin haploinsufficiency may drive frontotemporal dementia pathogenesis by disrupting lysosomal function, as patients with GRN mutations on both alleles develop the lysosomal storage disorder neuronal ceroid lipofuscinosis, and frontotemporal dementia patients with GRN mutations (FTD-GRN) also accumulate lipofuscin. The specific lysosomal deficits caused by progranulin insufficiency remain unclear, but emerging data indicate that progranulin insufficiency may impair lysosomal sphingolipid-metabolizing enzymes. We investigated the effects of progranulin insufficiency on sphingolipid-metabolizing enzymes in the inferior frontal gyrus of FTD-GRN patients using fluorogenic activity assays, biochemical profiling of enzyme levels and posttranslational modifications, and quantitative neuropathology. Of the enzymes studied, only β-glucocerebrosidase exhibited impairment in FTD-GRN patients. Brains from FTD-GRN patients had lower activity than controls, which was associated with lower levels of mature β-glucocerebrosidase protein and accumulation of insoluble, incompletely glycosylated β-glucocerebrosidase. Immunostaining revealed loss of neuronal β-glucocerebrosidase in FTD-GRN patients. To investigate the effects of progranulin insufficiency on β-glucocerebrosidase outside of the context of neurodegeneration, we investigated β-glucocerebrosidase activity in progranulin-insufficient mice. Brains from Grn-/- mice had lower β-glucocerebrosidase activity than wild-type littermates, which was corrected by AAV-progranulin gene therapy. These data show that progranulin insufficiency impairs β-glucocerebrosidase activity in the brain. This effect is strongest in neurons and may be caused by impaired β-glucocerebrosidase processing

Effects of Heterogeneous and Clustered Contact Patterns on Infectious Disease Dynamics

The spread of infectious diseases fundamentally depends on the pattern of contacts between individuals. Although studies of contact networks have shown that heterogeneity in the number of contacts and the duration of contacts can have far-reaching epidemiological consequences, models often assume that contacts are chosen at random and thereby ignore the sociological, temporal and/or spatial clustering of contacts. Here we investigate the simultaneous effects of heterogeneous and clustered contact patterns on epidemic dynamics. To model population structure, we generalize the configuration model which has a tunable degree distribution (number of contacts per node) and level of clustering (number of three cliques). To model epidemic dynamics for this class of random graph, we derive a tractable, low-dimensional system of ordinary differential equations that accounts for the effects of network structure on the course of the epidemic. We find that the interaction between clustering and the degree distribution is complex. Clustering always slows an epidemic, but simultaneously increasing clustering and the variance of the degree distribution can increase final epidemic size. We also show that bond percolation-based approximations can be highly biased if one incorrectly assumes that infectious periods are homogeneous, and the magnitude of this bias increases with the amount of clustering in the network. We apply this approach to model the high clustering of contacts within households, using contact parameters estimated from survey data of social interactions, and we identify conditions under which network models that do not account for household structure will be biased