  Many cognitive processes, including working memory, recruit multiple
distributed interacting brain regions to encode information. How to understand
the underlying cognition function mechanism of working memory is a challenging
problem, which involves neural circuit configuration from multiple brain
regions as well as stochastic transition dynamics between brain states. The
energy landscape idea provides a tool to study the global stability and
stochastic transition dynamics in the distributed cognitive function system.
However, how to quantify the energy landscape in a realistic large-scale brain
network remains unclear. Here, based on an anatomically constrained
computational model of large-scale macaque cortex, we quantified the underlying
multistable attractor landscape of distributed working memory. In the absence
of external stimulation, the landscape exhibits three stable attractors, a
spontaneous state, and two memory states. In the attractor landscape framework,
the working memory function is governed by the change of landscape topography
and the switch of system state according to the task requirement. The barrier
height inferred from landscape topography quantifies the global stability of
memory state and robustness to non-selective random fluctuations and distractor
stimuli. The kinetic transition path identified by the minimum action path
approach reveals that the spontaneous state serves as an intermediate state
during the switch between the two memory states, the memory stored in the
cortical area with higher hierarchy is more stable, and information flow
follows the direction of hierarchical structure. These results provide new
insights into the underlying mechanism of distributed working memory function,
and the landscape and kinetic path approach can be applied to other cognitive
function-related problems in brain networks.
