Too hot to handle? Analytic solutions for massive neutrino or warm dark matter cosmologies

We obtain novel closed-form solutions to the Friedmann equation for cosmological models containing a component whose equation of state is that of radiation (w =1/3) at early times and that of cold pressureless matter (w=0) at late times. The equation of state smoothly transitions from the early- to late-time behaviour and exactly describes the evolution of a species with a Dirac delta function distribution inmomentum magnitudes | p0| (i.e. all particles have the same | p0|). Such a component, here termed 'hot matter', is an approximate model for both neutrinos and warm dark matter. We consider it alone and in combination with cold matter and with radiation, also obtaining closed-form solutions for the growth of superhorizon perturbations in each case. The idealized model recovers t(a) to better than 1.5 per cent accuracy for all a relative to a Fermi-Dirac distribution (as describes neutrinos). We conclude by adding the second moment of the distribution to our exact solution and then generalizing to include all moments of an arbitrary momentum distribution in a closed-form solution

Automatic Kalman-filter-based wavelet shrinkage denoising of 1D stellar spectra

We propose a non-parametric method to denoise 1D stellar spectra based on wavelet shrinkage followed by adaptive Kalman thresholding. Wavelet shrinkage denoising involves applying the discrete wavelet transform (DWT) to the input signal, 'shrinking' certain frequency components in the transform domain, and then applying inverse DWT to the reduced components. The performance of this procedure is influenced by the choice of base wavelet, the number of decomposition levels, and the thresholding function. Typically, these parameters are chosen by 'trial and error', which can be strongly dependent on the properties of the data being denoised. We here introduce an adaptive Kalman-filter-based thresholding method that eliminates the need for choosing the number of decomposition levels. We use the 'Haar' wavelet basis, which we found to provide excellent filtering for 1D stellar spectra, at a low computational cost. We introduce various levels of Poisson noise into synthetic PHOENIX spectra, and test the performance of several common denoising methods against our own. It proves superior in terms of noise suppression and peak shape preservation. We expect it may also be of use in automatically and accurately filtering low signal-to-noise galaxy and quasar spectra obtained from surveys such as SDSS, Gaia, LSST, PESSTO, VANDELS, LEGA-C, and DESI

Quasi-Static Fractures in Disordered Media and Iterated Conformal Maps

We study the geometrical characteristic of quasi-static fractures in disordered media, using iterated conformal maps to determine the evolution of the fracture pattern. This method allows an efficient and accurate solution of the Lam\'e equations without resorting to lattice models. Typical fracture patterns exhibit increased ramification due to the increase of the stress at the tips. We find the roughness exponent of the experimentally relevant backbone of the fracture pattern; it crosses over from about 0.5 for small scales to about 0.75 for large scales, in excellent agreement with experiments. We propose that this cross-over reflects the increased ramification of the fracture pattern.Comment: submitted to Physical Review Letter

Scalar and vector Slepian functions, spherical signal estimation and spectral analysis

It is a well-known fact that mathematical functions that are timelimited (or spacelimited) cannot be simultaneously bandlimited (in frequency). Yet the finite precision of measurement and computation unavoidably bandlimits our observation and modeling scientific data, and we often only have access to, or are only interested in, a study area that is temporally or spatially bounded. In the geosciences we may be interested in spectrally modeling a time series defined only on a certain interval, or we may want to characterize a specific geographical area observed using an effectively bandlimited measurement device. It is clear that analyzing and representing scientific data of this kind will be facilitated if a basis of functions can be found that are "spatiospectrally" concentrated, i.e. "localized" in both domains at the same time. Here, we give a theoretical overview of one particular approach to this "concentration" problem, as originally proposed for time series by Slepian and coworkers, in the 1960s. We show how this framework leads to practical algorithms and statistically performant methods for the analysis of signals and their power spectra in one and two dimensions, and, particularly for applications in the geosciences, for scalar and vectorial signals defined on the surface of a unit sphere.Comment: Submitted to the 2nd Edition of the Handbook of Geomathematics, edited by Willi Freeden, Zuhair M. Nashed and Thomas Sonar, and to be published by Springer Verlag. This is a slightly modified but expanded version of the paper arxiv:0909.5368 that appeared in the 1st Edition of the Handbook, when it was called: Slepian functions and their use in signal estimation and spectral analysi

Integrated information increases with fitness in the evolution of animats

One of the hallmarks of biological organisms is their ability to integrate disparate information sources to optimize their behavior in complex environments. How this capability can be quantified and related to the functional complexity of an organism remains a challenging problem, in particular since organismal functional complexity is not well-defined. We present here several candidate measures that quantify information and integration, and study their dependence on fitness as an artificial agent ("animat") evolves over thousands of generations to solve a navigation task in a simple, simulated environment. We compare the ability of these measures to predict high fitness with more conventional information-theoretic processing measures. As the animat adapts by increasing its "fit" to the world, information integration and processing increase commensurately along the evolutionary line of descent. We suggest that the correlation of fitness with information integration and with processing measures implies that high fitness requires both information processing as well as integration, but that information integration may be a better measure when the task requires memory. A correlation of measures of information integration (but also information processing) and fitness strongly suggests that these measures reflect the functional complexity of the animat, and that such measures can be used to quantify functional complexity even in the absence of fitness data.Comment: 27 pages, 8 figures, one supplementary figure. Three supplementary video files available on request. Version commensurate with published text in PLoS Comput. Bio

The Catalogue for Astrophysical Turbulence Simulations (CATS)

Turbulence is a key process in many fields of astrophysics. Advances in numerical simulations of fluids over the last several decades have revolutionized our understanding of turbulence and related processes such as star formation and cosmic ray propagation. However, data from numerical simulations of astrophysical turbulence are often not made public. We introduce a new simulation-oriented database for the astronomical community: the Catalogue for Astrophysical Turbulence Simulations (CATS), located at www.mhdturbulence.com. CATS includes magnetohydrodynamic (MHD) turbulent box simulation data products generated by the public codes athena++, arepo, enzo, and flash. CATS also includes several synthetic observational data sets, such as turbulent HI data cubes. We also include measured power spectra and three-point correlation functions from some of these data. We discuss the importance of open-source statistical and visualization tools for the analysis of turbulence simulations such as those found in CATS.D.C. acknowledges compute resources provided by NSF TRAC allocations TG-AST090110, TG-MCA07S014, and TG-AST140008. C.F. acknowledges funding provided by the Australian Research Council (Discovery Project DP170100603 and Future Fellowship FT180100495), the Australia-Germany Joint Research Cooperation Scheme (UA-DAAD), and high-performance computing resources provided by the Leibniz Rechenzentrum, the Gauss Centre for Supercomputing (grants pr32lo), and the Australian National Computational Infrastructure (grant ek9) in the framework of the National Computational Merit Allocation Scheme and the ANU Merit Allocation Scheme.and from the National Computational Infrastructure (NCI), which is supported by the Australian Government (award jh2). A.L. acknowledges the support of the NSF AST 1816234, NASA TCAN 144AAG1967, NASA ATP AAH7546 and the Flatiron Institute. J.C.'s work is supported by the National R&D Program through the National Research Foundation of Korea Grants funded by the Korean Government (NRF2016R1A5A1013277 and NRF-2016R1D1A1B02015014). P.M. acknowledges support for this work provided by NASA through Einstein Postdoctoral Fellowship grant number PF7-180164 awarded by the Chandra X-ray Center, which is operated by the Smithsonian Astrophysical Observatory for NASA under contract NAS8-0306

A Flow Induced Autoimmune Response and Accelerated Senescence of Red Blood Cells in Cardiovascular Devices

Red blood cells (RBCs) passing through heart pumps, prosthetic heart valves and other cardiovascular devices undergo early senescence attributed to non-physiologic forces. We hypothesized that mechanical trauma accelerates aging by deformation of membrane proteins to cause binding of naturally occurring IgG. RBCs isolated from blood of healthy volunteers were exposed to high shear stress in a viscometer or microfluidics channel to mimic mechanical trauma and then incubated with autologous plasma. Increased binding of IgG was observed indicating forces caused conformational changes in a membrane protein exposing an epitope(s), probably the senescent cell antigen of band 3. The binding of immunoglobulin suggests it plays a role in the premature sequestration and phagocytosis of RBCs in the spleen. Measurement of IgG holds promise as a marker foreshadowing complications in cardiovascular patients and as a means to improve the design of medical devices in which RBCs are susceptible to sublethal trauma.Research in this publication was supported by the National Institutes of Health Small Business Innovation Research program under award number R44HL114246 as a subcontract to the University of Oklahoma from VADovations and NIH grant R21HL132286 to DWS and TAS. Open Access fees paid for in whole or in part by the University of Oklahoma Libraries.Ye