A computational model of cerebral cortex folding

The geometric complexity and variability of the human cerebral cortex has long intrigued the scientific community. As a result, quantitative description of cortical folding patterns and the understanding of underlying folding mechanisms have emerged as important research goals. This paper presents a computational 3-dimensional geometric model of cerebral cortex folding initialized by MRI data of a human fetal brain and deformed under the governance of a partial differential equation modeling cortical growth. By applying different simulation parameters, our model is able to generate folding convolutions and shape dynamics of the cerebral cortex. The simulations of this 3D geometric model provide computational experimental support to the following hypotheses: 1) Mechanical constraints of the skull regulate the cortical folding process. 2) The cortical folding pattern is dependent on the global cell growth rate of the whole cortex. 3) The cortical folding pattern is dependent on relative rates of cell growth in different cortical areas. 4) The cortical folding pattern is dependent on the initial geometry of the cortex.