Optimal foraging strategies: Lévy walks balance searching and patch exploitation under a very broad range of conditions

While evidence for optimal random search patterns, known as Lévy walks, in empirical movement data is mounting for a growing list of taxa spanning motile cells to humans, there is still much debate concerning the theoretical generality of Lévy walk optimisation. Here, using a new and robust simulation environment, we investigate in the most detailed study to date (24×10(6) simulations) the foraging and search efficiencies of 2-D Lévy walks with a range of exponents, target resource distributions and several competing models. We find strong and comprehensive support for the predictions of the Lévy flight foraging hypothesis and in particular for the optimality of inverse square distributions of move step-lengths across a much broader range of resource densities and distributions than previously realised. Further support for the evolutionary advantage of Lévy walk movement patterns is provided by an investigation into the 'feast and famine' effect, with Lévy foragers in heterogeneous environments experiencing fewer long 'famines' than other types of searchers. Therefore overall, optimal Lévy foraging results in more predictable resources in unpredictable environments.

A new approach for objective identification of turns and steps in organism movement data relevant to random walk modelling

1. A first step in the analysis of complex movement data often involves discretisation of the path into a series of step-lengths and turns, for example in the analysis of specialised random walks, such as Lévy flights. However, the identification of turning points, and therefore step-lengths, in a tortuous path is dependent on ad-hoc parameter choices. Consequently, studies testing for movement patterns in these data, such as Lévy flights, have generated debate. However, studies focusing on one-dimensional (1D) data, as in the vertical displacements of marine pelagic predators, where turning points can be identified unambiguously have provided strong support for Lévy flight movement patterns. 2. Here, we investigate how step-length distributions in 3D movement patterns would be interpreted by tags recording in 1D (i.e. depth) and demonstrate the dimensional symmetry previously shown mathematically for Lévy-flight movements. We test the veracity of this symmetry by simulating several measurement errors common in empirical datasets and find Lévy patterns and exponents to be robust to low-quality movement data. 3. We then consider exponential and composite Brownian random walks and show that these also project into 1D with sufficient symmetry to be clearly identifiable as such. 4. By extending the symmetry paradigm, we propose a new methodology for step-length identification in 2D or 3D movement data. The methodology is successfully demonstrated in a re-analysis of wandering albatross Global Positioning System (GPS) location data previously analysed using a complex methodology to determine bird-landing locations as turning points in a Lévy walk. For this high-resolution GPS data, we show that there is strong evidence for albatross foraging patterns approximated by truncated Lévy flights spanning over 3·5 orders of magnitude. 5. Our simple methodology and freely available software can be used with any 2D or 3D movement data at any scale or resolution and are robust to common empirical measurement errors. The method should find wide applicability in the field of movement ecology spanning the study of motile cells to humans.

Scaling laws of ambush predator 'waiting' behaviour are tuned to a common ecology

The decisions animals make about how long to wait between activities can determine the success of diverse behaviours such as foraging, group formation or risk avoidance. Remarkably, for diverse animal species, including humans, spontaneous patterns of waiting times show random ‘burstiness’ that appears scale-invariant across a broad set of scales. However, a general theory linking this phenomenon across the animal kingdom currently lacks an ecological basis. Here, we demonstrate from tracking the activities of 15 sympatric predator species (cephalopods, sharks, skates and teleosts) under natural and controlled conditions that bursty waiting times are an intrinsic spontaneous behaviour well approximated by heavy-tailed (power-law) models over data ranges up to four orders of magnitude. Scaling exponents quantifying ratios of frequent short to rare very long waits are species-specific, being determined by traits such as foraging mode (active versus ambush predation), body size and prey preference. A stochastic–deterministic decision model reproduced the empirical waiting time scaling and species-specific exponents, indicating that apparently complex scaling can emerge from simple decisions. Results indicate temporal power-law scaling is a behavioural ‘rule of thumb’ that is tuned to species’ ecological traits, implying a common pattern may have naturally evolved that optimizes move–wait decisions in less predictable natural environments.

Semiclassical treatment of Wannier's theory when the exponent diverges

We consider a non-standard application of the Wannier model. A physical example is the single ionization of a hydrogenic beryllium ion with a fully stripped beryllium ion, where the ratio of the charge of the third particle to the charges of the escaping particles is 1/4; we investigate the single ionization by an electron of an atom comprising an electron and a nucleus of charge 1/4. An infinite exponent is obtained suggesting that this process is not tractable within the Wannier model. A modified version of Crothers' uniform semiclassical wavefunction for the outgoing particles has been adopted, since the Wannier exponents and are infinite for an effective charge of Z = 1/4. We use Bessel functions to describe the Peterkop functions u and u and derive a new turning point ?. Since u is well behaved at infinity, there exists only the singularity in u at infinity, thus we employ a one- (rather than two-) dimensional change of dependent variable, ensuring that a uniform solution is obtained that avoids semiclassical breakdown on the Wannier ridge. The regularized final-state asymptotic wavefunction is employed, along with a continuum-distorted-wave approximation for the initial-state wavefunction to obtain total cross sections on an absolute scale. © 2006 IOP Publishing Ltd.

Temperature-Driven Mixing-Demixing Behavior of Binary Mixtures of the Ionic Liquid Choline Bis(trifluoromethylsulfonyl)imide and Water

The ionic liquid (2-hydroxyethylammonium)trimethylammonium) bis(trifluoromethylsulfonyl)imide (choline bistriflimide) was obtained as a supercooled liquid at room temperature (melting point = 30 degrees C). Crystals of choline bistriflimide suitable for structure determination were grown from the melt in situ on the X-ray diffractometer. The choline cation adopts a folded conformation, whereas the bistriflimide anion exhibits a transoid conformation. The choline cation and the bistriflimide anion are held together by hydrogen bonds between the hydroxyl proton and a sulfonyl oxygen atom. This hydrogen bonding is of importance for the temperature-dependent solubility proper-ties of the ionic liquid. Choline bistriflimide is not miscible with water at room temperature, but forms one phase with water at temperatures above 72 degrees C (equals upper critical solution temperature). H-1 NMR studies show that the hydrogen bonds between the choline cation and the bistriflimide anion are substantially weakened above this temperature. The thermophysical properties of water-choline bistriflimide binary mixtures were furthermore studied by a photopyroelectric technique and by adiabatic scanning calorimetry (ASC). By photothermal analysis, besides highly accurate values for the thermal conductivity and effusivity of choline bistriflimide at 30 degrees C, the detailed temperature dependence of both the thermal conductivity and effusivity of the upper and lower part of a critical water-choline bistriflimide mixture in the neighborhood of the mixing-demixing phase transition could be determined with high resolution and accuracy. Together with high resolution ASC data for the heat capacity, experimental values were obtained for the critical exponents alpha and beta, and for the critical amplitude ratio G(+)/G(-). These three values were found to be consistent with theoretical expectations for a three dimensional Ising-type of critical behavior of binary liquid mixtures.

Universal power-law diet partitioning by marine fish and squid with surprising stability–diversity implications

A central question in community ecology is how the number of trophic links relates to community species richness. For simple dynamical food-web models, link density (the ratio of links to species) is bounded from above as the number of species increases; but empirical data suggest that it increases without bounds. We found a new empirical upper bound on link density in large marine communities with emphasis on fish and squid, using novel methods that avoid known sources of bias in traditional approaches. Bounds are expressed in terms of the diet-partitioning function (DPF): the average number of resources contributing more than a fraction f to a consumer's diet, as a function of f. All observed DPF follow a functional form closely related to a power law, with power-law exponents indepen- dent of species richness at the measurement accuracy. Results imply universal upper bounds on link density across the oceans. However, the inherently scale-free nature of power-law diet partitioning suggests that the DPF itself is a better defined characterization of network structure than link density.

NON-GAUSSIAN STATISTICS OF OIL PRICING TIME-SERIES: A CASE STUDY

The stochastic nature of oil price fluctuations is investigated over a twelve-year period, borrowing feedback from an existing database (USA Energy Information Administration database, available online). We evaluate the scaling exponents of the fluctuations by employing different statistical analysis methods, namely rescaled range analysis (R/S), scale windowed variance analysis (SWV) and the generalized Hurst exponent (GH) method. Relying on the scaling exponents obtained, we apply a rescaling procedure to investigate the complex characteristics of the probability density functions (PDFs) dominating oil price fluctuations. It is found that PDFs exhibit scale invariance, and in fact collapse onto a single curve when increments are measured over microscales (typically less than 30 days). The time evolution of the distributions is well fitted by a Levy-type stable distribution. The relevance of a Levy distribution is made plausible by a simple model of nonlinear transfer. Our results also exhibit a degree of multifractality as the PDFs change and converge toward to a Gaussian distribution at the macroscales.

Random neighbor sandpile model with constrained total energy

We propose as energy-constrained sandpile model with random neighbors. The critical behavior of the model is in the same universality class as the mean-field self-organized criticality sandpile. The critical energy E-c depends on the number of neighbors n of each site, but the various exponents do not. For n = 6, we got that E-c = 0.4545; and a self-similar structure of the energy distribution function with five major peaks is also observed. This is a natural result of system dynamics and the way the system is disturbed.

Universality and self-similarity of an energy-constrained sandpile model with random neighbors

We study an energy-constrained sandpile model with random neighbors. The critical behavior of the model is in the same universality class as the mean-field self-organized criticality sandpile. The critical energy E-c depends on the number of neighbors n for each site, but the various exponents are independent of n. A self-similar structure with n-1 major peaks is developed for the energy distribution p(E) when the system approaches its stationary state. The avalanche dynamics contributes to the major peaks appearing at E-Pk = 2k/(2n - 1) with k = 1,2,...,n-1, while the fine self-similar structure is a natural result of the way the system is disturbed. [S1063-651X(99)10307-6].

On the universality of a one-dimensional model of a rice pile

A one-dimensional model of a rice pile is numerically studied for different driving mechanisms and different levels of medium disorder. The universality of the scaling exponents for the transit time distribution and avalanche size distribution is discussed. (C) 1997 Published by Elsevier Science B.V.

ON THE UNIVERSALITY OF 2-DIMENSIONAL DIFFUSION-LIMITED DEPOSITION

A history dependent stick probability is introduced to the diffusion-limited deposition model. The exponents in the scaling laws are calculated. The universality class is also discussed.

STO and GTO field-induced polarization functions for H to Kr

Field-induced polarization (FIP) functions were proposed over two decades ago to improve the accuracy of calculated response properties, and the FIP functions in GTO form for H and C to F were tested on small molecules, with encouraging results. The concept of FIP,is now extended to all atoms up to Kr. New simplifying approximations for the description of asymptotic highest occupied atomic orbitals. (HOAOs) are introduced in this study. They provide the basis for STO and GTO exponents of a complete set of FIP functions from H to Kr, which are both listed for the convenience of the users. Tests on the polarizabilities of a series of atoms and molecules demonstrate that addition of the FIP basis functions to a series' of standard basis sets drastically improves the performance of all these basis sets compared to converged results. Moreover, the byproduct of this study (approximate asymptotic HOAOs) provides information for the construction of accurate basis sets for long-range ground state properties. (C) 2003 Wiley Periodicals, Inc.

Testing metabolic scaling theory using intraspecific allometries in Antarctic microarthropods

Quantitative scaling relationships among body mass, temperature and metabolic rate of organisms are still controversial, while resolution may be further complicated through the use of different and possibly inappropriate approaches to statistical analysis. We propose the application of a modelling strategy based on the theoretical approach of Akaike's information criteria and non-linear model fitting (nlm). Accordingly, we collated and modelled available data at intraspecific level on the individual standard metabolic rate of Antarctic microarthropods as a function of body mass (M), temperature (T), species identity (S) and high rank taxa to which species belong (G) and tested predictions from metabolic scaling theory (mass-metabolism allometric exponent b = 0.75, activation energy range 0.2-1.2 eV). We also performed allometric analysis based on logarithmic transformations (lm). Conclusions from lm and nlm approaches were different. Best-supported models from lm incorporated T, M and S. The estimates of the allometric scaling exponent linking body mass and metabolic rate resulted in a value of 0.696 +/- 0.105 (mean +/- 95% CI). In contrast, the four best-supported nlm models suggested that both the scaling exponent and activation energy significantly vary across the high rank taxa (Collembola, Cryptostigmata, Mesostigmata and Prostigmata) to which species belong, with mean values of b ranging from about 0.6 to 0.8. We therefore reached two conclusions: 1, published analyses of arthropod metabolism based on logarithmic data may be biased by data transformation; 2, non-linear models applied to Antarctic microarthropod metabolic rate suggest that intraspecific scaling of standard metabolic rate in Antarctic microarthropods is highly variable and can be characterised by scaling exponents that greatly vary within taxa, which may have biased previous interspecific comparisons that neglected intraspecific variability.

Global quantum correlations in finite-size spin chains

We perform an extensive study of the properties of global quantum correlations in finite-size one-dimensional quantum spin models at finite temperature. By adopting a recently proposed measure for global quantum correlations (Rulli and Sarandy 2011 Phys. Rev. A 84 042109), called global discord, we show that critical points can be neatly detected even for many-body systems that are not in their ground state. We consider the transverse Ising model, the cluster-Ising model where three-body couplings compete with an Ising-like interaction, and the nearest-neighbor XX Hamiltonian in transverse magnetic field. These models embody our canonical examples showing the sensitivity of global quantum discord close to criticality. For the Ising model, we find a universal scaling of global discord with the critical exponents pertaining to the Ising universality class.

Scaling of the entanglement spectrum near quantum phase transitions

The entanglement spectrum describing quantum correlations in many-body systems has been recently recognized as a key tool to characterize different quantum phases, including topological ones. Here we derive its analytically scaling properties in the vicinity of some integrable quantum phase transitions and extend our studies also to nonintegrable quantum phase transitions in one-dimensional spin models numerically. Our analysis shows that, in all studied cases, the scaling of the difference between the two largest nondegenerate Schmidt eigenvalues yields with good accuracy critical points and mass scaling exponents.