-
Critical behavior of the stochastic SIR model on random bond-diluted lattices
Authors:
Carlos Handrey A. Ferraz,
José Luiz S. Lima
Abstract:
In this paper, we investigate the impact of bond-dilution disorder on the critical behavior of the stochastic SIR model. Monte Carlo simulations were conducted using square lattices with first- and second-nearest neighbor interactions. Quenched bond-diluted lattice disorder was introduced into the systems, allowing them to evolve over time. By employing percolation theory and finite-size scaling a…
▽ More
In this paper, we investigate the impact of bond-dilution disorder on the critical behavior of the stochastic SIR model. Monte Carlo simulations were conducted using square lattices with first- and second-nearest neighbor interactions. Quenched bond-diluted lattice disorder was introduced into the systems, allowing them to evolve over time. By employing percolation theory and finite-size scaling analysis, we estimate both the critical threshold and leading critical exponent ratios of the model for different bond-dilution rates ($p$). An examination of the average size of the percolating cluster and the size distribution of non-percolating clusters of recovered individuals was performed to ascertain the universality class of the model. The simulation results strongly indicate that the present model belongs to a new universality class distinct from that of 2D dynamical percolation, depending on the specific $p$ value under consideration.
△ Less
Submitted 20 March, 2024;
originally announced March 2024.
-
Random packing dynamics of $Σ_{2v}(2π/3)$-triplets
Authors:
Carlos Handrey Araujo Ferraz
Abstract:
In this letter, we used a combination of DEM and the multi-sphere method to investigate the random packing dynamics of $Σ_{2v}(2π/3)$-triplets. These triplets consist of three overlapping primary spheres, forming a bent structure with a bond angle of $2π/3$ and belonging to the $C_{2v}$ symmetry group. The selection of this specific structure was motivated by its similarity to molecules such as wa…
▽ More
In this letter, we used a combination of DEM and the multi-sphere method to investigate the random packing dynamics of $Σ_{2v}(2π/3)$-triplets. These triplets consist of three overlapping primary spheres, forming a bent structure with a bond angle of $2π/3$ and belonging to the $C_{2v}$ symmetry group. The selection of this specific structure was motivated by its similarity to molecules such as water, which displays crucial physicochemical properties and finds extensive application in various fields. To ensure non-overlapping particles at the beginning of the simulations, the rectangular confinement box was divided into basic cells. Each triplet was then inserted into a basic cell with a random orientation. After that, the system is allowed to settle under gravity towards the bottom of the box. An implicit leapfrog algorithm with quaternion acceleration was used to numerically integrate the rotational motion equations. By assuming a molecular approach, we account for the long-range cohesive forces using a Lennard-Jones (LJ)-like potential. The packing processes are studied assuming different long-range interaction strengths. We performed statistical calculations of the different quantities studied including packing density, radial distribution function, and orientation pair correlation function. In addition, the force probability distributions in the random packing structures have been analyzed.
△ Less
Submitted 8 August, 2023;
originally announced August 2023.
-
A non-absorbing SIR stochastic lattice gas model on hybrid lattices
Authors:
Carlos Handrey Araujo Ferraz
Abstract:
In this paper, we perform Monte Carlo calculations to study the critical behavior of the spread of infectious diseases through a novel approach to the SIR epidemiological model. A stochastic lattice gas version of the model was applied on hybrid lattices which, in turn, are generated from typical square lattices when inserting a connection probability $p$ that a given lattice site has both first-…
▽ More
In this paper, we perform Monte Carlo calculations to study the critical behavior of the spread of infectious diseases through a novel approach to the SIR epidemiological model. A stochastic lattice gas version of the model was applied on hybrid lattices which, in turn, are generated from typical square lattices when inserting a connection probability $p$ that a given lattice site has both first- and second-nearest neighbor interactions. By combining percolation theory and finite-size scaling analysis, we estimate both the critical threshold and leading critical exponent ratios of the non-absorbing SIR model in different cases of hybrid lattices. An analysis of the average size of the percolating cluster and the size distribution of non-percolating clusters of recovered individuals was carried out to determine the universality class of the model.
△ Less
Submitted 30 November, 2021; v1 submitted 3 October, 2021;
originally announced October 2021.
-
Effects of the long-range cohesive forces in binary particle packing dynamics
Authors:
Carlos Handrey Araujo Ferraz
Abstract:
Studies on random packing of bidispersive particles have shown that such systems can capture the underlying behavior of more complex phenomena found in physics and materials engineering. In industry, bidispersive particles are used to allow the increase of density and fluidity of the formed compounds. The understanding of the dynamics of these processes is therefore of great theoretical and practi…
▽ More
Studies on random packing of bidispersive particles have shown that such systems can capture the underlying behavior of more complex phenomena found in physics and materials engineering. In industry, bidispersive particles are used to allow the increase of density and fluidity of the formed compounds. The understanding of the dynamics of these processes is therefore of great theoretical and practical interest. In this paper, we perform molecular dynamics (MD) simulations to study the packing process of particles with binary size distribution. Samples with different particle population densities ($p$) as well as particle size ratios ($λ$) have been assessed. The initial positions of five thousand non-overlapping particles are assigned inside a confining rectangular box. After that, the system is allowed to settle under gravity towards the bottom of the box. Both the translational and rotational movements of each particle are considered in the simulations. In order to deal with interacting particles, we take into account both the contact and long-range cohesive forces. The normal viscoelastic force is calculated according to the nonlinear Hertz model, whereas the tangential force is calculated through an accurate nonlinear-spring model. Assuming a molecular approach, we account for the long-range cohesive forces using a Lennard-Jones(LJ)-like potential. The packing processes are studied assuming different long-range interaction strengths. We carry out statistical calculations of the different quantities studied including packing density, radial distribution function and orientation pair correlation function.
△ Less
Submitted 19 October, 2020; v1 submitted 28 April, 2020;
originally announced April 2020.
-
Influence of the Long-Range Forces in Non-Gaussian Random-Packing Dynamics
Authors:
Carlos Handrey Araujo Ferraz
Abstract:
In this paper, we perform molecular dynamics (MD) simulations to study the random packing of spheres with different particle size distributions. In particular, we deal with non-Gaussian distributions by means of the Lévy distributions. The initial positions, as well as the radii of five thousand non-overlapping particles, are assigned inside a confining rectangular box. After that, the system is a…
▽ More
In this paper, we perform molecular dynamics (MD) simulations to study the random packing of spheres with different particle size distributions. In particular, we deal with non-Gaussian distributions by means of the Lévy distributions. The initial positions, as well as the radii of five thousand non-overlapping particles, are assigned inside a confining rectangular box. After that, the system is allowed to settle under gravity towards the bottom of the box. Both the translational and rotational movements of each particle are considered in the simulations. In order to deal with interacting particles, we take into account both the contact and long-range cohesive forces. The normal viscoelastic force is calculated according to the nonlinear Hertz model, whereas the tangential force is calculated through an accurate nonlinear-spring model. Assuming a molecular approach, we account for the long-range cohesive forces using a Lennard-Jones(LJ)-like potential. The packing processes are studied assuming different long-range interaction strengths.
△ Less
Submitted 14 August, 2019; v1 submitted 26 November, 2018;
originally announced November 2018.
-
Kauffman cellular automata on quasicrystal topology
Authors:
Carlos Handrey A. Ferraz,
José Luiz S. Lima
Abstract:
In this paper we perform numerical simulations to study Kauffman cellular automata (KCA) on quasiperiod lattices. In particular, we investigate phase transition, magnetic entropy and propagation speed of the damage on these lattices. Both the critical threshold parameter $p_{c}$ and the critical exponents are estimated with good precision. In order to investigate the increase of statistical fluctu…
▽ More
In this paper we perform numerical simulations to study Kauffman cellular automata (KCA) on quasiperiod lattices. In particular, we investigate phase transition, magnetic entropy and propagation speed of the damage on these lattices. Both the critical threshold parameter $p_{c}$ and the critical exponents are estimated with good precision. In order to investigate the increase of statistical fluctuations and the onset of chaos in the critical region of the model, we have also defined a magnetic entropy to these systems. It is seen that the magnetic entropy behaves in a different way when one passes from the frozen regime ($p<p_{c}$) to the chaotic regime ($p>p_{c}$). For a further analysis, the robustness of the propagation of failures is checked by introducing a quenched site dilution probability $q$ on the lattices. It is seen that the damage spreading is quite sensitive when a small fraction of the lattice sites are disconnected. A finite-size scaling analysis is employed to estimate the critical exponents. From these numerical estimates, we claim that on both pure ($q=0$) and diluted ($q=0.05$) quasiperiodic lattices, the KCA model belongs to the same universality class than on square lattices. Furthermore, with the aim of comparing the dynamical behavior between periodic and quasiperiodic systems, the propagation speed of the damage is also calculated for the square lattice assuming the same conditions. It is found that on square lattices the propagation speed of the damage obeys a power law as $v\sim (p-p_{c})^α$, whereas on quasiperiod lattices it follows a logarithmic law as $v \sim \ln(p-p_{c})^α$.
△ Less
Submitted 27 June, 2018; v1 submitted 26 March, 2018;
originally announced March 2018.
-
Dynamic Simulation of Random Packing of Polydispersive Fine Particles
Authors:
Carlos Handrey Araujo Ferraz,
Samuel Apolinário Marques
Abstract:
In this paper, we perform molecular dynamics (MD) simulations to study the two-dimensional packing process of both monosized and random size particles with radii ranging from $1.0 \, μm$ to $7.0 \, μm$. The system was allowed to settle under gravity towards the bottom of a $300 \, μm \times 500 \, μm$ rectangular box. The initial positions as well as the radii of five thousand fine particles were…
▽ More
In this paper, we perform molecular dynamics (MD) simulations to study the two-dimensional packing process of both monosized and random size particles with radii ranging from $1.0 \, μm$ to $7.0 \, μm$. The system was allowed to settle under gravity towards the bottom of a $300 \, μm \times 500 \, μm$ rectangular box. The initial positions as well as the radii of five thousand fine particles were defined along the box by using a random number generator. Both the translational and the rotational movement of each particle were considered in the simulations. In order to deal with interacting fine particles, we take into account both the contact forces and the long-range dispersive forces. We account for normal and static/sliding tangential friction forces between particles and between particle and wall by means of a linear model approach, while the long-range dispersive forces are computed by using a Lennard-Jones like potential. The packing processes were studied assuming different long- range interaction strengths. We carry out statistical calculations of the different quantities studied such as packing density, mean coordination number and time derivative of the kinetic energy as the system evolves over time. A size spectral analysis was employed to obtain the radial distribution function (RDF) of the random close-packed structures (RCPS) for the case of random size particles. We find that the long-range dispersive forces can strongly influence the packing process dynamics as they might form large particle clusters, depending on the intensity of the long-range interaction strength. However, the general shape of the RDFs for the RCPS is seen to be more influenced by the hardness of the particles than by the long-range dispersive forces.
△ Less
Submitted 30 November, 2017; v1 submitted 3 December, 2016;
originally announced December 2016.
-
Three-state Potts model on Non-local Directed Small-World Lattices
Authors:
Carlos Handrey Araujo Ferraz,
José Luiz Sousa Lima
Abstract:
In this paper, we study the effects of non-local directed Small-World (NLDSW) disorder in the three-state Potts model as a form to capture the essential features shared by real complex systems where non-locality effects play a important role in the behavior of these systems. Using Monte Carlo techniques and finite-size scaling analysis, we estimate the infinite lattice critical temperatures and th…
▽ More
In this paper, we study the effects of non-local directed Small-World (NLDSW) disorder in the three-state Potts model as a form to capture the essential features shared by real complex systems where non-locality effects play a important role in the behavior of these systems. Using Monte Carlo techniques and finite-size scaling analysis, we estimate the infinite lattice critical temperatures and the leading critical exponents in this model. In particular, we investigate the first- to second-order phase transition crossover when NLDSW links are inserted. A cluster-flip algorithm was used to reduce the critical slowing down effect in our simulations. We find that for a NLDSW disorder densities $p<p^{*}=0.05(4)$, the model exhibits a continuous phase transition falling into a new universality class, which continuously depends on $p$, while for $p^{*}\leqslant p \leqslant 1.0$, the model presents a weak first-order phase transition.
△ Less
Submitted 27 May, 2017; v1 submitted 5 September, 2016;
originally announced September 2016.
-
Critical Behavior of the q = 3, 4-Potts model on Quasiperiodic Decagonal Lattices
Authors:
Carlos Handrey Araujo Ferraz
Abstract:
In this study, we performed Monte Carlo simulations of the $q=3,4$-Potts model on quasiperiodic decagonal lattices (QDL) to assess the critical behavior of these systems. Using the single histogram technique in conjunction with the finite-size scaling analysis, we estimate the infinite lattice critical temperatures and the leading critical exponents for $q=3$ and $q=4$ states. Our estimates for th…
▽ More
In this study, we performed Monte Carlo simulations of the $q=3,4$-Potts model on quasiperiodic decagonal lattices (QDL) to assess the critical behavior of these systems. Using the single histogram technique in conjunction with the finite-size scaling analysis, we estimate the infinite lattice critical temperatures and the leading critical exponents for $q=3$ and $q=4$ states. Our estimates for the critical exponents on QDL are in good agreement with the exact values on 2D periodic lattices, supporting the claim that both the $q=3$ and $q=4$ Potts model on quasiperiodic lattices belong to the same universality class as those on 2D periodic lattices.
△ Less
Submitted 1 October, 2015; v1 submitted 28 May, 2014;
originally announced May 2014.
-
The Strange Man in Random Networks of Automata
Authors:
Carlos Handrey A. Ferraz,
Hans J. Herrmann
Abstract:
We have performed computer simulations of Kauffman's automata on several graphs such as the regular square lattice and invasion percolation clusters in order to investigate phase transitions, radial distributions of the mean total damage (dynamical exponent $z$) and propagation speeds of the damage when one adds a damaging agent, nicknamed "strange man". Despite the increase in the damaging effi…
▽ More
We have performed computer simulations of Kauffman's automata on several graphs such as the regular square lattice and invasion percolation clusters in order to investigate phase transitions, radial distributions of the mean total damage (dynamical exponent $z$) and propagation speeds of the damage when one adds a damaging agent, nicknamed "strange man". Despite the increase in the damaging efficiency, we have not observed any appreciable change at the transition threshold to chaos neither for the short-range nor for the small-world case on the square lattices when the strange man is added in comparison to when small initial damages are inserted in the system. The propagation speed of the damage cloud until touching the border of the system in both cases obeys a power law with a critical exponent $α$ that strongly depends on the lattice. Particularly, we have ckecked the damage spreading when some connections are removed on the square lattice and when one considers special invasion percolation clusters (high boundary-saturation clusters). It is seen that the propagation speed in these systems is quite sensible to the degree of dilution.
△ Less
Submitted 7 February, 2007;
originally announced February 2007.
-
The Kauffman model on Small-World Topology
Authors:
Carlos Handrey A. Ferraz,
Hans J. Herrmann
Abstract:
We apply Kauffman's automata on small-world networks to study the crossover between the short-range and the infinite-range case. We perform accurate calculations on square lattices to obtain both critical exponents and fractal dimensions. Particularly, we find an increase of the damage propagation and a decrease in the fractal dimensions when adding long-range connections.
We apply Kauffman's automata on small-world networks to study the crossover between the short-range and the infinite-range case. We perform accurate calculations on square lattices to obtain both critical exponents and fractal dimensions. Particularly, we find an increase of the damage propagation and a decrease in the fractal dimensions when adding long-range connections.
△ Less
Submitted 18 April, 2006; v1 submitted 30 March, 2006;
originally announced March 2006.