Photon Stars

We discuss numerical solutions of Einstein's field equation describing static, spherically symmetric conglomerations of a photon gas. These equations imply a back reaction of the metric on the energy density of the photon gas according to Tolman's equation. The 3-fold of solutions corresponds to a class of physically different solutions which is parameterized by only two quantities, e.g. mass and surface temperature. The energy density is typically concentrated on a shell because the center contains a repelling singularity, which can, however, not be reached by timelike or null geodesics. The physical relevance of these solutions is completely open, although their existence may raise some doubts w.r. to the stability of black holes.Comment: 10 pages, 5 figures, talk at the DPG Spring Meeting 199

A Lagrangian model for the evolution of turbulent magnetic and passive scalar fields

In this paper we present an extension of the \emph{Recent Fluid Deformation (RFD)} closure introduced by Chevillard and Meneveau (2006) which was developed for modeling the time evolution of Lagrangian fluctuations in incompressible Navier-Stokes turbulence. We apply the RFD closure to study the evolution of magnetic and passive scalar fluctuations. This comparison is especially interesting since the stretching term for the magnetic field and for the gradient of the passive scalar are similar but differ by a sign such that the effect of stretching and compression by the turbulent velocity field is reversed. Probability density functions (PDFs) of magnetic fluctuations and fluctuations of the gradient of the passive scalar obtained from the RFD closure are compared against PDFs obtained from direct numerical simulations

Lagrangian Statistics of Navier-Stokes- and MHD-Turbulence

We report on a comparison of high-resolution numerical simulations of Lagrangian particles advected by incompressible turbulent hydro- and magnetohydrodynamic (MHD) flows. Numerical simulations were performed with up to $1024^3$ collocation points and 10 million particles in the Navier-Stokes case and $512^3$ collocation points and 1 million particles in the MHD case. In the hydrodynamics case our findings compare with recent experiments from Mordant et al. [1] and Xu et al. [2]. They differ from the simulations of Biferale et al. [3] due to differences of the ranges choosen for evaluating the structure functions. In Navier-Stokes turbulence intermittency is stronger than predicted by a multifractal approach of [3] whereas in MHD turbulence the predictions from the multifractal approach are more intermittent than observed in our simulations. In addition, our simulations reveal that Lagrangian Navier-Stokes turbulence is more intermittent than MHD turbulence, whereas the situation is reversed in the Eulerian case. Those findings can not consistently be described by the multifractal modeling. The crucial point is that the geometry of the dissipative structures have different implications for Lagrangian and Eulerian intermittency. Application of the multifractal approach for the modeling of the acceleration PDFs works well for the Navier-Stokes case but in the MHD case just the tails are well described.Comment: to appear in J. Plasma Phy

Lagrangian statistics in forced two-dimensional turbulence

We report on simulations of two-dimensional turbulence in the inverse energy cascade regime. Focusing on the statistics of Lagrangian tracer particles, scaling behavior of the probability density functions of velocity fluctuations is investigated. The results are compared to the three-dimensional case. In particular an analysis in terms of compensated cumulants reveals the transition from a strong non-Gaussian behavior with large tails to Gaussianity. The reported computation of correlation functions for the acceleration components sheds light on the underlying dynamics of the tracer particles.Comment: 8 figures, 1 tabl

Effect of turbulence on collisions of dust particles with planetesimals in protoplanetary disks

Planetesimals in gaseous protoplanetary disks may grow by collecting dust particles. Hydrodynamical studies show that small particles generally avoid collisions with the planetesimals because they are entrained by the flow around them. This occurs when $St$, the Stokes number, defined as the ratio of the dust stopping time to the planetesimal crossing time, becomes much smaller than unity. However, these studies have been limited to the laminar case, whereas these disks are believed to be turbulent. We want to estimate the influence of gas turbulence on the dust-planetesimal collision rate and on the impact speeds. We used three-dimensional direct numerical simulations of a fixed sphere (planetesimal) facing a laminar and turbulent flow seeded with small inertial particles (dust) subject to a Stokes drag. A no-slip boundary condition on the planetesimal surface is modeled via a penalty method. We find that turbulence can significantly increase the collision rate of dust particles with planetesimals. For a high turbulence case (when the amplitude of turbulent fluctuations is similar to the headwind velocity), we find that the collision probability remains equal to the geometrical rate or even higher for $St\geq 0.1$, i.e., for dust sizes an order of magnitude smaller than in the laminar case. We derive expressions to calculate impact probabilities as a function of dust and planetesimal size and turbulent intensity

Clustering of passive impurities in MHD turbulence

The transport of heavy, neutral or charged, point-like particles by incompressible, resistive magnetohydrodynamic (MHD) turbulence is investigated by means of high-resolution numerical simulations. The spatial distribution of such impurities is observed to display strong deviations from homogeneity, both at dissipative and inertial range scales. Neutral particles tend to cluster in the vicinity of coherent vortex sheets due to their viscous drag with the flow, leading to the simultaneous presence of very concentrated and almost empty regions. The signature of clustering is different for charged particles. These exhibit in addition to the drag the Lorentz-force. The regions of spatial inhomogeneities change due to attractive and repulsive vortex sheets. While small charges increase clustering, larger charges have a reverse effect.Comment: 9 pages, 13 figure

Dispersion Relations for Thermally Excited Waves in Plasma Crystals

Thermally excited waves in a Plasma crystal were numerically simulated using a Box_Tree code. The code is a Barnes_Hut tree code proven effective in modeling systems composed of large numbers of particles. Interaction between individual particles was assumed to conform to a Yukawa potential. Particle charge, mass, density, Debye length and output data intervals are all adjustable parameters in the code. Employing a Fourier transform on the output data, dispersion relations for both longitudinal and transverse wave modes were determined. These were compared with the dispersion relations obtained from experiment as well as a theory based on a harmonic approximation to the potential. They were found to agree over a range of 0.9<k<5, where k is the shielding parameter, defined by the ratio between interparticle distance a and dust Debye length lD. This is an improvement over experimental data as current experiments can only verify the theory up to k = 1.5.Comment: 8 pages, Presented at COSPAR '0

Automatic Neuron Detection in Calcium Imaging Data Using Convolutional Networks

Calcium imaging is an important technique for monitoring the activity of thousands of neurons simultaneously. As calcium imaging datasets grow in size, automated detection of individual neurons is becoming important. Here we apply a supervised learning approach to this problem and show that convolutional networks can achieve near-human accuracy and superhuman speed. Accuracy is superior to the popular PCA/ICA method based on precision and recall relative to ground truth annotation by a human expert. These results suggest that convolutional networks are an efficient and flexible tool for the analysis of large-scale calcium imaging data.Comment: 9 pages, 5 figures, 2 ancillary files; minor changes for camera-ready version. appears in Advances in Neural Information Processing Systems 29 (NIPS 2016

Transformation of stimulus correlations by the retina

Redundancies and correlations in the responses of sensory neurons seem to waste neural resources but can carry cues about structured stimuli and may help the brain to correct for response errors. To assess how the retina negotiates this tradeoff, we measured simultaneous responses from populations of ganglion cells presented with natural and artificial stimuli that varied greatly in correlation structure. We found that pairwise correlations in the retinal output remained similar across stimuli with widely different spatio-temporal correlations including white noise and natural movies. Meanwhile, purely spatial correlations tended to increase correlations in the retinal response. Responding to more correlated stimuli, ganglion cells had faster temporal kernels and tended to have stronger surrounds. These properties of individual cells, along with gain changes that opposed changes in effective contrast at the ganglion cell input, largely explained the similarity of pairwise correlations across stimuli where receptive field measurements were possible.Comment: author list corrected in metadat