Day and Night Changes of Cardiovascular Complexity: A Multi-Fractal Multi-Scale Analysis

Recently, a multifractal-multiscale approach to detrended fluctuation analysis (DFA) was proposed to evaluate the cardiovascular fractal dynamics providing a surface of self-similarity coefficients α(q,τ), function of the scale τ, and moment order q. We hypothesize that this versatile DFA approach may reflect the cardiocirculatory adaptations in complexity and nonlinearity occurring during the day/night cycle. Our aim is, therefore, to quantify how α(q, τ) surfaces of cardiovascular series differ between daytime and night-time. We estimated α(q,τ) with -5 ≤ q ≤ 5 and 8 ≤ τ ≤ 2048 s for heart rate and blood pressure beat-to-beat series over periods of few hours during daytime wake and night-time sleep in 14 healthy participants. From the α(q,τ) surfaces, we estimated short-term (<16 s) and long-term (from 16 to 512 s) multifractal coefficients. Generating phase-shuffled surrogate series, we evaluated short-term and long-term indices of nonlinearity for each q. We found a long-term night/day modulation of α(q,τ) between 128 and 256 s affecting heart rate and blood pressure similarly, and multifractal short-term modulations at q < 0 for the heart rate and at q > 0 for the blood pressure. Consistent nonlinearity appeared at the shorter scales at night excluding q = 2. Long-term circadian modulations of the heart rate DFA were previously associated with the cardiac vulnerability period and our results may improve the risk stratification indicating the more relevant α(q,τ) area reflecting this rhythm. Furthermore, nonlinear components in the nocturnal α(q,τ) at q ≠ 2 suggest that DFA may effectively integrate the linear spectral information with complexity-domain information, possibly improving the monitoring of cardiac interventions and protocols of rehabilitation medicine.

Bubble pinch-off in turbulence

Although bubble pinch-off is an archetype of a dynamical system evolving toward a singularity, it has always been described in idealized theoretical and experimental conditions. Here, we consider bubble pinch-off in a turbulent flow representative of natural conditions in the presence of strong and random perturbations, combining laboratory experiments, numerical simulations, and theoretical modeling. We show that the turbulence sets the initial conditions for pinch-off, namely the initial bubble shape and flow field, but after the pinch-off starts, the turbulent time at the neck scale becomes much slower than the pinching dynamics: The turbulence freezes. We show that the average neck size, [Formula: see text], can be described by [Formula: see text], where [Formula: see text] is the pinch-off or singularity time and [Formula: see text], in close agreement with the axisymmetric theory with no initial flow. While frozen, the turbulence can influence the pinch-off through the initial conditions. Neck shape oscillations described by a quasi-2-dimensional (quasi-2D) linear perturbation model are observed as are persistent eccentricities of the neck, which are related to the complex flow field induced by the deformed bubble shape. When turbulent stresses are less able to be counteracted by surface tension, a 3-dimensional (3D) kink-like structure develops in the neck, causing [Formula: see text] to escape its self-similar decrease. We identify the geometric controlling parameter that governs the appearance of these kink-like interfacial structures, which drive the collapse out of the self-similar route, governing both the likelihood of escaping the self-similar process and the time and length scale at which it occurs.

Tree crowns grow into self-similar shapes controlled by gravity and light sensing

Plants have developed different tropisms: in particular, they reorient the growth of their branches towards the light (phototropism) or upwards (gravitropism). How these tropisms affect the shape of a tree crown remains unanswered. We address this question by developing a propagating front model of tree growth. Being length-free, this model leads to self-similar solutions after a long period of time, which are independent of the initial conditions. Varying the intensities of each tropism, different self-similar shapes emerge, including singular ones. Interestingly, these shapes bear similarities to existing tree species. It is concluded that the core of specific crown shapes in trees relies on the balance between tropisms.

Modeling and Optimization of Time-Of-Use Electricity Pricing Systems

Time-of-Use (TOU) pricing is an important strategy for electricity providers to manage supply and hence making the grid more efficient and for consumers to manage their costs. In this paper, we discuss a general stochastic modeling framework for consumer's power demand based on which the TOU contract characteristics can be selected, so as to minimize the mean electricity price paid by the customer. We exploit the characteristics of power demand observed in real grids to propose to model it during homogeneous peak periods as a constant level with fluctuations described by a scaled fractional Brownian motion. We analyze the exceedance process over pre-specified thresholds and use this information for formulating an optimization problem to determine the key features of the TOU contract. Due to the analytical intractability of certain expressions with the exception of short-range dependence fluctuations, the solution of the posited optimization problem requires using techniques such as Monte Carlo simulation and numerical search. The methodology for two pricing schemes is illustrated using real data.

Non-equilibrium turbulence scalings and self-similarity in turbulent planar jets

We study the self-similarity and dissipation scalings of a turbulent planar jet and the theoretically implied mean flow scalings. Unlike turbulent wakes where such studies have already been carried out (Dairay et al. 2015 J. Fluid Mech. 781, 166-198. (doi:10.1017/jfm.2015.493); Obligado et al. 2016 Phys. Rev. Fluids 1, 044409. (doi:10.1103/PhysRevFluids.1.044409)), this is a boundary-free turbulent shear flow where the local Reynolds number increases with distance from inlet. The Townsend-George theory revised by (Dairay et al. 2015 J. Fluid Mech. 781, 166-198. (doi:10.1017/jfm.2015.493)) is applied to turbulent planar jets. Only a few profiles need to be self-similar in this theory. The self-similarity of mean flow, turbulence dissipation, turbulent kinetic energy and Reynolds stress profiles is supported by our experimental results from 18 to at least 54 nozzle sizes, the furthermost location investigated in this work. Furthermore, the non-equilibrium dissipation scaling found in turbulent wakes, decaying grid-generated turbulence, various instances of periodic turbulence and turbulent boundary layers (Dairay et al. 2015 J. Fluid Mech. 781, 166-198. (doi:10.1017/jfm.2015.493); Vassilicos 2015 Annu. Rev. Fluid Mech. 95, 114. (doi:10.1146/annurev-fluid-010814-014637); Goto & Vassilicos 2015 Phys. Lett. A 3790, 1144-1148. (doi:10.1016/j.physleta.2015.02.025); Nedic et al. 2017 Phys. Rev. Fluids 2, 032601. (doi:10.1103/PhysRevFluids.2.032601)) is also observed in the present turbulent planar jet and in the turbulent planar jet of (Antonia et al. 1980 Phys. Fluids 23, 863055. (doi:10.1063/1.863055)). Given these observations, the theory implies new mean flow and jet width scalings which are found to be consistent with our data and the data of (Antonia et al. 1980 Phys. Fluids 23, 863055. (doi:10.1063/1.863055)). In particular, it implies a hitherto unknown entrainment behaviour: the ratio of characteristic cross-stream to centreline streamwise mean flow velocities decays as the -1/3 power of streamwise distance in the region, where the non-equilibrium dissipation scaling holds.

Trait Mindfulness Is Associated With the Self-Similarity of Heart Rate Variability

Previous studies have linked trait mindfulness with better self-regulation and adaptation. Heart rate variability (HRV) is a good physiological indicator of the capacity for self-regulation and adaptation. The present study explored the relationship between trait mindfulness and HRV from the viewpoint of crosstalking between different HRV parameter pairs, which would reflect the dynamic interactions between each pair of HRV parameters in different processes. We measured the trait mindfulness of seventy-four undergraduate students and recorded nine HRV parameters during the following four consecutive experimental phases: (1) calming phase, (2) mental arithmetic task phase, (3) recovery phase, and (4) mindfulness practice phase. The relationship between trait mindfulness and HRV was explored at the following three levels: (1) the absolute level, i.e., HRV parameters in four different states, (2) the difference-change level, i.e., differences in HRV parameters between different states, and (3) the crosstalking level, i.e., self-similarity of crosstalking HRV parameter pairs. The results supported the following hypothesis: trait mindfulness, as measured by the Mindful Attention Awareness Scale (MAAS), was significantly and positively correlated with the self-similarity of crosstalking HRV parameter pairs but was not significantly correlated with the HRV parameters at the difference-change and absolute levels. These findings indicate that as trait mindfulness increases, the ability to maintain ANS function homeostasis improves. HIGHLIGHTS -Trait mindfulness is associated with better self-regulation and adaptation.-Heart rate variability (HRV) is a good physiological indicator of the capacity for self-regulation and adaptation.-Trait mindfulness is significantly correlated with self-similarity of crosstalking HRV parameter pairs but not with the HRV parameters at the difference-change or absolute levels.

Assessment of Driver's Drowsiness Based on Fractal Dimensional Analysis of Sitting and Back Pressure Measurements

The most effective way of preventing motor vehicle accidents caused by drowsy driving is through a better understanding of drowsiness itself. Prior research on the detection of symptoms of drowsy driving has offered insights on providing drivers with advance warning of an elevated risk of crash. The present study measured back and sitting pressures during a simulated driving task under both high and low arousal conditions. Fluctuation of time series of center of pressure (COP) movement of back and sitting pressure was observed to possess a fractal property. The fractal dimensions were calculated to compare the high and low arousal conditions. The results showed that under low arousal (the drowsiness state) the fractal dimension was significantly lower than what was calculated with high arousal. Accumulated drowsiness thus contributed to the loss of self-similarity and unpredictability of time series of back and sitting pressure measurement. Drowsiness further reduces the complexity of the posture control system as viewed from back and sitting pressure. Thus, fractal dimension is a necessary and sufficient condition of a decreased arousal level. It further is a necessary condition for detecting the interval or point in time with high risk of crash.

Oscillating path between self-similarities in liquid pinch-off

Many differential equations involved in natural sciences show singular behaviors; i.e., quantities in the model diverge as the solution goes to zero. Nonetheless, the evolution of the singularity can be captured with self-similar solutions, several of which may exist for a given system. How to characterize the transition from one self-similar regime to another remains an open question. By studying the classic example of the pinch-off of a viscous liquid thread, we show experimentally that the geometry of the system and external perturbations play an essential role in the transition from a symmetric to an asymmetric solution. Moreover, this transient regime undergoes unexpected log-scale oscillations that delay dramatically the onset of the final self-similar solution. This result sheds light on the strong impact external constraints can have on predictions established to explain the formation of satellite droplets or on the rheological tests applied on a fluid, for example.

Symmetric coalescence of two hydraulic fractures

The formation of a fracture network is a key process for many geophysical and industrial practices from energy resource recovery to induced seismic management. We focus on the initial stage of a fracture network formation using experiments on the symmetric coalescence of two equal coplanar, fluid-driven, penny-shaped fractures in a brittle elastic medium. Initially, the fractures propagate independently of each other. The fractures then begin to interact and coalesce, forming a bridge between them. Within an intermediate period after the initial contact, most of the fracture growth is localized along this bridge, perpendicular to the line connecting the injection sources. Using light attenuation and particle image velocimetry to measure both the fracture aperture and velocity field, we characterize the growth of this bridge. We model this behavior using a geometric volume conservation argument dependent on the symmetry of the interaction, with a 2D approximation for the bridge. We also verify experimentally the scaling for the bridge growth and the shape of the thickness profile along the bridge. The influence of elasticity and toughness of the solid, injection rate of the fluid, and initial location of the fractures are captured by our scaling.

Fractal Structure and Non-Extensive Statistics

The role played by non-extensive thermodynamics in physical systems has been under intense debate for the last decades. With many applications in several areas, the Tsallis statistics have been discussed in detail in many works and triggered an interesting discussion on the most deep meaning of entropy and its role in complex systems. Some possible mechanisms that could give rise to non-extensive statistics have been formulated over the last several years, in particular a fractal structure in thermodynamic functions was recently proposed as a possible origin for non-extensive statistics in physical systems. In the present work, we investigate the properties of such fractal thermodynamical system and propose a diagrammatic method for calculations of relevant quantities related to such a system. It is shown that a system with the fractal structure described here presents temperature fluctuation following an Euler Gamma Function, in accordance with previous works that provided evidence of the connections between those fluctuations and Tsallis statistics. Finally, the scale invariance of the fractal thermodynamical system is discussed in terms of the Callan–Symanzik equation.

Assessing and Interpreting the Within-Body Biogeography of Human Microbiome Diversity

A human body hosts a relatively independent microbiome including five major regional biomes (i.e., airway, oral, gut, skin, and urogenital). Each of them may possess different regional characteristics with important implications to our health and diseases (i.e., so-termed microbiome associated diseases). Nevertheless, these regional microbiomes are connected with each other through diffusions and migrations. Here, we investigate the within-body (intra-individual) distribution feature of microbiome diversity via diversity area relationship (DAR) modeling, which, to the best of our knowledge, has not been systematically studied previously. We utilized the Hill numbers for measuring alpha and beta-diversities and built 1,200 within-body DAR models with to date the most comprehensive human microbiome datasets of 18 sites from the human microbiome project (HMP) cohort. We established the intra-DAR profile (z-q pattern: the diversity scaling parameter z of the power law (PL) at diversity order q = 0–3), intra-PDO (pair-wise diversity overlap) profile (g-q), and intra-MAD (maximal accrual diversity) profile (D_max-q) for the within-body biogeography of the human microbiome. These profiles constitute the “maps” of the within-body biogeography, and offer important insights on the within-body distribution of the human microbiome. Furthermore, we investigated the heterogeneity among individuals in their biogeography parameters and found that there is not an “average Joe” that can represent majority of individuals in a cohort or population. For example, we found that most individuals in the HMP cohort have relatively lower maximal accrual diversity (MAD) or in the “long tail” of the so-termed power law distribution. In the meantime, there are a small number of individuals in the cohort who possess disproportionally higher MAD values. These findings may have important implications for personalized medicine of the human microbiome associated diseases in practice, besides their theoretical significance in microbiome research such as establishing the baseline for the conservation of human microbiome.

Branch Thinning and the Large-Scale, Self-Similar Structure of Trees

Branch formation in trees has an inherent tendency toward exponential growth, but exponential growth in the number of branches cannot continue indefinitely. It has been suggested that trees balance this tendency toward expansion by also losing branches grown in previous growth cycles. Here, we present a model for branch formation and branch loss during ontogeny that builds on the phenomenological assumption of a branch carrying capacity. The model allows us to derive approximate analytical expressions for the number of tips on a branch, the distribution of growth modules within a branch, and the rate and size distribution of tree wood litter produced. Although limited availability of data makes empirical corroboration challenging, we show that our model can fit field observations of red maple (Acer rubrum) and note that the age distribution of discarded branches predicted by our model is qualitatively similar to an empirically observed distribution of dead and abscised branches of balsam poplar (Populus balsamifera). By showing how a simple phenomenological assumption-that the number of branches a tree can maintain is limited-leads directly to predictions on branching structure and the rate and size distribution of branch loss, these results potentially enable more explicit modeling of woody tissues in ecosystems worldwide, with implications for the buildup of flammable fuel, nutrient cycling, and understanding of plant growth.

The Fractality of Polar and Reed-Muller Codes

The generator matrices of polar codes and Reed-Muller codes are submatrices of the Kronecker product of a lower-triangular binary square matrix. For polar codes, the submatrix is generated by selecting rows according to their Bhattacharyya parameter, which is related to the error probability of sequential decoding. For Reed-Muller codes, the submatrix is generated by selecting rows according to their Hamming weight. In this work, we investigate the properties of the index sets selecting those rows, in the limit as the blocklength tends to infinity. We compute the Lebesgue measure and the Hausdorff dimension of these sets. We furthermore show that these sets are finely structured and self-similar in a well-defined sense, i.e., they have properties that are common to fractals.

Statistical Image Properties in Large Subsets of Traditional Art, Bad Art, and Abstract Art

Several statistical image properties have been associated with large subsets of traditional visual artworks. Here, we investigate some of these properties in three categories of art that differ in artistic claim and prestige: (1) Traditional art of different cultural origin from established museums and art collections (oil paintings and graphic art of Western provenance, Islamic book illustration and Chinese paintings), (2) Bad Art from two museums that collect contemporary artworks of lesser importance (© Museum Of Bad Art [MOBA], Somerville, and Official Bad Art Museum of Art [OBAMA], Seattle), and (3) twentieth century abstract art of Western provenance from two prestigious museums (Tate Gallery and Kunstsammlung Nordrhein-Westfalen). We measured the following four statistical image properties: the fractal dimension (a measure relating to subjective complexity); self-similarity (a measure of how much the sections of an image resemble the image as a whole), 1st-order entropy of edge orientations (a measure of how uniformly different orientations are represented in an image); and 2nd-order entropy of edge orientations (a measure of how independent edge orientations are across an image). As shown previously, traditional artworks of different styles share similar values for these measures. The values for Bad Art and twentieth century abstract art show a considerable overlap with those of traditional art, but we also identified numerous examples of Bad Art and abstract art that deviate from traditional art. By measuring statistical image properties, we quantify such differences in image composition for the first time.

A Statistical Description of Plant Shoot Architecture

Plant architectures can be characterized statistically by their spatial density function, which specifies the probability of finding a branch at each location in the territory occupied by a plant. Using high-precision 3D scanning, we analyzed 557 plant shoot architectures, representing three species, grown across three to five environmental conditions, and through 20-30 developmental time points. We found two elegant properties in the spatial density functions of these architectures: all functions could be nearly modified in one direction without affecting the density in orthogonal directions (called "separability"), and all functions shared the same underlying shape, aside from stretching and compression (called "self-similarity"). Surprisingly, despite their striking visual diversity, we discovered that all architectures could be described as variations on a single underlying function: a Gaussian density function truncated at roughly two SDs. We also observed systematic variation in the spatial density functions across species, growth conditions, and time, which suggests functional specialization despite following the same general design form.

Using CNN Features to Better Understand What Makes Visual Artworks Special

One of the goal of computational aesthetics is to understand what is special about visual artworks. By analyzing image statistics, contemporary methods in computer vision enable researchers to identify properties that distinguish artworks from other (non-art) types of images. Such knowledge will eventually allow inferences with regard to the possible neural mechanisms that underlie aesthetic perception in the human visual system. In the present study, we define measures that capture variances of features of a well-established Convolutional Neural Network (CNN), which was trained on millions of images to recognize objects. Using an image dataset that represents traditional Western, Islamic and Chinese art, as well as various types of non-art images, we show that we need only two variance measures to distinguish between the artworks and non-art images with a high classification accuracy of 93.0%. Results for the first variance measure imply that, in the artworks, the subregions of an image tend to be filled with pictorial elements, to which many diverse CNN features respond (richness of feature responses). Results for the second measure imply that this diversity is tied to a relatively large variability of the responses of individual CNN feature across the subregions of an image. We hypothesize that this combination of richness and variability of CNN feature responses is one of properties that makes traditional visual artworks special. We discuss the possible neural underpinnings of this perceptual quality of artworks and propose to study the same quality also in other types of aesthetic stimuli, such as music and literature.

Statistical evidence of anasymptotic geometric structure to the momentum transporting motions in turbulent boundary layers

The turbulence contribution to the mean flow is reflected by the motions producing the Reynolds shear stress (〈-uv〉) and its gradient. Recent analyses of the mean dynamical equation, along with data, evidence that these motions asymptotically exhibit self-similar geometric properties. This study discerns additional properties associated with the uv signal, with an emphasis on the magnitudes and length scales of its negative contributions. The signals analysed derive from high-resolution multi-wire hot-wire sensor data acquired in flat-plate turbulent boundary layers. Space-filling properties of the present signals are shown to reinforce previous observations, while the skewness of uv suggests a connection between the size and magnitude of the negative excursions on the inertial domain. Here, the size and length scales of the negative uv motions are shown to increase with distance from the wall, whereas their occurrences decrease. A joint analysis of the signal magnitudes and their corresponding lengths reveals that the length scales that contribute most to 〈-uv〉 are distinctly larger than the average geometric size of the negative uv motions. Co-spectra of the streamwise and wall-normal velocities, however, are shown to exhibit invariance across the inertial region when their wavelengths are normalized by the width distribution, W(y), of the scaling layer hierarchy, which renders the mean momentum equation invariant on the inertial domain.This article is part of the themed issue 'Toward the development of high-fidelity models of wall turbulence at large Reynolds number'.

On the shape of giant soap bubbles

We study the effect of gravity on giant soap bubbles and show that it becomes dominant above the critical size [Formula: see text], where [Formula: see text] is the mean thickness of the soap film and [Formula: see text] is the capillary length ([Formula: see text] stands for vapor-liquid surface tension, and [Formula: see text] stands for the liquid density). We first show experimentally that large soap bubbles do not retain a spherical shape but flatten when increasing their size. A theoretical model is then developed to account for this effect, predicting the shape based on mechanical equilibrium. In stark contrast to liquid drops, we show that there is no mechanical limit of the height of giant bubble shapes. In practice, the physicochemical constraints imposed by surfactant molecules limit the access to this large asymptotic domain. However, by an exact analogy, it is shown how the giant bubble shapes can be realized by large inflatable structures.

Statistical Image Properties in Works from the Prinzhorn Collection of Artists with Schizophrenia

The Prinzhorn Collection preserves and exhibits thousands of visual artworks by patients who were diagnosed to suffer from mental disease. From this collection, we analyzed 1,256 images by 14 artists who were diagnosed with dementia praecox or schizophrenia. Six objective statistical properties that have been used previously to characterize visually aesthetic images were calculated. These properties reflect features of formal image composition, such as the complexity and distribution of oriented luminance gradients and edges, as well as Fourier spectral properties. Results for the artists with schizophrenia were compared to artworks from three public art collections of paintings and drawings that include highly acclaimed artworks as well as artworks of lesser artistic claim (control artworks). Many of the patients' works did not differ from these control images. However, the artworks of 6 of the 14 artists with schizophrenia possess image properties that deviate from the range of values obtained for the control artworks. For example, the artworks of four of the patients are characterized by a relative dominance of specific edge orientations in their images (low first-order entropy of edge orientations). Three patients created artworks with a relatively high ratio of fine detail to coarse structure (high slope of the Fourier spectrum). In conclusion, the present exploratory study opens novel perspectives for the objective scientific investigation of visual artworks that were created by persons who suffer from schizophrenia.

On the organization of human T-cell receptor loci: log-periodic distribution of T-cell receptor gene segments

The human T-cell repertoire is complex and is generated by the rearrangement of variable (V), diversity (D) and joining (J) segments on the T-cell receptor (TCR) loci. The T-cell repertoire demonstrates self-similarity in terms clonal frequencies when defined by V, D and J gene segment usage; therefore to determine whether the structural ordering of these gene segments on the TCR loci contributes to the observed clonal frequencies, the TCR loci were examined for self-similarity and periodicity in terms of gene segment organization. Logarithmic transformation of numeric sequence order demonstrated that the V and J gene segments for both T-cell receptor α (TRA) and β (TRB) loci are arranged in a self-similar manner when the spacing between the adjacent segments was considered as a function of the size of the neighbouring gene segment, with an average fractal dimension of approximately 1.5. Accounting for the gene segments occurring on helical DNA molecules with a logarithmic distribution, sine and cosine functions of the log-transformed angular coordinates of the start and stop nucleotides of successive TCR gene segments showed an ordered progression from the 5' to the 3' end of the locus, supporting a log-periodic organization. T-cell clonal frequency estimates, based on V and J segment usage, from normal stem cell donors were plotted against the V and J segment on TRB locus and demonstrated a periodic distribution. We hypothesize that this quasi-periodic variation in gene-segment representation in the T-cell clonal repertoire may be influenced by the location of the gene segments on the periodic-logarithmically scaled TCR loci. Interactions between the two strands of DNA in the double helix may influence the probability of gene segment usage by means of either constructive or destructive interference resulting from the superposition of the two helices.

The self-similarity theory of high pressure torsion

By analyzing the problem of high pressure torsion (HPT) in the rigid plastic formulation, we show that the power hardening law of plastically deformed materials leads to self-similarity of HPT, admitting a simple mathematical description of the process. The analysis shows that the main parameters of HPT are proportional to β ^q , with β being the angle of the anvil rotation. The meaning of the parameter q is: q = 0 for velocity and strain rate, q = 1 for shear strain and von Mises strain, q = n for stress, pressure and torque (n is the exponent of a power hardening law). We conclude that if the hardening law is a power law in a rotation interval β, self-similar regimes can emerge in HPT if the friction with the lateral wall of the die is not too high. In these intervals a simple mathematical description can be applied based on self-similarity. Outside these ranges, the plasticity problem still has to be solved for each value of β. The results obtained have important practical implications for the proper design and analysis of HPT experiments.

The dynamic generalization of the Eshelby inclusion problem and its static limit

The dynamic generalization of the celebrated Eshelby inclusion with transformation strain is the (subsonically) self-similarly expanding ellipsoidal inclusion starting from the zero dimension. The solution of the governing system of partial differential equations was obtained recently by Ni & Markenscoff (In press. J. Mech. Phys. Solids (doi:10.1016/j.jmps.2016.02.025)) on the basis of the Radon transformation, while here an alternative method is presented. In the self-similarly expanding motion, the Eshelby property of constant constrained strain is valid in the interior domain of the expanding ellipsoid where the particle velocity vanishes (lacuna). The dynamic Eshelby tensor is obtained in integral form. From it, the static Eshelby tensor is obtained by a limiting procedure, as the axes' expansion velocities tend to zero and time to infinity, while their product is equal to the length of the static axis. This makes the Eshelby problem the limit of its dynamic generalization.

Box-Counting Dimension Revisited: Presenting an Efficient Method of Minimizing Quantization Error and an Assessment of the Self-Similarity of Structural Root Systems

Fractal dimension (FD), estimated by box-counting, is a metric used to characterize plant anatomical complexity or space-filling characteristic for a variety of purposes. The vast majority of published studies fail to evaluate the assumption of statistical self-similarity, which underpins the validity of the procedure. The box-counting procedure is also subject to error arising from arbitrary grid placement, known as quantization error (QE), which is strictly positive and varies as a function of scale, making it problematic for the procedure's slope estimation step. Previous studies either ignore QE or employ inefficient brute-force grid translations to reduce it. The goals of this study were to characterize the effect of QE due to translation and rotation on FD estimates, to provide an efficient method of reducing QE, and to evaluate the assumption of statistical self-similarity of coarse root datasets typical of those used in recent trait studies. Coarse root systems of 36 shrubs were digitized in 3D and subjected to box-counts. A pattern search algorithm was used to minimize QE by optimizing grid placement and its efficiency was compared to the brute force method. The degree of statistical self-similarity was evaluated using linear regression residuals and local slope estimates. QE, due to both grid position and orientation, was a significant source of error in FD estimates, but pattern search provided an efficient means of minimizing it. Pattern search had higher initial computational cost but converged on lower error values more efficiently than the commonly employed brute force method. Our representations of coarse root system digitizations did not exhibit details over a sufficient range of scales to be considered statistically self-similar and informatively approximated as fractals, suggesting a lack of sufficient ramification of the coarse root systems for reiteration to be thought of as a dominant force in their development. FD estimates did not characterize the scaling of our digitizations well: the scaling exponent was a function of scale. Our findings serve as a caution against applying FD under the assumption of statistical self-similarity without rigorously evaluating it first.

Ants determine their next move at rest: motor planning and causality in complex systems

To find useful work to do for their colony, individual eusocial animals have to move, somehow staying attentive to relevant social information. Recent research on individual Temnothorax albipennis ants moving inside their colony's nest found a power-law relationship between a movement's duration and its average speed; and a universal speed profile for movements showing that they mostly fluctuate around a constant average speed. From this predictability it was inferred that movement durations are somehow determined before the movement itself. Here, we find similar results in lone T. albipennis ants exploring a large arena outside the nest, both when the arena is clean and when it contains chemical information left by previous nest-mates. This implies that these movement characteristics originate from the same individual neural and/or physiological mechanism(s), operating without immediate regard to social influences. However, the presence of pheromones and/or other cues was found to affect the inter-event speed correlations. Hence we suggest that ants' motor planning results in intermittent response to the social environment: movement duration is adjusted in response to social information only between movements, not during them. This environmentally flexible, intermittently responsive movement behaviour points towards a spatially allocated division of labour in this species. It also prompts more general questions on collective animal movement and the role of intermittent causation from higher to lower organizational levels in the stability of complex systems.