Si Tang
Current Organizations / Institutes
The organizations/institutes from the scholar's latest article
Research Subjects
Si Tang has a total of 89 co-authors. Here are the top co-authors by number of shared publications:
- Stefano Allesina
University of Chicago
6 shared publications - Wei-Kuo Chen
University of Minnesota, Minneapolis, USA
5 shared publications Antonio Auffinger
2 shared publications- Antonio Auffinger
Department of Mathematics, Northwestern University, Evanston, USA
2 shared publications Dylan Hafner
Lehigh University
2 shared publications- Jonathan Touboul
Brandeis University
2 shared publications
Latest Publications
- Open Access
- 26 August 2026
- Published by American Physical Society (APS) in Physical Review E
- Vol. 114 (2) , 024408
Complexity and dynamics of partially symmetric random neural networks
Neural circuits exhibit structured connectivity, including an overrepresentation of reciprocal connections between neuron pairs. Despite important advances, a full understanding of how such partial symmetry in connectivity shapes neural dynamics remains elusive. Here we ask how correlations between reciprocal connections in a random, recurrent neural network affect phase-space complexity, defined as the exponential proliferation rate (with network size) of the number of fixed points that accompanies the transition to chaotic dynamics. We find a striking pattern: partial antisymmetry strongly amplifies complexity, while partial symmetry suppresses it. These opposing trends closely track changes in other measures of dynamical behavior, such as dimensionality, Lyapunov exponents, and transient path length, supporting the view that fixed-point structure is a key determinant of network dynamics. Thus, positive reciprocal correlations favor low-dimensional, slowly varying activity, whereas negative correlations promote high-dimensional, rapidly fluctuating chaotic activity. These results yield testable predictions about how connection reciprocity may shape neural dynamics and function.
- 26 August 2026
- Published by American Physical Society (APS) in Physical Review E
- Vol. 114 (2) , 024408
Open AccessComplexity and dynamics of partially symmetric random neural networks
Neural circuits exhibit structured connectivity, including an overrepresentation of reciprocal connections between neuron pairs. Despite important advances, a full understanding of how such partial symmetry in connectivity shapes neural dynamics remains elusive. Here we ask how correlations between reciprocal connections in a random, recurrent neural network affect phase-space complexity, defined as the exponential proliferation rate (with network size) of the number of fixed points that accompanies the transition to chaotic dynamics. We find a striking pattern: partial antisymmetry strongly amplifies complexity, while partial symmetry suppresses it. These opposing trends closely track changes in other measures of dynamical behavior, such as dimensionality, Lyapunov exponents, and transient path length, supporting the view that fixed-point structure is a key determinant of network dynamics. Thus, positive reciprocal correlations favor low-dimensional, slowly varying activity, whereas negative correlations promote high-dimensional, rapidly fluctuating chaotic activity. These results yield testable predictions about how connection reciprocity may shape neural dynamics and function.
- Open Access
Complexity and dynamics of partially symmetric random neural networks
Neural circuits exhibit structured connectivity, including an overrepresentation of reciprocal connections between neuron pairs. Despite important advances, a full understanding of how such partial symmetry in connectivity shapes neural dynamics remains elusive. Here we ask how correlations between reciprocal connections in a random, recurrent neural network affect phase-space complexity, defined as the exponential proliferation rate (with network size) of the number of fixed points that accompanies the transition to chaotic dynamics. We find a striking pattern: partial anti-symmetry strongly amplifies complexity, while partial symmetry suppresses it. These opposing trends closely track changes in other measures of dynamical behavior, such as dimensionality, Lyapunov exponents, and transient path length, supporting the view that fixed-point structure is a key determinant of network dynamics. Thus, positive reciprocal correlations favor low-dimensional, slowly varying activity, whereas negative correlations promote high-dimensional, rapidly fluctuating chaotic activity. These results yield testable predictions about the link between connection reciprocity, neural dynamics and function.
Open AccessComplexity and dynamics of partially symmetric random neural networks
Neural circuits exhibit structured connectivity, including an overrepresentation of reciprocal connections between neuron pairs. Despite important advances, a full understanding of how such partial symmetry in connectivity shapes neural dynamics remains elusive. Here we ask how correlations between reciprocal connections in a random, recurrent neural network affect phase-space complexity, defined as the exponential proliferation rate (with network size) of the number of fixed points that accompanies the transition to chaotic dynamics. We find a striking pattern: partial anti-symmetry strongly amplifies complexity, while partial symmetry suppresses it. These opposing trends closely track changes in other measures of dynamical behavior, such as dimensionality, Lyapunov exponents, and transient path length, supporting the view that fixed-point structure is a key determinant of network dynamics. Thus, positive reciprocal correlations favor low-dimensional, slowly varying activity, whereas negative correlations promote high-dimensional, rapidly fluctuating chaotic activity. These results yield testable predictions about the link between connection reciprocity, neural dynamics and function.
- Open AccessOn the Fitness Landscape in the $NK$ ModelThe $NK$ model, introduced by Kauffman, Levin, and Weinberger, is a random field used to describe the fitness landscape of certain species with $N$ genetic loci, each interacting with $K$ others. The model has wide applications in understanding evolutionary and natural selection as it captures ruggedness feature of the fitness landscape. Earlier literature has been focused on the case $K$ being a fixed positive integer and used tools from Ergodic and Markov theory. In this paper, by viewing it as a statistical physics object, we investigate the $NK$ model in the regime $K/N\toα\in(0,1]$ via the spin glass methodologies. Our main result identifies the exact limits for the free energy at any temperature and the maximum fitness. Moreover, we show that the $NK$ model exhibits a multiple-peak structure, namely, the number of near-fittest genomes that are asymptotically orthogonal to each other is exponentially large. Based on establishing the overlap gap properties, we obtain quantitative descriptions for the geometry of the fitness landscape and deduce that, in particular, near-fittest evolutionary paths become impossible as the fitness levels of the genomes approach the global maximum for any $α\in (0,1]$. Nevertheless, we also show that by choosing $α$ sufficiently small, an evolutionary path maintained at a given fitness level can be constructed with high probability.Open AccessOn the Fitness Landscape in the $NK$ ModelThe $NK$ model, introduced by Kauffman, Levin, and Weinberger, is a random field used to describe the fitness landscape of certain species with $N$ genetic loci, each interacting with $K$ others. The model has wide applications in understanding evolutionary and natural selection as it captures ruggedness feature of the fitness landscape. Earlier literature has been focused on the case $K$ being a fixed positive integer and used tools from Ergodic and Markov theory. In this paper, by viewing it as a statistical physics object, we investigate the $NK$ model in the regime $K/N\toα\in(0,1]$ via the spin glass methodologies. Our main result identifies the exact limits for the free energy at any temperature and the maximum fitness. Moreover, we show that the $NK$ model exhibits a multiple-peak structure, namely, the number of near-fittest genomes that are asymptotically orthogonal to each other is exponentially large. Based on establishing the overlap gap properties, we obtain quantitative descriptions for the geometry of the fitness landscape and deduce that, in particular, near-fittest evolutionary paths become impossible as the fitness levels of the genomes approach the global maximum for any $α\in (0,1]$. Nevertheless, we also show that by choosing $α$ sufficiently small, an evolutionary path maintained at a given fitness level can be constructed with high probability.
- Open Access
On the time constant of high dimensional first passage percolation, revisited
In [2], it was claimed that the time constant $μ_{d}(e_{1})$ for the first-passage percolation model on $\mathbb Z^{d}$ is $μ_{d}(e_{1}) \sim \log d/(2ad)$ as $d\to \infty$, if the passage times $(τ_{e})_{e\in \mathbb E^{d}}$ are i.i.d., with a common c.d.f. $F$ satisfying $\left|\frac{F(x)}{x}-a\right| \le \frac{C}{|\log x|}$ for some constants $a, C$ and sufficiently small $x$. However, the proof of the upper bound, namely, Equation (2.1) in [2] \begin{align} \limsup_{d\to\infty} \frac{μ_{d}(e_{1})ad}{\log d} \le \frac{1}{2} \end{align} is incorrect. In this article, we provide a different approach that establishes this inequality. As a side product of this new method, we also show that the variance of the non-backtracking passage time to the first hyperplane is of order $o\big((\log d/d)^{2}\big)$ as $d\to \infty$ in the case of the when the edge weights are exponentially distributed.Open AccessOn the time constant of high dimensional first passage percolation, revisited
In [2], it was claimed that the time constant $μ_{d}(e_{1})$ for the first-passage percolation model on $\mathbb Z^{d}$ is $μ_{d}(e_{1}) \sim \log d/(2ad)$ as $d\to \infty$, if the passage times $(τ_{e})_{e\in \mathbb E^{d}}$ are i.i.d., with a common c.d.f. $F$ satisfying $\left|\frac{F(x)}{x}-a\right| \le \frac{C}{|\log x|}$ for some constants $a, C$ and sufficiently small $x$. However, the proof of the upper bound, namely, Equation (2.1) in [2] \begin{align} \limsup_{d\to\infty} \frac{μ_{d}(e_{1})ad}{\log d} \le \frac{1}{2} \end{align} is incorrect. In this article, we provide a different approach that establishes this inequality. As a side product of this new method, we also show that the variance of the non-backtracking passage time to the first hyperplane is of order $o\big((\log d/d)^{2}\big)$ as $d\to \infty$ in the case of the when the edge weights are exponentially distributed. - Open Access
- 1 January 2025
- Published by Institute of Mathematical Statistics in Electronic Journal of Probability
- Vol. 30 (none) , 1-13
On the time constant of high dimensional first passage percolation, revisited
In [2], it was claimed that the time constant μd(e1) for the first-passage percolation model on Zd is μd(e1)∼logd∕(2ad) as d→∞, if the passage times (τe)e∈Ed are i.i.d., with a common c.d.f. F satisfying F(x) x−a≤C |logx| for some constants a,C and sufficiently small x. However, the proof of the upper bound, namely, Equation (2.1) in [2] lim supd→∞μd(e1)ad logd≤1 2(0.1) is incorrect. In this article, we provide a different approach that establishes (0.1). As a side product of this new method, we also show that the variance of the non-backtracking passage time to the first hyperplane is of order o((logd∕d)2) as d→∞ in the case of the when the edge weights are exponentially distributed.
- 1 January 2025
- Published by Institute of Mathematical Statistics in Electronic Journal of Probability
- Vol. 30 (none) , 1-13
Open AccessOn the time constant of high dimensional first passage percolation, revisited
In [2], it was claimed that the time constant μd(e1) for the first-passage percolation model on Zd is μd(e1)∼logd∕(2ad) as d→∞, if the passage times (τe)e∈Ed are i.i.d., with a common c.d.f. F satisfying F(x) x−a≤C |logx| for some constants a,C and sufficiently small x. However, the proof of the upper bound, namely, Equation (2.1) in [2] lim supd→∞μd(e1)ad logd≤1 2(0.1) is incorrect. In this article, we provide a different approach that establishes (0.1). As a side product of this new method, we also show that the variance of the non-backtracking passage time to the first hyperplane is of order o((logd∕d)2) as d→∞ in the case of the when the edge weights are exponentially distributed.
- Open Access
On the TAP equations via the cavity approach in the generic mixed $p$-spin models
In 1977, Thouless, Anderson, and Palmer (TAP) derived a system of consistent equations in terms of the effective magnetization in order to study the free energy in the Sherrington-Kirkpatrick (SK) spin glass model. The solutions to their equations were predicted to contain vital information about the landscapes in the SK Hamiltonian and the TAP free energy and moreover have direct connections to Parisi's replica ansatz. In this work, we aim to investigate the validity of the TAP equations in the generic mixed $p$-spin model. By utilizing the ultrametricity of the overlaps, we show that the TAP equations are asymptotically satisfied by the conditional local magnetizations on the asymptotic pure states.
Open AccessOn the TAP equations via the cavity approach in the generic mixed $p$-spin models
In 1977, Thouless, Anderson, and Palmer (TAP) derived a system of consistent equations in terms of the effective magnetization in order to study the free energy in the Sherrington-Kirkpatrick (SK) spin glass model. The solutions to their equations were predicted to contain vital information about the landscapes in the SK Hamiltonian and the TAP free energy and moreover have direct connections to Parisi's replica ansatz. In this work, we aim to investigate the validity of the TAP equations in the generic mixed $p$-spin model. By utilizing the ultrametricity of the overlaps, we show that the TAP equations are asymptotically satisfied by the conditional local magnetizations on the asymptotic pure states.
- 18 March 2024
- Published by Springer Nature in Communications in Mathematical Physics
- Vol. 405 (4) , 1-43
On the TAP Equations via the Cavity Approach in the Generic Mixed p-Spin Models
- 18 March 2024
- Published by Springer Nature in Communications in Mathematical Physics
- Vol. 405 (4) , 1-43
On the TAP Equations via the Cavity Approach in the Generic Mixed p-Spin Models
- 13 May 2021
- Published by Springer Nature in Communications in Mathematical Physics
- Vol. 386 (2) , 1209-1242
On Convergence of the Cavity and Bolthausen’s TAP Iterations to the Local Magnetization
- 13 May 2021
- Published by Springer Nature in Communications in Mathematical Physics
- Vol. 386 (2) , 1209-1242
On Convergence of the Cavity and Bolthausen’s TAP Iterations to the Local Magnetization
- Open Access
On convergence of the cavity and Bolthausen's TAP iterations to the local magnetization
The cavity and TAP equations are high-dimensional systems of nonlinear equations of the local magnetization in the Sherrington-Kirkpatrick model. In the seminal work [Comm. Math. Phys., 325(1):333-366, 2014], Bolthausen introduced an iterative scheme that produces an asymptotic solution to the TAP equations if the model lies inside the Almeida-Thouless transition line. However, it was unclear if this asymptotic solution coincides with the local magnetization. In this work, motivated by the cavity equations, we introduce a new iterative scheme and establish a weak law of large numbers. We show that our new scheme is asymptotically the same as the so-called Approximate Message Passing algorithm, a generalization of Bolthausen's iteration, that has been popularly adapted in compressed sensing, Bayesian inferences, etc. Based on this, we confirm that our cavity iteration and Bolthausen's scheme both converge to the local magnetization as long as the overlap is locally uniformly concentrated.
Open AccessOn convergence of the cavity and Bolthausen's TAP iterations to the local magnetization
The cavity and TAP equations are high-dimensional systems of nonlinear equations of the local magnetization in the Sherrington-Kirkpatrick model. In the seminal work [Comm. Math. Phys., 325(1):333-366, 2014], Bolthausen introduced an iterative scheme that produces an asymptotic solution to the TAP equations if the model lies inside the Almeida-Thouless transition line. However, it was unclear if this asymptotic solution coincides with the local magnetization. In this work, motivated by the cavity equations, we introduce a new iterative scheme and establish a weak law of large numbers. We show that our new scheme is asymptotically the same as the so-called Approximate Message Passing algorithm, a generalization of Bolthausen's iteration, that has been popularly adapted in compressed sensing, Bayesian inferences, etc. Based on this, we confirm that our cavity iteration and Bolthausen's scheme both converge to the local magnetization as long as the overlap is locally uniformly concentrated.
- Open Access
On the minimal drift for recurrence in the frog model on $d$-ary trees
We study the recurrence of one-per-site frog model $\text{FM}(d, p)$ on a $d$-ary tree with drift parameter $p\in [0,1]$, which determines the bias of frogs' random walks. We are interested in the minimal drift $p_{d}$ so that the frog model is recurrent. Using a coupling argument together with a generating function technique, we prove that for all $d \ge 2$, $p_{d}\le 1/3$, which is the optimal universal upper bound.Open AccessOn the minimal drift for recurrence in the frog model on $d$-ary trees
We study the recurrence of one-per-site frog model $\text{FM}(d, p)$ on a $d$-ary tree with drift parameter $p\in [0,1]$, which determines the bias of frogs' random walks. We are interested in the minimal drift $p_{d}$ so that the frog model is recurrent. Using a coupling argument together with a generating function technique, we prove that for all $d \ge 2$, $p_{d}\le 1/3$, which is the optimal universal upper bound.