key: cord-0692309-fa7jy8mp authors: Jain, Apurva; Rayal, Swapnil title: Managing medical equipment capacity with early spread of infection in a region date: 2022-03-03 journal: Prod Oper Manag DOI: 10.1111/poms.13684 sha: eeee31afcfdc0f03e9478c86b4a6fe9fcffa20a2 doc_id: 692309 cord_uid: fa7jy8mp We develop a model for a regional decision‐maker to analyze the requirement of medical equipment capacity in the early stages of a spread of infections. We use the model to propose and evaluate ways to manage limited equipment capacity. Early‐stage infection growth is captured by a stochastic differential equation (SDE) and is part of a two‐period community spread and shutdown model. We use the running‐maximum process of a geometric Brownian motion to develop a performance metric, probability of breach, for a given capacity level. Decision‐maker estimates costs of economy versus health and the time till the availability of a cure; we develop a heuristic rule and an optimal formulation that use these estimates to determine the required medical equipment capacity. We connect the level of capacity to a menu of actions, including the level and timing of shutdown, shutdown effectiveness, and enforcement. Our results show how these actions can compensate for the limited medical equipment capacity in a region. We next address the sharing of medical equipment capacity across regions and its impact on the breach probability. In addition to traditional risk‐pooling, we identify a peak‐timing effect depending on when infections peak in different regions. We show that equipment sharing may not benefit the regions when capacity is tight. A coupled SDE model captures the messaging coordination and movement across regional borders. Numerical experiments on this model show that under certain conditions, such movement and coordination can synchronize the infection trajectories and bring the peaks closer, reducing the benefit of sharing capacity. During the early spread of COVID-19, there has been widespread reporting about shortages of medical equipment like ventilators, vital sign monitors, and intensive-care-unit beds (Halpern & Tan, 2020) . Institutions like the FDA and WHO have provided information about monitoring and tracking such shortages (US FDA, 2020) . The demand for such equipment is driven by the evolution of infections in different regions. The supply is severely constrained. The regional administrators and healthcare authorities must determine the level of capacity they need in their regions. Given the limited availability, they must also identify ways to use this capacity to improve health outcomes in their regions. Our goal in this paper is to develop an analytical model to guide such decisions. Accepted by Special Issue Pandemic Editors, after 4 revisions. The determination of the right capacity and its management to improve performance is a core concern of the production and operations management (POM) research. The literature, however, has paid limited attention to managing demands that are driven by the spread of a contagious infection. The growth pattern of infections, an early exponential rise followed by a decline or a plateau, determines how much capacity will be needed and when. The literature is also relatively silent about situations where decision-makers, facing considerable uncertainty about infection-growth patterns and other model inputs, exhibit a wide range in the decisions they make (KFF, 2021) . We aim to develop a model that can capture the essential new features of this context: uncertainty in how infections grow and drive demand, a wide range of decision-making choices in different regions, and the need to compensate for the limited capacity. We take the perspective of a regional decision-maker (she) to address the following questions. Based on her region's infection growth and her estimates of its impact on the region's health and economy, how should she determine the medical equipment capacity? Given that such capacity may not be fully available, what is the impact of compensating actions-from early shutdowns to consistent enforcement messaging and regional capacity sharing-on her region's health outcomes? Our approach to answering these questions is analytical, but the model provides an understanding of trade-offs and insights into the value of collaboration that is practically useful. On the theoretical side, our contributions include developing a model that (i) derives medical equipment demand from the stochastic spread of infections, (ii) accounts for the variety in decision-makers ' actions, and (iii) proposes and evaluates ideas for improving health outcomes with constrained supply. Our focus is on the demand side; we do not model the acquisition of medical equipment supply. Specifically, we analyze a stochastic differential equation (SDE) of infection growth in the early stages. Based on this growth, we develop a measure of the chance that the need for medical equipment over a period will exceed a given capacity, the breach probability. We then describe a formulation of a regional decision maker's beliefs about costs related to economy and health and the expectation of a cure. We show how her preferred estimates interact with the spread of infections to determine the infection level at which a region should shut down to reduce its infection rate. We also analyze the capacity needed at the shutdown time. We present a heuristic rule for this decision and an optimal stopping formulation. We show that while a lack of capacity can drive up the breach probability, actions like reusability, early or high-intensity shutdowns, and consistent messaging can reduce it. We then develop a coupled SDE model to analyze the impact of capacity coordination across two regions. The model allows for travel between the regions and for the coordination of messaging and shutdown policy. We analyze the impact of such coordination on the regions' breach probability. We develop the following results and insights: • Modeling uncertainty in infection rates in an SDE formulation induces the needed realism in the model while still maintaining tractability for analyzing outcomes like the breach probability. • A region can shut down to reduce infection rate; should it do so and when? The model clarifies the trade-off between incurring a cost to shutting down the economy and reducing the probability of a breach. A heuristic decision rule determines the level of infections at which to shut down. • At the shutdown level of infections, capacitydetermination is a function of the decision maker's two estimates: the relative cost impacts of an economic shutdown versus a breach of health capacity and the time until the availability of a cure or a vaccine. A simple visual guide lets her determine the capacity she needs corresponding to her estimates. The guide maps the variety of decision-makers' beliefs and estimates to the range of preparatory actions we observe across various US states. • An optimal stopping formulation of the problem leads to observations similar to those based on the heuristic rule. • If capacity is limited, what actions can a region take to reduce its breach probability? The toolbox includes reusability, early shut down or high-intensity shut down, and consistent enforcement messaging. The model can determine the reduction in breach probability that these actions can achieve. Consistency in messaging reduces noise by limiting the variation in people's behavior; its impact on the breach probability is not always beneficial. • Looking beyond a region's boundaries, what is the impact of sharing capacity with a neighboring region? We model the dynamics as a coupled SDE and prove a technical result for the solution's existence and positivity. The model captures travel between regions and the correlation between the random processes underlying the two regions' infections. • Without any travel and correlation, capacity-sharing between the two regions will reduce the breach probability. Interestingly, this holds only if the capacity is not too tight. We identify two separate drivers for this effect: a pooling effect and a peak-timing effect. • Movement across the border of the two regions can increase the breach probability with shared capacity. • Consistency in pandemic-related messaging and joint efforts to maintain data accuracy correlates with the random processes underlying the two regions' infections. An increase in correlation can increase the breach probability with shared capacity. In addition to these practical insights, our theoretical contribution to the POM literature includes offering a tractable model to determine the equipment capacity as a function of the running maximum of infections evolving as a geometric Brownian motion. We believe that the paper contributes to the POM literature as an early effort to introduce uncertainty in the pandemic-driven demand model while keeping it tractable for capacity questions. The coupled SDE formulation for regional coordination is new to the literature as well; the proof of the condition for its solution's existence and positivity is a technical contribution to the literature. The insights into the drivers of capacity-sharing benefits introduce the idea of peak-timing, which is new to the literature. Section 2 describes how our model connects to three different streams of literature. Section 3 presents the model, and Section 4 offers an analysis of the shutdown level and the capacity. Section 5 presents the impact of compensatory actions. Section 6 presents the regional collaboration coupled SDE model. Section 7 discusses possible extensions and concludes the paper. We draw on several different streams of literature: epidemiology literature, stochastic finance literature, and capacity models with diffusion-driven demand in operations management. This section highlights how our model connects to these streams and how we differentiate ourselves from the existing models. Some of the recent and rapidly growing literature on pandemic modeling is also discussed. Our starting point is a model of how infections spread in a pandemic. There is a vast body of literature in epidemiology that develops and analyzes such models. The bestknown class of such models is compartment models, with a well-known example being the susceptible, infectious, recovered (SIR) model. These types of models divide the population into separate compartments and estimate transition rates between these compartments. If we drop the recovered compartment, it is called an susceptible, infectious, susceptible (SIS) model. If we add an exposed compartment, it is a susceptible, exposed, infectious, recovered model. If postinfection immunity fades away, the model is called susceptible, infectious, recovered, susceptible. Demographic and geographic details can be added to population transition rate inputs. The resulting ordinary differential equations can be solved numerically or simulated to project the number of infections over time. For a good introduction, see Martcheva (2015) . Our goal is to analyze capacity questions in a setting that captures the uncertainty inherent in the spread of infections. Most compartment models, however, do not offer analytical solutions even without any uncertainty. Our approach is to present a simplified, early-stage version of the SIR model before we add uncertainty to it. The early-stage approach, discussed in Britton (2009) among others, relies on the observation that at the beginning of the pandemic, the number of susceptible (almost the whole population) is extremely large, compared to the number of infections. The next section discusses the early-stage model in detail and offers other supporting references for this approach. We add uncertainty by explicitly modeling the probabilistic spread using an SDE. Some recent reporting (see, e.g., Roberts, 2020) has argued that classic epidemiological models tend to use constant reproduction rate numbers even though, in reality, the spread depends on uncertain factors. There is, however, plenty of discussion of stochastic modeling in epidemiology. For basic computational methods for analyzing probabilistic spread, see Vynnycky and White (2010) . The direct approach to introduce stochasticity in a compartment model is to use continuous-time Markov chains (CTMC). For the CTMC approach, now known as the stochastic general epidemic model, see Bailey (1975) . In general, these models yield only to computational analysis. The use of SDEs to model the spread is relatively new. The literature justifies the use of SDE models for epidemic modeling in several ways. Allen (2015) presents discrete-state CTMC models and then shows how they can be approximated by SDEs. In some cases, SDEs are introduced as a reasonable way to model growth without any need for further justification; see, for example, Gard (1988) . Another approach is to perturb the infection rate parameter by adding noise to it. This is the approach we take; the next section describes it in detail with supporting references. Our contribution is to combine the SDE method with the early-stage modeling discussed above to present an approach that is tractable enough to analyze questions about capacity, the area of interest in the production and operations management literature. In contrast to population-based compartment models, an alternative approach is to treat each individual as an agent and develop the epidemic dynamics based on interactions between agents. These individual-based models allow for the incorporation of agents' choices and reactions to testing and shutdown policies. See, for example, Acemoglu et al. (2020) and Drakopoulos and Randhawa (2021) . Given our focus on the aggregated need for capacity and on how it evolves over time, we do not take this individual agent-based approach. Our paper fits a rapidly growing literature on modeling the impact of economic, social, and operational policies in the context of the recent pandemic. Some examples of papers that focus on operations issues include Mehrotra et al. (2020), Alvarez et al. (2020) , and Yang et al. (2022) . The second stream of literature that we draw upon is stochastic finance. While we do not interpret our model using finance terminology, the analytical techniques we use are similar to those used in finance. Analysis and solution of SDEs is a foundation of stochastic finance. See, for example, Shreve (2004) . The pricing of exotic options requires an analysis of running maximum; we use a similar approach, see, for example, Buchen (2012) . Finally, we will model the decision to shut down as a stopping problem. Stopping problems are rigorously presented in Shiryaev (2007) and Oksendal (2003) , but they are also used in finance to make investment decisions, see, for example, Dixit and Pindyck (1994) . We extend this literature by presenting a new application of SDEs to model the spread of infections in two neighboring regions. The regions influence each other due to the possibility of people traveling across the border. They are also related through the correlation of Wiener processes that underlie the evolution of infections. This leads to the creation of a set of two coupled SDEs. Our contribution to this literature is a new result to show the existence and positivity of a solution to these equations. Finally, the third stream of literature we relate to is the analysis of capacity or inventory decisions in the presence of demands that follow the increasing-then-decreasing pattern discussed above. While we are not aware of any work that directly addresses demands influenced by the spread of infections, the first-increasing-then-decreasing pattern is similar to the new product diffusion. There has been some work in the operations management literature that relates to diffusion demand. Ho et al. (2002) and Kumar and Swaminathan (2003) compare the option of delaying the launch with the option to reject customers when capacity is limited. Shen et al. (2011) include pricing as a decision variable in addition to capacity and sales policy. These papers follow a new Production and Operations Management product demand model based on Bass (1969) diffusion, a deterministic approach. The use of SDEs and optimal stopping formulation for other POM problems is comparatively limited; see Kwon (2010) and Angelus (2020) for two recent examples. We contribute to this literature by offering a stochastic diffusion model to consider capacity questions in a new context, the spread of infections in a pandemic. We model a region with a decision-maker. The region is observing the beginning of the spread of infections. A continuous variable I t represents the number of infections in the region at time t. We will first present the dynamics of the spread and then discuss its underlying assumptions based on epidemiology literature. We model the growth of infections by focusing on the new infections generated by an infected individual. We define the transmission rate b to be equal to the new infections as the instantaneous fraction of the infected population. To consider the recovery of infected individuals, we model the recovery rate d as the instantaneous fraction of the infected population that leaves the infected compartment due to recovery or death. For notational brevity, we define the effective infection spread rate as r = b − d, a constant estimated based on actual observations. The rate of change in the number of infections is as follows: Our approach is a simplification of the SIR model as applied to the early spread stage. In the standard SIR model, the infection dynamics are presented as −dI t ∕dt = ( ) − I t where S t represents the number of susceptible at time t, N is the fixed total population, and , are rate parameters. For detailed explanations and parameter estimation, please see the references available in the previous section. The first term captures the increase in the number of infections as a fraction of interactions between the susceptible and infected populations. In the early stage spread of infections, the whole population N of the region, an extremely large number, is susceptible to infection. That is, there is an early stage in which S∕N ≈ 1 is a reasonable assumption. Substituting this in the above SIR equation leads to our model. In other words, when the number of susceptible (almost the whole population) is extremely large, compared to the number of infections, it is reasonable to assume that each infection generates a certain number of new infections irrespective of how many sus-ceptible remain. Therefore, it is sufficient to track only the number of infections to model early-stage spread. This earlystage approach is well supported in the literature. Employ this approach to model the spread of COVID-19 in stages. For a unit-size population, they assume S ≈ 1 in the first stage of the epidemic. Britton (2009) suggests that "The key reason for having an approximation during the early stages of an outbreak when there are many initial susceptible is that it is then very unlikely that any of the first numbers of infectious contacts happen to be with the same susceptible individual. Conversely, it is very likely that all of the first set of infectious contacts happen with distinct individuals." This formulation is equivalent to our model. Sazonov et al. (2011) model the spread in two stages where the first (early) stage is analyzed without the recovered compartment. The modeling advantage of this early-stage approach is that it allows us to introduce uncertainty in the model while keeping it tractable for capacity considerations. The disadvantage is that in the absence of an intervention, the number of infections in the model will continue to grow. We address this later by adding a second stage where a shutdown decreases the infection rate. The above dynamics are deterministic, but there is a great deal of variation in how infected individuals interact with others and spread the infection. At the beginning of the pandemic, the drivers of this variation, such as people's behavior, the accuracy of information, penetration of testing, and so forth, are poorly understood. The decision-maker perceives them as unpredictable variability in the spread. We model the infection spread rate as a stochastic processr where dr = rdt + s dW t , s is the volatility, and dW t is the increment of the standard Brownian motion (Wiener process) {W t } t≥0 . We can now rewrite the dynamics of the infection spread as follows: The stochastic process I t , which satisfies the above SDE is defined on an underlying probability space (Ω,  , P), which is equipped with a filtration F = { t } t≥0 and the stochastic process is adapted to F. Parameters r and s are constants. We have briefly reviewed the literature on stochastic epidemic models in the previous section. The parameter perturbation approach (addition of a time-dependent noise in the infection rate) that we have used to create the SDE model is directly supported by Gray et al. (2011 ), Cai et al. (2019 , Tornatore et al. (2005) , and references therein. It is useful to consider the interpretation of the volatility parameter s and the underlying Wiener process W t because these notations are not usually part of infection-spread models. The SDE is interpreted to mean that in successive small the numbers of new infections created by an infected individual are independent and identically distributed random variables that follow a normal distribution with mean rdt and variance s 2 dt. Thus, the Wiener process captures the timedependent variability of the number of contacts an infected individual will make in a small period, and s can be adjusted to scale its impact on the spread. This variability is driven by the inconsistency in the behavior of the infected individual and the behavior of others, a natural part of human interaction. A lack of standard messaging will add to this variability. Incomplete and non-uniform testing leads to noise in the observed number of infections and that too will add to this variability. A possible way to estimate s is to collect estimates of infection rates from various sources like hospitals, neighborhood health centers and use this sample to estimate it. For the interpretation and estimation of the infection rate, we refer to Fisher (2020) . In general, for the entire timeline defined below, we assume r − s 2 ∕2 ≥ 0 in order to avoid scenarios in which I t → 0 a.s. when t → ∞. We frame our model in terms of two periods, beginning at two discrete points in time. At timet 0 , when a very small number of initial infections I t 0 are observed, the regional decisionmaker starts tracking the number of infections. The infection rate is r 0 and the volatility is s 0 , both given constants. At a time t 1 of her choosing, she observes the number of infections I t 1 and decides if she will take action (e.g., mask mandates, shutdowns) to decrease the infections' spread. If she chooses to shut down, we change the infection rate to r 1 < r 0 and the volatility to s 1 , assumed to be equal to s 0 , s 1 = s 0 = s, unless otherwise mentioned. If she chooses not to shut down, the rate and volatility parameters remain unchanged. Even though it is a shutdown-or-not decision, for convenience, we refer to it as a shutdown decision. t 0 : A given constant, representing the start of the decision-maker's observations. I t 0 : A given constant, representing the very small number of infections observed at the time t 0 . t 1 : Decision time at which the decision-maker observes the number of infections I t 1 and determines whether to shut down. For now, it is assumed to be a given constant t 1 > t 0 . Later, we develop a decision rule to determine t 1 . I t orI t (I u , r, s): Stochastic process I t at time t > u, starting at time u with a known constant value I u , and following the early-stage SDE described above with known constant parameters r and s. For t 0 ≤ t < t 1 , we have I u = I t 0 , r = r 0 , s = s 0 . For t > t 1 , the process evolves conditional on the filtration up to the time t 1 and the decision taken at that time: If the decision at t 1 is to shut down, we have I u = I t 1 , r = r 1 < r 0 , s = s 1 ; and if the decision is to not shut down, we have I u = I t 1 , F I G U R E 1 Spread of infection in a region with and without shutdown during different time periods [Color figure can be viewed at wileyonlinelibrary.com] r = r 0 , s = s 0 . Depending on the context, we will drop the arguments and simply refer to the process as I t . We refer to the interval [t 0 , t 1 ) as the early community spread period. During this period, some medical understanding of the pandemic is starting to build, but there is no clear indication about if and when a therapeutic cure or vaccine may be discovered. There is no common understanding of when a medical solution may be available to halt the spread of infections. Nevertheless, a regional decision-maker must look ahead toward such a possibility when making a shutdown decision. That is, at the decision time t 1 , she needs to have an estimate of the horizon for her decision problem. It may be possible to argue that given the state of medical knowledge about the pandemic at the decision time, one can estimate a random variable to represent the horizon she should consider. A more likely scenario is that, given the extreme uncertainty, decision-makers let their preferences guide the level of optimism they feel about the pandemics' end. See, for example, Breuninger (2020) and Goldberg (2020) for differences in decision-makers' estimates of how long the pandemic will last. Therefore, we chose to model the decision horizon T as a preferred estimate of the decision-maker, which would vary from one decision-maker to the other. T : The decision-horizon, a constant; the decision-maker who is considering a shutdown decision at a time t 1 > t 0 will consider the impact of her decisions in [t 1 , If a shutdown is implemented at a time t 1 , we refer to [t 1 , t 1 + T] as the shutdown period, but we emphasize that t 1 + T does not model the actual reopening time. It only represents the decision maker's estimate of how far she prefers to look ahead. Given our focus on the capacity question during the early spread, we do not specifically model the reopening decision. In the extensions section, we provide a brief discussion of the reopening decision. Figure 1 shows two sample trajectories with and without shutdown. We believe that this two-period framing reflects how the recent pandemic unfolded in various regions, U.S. states, and other countries. At the same time, it is flexible enough to Production and Operations Management capture the variety of choices made by the regional decisionmakers about the timing and intensity of shutdowns. Our focus in this paper is on the decision time t 1 and the shutdown period. This is the period in which the region experiences a peak in infections and in which the questions about medical equipment capacity have the most relevance. As the decision-maker watches the infections increase in the early community spread period, it becomes clear that the disease will continue to grow and impose adverse health consequences. The degree of this adverse impact will depend on the region's level of preparation (medical equipment capacity). In the absence of therapeutics or vaccines, the decision-maker must decide if she should take any preventive measures and, if so, when. These measures can range from a ban on large gatherings to a wholesale shutdown of businesses. Different measures may bring different levels of improvements in the infection rates, but they will all impose an adverse economic consequence on the region. In our basic model, we capture these adverse health and economic consequences by the following parameters. m: Preparation level or medical equipment capacity, a known constant, expressed in terms of the maximum number of infections that the regional health system can simultaneously manage. Note that there are two types of equipment, those that cannot be reused and those that can be. The former's required capacity depends on the cumulative number of infections and the latter on the maximum number of infections. The required capacity for equipment like ventilators is of the latter type; this is how we define m here. PrB(I t 1 , m; r i, s i , T): Breach probability; the probability that starting with I t 1 number of infections at the time t 1 and evolving with the rate r i and volatility s i , i = 0, 1, the maximum number of infections in [t 1 , t 1 + T] will exceed the capacity m. C h : A fixed health cost incurred if the maximum number of infections exceeds the capacity m. C e : A fixed economy cost that the decision-maker believes will be incurred if she chooses to implement a regional shutdown. The primary performance measure in our model is the probability of a breach. If a breach occurs, the region incurs a fixed health cost. The practical relevance of these measures is evident from reporting such as Sanger-Katz et al. (2020) that discusses the risk of shortage and efforts to prepare for a breach. The region can influence this probability by incurring a fixed economy cost. Once again, the extreme uncertainty surrounding these events makes it unlikely that regional decision-makers have precise estimates for such costs. We believe that a model built on the decision-makers' preferred estimates or beliefs regarding such costs better suits the reality. See, for example, Leachman and Sullivan (2020) to support our modeling of differences in decision-makers' estimates of these costs. Our performance measures capture this situation's essential trade-off while directly connecting it to the paper's main new feature, infection rate uncertainty. In the extensions section, we briefly discuss other, more detailed metrics to capture the impact on health and economic performance. Decision-makers in different regions have different preferences about how they view the health versus economy debate and how optimistic they are about the availability of a cure. These preferences guide their estimates of C h , C e , and T. Given these preferences, a decision-maker can compare the costs of shutdown versus not shutting down as we analyze below. In this section, we analyze the expression for the number of infections, develop a formula for the breach probability, and then analyze the decision to shut down and how it relates to the capacity level. The first step in the analysis is to obtain an expression for the number of infections at time t, I t . We use standard methods (Shreve, 2004) from the analysis of SDEs to analyze the infection dynamics to obtain the following solution. Proposition 1. (3) It is straightforward to recognize the dynamics of infections as a geometric Brownian motion, leading to the above result. Proofs are available in the Appendix in Supporting Information. The equations describe how the number of infections evolves, first with community spread in the period t 0 ≤ t < t 1 and then with shutdown interventions in the period t ≥ t 1 . We next develop an expression for the probability that the number of infections exceeds the capacity, a measure we have defined as the breach probability. The explicit expression of I t in the above result allows us to determine the probability that the maximum number of simultaneous infections will exceed a given level of capacity over a time horizon. We evaluate this probability at the decision time t 1 , conditional on the filtration up to the time t 1 . the breach probability is trivially 1. For 0 < I t 1 < m, we have the following result: where ∅(⋅) is the standard normal cumulative cistribution function. We briefly discuss the plan for the rest of this section. As the time progresses forward from t 0 , we arrive at the given decision time t 1 > t 0 and observe I t 1 . We first propose a heuristic shutdown rule to determine if the decision-maker should shut down at time t 1 . Next, rather than using a given time t 1 , we use the rule to determine t 1 when the region should shut down and, at this time, what should be the region's shutdown level capacity. We then extend the basic model and present an optimal stopping formulation of the problem. The decision-maker has her preferred estimates C h , C e , and T, knows the current rate and volatility parameters r 0, s 0 , and estimates that a shutdown will change them to r 1 < r 0, s 1 . At time t 1 with I t 1 infections, should she shut down or not? We note that she is making this decision in a context characterized by a great deal of uncertainty about the pandemic's evolution and difficulty estimating precise cost parameters. Therefore, our first approach is to present a simple heuristic rule that can capture the major drivers of this decision and analyze their relationships and trade-offs. The rule is as follows. At a given time t 1 , with I t 1 infections, a decision-maker will choose to remain open if for all i ≤ I t 1 that is if, ΔPrB(i, m; r 0, s 0 , r 1, s 1 , T) = PrB(i, m; r 0, s 0 , T) − PrB(i, m; r 1, s 1 , T) ≤ C e ∕C h and will choose to shut down otherwise. The heuristic is directly motivated by the underlying tradeoff. If a regional decision-maker chooses not to shut down, the higher infection rate r 0 will continue, and the region will incur the health cost C h with a probability that, given her decision horizon, she will estimate as PrB(I t 1 , m; r 0, s 0 , T). If she chooses to shut down, the infection rate will decrease to r 1 . The region will incur the economy cost C e and will also incur the health cost C h with probability PrB(I t 1 , m; r 1, s 1 , T). The heuristic decision rule suggests that the decision-maker will remain open as long as the shutdown has a higher cost. The rule relies on the assumption that the region starts its observation at a very small number of infections, I t 0 . As it observes the number of infections grow from a small number, the rule recommends a shutdown at the first time when the estimated economic cost of shutdown exceeds the health cost of remaining open. We will highlight the limitations of this rule in the next section, but its advantage is clear. It is a simple connection between the decision-maker's economy and health beliefs, and it roots the shutdown decision in the dynamics of how infections evolve. Next, to make sure that the rule is easily implementable, it will be useful to show that the decision can simply be based on the observed number of infections at the decision time. The following result shows that there exist thresholds on different input parameters (costs, decision horizon, capacity level) such that if they are crossed, the region will not shut down. Otherwise, the rule is well-formed in the sense that there is a unique threshold on the number of infections that would determine the decision. Proposition 3. At the decision time t 1 , there exist thresh-oldsĈ e ∕C h (m, T, I t 1 ), andT(C e ∕C h , m, I t 1 ) such that if either C e ∕C h >Ĉ e ∕C h , or T 0 : I t ∉ D}. The proof of Proposition 4 holds for any initial infection level and defines the threshold I o as the supremum of values where smooth pasting condition is valid. If the infections exceed I o , the cost of infections will exceed the cost of imposing a shutdown. The role of the assumption c > threshold I o such that if I t > I o , it would always be optimal to shut down. This is so because the continuing infection cost increases in the number of infections, but in contrast, the terminal economy and health costs are bounded and also the underlying geometric Brownian motion has a positive drift. Though there would be no guarantee that I o < m, making the result less insightful in a practical setting. It is useful to consider the difference between the decision rule and the optimal stopping approaches. Without a perinfection cost, the decision rule looks for the first time that the trade-off between fixed economy and health costs justifies shutdown. With the per-infection cost, the optimal stopping formulation shifts focus to how that cost determines a level at which the influence of the per-infection cost justifies a shutdown. The qualitative behaviors of the optimal stopping problem and the decision rule strategy are fundamentally different. For the decision rule strategy, without the perinfection cost, it is possible that it is never optimal to shut down. But with the optimal stopping problem, with a nonzero per-infection cost, the shutdown should eventually occur with Probability 1. This is because the total per-unit infection cost will continue to increase and exceed the combined impact of health and economy costs as long as the latter costs are bounded. Starting with a decision maker's preferred estimates of costs of health and economy (C h C e ), and the level of her optimism about the availability of a cure or vaccine (T), the previous section analyzed the shutdown level of capacity. In a capacity-constrained capacity environment, the desired level of capacity may not be available. The objective would then be to reduce the probability of breach. Therefore, in this section, we analyze other actions a decision-maker can take to reduce the breach probability if the medical equipment capacity is constrained to a given value. We start with a simple observation regarding the impact of reusability of the capacity on the breach probability. See, for example, Srinivasan et al. (2020) for ventilator sharing. Recall that the capacity m was introduced in terms of the maximum number of infections that can be simultaneously treated. For example, let V represent the ventilator capacity and let the fraction ∝ (Gelles & Petras, 2020) represent the fraction of infected persons who will need the equipment. Then, m = V∕ ∝. Any reuse decreases ∝ and, for a fixed V, increases m. As a result, the breach probability decreases. Next, we discuss the impact of r 1 , the infection spread rate in the shutdown period. The decision-maker may have a range of mandated behaviors as possible options, from stayat-home orders to mask orders, which translate into different degrees of reduction in the infection rate post shutdown. For example, Flaxman (2020) lists recommended strategies ranging from self-isolation, ban on public events, social distancing, school closures, and lockdown with the corresponding percentage decrease in transmissibility from 5% to 45%. Our model can predict the impact of such reduction on the breach probability. Part (a) of the following result formalizes this observation. a. The previous section proposed a heuristic rule and an optimal stopping formulation to determine the shutdown level, but if the region does not have the capacity it needs, it can choose to deviate from the prescribed shutdown level. To formalize this argument, we first look at the impact of the shutdown level on the breach probability. Part (b) of the result shows that the breach probability increases as the shutdown level infections increase. These observations allow us to point out substitutability between the capacity and the shutdown infection rate and shutdown level. As Figure 4 shows, a lack of capacity increases the breach probability, but it can be compensated by taking action to increase shutdown intensity. The figure also shows that any decrease in the breach probability due to lack of capacity can be compensated by starting the shutdown early at a smaller number of infections. These strategies can be observed in real contexts as well. For example, by some measures of available bed capacity (adjusted for size), states like Oregon and Washington are the ones with the tightest capacity, while Texas and Minnesota come out at the top (Jones, 2020) . Both Washington and Oregon were the earliest states to issue stay-at-home orders, and Texas was one of the last ones (see Mervosh et al., 2020) . Shutdown enforcement s 1 Recalling our discussion for the introduction of volatility parameter in an earlier section, we note that volatility exists due to variations in how people behave in their interactions with each other. Before the shutdown, people follow their normal behavior patterns; some interact more and others less, with each interaction a possible source of infection. After the shutdown, the variation in people's behavior may be influenced by the strength of shutdown messaging and enforcement of its rules. Improvements in testing capacity and accuracy may also reduce noise in the system. We interpret these as reductions in volatility. Its impact on the breach probability can be complex. The intuition suggests that reduction in variability should benefit the performance. As Figure 5 shows, this is indeed the case. But, in cases with long time horizons, it can also increase the breach probability. This section considered strategies that a region can implement on its own. In the next section, we consider strategies that require coordination across regions. A major theme in reporting about the early spread of the disease in the United States is the lack of coordination between different states and regions. In the absence of coordination, many states and regions had closed down their borders, charted their own communication strategies, and relied on their own medical equipment capacity (Jancowicz et al., 2020) . This description fits closely with the model we have developed in the previous section. We will refer to the previous model as the individual-region model. There are different ways we can conceive coordination between individual regions in order to improve outcomes. In the United States, a lack of coordination on acquiring and managing the ventilator capacity came in for early criticism. Later, several states changed their policies to make their beds and ventilators available to both neighboring and non-adjacent states (Mulvihill, 2020; Ronayne, 2020) . A natural question to ask is about the improvement in outcomes that can be achieved by combining two individual regions' capacity. When two regions agree to make their total medical equipment capacity available to treat patients from either region, we refer to it as medical equipment capacity-sharing. It appears intuitively clear that such a sharing of capacity should be better for the regions. Our model offers a direct way to check this intuition; does this always hold true? We address this question as the first step in this section. Our interest, however, is not just in showing the existence of such capacity-sharing improvement. We are also interested in identifying and investigating other features of regional coordination that can influence the benefit of capacity sharing. As we better understand the drivers of this value, we aim to show that other features of regional coordination, such as information and movement coordination, can have either a positive or negative influence on the benefit of capacity sharing. Later sub-sections address such features. Let us identify two regions by assigning a superscript index k to each relevant parameter where k = a, b, with each index value representing one region. The regions are assumed to be independent in the sense that the infections evolve independently in each region. The next section considers the Production and Operations Management case of interdependent regions. We will focus on the dynamics of infections in the shutdown periods t ≥ t 1 in each region, assuming a shutdown is implemented at t 1 . Recall that in the individual region model discussed earlier,Î k t 1 represents the maximum number of infections in [t 1 , t 1 + T] and the capacity m k represents the number of simultaneous infections that the existing medical equipment capacity can treat. The health cost is incurred if the maximum number of infections exceeds the capacity, an event we refer to as a breach. To keep the focus on the operational performance, we consider the impact of capacity sharing on the breach probability represented by Pr[Î a t 1 + I b t 1 > m a + m b ] in the capacity-sharing case and by Pr[Î k t 1 > m k ], k = a, b in the individual-region case. The result states that the maximum of the sum of the two regions' infections is less likely to exceed the combined capacity than the sum of individual regions' probabilities of exceeding their own capacities. The result confirms intuition and is also general in the sense that it does not depend on the relative sizes, infection rates, and capacities of the two regions. There is less risk of breaching capacity from the combined region's perspective by combining medical equipment capacity. However, if we focus on the health cost as our metric, the benefit is less clear. It would depend on how we estimate each region's health cost in the capacity-sharing case. From an individual region's perspective, the benefit is less clear as it depends on the other region's size and capacity. To investigate this further, we consider two identical and independent regions. That is, they shut down at the same time t 1 with an equal number of infections I k t 1 and equal preparation levels m k , and that they have the same infection rate and volatility after shutdown. This allows us to focus only on the impact of sharing capacity, unmediated by any other differences. From each region's perspective, do we expect to see the breach probability always decrease if the two regions share a capacity? Interestingly, the next result suggests not always. still be below the combined capacity. Both risk-pooling and peak-timing play a role in this phenomenon. We are now ready to consider other aspects of coordination across regions that can promote interdependence between the two regions' infection evolutions. Specifically, we analyze two different features of regional coordination. First, interregional coordination may allow movement across the physical boundaries. We label such physical coordination leakage and seek to understand its impact on capacity sharing. Second, coordination may allow regions to exercise uniformity in communication strategies. We label this practice information coordination. We propose a way to incorporate it into our model and ask how this may influence capacity sharing benefits. Later, we briefly discuss another feature, shutdown coordination, where regions may synchronize their shutdown policies. This section models the evolution of infections with these new features, and the next section considers the impact of such interdependent evolution on the benefit of capacity sharing. Unlike our discussion of capacity sharing in the previous section, where each region's infections evolved independently, the coordination features we describe in this section require us to model the interdependence between how regional infections evolve. To that end, we first propose a new dynamics of infection evolution that links the two regions. We introduce two new parameters: a correlation index , and a leakage parameter l. They are explained in detail right after the dynamics presented below. Assuming a shutdown at t 1 , for t ≥ t 1 : where W a t and W b t are correlated Brownian motions. That is, To model physical coordination between different regions such that the travel between the regions is permitted but the travel from outside to either region is not, we use the leakage parameter l. Physical coordination across the two regions allows some movement across the borders, thus enlarging the area that can be insulated from outside influences. Such a policy has its advantages in providing the residents with a larger area for movement, and thus lessening the pressure brought on by the shutdown orders. For example, such physical coordination existed between states like New York and New Jersey in the early days of the pandemic. As seen in the equations above, a fraction l of the current infected population in one region may leak to the other region. Once there, this infected fraction will contribute to the other regions'' spread. To model information coordination, we focus on the sources of noise in the infection rate. Each region has its idiosyncratic character that drives the variation in the range of people's social behavior. The lack of testing capacity and inconsistent data reporting also introduce uncertainty. The Wiener process captures this variation. With information coordination, regions implement coordinated messaging and uniform standards of data capture and reporting. Initially, the lack of such coordination was heavily criticized. Later, multi-state collaborations were launched to coordinate data accuracy; see, for example, Battaglia (2020) . The evidence of heterogeneity in health-related social behavior is available in Allcott (2020) . The connection between consistent messaging and uniformity in social behavior is discussed in Gadarian et al. (2021) . In this framing, coordination of information and messaging across two regions will induce consistency in people's social behavior. We model such coordination as an increase in the correlation between the two regions' Wiener processes. Since we are proposing a new system of SDEs to model the spread of infections due to interactions between two regions, it is important to derive the conditions under which there exists a unique positive solution to the coupled equations. We proceed to do so in the following result. We show the existence of a unique solution using the atmost linear growth and Lipschitz continuity condition. The proof of the positivity of the solution is more involved. We create a sequence of stopping times, which are the first hitting times of a sequence of boundaries very close to 0. We then show that as the boundaries tend toward 0, the stopping time tends to infinity. The detailed proof is presented in the Appendix in Supporting Information. The literature on SDEs to model the spread of infections is limited. We believe that with good parameter estimation, the above equations can be used to model the early spread of infections. Because the set of equations above does not yield an analytical solution, we present further insights based on the above system's simulation analysis. The analysis is carried out on a testbed of parameter ranges. The observations reported in the next section hold true for this testbed. For clarity, we present graphs for specific parameter combinations, which are mentioned with each graph. Before analyzing the impact of regional interdependence, driven by information and movement coordination, on capacity sharing, we first use the simulation to set a baseline by revisiting Section 6.1 where there was no movement coordination (zero leakage) and no information coordination (zero correlation). Figure 6a shows that as capacity increases, the breach probability for both the individual region case, Pr[Î a t 1 > m], and the shared capacity case, Pr[Î a t 1 + I b t 1 > 2m], decreases as we expect them to. There is, however, a unique threshold capacity level below in which sharing capacity does not reduce breach probability. Benefits of capacitysharing occur only if the capacity is above a threshold; this is in line with Proposition 7. Figure 6b digs deeper into this by adding the breach probability for the intermediate system, Pr[Î a t 1 +Î b t 1 > 2m]. When capacity is above the threshold, the graph can differentiate between the riskpooling effect and the peak timing effect as discussed below Proposition 7. Physical coordination across the two regions allows movement across the borders increasing the leakage l and inducing interdependence between two regions' evolutions. Our focus is on its impact on the breach probability. We note that unlike Figure 6 and Section 6.1, the following discussion computes, I a t 1 ,Î b t 1 , andÎ a t 1 + I b t 1 , based on simulation of the coupled SDEs presented in Section 6.2. While the existence of a capacity threshold mentioned in Proposition 7 was proved for the independent region case, it is notable that a similar notion holds in the case of interdependent regions. Figure 7 shows that capacity sharing improves breach probability for the high-capacity case but not for the low-capacity case. When capacity is high, the figure shows that the breach probability in both individual-region and capacity-sharing cases is relatively indifferent to increasing leakage, but when capacity is tight, increasing leakage increases breach probability in the capacity-sharing case. At any given time, the SDEs with leakage replace the individual regions' number of infections with a combination of both regions' infections. This synchronizes the two trajectories bringing their peaks closer. This peak-timing effect then leads to an increase in breach probability when capacity is shared. The sample trajectories in Figure 8 show a visual example of how an increase in leakage can bring the two regions' trajectories together. If we assume increased travel to be a proxy for increased leakage then this suggests inverse outcomes for an increase in travel, an observation similar to that made in some empirical work (see, e.g., Knittel & Ozaltun, 2020) . Information coordination across the two regions means coordination on messaging and communications transmitted to populations in two regions. We interpret this as an increase in correlation between the two Wiener processes underlying the spread of infections in the regions. We consider the impact of the increasing correlation on the breach probability. Once again, Figure 9 supports the existence of a capacity threshold even with interdependent regions; capacity-sharing reduces breach probability only with high capacity. The figure shows that increasing correlation does not have a measurable influence on the individual region's breach probability. As the SDEs show, combining two Wiener processes still The figure further shows that with high capacity, increasing correlation brings the individual-region and capacity-sharing breach probabilities closer, thus reducing the benefit derived by capacity-sharing. This too can be explained by the peak timing effect. The larger the correlation, the more likely it is that peaks occur close to each other, driving up the breach probability for capacity-sharing. We close this section with an observation related to the coordination of shutdown timing. In our analysis till now, we have assumed that the two regions' shutdown levels are the same and the shutdowns at the same time. However, given the random evolution of infection in the early spread period, the two shutdown times will not be the same. This would be a further reason for the peaks to occur at different times. Therefore, the peak timing effect will drive the breach probability down even further. Conversely, we expect coordination on shutdown timing to decrease the breach probability under capacity sharing. To summarize, this section shows that medical equipment capacity sharing across two regions can be beneficial but not when capacity is very tight. Allowing movement across the border, interestingly, can smooth the infection trajectories and offer further benefits. Any coordination that increases the correlation between the underlying random processes can bring peaks in sync and reduce the benefits. We note that these conclusions say nothing about other benefits regional coordination can bring. However, our results suggest that regions should assess the benefit of capacity sharing in light of other elements, such as leakage and correlation, included in the coordination agreement. We believe that a significant part of the paper's contribution stems from being the first to apply SDE-based modeling to pandemic-related capacity. The model can capture the uncertainty in infection spread and yet remain tractable enough to address capacity questions. Being early in this research stream, we have mainly focused on building a model that captures this situation's essential features and asks fundamental questions. However, other questions remain, raising the possibility of a variety of extensions. We discuss them in this section. Our focus on the early-stage dynamics led to the SDE model we used here. The main modeling goal is to produce a trajectory that can capture an early exponential increase and then follow it up with a variety of possibilities leading to a plateau, long-term decrease, or resurgence. One extension would be to use time-dependent rate and volatility parameters in the geometric Brownian motion. Depending on how these parameters change, it is possible to produce a variety of trajectories. Such a model can better fit the wide range of long-term infection-spread patterns we have observed across US states and other regions. Another extension would be to bring the infection-spread SDE closer to the epidemiological literature by modeling additional compartments. Specifically, we can include the susceptible compartment in the model where the sum of the numbers of susceptible people and infected people equals a fixed population size. We have focused on the breach probability as our primary performance measure. There are, however, other performance measures that directly depend on medical equipment capacity. The total number of infections over the horizon or the cumulative number of infections for which capacity is unavailable are important performance measures to compute. On the economy side, the region's size and the shutdown length will be important determinants of the total cost. An extension can include these considerations in the optimal stopping problem. Combining the previous two extensions, modeling longterm trajectories and additional performance metrics, would also allow us to include other decisions in the optimal stopping problem. The most important example is the reopening decision. The reopening decision will depend on how the infections are likely to spread after reopening and its implications for health and the economy. We limited the scope of our model to modeling the demandside impact of capacity. In our judgment, that is the most relevant question in a constrained capacity environment. However, this leaves open the question of how this capacity is acquired. The most straightforward extension would be to include a per-unit capacity cost in the model. The reality, however, is more complex. Given the severe shortages, each region is competing against others to acquire limited capacity. Thus, an extension with a competitive supply model would be more appropriate. The previous point also offers a reason to extend our analysis of regional coordination. While we have focused on the Production and Operations Management coordination of existing capacity, the coordination on acquiring new capacity can be a worthy extension. This extension would be akin to the development of cooperative game models for distributing coordination benefits in newsvendormodel setting; such models followed the earlier literature showing the benefits of pooling coordination. Other extensions for the coordination model would include consideration of non-identical regions and shutdown policies. Our theoretical model delivers conclusions that relate differences in decision-makers' preferences (about the relative costs of economy and health and expectation of a cure) to differences in their shutdown policies. Another conclusion we draw relates the value of capacity sharing to the tightness of capacity and the extent of movement and messaging coordination across neighboring regions. An extension that tests these theoretical conclusions against actual observations in the US states and other countries would be a useful exercise. We conclude by reciting some of the main themes in the paper. This paper is a first attempt to explicitly include infection rate uncertainty in the pandemic modeling while maintaining tractability to analyze the impact of limited capacity on the region's health performance. We develop an SDE model of the spread of infection in the early stages of a pandemic, analyze its running maximum process, and use it to compute the required medical equipment capacity. The performance measure of our choice is the breach probability. Our results relate regional decision-makers' preferred estimates of economy and health costs and time-till-cure to the shutdown policy, the level of infections at which the region should shut down, and the corresponding medical equipment capacity. Given the shortage of capacity, we analyze a list of options to compensate for it. These options include reusability of the equipment, shutdown intensity, earlier shut down, and reduction in the messaging noise. The model analyzes the impact of these actions on the breach probability. We then consider sharing capacity across regions as a way to influence the breach probability. We show that shared capacity may not benefit individual regions if the capacity is too tight. Movement across borders decreases the breach probability, and an increasing correlation between the random processes underlying two regions' infections increases the breach probability. As discussed above, we offer several ways to extend the model. Testing, voluntary social distancing and the spread of an infection Polarization and public health: Partisan differences in social distancing during the coronavirus pandemic Stochastic population and epidemic models: Persistence and extinction A simple planning problem for COVID-19 lockdown Distributed renewable power generation and implications for capacity investment and electricity prices. Production and Operations Management The mathematical theory of infectious diseases and its applications A new product growth model for consumer durables State auditor unveils new multi-state collaboration to improve COVID-19 data accuracy Cuomo says the 'ultimate resolution' of coronavirus will come with vaccine in 18 months Stochastic epidemic models: A survey An introduction to exotic option pricing A stochastic differential equation SIS epidemic model with two independent Brownian motions Investment under uncertainty Demand control of information products: Why perfect tests may not be worth waiting for R0, the messy metric that may soon shape our lives, explained Estimating the number of infections and the impact of non-pharmaceutical interventions on COVID-19 in 11 European countries. Imperial College COVID-19 Response Team Partisanship, health behavior, and policy attitudes in the early stages of the COVID-19 pandemic Introduction to stochastic differential equations How ventilators work and why COVID-19 patients need them to survive coronavirus It's going to disappear': Trump's changing tone on coronavirus A stochastic differential equation SIS epidemic model United States resource availability for COVID-19 Brownian models of performance and control Managing demand and sales dynamics in new product diffusion under supply constraint Why the U.S. still has a severe shortage of medical supplies Would the United States have had too few beds for universal emergency care in the event of a more widespread COVID-19 epidemic? State COVID-19 data and policy actions What does and does not correlate with COVID-19 death rates Diffusion of innovations under supply constraints Invest or exit? Optimal decisions in the face of a declining profit stream Some states much better prepared than others for recession A model of supply-chain decisions for resource sharing with an application to ventilator allocation to combat COVID-19 See which states and cities have told residents to stay at home U.S. states share, get creative in hunt for medical supplies Stochastic differential equations: An Introduction with applications How to think like an epidemiologist California lends 500 ventilators to 4 states, 2 territories These places could run out of hospital beds as coronavirus spreads A two-stage model for the SIR outbreak: Accounting for the discrete and stochastic nature of the epidemic at the initial contamination stage Stochastic orders New product diffusion decisions under supply constraints Optimal stopping rules Stochastic calculus for finance II A rapidly deployable individualized system for augmenting ventilator capacity Stability of a stochastic SIR system Medical device shortages during the COVID-19 public health emergency Stochastic calculus, filtering, and stochastic control lecture notes An Introduction to infectious disease modelling Design of Covid-19 testing queues U.S. state and local government responses to the COVID-19 pandemic