abstract: As an interdisciplinary discipline, data mining (DM) is popular in education area especially when examining students' learning performances. It focuses on analyzing educational related data to develop models for improving learners' learning experiences and enhancing institutional effectiveness. Therefore, DM does help education institutions provide high-quality education for its learners. Applying data mining in education also known as educational data mining (EDM), which enables to better understand how students learn and identify how improve educational outcomes. Present paper is designed to justify the capabilities of data mining approaches in the field of education. The latest trends on EDM research are introduced in this review. Several specific algorithms, methods, applications and gaps in the current literature and future insights are discussed here.
author:
title: Data-Mining Research in Education

# INTRODUNCTION

One of the biggest challenges that educational institutions facing today is the exponential growth of educational data and how to apply this data to improve the quality of managerial decisions \[1\]. Education Institutions would like to know, for instance, which students will enroll in particular course programs, and which students will need assistance for graduation. Through the analysis and presentation of data they collected, or data mining, the challenges of these student or learners are able be effectively addressed. Data mining enables organizations to uncover and understand hidden patterns in vast databases by using their current reporting capabilities. And these patterns are then built into data mining models and applied to predict individual behavior and performance with high accuracy. In this way, resources and staff can be allocated by institutions more effectively. Data mining may also, for example, efficiently allocate resources with an accurate estimate of how many students will take action before he or she drops out.

Educational data mining (EDM) is an emerging discipline including but not limited to information retrieval, recommender systems, visual data analytics, social network analysis (SNA), cognitive psychology, psychometrics, and so on. Its methods is often different from those methods from the broader data mining literature. What's more, EDM draws from several reference disciplines including data mining, learning theory, data visualization, machine learning and psychometrics \[2\]. And this emerging field of EDM examines the unique ways of applying data mining techniques to solve educationally related issues.

This paper is to synthesize and share various examples by using data mining in education, to support reflection on teaching and learning. The background of EDM is described, then various algorithms that frequently used are briefly presented. Some specific EDM methods found are described. Subsequently, several examples of applications demonstrate how data mining are used to save resources and help teachers and learners. Finally, we conclude the paper.

# BACKGROUND OF DATA MINING

Data mining is an interdisciplinary subfield of computer science \[3-5\]. Data mining is the analysis step of the "knowledge discovery in databases" process, or KDD \[6\]. Data mining techniques have their roots in machine learning, artificial intelligence, computer science, and statistics etc. \[7\]. And data mining is an exploratory process, but it can be used for confirmatory investigations \[8\]. It is different from other searching and analysis techniques because data mining is highly exploratory, where other analyses are typically problem-driven and confirmatory. Through the combination of an explicit knowledge base, sophisticated analytical skills, and domain knowledge, hidden trends and patterns are able to be uncovered. These trends and patterns form the predictive models that enable to assist organizations with uncovering useful information then guide decision-making \[9\].

The Cross Industry Standard Process for Data Mining (CRISP-DM) is a cycle process for development and analysis of data mining models \[10\]. As the demand for data mining increases and more algorithms are created, CRISP-DM ensures practices that everyone can follow, and it gives specific tips and techniques on how to understand business data by deploying a data-mining model. CRISP-DM has six phases including business understanding, data understanding, data preparation, modeling, evaluation, and deployment \[10\].

# BACKGROUND OF EDUCATIONAL DATA MINING

Educational data mining as a field for solving educationally-related problems, at a high level, it seeks solutions to improve methods for exploring the data, which usually has meaningful hierarchy at multiple levels, in order to discover new insights of how people learn in the context of these settings \[11\]. For instance, a student's college transcript may contain a temporally ordered list of courses taken by him or her, the grade that the student earned in each course, and information about when the student selected or changed his or her academic major or minor. We might also understand how different individuals engage with or potentially 'game' the EDM system. Taken together, these learning analytics provide much useful information for the design of learning environments.

EDM applies lots of techniques such as Decision Trees, Neural Networks, Naïve Bayes, K-Nearest neighbor, etc. into its examples. Qualitative techniques such as interviews and document analysis are frequently used to support case studies in EDM. EDM impact students directly within areas of course content, course selections, recommender systems and admissions. In addition, applications of specific data mining methods like web mining, classification and multivariate statistics are key techniques used in educationally related data \[12\]. These approaches can be applied into modeling students' individual difference and respond to those differences by providing a way, which help improve students' performance \[13\].

# ALGORITHMS OF DATA MINING

Data mining relies on disciplines like classification, categorization, estimation, and visualization. Classification assists with identifying associations and clusters, and separates subjects under study. E.g., education institutions can use classification comprehensively to analyze students' characteristics. Categorization applies rule induction algorithms to handle categorical outcomes. Estimation includes predictive functions or likelihood deals with continuous outcome variables. Estimation and classification use unsupervised or supervised modeling techniques. Visualization uses interactive graphs to demonstrate mathematically induced data and scores, and is much more sophisticated than traditional bar charts or pie charts.

An algorithm is a specific, mathematically driven data mining function, such as a neural network, classification and Data mining relies on disciplines like classification, categorization, estimation, and visualization. Classification assists with identifying associations and clusters, and separates subjects under study. E.g., education institutions can use classification comprehensively to analyze students' characteristics. Categorization applies rule induction algorithms to handle categorical outcomes. Estimation includes predictive functions or likelihood deals with continuous outcome variables. Estimation and classification use unsupervised or supervised modeling techniques. Visualization uses interactive graphs to demonstrate mathematically induced data and scores, and is much more sophisticated than traditional bar charts or pie charts.

An algorithm is a specific, mathematically driven data mining function, such as a neural network, classification and regression tree (C& RT), or K-means. Data mining techniques including algorithms such as clustering, classification, regression, neural networks, association rules, decision trees, some of them have been applied successfully in the educational area \[14\]. E.g., methods for hierarchical data mining and longitudinal data modeling have been applied into EDM.

## Clustering

The goal of clustering is to find data points that naturally group together, splitting the full data set into clusters sets. By using clustering techniques we are able to further identify dense and sparse regions in object space, and discover overall distribution pattern and correlations among data attributes. Here are some researches about clustering techniques \[15-17, 18, 19\].

## Classification

Classification is used to predict values for some variables. This algorithm frequently employs the decision tree or neural network-based classification algorithms. In classification test, data are used to estimate the accuracy of the classification rules. If the accuracy is acceptable then rules are able to be applied into new data. Some popular classification methods include decision trees, logistic regression (for binary predictions) and support vector machines.

## Association rule

Association rules are used to find relations between different items \[15\]. Back to 1995, the analysis method of association rule was frequently utilized in most studies on educational data mining because of its less extensive expertise while comparing with other methods \[20, 21\]. However, after the year of 2005, as researchers frequently adopting clustering and classification methods, the trend changed.

## Regression

Regression analysis can be applied to model the relationship between independent variables and dependent variables. Independent variables are attributes we already known and response variables are able to predict what we want. But a number of real-world problems are not simply prediction. Therefore, more complex techniques (e.g., decision trees, or neural nets) may be necessary for future prediction. Some popular regression methods within educational data mining include linear regression, neural networks, and support vector machine regression.

## Neural Networks

Neural networks have the remarkable ability to obtain meaning from complicated or imprecise data and can be used to extract patterns and detect trends that are too complex to be noticed by humans or computer techniques, they are good at identifying patterns or trends for future forecasting needs.

## Decision Trees

Decision tree is tree-shaped structures, which represents decisions sets. Specific decision tree methods include classification and regression trees and Chi Square Automatic Interaction Detection.

## Nearest Neighbor Method

Nearest neighbor method, also called the k-nearest neighbor technique sometimes, classifies each record in the dataset based on a combination of different classes of the k record(s), which is similar to that in a historical dataset (where k \>=1).

To choose the appropriate algorithms, researchers need design the data and align it with the desired output. As for small-scale, they can opt for clustering approach since it does not require necessary splitting data in classification.

# METHODS OF EDUCATIONAL DATA MINING

EDM has a large number of powerful methods \[2, 22\], some of which are widely acknowledged to be universal data mining types (E.g., clustering, prediction, outlier detecting, relationship mining, etc.). However, Discovery with Models and Distillation of Data for Human Judgment are considered more prominent methods recently in educational data mining \[14\].

## Clustering

Clustering approaches had been applied to obtain clear distinction between the clusters. Once a set of clusters has been determined, new instances can be classified by determining the closest cluster. Clustering is able to be applied into grouping similar course materials or grouping students based on their learning behavior patterns \[24\]. There was an e-learning report was design to forecast the students' behavior pattern within the data mining approaches \[25\].

### Prediction

Prediction can deduce single aspect of the data from combinations of data in other aspects. Classification, regression and density estimation are main types of prediction methods. Prediction has already applied into predicting students' performance \[26\].

### Relationship Mining

Relationship mining can identify relationships among variables and encode them in rules for later application. Broadly, there are four different types of relationship mining: association rule mining, sequential pattern mining, correlation mining, and causal data mining, some of them has already utilized into identifying the relationship in patterns of students' behavior, difficulties or mistakes that learners usually encountered with \[27\].

### Process Mining

Process mining has been used to extract data from event logs in an information system to form a clear presentation in the overall activities. There are three different subfields of process mining: model discovery, model extension and conformance checking. Process mining is reported able to reflect students' behavior in the sequence of course, grade etc. \[28\].

## Text Mining

Text mining can be used for deriving information with high accuracy from text data and resources. The main contents of the text mining contain text categorization, text clustering, concept/entity extraction sentiment analysis and document summarization etc. Text mining is applied to analyze contents from forums, chats, Webpage or text resources etc. \[29\].

## Social Network Analysis

Social Network Analysis (SNA) is used to measure the relationships among different entities of information, and it is able to analyze the relationships in various tasks \[30\].

## Discovery with Models

In discovery with models, a model is validated via prediction, clustering, or manual knowledge engineering. Key applications of this method include discovering relationships among student behaviors, characteristics and contextual variables \[14, 31\].

## Distillation of Data for Human Judgment

Distilled data enables humans to identify well-known patterns for identification, it is also applicable for classifying data features. The goal of this method is to summarize and present the information in a useful, interactive and visually appealing way for understanding the large amounts of education data and supporting decision making \[32\]. Data is distilled for human judgment in educational data mining for two main purposes: identification and classification.

# APPLICATIONS IN EDUCATION

There are three applications of educational data mining with having received particular attention are discussed here.

## Predicting Student Performance

Lin applied classification and regression trees to predict what types of students would drop out from school, and then return to school later on \[33\]. The models were able to provide short-term accuracy for predicting which types of students would benefit from student retention programs \[33\]. Chacon and Spicer et al. developed a system based on data mining helps the institution identify and respond to students at-risk \[34\]. Their work is highly representative of the discipline, because it follows with a strict data mining process, which is quantitative. The research by Yeats, Reddy and Wheeler found that students who attend writing centers tend to do a good job in their classes \[35\]. Yu and DiGangi, et al. discovered that east coast students in USA tend to keep enrolled longer than their west coast counterparts do \[36\].

## Course Management System

EDM is often used in course management systems, like Moodle, which contains usage data that includes different activities. García, Romero, Ventura, and de Castro developed a simplified data mining toolkit that operates within the course management system and allows students and their learning users to get data mining information for their courses \[37\]. This research and application contributions will allow non-technical faculty to engage in educational data mining activities. Instead of traditional static course patterns, data mining can be applied to customize learning activities and adapt the pace for learners to complete courses \[38, 39\]. It will create significant and optimal learning experiences for each student. Also, Blikstein found different types of programming behaviors in an online course \[23\].

## Planning and Scheduling

Researches on mobile learning environments recently suggest that data mining can be applied to help provide personalized contents to different mobile users, despite the differences between mobile devices and conventional PCs. EDM applications will allow non-technical users engage in data mining tools and activities making processing more accessible for all EDM users \[39\]. There are some examples, including statistical and visualization tools, analyzing social networks and related influence on learning outcomes \[27\].

# CHALLENGES

As the related technology developed, costs and challenges associated with implementing EDM applications, like storing logged data and managing data systems \[19\]. Moreover, choosing which data to mine and analyze may also be a challenge. In addition, individual privacy is a continued concern for the application of educational data mining tools. With free, accessible tools in the market, students and learners may be at risk providing information to the learning system. Protecting individual privacy should be considered for the long-tern development of EDM. Moreover, it's unclear what data displays, visualizations, and visual analysis are most informative and support effective decision making for different stakeholders.

# CONCLUSION AND FUTURE INSIGHTS

Data mining is a powerful analytical tool to enhance decision making and analyzing new patterns and relationships for organizations. And EDM contains techniques including data mining, statistics, machine learning. DM need to analyze data coming from teaching and learning, tests learning theories, and policy decision-making etc. There are a number of opportunities exist in EDM, from an analysis at organizational level to the analysis at individual level. What's more, EDM is widely used and applied by learners, researchers and teachers, even institutions.

Recently, there are several studies focus on applying EDM into admissions and enrollment, but we don't know exactly how institutions using data mining to enhance student learning or improving related educational processes. And results from EDM research are typically obtained from the narrow context of specific educational settings. Therefore, the need for studies to examine in the broader context is necessary. For the overall EDM work to be completed, the urgent need of examining how to widespread the adoption of educational data mining is necessary. Furthermore, research indicates the area of educational data mining is concentrated in western cultures and subsequently, other countries like Asians may not be represented in the related researches and studies. Therefore, applications across multiple contexts should be considered in development of future models \[39\].

# Acknowledgment

This work was supported by International School of Software at Wuhan University. The authors want to thank teachers and staffs at International School of Software for the provided computational resources and support, as well as the anonymous reviewers for their insightful comments.


abstract: The Advanced Technology Large Aperture Space Telescope (ATLAST) is a set of mission concepts for the next generation UV-Optical-Near Infrared space telescope with an aperture size of 8 to 16 meters. ATLAST, using an internal coronagraph or an external occulter, can characterize the atmosphere and surface of an Earth-sized exoplanet in the Habitable Zone of long-lived stars at distances up to $\sim$`<!-- -->`{=html}45 pc, including its rotation rate, climate, and habitability. ATLAST will also allow us to glean information on the nature of the dominant surface features, changes in cloud cover and climate, and, potentially, seasonal variations in surface vegetation. ATLAST will be able to visit up to 200 stars in 5 years, at least three times each, depending on the technique used for starlight suppression and the telescope aperture. More frequent visits can be made for interesting systems.
title: Advanced Technology Large-Aperture Space Telescope (ATLAST): Characterizing Habitable Worlds

> M. Postman$^1$, W. Traub$^2$, J. Krist$^2$, K. Stapelfeldt$^2$, R. Brown$^1$,
> W. Oegerle$^3$, A. Lo$^4$, M. Clampin$^3$, R. Soummer$^1$, J. Wiseman$^3$, and M. Mountain$^1$
> *$^1$Space Telescope Science Institute, Baltimore, MD, USA*
> *$^2$Jet Propulsion Laboratory, California Inst. of Technology, Pasadena, CA, USA*
> *$^3$NASA/Goddard Space Flight Center, Greenbelt, MD, USA*
> *$^4$Northrop Grumman Aerospace Systems, Redondo Beach, CA, USA*

# Overview of the ATLAST Concept

We are at the brink of answering two paradigm-changing questions: Do Earth-sized planets exist in the Habitable Zones of their host stars? Do any of them harbor life? The tools for answering the first question already exist (e.g., Kepler, CoRoT); those that can address the second can be developed within the next 10-20 years (Kasting, Traub et al. 2009). ATLAST is our best option for an extrasolar life-finding facility and is consistent with the long-range strategy for space-based initiatives recommended by the AAAC Exoplanet Task Force . ATLAST is a NASA Astrophysics Strategic Mission Concept for the next generation flagship UVOIR space observatory (wavelength coverage: 110 nm – 2400 nm), designed to answer some of the most compelling astronomical questions, including "Is there life elsewhere in the Galaxy?" The ATLAST team investigated two different observatory architectures: a telescope with an 8-m monolithic primary mirror and two variations of a telescope with a large segmented primary mirror (9.2-m and 16.8-m). The two architectures span the range in viable technologies: e.g., monolithic vs. segmented apertures, Ares V launch vehicle vs. EELV, and passive vs. fully active wavefront control. This approach provides several pathways to realize the mission that would be narrowed to one as the needed technology development progresses and the availability of launch vehicles is clarified.

While ATLAST requires some technology development, all the observatory concepts take full advantage of heritage from previous NASA missions, as well as technologies available for missions in development. The 8 m monolith architecture is similar to the Hubble Space Telescope (HST), although the optical design is different. The 8 m mirror has an areal density exceeding that of the HST mirror, providing superb stiffness and thermal inertia. The use of such a massive mirror provides excellent wavefront quality and is made possible given a launch vehicle with the capabilities of the proposed Ares V. The 9.2 m and 16.8 m segmented mirror concepts rely heavily on design heritage from the James Webb Space Telescope (JWST), in development of lightweight segmented optics, the OTA deployment mechanics, and wavefront sensing and control. The 9.2 m segmented concept can be launched on a Delta IV Heavy launch vehicle with a modified 6.5 m fairing; the 16.8 m concept requires an Ares V. The non-cryogenic nature of ATLAST makes the construction and testing of the observatory much simpler than for JWST. We have also identified departures from existing NASA mission designs to capitalize on newer technologies, minimize complexity, and enable the required improvements in performance.

Four significant drivers dictate the need for a large space-based telescope if one wishes to conduct a successful search for biosignatures on exoplanets. First, and foremost, Earth-mass planets are faint – an Earth twin at 10 pc, seen at maximum elongation around a G-dwarf solar star, will have V $\sim$ 29.8 AB mag. Detecting a biosignature, such as the presence of molecular oxygen in the exoplanet's atmosphere, will require the ability to obtain direct low-resolution spectroscopy of such extremely faint sources. Second, the average projected angular radius of the Habitable Zone (HZ) around nearby F,G,K stars is less than 100 milli-arcseconds (mas). One thus needs an imaging system capable of angular resolutions of $\sim$`<!-- -->`{=html}10 to 25 mas to adequately sample the HZ and isolate the exoplanet point source in the presence of an exo-zodiacal background. Third, direct detection of an Earth-sized planet in the HZ requires high contrast imaging, typically requiring starlight suppression factors of 10$^{-9}$ to 10$^{-10}$. Several techniques are, in principle, capable of delivering such high contrast levels but all require levels of wavefront stability not possible with ground-based telescopes because the timescale of wavefront variations induced by the Earth's atmosphere is shorter than the time to measure the wavefront error to the required precision. A space-based platform is required to achieve the wavefront stability that is needed for such high contrast imaging. [^1] Lastly, biosignature-bearing planets may well be rare, requiring one to search tens or even several hundred stars to find even a handful with compelling signs of life. The number of stars for which one can obtain an exoplanet's spectrum at a given SNR and in less than a given exposure time scales approximately as D$^3$, where D is the telescope aperture diameter. This is demonstrated in Figure 1 where we have averaged over different simulations done using various starlight suppression options (internal coronagraphs of various kinds as well as an external occulter). To estimate the number of potentially habitable worlds detected, one must multiply the numbers in Figure 1 by the fraction of the stars that have an exoplanet with detectable biosignatures in their HZ ($\eta_{\oplus}$). The value of $\eta_{\oplus}$ is currently not constrained but it is not likely to be close to unity. One must conclude that to maximize the chance for a successful search for life in the solar neighborhood requires a space telescope with an aperture size of at least 8 meters.

All ATLAST concepts require many of the same key technologies. We believe these designs compose a robust set of options for achieving the next generation of UVOIR space observatory in the 2020 era. There are several fundamental features common to all our designs. All ATLAST concepts are designed to operate in a halo orbit at the Sun-Earth L2 point. The optical designs are diffraction limited at 500 nm (36 nm RMS WFE) and the optical telescope assembly (OTA) operates near room temperature (280 K – 290 K). All OTAs employ two, simultaneously usable foci: a three-mirror anastigmat channel for multiple, wide field of view instruments, and a Cassegrain channel for high-throughput UV instruments and instruments for imaging and spectroscopy of exoplanets (all designs have an RMS WFE of $<5$ nm at $<2$ arcsec radial offset from Cass optical axis). The ATLAST concept study is available on line and on the astro-ph archives .

# Simulated Data and Exoplanet Characterization Performance

Estimates of the SNR of habitability and biosignature features in an Earth-twin spectrum, achievable with ATLAST, are shown in Table 1. For these calculations we use a fully validated model of the Earth's spectrum , in combination with the observed visible reflection spectrum of the present Earth. We assume that the exoplanet is at maximum elongation and that the planet is observed for a length of time sufficient to achieve an SNR of 10, at a spectral resolution R = 70, in the red continuum. The vegetation signal is the enhanced albedo of the Earth from land plants, for wavelengths longer than  720 nm , with a modest SNR. Column 3 gives the width of the spectral feature. All of these SNR values can easily be improved with re-visits. In addition, ATLAST will allow us to glean substantial information about an exo-Earth from temporal variations in its features. Such variations inform us about the nature of the dominant surface features, changes in climate, changes in cloud cover, and, potentially, seasonal variations in surface vegetation . Constraints on variability require multiple visits to each target. The 8-m ATLAST (with internal coronagraph) will be able to observe $\sim$`<!-- -->`{=html}100 different star systems 3 times each in a 5-year interval and not exceed 20% of the total observing time available to the community. The 16-m version (with internal coronagraph) could visit up to $\sim$`<!-- -->`{=html}250 stars three times each in a five-year period. The 8-m or 16-m ATLAST (with a single external occulter) can observe $\sim$`<!-- -->`{=html}85 stars 3 times each in a 5-year period, limited by the transit times of the occulter. Use of two occulters would remove this limitation.

Habitability and Bio-Signature Characteristics

Feature

λ (nm)
Δλ (nm)
SNR
Significance

Reference Continuum

750
11
10

Rayleigh Scattering
<500
100
4
Protective atmosphere

Ozone (O3)
580
100
5
Source is O2; UV Shield

Oxygen (O2)
760
11
5
Plants emit; Animals inhale

Cloud/Surface reflection
750
100
30
Rotation indicator

Land Plant reflection
770
100
2
Vegetated land area

Water vapor (H2O)
940
60
16
Needed for life

Figure 2 shows two simulated ATLAST spectra for an Earth-twin at 10 pc, one at R=100 and one at R=500, taken with sufficient exposure to reach SNR=10 at 750 nm in the continuum. A 3-zodi background was used (local plus exosolar). The R=100 exposure times are 46 ksec and 8 ksec, respectively, for an 8-m and 16-m space telescope. The corresponding exposure times for the R=500 spectrum are 500 ksec and 56 ksec, respectively, for the 8-m and 16-m telescopes. The reflected flux from an Earth-like rocky planet increases as M$^{2/3}$, where M is the exoplanet mass. Hence, the exposure times for a 5M$_{\oplus}$ super-Earth would be $\sim$`<!-- -->`{=html}3 times shorter. At both resolutions, the O$_2$ features at 680 nm and 760 nm are detected, as are the H$_2$O features at 720, 820, 940, 1130, 1410, and 1880 nm. Rayleigh scattering is detected as an increase in reflectivity bluewards of 550 nm. The higher spectral resolutions enabled by large-aperture space telescopes enable the detection of molecular oxygen in exoplanets with lower abundances than those on Earth and provide constraints on the kinematics and thermal structure of the atmosphere that are not accessible at lower resolution.

For a 16-m class space telescope, time-resolved spectroscopy over intervals of a few hours may reveal surface composition variations, if the planet is not cloud dominated, as the exoplanet rotates. However, even broadband photometry can be used to detect short-term variations in albedo that can determine the rotation period and constrain the amount of cloud cover. Ford et al. (2003) generated model light curves for the Earth over 6 consecutive days using data from real satellite observations. Photometric variations of 20 – 30% on timescales of 6 hours were typical in the B,V,R,I passbands. Their models are shown in Figure 3. On top of these light curves we show grey bands whose height represents the $\pm$`<!-- -->`{=html}5% uncertainty for SNR=20 broadband (R=4) photometry and whose length represents the time it would take to perform such an observation with a 4-m, 8-m, and 16-m space telescope. We show the results of such calculations for an Earth twin at 10 pc and 20 pc. The 4-m has marginal capability to study such photometric variations at 10 pc but telescopes with apertures of 8-m or larger would be able to perform the measurements well as the integration time is less than the typical period between significant albedo changes. At 20 pc, even an 8-m telescope reaches its limits but a 16-m telescope is still able to acquire the needed accuracy in photometry in less than 4 hours.

Transit spectroscopy with ATLAST will permit characterization of super-Earth mass exoplanets. Figure 4 shows two simulated ATLAST transit spectra for planets around a 6th magnitude M2 star where the orbital period in the HZ is $\sim$`<!-- -->`{=html}20 – 30 days. Such observations are time consuming but do not require the use of a coronagraph or occulter.

In summary, with ATLAST, we will be able to determine if HZ exoplanets are indeed habitable, and if they show signs of life as evidenced by the presence of oxygen, water, and ozone. ATLAST also will provide useful information on the column abundance of the atmosphere, the presence of continents and oceans, the rotation period, and the degree of daily large-scale weather variations.

[^1]: Such stability is achievable regardless of the space telescope's aperture (within reasonable limits) – so this driver is primarily for the space locale of the telescope rather than its aperture.


abstract: The Golgi apparatus has intrigued researchers since its discovery and despite the advances, there are still many open questions in regards to its shape and function. Here, we propose a mechanical model of Golgi apparatus stack and explain its most elementary geometrical properties: the equilibrium number of cisternae, the stack size, and its general equilibrium shape. In order to gain comprehensive insight into the stack morphology, we combine both analytical and numerical methods. We successfully reconstruct the stack morphology using minimal assumptions, thus identifying the most dominant factors regulating the stack shape. We demonstrate that energy-wise the stack prefers an overall bent shape and show strong evidence that the adhesion strength determines the equilibrium number of cisternae per stack. We explore the morphological role of fenestrations and discuss their impact on the overall stack structure. We also comment on the effects of the membrane asymmetry on the shape of the cisternae and thus offer a broad steady-state study of the stack morphology and present a method that can be used also for other membrane-bound organelles.
author: U. Jelerčič
title: Mechanical model of Golgi apparatus

# Introduction

The Golgi apparatus is an essential cell organelle found in most eukaryotic cells. It is a highly ordered structure, typically composed of flattened stacks which are surrounded and interconnected by tubular networks and small spherical vesicles . Each stack consists of smooth, flattened membrane-bound units referred to as cisternae which are aligned roughly on top of one another. It is typically composed of 4 to 8 cisternae  although this number can vary between different species. The number of cisternae per stack is a very well-defined and robust quantity – it is generally conserved within a specific cell line  and is impervious to disassembly and subsequent assembly of the whole organelle .

The stack has three basic regions which differ in function, structure, and in microscopic membrane composition: the *cis* Golgi, the *mid* Golgi, and the *trans* Golgi region . The *cis* side of the organelle is normally convex and represents the part of the stack where matter enters the organelle. The *trans* side is on the other side of the stack and is usually concave, and it serves as the exit site for the modified matter. Both *cis* and *trans* cisternae exhibit buds and dynamic tubules with budding profiles which stretch from their margins . Furthermore, a single cisterna is structurally, chemically, and functionally unique . It is usually large (up to several micrometres in diameter ) and narrow (10 – 20 nm). The cisternae are stacked close to one another, with typical intercisternal distance of 7 – 15 nm  and although the cisternae do not touch physically, they maintain a very uniform distance.

All cisternae are perforated, albeit not homogeneously, with the holes referred to as fenestrations. These can be divided into three classes according to size : i) large well-like openings (\< 100 nm) are aligned over several consecutive cisternae and often contain budding profiles and free vesicles, ii) medium (65 – 100 nm) and iii) small holes (\< 65 nm) are unaligned and are separated from one another by 60 – 70 nm . The openings are abundant on the *cis* side, decrease in number towards the *mid* cisternae and then in general increase again on the *trans* side  (note that some older works report a different distribution of fenestrations ). Lateral distribution of fenestrations within a single cisterna is also specific with an average cisterna being characterised by a non-fenestrated central region and fenestrated outer parts . The shape of the fenestrations is less regular close to the rims of the cisternae and the fenestrations themselves can be, in the extreme case, big enough to create tubules extending from the stack .

Even though fenestrations are very abundant, their function is not clear. Ladinsky et al.  suggest that they act as spot welds to prevent cisternal swelling or they might help in the budding processes by creating regions of high local curvature . The former argument seems to agree with the experimental observations showing a significant difference in the thickness of perforated and unperforated cisternae. The cisternae with fenestrations were found to have a thickness of 50 – 100 nm, whereas the unfenestrated ones had a thickness of only $\approx$ 20 nm . Furthermore, it has been speculated early on that the openings may enable the cisternae to have a greater functional variability due to increased shape pliancy  but the exact connection is yet to be established.

Current theoretical work and *in vitro* experiments are greatly limited to simplified systems of isolated vesicles and vesicle stacks. They are not always directed towards the explanation of the actual organelle shape but rather touch on various phenomena otherwise observed in Golgi apparatus stacks. In the case of unfenestrated structures, most of the research focused on accurately reproducing various stacked erythrocyte morphologies  finding that the number of vesicles in the stack depends on the strength of adhesion between them . Higher genus vesicles have also been thoroughly investigated (both theoretically and experimentally), albeit only isolated – within the family of toroidal vesicles  and in the case of vesicles with 2 or 3 holes . Further theoretical generalisation to vesicles with an even higher genus is significantly more limited . Finally, we point out a recent study by Tachikawa *et al.*  who focused on the morphogenesis of the unfenestrated Golgi apparatus stack and presented an interesting dynamic pathway for the reassembly of the stack structure.

To our knowledge no research successfully addressed highly fenestrated structures in a stacked formation, let alone treated them numerically and with a great degree of generality. The subject of the equilibrium shape of Golgi apparatus stack thus remains open and represents a great challenge both for theoreticians and experimentalists.

In this paper we address fenestrated Golgi apparatus stacks and present a computational framework (analytical and numerical) based on the Helfrich theory of bending elasticity used to determine their equilibrium shape. Qualitatively, we answer three main questions; what is the equilibrium shape of the stack as predicted by bending energy minimisation; does this shape resemble the native morphology of the Golgi apparatus stack thus justifying the use of elasticity theory as a method of calculating its shape; and is it possible to find a simple mechanism which can determine the equilibrium number of cisternae per stack and their thickness? All in all, we successfully reproduce native-like stack shapes and identify adhesion strength and membrane composition as possible factors determining the stack geometry. We furthermore discuss the role of the fenestrations, membrane asymmetry, and cell environment on the stack shape, thus presenting an exhaustive study of the main aspects of the Golgi apparatus morphology.

# Model

In order to calculate the free energy of Golgi-like vesicle stack the bending energy of each cisterna has to be determined. In our model, each cisterna is approximated as a flattened pancake-like vesicle with volume $V_i$ and surface area $A_i$. The vesicles are then stacked one on top of another and are assumed to adhere *via* contact adhesion mechanism. If there are $n$ cisternae in the stack, the total volume and surface area of the stack are then $V=\sum_i^n V_i$ and $A=\sum_i^n A_i$, respectively and the reduced volume of the stack can be defined as $v=6\sqrt{\pi}V/A^{3/2}$.

Since the Golgi stack is observed to preserve its size over a long period of time, the total volume and surface area are assumed to be fixed. However, volumes and surface areas of individual cisternae can vary and thus reflect the dynamic nature of the organelle (*i.e.* vesicular transport through the stack and the tubular sections connecting the cistarnae). The free energy of the stack is in our model governed by the bending energy of the cisternal membrane and the adhesion energy due to the association of the cisternae (thermal energy is omitted because we assume that it is suppressed by adhesion). The bending energy of a single cisterna itself is determined using the basic Helfrich bending energy term $$F_i=2\kappa\int H^2\mbox{d}A+4\pi\kappa_G(1-g), \quad \mbox{for } i=0,...,n
\label{eq:Fi}$$ where $H$, $\kappa$, $\kappa_G$ and $g$ denote the mean curvature, mean bending modulus, Gaussian bending modulus and the genus of the vesicle, respectively. The second term in Eq. (), *i.e.* the Gaussian term, is a topological invariant and is explicitly taken into account when the cisternae feature highly non-trivial topology (*i.e.* many fenestrations) – it can be omitted when considering structures of the same topological order but should be included in all other general cases, especially if the genus is very high. Finally, the total free energy of the stack is the sum of the bending energies of individual cisternae and the adhesion energy $$F=\sum_i^n F_i-\Gamma A_c,
\label{eq:F}$$ where $\Gamma$ and $A_c$ represent the adhesion strength and the total surface area of the contact zones between vesicles, respectively. The adhesion strength can be alternatively expressed as a reduced quantity $\gamma=\Gamma R_s^2/2\kappa$, where $R_s$ refers to the radius of the sphere with the surface area $A$ and the free energy is measured in the units of $8\pi\kappa$.

We do not include the area difference elasticity term because we assume that the lipid membrane continuously relaxes due to the steady transport along the organelle. The membrane area difference is in our model therefore not fixed albeit non-vanishing – it amounts to the value at which the total free energy is minimal and differs between stacks with different equilibrium shapes. Furthermore, spontaneous curvature term is also omitted thus neglecting specific composition of the membrane (*e.g.* highly conical lipid molecules) and the membrane is assumed to be homogeneous.

In each particular case, the equilibrium configuration is determined by analytically and numerically minimising the free energy \[Eq. ()\] while keeping the total volume $V$, total surface area $A$ and adhesion strength $\Gamma$ fixed. The number of cisternae in the stack $n$ and the topological order of cisternae are also fixed during an individual minimisation.

# Equilibrium shapes of Golgi apparatus stack

Golgi apparatus stack is topologically extremely complex (involving hundreds of differently sized fenestrations) and thus mathematically very difficult to treat. In general, the analytical models give an excellent insight into the dominant characteristics of the stack, but are limited by the assumptions necessary for the algebra to be manageable. On contrary, numerical solutions offer an exact representation of the actual morphology of the stack within the model, but may be computationally very demanding. A joint approach, used in this paper, overcomes limitations of each individual method and provides us with a comprehensive view of the stack morphology.

***Analytical results.*** In its most basic form, the Golgi stack can be represented as a system of $n$ identical flat discs as first proposed by Derganc *et al.*  (Fig. a; termed DMS model). The stack has a fixed total volume and surface area, while the discs themselves can exchange both membrane and lumen. The free energy of such stack combines the bending energy associated with the rims and the adhesion energy (see Supplementary Information). Equilibrium shapes of the stack are determined by minimising the free energy at fixed $V$ and $A$ while varying the adhesion strength $\Gamma$ and the number of cisternae $n$. Derganc *et al.* found that at fixed reduced volume $v$ the equilibrium number of cisternae (*i.e.* the size of the stack) depends on the reduced adhesion strength. We generalise this result by considering stacks with alternative geometries – truncated pyramid/truncated double pyramid (where the radial size of each subsequent cisterna differs by the same amount; Fig. b) and curved (with cisternae in the shape of spherical caps; Fig. c).

When comparing the energies of different stack variants shown in Fig.  (see Supplementary Information), we find that in the regime of relevant reduced adhesion strengths (*i.e.* large enough to guarantee that the cisternae form a stack and not a cluster of isolated vesicles) and at small reduced volumes, the DMS model always energy-wise outperforms both pyramidal variants while the curved stack alternative proves to be superior to all three. The actual equilibrium degree of stack curvature depends on several parameters ($n$, $v$, and $\gamma$), but it is always the maximal possible in a given case. Since we artificially limit the maximal curvature, the shapes at biological reduced volumes and within our model conform to the hemispherical shape (Fig. ). We can generalise this result by hypothesising that without the limitation, the stack prefers to be fully curved, *i.e.* resembling nested stomatocytes with only a small opening. We therefore analytically show not only that the curved morphology is optimal but also that virtually any equilibrium number of cisternae in the stack is realisable by varying the reduced adhesion strength (Fig. ). Unless too large ($v\gg0.1$), our results qualitatively do not depend on reduced volume.

***Numerical results.*** The stack morphology is evaluated numerically using C-based Surface Evolver program package . Within its framework, a surface is represented as a simplicial complex and the equilibrium shape is determined by a gradual minimisation of the predefined free energy – in our case free energy defined by Eq. (). As before, we additionally constrain the total volume and total surface area and assume that, in equilibrium, the membrane is relaxed. We search for the deepest local energy minima within the space of physically reasonable shapes set by initial conditions and conjecture the equilibrium morphology by analysing the stack shape and total energy of thus found structures.

As a proof of principle, and in connection with the analytical results, we first focus on simple systems with a small number of cisternae and start with non-fenestrated stacks. We confirm the analytical prediction that the flat stack does not represent a stable shape and demonstrate that the free energy decreases as the overall bending of the stack increases (Fig. a, b). Apart from being energetically optimal due to the increased contact area, the bent stack shape also features a stable cisternal thicknesses as observed in nature. We note that we find a third group of shapes (besides flat and bent) that has even smaller free energy and resembles a stack of nested stomatocytes where the inner-most cisterna adapts a spherical shape (cup-like structure, Fig. c). We do not include this solution in further analysis because it proves irrelevant (*i.e.* geometrically impossible) in the general case of fenestrated stacks and cannot exist at any given small reduced volume.

We extend the analysis to stacks with larger number of cisternae and compute their equilibrium shape at different values of $\gamma$. We find that, in agreement with the analytical solutions, the reduced adhesion strength represents a robust and simple parameter able to determine different equilibrium numbers of cisternae $n_{eq}$ (Fig. ). Even though we illustrate this behaviour on a simple system and with agreement with the analytical results, we can assume that the trend numerically observed for small $n_{eq}$ continues into the regime of more realistic stack parameters ($n_{eq}\simeq 5$). We thus postulate that a stack with any number of cisternae can be theoretically stabilised by a large enough reduced adhesion strength, although not any value of $\gamma$ may be biologically relevant.

Treating systems of interconnected vesicles with fenestrations is very demanding, albeit unavoidable when discussing more realistic Golgi apparatus models. Apart from purely technical issues connected to the surface triangulation there is also a conceptual problem regarding the accurate description of Gaussian bending energy term. Contrary to the mean bending modulus $\kappa$, the Gaussian bending modulus $\kappa_G$ cannot be easily measured and has been evaluated only for selected systems – in the case of simple lipid membranes, the majority of experimental and theoretical works find $\kappa_G$ to be negative and close to the value of $\kappa$ . On contrary, Golgi apparatus membrane is highly complex and experiments have shown that the regions with high curvature (*e.g.* rims and fenestrations) contain specialised proteins  thought to act as shape stabilisers . The composition of such membrane can thus in principle effect $\kappa_G$, making it either smaller in absolute value or even positive . The sign of $\kappa_G$ in turn strongly effects the energy landscape in heavily fenestrated stacks – a large negative $\kappa_G$ discourages formation of fenestrations, while positive $\kappa_G$ promotes them. Due to the lack of quantitative estimates, we perform the calculations assuming $\kappa_G=\pm\kappa$ thus illustrating the general impact of fenestrations on stack morphology. We note however, that the value of $\kappa_G$ has no influence on the final shape within an individual minimisation (where $A$, $V$, $\gamma$, and $n$ are fixed). We therefore numerically determine the equilibrium shape by setting $\kappa_G=0$ and consider the Gaussian contribution by adding $4\pi\kappa_G(1-g)$ only when we are interested in the actual value of the free energy (*e.g.* when we determine equilibrium number of cistarnae per stack).

As already mentioned before, the cup-like solution found in the case of trivial topology does not survive the generalisation to a given high genus, because it is physically impossible to have a bulk spherical shape riddled with fenestrations. In the case of flat stacks, we see that the stacks with fenestrations maintain a stable lumen even at very small, biologically relevant reduced volumes, where non-fenestrated stacks collapse (Fig. ). This is an important structural feature of the stack consistent with the experimental observations and it is a direct consequence of the introduction of fenestrations. Bent stacks, on the other hand, do not experience such drastic changes in cisternal morphology as flat ones when the genus is increased. At small reduced volumes, the lumen of the cisternae is stable even without fenestrations, although their presence in general guarantees a more uniform thickness along the stack.

Within our approach, the fenestrations are free to move, although they do not move significantly[^1]. Initially, we distribute them randomly and in the same fashion in all the cisternae. Based on our solutions it seems that even though the fenestrations adjust their initial positions, they do not tend to cluster, as already established by Michalet *et al.*  in the case of necks in an infinite membrane sheet. We furthermore find that in cases where the fenestrations are scarce, their radii increase towards the rims of the cisternae which agrees with experimental observations remarkably well. This effect is less prominent in strongly fenestrated stacks where the fenestrations express more uniform radii. The radii themselves are found to depend mostly on the relative ratio between their total number and the reduced volume. As expected energy-wise, the fenestrations are roughly catenoidal in shape and thus contribute very little to the total bending energy of the stack, but due to their specific geometry (neck radius *vs.* height ratio) efficiently determine the cisterna thickness. Since the details of fenestrations' sizes and their distribution depend on many parameters, we refrained from further analysis in this direction. It, however, remains a most interesting topic for future work.

As we consider overall bent stacks, we find a similar trend as in the topologically trivial case – the free energy of the stack decreases as the bending increases while maintaining a remarkably uniform cisternal thickness at any given overall bending (Fig. ). The lower energy of curved shapes is also observed in stacks with an even larger number of cisternae, although the numerical analysis is then much more strenuous. It is not entirely clear, however, whether the fenestrated stacks prefer a completely closed stomatocyte morphology seen in stacks without fenestrations, or whether a curved shape with a finite opening is favoured instead. We find that in general, the free energy of a fenestrated stack decreases with the overall curvature of the cisternae at all numbers of cisternae and all numbers of fenestrations per cisterna that we have studied. Based on these results, it is reasonable to expect that the minimal-energy stack is characterised by a small opening (as seen in Fig. ), but this limit is technically difficult to reach.

Finally, we test whether the fenestrated stacks exhibit the same relationship between the equilibrium number of cisternae and reduced adhesion strength that we report in non-fenestrated stacks. Since the curved stomatocyte solutions cannot be easily computed, we perform the analysis using the flat stack solution (Fig. a) and assume that the results apply also to the general case. We find again that we can qualitatively regulate the equilibrium number of cisternae per stack simply by varying $\gamma$ – however, the actual value of $n_{eq}$ now depends on the choice of Gaussian bending modulus. For example, if we compare only the bending parts of the free energy ($\kappa_G=0$) we find that in the case of $v=0.09$ and 30 fenestrations per cisterna, the optimal number of cisternae in the stack changes from 3 (at $\gamma=100$) to 4 (at $\gamma=1000$). Taking into account the Gaussian term causes the global energy minimum to shift – if $\kappa_G<0$ there is an energy penalty per each fenestration and the stack with the smaller number of cisternae (and consequently fenestrations) is preferred. On the other hand, if $\kappa_G>0$ stacks with many cisternae and many fenestrations are stable and as such agree with the native shapes of the organelle (Fig. b).

In order to make our model more realistic, we briefly discuss the effect of membrane asymmetry on the shape of the stack. Even though we still consider that the majority of the stack relaxes the asymmetry *via* inter-organelle transport, we assume that outmost *cis* cisterna[^2] maintains a fixed preferred membrane area difference $\Delta A_0$. This simulates the matter entering the organelle *via* small spherical vesicles carrying large amounts of $\Delta A_0$ which is absorbed into the stack membrane after the vesicles fuse with the cisterna. A portion of the initial membrane area difference relaxes through membrane rearrangement (*e.g.* flip-flop mechanism) and the rest provokes a shape modification of the cisterna. In this case we use a modified bilayer-couple model and find that even a small increase of membrane area difference in comparison to the relaxed value in the rest of the stack noticeably changes the appearance of the affected cisterna (Fig. a). The radius of the cisterna increases and, even more dramatically, the radii of the fenestrations greatly increase in order to maximise the total surface area of curved regions capable of accommodating incoming $\Delta A_0$ (Fig. b). We expect a similar behaviour also on the *trans* side, although this is more difficult to explain – contrary to the *cis* side where $\Delta A_0$ increases passively via inter-organelle transport, on the *trans* side the membrane asymmetry has to be first actively generated in order for the vesiculation and secretion to occur.

We hypothesise that en even greater increase of $\Delta A_0$ results in such an enlargement of fenestrations that the cisterna transitions from the flattened perforated shape towards the form of a reticular network. This is in agreement with the observations of the native shape of Golgi apparatus, where outmost cisternae feature large openings and often continue into reticular/tubular networks.

# Discussion

***Overall shape of the stack.*** Using simple models, we show that a bent stack (global solution) is in general always energetically preferred over a flat one (local solution) which is consistent with observations and reflects the native morphology of the majority of the mammalian Golgi apparatus stacks. As expected, based on the model assumptions, the exact degree of curvature is not reproduced accurately in our solutions and there are several conceptual improvements that can be implemented easily, albeit are very time-consuming (*e.g.* dissimilar $\kappa$ along the stack mimicking chemical uniqueness of cisternae). These can provide us with quantitatively more exact results, but do not introduce qualitatively new findings.

Nevertheless, we stress that the whole energy landscape, as presented in this paper, is not necessarily reachable in every biological system. This is consistent with our minimal energy approach and merely takes into account that the solution found in nature is not a global minimum, but a minimum within a set of additional constraints. Together with factors determining the specific composition of membrane, these constraints can be also determined by the biochemical environment (*e.g.* $\Delta A_0$, pressure, surface tension, confinement). Surface tension or pressure are especially intriguing because they can be both manipulated *via* different cellular processes (*e.g.* cytoskeleton stresses or pH) and if we compare their values in stacks with different degrees of bending, significant differences are observed – compared to the flat solution, surface tension and pressure visibly decrease as the overall curvature of the cisternae increases (Fig. ). If, for example, there is a specific range of biological values for surface tension, the equilibrium shape has to be chosen as a minimum energy solution within this range and, depending on the values, less bent stack can be preferred. Given that there are differences in the composition of the Golgi apparatus membrane and the cytosol among species, this argument also offers a possible mechanism for the observed biological variability of stack shapes.

Apart from general insights into stack morphology, this work presents important new understandings of the role of fenestrations in Golgi apparatus. We see that they are essential for the lumen stability in flat and mildly bend stacks and affect the cisterna morphology in two distinct ways: i) in the central region the fenestrations act as pillars and keep the structure from collapsing, while ii) the fenestrations that are closer to the rims seem to pull the outmost membrane leaflets slightly together and thus have the role of spot welds as hypothesised by Ladinsky *et al.* . Both mechanisms complement each other in order to create a large region of uniform thickness where the biochemical reactions can take place thus supporting the organelle function. The increase in the cisterna thickness as seen in our analysis is in agreement with the experimental results, where fenestrated cisternae were reported to be up to 5 times thicker than the unfenestrated ones . Moreover, in a different experiment, temperature blocked cells with an eliminated outward transport were found to have a strongly reduced function and also a significantly decreased number of fenestrations  – combining these two observations, we can speculate that the fenestrations ensure well regulated thickness of cisternae, which in turn positively affects the function of the organelle. The size of the fenestrations themselves is found to be transport-mediated and depends on the details of the membrane rearrangement due to the increased membrane area difference in the outer cisternae. Varying only $\Delta A_0$, we can thus control the extent of reticulation, which is consistent with the dynamic nature of the Golgi apparatus peripheral regions.

***Equilibrium number of cisternae.*** The relationship between the adhesion strength and equilibrium number of cisternae is found to be universal and is valid regardless of the number of fenestrations or the degree of curvature in the stack. Even though the adhesion strength represents a phenomenological extensive quantity and we cannot determine its value directly from experiments, it is reasonable to expect that it is closely related to the composition of the membrane and of the matrix surrounding the cisternae. Similar is true also for the Gaussian bending modulus of the membrane, which in our model implicitly determines the equilibrium number of cisternae by maximising or minimising the total number of fenestrations. We can assume that the specific protein and lipid configuration of the membrane remains roughly the same within the same cell line but can significantly vary for cells with different roles in the organism (and even more between different species) – by depending on membrane composition, both parameters allow for great robustness and significant variability.

***Bending elasticity approach.*** Based on our results we conclude that we can use the bending energy minimisation approach to successfully determine the equilibrium shape of a Golgi apparatus stack even by employing the most basic Helfrich model. We show that many of the stack properties derive from a mechanical equilibrium and exhibit a great level of robustness and generality. We thus believe that any other purely function-related processes (*e.g.* ) that effect the stack morphology serve only to stabilise the shape and consequently have a supporting role in the Golgi apparatus shape determination.

Presented theoretical model introduces a new step in the understanding of the Golgi apparatus shape and shows that significant insights can be gained by using conceptually simple framework and only a few limiting assumptions. Our model thus represents a solid basis for further generalisation and development, especially by taking into account the details of membrane composition (*e.g.* variation of $\kappa$ and $\Delta A_0$ along the stack, cytoskeleton, localisation of proteins). The method can be further extended beyond Golgi apparatus morphology and can be applied to other membrane-bound cell organelles, giving a deeper insight into the biomechanics of the cell interior.

**Acknowledgments.** The author thanks P. Ziherl for his significant contribution to this work and N. Gov, A. Bernheim, A. Lošdorfer-Božič, G. Orly, and S. Čopar for critically reading the manuscript. This work was supported by the Slovenian Research Agency through Grant No. P1-0055.

[^1]: This is at least partially contributed to the complexity of the system, where a move by one fenestration effects all the rest through coupling with adhesion surface.

[^2]: We can easily generalise the approach to include other combinations of fixed $\Delta A_0$ or even a gradient along the stack.