25 Years of Self-Organized Criticality: Numerical Detection Methods

The detection and characterization of self-organized criticality (SOC), in both real and simulated data, has undergone many significant revisions over the past 25 years. The explosive advances in the many numerical methods available for detecting, discriminating, and ultimately testing, SOC have played a critical role in developing our understanding of how systems experience and exhibit SOC. In this article, methods of detecting SOC are reviewed; from correlations to complexity to critical quantities. A description of the basic autocorrelation method leads into a detailed analysis of application-oriented methods developed in the last 25 years. In the second half of this manuscript space-based, time-based and spatial-temporal methods are reviewed and the prevalence of power laws in nature is described, with an emphasis on event detection and characterization. The search for numerical methods to clearly and unambiguously detect SOC in data often leads us outside the comfort zone of our own disciplines - the answers to these questions are often obtained by studying the advances made in other fields of study. In addition, numerical detection methods often provide the optimum link between simulations and experiments in scientific research. We seek to explore this boundary where the rubber meets the road, to review this expanding field of research of numerical detection of SOC systems over the past 25 years, and to iterate forwards so as to provide some foresight and guidance into developing breakthroughs in this subject over the next quarter of a century.Comment: Space Science Review series on SO

Three-dimensional random Voronoi tessellations: From cubic crystal lattices to Poisson point processes

We perturb the SC, BCC, and FCC crystal structures with a spatial Gaussian noise whose adimensional strength is controlled by the parameter a, and analyze the topological and metrical properties of the resulting Voronoi Tessellations (VT). The topological properties of the VT of the SC and FCC crystals are unstable with respect to the introduction of noise, because the corresponding polyhedra are geometrically degenerate, whereas the tessellation of the BCC crystal is topologically stable even against noise of small but finite intensity. For weak noise, the mean area of the perturbed BCC and FCC crystals VT increases quadratically with a. In the case of perturbed SCC crystals, there is an optimal amount of noise that minimizes the mean area of the cells. Already for a moderate noise (a>0.5), the properties of the three perturbed VT are indistinguishable, and for intense noise (a>2), results converge to the Poisson-VT limit. Notably, 2-parameter gamma distributions are an excellent model for the empirical of of all considered properties. The VT of the perturbed BCC and FCC structures are local maxima for the isoperimetric quotient, which measures the degre of sphericity of the cells, among space filling VT. In the BCC case, this suggests a weaker form of the recentluy disproved Kelvin conjecture. Due to the fluctuations of the shape of the cells, anomalous scalings with exponents >3/2 is observed between the area and the volumes of the cells, and, except for the FCC case, also for a->0. In the Poisson-VT limit, the exponent is about 1.67. As the number of faces is positively correlated with the sphericity of the cells, the anomalous scaling is heavily reduced when we perform powerlaw fits separately on cells with a specific number of faces

Symmetry-break in Voronoi tessellations

We analyse in a common framework the properties of the Voronoi tessellations resulting from regular 2D and 3D crystals and those of tessellations generated by Poisson distributions of points, thus joining on symmetry breaking processes and the approach to uniform random distributions of seeds. We perturb crystalline structures in 2D and 3D with a spatial Gaussian noise whose adimensional strength is α and analyse the statistical properties of the cells of the resulting Voronoi tessellations using an ensemble approach. In 2D we consider triangular, square and hexagonal regular lattices, resulting into hexagonal, square and triangular tessellations, respectively. In 3D we consider the simple cubic (SC), body-centred cubic (BCC), and face-centred cubic (FCC) crystals, whose corresponding Voronoi cells are the cube, the truncated octahedron, and the rhombic dodecahedron, respectively. In 2D, for all values α>0, hexagons constitute the most common class of cells. Noise destroys the triangular and square tessellations, which are structurally unstable, as their topological properties are discontinuous in α=0. On the contrary, the honeycomb hexagonal tessellation is topologically stable and, experimentally, all Voronoi cells are hexagonal for small but finite noise with α0.5), memory of the specific initial unperturbed state is lost, because the statistical properties of the three perturbed regular tessellations are indistinguishable. When α>2, results converge to those of Poisson-Voronoi tessellations. In 2D, while the isoperimetric ratio increases with noise for the perturbed hexagonal tessellation, for the perturbed triangular and square tessellations it is optimised for specific value of noise intensity. The same applies in 3D, where noise degrades the isoperimetric ratio for perturbed FCC and BCC lattices, whereas the opposite holds for perturbed SCC lattices. This allows for formulating a weaker form of the Kelvin conjecture. By analysing jointly the statistical properties of the area and of the volume of the cells, we discover that also the cells shape heavily fluctuates when noise is introduced in the system. In 2D, the geometrical properties of n-sided cells change with α until the Poisson-Voronoi limit is reached for α>2; in this limit the Desch law for perimeters is shown to be not valid and a square root dependence on n is established, which agrees with exact asymptotic results. Anomalous scaling relations are observed between the perimeter and the area in the 2D and between the areas and the volumes of the cells in 3D: except for the hexagonal (2D) and FCC structure (3D), this applies also for infinitesimal noise. In the Poisson-Voronoi limit, the anomalous exponent is about 0.17 in both the 2D and 3D case. A positive anomaly in the scaling indicates that large cells preferentially feature large isoperimetric quotients. As the number of faces is strongly correlated with the sphericity (cells with more faces are bulkier), in 3D it is shown that the anomalous scaling is heavily reduced when we perform power law fits separately on cells with a specific number of faces

Speckle from phase ordering systems

The statistical properties of coherent radiation scattered from phase-ordering materials are studied in detail using large-scale computer simulations and analytic arguments. Specifically, we consider a two-dimensional model with a nonconserved, scalar order parameter (Model A), quenched through an order-disorder transition into the two-phase regime. For such systems it is well established that the standard scaling hypothesis applies, consequently the average scattering intensity at wavevector _k and time t' is proportional to a scaling function which depends only on a rescaled time, t ~ |_k|^2 t'. We find that the simulated intensities are exponentially distributed, with the time-dependent average well approximated using a scaling function due to Ohta, Jasnow, and Kawasaki. Considering fluctuations around the average behavior, we find that the covariance of the scattering intensity for a single wavevector at two different times is proportional to a scaling function with natural variables mt = |t_1 - t_2| and pt = (t_1 + t_2)/2. In the asymptotic large-pt limit this scaling function depends only on z = mt / pt^(1/2). For small values of z, the scaling function is quadratic, corresponding to highly persistent behavior of the intensity fluctuations. We empirically establish a connection between the intensity covariance and the two-time, two-point correlation function of the order parameter. This connection allows sensitive testing, either experimental or numerical, of existing theories for two-time correlations in systems undergoing order-disorder phase transitions. Comparison between theory and our numerical results requires no adjustable parameters.Comment: 18 pgs RevTeX, to appear in PR

Probing vacuum birefringence by phase-contrast Fourier imaging under fields of high-intensity lasers

In vacuum high-intensity lasers can cause photon-photon interaction via the process of virtual vacuum polarization which may be measured by the phase velocity shift of photons across intense fields. In the optical frequency domain, the photon-photon interaction is polarization-mediated described by the Euler-Heisenberg effective action. This theory predicts the vacuum birefringence or polarization dependence of the phase velocity shift arising from nonlinear properties in quantum electrodynamics (QED). We suggest a method to measure the vacuum birefringence under intense optical laser fields based on the absolute phase velocity shift by phase-contrast Fourier imaging. The method may serve for observing effects even beyond the QED vacuum polarization.Comment: 14 pages, 9 figures. Accepted by Applied Physics

The Small World of Osteocytes: Connectomics of the Lacuno-Canalicular Network in Bone

Osteocytes and their cell processes reside in a large, interconnected network of voids pervading the mineralized bone matrix of most vertebrates. This osteocyte lacuno-canalicular network (OLCN) is believed to play important roles in mechanosensing, mineral homeostasis, and for the mechanical properties of bone. While the extracellular matrix structure of bone is extensively studied on ultrastructural and macroscopic scales, there is a lack of quantitative knowledge on how the cellular network is organized. Using a recently introduced imaging and quantification approach, we analyze the OLCN in different bone types from mouse and sheep that exhibit different degrees of structural organization not only of the cell network but also of the fibrous matrix deposited by the cells. We define a number of robust, quantitative measures that are derived from the theory of complex networks. These measures enable us to gain insights into how efficient the network is organized with regard to intercellular transport and communication. Our analysis shows that the cell network in regularly organized, slow-growing bone tissue from sheep is less connected, but more efficiently organized compared to irregular and fast-growing bone tissue from mice. On the level of statistical topological properties (edges per node, edge length and degree distribution), both network types are indistinguishable, highlighting that despite pronounced differences at the tissue level, the topological architecture of the osteocyte canalicular network at the subcellular level may be independent of species and bone type. Our results suggest a universal mechanism underlying the self-organization of individual cells into a large, interconnected network during bone formation and mineralization

Generation of spatially pure photon pairs in a multimode nonlinear waveguide using intermodal dispersion

We present experimental realization of type-II spontaneous parametric down-conversion in a periodically poled potassium titanyl phosphate (KTiOPO4) nonlinear waveguide. We demonstrate that by careful exploitation of intermodal dispersion in the waveguide it is feasible to produce photon pairs in well defined transverse modes without any additional spatial filtering at the output. Spatial characteristics is verified by measurements of the M2 beam quality factors. We also prepared a postselected polarization-entangled two-photon state shown to violate Bell's inequality. Similar techniques based on intermodal dispersion can be used to generate spatial entanglement and hyperentanglement.Comment: 11 pages, 5 figures, submitted to Proceedings of Quantum Communications and Quantum Imaging X Conference at SPIE 2012 Optics + Photonics, San Diego, 12-16 August 201