Search

One way to study the spread of disease is through mathematical models. The most successful models compartmentalize the host population according to their infectious stage, e.g., susceptible (S), infected (I), exposed (E), and recovered (R). The composition of these compartments leads to the SI, SIS, SIR, and SEIR models. In this Chapter, we present and compare three formulations of SI, SIS, SIR, and SEIR models in the framework of standard (integer operators), fractional (Caputo sense), and fractal derivatives (Hausdorff sense). As an application of the SI model, we study the evolution of AIDS cases in Bangladesh from 2001 to 2021. For this case, our simulations suggest that fractal formulation describes the data well. For the SIS model, we consider syphilis data from Brazil from 2006 to 2017. In this case, the three frameworks describe the data with good accuracy. We used data from Influenza A to adjust the SIR model in previous approaches and observed that the fractional formulation was better. The last application considers the COVID-19 data from India in the range 2020-04-10 to 2020-12-31 to adjust the parameters of the SEIR model. The standard formulation fits the data better than the other approaches. As a common result, all models exhibit steady solutions in the different formulations. The time to reach a steady solution is correlated to the considered approach. The standard and fractal formulations reach the steady state earlier when compared with the fractional formulation.

Permutation entropy and its associated frameworks are remarkable examples of physics-inspired techniques adept at processing complex and extensive datasets. Despite substantial progress in developing and applying these tools, their use has been predominantly limited to structured datasets such as time series or images. Here, we introduce the k-nearest neighbor permutation entropy, an innovative extension of the permutation entropy tailored for unstructured data, irrespective of their spatial or temporal configuration and dimensionality. Our approach builds upon nearest neighbor graphs to establish neighborhood relations and uses random walks to extract ordinal patterns and their distribution, thereby defining the k-nearest neighbor permutation entropy. This tool not only adeptly identifies variations in patterns of unstructured data, but also does so with a precision that significantly surpasses conventional measures such as spatial autocorrelation. Additionally, it provides a natural approach for incorporating amplitude information and time gaps when analyzing time series or images, thus significantly enhancing its noise resilience and predictive capabilities compared to the usual permutation entropy. Our research substantially expands the applicability of ordinal methods to more general data types, opening promising research avenues for extending the permutation entropy toolkit for unstructured data., Comment: 14 two-column pages, 6 figures

The description of neuronal activity has been of great importance in neuroscience. In this field, mathematical models are useful to describe the electrophysical behaviour of neurons. One successful model used for this purpose is the Adaptive Exponential Integrate-and-Fire (Adex), which is composed of two ordinary differential equations. Usually, this model is considered in his standard formulation, i.e., with integer order derivatives. In this work, we propose and study the fractal extension of Adex model, which in simple terms corresponds to replacing the integer derivative by non-integer. As non-integer operators, we choose the fractal derivatives. We explore the effects of equal and different orders of fractal derivatives in the firing patterns and mean frequency of the neuron described by the Adex model. Previous results suggest that fractal derivatives can provide a more realistic representation due to the fact that the standard operators are generalized. Our findings show that the fractal order influences the inter-spike intervals and changes the mean firing frequency. In addition, the firing patterns depend not only on the neuronal parameters but also on the order of respective fractal operators. As our main conclusion, the fractal order below the unit value increases the influence of the adaptation mechanism in the spike firing patterns.

This work studies the SIS model extended by fractional and fractal derivatives. We obtain explicit solutions for the standard and fractal formulations; for the fractional case, we study numerical solutions. As a real data example, we consider the Brazilian syphilis data from 2011 to 2021. We fit the data by considering the three variations of the model. Our fit suggests a recovery period of 11.6 days and a reproduction ratio ($R_0$) equal to 6.5. By calculating the correlation coefficient ($r$) between the real data and the theoretical points, our results suggest that the fractal model presents a higher $r$ compared to the standard or fractional case. The fractal formulation is improved when two different fractal orders with distinguishing weights are considered. This modification in the model provides a better description of the data and improves the correlation coefficient.

In this work, we analyze the effects of fractional derivatives in the chaotic dynamics of a cancer model. We begin by studying the dynamics of a standard model, {\it i.e.}, with integer derivatives. We study the dynamical behavior by means of the bifurcation diagram, Lyapunov exponents, and recurrence quantification analysis (RQA), such as the recurrence rate (RR), the determinism (DET), and the recurrence time entropy (RTE). We find a high correlation coefficient between the Lyapunov exponents and RTE. Our simulations suggest that the tumor growth parameter ($\rho_1$) is associated with a chaotic regime. Our results suggest a high correlation between the largest Lyapunov exponents and RTE. After understanding the dynamics of the model in the standard formulation, we extend our results by considering fractional operators. We fix the parameters in the chaotic regime and investigate the effects of the fractional order. We demonstrate how fractional dynamics can be properly characterized using RQA measures, which offer the advantage of not requiring knowledge of the fractional Jacobian matrix. We find that the chaotic motion is suppressed as $\alpha$ decreases, and the system becomes periodic for $\alpha \lessapprox 0.9966$. We observe limit cycles for $\alpha \in (0.9966,0.899)$ and fixed points for $\alpha<0.899$. The fixed point is determined analytically for the considered parameters. Finally, we discover that these dynamics are separated by an exponential relationship between $\alpha$ and $\rho_1$. Also, the transition depends on a supper transient which obeys the same relationship.

Heterogeneous media diffusion is often described using position-dependent diffusion coefficients and estimated indirectly through mean squared displacement in experiments. This approach may overlook other mechanisms and their interaction with position-dependent diffusion, potentially leading to erroneous conclusions. Here, we introduce a hybrid diffusion model that merges a position-dependent diffusion coefficient with the trapping mechanism of the comb model. We derive exact solutions for position distributions and mean squared displacements, validated through simulations of Langevin equations. Our model shows that the trapping mechanism attenuates the impact of media heterogeneity. Superdiffusion occurs when the position-dependent coefficient increases superlinearly, while subdiffusion occurs for sublinear and inverse power-law relations. This nontrivial interplay between heterogeneity and state-independent mechanisms also leads to anomalous yet Brownian and non-Brownian yet Gaussian regimes. These findings emphasize the need for cautious interpretations of experiments and highlight the limitations of relying solely on mean squared displacements or position distributions for diffusion characterization., Comment: 13 two-column pages, 6 figures; accepted for publication in Communications Physics

Recent advances in deep learning methods have enabled researchers to develop and apply algorithms for the analysis and modeling of complex networks. These advances have sparked a surge of interest at the interface between network science and machine learning. Despite this, the use of machine learning methods to investigate criminal networks remains surprisingly scarce. Here, we explore the potential of graph convolutional networks to learn patterns among networked criminals and to predict various properties of criminal networks. Using empirical data from political corruption, criminal police intelligence, and criminal financial networks, we develop a series of deep learning models based on the GraphSAGE framework that are able to recover missing criminal partnerships, distinguish among types of associations, predict the amount of money exchanged among criminal agents, and even anticipate partnerships and recidivism of criminals during the growth dynamics of corruption networks, all with impressive accuracy. Our deep learning models significantly outperform previous shallow learning approaches and produce high-quality embeddings for node and edge properties. Moreover, these models inherit all the advantages of the GraphSAGE framework, including the generalization to unseen nodes and scaling up to large graph structures., Comment: 14 two-column pages, 5 figures

Recent research has shown that criminal networks have complex organizational structures, but whether this can be used to predict static and dynamic properties of criminal networks remains little explored. Here, by combining graph representation learning and machine learning methods, we show that structural properties of political corruption, police intelligence, and money laundering networks can be used to recover missing criminal partnerships, distinguish among different types of criminal and legal associations, as well as predict the total amount of money exchanged among criminal agents, all with outstanding accuracy. We also show that our approach can anticipate future criminal associations during the dynamic growth of corruption networks with significant accuracy. Thus, similar to evidence found at crime scenes, we conclude that structural patterns of criminal networks carry crucial information about illegal activities, which allows machine learning methods to predict missing information and even anticipate future criminal behavior., Comment: 10 pages, 4 figures, supplementary information; accepted for publication in Scientific Reports

It is well-known that the theories of semi-vector spaces and semi-algebras -- which were not much studied over time -- are utilized/applied in Fuzzy Set Theory in order to obtain extensions of the concept of fuzzy numbers as well as to provide new mathematical tools to investigate properties and new results on fuzzy systems. In this paper we investigate the theory of semi-vector spaces over the semi-field of nonnegative real numbers ${\mathbb R}_{0}^{+}$. We prove several results concerning semi-vector spaces and semi-linear transformations. Moreover, we introduce in the literature the concept of eigenvalues and eigenvectors of a semi-linear operator, describing in some cases how to compute them. Topological properties of semi-vector spaces such as completeness and separability are also investigated. New families of semi-vector spaces derived from semi-metric, semi-norm, semi-inner product, among others are exhibited. Additionally, some results on semi-algebras are presented.

Machine learning algorithms have been available since the 1990s, but it is much more recently that they have come into use also in the physical sciences. While these algorithms have already proven to be useful in uncovering new properties of materials and in simplifying experimental protocols, their usage in liquid crystals research is still limited. This is surprising because optical imaging techniques are often applied in this line of research, and it is precisely with images that machine learning algorithms have achieved major breakthroughs in recent years. Here we use convolutional neural networks to probe several properties of liquid crystals directly from their optical images and without using manual feature engineering. By optimizing simple architectures, we find that convolutional neural networks can predict physical properties of liquid crystals with exceptional accuracy. We show that these deep neural networks identify liquid crystal phases and predict the order parameter of simulated nematic liquid crystals almost perfectly. We also show that convolutional neural networks identify the pitch length of simulated samples of cholesteric liquid crystals and the sample temperature of an experimental liquid crystal with very high precision., Comment: 11 pages, 4 figures, supplementary information; accepted for publication in Scientific Reports

In this note, a numerical method based on finite differences to solve a class of nonlinear advection-diffusion fractional differential equation is proposed. The fractional operator considered here is the fractional Riemann-Liouville derivative or the fractional Riesz derivative of order $\alpha$. The consistency and unconditionally stability of the method are shown. Finally, an example of application of this method is presented., Comment: 9 pages, 2 figures

One of the most useful tools for distinguishing between chaotic and stochastic time series is the so-called complexity-entropy causality plane. This diagram involves two complexity measures: the Shannon entropy and the statistical complexity. Recently, this idea has been generalized by considering the Tsallis monoparametric generalization of the Shannon entropy, yielding complexity-entropy curves. These curves have proven to enhance the discrimination among different time series related to stochastic and chaotic processes of numerical and experimental nature. Here we further explore these complexity-entropy curves in the context of the R\'enyi entropy, which is another monoparametric generalization of the Shannon entropy. By combining the R\'enyi entropy with the proper generalization of the statistical complexity, we associate a parametric curve (the R\'enyi complexity-entropy curve) with a given time series. We explore this approach in a series of numerical and experimental applications, demonstrating the usefulness of this new technique for time series analysis. We show that the R\'enyi complexity-entropy curves enable the differentiation among time series of chaotic, stochastic, and periodic nature. In particular, time series of stochastic nature are associated with curves displaying positive curvature in a neighborhood of their initial points, whereas curves related to chaotic phenomena have a negative curvature; finally, periodic time series are represented by vertical straight lines., Comment: 19 pages, 6 figures; accepted for publication in Physica A

Corruptive behaviour in politics limits economic growth, embezzles public funds, and promotes socio-economic inequality in modern democracies. We analyse well-documented political corruption scandals in Brazil over the past 27 years, focusing on the dynamical structure of networks where two individuals are connected if they were involved in the same scandal. Our research reveals that corruption runs in small groups that rarely comprise more than eight people, in networks that have hubs and a modular structure that encompasses more than one corruption scandal. We observe abrupt changes in the size of the largest connected component and in the degree distribution, which are due to the coalescence of different modules when new scandals come to light or when governments change. We show further that the dynamical structure of political corruption networks can be used for successfully predicting partners in future scandals. We discuss the important role of network science in detecting and mitigating political corruption., Comment: 20 pages, 5 figures; accepted for publication in Journal of Complex Networks

The generalized diffusion equations with fractional order derivatives have shown be quite efficient to describe the diffusion in complex systems, with the advantage of producing exact expressions for the underlying diffusive properties. Recently, researchers have proposed different fractional-time operators (namely: the Caputo-Fabrizio and Atangana-Baleanu) which, differently from the well-known Riemann-Liouville operator, are defined by non-singular memory kernels. Here we proposed to use these new operators to generalize the usual diffusion equation. By analyzing the corresponding fractional diffusion equations within the continuous time random walk framework, we obtained waiting time distributions characterized by exponential, stretched exponential, and power-law functions, as well as a crossover between two behaviors. For the mean square displacement, we found crossovers between usual and confined diffusion, and between usual and sub-diffusion. We obtained the exact expressions for the probability distributions, where non-Gaussian and stationary distributions emerged. This former feature is remarkable because the fractional diffusion equation is solved without external forces and subjected to the free diffusion boundary conditions. We have further shown that these new fractional diffusion equations are related to diffusive processes with stochastic resetting, and to fractional diffusion equations with derivatives of distributed order. Thus, our results show that these new operators are a simple and efficient way for incorporating different structural aspects into the system, opening new possibilities for modeling and investigating anomalous diffusive processes.

The search for patterns in time series is a very common task when dealing with complex systems. This is usually accomplished by employing a complexity measure such as entropies and fractal dimensions. However, such measures usually only capture a single aspect of the system dynamics. Here, we propose a family of complexity measures for time series based on a generalization of the complexity-entropy causality plane. By replacing the Shannon entropy by a mono-parametric entropy (Tsallis $q$-entropy) and after considering the proper generalization of the statistical complexity ($q$-complexity), we build up a parametric curve (the $q$-complexity-entropy curve) that is used for characterizing/classifying time series. Based on simple exact results and numerical simulations of stochastic processes, we show that these curves can distinguish among different long-range, short-range and oscillating correlated behaviors. Also, we verify that simulated chaotic and stochastic time series can be distinguished based on whether these curves are open or closed. We further test this technique in experimental scenarios related to chaotic laser intensity, stock price, sunspot, and geomagnetic dynamics, confirming its usefulness. Finally, we prove that these curves enhance the automatic classification of time series with long-range correlations and interbeat intervals of healthy subjects and patients with heart disease., Comment: Accepted for publication in PRE

Statistical similarities between earthquakes and other systems that emit cracking noises have been explored in diverse contexts, ranging from materials science to financial and social systems. Such analogies give promise of a unified and universal theory for describing the complex responses of those systems. There are, however, very few attempts to simultaneously characterize the most fundamental seismic laws in such systems. Here we present a complete description of the Gutenberg-Richter law, the recurrence times, Omori's law, the productivity law, and Bath's law for the acoustic emissions that happen in the relaxation process of uncrumpling thin plastic sheets. Our results show that these laws also appear in this phenomenon, but (for most cases) with different parameters from those reported for earthquakes and fracture experiments. This study thus contributes to elucidate the parallel between seismic laws and cracking noises in uncrumpling processes, revealing striking qualitative similarities but also showing that these processes display unique features., Comment: Accepted for publication in EPL

More than a half of world population is now living in cities and this number is expected to be two-thirds by 2050. Fostered by the relevancy of a scientific characterization of cities and for the availability of an unprecedented amount of data, academics have recently immersed in this topic and one of the most striking and universal finding was the discovery of robust allometric scaling laws between several urban indicators and the population size. Despite that, most governmental reports and several academic works still ignore these nonlinearities by often analyzing the raw or the per capita value of urban indicators, a practice that actually makes the urban metrics biased towards small or large cities depending on whether we have super or sublinear allometries. By following the ideas of Bettencourt et al., we account for this bias by evaluating the difference between the actual value of an urban indicator and the value expected by the allometry with the population size. We show that this scale-adjusted metric provides a more appropriate/informative summary of the evolution of urban indicators and reveals patterns that do not appear in the evolution of per capita values of indicators obtained from Brazilian cities. We also show that these scale-adjusted metrics are strongly correlated with their past values by a linear correspondence and that they also display crosscorrelations among themselves. Simple linear models account for 31%-97% of the observed variance in data and correctly reproduce the average of the scale-adjusted metric when grouping the cities in above and below the allometric laws. We further employ these models to forecast future values of urban indicators and, by visualizing the predicted changes, we verify the emergence of spatial clusters characterized by regions of the Brazilian territory where we expect an increase or a decrease in the values of urban indicators., Comment: Accepted for publication in PLoS ONE

We report on the time dependent solutions of the $q-$generalized Schr\"odinger equation proposed by Nobre et al. [Phys. Rev. Lett. 106, 140601 (2011)]. Here we investigate the case of two free particles and also the case where two particles were subjected to a Moshinsky-like potential with time dependent coefficients. We work out analytical and numerical solutions for different values of the parameter $q$ and also show that the usual Schr\"odinger equation is recovered in the limit of $q\rightarrow 1$. An intriguing behavior was observed for $q=2$, where the wave function displays a ring-like shape, indicating a bind behavior of the particles. Differently from the results previously reported for the case of one particle, frozen states appear only for special combinations of the wave function parameters in case of $q=3$., Comment: Accepted for publication in Physica A

We report on the existing connection between power-law distributions and allometries. As it was first reported in [PLoS ONE 7, e40393 (2012)] for the relationship between homicides and population, when these urban indicators present asymptotic power-law distributions, they can also display specific allometries among themselves. Here, we present an extensive characterization of this connection when considering all possible pairs of relationships from twelve urban indicators of Brazilian cities (such as child labor, illiteracy, income, sanitation and unemployment). Our analysis reveals that all our urban indicators are asymptotically distributed as power laws and that the proposed connection also holds for our data when the allometric relationship displays enough correlations. We have also found that not all allometric relationships are independent and that they can be understood as a consequence of the allometric relationship between the urban indicator and the population size. We further show that the residuals fluctuations surrounding the allometries are characterized by an almost constant variance and log-normal distributions., Comment: Accepted for publication in Physica A

We report on a quantitative analysis of relationships between the number of homicides, population size and other ten urban metrics. By using data from Brazilian cities, we show that well defined average scaling laws with the population size emerge when investigating the relations between population and number of homicides as well as population and urban metrics. We also show that the fluctuations around the scaling laws are log-normally distributed, which enabled us to model these scaling laws by a stochastic-like equation driven by a multiplicative and log-normally distributed noise. Because of the scaling laws, we argue that it is better to employ logarithms in order to describe the number of homicides in function of the urban metrics via regression analysis. In addition to the regression analysis, we propose an approach to correlate crime and urban metrics via the evaluation of the distance between the actual value of the number of homicides (as well as the value of the urban metrics) and the value that is expected by the scaling law with the population size. This approach have proved to be robust and useful for unveiling relationships/behaviors that were not properly carried out by the regression analysis, such as i) the non-explanatory potential of the elderly population when the number of homicides is much above or much below the scaling law, ii) the fact that unemployment has explanatory potential only when the number of homicides is considerably larger than the expected by the power law, and iii) a gender difference in number of homicides, where cities with female population below the scaling law are characterized by a number of homicides above the power law., Comment: Accepted for publication in PLoS ONE

The Poisson-Nernst-Planck (PNP) diffusional model for the immittance or impedance spectroscopy response of an electrolytic cell in a finite-length situation is extended to a general framework. In this new formalism, the bulk behavior of the mobile charges is governed by a fractional diffusion equation in the presence of a reaction term. The solutions have to satisfy a general boundary condition embodying, in a single expression, most of the surface effects commonly encountered in experimental situations. Among these effects, we specifically consider the charge transfer process from an electrolytic cell to the external circuit and the adsorption-desorption phenomenon at the interfaces. The equations are exactly solved in the small AC signal approximation and are used to obtain an exact expression for the electrical impedance as a funcion of the frequency. The predictions of the model are compared to and found to be in good agreement with the experimental data obtained for an electrolytic solution of CdCl2H20., Comment: 5 figures