Henrique Almeida Fernandes
Current Organizations / Institutes
The organizations/institutes from the scholar's latest article
Previous Organizations / Institutes
The organizations / institutes from the scholar's other articles in the past ten years
Research Subjects
Henrique Almeida Fernandes has a total of 26 co-authors. Here are the top co-authors by number of shared publications:
- Roberto da Silva
Universidade Federal do Rio Grande do Sul
20 shared publications - Jose Roberto Drugowich de Felicio
Departamento de Física e Matemática FFCLRP Universidade de São Paulo Avenida Bandeirantes 3900 014040-901 Ribeirão Preto São Paulo CEP Brazil
11 shared publications - Paulo Freitas Gomes
Grupo de Redes Complexas Aplicadas de Jataí Universidade Federal de Jataí Jataí Brazil
7 shared publications - Everaldo Arashiro
Universidade Federal do Rio Grande
3 shared publications - Alvaro de Almeida Caparica
Instituto de Física, Universidade Federal de Goiás, Goiânia, Brazil
2 shared publications Eduardo Velasco Stock
2 shared publications
Latest Publications
- Open Access
- 15 May 2026
- Published by Elsevier in Physica A: Statistical Mechanics and its Applications
- Vol. 695
Phase diagrams of the YK surface-reaction model on 2D lattices with exchange diffusion
- 15 May 2026
- Published by Elsevier in Physica A: Statistical Mechanics and its Applications
- Vol. 695
Open AccessPhase diagrams of the YK surface-reaction model on 2D lattices with exchange diffusion
- 2 April 2026
- Published by American Physical Society (APS) in Physical Review E
- Vol. 113 (4) , 044301
Number of spanning trees as an indicator of critical phenomena: When Kirchhoff meets Ising
Visibility graphs are spatial interpretations of time series. When derived from the time evolution of physical systems, the graphs associated with such series may exhibit properties that can reflect aspects such as ergodicity, criticality, or other dynamical behaviors. It is important to describe how the criticality of a system is manifested in the structure of the corresponding graphs or, in a particular way, in the spectra of certain matrices constructed from them. In this paper, we show how the critical behavior of an Ising spin system manifests in the spectra of the adjacency and Laplacian matrices constructed from an ensemble of time evolutions simulated via Monte Carlo (MC) Markov chains, even for small systems and short MC steps. In particular, we show that the number of spanning trees-or its logarithm-which represents a kind of structural entropy or topological complexity here obtained from Kirchhoff's theorem, can, in an alternative way, describe the criticality of the spin system. These findings parallel those obtained from the spectra of correlation matrices, which similarly encode signatures of critical and chaotic behavior.
- 2 April 2026
- Published by American Physical Society (APS) in Physical Review E
- Vol. 113 (4) , 044301
Number of spanning trees as an indicator of critical phenomena: When Kirchhoff meets Ising
Visibility graphs are spatial interpretations of time series. When derived from the time evolution of physical systems, the graphs associated with such series may exhibit properties that can reflect aspects such as ergodicity, criticality, or other dynamical behaviors. It is important to describe how the criticality of a system is manifested in the structure of the corresponding graphs or, in a particular way, in the spectra of certain matrices constructed from them. In this paper, we show how the critical behavior of an Ising spin system manifests in the spectra of the adjacency and Laplacian matrices constructed from an ensemble of time evolutions simulated via Monte Carlo (MC) Markov chains, even for small systems and short MC steps. In particular, we show that the number of spanning trees-or its logarithm-which represents a kind of structural entropy or topological complexity here obtained from Kirchhoff's theorem, can, in an alternative way, describe the criticality of the spin system. These findings parallel those obtained from the spectra of correlation matrices, which similarly encode signatures of critical and chaotic behavior.
- Open Access
Phase Diagrams of the YK Surface-Reaction Model on 2D lattices with Exchange Diffusion
In this work, we investigate the phase diagrams of the Yaldram and Khan catalytic surface model on square and hexagonal lattices when exchange diffusion is allowed for carbon monoxide (CO) and nitrogen (N) atoms. To reach our goal, we carried out steady-state Monte Carlo (MC) simulations over $4\times 10^5$ points, for both lattices, in order to obtain a framework of the steady reactive state of the model for different values of the nitric oxide dissociation rate, $r_{NO}$. The results show the emergence of steady reactive state for certain values of $r_{NO}$ and of exchange diffusion rate $x$ when the simulations take place on square lattices. Our findings on the hexagonal lattice also show that the diffusion of these species favors the appearance of the active phase for values of $r_{NO}$ lower than that found for the standard model. In addition, we observed that the system possesses a continuous phase transition and a discontinuous one separating the active phase from absorbing states for both lattices, except for $r_{NO}=1$ in which the continuous phase transition is destroyed and a steady reactive state emerges from the beginning since very small values of $x$.Open AccessPhase Diagrams of the YK Surface-Reaction Model on 2D lattices with Exchange Diffusion
In this work, we investigate the phase diagrams of the Yaldram and Khan catalytic surface model on square and hexagonal lattices when exchange diffusion is allowed for carbon monoxide (CO) and nitrogen (N) atoms. To reach our goal, we carried out steady-state Monte Carlo (MC) simulations over $4\times 10^5$ points, for both lattices, in order to obtain a framework of the steady reactive state of the model for different values of the nitric oxide dissociation rate, $r_{NO}$. The results show the emergence of steady reactive state for certain values of $r_{NO}$ and of exchange diffusion rate $x$ when the simulations take place on square lattices. Our findings on the hexagonal lattice also show that the diffusion of these species favors the appearance of the active phase for values of $r_{NO}$ lower than that found for the standard model. In addition, we observed that the system possesses a continuous phase transition and a discontinuous one separating the active phase from absorbing states for both lattices, except for $r_{NO}=1$ in which the continuous phase transition is destroyed and a steady reactive state emerges from the beginning since very small values of $x$. - 18 July 2025
- Published by Wiley in Journal of Computational Chemistry
- Vol. 46 (20) , e70182
Topology‐Induced Modifications in the Critical Behavior of the Yaldram–Khan Catalytic Reaction Model
In this work, we investigated how the use of complex networks as catalytic surfaces can affect the phase diagram of the Yaldram–Khan model, as well as how the order of the phase transitions present in the seminal work behaves when randomness is added to the model. The study was conducted by taking into consideration two well‐known random networks, the Erdős‐Rényi network (ERN), with its long‐range randomness, and the random geometric graph (RGG), with its spatially constrained randomness. We perform extensive steady‐state Monte Carlo simulations for , the NO dissociation rate, and show the behavior of the reactive window as a function of the average degree of the networks. Our results also show that, different from the ERN, which preserves the nature of the phase transitions of the original model for all considered average degrees, the RGG seems to have two second‐order phase transitions for small values of average degree.
- 18 July 2025
- Published by Wiley in Journal of Computational Chemistry
- Vol. 46 (20) , e70182
Topology‐Induced Modifications in the Critical Behavior of the Yaldram–Khan Catalytic Reaction Model
In this work, we investigated how the use of complex networks as catalytic surfaces can affect the phase diagram of the Yaldram–Khan model, as well as how the order of the phase transitions present in the seminal work behaves when randomness is added to the model. The study was conducted by taking into consideration two well‐known random networks, the Erdős‐Rényi network (ERN), with its long‐range randomness, and the random geometric graph (RGG), with its spatially constrained randomness. We perform extensive steady‐state Monte Carlo simulations for , the NO dissociation rate, and show the behavior of the reactive window as a function of the average degree of the networks. Our results also show that, different from the ERN, which preserves the nature of the phase transitions of the original model for all considered average degrees, the RGG seems to have two second‐order phase transitions for small values of average degree.
- Open Access
Revisiting the Contact Model with Diffusion Beyond the Conventional Methods
The contact process is a nonequilibrium Hamiltonian model that, even in one dimension, lacks an exact solution and has been extensively studied via Monte Carlo simulations, both in steady-state and time-dependent scenarios. Although the effects of particle mobility and diffusion on criticality have been preliminarily explored, they remain poorly understood in many aspects. In this work, we examine how the critical rate of the model varies with the probability of particle mobility. By analyzing different stochastic evolutions of the system, we employ two modern approaches: (1) Random Matrix Theory (RMT): By building on the success of RMT, particularly Wishart-like matrices, in studying statistical physics of systems with up-down symmetry via magnetization dynamics [R. da Silva, IJMPC 2022], we demonstrate its applicability to models with an absorbing state; (2) Optimized Temporal Power Laws: By using short-time dynamics, we optimize power laws derived from ensemble-averaged evolutions of the system. Both methods consistently reveal that the critical rate decays with mobility according to a simple Belehradek function. Additionally, a straightforward mean-field analysis supports the decay of the critical parameter with mobility, although it predicts a simpler linear dependence. We also demonstrate that the more sophisticated pair approximation mean-field model developed by ben-Avraham and Köhler aligns closely with the Belehradek function, precisely matching our lattice simulation results.
Open AccessRevisiting the Contact Model with Diffusion Beyond the Conventional Methods
The contact process is a nonequilibrium Hamiltonian model that, even in one dimension, lacks an exact solution and has been extensively studied via Monte Carlo simulations, both in steady-state and time-dependent scenarios. Although the effects of particle mobility and diffusion on criticality have been preliminarily explored, they remain poorly understood in many aspects. In this work, we examine how the critical rate of the model varies with the probability of particle mobility. By analyzing different stochastic evolutions of the system, we employ two modern approaches: (1) Random Matrix Theory (RMT): By building on the success of RMT, particularly Wishart-like matrices, in studying statistical physics of systems with up-down symmetry via magnetization dynamics [R. da Silva, IJMPC 2022], we demonstrate its applicability to models with an absorbing state; (2) Optimized Temporal Power Laws: By using short-time dynamics, we optimize power laws derived from ensemble-averaged evolutions of the system. Both methods consistently reveal that the critical rate decays with mobility according to a simple Belehradek function. Additionally, a straightforward mean-field analysis supports the decay of the critical parameter with mobility, although it predicts a simpler linear dependence. We also demonstrate that the more sophisticated pair approximation mean-field model developed by ben-Avraham and Köhler aligns closely with the Belehradek function, precisely matching our lattice simulation results.
- 1 April 2025
- Published by American Physical Society (APS) in Physical Review E
- Vol. 111 (4) , 044101
Unveiling the hidden weak universality of the Ziff-Gulari-Barshad model
In this work, we revisited the Ziff-Gullari-Barshad (ZGB) model in order to investigate its critical behavior when carbon monoxide (CO) molecules are allowed to desorb from the catalytic surface. As shown by several authors, when this kind of desorption takes place, the first-order phase transition of the standard model disappears, and an Ising-like critical point is found for a very small value of the desorption rate. However, our time-dependent Monte Carlo simulations reveal that, instead of a single critical point, there exists a critical line that appears to encompass multiple universality classes, passing through the Ising and three-state Potts points, as well as providing reasonable evidence for the four-state Potts point, resulting in an unprecedented critical line exhibiting weak universality.
- 1 April 2025
- Published by American Physical Society (APS) in Physical Review E
- Vol. 111 (4) , 044101
Unveiling the hidden weak universality of the Ziff-Gulari-Barshad model
In this work, we revisited the Ziff-Gullari-Barshad (ZGB) model in order to investigate its critical behavior when carbon monoxide (CO) molecules are allowed to desorb from the catalytic surface. As shown by several authors, when this kind of desorption takes place, the first-order phase transition of the standard model disappears, and an Ising-like critical point is found for a very small value of the desorption rate. However, our time-dependent Monte Carlo simulations reveal that, instead of a single critical point, there exists a critical line that appears to encompass multiple universality classes, passing through the Ising and three-state Potts points, as well as providing reasonable evidence for the four-state Potts point, resulting in an unprecedented critical line exhibiting weak universality.
- 11 May 2022
- Published by Elsevier in Physica A: Statistical Mechanics and its Applications
- Vol. 600
A physics-based algorithm to perform predictions in football leagues
- 11 May 2022
- Published by Elsevier in Physica A: Statistical Mechanics and its Applications
- Vol. 600
A physics-based algorithm to perform predictions in football leagues
- 23 March 2022
- Published by Elsevier in Physica A: Statistical Mechanics and its Applications
- Vol. 597
Topological transition in a coupled dynamics in random networks
- 23 March 2022
- Published by Elsevier in Physica A: Statistical Mechanics and its Applications
- Vol. 597
Topological transition in a coupled dynamics in random networks
- 1 May 2019
- Published by IOP Publishing in Journal of Statistical Mechanics: Theory and Experiment
- Vol. 2019 (5) , 053205
Mapping the phase transitions of the ZGB model with inert sites via nonequilibrium optimizing methods
In this paper, we revisit the ZGB model and explore the effects of the presence of inert sites on the catalytic surface. The continuous and discontinuous phase transitions of the model are studied via time-dependent Monte Carlo simulations. In our study, we are concerned with building a refinement procedure, based on a simple concept known as coefficient of determination r , in order to find the possible phase transition points given by the two parameters of the model: the adsorption rates of carbon monoxide, y , and the density of inert sites, . First, we obtain 10 6 values of r by sweeping the whole set of possible values of the parameters with an increment 10 −3 , i.e. and with . Then, with the possible phase transition points in hand, we turn our attention to some fixed values of and perform a more detailed refinement considering larger lattices and increasing the increment by one order of magnitude to estimate the critical points with higher precision. Finally, we estimate the static critical exponents , , and , as well as the dynamic critical exponents z and .
- 1 May 2019
- Published by IOP Publishing in Journal of Statistical Mechanics: Theory and Experiment
- Vol. 2019 (5) , 053205
Mapping the phase transitions of the ZGB model with inert sites via nonequilibrium optimizing methods
In this paper, we revisit the ZGB model and explore the effects of the presence of inert sites on the catalytic surface. The continuous and discontinuous phase transitions of the model are studied via time-dependent Monte Carlo simulations. In our study, we are concerned with building a refinement procedure, based on a simple concept known as coefficient of determination r , in order to find the possible phase transition points given by the two parameters of the model: the adsorption rates of carbon monoxide, y , and the density of inert sites, . First, we obtain 10 6 values of r by sweeping the whole set of possible values of the parameters with an increment 10 −3 , i.e. and with . Then, with the possible phase transition points in hand, we turn our attention to some fixed values of and perform a more detailed refinement considering larger lattices and increasing the increment by one order of magnitude to estimate the critical points with higher precision. Finally, we estimate the static critical exponents , , and , as well as the dynamic critical exponents z and .
- 1 April 2019
- Published by Elsevier in Physics Letters A
- Vol. 383 (12) , 1235-1241
Investigating the nonequilibrium aspects of long-range Potts model: Refinement of critical temperatures and raw exponents
- 1 April 2019
- Published by Elsevier in Physics Letters A
- Vol. 383 (12) , 1235-1241
Investigating the nonequilibrium aspects of long-range Potts model: Refinement of critical temperatures and raw exponents