Postdoctoral Researcher · Centro de Ciencias de la Complejidad, UNAM

Saúl Huitzil

Mexican physicist studying complex systems — how criticality shows up in time series, how modularity emerges in living systems, and how collective motion arises.

  • Complex systems
  • Criticality & time series
  • Complex networks & modularity
  • Microbiome & evolution
  • Artificial intelligence
Portrait of Saúl Huitzil outdoors in the snow
Ph.D. in Physics, UNAM
Ph.D.Physics, UNAM · 2021, with honorable mention
5publications & book chapters, 2 more in review or in preparation
4interactive models you can run below

01 — About

Modelling complex systems across scales

I am a Mexican physicist with interdisciplinary training in the study of complex networks, evolution and critical dynamics. I did my B.Sc. in Physics at the Benemérita Universidad Autónoma de Puebla, followed by an M.Sc. and a Ph.D. in Sciences (Physics) at UNAM, where I developed models to understand how the interactions between organisms and their microbiomes drive biological adaptability and complexity.

As a postdoctoral researcher at Northwestern University (2022–2024) I widened my focus to the emergence of modularity, hierarchical organisation and collective motion, going deeper into how living systems structure their complexity across scales. Today, at the Centro de Ciencias de la Complejidad (UNAM), my research centres on time-series analysis to detect and characterise critical states in complex systems, and on their relation to optimisation processes.

Along the way I have built skills in mathematical modelling, artificial intelligence, computational simulation and the analysis of complex data.

02 — Research in the browser

Interactive models

Minimal models of collective behaviour, network structure and gene regulation, written in p5.js and running live in this page. They load as you scroll, so the text stays fast on slow connections.

Boolean network on its way to an attractor

Criticality Gene networks

A random Boolean network is a minimal model of a gene regulatory circuit. Each node is a gene that is either on or off (expressed or not), with k randomly chosen regulators and a random Boolean rule that decides its next state from theirs. Every gene updates at the same time, one step per frame.

The graph on the left is the network: bright nodes are genes that are on, faint ones are off, and the edges are regulatory interactions. The grid on the right is the state trajectory — one row per time step, one column per gene, green for on. After an irregular transient the same block of rows starts repeating: that cycle is an attractor, the set of states the network settles into.

The numbers D(s0, s1) … D(s3, s4) average, over many random starting states, the fraction of genes that flip at each step. They are a quick read on the dynamical regime: values near zero describe a frozen, ordered network, values near 0.5 mean that half the genes change at every step, as in the chaotic regime. The k slider moves you between those regimes — k = 1 tends to freeze, k = 2 is the classic critical case, and larger k becomes chaotic.

Runs live in this page
Open full screen Controls: Make Network builds a new random network, Change state picks a new initial state, and the k and N sliders set the number of regulators per gene and the network size (5 × the slider value).

The Vicsek model

Collective motion

The Vicsek model simulates the collective behaviour of self-propelled particles — birds or fish in a flock, or cells in a tissue. Each particle updates its velocity to the average velocity of its neighbours within a certain distance, and then moves in the updated direction.

It matters for living systems because it shows how coordinated behaviour can emerge in a group without a leader or central control, and it helps explain how cells in a tissue coordinate to form patterns during development.

Runs live from the p5.js editor

Force-directed graph drawing

Network science

Force-directed graph drawing visualises complex network structures — social, transport or biological. Nodes are treated as particles and edges as springs, and the algorithm simulates the forces between them (Hooke's law among others) until the energy of the system is minimised, producing a layout that is easy to read.

It is useful because it uncovers hidden patterns in large networks: groups of highly connected nodes (communities or clusters), network hierarchies, and the key players that hold a structure together — often in a way non-experts can interpret at a glance.

Runs live from the p5.js editor

Active elastic model

Active matter

The active elastic model is a minimal description of self-propelled agents that interact through attraction, repulsion and alignment using elastic interactions. It has a simple mechanical realisation and applies to real systems such as active cell membranes, robotic swarms or animal groups.

Agents are connected to their neighbours by linear springs placed a distance R in front of their centres of rotation. Depending on R, the elastic interactions produce mainly attraction–repulsion or mainly alignment.

Runs live from the p5.js editor

03 — Projects & reports

Longer pieces of work

Standalone reports and collaborative projects, hosted here or on their own sites.

Work in progress Network science

A citation map of research on modularity in biology

207 hand-picked papers on modularity, each scored on five ways of thinking about modules, together with the thousands of papers they cite or are cited by. The report sets out how the corpus was built and includes an interactive pentagon map where every dot is one paper.

Work in progress with Joseph L. Hellerstein (University of Washington), Brandilyn Stigler (Southern Methodist University), Cristián Huepe (Northwestern University), Elena Dimitrova (Cal Poly) and Virginia Pasour, among many others.

Interactive report · 207 seed papers · 14,315 papers · 21,276 links

Open the report
Active matter

Multiscale Organizational Principles in Active Matter and Living Systems

A project led by Cristián Huepe at Northwestern and supported by the John Templeton Foundation, on how multiscale structures and dynamics produce the self-organising complexity of active matter and living systems — combining agent-based simulations, data and theory with evolutionary models of Boolean and adaptive networks.

I was part of the team during my postdoc at Northwestern (2022–2024). The project site gathers its animations and news.

Visit the project

04 — Publications

Selected work

Five published papers and book chapters, plus two manuscripts under review or in preparation.

Under review & in preparation

In review

Criticality detection in multi-variable systems through matrix scale invariance

Manuscript under review.

In prep.

Emergence of modularity in hierarchical living systems

Manuscript in preparation.

Full and up-to-date list on Google Scholar

05 — Background

Career, education & skills

Education

  • 2021 Ph.D. in Science (Physics) National Autonomous University of Mexico (UNAM) · passed with honorable mention
  • 2016 M.Sc. in Physics National Autonomous University of Mexico (UNAM)
  • 2013 B.Sc. in Physics Benemérita Universidad Autónoma de Puebla (BUAP)

Positions

  • 2025–2026 Postdoctoral Researcher Centro de Ciencias de la Complejidad (C3), UNAM · criticality and time series in complex systems
  • 2022–2024 Postdoctoral Researcher Northwestern University, USA · modularity, hierarchical organisation and collective motion

Research interests

  • Complex networks
  • Modularity
  • Criticality
  • Time-series analysis
  • Microbiome
  • Gene regulation
  • Heritability
  • Artificial intelligence

Teaching

  • Teaching Assistant
    • 2021–2022 Introduction to complex systems, UNAM
    • 2017–2018 Complex variables, UAEM
    • 2017–2018 Differential equations, UAEM
    • 2016–2017 Statistical physics, UNAM

Talks & conferences

  • 2024Biology Across Scales symposiumPortugal
  • 2023Complex NetworksFrance
  • 2023New Perspectives in Active SystemsGermany
  • 2022APS March MeetingUnited States
  • 20192nd International Conference on HolobiontsCanada
  • 2017Conference on Complex SystemsMexico

Peer review & recognition

  • Reviewer for Frontiers Peer review for Frontiers journals.
  • 2022/2023 Rising Star in Coevolution Publication selected by the Frontiers editors for Rising Stars in Coevolution, which showcases early-career scientists.
  • Co-organiser, Student Seminar Institute of Physical Sciences, UNAM.

Programming

  • Python
  • JavaScript p5.js
  • C#
  • Processing
  • HTML & CSS

Languages

07 — Contact

Let's talk about complex systems, data or collaborations

Postdoctoral Researcher, Centro de Ciencias de la Complejidad (C3), UNAM · the fastest way to reach me is email.

saulhuitzil@gmail.com