1
|
Rahimi-Majd M, Shirzad T, Najafi MN. A self-organized critical model and multifractal analysis for earthquakes in Central Alborz, Iran. Sci Rep 2022; 12:8364. [PMID: 35589782 PMCID: PMC9120491 DOI: 10.1038/s41598-022-12362-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2021] [Accepted: 04/28/2022] [Indexed: 11/25/2022] Open
Abstract
This paper is devoted to a phenomenological study of the earthquakes in central Alborz, Iran. Using three observational quantities, namely the weight function, the quality factor, and the velocity model in this region, we develop a modified dissipative sandpile model which captures the main features of the system, especially the average activity field over the region of study. The model is based on external stimuli, the location of which is chosen (I) randomly, (II) on the faults, (III) on the low active points, (IV) on the moderately active points, and (V) on the highly active points in the region. We uncover some universal behaviors depending slightly on the method of external stimuli. A multi-fractal detrended fluctuation analysis is exploited to extract the spectrum of the Hurst exponent of the time series obtained by each of these schemes. Although the average Hurst exponent depends slightly on the method of stimuli, we numerically show that in all cases it is lower than 0.5, reflecting the anti-correlated nature of the system. The lowest average Hurst exponent is found to be associated with the case (V), in such a way that the more active the stimulated sites are, the lower the average Hurst exponent is obtained, i.e. the large earthquakes are more anticorrelated. Moreover, we find that the activity field achieved in this study provide information about the depth and topography of the basement, and also the area that can potentially be the location of the future large events. We successfully determine a high activity zone on the Mosha Fault, where the mainshock occurred on May 7th, 2020 (M\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$_W$$\end{document}W 4.9).
Collapse
Affiliation(s)
- M Rahimi-Majd
- Department of Physics, Shahid Beheshti University, 1983969411, Tehran, Iran
| | - T Shirzad
- Institute of Geophysics, Polish Academy of Sciences - 01-452, Warsaw, Poland
| | - M N Najafi
- Department of Physics, University of Mohaghegh Ardabili, P.O. Box 179, Ardabil, Iran.
| |
Collapse
|
2
|
Valizadeh N, Samadpour M, Hamzehpour H, Najafi MN. Edwards-Wilkinson depinning transition in random Coulomb potential background. Phys Rev E 2021; 104:064140. [PMID: 35030907 DOI: 10.1103/physreve.104.064140] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2021] [Accepted: 12/01/2021] [Indexed: 06/14/2023]
Abstract
The quenched Edwards-Wilkinson growth of the 1+1 interface is considered in the background of the correlated random noise. We use random Coulomb potential as the background long-range correlated noise. A depinning transition is observed in a critical driving force F[over ̃]_{c}≈0.037 (in terms of disorder strength unit) in the vicinity of which the final velocity of the interface varies linearly with time. Our data collapse analysis for the velocity shows a crossover time t^{*} at which the velocity is size independent. Based on a two-variable scaling analysis, we extract the exponents, which are different from all universality classes we are aware of. Especially noting that the dynamic and roughness exponents are z_{w}=1.55±0.05, and α_{w}=1.05±0.05 at the criticality, we conclude that the system is different from both Edwards-Wilkinson (EW) and Kardar-Parisi-Zhang (KPZ) universality classes. Our analysis shows therefore that making the noise long-range correlated, drives the system out of the EW universality class. The simulations on the tilted lattice show that the nonlinearity term (λ term in the KPZ equations) goes to zero in the thermodynamic limit.
Collapse
Affiliation(s)
- N Valizadeh
- Department of Physics, K.N. Toosi University of Technology, Tehran 15875-4416, Iran
| | - M Samadpour
- Department of Physics, K.N. Toosi University of Technology, Tehran 15875-4416, Iran
| | - H Hamzehpour
- Department of Physics, K.N. Toosi University of Technology, Tehran 15875-4416, Iran
| | - M N Najafi
- Department of Physics, University of Mohaghegh Ardabili, P.O. Box 179, Ardabil, Iran
| |
Collapse
|
3
|
Najafi MN, Tizdast S, Cheraghalizadeh J, N HD. Invasion percolation in short-range and long-range disorder background. Phys Rev E 2021; 104:064119. [PMID: 35030889 DOI: 10.1103/physreve.104.064119] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/10/2021] [Accepted: 11/23/2021] [Indexed: 06/14/2023]
Abstract
In the original invasion percolation model, a random number quantifies the role of necks, or generally the quality of pores, ignoring the structure of pores and impermeable regions (to which the invader cannot enter). In this paper, we investigate invasion percolation (IP), taking into account the impermeable regions, the configuration of which is modeled by ordinary and Ising-correlated site percolation (with short-range interactions, SRI), on top of which the IP dynamics is defined. We model the long-ranged correlations of pores by a random Coulomb potential (RCP). By examining various dynamical observables, we suggest that the critical exponents of Ising-correlated cases change considerably only in the vicinity of the critical point (critical temperature), while for the ordinary percolation case the exponents are robust against the occupancy parameter p. The properties of the model for the long-range interactions [LRI (RCP)] are completely different from the normal IP. In particular, the fractal dimension of the external frontier of the largest hole is nearly 4/3 for SRI far from the critical points, which is compatible with normal IP, while it converges to 1.099±0.04 for RCP. For the latter case, the time dependence of our observables is divided into three parts: the power law (short time), the logarithmic (moderate time), and the linear (long time) regimes. The second crossover time is shown to go to infinity in the thermodynamic limit, whereas the first crossover time is nearly unchanged, signaling the dominance of the logarithmic regime. The average gyration radius of the growing clusters, the length of their external perimeter, and the corresponding roughness are shown to be nearly constant for the long-time regime in the thermodynamic limit.
Collapse
Affiliation(s)
- M N Najafi
- Department of Physics, University of Mohaghegh Ardabili, P.O. Box 179, Ardabil, Iran
| | - S Tizdast
- Department of Physics, University of Mohaghegh Ardabili, P.O. Box 179, Ardabil, Iran
| | - J Cheraghalizadeh
- Department of Physics, University of Mohaghegh Ardabili, P.O. Box 179, Ardabil, Iran
| | - H Dashti N
- School of Physics, Korea Institute for Advanced Study, Seoul 02455, Korea
| |
Collapse
|
4
|
Cheraghalizadeh J, Seifi M, Ebadi Z, Mohammadzadeh H, Najafi MN. Superstatistical two-temperature Ising model. Phys Rev E 2021; 103:032104. [PMID: 33862766 DOI: 10.1103/physreve.103.032104] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2020] [Accepted: 02/12/2021] [Indexed: 11/07/2022]
Abstract
The previous approach of the nonequilibrium Ising model was based on the local temperature in which each site or part of the system has its own specific temperature. We introduce an approach of the two-temperature Ising model as a prototype of the superstatistic critical phenomena. The model is described by two temperatures (T_{1},T_{2}) in a zero magnetic field. To predict the phase diagram and numerically estimate the exponents, we develop the Metropolis and Swendsen-Wang Monte Carlo method. We observe that there is a nontrivial critical line, separating ordered and disordered phases. We propose an analytic equation for the critical line in the phase diagram. Our numerical estimation of the critical exponents illustrates that all points on the critical line belong to the ordinary Ising universality class.
Collapse
Affiliation(s)
- J Cheraghalizadeh
- Department of Physics, University of Mohaghegh Ardabili, P.O. Box 179, Ardabil, Iran
| | - M Seifi
- Department of Physics, University of Mohaghegh Ardabili, P.O. Box 179, Ardabil, Iran
| | - Z Ebadi
- Department of Physics, University of Mohaghegh Ardabili, P.O. Box 179, Ardabil, Iran
| | - H Mohammadzadeh
- Department of Physics, University of Mohaghegh Ardabili, P.O. Box 179, Ardabil, Iran
| | - M N Najafi
- Department of Physics, University of Mohaghegh Ardabili, P.O. Box 179, Ardabil, Iran
| |
Collapse
|
5
|
Najafi MN, Cheraghalizadeh J, Herrmann HJ. Elastic backbone phase transition in the Ising model. Phys Rev E 2019; 100:042132. [PMID: 31770915 DOI: 10.1103/physreve.100.042132] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2019] [Indexed: 11/07/2022]
Abstract
The two-dimensional (zero magnetic field) Ising model is known to undergo a second-order paraferromagnetic phase transition, which is accompanied by a correlated percolation transition for the Fortuin-Kasteleyn (FK) clusters. In this paper we uncover that there exists also a second temperature T_{eb}<T_{c} at which the elastic backbone of FK clusters undergoes a second-order phase transition to a dense phase. The corresponding universality class, which is characterized by determining various percolation exponents, is shown to be completely different from directed percolation, which leads us to propose a new anisotropic universality class with β=0.54±0.02, ν_{||}=1.86±0.01, ν_{⊥}=1.21±0.04, and d_{f}=1.53±0.03. All tested hyperscaling relations are shown to be valid.
Collapse
Affiliation(s)
- M N Najafi
- Department of Physics, University of Mohaghegh Ardabili, P.O. Box 179, Ardabil, Iran.,Computational Physics, IfB, ETH Zurich, Stefano-Franscini-Platz 3, CH-8093 Zurich, Switzerland
| | - J Cheraghalizadeh
- Department of Physics, University of Mohaghegh Ardabili, P.O. Box 179, Ardabil, Iran
| | - H J Herrmann
- Computational Physics, IfB, ETH Zurich, Stefano-Franscini-Platz 3, CH-8093 Zurich, Switzerland.,Departamento de Física, Universidade Federal do Ceara, 60451-970 Fortaleza, Brazil.,ESPCI, CNRS UMR 7636, Laboratoire PMMH, 75005 Paris, France
| |
Collapse
|
6
|
Cheraghalizadeh J, Najafi MN, Mohammadzadeh H, Saber A. Self-avoiding walk on a square lattice with correlated vacancies. Phys Rev E 2018; 97:042128. [PMID: 29758691 DOI: 10.1103/physreve.97.042128] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/30/2017] [Indexed: 11/07/2022]
Abstract
The self-avoiding walk on the square site-diluted correlated percolation lattice is considered. The Ising model is employed to realize the spatial correlations of the metric space. As a well-accepted result, the (generalized) Flory's mean-field relation is tested to measure the effect of correlation. After exploring a perturbative Fokker-Planck-like equation, we apply an enriched Rosenbluth Monte Carlo method to study the problem. To be more precise, the winding angle analysis is also performed from which the diffusivity parameter of Schramm-Loewner evolution theory (κ) is extracted. We find that at the critical Ising (host) system, the exponents are in agreement with Flory's approximation. For the off-critical Ising system, we find also a behavior for the fractal dimension of the walker trace in terms of the correlation length of the Ising system ξ(T), i.e., D_{F}^{SAW}(T)-D_{F}^{SAW}(T_{c})∼1/sqrt[ξ(T)].
Collapse
Affiliation(s)
- J Cheraghalizadeh
- Department of Physics, University of Mohaghegh Ardabili, P.O. Box 179, Ardabil, Iran
| | - M N Najafi
- Department of Physics, University of Mohaghegh Ardabili, P.O. Box 179, Ardabil, Iran
| | - H Mohammadzadeh
- Department of Physics, University of Mohaghegh Ardabili, P.O. Box 179, Ardabil, Iran
| | - A Saber
- Department of Physics, University of Mohaghegh Ardabili, P.O. Box 179, Ardabil, Iran
| |
Collapse
|