Projects

A collection of my research projects.

Physical networks

Physical networks are networks of tangible components, in which nodes and edges are embedded in physical space and subject to constraints such as volume exclusion [(Dehmamy et al. 2018)]. Examples include biological neural networks, vascular systems, porous media, and granular materials. Advances in imaging technologies and accelerated 3D reconstruction by ML/AI have enabled to obtain detailed three-dimensional maps of such complex systems, providing new opportunities to study general organizing principles for tangible complex systems.

While the physical layout of networks has been shown to give rise to emergent structural properties, such as entanglement [Glover and Barabási (2024)], bundling (Bonamassa et al. 2025), jamming (Pósfai et al. 2024), and correlations between node degree and node volume (Pete et al. 2024), it remains elusive how these would impact dynamical processes that unfold on physical networks.

For my PhD thesis, I am exploring the spectral and topological properties of physical networks and their dynamical implications through physical Laplacian, a vertex-weighted graph Laplacian for physical networks derived in (Pete et al. 2024) as \begin{align} \mathbf{Q}_{\textrm{P}} = \mathbf{V}^{-1/2} \mathbf{Q}_{\textrm{G}} \mathbf{V}^{-1/2} \end{align} where \mathbf{G_{\textrm{G}}} is the combinatorial graph Laplacian and \mathbf{V}=\mathrm{diag}(v_1, \dots, v_{N}) is the diagonal volume matrix.

NoteLocalization

Coming soon…

TipSynchronization

Coming soon…

Higher-order networks with triadic interactions

Triadic interactions are higher-order, three-body interactions in which a node modulates edge between two other nodes. Such interactions are reported in various real-world systems; for example, glia cells in the brain regulate synaptic connectivity between two neurons, either inhibiting or exciting synaptic signaling. To model such interactions, we defined triadic interactions on network G=(V,E) as an interaction between a node i \in V and an edge \ell \in E encoded by incidence matrix \mathbf{K} \in \{-1, 0, 1\}^{|E| \times |V|} with entries \begin{align} K_{\ell i} = \begin{cases} -1 & \text{if node $i$ inhibits edge $\ell$}, \\ 1 & \text{if node $i$ excites edge $\ell$}, \\ 0 & \text{otherwise}. \end{cases} \end{align} Signed triadic interactions turn percolation processes into dynamical ones, with phenomena such as a period doubling and chaos [(Sun et al. 2023), (Sun and Bianconi 2024), (Millán et al. 2024)], and dynamical implications of triadic interactions remain an important but open question.

NoteLangevin dynamics on networks with triadic interactions

To investigate how triadic interactions affect dynamics, we studied Langevin dynamics on networks incorporating triadic interactions \begin{align} d\mathbf{X} = - (\mathbf{L}^{\textrm{(T)}} + \alpha \mathbf{I}) \mathbf{X} dt + \mathbf{\Gamma} d\mathbf{W}_{t} \end{align} where d\mathbf{W}_{t} is the |V|-dimensional Wiener process and \mathbf{L}^{\textrm{(T)}} is the triadic Laplacian \begin{align} L_{ij}^{\textrm{(T)}} (\mathbf{X}) &= \left(\sum_{k} J_{i k} (\mathbf{X}) \right) \delta_{ij} - J_{ij} (\mathbf{X})\\ J_{ij}(\mathbf{X}; \hat{T}) &= w_{-} + (w_{+} - w_{-}) \theta\left(\sum_{k} K_{[i,j] k} X_{k} - \hat{T} \right) \end{align} that coevolves with \mathbf{X}. Excitation and inhibitation of triadic interactions are triggered when the states of triadically incident nodes collectively exceed a threshold \hat{T}. Our results show that triadic interactions in the model are observable in conditional correlations and conditional mutual information among the three participating nodes, revealing both the strength and direction of the triadic interactions.

We applied this framework as a null model to identify potential triadic interactions in gene regulatory networks.

For details, see (Niedostatek et al. 2025).

Fractal scale-free networks

A network is fractal with respect to shortest-path lengths, if the minimum number of subgraphs (boxes) of diameter l_{\textup{B}} required to cover the network scales as \begin{equation} N_{\textup{B}}(l_{\textup{B}}) \sim l_{\textup{B}}^{-D_{\textup{f}}} \end{equation} where D_{\textup{f}} is the fractal dimension of the network (Song et al. 2005). In fractal networks, the average path lengths also scales as the power of the network size, \begin{equation} \langle l \rangle \sim N^{1/D_{\textup{f}}}. \end{equation} While many empirical networks exhibit the small-world property, where the average path lengths scales logarithmically with the number of nodes, the fractal property can still hold at length scales shorter than the network’s characteristic size. Detecting this scaling through shortest paths, however, is often challenging as most real networks do not have a sufficiently large diameter.

Fractal networks are of particular interest due to their well-defined structural scale invariance, which gives rise to distinctive features such as long-range degree correlations manifested as hub repulsion (Rozenfeld et al. 2008). Fractal network models therefore provide a natural analytical framework for studying network coarse-graining and renormalization.

NoteBifractality

When scale-free and fractal properties coexist, a broad class of networks exhibits bifractality, characterized by two distinct local fractal dimensions, d_{\textup{f}}^{\min} and d_{\textup{f}}^{\max} that depend on position within the network (Yamamoto and Yakubo 2023). In particular, consider a fractal scale-free network in which, under renormalization, the number of nodes \nu_{\textup{B}} in a supernode (subgraph) of diameter l_{\textup{B}} is proportional to the degree k_{\textup{B}} of that supernode (i.e., the number of adjacent subgraphs): \begin{equation} \nu_{\textup{B}} \propto k_{\textup{B}}. \end{equation} Under this condition, the network exhibits two local fractal dimensions (equivalently, Hölder exponents): \begin{equation} d_{\textup{f}}^{\max} = D_{\textup{f}}, \quad d_{\textup{f}}^{\min} = D_{\textup{f}} \left(\frac{\gamma-2}{\gamma-1}\right), \end{equation} where \gamma is the exponent of the degree distribution and D_{\textup{f}} is the global fractal dimension of the network.

Extensive analytical and numerical calculations for both deterministic and stochastic models suggest that this proportionality condition holds broadly across fractal scale-free networks. We therefore conjecture that bifractality is a general property of fractal scale-free networks.

For details, see (Yamamoto and Yakubo 2023).

TipRandom walks on bifractal networks

Building on the bifractality conjecture, we studied the walk and spectral dimensions of random walks on fractal scale-free networks. We found that the walk dimension d_{\textup{w}}, which characterizes how the mean topological displacement of random walkers scales, remains constant across the network and does not depend on position.

In contrast, the spectral dimension, which governs the scaling of the return probability of a random walker, splits into two distinct values, d_{\textup{s}}^{\min} and d_{\textup{s}}^{\max}. This positional dependence arises directly from the network’s bifractal structure.

For details, see (Yakubo et al. 2024).

Back to top

References

Bonamassa, Ivan, Balázs Ráth, Márton Pósfai, et al. 2025. “Logarithmic Kinetics and Bundling in Random Packings of Elongated 3D Physical Links.” Proceedings of the National Academy of Sciences 122 (32): e2427145122.
Dehmamy, Nima, Soodabeh Milanlouei, and Albert-László Barabási. 2018. “A Structural Transition in Physical Networks.” Nature 563 (7733): 676–80.
Glover, Cory, and Albert-László Barabási. 2024. “Measuring Entanglement in Physical Networks.” Physical Review Letters 133 (7): 077401.
Liu, Yanchen, Nima Dehmamy, and Albert-László Barabási. 2021. “Isotopy and Energy of Physical Networks.” Nature Physics 17 (2): 216–22.
Millán, Ana P, Hanlin Sun, Joaquı́n J Torres, and Ginestra Bianconi. 2024. “Triadic Percolation Induces Dynamical Topological Patterns in Higher-Order Networks.” PNAS Nexus 3 (7): pgae270.
Niedostatek, Marta, Anthony Baptista, Jun Yamamoto, et al. 2025. “Mining Higher-Order Triadic Interactions.” Nature Communications.
Pete, Gábor, Ádám Timár, Sigurdur Örn Stefánsson, Ivan Bonamassa, and Márton Pósfai. 2024. “Physical Networks as Network-of-Networks.” Nature Communications 15 (1): 4882.
Pósfai, Márton, Balázs Szegedy, Iva Bačić, et al. 2024. “Impact of Physicality on Network Structure.” Nature Physics 20 (1): 142–49.
Rozenfeld, Hernán D, Lazaros K Gallos, Chaoming Song, and Hernán A Makse. 2008. “Fractal and Transfractal Scale-Free Networks.” arXiv Preprint arXiv:0808.2206.
Song, Chaoming, Shlomo Havlin, and Hernan A Makse. 2005. “Self-Similarity of Complex Networks.” Nature 433 (7024): 392–95.
Sun, Hanlin, and Ginestra Bianconi. 2024. “Higher-Order Triadic Percolation on Random Hypergraphs.” Physical Review E 110 (6): 064315.
Sun, Hanlin, Filippo Radicchi, Jürgen Kurths, and Ginestra Bianconi. 2023. “The Dynamic Nature of Percolation on Networks with Triadic Interactions.” Nature Communications 14 (1): 1308.
Yakubo, Kousuke, Gentaro Shimojo, and Jun Yamamoto. 2024. “Random Walks on Bifractal Networks.” Phys. Rev. E 110 (December): 064318. https://doi.org/10.1103/PhysRevE.110.064318.
Yamamoto, Jun, and Kousuke Yakubo. 2023. “Bifractality of Fractal Scale-Free Networks.” Phys. Rev. E 108: 024302. https://doi.org/10.1103/PhysRevE.108.024302.