Local starburst galaxies and their descendants

Despite strong interest in the starburst (hereafter SB) phenomenon, the concept remains ill-defined. We use a strict definition of SB to examine the statistical properties of local SB and post-starburst (hereafter PB) galaxies. We also seek relationships to active galaxies. Potential SB galaxies are selected from the SDSS DR7 and their stellar content is analysed. We apply an age dependent dust attenuation correction and derive star formation rates (SFR), ages and masses of the young and old populations. The photometric masses nicely agree with dynamical masses derived from the H-alpha emission line width. To select SB galaxies, we use the birthrate parameter b=SFR/, requiring b>=3. The PB sample is selected from the citerion EW(Hdelta_abs)>=6 A. Only 1% of star-forming galaxies are found to be SB galaxies. They contribute 3-6% to the stellar production and are therefore unimportant for the local star formation activity. The median SB age is 70 Myr, roughly independent of mass. The b-parameter strongly depends on burst age. Values close to b=60 are found at ages ~10 Myr, while almost no SBs are found at ages >1 Gyr. The median baryonic burst mass fraction of sub-L* galaxies is 5%, decreasing slowly with mass. The median mass fraction of the recent burst in the PB sample is 5-10%. The age-mass distribution of the progenitors of the PBs is bimodal with a break at log(M)~10.6 above which the ages are doubled. The SB and PB luminosity functions (hereafter LFs) follow each other closely until M_r~-21, when AGNs begin to dominate. The PB LF continues to follow the AGN LF while SB loose significance. This suggests that the number of luminous SBs is underestimated by about one dex at high luminosities, due to large amounts of dust and/or AGN blending. It also indicates that the SB phase preceded the AGN phase. We also discuss the conditions for global gas outflow caused by stellar feedback.Comment: Accepted for publication in Astronomy and Astrophysics. This is an extended, substantially revised and corrected version with partly modified conclusion

Local Finite Element Approximation of Sobolev Differential Forms

We address fundamental aspects in the approximation theory of vector-valued finite element methods, using finite element exterior calculus as a unifying framework. We generalize the Cl\'ement interpolant and the Scott-Zhang interpolant to finite element differential forms, and we derive a broken Bramble-Hilbert Lemma. Our interpolants require only minimal smoothness assumptions and respect partial boundary conditions. This permits us to state local error estimates in terms of the mesh size. Our theoretical results apply to curl-conforming and divergence-conforming finite element methods over simplicial triangulations.Comment: 22 pages. Comments welcom

3D mesh processing using GAMer 2 to enable reaction-diffusion simulations in realistic cellular geometries

Recent advances in electron microscopy have enabled the imaging of single cells in 3D at nanometer length scale resolutions. An uncharted frontier for in silico biology is the ability to simulate cellular processes using these observed geometries. Enabling such simulations requires watertight meshing of electron micrograph images into 3D volume meshes, which can then form the basis of computer simulations of such processes using numerical techniques such as the Finite Element Method. In this paper, we describe the use of our recently rewritten mesh processing software, GAMer 2, to bridge the gap between poorly conditioned meshes generated from segmented micrographs and boundary marked tetrahedral meshes which are compatible with simulation. We demonstrate the application of a workflow using GAMer 2 to a series of electron micrographs of neuronal dendrite morphology explored at three different length scales and show that the resulting meshes are suitable for finite element simulations. This work is an important step towards making physical simulations of biological processes in realistic geometries routine. Innovations in algorithms to reconstruct and simulate cellular length scale phenomena based on emerging structural data will enable realistic physical models and advance discovery at the interface of geometry and cellular processes. We posit that a new frontier at the intersection of computational technologies and single cell biology is now open.Comment: 39 pages, 14 figures. High resolution figures and supplemental movies available upon reques

Semilinear mixed problems on Hilbert complexes and their numerical approximation

Arnold, Falk, and Winther recently showed [Bull. Amer. Math. Soc. 47 (2010), 281-354] that linear, mixed variational problems, and their numerical approximation by mixed finite element methods, can be studied using the powerful, abstract language of Hilbert complexes. In another recent article [arXiv:1005.4455], we extended the Arnold-Falk-Winther framework by analyzing variational crimes (a la Strang) on Hilbert complexes. In particular, this gave a treatment of finite element exterior calculus on manifolds, generalizing techniques from surface finite element methods and recovering earlier a priori estimates for the Laplace-Beltrami operator on 2- and 3-surfaces, due to Dziuk [Lecture Notes in Math., vol. 1357 (1988), 142-155] and later Demlow [SIAM J. Numer. Anal., 47 (2009), 805-827], as special cases. In the present article, we extend the Hilbert complex framework in a second distinct direction: to the study of semilinear mixed problems. We do this, first, by introducing an operator-theoretic reformulation of the linear mixed problem, so that the semilinear problem can be expressed as an abstract Hammerstein equation. This allows us to obtain, for semilinear problems, a priori solution estimates and error estimates that reduce to the Arnold-Falk-Winther results in the linear case. We also consider the impact of variational crimes, extending the results of our previous article to these semilinear problems. As an immediate application, this new framework allows for mixed finite element methods to be applied to semilinear problems on surfaces.Comment: 22 pages; v2: major revision, particularly sharpening of error estimates in Section

Emotional and Adrenocortical Responses of Infants to the Strange Situation: The Differential Function of Emotional Expression

The aim of the study was to investigate biobehavioural organisation in infants with different qualities of attachment. Quality of attachment (security and disorganisation), emotional expression, and adrenocortical stress reactivity were investigated in a sample of 106 infants observed during Ainsworth’s Strange Situation at the age of 12 months. In addition, behavioural inhibition was assessed from maternal reports. As expected, securely attached infants did not show an adrenocortical response. Regarding the traditionally defined insecurely attached groups, adrenocortical activation during the strange situation was found for the ambivalent group, but not for the avoidant one. Previous ndings of increased adrenocortical activity in disorganised infants could not be replicated. In line with previous ndings, adrenocortical activation was most prominent in insecure infants with high behavioural inhibition indicating the function of a secure attachment relationship as a social buffer against less adaptive temperamental dispositions. Additional analyses indicated that adrenocortical reactivity and behavioural distress were not based on common activation processes. Biobehavioural associations within the different attachment groups suggest that biobehavioural processes in securely attached infants may be different from those in insecurely attached and disorganised groups. Whereas a coping model may be applied to describe the biobehavioural organisation of secure infants, an arousal model explanation may be more appropriate for the other groups

Modeling effects of L-type ca(2+) current and na(+)-ca(2+) exchanger on ca(2+) trigger flux in rabbit myocytes with realistic T-tubule geometries.

The transverse tubular system of rabbit ventricular myocytes consists of cell membrane invaginations (t-tubules) that are essential for efficient cardiac excitation-contraction coupling. In this study, we investigate how t-tubule micro-anatomy, L-type Ca(2+) channel (LCC) clustering, and allosteric activation of Na(+)/Ca(2+) exchanger by L-type Ca(2+) current affects intracellular Ca(2+) dynamics. Our model includes a realistic 3D geometry of a single t-tubule and its surrounding half-sarcomeres for rabbit ventricular myocytes. The effects of spatially distributed membrane ion-transporters (LCC, Na(+)/Ca(2+) exchanger, sarcolemmal Ca(2+) pump, and sarcolemmal Ca(2+) leak), and stationary and mobile Ca(2+) buffers (troponin C, ATP, calmodulin, and Fluo-3) are also considered. We used a coupled reaction-diffusion system to describe the spatio-temporal concentration profiles of free and buffered intracellular Ca(2+). We obtained parameters from voltage-clamp protocols of L-type Ca(2+) current and line-scan recordings of Ca(2+) concentration profiles in rabbit cells, in which the sarcoplasmic reticulum is disabled. Our model results agree with experimental measurements of global Ca(2+) transient in myocytes loaded with 50 μM Fluo-3. We found that local Ca(2+) concentrations within the cytosol and sub-sarcolemma, as well as the local trigger fluxes of Ca(2+) crossing the cell membrane, are sensitive to details of t-tubule micro-structure and membrane Ca(2+) flux distribution. The model additionally predicts that local Ca(2+) trigger fluxes are at least threefold to eightfold higher than the whole-cell Ca(2+) trigger flux. We found also that the activation of allosteric Ca(2+)-binding sites on the Na(+)/Ca(2+) exchanger could provide a mechanism for regulating global and local Ca(2+) trigger fluxes in vivo. Our studies indicate that improved structural and functional models could improve our understanding of the contributions of L-type and Na(+)/Ca(2+) exchanger fluxes to intracellular Ca(2+) dynamics

Resistant starch and protein intake enhances fat oxidation and feelings of fullness in lean and overweight/obese women

BACKGROUND: Diets high in either resistant starch or protein have been shown to aid in weight management. We examined the effects of meals high in non-resistant or resistant starch with and without elevated protein intake on substrate utilization, energy expenditure, and satiety in lean and overweight/obese women. METHODS: Women of varying levels of adiposity consumed one of four pancake test meals in a single-blind, randomized crossover design: 1) waxy maize (control) starch (WMS); 2) waxy maize starch and whey protein (WMS+WP); 3) resistant starch (RS); or 4) RS and whey protein (RS+WP). RESULTS: Total post-prandial energy expenditure did not differ following any of the four test meals (WMS = 197.9 ± 8.9; WMS+WP = 188 ± 8.1; RS = 191.9 ± 8.9; RS+WP = 195.8 ± 8.7, kcals/180 min), although the combination of RS+WP, but not either intervention alone, significantly increased (P <0.01) fat oxidation (WMS = 89.5 ± 5.4; WMS+WP = 84.5 ± 7.2; RS = 97.4 ± 5.4; RS+WP = 107.8 ± 5.4, kcals/180 min). Measures of fullness increased (125 % vs. 45 %) and hunger decreased (55 % vs. 16 %) following WP supplemented versus non-whey conditions (WMS+WP, RS+WP vs. WMS, RS), whereas circulating hunger and satiety factors were not different among any of the test meals. However, peptide YY (PYY) was significantly elevated at 180 min following RS+WP meal. CONCLUSIONS: The combined consumption of dietary resistant starch and protein increases fat oxidation, PYY, and enhances feelings of satiety and fullness to levels that may be clinically relevant if maintained under chronic conditions. This trial was registered at clinicaltrials.gov as NCT02418429

An Open Source Mesh Generation Platform for Biophysical Modeling Using Realistic Cellular Geometries

Advances in imaging methods such as electron microscopy, tomography and other modalities are enabling high-resolution reconstructions of cellular and organelle geometries. Such advances pave the way for using these geometries for biophysical and mathematical modeling once these data can be represented as a geometric mesh, which, when carefully conditioned, enables the discretization and solution of partial differential equations. In this study, we outline the steps for a na\"ive user to approach GAMer 2, a mesh generation code written in C++ designed to convert structural datasets to realistic geometric meshes, while preserving the underlying shapes. We present two example cases, 1) mesh generation at the subcellular scale as informed by electron tomography, and 2) meshing a protein with structure from x-ray crystallography. We further demonstrate that the meshes generated by GAMer are suitable for use with numerical methods. Together, this collection of libraries and tools simplifies the process of constructing realistic geometric meshes from structural biology data.Comment: 6 pages and 4 figures. Supplemental Movie available upon reques

Barrett-Crane spin foam model from generalized BF-type action for gravity

We study a generalized action for gravity as a constrained BF theory, and its relationship with the Plebanski action. We analyse the discretization of the constraints and the spin foam quantization of the theory, showing that it leads naturally to the Barrett-Crane spin foam model for quantum gravity. Our analysis holds true in both the Euclidean and Lorentzian formulation.Comment: 15 pages, revtex; a sign corrected (area spectrum); some of these results were presented in a preliminary form in gr-qc/0103081; v2: improved presentation of the results, some changes in the text; to appear in Phys. Rev.