abstract: Agent-based models have been employed to describe numerous processes in immunology. Simulations based on these types of models have been used to enhance our understanding of immunology and disease pathology. We review various agent-based models relevant to host-pathogen systems and discuss their contributions to our understanding of biological processes. We then point out some limitations and challenges of agent-based models and encourage efforts towards reproducibility and model validation.
author: Amy L. Bauer[^1]; Catherine A.A. Beauchemin; Alan S. Perelson[^2]
bibliography: allbibliographies.bib; abm-review.bib
date: 2024-10-02
title: Agent-Based Modeling of Host-Pathogen Systems: The Successes and Challenges

# Introduction

The ideas gleaned from studying immunology and host-pathogen systems may be relevant not only to human health but also to a wide array of other systems. A pathogen is any infectious agent that can lead to illness or disease of a host. Examples of pathogens include human immunodeficiency virus (HIV), *Mycobacterium tuberculosis*, the etiological agent for tuberculosis, the SARS coronavirus, and the influenza virus. In these cases, the host is usually a human being or an animal. However, in general terms, a host could just as easily be a computer network and the pathogen a computer virus. The immune system has memory and learns about the pathogens it encounters. It also must discriminate between self and non-self. Consequently, applications of the ideas intuited from immune system dynamics can be translated into algorithms relevant to learning, pattern discrimination, artificial intelligence adaptive behavior, and applied towards goals such as the development of new computer virus security applications . The interested readers should see .

Because of the difficulty in reasoning about large numbers of interacting components with non-linear interactions, mathematical modeling and simulation are becoming important research tools. Depending on what aspect of the host-pathogen system is being investigated, different mathematical modeling tools are employed. Below we discuss several modeling techniques commonly used to describe such systems.

Ordinary differential equation (ODE) models are often used as a starting point to describe host-pathogen systems. One of the advantages of using ODEs is that a lot is known about their behavior. ODE models are simple and elegant and require fewer parameters than their spatial counterparts (e.g., agent-based models or partial differential equations). This is an important consideration when experimental data is obtained from a well-mixed compartment, such as blood, or from a homogenate of a tissue, such as the spleen or lymph node, where spatial information is lacking. For example, in a recent model of influenza infection in humans, Baccam et al. used an ODE model not because the body is well-mixed but due to the fact that the data being analyzed was the amount of virus in nasal wash. When sufficiently simple, ODEs have the added benefit of being analytically tractable, and techniques such as bifurcation theory can be used to predict, for example, parameter values that switch a system from pathogen clearance to a chronically infected state. There are, however, limitations to using ODE models. These models assume that the populations (dependent variables) are homogeneous and uniformly distributed over the simulation space for all time. This is an assumption that may not be realistic, and that may significantly influence the resulting dynamics. To examine the effects of the assumption of spatial homogeneity, there is a growing body of research probing the effects of spatial distribution on systems in ecology , epidemiology , and immunology . One option to address this limitation is to use partial differential equations (PDEs), which capture changes in both time and space, but, in general, as these equations get more complicated, and consequently more computationally challenging, the advantages to using PDE models wane. For both ODE and PDE models, one must also consider that solutions to these equations only provide an average or mean field description of the system behavior with little or no information about the possible deviations from this aggregated behavior.

An alternative to differential equation models are agent-based models (ABMs). ABMs are stochastic models used to describe populations of interacting agents, such as insects and people, using simple rules that dictate their behaviors. These models were originally introduced by John von Neumann and Stanislaw Ulam under the name of "cellular spaces" as a possible idealization of biological systems. They sought to show that biological processes such as the reproduction and evolution of organized forms could be modeled by simple cells following local rules . A well known feature of ABMs is their ability to generate surprisingly complex and emergent behavior from very simple rules, including periodic behaviors or intricate spatial and temporal patterns . Agent actions are asynchronous, that is, they do not evolve at constant time steps. Instead, agents respond dynamically and independently to changing environmental or discrete event cues. Consequently, nonlinearities and time-delays are not difficult to treat empirically since they can be incorporated into the agent's rules or they may even emerge naturally as a consequence of the system's collective dynamics. Another advantage to ABMs is that their computational structure is inherently parallel and therefore can be implemented on parallel computers very efficiently.

In this article, we review a variety of agent-based modeling approaches and their contributions to our understanding of host-pathogen interactions and disease dynamics.

# Applications of Agent-Based Models

Agent-based models are quickly gaining in popularity. As experimental assays are developed that increase our understanding of host-pathogen interactions, the level of description desired for realistic and relevant models is increasingly more complex. Computer processing is becoming faster and more efficient, expanding the computational ability of computers and making possible the use of ABMs for complex systems. Moreover, because ABM implementation is achieved at the agent level, the description of the agents and rules tends to mimic the language used to describe the real system, that is, the description is more physical in character than mathematical. For example, models using differential equations consider rates of creation, rates of death, rates of binding, or rates of diffusion for whole populations of agents, whereas an ABM considers the rules guiding the actions of the agent. The familiar and natural modeling language used in agent based models not only makes ABMs approachable and useful to experimentalists and clinicians, but also engages them in the modeling process.

## ABMs as Immune System and Disease Simulators

A number of immune simulators have been constructed that provide a programming framework that in principle could incorporate all current knowledge of immunology and could be used to model any aspect of immune dynamics. Platforms like IMMSIM , SIMMUNE , Reactive Animation , and SIS are true immune simulators and allow users to define the rules of interactions and simulate an immune reaction. Most of the immune simulators are developed to make the interaction rules simple to define and easy to change in order to facilitate the exploration and impact of different rules on the development and outcome of an immune response. Some of these simulators are more flexible than others. With IMMSIM, for example, one would typically tune interactions by changing parameter values and one could, for example, completely turn off the humoral response by setting to zero the right set of parameters . With a system like Reactive Animation, users can potentially go even further and choose between what the authors refer to as "running theories" . When the rules of interactions between, for example, an epithelial cell and a T cell are not known there typically exist various hypotheses for the way in which such interactions proceed. The Reactive Animation system would allow the user to choose between different hypotheses and observe the impact of that choice on the simulation dynamics.

The usefulness and applicability of these simulators vary, but some have been applied to important immunological problems and their findings published in experimental journals. For example, IMMSIM was used to model affinity maturation and hypermutation in the humoral immune system , to test approaches to vaccine design , and to investigate mechanisms for tolerance to pathologic rheumatoid factors . IMMSIM has also been used as a pedagogical tool: Prof. Martin Weigert used IMMSIM in his seminar "Why Immune Systems Fail: Autoimmunity, Influenza Pandemics and HIV" at Princeton to demonstrate rather than just describe immune cell interactions . SIMMUNE became the cornerstone of the Program in Systems Immunology and Infectious Disease Modeling (PSIIM) at the National Institute of Allergy and Infectious Diseases , a program dedicated to the use of computational approaches to problems in immunology. Recently, immunologist Ronald N. Germain, in collaboration with SIMMUNE's creator, Martin Meier-Schellersheim, and others used the simulator to investigate the mechanism of chemo-sensing .

Another type of simulator, which we refer to as a *disease simulator*, is a general programing framework that can be tuned to model a specific disease including tumor growth, tuberculosis, or influenza. By changing parameters, such as the rate of spread of the infection, the lifespan of infected cells, or the binding rates of cytokines, users can calibrate these frameworks to model a variety of diseases. Three such simulators are CyCells , PathSim and the MASyV modules ma_immune and ma_virions . These simulators reproduce a variety of host-pathogen interactions, and are typically easier to use and calibrate than the immune simulators.

In contrast to general disease simulators, there is another class of ABMs that concentrates instead on a single disease or on specific aspects of immune mechanisms. Among these are models of: HIV , *M. tuberculosis* , Epstein-Barr , influenza , cancer vaccination , tumor growth and invasion , tumor-induced angiogenesis , and acute inflammation . Some of these simulators have greatly contributed to improving our understanding of the dynamics of a disease and represent novel and original approaches to modeling.

## The Use of Agent-Based Models by Experimentalists

The appeal of ABMs is such that their use goes beyond the realm of traditional modelers. While the involvement of experimentalists in mathematical modelling remains somewhat limited, ABMs in immunology/disease modeling have engaged experimentalists since their beginning. Indeed, one of the first ABM of immune interactions, IMMSIM, was the result of a close collaboration between Philip E. Seiden, a computer scientist, and Franco Celada, a well-respected experimentalist. This is a testament to the appeal of ABMs, where the description and representation of disease systems are close to that of the true biological system. This attraction for experimentalists is helping to bridge the gap between theoretical models and experiments.

Marc Jenkins of the University of Minnesota is one such experimentalist. New imaging methods have, for the first time, permitted the visualization of individual T cells and their interactions with antigen-presenting cells *in vivo* . Jenkins and his team assembled the emerging literature into a coherent picture of the first 50 hours of a primary immune response to an antigen in the form of an animation . Their model, or rather their animation, was implemented as a movie simulation using Macromedia Flash MX. It consists of a two-dimensional plane, which corresponds to a 10 slice (approximately one cell diameter) through a hypothetical spherical lymph node 2 in diameter, and includes T cells, B cells, and dendritic cells. Random motion paths for each cell in the simulation are pre-scripted individually and scaled to their known approximate speeds taken from the two-photon microscopy literature . Restrictions are put on the paths, such that B cells are confined to the follicles (except 5 hours after exposure to antigen when their movement is restricted to the outer edge of the follicles near the T cell area), dendritic cells to the T cell area, and T cells to the T cell area and outer edges of the follicles for 90% and 10% of their paths, respectively. When a collision occurs between two cells, the cells' motion along their respective path is halted for a defined amount of time to mimic real-time interactions. The time alloted to each interaction depends on the circumstances. The restrictions on the cell paths and the scripted interaction times were all taken from published experimental results.

Jenkins' approach is novel; rather than writing what they feel would be a good model at an appropriate level of detail and then seeking empirical data to calibrate the model, their model was built piece by piece from the experimental literature. The resulting simulation constitutes a very complete and convincing visual summary of the available knowledge at the time and is exemplary of the capability of ABMs to translate and synthesize a large body of compartmentalized research on a complex biological system. Moreover, movies of Jenkins' simulations make other important biological contributions. For example, these movies showed that it is possible to have B-T cell interactions in the absence of directed cell chemotaxis. Indeed, random mechanisms based on restriction of B cells to a specific spatial area rarely patrolled by T cells can also result in T-B cell interactions if T cells proliferate first, thus increasing the probability of a T-B cell encounter .

Another experimentalist seduced by the appeal of ABMs is Gary An, an M.D. in Trauma and Critical Care at Northwestern University Feinberg School of Medicine. Dr. An employed agent-based modeling to describe systemic inflammatory response syndrome and multiple organ failure , the first application of agent-based modeling to critical care. An important aspect of this ongoing work is the utilization of agent-based modeling to simulate clinical trials, providing a crucial link between theoretical modeling and the clinical setting . This model assumes that the primary immune response is initiated by endothelial cell injury at the cell-blood boundary and that the inflammatory response is translated from a local response to a systemic one. Model agents are endothelial cells, neutrophils, monocytes, several populations of T helper cells, and neutrophil/monocyte/T lymphocyte progenitor cells. The model also considers the role of a number of cell receptors, including the ICAM, IL-1, and TNF-$\alpha$ receptors, and relevant chemical mediators. Compared with results from other clinical trials and animal studies, this model reproduced the general behavior of the innate inflammatory response as measured by patient outcome and cause of death. *In silico* trials of anti-cytokine therapy were simulated and compared qualitatively well with those in published phase III anti-cytokine clinical trials. Dr. An then formulated a series of hypothetical treatment strategies and used his model to evaluate the efficacy of the proposed interventions. Even this relatively simple ABM of the immune response displays complex and counter-intuitive system dynamics. Conceptually, the use of ABMs to simulate clinical trials holds a great deal of promise and sets the stage for other theoreticians and experimentalists to use models to formulate and test treatment strategies before commencing actual animal or human trials.

## Studying Localized Spatial Effects

ABMs are stochastic models; as such they can reveal unique dynamics resulting from very specific spatial configurations or from rare localized events that would be missed using a mean field approach. This is a particularly valuable feature of these types of models, especially since similar features are responsible for variations seen in the development and outcome of infection in different individuals or in a given individual at different times. Many papers have investigated the effects of varying spatial configurations or the occurrence of rare localized events on infections .

Louzoun et al.  investigated the effects of spatial heterogeneity of antigen concentration on lymphocyte proliferation. First, the authors define a simple system with two species: antigen and lymphocytes. The antigen population grows at a constant rate and decays at a rate proportional to its concentration. Lymphocytes grow at a rate proportional to both their concentration and to the concentration of antigen, and die at a rate proportional to their concentration. For this model, the authors showed that if the concentration of antigen was above some threshold, the lymphocytes would proliferate, otherwise they would not. The authors then discretized the antigen and lymphocyte populations in both space and time, and incorporated the discrete addition or removal of agents (lymphocytes or antigen) with a given probability rather than at a given rate. This change resulted in the emergence of local *hot spots* of lymphocyte proliferation. So while the average antigen concentration was not above the threshold required for lymphocyte proliferation, regions where the local antigen concentration was above the threshold quickly dominated the system average with a growth rate proportional to the peaks of local concentration of antigen, and not to the average concentration throughout the system .

Strain et al. examined the specific case of an HIV infection and compared the behavior of a commonly used ODE model to an analogous cellular automaton model where each lattice site represents a T cell. From the diffusion rate of virions and the rate of encounter and fusion between virus and T cells, the authors derived an expression for the probability that a T cell located a certain distance away from a bursting infected T cell will become infected. The authors initiate their simulation with a single infected cell which in turn infects other cells according to the probability they derived. Subsequent infections result in a radial wave of infection moving outward. This wave leaves behind a pool of infected cells. If the replacement rate of dead cells with new target cells is slow, the wave dissipates and uninfected cells replenish the area. If the replacement rate is fast, the newly regenerated target cells cause the infection wave to recede and a chaotic steady state pattern forms . Most importantly, because infectivity depends on the concentration of T cells (neighbors can be infected more readily when cells are tightly packed), the authors found that the ODE model overestimated viral infectiousness by more than an order of magnitude when compared with the spatial model .

Funk et al.  used an approach similar to that of Louzoun et al. , but proposed a more complete and realistic model of three species: uninfected target cells, virus-infected cells, and virus. The authors first present the dynamics of the non-spatial homogeneous model, and then extend the model to a 2-D spatial model, and finally they randomly vary the parameters from site to site according to a uniform distribution. Comparing results from the non-spatial homogeneous model with the spatial model, the authors found that as the infection spreads outwards from the initial infection site, the viral titer curves for both models are initially similar. However, as the initially infected sites approach an equilibrium state and outer cells become affected by the infection wave, the increase in virus slows. This phenomenon is not seen in the non-spatial homogeneous model. As in the Strain et al. model, the Funk et al. spatial models also resulted in lower average viral concentrations when compared to non-spatial homogeneous models .

When the parameters of the Funk et al. model were varied from site to site, a new dynamic was revealed. This new system formed a series of source and sink sites and the ability of the virus to migrate from site to site (its diffusion rate) had a significant impact on the dynamics. Finally, the authors added immune cells to their heterogeneous spatial model. For this case, they found that a spatial model greatly improved the stability of the infection. Stability was enhanced because the equivalent non-spatial homogeneous model gave rise to large oscillations, analogous to what is observed in predator-prey models in ecology. In contrast, the spatial coupling by local dispersal of virus and immune cells in the non-homogeneous spatial model resulted in oscillations at different sites that were out of phase. These oscillations were damped and had the effect of spatially averaging the dynamics. This is a very interesting finding as the large long-lasting periods of oscillations commonly observed in non-spatial models are only very rarely observed experimentally or clinically, while the equalized dynamics of spatial models are more in line with experimental and clinical observations.

Finally, Beauchemin investigated the effect of viral infection spread in tissue and the particular effect on influenza dynamics. The model is a 2-D square lattice where each site represents a target cell and the grid, which represents the tissue, is patrolled by generic immune cells which can kill infected cells. While the models mentioned above all initiate infection with a single infected cell, this work explored the effect of altering the initial distribution of infected cells. The cases examined included a given number of initial cells distributed either randomly on the grid, in small isolated groups of several cells, or as a single lump. Beauchemin found that infection dynamics was sensitive to the initial spatial distribution. This is due to the fact that when an infected cell is part of an infected patch, it already has infected neighbors and therefore has fewer infectable neighbors than a single infected cell whose neighbors are all infectable. This results in a smaller effective infection rate in simulations where initially infected cells are distributed in larger lumps. The author then explored the effect of target cell replenishment. In ODE models, target cells are typically replenished at a fixed rate or at a rate proportional to the number of target cells as a function of cell division. Beauchemin explored the effects of changing the replacement rule such that an empty site can only be replenished by a new target cell if it is in contact with a dividing uninfected target cell. Screenshots resulting from both regeneration rules are presented in Figure .

As the infection wave spread outward, it left behind a pool of dead cells, which could only be replenished if in contact with uninfected cells. This could only happen once immune cells had detected the infection and began breaching the infection wave that segregates the dead cells from the uninfected cells (see Figure ). This modification to the model yielded a curve for the fraction of dead cells over time which closely resembled experimentally observed curves and contrasted with the sharp and narrow curves obtained with a model where replacement rate only depends on the number of uninfected cells but not on their location.

## Infection of Cell Layers

Because emerging spatial patterns can have a significant effect on the development and outcome of a simulated infection, it is important to define the model in a manner that is consistent with the true spatial nature of the real system. For example, two existing ABMs of HIV are implemented on 2-D square lattices where each site represents a T cell and the infection spreads from infected cells to their neighbors. The infection spreads as a wave and under certain conditions, intricate wave patterns can form. Unfortunately, while it is true that T cells are tightly packed (as in a lattice) , they are quite motile with velocities of approximately $\unit{11}{\micro\meter/\minute}$ . It is likely that if the T cells in these HIV models were motile, these intricate infection wave patterns would not form.

There are, however, models where a 2-D lattice of static cells is the best and most accurate representation. An example of this are *in vitro* experiments conducted on 2-D cell monolayers . Howat et al.  used data from 2-D Madine-Darby bovine kidney (MDBK) and Vero cell cultures infected with Herpes simplex virus 1 (HSV-1) to study the antiviral effect of interferon $\beta$ (IFN-$\beta$). Their ABM consists of a 2-D hexagonal lattice where each site represents a target cell that can either be susceptible, infected, in an antiviral state, or dead. Because IFN-$\beta$ diffusion is very fast, it is taken to be instantaneous and IFN-$\beta$ is tracked as a single global variable. Because the experimental cellular monolayer is overlaid with an agar that impedes viral diffusion, infection is assumed to spread locally: Virions are released by an infected cell as a single burst upon cell death and distributed to the 6 immediate neighbors and 12 next neighbors weighted according to distance. By comparing the fraction of dead cells over time obtained experimentally to what their ABM produces, the authors were able to calibrate their model and determine key parameters of HSV-1 infection and IFN-$\beta$ protection. For example, they found that it takes about 95 hours for the MDBK monolayer to reach an antiviral state . Once calibrated, Howat et al. used their model to determine the dosage of HSV-1 needed to obtain total monolayer destruction, maximum IFN-$\beta$ production, or minimum recovery time. Notably, the authors determined that the minimum recovery time is obtained for the maximum boundary length, i.e., when the dead cell plaque geometry is such that it maximizes the number of cells in contact with an empty area (a pool of dead cells).

Beyond *in vitro* experimental systems, there are also *in vivo* diseases which can be accurately represented by a 2-D lattice of non-motile cells. An uncomplicated influenza infection, i.e., one which infects only the upper lung and does not degenerate to pneumonia, is a good example. Beauchemin et al.  have taken advantage of this fact and developed MASyV's `ma_immune`, a 2-D square lattice model of influenza A infection . This representation is a reasonable description of the true system because the target of influenza A is the tightly packed monolayer of cells that line the airways. Additionally, since influenza virions only bud apically (from the top of the infected cells) the infection effectively occurs in a 2-D plane. Interestingly, when this model was published in 2005, it was only the second model ever proposed for influenza infection in humans; the other one being the 60 parameter ODE model with delays proposed by Bocharov and Romanyukha 20 years earlier . At the time, very little experimental data on the dynamics of influenza was available and a spatial model was especially difficult to justify. Nevertheless, Beauchemin et al. were able to fix 3 of their parameters and set realistic biological bounds on 6 of the 8 other parameters from information gathered in the literature. Then, picking parameter values within these bounds, the authors were able to reproduce a number of features of influenza infection, for example, that the peak of infected cells occurs at 48 and that infection resolves by day $5\pm2$. This work was originally intended to show that the ABM is complex enough to reproduce the general shape of a response to an uncomplicated viral infection, and that it yields quantitatively reasonable results when parameterized for a particular viral infection. The ABM was later used to study the influence of spatial heterogeneities on the dynamical evolution of a viral infection , as discussed in Section 2.3.

This model also improved our understanding of flu kinetics. The model was used to derive values for some parameters (e.g. replication rate of cytotoxic T lymphocytes per infected cell encounter) that were subsequently employed in other influenza models . The model's agreement with the limited available data motivated its expansion into a more detailed and accurate model of influenza infection in vitro. In collaboration with experimentalist Dr. Frederick T. Koster, M.D., and University of New Mexico computer scientist Stephanie Forrest, Beauchemin and others developed an improved version of the model, MASyV's `ma_virions`, which accounts for viral release, dispersion, and clearance . The ABM consists of a 2-D hexagonal grid representing lung tissue and on which viral concentrations evolve based on a discretized version of the diffusion equation. Calibration of the ABM was accomplished by quantitatively comparing the ABM's results with experimental data for fraction of infected cells, patchiness, and viral concentrations versus time produced by an in vitro human lung cell monolayer . The calibration to preliminary data provided an estimate for the diffusion rate of virions that was 1,000-fold lower than the natural diffusion expected from the Stokes-Einstein equation. The calibrated virion diffusion rate is consistent with virions getting trapped by cell receptors and mucus .

## Agent-Based Models in Shape Space

The space in which ABMs evolve does not always need to correspond to 3-D Euclidean space. Another variety of ABMs of host-pathogen/immune systems are called *shape space* models . In these ABMs, the spatial dimensions represent the generalized shapes of cell receptors and epitopes, the parts of antigens recognized by the immune system. The affinity between a receptor and an epitope is a function of their distance in shape space . The evolution of the system in space and time then can be used to represent the development of an immune response including receptor variability due to somatic hypermutation, and the dynamics of clonal expansion and contraction. Typically, an epitope corresponds to a given location in shape space and its presence stimulates the activation and proliferation of lymphocytes located at and around that site, that is, those with a high receptor affinity for that epitope. In modeling viral infections it is also possible to allow the virus to mutate and evolve its epitopes so as to avoid immune elimination. Thus, the shape space framework is very flexible and has obvious use in immunological modeling.

One shape space model that has had an impact on influenza research was proposed by Smith et al.  in order to study the efficacy of repeated annual influenza vaccination. Although annual influenza vaccination is recommended for high risk groups, such as infants and the elderly, clinical studies of the efficacy of annual vaccination have lead to conflicting results. One empirical study by Hoskins et al. suggested that during some flu seasons patients given repeated annual vaccinations had less protection that first-time vaccinees, while another study suggested that repeated annual vaccination did in fact offer long-term continual protection . Smith et al.  showed that differences in the antigenic distance (i.e. the distance in shape space) between the vaccine strain and the circulating influenza strain in both studies, which were conducted in different years, could account for the observed discrepancies in the results. In particular, they showed that a first vaccine will negatively interfere with the protection potentially afforded by a second vaccine when the two vaccine strains are close in shape space. This is because the first vaccination raises antibodies that by cross-reactivity can partially eliminate the second vaccine. Using an ABM that followed the dynamics of a broad array of B cells and antibodies for the vaccine epitopes, Smith et al.  were able to show in a quantitative manner how the model could explain the available experimental data on the efficacy of repeated vaccination. Further, the idea of monitoring the antigenic distance in shape space between flu strains has proved to be a valuable approach in studying influenza evolution and is now a factor in the World Health Organization's biannual decision making process for selecting which vaccine strains to include in the seasonal influenza vaccine .

## The Role of Agent-Based Models in Multiscale Systems

Biological systems naturally span multiple time and length scale hierarchies. Time and spatial scales can range from $10^{-2}$ seconds and $10^{-9}$ meters during gene transcription and protein production at the molecular level to $10^6$ seconds and up to several meters at the level of the whole organism . However, experimental investigations and the data generated focus on specific mechanisms relegated to isolated time or length scales and must be translated between biological scales. A multiscale approach to modeling biological systems is crucial to developing realistic and relevant models capable of predicting complex biological phenomena.

Bridging multiple scales is a formidable modeling challenge. Because ABMs are constructed by considering the behaviors of individual system components and can be developed for each biological subsystem or hierarchy, they are also naturally well suited for linking these different models together. For example, merging ABMs focusing on different biological scales has been used to translate the results of and facilitate collaboration between experimental groups working on various aspects of the acute inflammatory response . In this work, a group of agent based models were developed to simulate intracellular signal transduction pathways that were then incorporated into ABMs for different cell types. Subsequently, these cellular models were integrated to simulate tissue and whole organ function during the acute inflammatory response.

Another modeling technique that has emerged in response to the computational challenge presented by multiple scales integrates differential equation models with ABMs to couple the dynamics occurring on different time and length scales . This technique has been applied to the development of models describing tumor growth, a complex biological system with important clinical applications . In Jiang et al. , a three-dimensional model of avascular tumor growth was developed that spans three distinct scales. A type of ABM, called a lattice Monte Carlo model, is used to describe tumor cell growth, adhesion, and viability. These cell dynamics are regulated at the intracellular level by a Boolean network for protein expression that controls the cell cycle. At the extracellular level, nutrients, metabolic waste, and growth factor and inhibitory chemical concentrations are described by a system of partial differential equations. Growth curve measurements and measurements of the size of the proliferating rim and necrotic core regions from simulated tumors agree with the quantitative results obtained from tumor spheroid experiments. This model predicts what environmental conditions tumor cells require for survival, and the molecular weight of potential growth promoters and inhibitors regulating tumor cell viability. As demonstrated by these results, validated models of biological systems provide new biological insights, can be used for prediction, and guide future experimental pursuits. Looking forward, coupling tumor growth models with agent-based and cell-based models of angiogenesis and the immune system response against the tumor will be a major step towards the development of an integrated systems approach to modeling tumor growth and cancer invasion and towards the ultimate goal of predicting the effects of novel cancer and anti-angiogenic therapies.

## Improved Experimental Data Fueling Advances in Modeling

Models of any type can only be as good as the data used to develop, calibrate and confirm them. Fortunately, new experimental techniques are allowing access to unprecedented amounts of data, but also to novel types of data. Because ABMs are spatial models they require spatial data: a rare commodity. High resolution spatial data is especially important given that localized spatial effects can have an important impact on the resulting dynamics, as mentioned above. The recent application of two-photon microscopy to *in situ* and *in vivo* imaging of mice lymph nodes has determined trajectories of individual T cells as they move within lymph nodes . This new technique has revealed that T cell motion in lymph nodes occurs largely as a result of the cells crawling along fibers of the fibroblastic reticular cell network . Additionally, examination of these trajectories has shown that the end-to-end distance that T cells travel increases approximately linearly with the square root of time, suggesting that the movement of T cell over long distances can be described by a random walk rather than motility guided by chemotaxis. These new data have facilitated the elaboration of models for the movement of T cells within lymph nodes: some simple and some more detailed .

The model by Beauchemin et al.  was a simple model where a T cell moves in a straight line at constant speed $v_\mathrm{free}$ for a time $t_\mathrm{free}$, then pauses for a time $t_\mathrm{pause}$ as it reorganizes its intracellular machinery allowing it to turn. It then picks a new direction, and undergoes another free run. Comparing simulations against experimental data yielded a mean free path of approximately $\unit{38}{\micro\meter}$, which is about twice the distance between intersections in the fibroblastic reticular cell network, which averages $\unit{17\pm7}{\micro\meter}$. From this, the authors concluded that when coming to a fiber intersection, a T cell will turn roughly $50\%$ of the time, and $50\%$ of the time it will continue along its original fiber.

The Meyer-Hermann and Maini model offers a more detailed description of the active crawling exhibited by both B and T cells. This model is a 2-D cellular Potts model where a B or T cell is represented by a group of neighboring lattice sites. The lattice sites that make up the lymphocyte are moved individually but their movements are constrained to maintain a constant cell volume and to respect the general direction chosen for active movement. The model yielded a distribution for the velocities of B and T cells that agreed impressively well with velocities obtained experimentally .

More recently, Beltman et al.  proposed a 3-D rather than a 2-D cellular Potts model, which also represented lymphocytes as a collection of neighboring lattice sites. The authors investigated the possibility that motility patterns observed for T cells and dendritic cells are a result of the anatomical structure of lymph nodes, rather than an intrinsic characteristic of the cells' motility. Indeed, they found that by restraining T cell movement using the constraints imposed by anatomical barriers in lymph nodes, the model can account for all characteristics of the experimental observations. In particular, they showed that lymph node anatomy can explain the large velocity fluctuations seen experimentally for T cells. Further, they found that as a result of the obstacles, T cells tend to organize into small, highly dynamic streams and provide experimental evidence to support these new findings. Finally, using their model, they were able to predict that dendritic cells would be able to contact approximately 2,000 different T cells per hour while T cells could contact roughly 100 dendritic cells per hour .

## Sensitivity and Uncertainty Analysis in Modeling and Prediction

Because ABMs consist of several thousand lines of code, it can often be just as complicated to understand the model as it is to understand the real system. In general, as the system grows in complexity, the number of parameters that must be estimated also increases. Parameter estimation introduces uncertainty into the system due to experimental error, differences in experimental assays, or error from data fitting techniques. Moreover, when the number of agents is large and their interactions numerous and complicated, it can become difficult to extract or isolate the key processes responsible for a given outcome. In their paper proposing an ABM of *M. tuberculosis* infection, Segovia-Juarez, Ganguli, and Kirschner use a simple and efficient way of analyzing their results by performing a sensitivity and uncertainty analysis of the model's parameters. The authors picked 12 of the 27 parameters of their model and investigated their effect on the growth of tuberculosis granulomas. Using the number of extracellular bacterium and the granuloma size as outcome variables, they mapped out the parameter space using the Latin hypercube sampling method and quantified the sensitivity of their model to these parameters using partial rank correlation coefficients on each parameter set. This type of analysis quantifies the importance and impact of each parameter on the outcome variables. Their analysis revealed that the intracellular growth rate of the bacteria is strongly and positively correlated with the number of extracellular bacteria at early times of the infection and at much later times post-infection, and negatively correlated at intermediate times. This suggests that at low bacterial levels, a large intracellular bacteria growth rate is critical to infection, while a small growth rate seems necessary to generate larger lasting granulomas. This is one example of how simulation and analysis can be coupled to provide insights into and enhance our understanding of the dynamics of disease (see for a more complete review of these analytic techniques).

# Discussion

Agent-based models have inspired significant interest through their visual appeal and because the languages used to describe the model and the natural system are very similar. Computational advances are making possible the use of ABMs to describe whole systems arising not only in the human immune system, which has been the focus of this article, but also in financial markets , the spread of epidemics , cancer dynamics , and the threat of bio-warfare to name just a few. In addition, ABMs have engaged experimentalists and have helped facilitate their increasing involvement in modeling. Here, we have presented a limited survey of ABMs in the context of host-pathogen dynamics. This survey is in no way exhaustive, but is meant to highlight the novel and original modeling approaches applied by researchers and the contributions to our knowledge of specific diseases and immune mechanisms.

ABMs are only one of a number of modeling techniques that could be applied to describe a particular system. The choice of model to employ is guided by multiple factors including the questions that one wishes to answer, the assumptions that can be made, parameter availability, and computational expense. Choosing agent-based modeling is not without cost or consequence. ABMs tend to require a larger number of parameters than both ODE models and their spatially-continuous partial differential equation analogs. Until experimental assays capable of measuring relevant experimental data catch up with the needs and complexity of model calibration, ODEs are sometimes the only viable modeling option. An important weakness of ABMs arises in the context of sharing and comparing results. At worst, a partial differential equation model with delays will consist of a large and complicated set of equations, but that can be written down in an unambiguous manner. In contrast, an ABM consists of many lines of codes which could be in any variety of programming languages, and its translation into words can overlook very important details of its implementation that are crucial to reproducing its results. This is a problem that researchers utilizing ABMs need to pay more attention to. Establishing a permanent home for any computer model developed and a static version corresponding to the exact state of the program used for published results are important first steps. Software repositories like SourceForge.net (<http://sourceforge.net>) and freshmeat (<http://freshmeat.net>) provide such services free of cost.

As with any modeling approach, ABMs need to adhere to the principles of scientific rigor. For an ABM to be useful to understanding the mechanisms of the system, it needs to convincingly relate to it. Models need to be calibrated and researchers should confirm that their parameter values are within physically relevant ranges. Toward this end, it is helpful to provide a sensitivity analysis of the parameters to identify key parameters and to fully grasp the consequences of parameter uncertainty and variability on the observed outcomes. Modelers should also make every effort possible to quantitatively validate the results of their simulations with independent experiments or with reports in the literature, using data different from those that were used for model calibration. Finally, and perhaps most importantly, models should strive to make experimentally verifiable predictions. Once verified, such predictions will establish that the model realistically captures some aspects of the system. Further, models validated in this way may offer new perspectives, provide a framework for formulating and testing hypotheses, and may suggest important new experiments.

[^1]: ALB and CAAB both contributed equally to this manuscript.

[^2]: Corresponding author. Mailing address: Los Alamos National Laboratory, MS-K710, T-10, Los Alamos, NM 87545, USA. Phone: +1 505 667-6829. Fax: +1 505 665-3493. E-mail: `email@example.com`.


abstract: We comment on the recent paper (hep-lat/0102003) by Horváth, Isgur, McCune, and Thacker, which concludes that the local chiral structure of fermionic eigenmodes is not consistent with instanton dominance. Our calculations, done with an overlap action, suggest the opposite conclusion.
address:  Department of Physics, University of Colorado, Boulder, CO 80309 USA
author: Thomas DeGrand, Anna Hasenfratz
date: 2024-10-03
title: Comment on "Evidence Against Instanton Dominance of Topological Charge Fluctuations in QCD" by I. Horváth *et al.*

# Introduction

In a recent paper, Ref. , Horváth, Isgur, McCune, and Thacker present evidence from lattice simulations claiming to show that the behavior of fermionic eigenmodes in the QCD vacuum is inconsistent with an instanton picture. They refer to a recent paper of ours, Ref. , and remark "Although their qualitative conclusions regarding the instantons differ from ours, there does not seem to be any direct conflict between their lattice results and ours." We disagree with that statement. To expose that disagreement, we reanalyzed a subset of our data along the lines of Ref. .

Why is it important to resolve this disagreement? There is an apparent inconsistency, first described by Witten in 1979 , between large number of colors (large-$N_c$) dynamics and instanton-based phenomenology. We do not wish to summarize arguments favoring one or the other approach, which are well-described in a variety of publications . In principle, lattice simulations might distinguish between the qualitative features of the two phenomenologies. However, the lattice calculations of Refs. and produce contradictory results.

The two simulations differ in several ways. Both are done in quenched approximation, using the Wilson gauge action. Ref. works at a lattice coupling $\beta=5.7$ and Ref. works at $\beta=5.9$. The lattice spacings $a$ are at nominal values of 0.17 and 0.11 fm, respectively, from the Sommer parameter, using the interpolating formula of Ref. . The fermion action of Ref. is the Wilson fermion action. This action has rather poor chiral properties, as evidenced by its large additive quark mass renormalization and the presence of exceptional configurations. Formally, it has $O(a)$ scaling violations. Its real eigenmodes are not eigenfunctions of $\gamma_5$, and the expectation value of $\gamma_5$ of a real eigenmode is a function of the eigenvalue.

The action of Ref. is an overlap action , which is constructed by using an approximate overlap action as its starting point. This action has exact chiral symmetry. There is no additive mass renormalization of the quark mass and there are no exceptional configurations. Zero-eigenvalue eigenmodes are chiral eigenstates. Nonzero eigenvalue eigenmodes have zero expectation value of $\gamma_5$. The action has only $O(a^2)$ discretization artifacts.

Our action uses "fat links." Fat link actions replace the usual one-link gauge connection with a combination of several gauge paths. This combination might or might not be projected to SU(3), but in any case it is gauge invariant and local. Locality assures that the fat link action is in the same universality class as the thin link one. Fat links improve the chiral behavior of nonchiral actions , and their use in our implementation of the overlap action is simply to make the calculation of the overlap operator more efficient (a factor of roughly 20 for the action of Ref. , compared to the use of the thin link Wilson action in the overlap.) The particular choice of fattening is not too important. In Ref. , we used an APE-blocked link as the gauge connection. In Ref. it is shown that the perturbative effect of N levels of APE-smearing with a smearing parameter $\alpha$ is to multiply the vertices by a form factor $(1 - \alpha {\hat q}^2 /6)^N$ with $q$ the gluon momentum and $\hat q^2 = 4/a\sum_\mu \sin^2 q_\mu a/2$. In coordinate space this corresponds to a Gaussian smearing with a spread $\langle x^2 \rangle = \alpha N/3 a^2$ of the quark gluon vertex. Some readers might be concerned that several levels of APE smearing might adversely affect our results. To address that point, we describe results from another fattening, one chosen to lie rigorously within a hypercube . We observe no qualitative change in our results.

# Analysis of Fermion Eigenmodes

The statement of the problem has been given by Ref. : "If an instanton-dominated picture of low eigenmodes is at all valid, we would expect the peaks of the wave function to closely resemble instanton zero modes. Thus, if instantons dominate, a local peak in the wave function for a low-lying eigenmode (a fermion lump) should be dominantly a lump in $\psi_L^\dagger\psi_L$ or $\psi_R^\dagger\psi_R$, but not both." (Here $\psi_L^\dagger\psi_L = \psi^\dagger (1-\gamma_5)\psi$ and $\psi_R^\dagger\psi_R = \psi^\dagger(1+\gamma_5)\psi$.) "On the other hand, fermion lumps without a definite chirality would be an indication of non-self-dual fluctuations, as expected from a confinement-related mechanism."

In Ref. we used an overlap action with $N=10$ levels of APE smearing with $\alpha=0.45$ APE parameter. We measured the autocorrelation function of the local chirality density $\omega(x) = \langle \psi(x) | \gamma_5 | \psi(x) \rangle$ and the correlator of $\omega(x)$ with a pure gauge observable sensitive to topological charge $Q(x)$. We saw that both kinds of correlators were strongly peaked at small operator separation indicating the the chiral density has spatially localized lumps and those lumps correlate with the lumps of the topological charge density. However, there is a logical possibility that a lump in $\omega(x)=\bar\psi_L\psi_L-\bar\psi_R\psi_R$ does not indicate a lump in either $\bar\psi_L\psi_L$ or $\bar\psi_R\psi_R$. To check that possibility, we follow Ref. and define a "chiral order parameter" $X(x)$ $$\tan({\pi \over 4}(1 + X(x)) = {{|\psi_L(x)|}\over{|\psi_R(x)|}}
 = ( {{\psi_L^\dagger(x)\psi_L(x)}\over{\psi_R^\dagger(x)\psi_R(x)}} )^{1/2}$$ and plot histograms of the density of the variable $X$, $\rho(X)$, for regions of the lattice where $\psi(x)^\dagger\psi(x)$ is large, a few per cent of the points of the lattice. We distinguish two cases:

\(a\) Zero modes are chiral. The analogs of Figs. 1 and 4 of Ref. for our action are just delta functions $\rho(X) \simeq \delta(X \pm 1)$. The chiral symmetry breaking inherent in the Wilson fermion action of Ref. broadens this distribution.

\(b\) Nonzero eigenvalue eigenmodes are nonchiral, but we find that the distribution $\rho(X)$ is strongly peaked near $X = \pm 1$. To show this, we take a set of six $12^4$ lattices from our data set and plot histograms of $\rho(X)$. Fig. shows the result for the lowest two nonchiral eigenmodes for the overlap action, with a cut keeping the top 2.5 per cent of $\psi^\dagger\psi$. The imaginary part of the eigenvalue of the Dirac operator for these modes is in the range $0.03/a$ to $0.09/a$. Keeping more lattice points in the histogram will fill in the dip. That simply implies that as we keep more and more lattice points we pick up more of the vacuum fluctuations around the chiral lumps and those eventually overwhelm the chiral structure. This happens very slowly with the two smallest non-chiral modes. Even with 30 per cent of the lattice points the two peak structure is clearly visible.

Figs. and shows the result for the largest modes we have recorded, the ninth and tenth eigenvalue modes of the squared Dirac operator. The imaginary part of the eigenvalue of the Dirac operator for these modes is around $0.20/a$-$0.26/a$. Figure , where about the top 0.3 per cent of $\psi^\dagger\psi$ is kept in the histogram, shows strong peaks around $X=\pm 1$. The peaks disappear as we add more points to the distribution. With 2.5 per cent of the lattice points kept the histogram is almost flat as Figure shows.

Lower eigenvalue eigenmodes show stronger peaking than higher eigenvalue modes, but all the modes we examined show peaking. We remind the reader that we showed in Ref. that hadronic correlators for light quark masses are well approximated by quark propagators containing only a few eigenmodes.

Next we consider results obtained with an overlap action that uses more local fat links. The fat links of the hypercubic blocking mix links only within a hypercube but achieve almost the same level of smoothness as the APE smearing considered previously.

The analog of Fig. from eigenmodes of the overlap action with hypercubic fat link is shown in Fig. . We used the same $12^4$ configurations and kept the same fraction of lattice points as in Fig. . The shape of the curve is qualitatively similar to that from APE-blocked links. Peaks of fermion density are again chiral.

# Conclusions

Our studies with a chiral fermion action at a lattice spacing near 0.11 fm showed that low eigenvalue fermionic eigenmodes have a structure strongly correlated with the locations of instantons and anti-instantons. This note demonstrates that in the locations where the fermionic modes are largest, the modes are also chiral. It seems to us that the simplest description of what we see is in terms of an instanton liquid model of the QCD vacuum.

We do not doubt that the numerical results of Ref. are correct, but we are concerned that the combination of a larger lattice spacing and the use of a nonchiral fermion action introduces lattice artifacts which compromise extrapolation to the continuum limit.

# Acknowledgments

We acknowledge conversations with I. Horvath and H. Thacker about the points we have raised in this paper. This work was supported by the U. S. Department of Energy.


author: Miguel Ángel Castro Tirado; Alberto J. Castro-Tirado
date: Received May 4, 2019; accepted June 11, 2019
title:  The Evolution of Astronomical Observatory Design

# Introduction

The theoretical approach to astronomical observatories has its origin in the roots of the modern observatory . However, historical-descriptive works that consider a particular project are the most common during the next centuries. In addition, some studies appeared focused on specific periods , others with a more general retrospective nature and only a few works focused on the architectural-constructive aspects of the observatories .

The concept of astronomical observatory itself lacks a precise and globally accepted definition given that historically there has not been a critical-intellectual movement that theorized about its nature or conception. Although it may seem obvious, what is understood by an observatory has changed over the years, and, in addition, it depends on the perspective from which it is valued. In this paper, the intention is to present an approach from the point of view of the architecture of the building itself, so that other aspects such as its astronomical equipment or its possible scientific usage will be set aside. That is to say, we will study the particularities related to the construction in which the instruments are located and whose purpose is the astronomical function. The analysis will also be limited to optical observatories.

Talking about evolution necessarily implies the extension of the study over time to analyze and assess how, when and why the changes that have defined the development of these astronomical centers have occurred. For this reason, following an earlier work , starting from the examination of plans, ruins and the buildings themselves, this study goes back to the genesis of the observatories trying to illustrate the reciprocal relationship that exists between the design strategies of the building and the practice of its astronomical function, with the goal of theorizing about their future development. We have focused on 32 observatories which we consider to be representative ones in the period 1259–2007, deeply analyzing the architectural drafts.

# Origins

Astronomy is one of the oldest scientific disciplines of mankind. Since man raised his gaze to the sky and contemplated the stars with his own eyes, a relationship of fascination and devotion was established. This would be reflected in one way or another in different cultures that would arise across the planet. And even without having tools of precision or a basic understanding of the physical laws that rule the functioning of the universe, the first civilizations already had certain notions of the main astronomical facts. All based, always, on mere observations.

In any case, it is not possible to determine the date on which astronomy arises, although its roots probably extend to the origins of the human societies. What does seem clear is that, regardless of a spiritual or pseudo-religious sense, the beginnings of this science point towards an instrumentalization of the sky as a system for measuring periods of time and as a mean of orientation. Beyond the numerous artistic representations of the Sun, the Moon and the stars that have been found , sometimes more realistic and others more allegorical or figurative, or of the written testimonies that have lasted over time, this astronomical predilection remains through the incidence of some celestial phenomena in the disposition of certain archaeological sets.

Although certain remains of ancient civilizations show an impact of astronomy on their culture, it will be the proto-observatories that stand out as points of interest to know the way in which these peoples approached the practice of this science. In these spaces the astronomical observatories have their deepest antecedents and roots. Although the astronomical principles that underlie the constitution of these assemblies are undoubted, there is no evidence to support the existence of research or scientific work in these places.

The oldest site that has been recorded so far is the Goseck Circle, in Germany, a solar proto-observatory with an annular plan that consists of a Neolithic structure dating back almost 7000 years . Of a similar period is the megalithic group of Nabta Playa, in Egypt. There are many other cases of lesser antiquity, among which the temple of Mnajdra, in Malta, or the famed Stonehenge, in England . All these sites belong to ancient societies, with a low degree of technology and even, frequently, of nomadic habits. These facts propitiate relatively simple identifications of their constructions. In addition, the scale and materials of these structures favor their durability up to the present, allowing them to be studied and identified, where relevant, as proto-observatories.

However, the evolution towards more advanced civilizations tend to create settlements of multiple buildings that do not reach the precise degree of specialization for the observation function as such. Since no scientific instruments had been developed that required demanding conditions, systematic surveying of the sky could be carried out with the naked eye in almost any construction or, even, in the vicinity of these out in the open.

The oldest building strictly dedicated to astronomy of which there is evidence is Cheomseongdae, in South Korea. It was designed as a small bottle-shaped hollow tower made of masonry. Its function was to enable the use of astronomical instruments from an elevated position to avoid surrounding obstacles. It dates from the seventh century and remains in a great state of preservation and its square top and round body reflect the astronomical concept of the time (i.e. the Earth was considered a square and the sky was considered to have a round shape) with the 29 layers of stone from the bottom to the top corresponding to the 29.5 days of full lunation (or synodic month, i.e. the time from one new moon to the next) and the 27 layers of the round body representing the 27.3-day sidereal month, amongst other astronomical symbolism .

Anyway, although references to the Tower of Babel or to some type of astronomical institution directed by Ptolemy in Alexandria can be found , there is no evidence to support the existence as such of pre-medieval astronomical observatories.

# The Medieval Islamic Observatory

Coinciding with the beginning of the Middle Ages, a period of cultural darkness began in the Western world from which astronomy would not escape. During this period, Islamic culture came to play an essential role in the history of astronomy by preserving the vast Greek knowledge of this subject and incorporating its own findings. All this would be recovered for the West by the School of Translators of Toledo of Alfonso X, during the thirteenth century .

The medieval Islamic observing sites can indeed be called observatories because of two main reasons: i) the fact that they were spaces dedicated specifically to the collective and prolonged study of celestial phenomena and ii) to the exchange of scientific knowledge. It can be said that the observatory, as such, appears for the first time in this civilization. This is not a coincidence; the features of the Islamic culture as a society link with the observatory as a scientific institution . However, although this permanence is implemented in the Islamic observatory, they would have a reduced life .

Two singularities that would ultimately be decisive for the relevance of the Islamic observatory as an institution were the size of the observing tools and its patronage. The fact that the tools of observation and study reached a dimension such as to prevent their portability led to a necessary settlement. Royal or state support meant the formation of institutions larger and longer lived than an individual scientist .

Although several observing posts, more or less temporary, had previously been established, it is difficult to speak of observatories as an institution before Al Mamûn (ninth century). The numerous observations collected by different authors allow us to recognize centers established in Shammasiya, Baghdad and Mount Qasiyun, Damascus . However, the fact that there are no remains, descriptions or illustrations of observatories could mean that no specific space was devoted to this use. This lack of material legacy makes it extremely difficult to distinguish between the observatories of this period and the specific observing posts established for an ephemeral study. The main differences amongst medieval islamic observatories were in the number of personnel dedicated to them, the quality of observing equipment and the life span of these centers. That is to say, at this initial moment the observatories were not a clearly defined institution and it is complicated to indicate which observation sites transcended the barrier between the most advanced temporary work post and the embryonic scientific center. It would be from this stage onward that the Islamic observatory began to develop as a specialized astronomical institution linked to a scientific team and a specific location. It then begins to be structured with its own administration around a work program and its instruments become more precise, larger and heavier.

All these circumstances suppose a change in the paradigm of the Islamic observatory that, in a hesitant process over several centuries, transcends from a semi-improvised observation post to the scientific reference center of its time. The construction of the observatory of Maragheh, in northwest Iran, was completed in 1264 and is part of the period that represents the climax of the evolution of the Islamic observatory . In fact, the main progress and features of this institution in the Islamic world appear with this observatory, which became the paradigm to imitate of its time.

Although only the ruins of the construction remain, these together with the existing testimonies and the subsequent observatories built in its image, allow a description of the building to be developed. The observatory was a circular wall erected around a large central space occupied by the large quadrant for meridian transit observations. On both sides of this area there would be spaces to house the library, workshops, stores and other functions distinct from mere observation. Even though the date of the fall of the observatory at Maragheh is not exactly known, it still became the model to imitate by the great Islamic observatories that would arise during the following centuries, such as those of Samarkand and Istanbul.

# The Modern Observatory

Although some references attribute the first European observatory to either Regiomontanus or Copernicus , there is no evidence of a construction or of fixed instruments that support these claims . There is a consensus to place the first modern observatory in the palace of Wilhelm IV of Hesse-Kassel (this is recognized by or ). However, a platform built as a lookout point in the palace where portable instruments are taken out does not represent an advance from the perspective of architecture. This does not invalidate its importance since it was a visit to these facilities that inspired Tycho Brahe, a few months later, to build the first two modern observatories in Hven characterized by having a certain degree of intention and architectural specialization for astronomy.

The first project of Tycho Brahe would be Uraniborg (1576), a palace in which the study of alchemy and astronomy would coexist with his own home. A dialogue between science and architecture in which the spatial configuration would adapt to the needs of research resulted in a solution that relegated alchemy to the basement and elevated astronomy to the upper floor, leaving the work and residence areas in the middle floor. Elevated platforms were built on slender poles from which it was possible to look at the sky without obstacles with covers that protected both the instruments and observers from poor weather conditions.

After several years Tycho Brahe noticed that the vibration of the platforms was being transmitted to the instruments and damaging the research unacceptably, so he planned a second center in which neither structural nor thermal instability would compromise astronomy. Following this idea Stjerneborg (1584) arose as a building excavated in the earth in which the observation tools would be fixed on pillars founded on the bedrock itself.

The appearance of the telescope in the early seventeenth century and its early connection with astronomy would condition the future design of observatories, in which the increase in accuracy would require improvements to avoid architecture from becoming a brake on scientific research. However, since the first telescopes were small and manageable, many astronomers could make their observations from the windows, balconies, terraces or gardens of their own residence, or from public spaces. Possibly, this is why, during this period, the development of astronomical centers tends towards simple structures or, directly, to adding platforms to existing buildings from which to get better views of the sky, as in the universities of Leiden (1632) and Ingolstadt (1637) or in the Danish Rundetårn (1637) . Therefore, many years would pass between the appearance of the first modern observatory and the consolidation of the idea that a center for scientific observation required a building expressly designed for that purpose.

The observatory of Paris (1667) was created with the support of Louis XIV as a national institution in a context of increasing interest for culture and science. In addition to its scientific purpose, the project was to accommodate the headquarters of the Academy of Sciences and to enhance the prestige of its promoter. The resulting building is a palace where the institutional function shares space with the residences and the areas dedicated to astronomy: observation rooms, a zenith study room in the basement and a flat walkway for instruments. For more than a century, observatories followed the example of the French, limited to more or less contrived buildings that combined residence, work area and spaces where telescopes could operate. This is how the high observation rooms with vertical windows were reproduced, such as at the Greenwich Observatory (1675) or the openings from north to south were implemented at the Radcliffe Observatory (1772).

It was not until the construction of the Dunsink Observatory (1785) that significant changes were assimilated, where the architectural configuration was defined by the astronomical instruments that were to be installed. The person responsible for this project defined the key elements of the building as "situation, foundation and soil" . So the location was selected for its good visibility and accessibility, and the importance of stability meant the construction of foundations for the instruments to become structurally independent of the rest of the building. Systems were also developed to promote thermal stability in observation rooms. The design even reserves a privileged position to the telescopes, placing them under movable domes. The modern observatory is designed and planned to respond to a specific astronomical use during the period between the seventeenth and eighteenth centuries. The new observatories implement changes derived from the development of astronomy itself as a science (especially in relation to its instruments) and improve its performance for the sake of an optimal performance of the research activity without underestimating the welfare of its occupants.

Sometimes for strategic purposes (because of its importance for maritime navigation), sometimes for cultural ambition and in others, perhaps, by imitation, these buildings began to emerge throughout Europe in the eighteenth century, from where they would end up spreading to the rest of the world. Spaces for observation, characterized by the instruments that occupied them, coexisted with spaces in which the activities derived from these observations would take place (calculation, data analysis, discussion of results ...), secondary ones (preparation of instruments, offices, classrooms, libraries ...), and others corresponding to service uses (residences, storage areas ...) that, on occasion, moved to one or more independent buildings. Likewise, the scientific and technological progress of this period had an impact on the observatories: on the one hand the buildings incorporated new functions such as physics or chemistry laboratories, photographic work rooms or optical and mechanical workshops, while on the other hand the work of astronomers was facilitated by mechanizing the movement of domes and telescopes. In addition, technical advances allowed a growth of the refractor that led to a race to having the largest telescope in the world. This struggle reached its peak at the Yerkes Observatory (1897) where the size of the refractor (1.02 m diameter) entailed an increase in the architectural scale of the building .

Although in this period the trend led to the observatories to move away from the city and, therefore, the town, some centers would still emerge, such as the Allegheny (1912) or the Griffith (1925), which would open to the public offering the use of some instruments and other activities to disseminate and promote astronomy.

```latex
\begin{table*}[t!]\caption{Percentages of surface areas of observatory buildings dedicated to different functions.\label{tab:jkastable1}}
\centering
\begin{tabular}{lrrrr}
\toprule
Observatory        & Area (m$^2$)  & Observation (\%) & Derivatives (\%) & Accessories (\%) \\
\midrule
Maragheh (1259)     &  1164  &  4  &  96  &  0   \\
Samarkand (1424)  &  3082  &  3  &  97  &  0   \\
                   &        &     &      &      \\
Uraniborg (1576)   &  1180  &  9  &  22  & 69   \\
Stjerneborg (1584) &    91  & 71  &  29  &  0   \\
Paris (1667)       &  3152  & 48  &   8  & 44   \\
Greenwich (1675)   &   501  & 25  &   4  & 71   \\
C\'adiz (1753)     &  1364  & 33  &  39  & 28   \\
Madrid (1790)      &   618  & 27  &  14  & 59   \\
Radcliffe (1772)   &   643  & 36  &  32  & 32   \\
Dusink (1783)      &   802  & 21  &  27  & 52   \\
                   &        &     &      &      \\
Seeberg (1790)     &   473  & 13  &  4   & 83   \\
G\"ottingen (1816) &  1092  & 16  & 16   & 68   \\
Hamburgo (1825)    &   627  & 18  & 19   & 63   \\
Cape of Good Hope (1828) & 355 & 26 & 3  & 61   \\
Helsinki (1834)    &  1097  & 23  & 16   & 61   \\
Harvard (1839)     &   744  & 29  & 20   & 51   \\
Pulkovo (1839)     &  1429  & 44  & 18   & 38   \\
U.S. Naval in Washington (1844) & 367 & 26 & 40 & 34 \\
Lisboa (1861)      &   893  & 40  & 10   & 50   \\
Yerkes (1897)      &  2396  & 28  & 40   & 32   \\
Allegheny (1912)   &   969  & 44  & 26   & 30   \\
                   &        &     &      &      \\
Lick (great refractor, 1881)  &  957   &  43 & 32   & 25   \\
Niza (large equatorial, 1881) & 1040 & 43 & 43 & 14 \\
Mount Wilson (60-inch, 1908)  &  975 & 52 & 48 &  0 \\
Mount Wilson (100-inch, 1917) & 2509 & 51 & 49 &  0 \\
Hamburg-Bergedorf (large refractor, 1912) & 562 & 32 & 56 & 12 \\
Palomar (Hale telescope, 1948)& 1443 & 28 & 63 & 9 \\
Lick (Shane reflector, 1959)  & 2108 & 33 & 60 & 7 \\
SAO (BTA-6m, 1975) &   4480   & 33   & 60 & 7 \\
Calar Alto (3.5m, 1975) & 3709 & 18 & 74 & 8 \\
Keck (1993)        &   3873   & 49  & 47 & 4 \\
Gran Telescopio Canarias (2007) & 3457 & 24 & 70 & 6 \\
\bottomrule
\end{tabular}
\end{table*}
```

# The Contemporary Observatory

The development of cities from the nineteenth century began to cause problems in some nineteenth-century observatories, which not only began to withdraw from large cities to smaller settlements but also to seek locations far from urban areas. This withdrawal became increasingly necessary in the following decades, underlining a physical distance from the ordinary citizen that also led to an intellectual disconnection from astronomers until the press/media started to recover the connection amongst them in the 1960s. Astronomical centers with spaces such as conference halls, classrooms or exhibitions, open to visitors, amateurs and experts, or even public observatories, lost prominence in favor of highly specialized scientific complexes almost exclusively for astronomy.

The search for the highest possible precision from the instruments resulted in optimization of the observation conditions, which led to the establishment of these institutions in areas not affected by anthropogenic (light) pollution and with more favorable atmospheric conditions. The Lick Observatory (1881) was the first continuously inhabited mountain observatory in the world . The location of the center was fundamentally based on its astronomical advantages, without regard of the distance to any city or university . The architectural characteristics of the building are more like those of modern observatories than other contemporary observatories, which is why the Lick Observatory represents an intermediate, transitional case between the two epochs.

Almost simultaneously, the Nice Observatory (1881) emerged as a set of independent structures scattered over the top of a mountain. Each building was designed and located for a specific function, including the first building devoted entirely to a telescope. In this way, it pioneered the contemporary observatory layout of a set of specialized and independent buildings .

The size of the new large reflectors, the differences in scale, the different needs intrinsic to the multiple uses or the orography of the mountainous terrain that makes it difficult to have wide level soils, are aspects that drive the segregation of the initially unified observatory building into a series of independent structures. This separation allows for a freedom of disposition and orientation in the territory that will favor the conditions of observation, minimizing interference (vibrations, smoke, light, heating) that would compromise the quality of the astronomical observation.

Although mirror-based telescopes had been installed in observatories for decades, their role had always been secondary to that of large refractors. However, their importance increased thanks to technical advances that enabled the manufacture of the required mirrors in parallel with the development of astrophysics. In addition, the dissociation of the observatory buildings allowed the adaptation of the constructions to the instruments they housed without external constraints. An example of this is the Mount Wilson Observatory (1904), initially harboring the large 60-inch reflector and, later, in the 100-inch Hooker telescope .

Ignoring the heterogeneous large-scale layout of the modern observatory, the contemporary design for an individual instrument is generally a base architectural structure scaled to the instrument it contains, topped by a movable cover, generally a dome. The size of telescopes reached such magnitudes that it rivaled the very scale of the architecture that had to contain it. This meant increasing the size of the buildings that housed the telescopes, producing a large amount of contained space available for other uses. Similarly, these large instruments require major structural supports, to such an extent that spaces inside the supports themselves can accommodate some functions.

Following the earlier Hooker telescope development, the Hale Telescope (1948) ended up consolidating different functions and hosting almost a thousand people during its opening ceremony. In this way, these large telescopes would incorporate or could incorporate enough secondary functions (offices, work canteens, meeting rooms, laboratories, warehouses, workshops, rest rooms, toilets ...) to be able to contain practically everything required for an entire observatory. Thus, the construction planned for each instrument would have sufficient autonomy to constitute an astronomical observatory by itself.

From this point on, the most important observatories were constituted as sets of independent structures where, habitually, they would include buildings for some accessory uses (residences, visitor areas ...), other structures allocated to secondary instruments, and the main buildings housing the great telescopes together with the primary functions (meeting rooms, work rooms, laboratories, offices ...), which would constitute the essence of the observatory.

The constant growth of the mirrors caused mechanical-structural problems in the BTA-6 (1975) (with a 6-meter primary mirror) that forced to introduce some fundamental changes in its design, adopting an alt-azimuth instead of an equatorial mount. This variation complicated the tracking of celestial objects, requiring a computerized control system. These two revolutionary advances were consolidated and replicated in later observatories . After the largest monolithic mirrors (8.2-meter diameter) were produced for the European Southern Observatory's Very Large Telescope (VLT, first light in 1998) and Subaru (8.2 m), and the Large Binocular Telescope (8.4 m) later on, the planners of the W. M. Keck Observatory (1993/1996) decided to design the primary mirror as a combination of smaller mirrors to make ever-increasing mirror sizes possible. Correctly aligning segmented mirrors at that time represented a very hard problem. Because of this, for VLT a single-mirror approach was selected. Keck did indeed have control problems for years, vindicating the VLT decision. However, the experience of fixing Keck's problems paved the way for the future larger segmented mirrors telescopes.

The alt-azimuth design allowed for a significant reduction of the structures, which could be light, and the size of the domes, which would also reduce their cost compared to equatorially mounted telescopes. The race to the largest reflector paved the way to the Gran Telescopio Canarias (GTC, first light in 2007), with effective aperture of 10.4 m (composed of 36 segmented mirrors of 1.8 m diameter each) and inspired by the Keck. However, other projects are already underway that will exceed it in size but without significant changes in its configuration: the GMT (24.5m), the TMT (30m) or the ELT (39m).

Despite the changes experienced during this period, or even from the very origin of this institution, the importance of technological development in terms of computerization and communication systems must be highlighted as the causes of the greatest observatory revolution. Both the remote control and the remote access to the observations mark a decisive turning point regarding these complexes, since it seems to offer the possibility of abandoning all physical presence in the observatories. And although it is not necessary for an astronomer to look through a lens to access an observation, there are still a series of tasks derived from it that require a space in which scientists and engineers can work.

# Space-function Analysis

Beyond this summary of the development of the observatory up to the present, anticipating the future evolution needs some analysis. Since it is difficult to evaluate all the variables that affect the constitution of these astronomical centers, this study focuses on the configuration of its architecture in relation to its uses, considering the evolution of astronomy and astronomical observing as a function of time. For this purpose, a classification is established for the different activities for which the different spaces of the buildings that house the instruments are intended. Observation spaces are considered to be those places occupied by active observation equipment, that is, they can be in operation if the environmental conditions make it possible. Research support spaces are those space that, though not directly participate in observations, provide functions necessary for astronomical work. Included in this category are work rooms, workshops, instrument warehouses, and so on. Finally, all of the dependencies that are not essential for astronomical research, such as conference rooms, rest areas, visitor galleries or exhibitions, are described as service spaces.

Based on this classification, we selected a representative sample (Figure ) of astronomical centers that had the greatest impact and influence on other observatories of their time and afterwards. They are grouped according to the periods described in the previous sections and ordered chronologically and summarized in Table , which is based on the analysis of the plans and drafts of the different observatories together with information on how various spaces were used.

Our work shows how the aspects defined in the different periods have their corresponding manifestation in this classification. This highlights the weak correlation between the Islamic observatory and the modern observatory. In addition, during the beginnings of the modern observatory (between the Uraniborg and Dunsink cases) there is a period of continuous changes in the distribution of uses, probably related to the different types of instrumentation and observing techniques used. However, since its consolidation (since Seeberg) a tendency has emerged regarding the development of its organization that extends to the current contemporary observatory.

On the one hand, Figure  shows a tendency toward an increase of the fraction occupied by research support spaces over time. That means a higher proportion of area in which astronomers and technicians can actually work (offices, laboratories, workshops ...), which suggests how the complexity of the instruments attached to the telescopes has increased as time went on. It underscores the growing importance of these functions for achieving results in research.

On the other hand, as shown in Figure , the fraction occupied by service spaces decreases with time. This can be due to a greater specialization of these buildings as they suppress some functions not strictly astronomical (rest areas, cafeterias, exhibition halls or visitors centers) which may be relegated to secondary structures.

An important reason for the scattered or segregated pattern of modern observatories is the addition of new projects after the initial foundation. Good sites are rare, getting water, electricity and roads to that site may be the single most expensive part of the project. By comparison, adding further instruments to an established site is often a minor expense. Also, in terms of project planning, existing observatories have a science-grade record of the local weather, eliminating the need for possibly several years of site monitoring to establish whether a new proposed site is good enough.

# The Forthcoming Observatory

Since its inception, the astronomical observatory has reflected a continuous architectural adaptation of the building to the needs of astronomy. Thus, aspects such as structural stability, thermal stability, orientation of the construction, the specific configuration of windows, the adaptation to the different astronomical instruments, the generation of service and comfort spaces for the users of the observatory, the incorporation of new functions derived from the advances of science, the change of scale according to the growth of the instruments, the adaptation to non-urban geographic environments or the conciliation with new technologies have been addressed.

Obviously the current observatories, of the beginning of the 21st century, are vastly different from their modern or Islamic-medieval ancestors. Nor are they the same as those of the mid-twentieth century, even though the influence of the Hale Telescope has spread through all the astronomical centers that succeeded it. In any case, this evolution is not taking place in an intentional or planned way, but the changes take place in an organic and almost improvised way. This, which is not a problem in itself, does denote a lack of deep understanding regarding observatories, understood as buildings in which astronomical research and human activity coexist and relate through architecture. All this means that occasionally some of the different functions of the observatory come into conflict with each other, indicating a certain lack of criticism and reflection. It is difficult to answer clearly and unequivocally where astronomical observatories are heading, since there does not seem to be a roadmap that points in a clear direction, due to the fact that technological advances are unpredictable.

# Conclusions

Although the selection of astronomical observatories analyzed in this work could include a larger sample, it would not significantly change the result since the architecture and configuration of the buildings was frequently replicated with few variations.

In any case, our investigation has found evidence for two opposite directions of evolution:

-   continuation of the dissociation of astronomical instruments from other functions related to research;

-   return to heterogeneous uses of space in the same building such that the needs of its users can match the astronomical activity.

The first of these possibilities would involve a split of what is now known as an observatory. The result is a pair of buildings, one of which houses all the spaces and dependencies necessary for the staff to do their research (laboratories, offices, etc.), all the elements that could result in better comfort for workers (cafeteria, dining room, rest areas, etc.). In addition, this property should also house all kinds of technical rooms (control rooms, workshops ...) or accessories (facilities, warehouses ...) that would serve the astronomical function. The second building, separate and independent, would be completely unrelated to human activity (for example, the electricity supply).

The second possibility would involve unifying in a single building both the strictly intrinsic functions of astronomical research and the rest of the uses and services associated with the workers of the center or even visitors. For this, it would be necessary to carefully consider and correct the possible interference and incompatibilities between the different uses. This would result in a rather complex single construction.

When planning a new observatory in an area where other research centers are concentrated, it would be reasonable to choose a structure for the new telescope that would be separated from a place for workers and other users that could even be shared with other facilities and neighbors. However, opting for a unitary building would be adequate if the new observatory considers aspects beyond pure research (such as teaching or public outreach) or if the project were to be implemented in a remote and isolated place, saving resources and optimizing supply lines.

Overall it seems reasonable that the two opposite directions of evolution can coexist in the near future where they could represent the most appropriate solutions in different scenarios; see also .

We have profited from valuable discussions with Ronan Cunniffe and Gemma Olmedo Paradas. We are grateful to Luis Alberto Rodríguez and Antonio Cabrera-Lavers (Gran Telescopio Canarias), Cristina Rodríguez (Observatorio de Sierra Nevada), Santos Pedraz (Calar Alto Observatory) and Tatyana Sokolova (Special Astrophysical Observatory) for providing the requested material and pictures of the astronomical observatories they are related to. We also acknowledge very useful comments from the two anonymous referees. M. A. C. T. has been supported by a Junta de Andalucía predoctoral grant under project TIC-2839. Finally, we also acknowledge financial support from the State Agency for Research of the Spanish MCIU through the "Center of Excellence Severo Ochoa" award for the Instituto de Astrofísica de Andalucía (SEV-2017-0709).