LLM-powered reasoning in agent-based modeling

arXiv cs.AI Papers

Summary

This paper introduces HALE, a scalable hybrid agent-based and language-driven epidemic modeling framework that leverages LLMs to predict human decision-making in ABM simulations, demonstrating improved accuracy in modeling COVID-19 spread in Salt Lake County.

arXiv:2607.06757v1 Announce Type: new Abstract: Agent-based modeling (ABM) has the capability to model millions of individuals and their interactions, which is useful for policy making. However, ABMs have traditionally relied on static prior, which prevents the models from adapting to real-time changes. Our research provides a novel approach to addressing this information gap. Large language models (LLMs) offer new opportunities to predict human decision-making. Here, we introduce a scalable Hybrid Agent-based and Language-driven Epidemic (HALE) modeling framework that leverages LLMs to predict human decision-making in an ABM simulation. As a proof-of-concept, we use HALE to simulate COVID-19 and its effects in Salt Lake County, UT.
Original Article
View Cached Full Text

Cached at: 07/09/26, 07:54 AM

# LLM-powered reasoning in agent-based modeling
Source: [https://arxiv.org/html/2607.06757](https://arxiv.org/html/2607.06757)
Dakotah Maguire1Adam Spannaus1Joe Tuccillo2Maksudul Alam3Sudip K\. Seal3John Gounley1and Heidi Hanson1

###### Abstract

Agent\-based modeling \(ABM\) has the capability to model millions of individual humans and their interactions, which is valuable for constructing a population\-scale digital twin to support policy decisions\. However, ABMs have traditionally relied on static and outdated data, which prevents the models from adapting to real\-time changes in a given scenario\. Complete real\-time information about the activities of a population is often dynamic, as mobility depends on epidemic conditions and human decision\-making\. The research presented here provides a novel approach to addressing this information gap\. Large language models \(LLMs\) offer new opportunities to predict human decision\-making as a scenario unfolds\. To this end, we introduce a scalable Hybrid Agent\-based and Language\-driven Epidemic \(HALE\) modeling framework that leverages LLMs to predict human decision\-making in an ABM simulation\. As a proof\-of\-concept, we use HALE to simulate COVID\-19 and its effects in Salt Lake County, UT, from September 2020 to February 2022\. The HALE framework better captures the observed epidemic peak and size, whereas simulations that use only ABM tend to overestimate the epidemic’s effects\.

Keywords:Large language models \(LLMs\), agent\-based modeling \(ABM\), activity\-based network, individual\-based network model, digital twin, COVID\-19\.

## 1Introduction

Agent\-based modeling \(ABM\) is a powerful tool for modeling complex dynamic systems at the individual level\. Understanding the impact of individual behavior on disease spread is critical for accurate disease forecasts\. ABM enables us to model individual human behavior \(e\.g\., mask\-wearing, vaccine acceptance\) based on a distribution of a population, their movement over time, and their social interactions\[[2](https://arxiv.org/html/2607.06757#bib.bib47)\]\. Typically, ABMs are initialized with these parameters, and the complex interactions of agents are modeled over a set time period\. However, human behavior is inherently uncertain and depends on various factors, including social environment, geographic location, climate, and weather\. Human preferences also change over time based on external environments and social pressures\. The reasoning abilities of large language models \(LLMs\) can assist in modeling change in human behavior during a disease outbreak\[[7](https://arxiv.org/html/2607.06757#bib.bib2),[39](https://arxiv.org/html/2607.06757#bib.bib3),[33](https://arxiv.org/html/2607.06757#bib.bib4),[15](https://arxiv.org/html/2607.06757#bib.bib5)\]\.

Although advancements in LLMs have improved knowledge across various domains\[[40](https://arxiv.org/html/2607.06757#bib.bib44)\], their application in epidemic modeling—especially individual\-based network modeling—has not yet been fully explored\.

LLMs can assist epidemic modeling and simulation in various ways, such as performing literature reviews, forecasting epidemics\[[16](https://arxiv.org/html/2607.06757#bib.bib42)\], and predicting human reactions to an epidemic\[[7](https://arxiv.org/html/2607.06757#bib.bib2)\]\. LLM agents can also exhibit realistic human behaviors\[[12](https://arxiv.org/html/2607.06757#bib.bib8)\]\. Interestingly, the collective intelligence of a generative agent\-based model can mimic real\-world responses to disease spreading\[[36](https://arxiv.org/html/2607.06757#bib.bib9)\]\. In this study, we investigate how LLMs can infer the change in human behavior during an epidemic and support ABMs in dynamically updating network structures\. These ABM\-based population\-scale digital twins use real\-world human activity data, which consist of different types of activities, including staying at home, going to school, being in the office, and shopping, along with their respective locations\. However, detailed activity data are often difficult to obtain, and this difficulty often requires researchers to generate network models from historical survey and demographic data\. Moreover, these models do not account for the change in an individual’s activities based on perceived threat or social pressures\. However, LLMs may be useful for modeling changes in behavior as the perceived threat level evolves over time\.

Network development is one of the pre\-processing steps in individual\-based network simulation\. In these networks, nodes can represent different entities \(e\.g\., people, locations, animals, plants\) depending on the research objective\. In epidemic models for diseases such as COVID\-19, nodes represent individuals, and links represent connections between two people\[[22](https://arxiv.org/html/2607.06757#bib.bib20)\]\. These links form when two people come into close proximity or visit the same location at the same time\. Information on individuals’ activities, including time and location, can be used to construct an activity\-based temporal network\. However, detailed activity data are often missing or only partially available\. Moreover, activity data may be available for a population during a specific time period, whereas inference is required for another period\. In an individual\-based temporal network model, the activities can be divided into two categories: deliberate and non\-deliberate actions\. LLMs can assist ABMs in modeling these deliberate actions when creditable data are not available\. For example, in an epidemic ABM simulation, spreading a virus from an infected person is a non\-deliberate action, whereas going out for outdoor activities is a deliberate action of an agent\. Deliberate actions indirectly influence non\-deliberate actions\. Additionally, agents in different locations may make different choices under similar conditions; they may perceive threats differently than their neighbors and react accordingly\. In this research, we address these data gaps by leveraging LLMs within ABM simulations\.

In this research, we have developed a spatiotemporal Hybrid Agent\-based and Language\-driven Epidemic \(HALE\) framework for near real–time forecasting\. LLM\-powered ABM simulation is a promising new research area\. However, most of the research on this topic is in other domains, is conceptual in nature, and has not yet reached the development phase\[[18](https://arxiv.org/html/2607.06757#bib.bib40)\]\. Researches have also replaced ABM agents with generative agents and modeled their behavior with generative AI\. Unfortunately, this approach is often computationally expensive and not intuitive or scalable\[[36](https://arxiv.org/html/2607.06757#bib.bib9)\]\. Therefore, we have developed a hybrid framework that leverages the strengths of both ABM and LLMs\. The HALE framework can simulate ABM models with millions of agents while enabling computationally efficient LLM inference\.

The major contributions of this research are as follows:

- •A scalable workflow to construct dynamic temporal networks models using a partially available real\-world activity dataset comprising millions of agents\.
- •A hybrid agent\-based and language\-driven epidemic modeling framework \(HALE\) to enable near real\-time forecasting\. Previous studies have primarily used ABM for epidemic analysis from retrospective data\. Although these methods are capable of forecasting, they do not adapt to real\-time changes in external environments and perceived threat level\. In reality, however, an individual’s mobility varies dynamically in response to an evolving epidemic situation\. In our HALE framework, we incorporate an LLM as a feedback loop to adjust human mobility in response to epidemic conditions\.
- •The performance analysis of the HALE framework by using an ablation study\. We also provide an analysis of the LLM’s inference across different spatial population groups\.

ABM is an important decision tool as it can model effects of different public health decisions \(e\.g\., lockdown, mandatory mask\-wearing\) at the individual level\. However, people often make different decisions depending on their surroundings and demographics features, which ABM cannot capture\. To address this challenge, we develop the HALE framework in which LLMs assist ABM model in individual\-level decision\-making\. The LLM supports the HALE framework to adjust mobility in response to an outbreak\. The research observes that different population groups react differently to an outbreak depending on their demographic characteristics\.

## 2Methods

### 2\.1Datasets

To develop the ABM, we use population data for Salt Lake County, UT, which represents a population of approximately 1\.1 million\. We then use Oak Ridge National Laboratory’s UrbanPop model to generate the synthetic population\[[31](https://arxiv.org/html/2607.06757#bib.bib19)\]\. Each individual in the synthetic population set has the following features assigned: age, race, sex, home ID, geographic location \(latitude, longitude\), living arrangement, occupation, and income level\. The statistical distribution of the synthetic population matches with the American Community Survey’s Public\-Use Microdata Sample \(PUMS\) at the census block group level\.

Each individual is assigned to a specific activity for a specific time by using data from the 2017 National Household Travel Survey \(NHTS\)\[[1](https://arxiv.org/html/2607.06757#bib.bib10)\]\. The Federal Highway Administration conducts the NHTS, which serves as an authoritative source of information on the travel behavior of the US public and enables analysis of trends in personal and household travel\. The activity schedule data include activities for each individual for every hour of the week, and an individual may engage in the same activity for multiple consecutive hours\. The activity schedule data contain information on time, activity type, home location, and activity location\. Locations are provided using Uber H3 hexagonal cells at a level\-8 resolution\[[6](https://arxiv.org/html/2607.06757#bib.bib11)\]\.

The UrbanPop synthetic population includes nighttime \(residential\) and daytime \(work or school\) locations at the census block group level\. Each individual is assigned a residential point location \(house\) from FEMA’s USA Structures dataset\[[38](https://arxiv.org/html/2607.06757#bib.bib23)\]by matching synthetic household dwelling characteristics with the properties of residential structures available within each block group\[[30](https://arxiv.org/html/2607.06757#bib.bib24)\]\.

Table[1](https://arxiv.org/html/2607.06757#S2.T1)presents a part of the activity schedule data for two persons\. In the synthetic population, each person has a unique identifier,P​I​DPID\. EachP​I​DPIDis linked to a household ID,H​I​DHID, which is a residence located in a hex H3 \(seeHome assignmentcolumn in Table[1](https://arxiv.org/html/2607.06757#S2.T1)\)\. AgentP​I​D01PID\_\{01\}in Table[1](https://arxiv.org/html/2607.06757#S2.T1)is at home on Monday from 0600 to 0700\. Therefore, the H3 hex IDs for the home location and the activity location are the same for this activity\. This individual \(agent\) goes to school after 0700, and the school location is in H3\_hex\_id 5\. Similarly,P​I​D02PID\_\{02\}has a work activity from 0700 to 0800, with the work location in H3\_hex\_id 8\. The activity schedule dataset contains approximately 189 million entries covering 168 hours \(7×247\\times 24\) for a population of 1\.1 million in Salt Lake County, UT\.

The activity schedule data include 22 distinct activities\. The activity listAAis \{home, school, work, exercise, services, retail, leisure, religious/community, meals, visiting friends and relatives, drop\-off/pickup, travel, changing transportation mode, errands, volunteering, medical, childcare, adult care, work home, work related, other, and unknown \(activity purpose not given by respondent\)\}\. The activities for home, school, and work are anchored activities in the schedule\. The activity schedule data contain activity location information only for the anchored activities \(last column of the Table[1](https://arxiv.org/html/2607.06757#S2.T1)\)\.

Table 1:Sample of activity schedule data at the individual level\.
### 2\.2Activity\-based temporal networks development

We develop a network,Ga​\(t\)G\_\{a\}\(t\), for each activity at each hour\. Therefore, for one week, we have 3,696Ga​\(t\)G\_\{a\}\(t\)networks that correspond to 22 activities\. Each network consists of 1,133,023 nodes \(agents\)\. The network development algorithms vary depending on the activity\.

#### 2\.2\.1Hourly networks for anchored activities

##### Home activity:

Thehomeactivity contains point locations for each residential place\. At hourtt, agents orP​I​DPIDs that share the sameH​I​DHIDand have the activityhomeorwork\-homeare connected in the networkGa​\(t\)G\_\{a\}\(t\)for thehomeactivity at hourtt\.

##### School activity:

Theschoolactivities also contain school locations\. From the activity schedule dataset, we first calculate the total number of students for each school\. Some schools have few students \(fewer than 19\), whereas others have many \(more than 500\)\. In the hourly school activity networkGa​\(t\)G\_\{a\}\(t\), students from smaller schools \(fewer than 19 students\) are fully connected if they participate in theschoolactivity at timett\.

We divide students from larger schools into classes and maintain the average public school class size based on data from the National Center for Education Statistics\[[13](https://arxiv.org/html/2607.06757#bib.bib12)\]\. We keep a record or map that links students to their respective classes to ensure that each student attends the same class every day\.

For theschoolactivity hourly network, we create two types of networks: the in\-class network and the out\-of\-class network\. In the in\-class network, students belonging to the same class and having theschoolactivity at timettform a fully connected network\. For the out\-of\-class network, we construct an Erdős–Rényi random network among students from the same school, with an average degree of 4 within the 8\-hour period\[[10](https://arxiv.org/html/2607.06757#bib.bib18)\]\. We assume that students have eight hours of school each day\. The out\-of\-class network represents contacts outside of the student’s own class \(e\.g\., gym, cafeteria, school bus\)\. Finally, we merge the in\-class and out\-of\-class networks to form the complete hourlyschoolactivity network\.

##### Office activity:

Theofficeactivities contain information about the office building location at the H3 level 8 index \(∼\\sim0\.74 km2\); the exact office is unknown\. Within each hex, some office buildings have many workers, whereas others have few, similar to schools\. However, the data do not include the exact office location\. They also do not include information about the layout, such as whether it is open\-plan or composed of cubicles\. Therefore, we do not divide offices into smaller units and utilize a configuration random network model to estimate contact\. In the hourlyofficeactivity network, we create a configuration random network in which the degree of each node is assigned from a Gaussian distribution:𝒩​\(μ,σ2\)\\mathcal\{N\}\(\\mu,\\sigma^\{2\}\)\. For the age range of 18–59, we consider hourlyμ=21\.1548\\mu=\\frac\{21\.154\}\{8\}andσ2=10\.588\\sigma^\{2\}=\\frac\{10\.58\}\{8\}\[[10](https://arxiv.org/html/2607.06757#bib.bib18),[22](https://arxiv.org/html/2607.06757#bib.bib20)\]\. We assume that the workday is eight hours\.

#### 2\.2\.2Hourly network for non\-anchored activities

Thenon\-anchored activities\(e\.g\., meals, retail, exercise\) do not have exact activity location information\. Therefore, we create contacts for them differently\. We use a spatial network model\[[21](https://arxiv.org/html/2607.06757#bib.bib22),[24](https://arxiv.org/html/2607.06757#bib.bib21)\], specifically an exponential distance dispersal kernel network, to construct hourly contact networks for eachnon\-anchored activity\. In this spatial network model, the probability of a connection between two nodes decays exponentially with distance\.

We select the spatial network model to capture the phenomenon in which people from nearby locations have a higher probability of meeting each other because they may go to the same restaurant, shop, or gym given their proximity\. In the exponential distance kernel, the probability,P​\(i,j\)P\(i,j\), forms a link between nodeiiandjj:

P​\(i,j\)=λ​e−λ​di​jP\(i,j\)=\\lambda e^\{\-\\lambda d\_\{ij\}\}\(1\)Here,λ\\lambdais the average travel distance for an activity, anddi​jd\_\{ij\}is the distance between the home locations ofiiandjj\. Average travel distanceλ\\lambdavaries across different activities\. In our spatial network model, we selectλ\\lambdavalues for the United States based on various previous surveys, research studies, and reports\. For example, the average travel distance for exercise is 4 miles\[[4](https://arxiv.org/html/2607.06757#bib.bib13)\]; average commuting time to church is around 15 minutes and roughly 11\.25 miles; the average distance between food establishments and homes is 2\.6 miles\[[19](https://arxiv.org/html/2607.06757#bib.bib15)\]; average travel distance for childcare is approximately 5 miles\[[5](https://arxiv.org/html/2607.06757#bib.bib16)\]; and a 2017 study reports the average distance for medical/dental care travel in the United States is 10\.04 miles\[[3](https://arxiv.org/html/2607.06757#bib.bib17)\]\. For activities for which average travel distances are unknown, we use the mean value of the knownλ\\lambdas\.

#### 2\.2\.3Hourly networks for all activities

For each hour, we develop 22 networks that correspond to 22 different activities\. To create the hourly network, we merge all the individual activity networks\. The hourly network contains no duplicate links because the activity schedule data ensure that a person or agent cannot be in two locations simultaneously\. The edge list of the hourly network records the activity associated with each edge as an attribute\. The average degree of the hourly networks differs between daytime hours and nighttime hours\.

#### 2\.2\.4Daily networks

We merge the hourly networks to create daily networks\. In the daily networks, we add an additional edge attribute calledweight\. If an edge appears inhhhourly networks for the same activity, then its weight in the daily network is set tohh\. Table[2](https://arxiv.org/html/2607.06757#S2.T2)shows a portion of the adjacency list of a daily network\. The first row in Table[2](https://arxiv.org/html/2607.06757#S2.T2)indicates thatP​I​D01PID\_\{01\}andP​I​D05PID\_\{05\}are connected in the daily network, with a weight corresponding to five hours of office activity in the same location\.

These daily networks are used for the ABM simulation\. The average number of connections per agent in the daily networks is approximately22\.4122\.41, which is consistent with previous literature\[[10](https://arxiv.org/html/2607.06757#bib.bib18)\]\.

Table 2:A portion of the adjacency list of a daily network\.

### 2\.3The HALE Framework

#### 2\.3\.1Architecture

HALE is an individual\-based framework that couples dynamic, AI\-generated mobility patterns with ABM\. In HALE,NNagents \(or individuals\) are connected via an activity\-based temporal network,G​\(t\)G\(t\)\. Each agent,nin\_\{i\}, is associated with a property vector,PP\. The properties can be demographic information, including sex, age, home location, school, workplace, average travel distance, and education\.

Figure[1](https://arxiv.org/html/2607.06757#S2.F1)presents the architecture of the HALE framework\. The ABM component of the framework models epidemic spreading by using a disease model such as susceptible\-infected\-susceptible \(SIS\), susceptible\-infected\-recovered \(SIR\), or susceptible\-exposed\-infected\-recovered \(SEIR\)\. A health status vector,H∈𝐑NH\\in\\mathbf\{R\}^\{N\}, tracks the health status of all agents\. ABM updates the health states of agents at discrete time steps,Δ​tA​B​M\\Delta t\_\{ABM\}\. For the SIR disease model, an agent’s health state ishi∈\{S,I,R\}h\_\{i\}\\in\\\{S,I,R\\\}\. At each time step,Δ​tA​B​M\\Delta t\_\{ABM\}, infected agents can spread the disease to their healthy neighbors in the networkG​\(t\)G\(t\), with a transmissibility,β\\beta\. Then, the state of a healthy agent changes to infected,hi:S→Ih\_\{i\}:S\\rightarrow I\. These events depend on the network or connectivity\. Infected agents recover,hi:I→Rh\_\{i\}:I\\rightarrow R, with a probability,γ\\gamma\. This event occurs independently of the network structure\.

![Refer to caption](https://arxiv.org/html/2607.06757v1/x1.png)Figure 1:The HALE framework—scalable agent\-based simulation with LLM reasoning—is applied to an individual\-level network model\. The LLM updates specific deliberate actions of agents\. In this model, agents are connected via temporal, activity\-based links\. Solid edges represent connections formed by fixed actions\. These connections are temporal but do not depend on LLM reasoning—for example, connections formed by living in the same home\. Dotted edges indicate connections resulting from deliberate actions that evolve with LLM reasoning\. The colors of individuals—green, orange, and blue—represent different health states\. ABM updates the health state of agents after each ABM time step,Δ​tA​B​M\\Delta t\_\{ABM\}\. After eachk​Δ​tA​B​Mk\\Delta t\_\{ABM\}\(orΔ​tL​L​M\\Delta t\_\{LLM\}\) time step, network updates the dotted edges from LLM reasoning\.For the LLM component, the agents are divided intoLLsubpopulations or groups based on their demographic properties and residential locations\. Here,L=\{l1,l2,l3,…​…\}L=\\\{l\_\{1\},l\_\{2\},l\_\{3\},\.\.\.\.\.\.\\\},lj=\{ni,ni\+1​…\}l\_\{j\}=\\\{n\_\{i\},n\_\{i\+1\}\.\.\.\\\}, and\|N\|≫\|L\|\|N\|\\gg\|L\|\. Att=0t=0, we initializeLLLLM agents – one for each group\. The groups or LLM agents have different sizes,\|li\|≠\|lj\|\|l\_\{i\}\|\\neq\|l\_\{j\}\|\. LLM agents are distinct and independent\. The outputs of the LLM agents are used in ABM simulation\. Therefore, we use structured LLM output formats\[[14](https://arxiv.org/html/2607.06757#bib.bib30)\]to make the LLM inference machine\-readable\.

The ABM simulation then models the disease dynamics ofNNagents on the networkG​\(t\)G\(t\)as time evolves\. The LLM agents generate reasoning for mobility patterns of theLLLLM agents based on the real and perceived environment over time\. Then, the HALE framework will update the network,G​\(t\)G\(t\), based on the LLM reasoning\. For example, if an LLM agentlil\_\{i\}decides to decline an outdoor activity and stay at home, then links formed from that activity for the subpopulation for that LLM agentlil\_\{i\}inG​\(t\)G\(t\)will be deactivated and will not spread the disease for that time\. ABM simulation updates the states of agents at each ABM time step,Δ​tA​B​M\\Delta t\_\{ABM\}, while the LLM assesses the current state of the system and provides inference on mobility at each LLM time step,Δ​tL​L​M\\Delta t\_\{LLM\}\. After eachΔ​tL​L​M\\Delta t\_\{LLM\}time step, ABM simulation waits for the LLM inference\. Because the ABM requires more frequent updates compared to the LLM component in the HALE framework, we haveΔ​tL​L​M=k​Δ​tA​B​M\\Delta t\_\{LLM\}=k\\Delta t\_\{ABM\}, wherekkis an integer greater than 1\.

#### 2\.3\.2Spatial LLM grouping

The ABM simulation contains millions of agents\. To add LLM reasoning to an ABM simulation while maintaining a reasonable compute time and accounting for Tobler’s First Law of Geography, we group the agents based on their spatial location and demographic properties\. We consider 23 municipalities in Salt Lake County, UT: Bluffdale, Cottonwood Heights, Draper, Herriman, Holladay, Kearns, Magna City, Midvale, Millcreek, Murray, Riverton, Salt Lake City, Sandy, South Jordan, South Salt Lake, Taylorsville, West Jordan, West Valley City, Alta, Brighton, Copperton, Emigration Canyon, and White City \(Figure[2](https://arxiv.org/html/2607.06757#S2.F2)\)\.

![Refer to caption](https://arxiv.org/html/2607.06757v1/x2.png)Figure 2:Spatial LLM grouping for Salt Lake County, UT\. We use race, sex, and age to further subdivide the people in cities with larger populations\.We divide the municipalities into two categories based on population size: larger cities and smaller cities\. Alta, Brighton, Copperton, Emigration Canyon, and White City are classified as smaller cities, whereas the remaining municipalities are classified as larger cities\.

We further divide the population of the larger cities by sex \{male, female\}, race \{White, Black or African American, Asian, Native American, Pacific Islander, Multiethnic, Other\}, and age group \{0–18, 19–24, 25–34, 35–44, 45–59, over 60\}\. However, we do not further divide the population of the smaller cities because this would yield very small, often empty, groups\. In total, there are 1,552 population groups, each represented by a corresponding LLM agent\.

#### 2\.3\.3ABM\-LLM structure

In this hybrid setup, the ABM and LLM run sequentially, passing their results back and forth to each other\. Here, ABM simulation and the LLM have different time scales: hourly time stepsΔ​tA​B​M\\Delta t\_\{ABM\}for ABM simulation and weekly time stepsΔ​tL​L​M\\Delta t\_\{LLM\}for the LLM\. These time steps are flexible and can be different for different research objectives\.

##### ABM simulation:

ABM will simulate the disease progression in the daily weighted network,G​\(t\)G\(t\), where,t∈\{d​a​y​0,d​a​y​1,d​a​y​2,…\}t\\in\\\{day0,day1,day2,\.\.\.\\\}\(Algorithm 1 in\[[28](https://arxiv.org/html/2607.06757#bib.bib46)\]\)\. For our experiment, we choose the SIR disease model and use incidence data for COVID\-19\.

We design the ABM simulation to model disease spreading using the Repast modeling framework\[[25](https://arxiv.org/html/2607.06757#bib.bib26),[28](https://arxiv.org/html/2607.06757#bib.bib46)\]\. Repast is a highly scalable distributed framework that can leverage multiple processing cores via MPI \(Message Passing Interface\)\.

##### LLM reasoning:

In this research, we use LLMs to explore an important question: whether a person would continue their normal daily activities \(e\.g\., school, work, shopping\) if they were aware that a disease is spreading in their community\.

Here, the LLM takes the incidence level from ABM and provides reasoning for different LLM groups\. After eachΔ​tL​L​M\\Delta t\_\{LLM\}, ABM sends the simulation or disease scenario to the LLM, which then provides the reasoning for each LLM agents\.

We use zero\-shot learning and assume that the LLM already have the baseline knowledge needed to infer typical human activities in a given area\[[8](https://arxiv.org/html/2607.06757#bib.bib28)\]\. For our experiments, we use pretrained and instruction\-tuned generative models with 8 billion parameters from Llama 3\.1\[[11](https://arxiv.org/html/2607.06757#bib.bib14)\]\. We also utilize vLLM, an open\-source library for fast inference of LLMs\[[17](https://arxiv.org/html/2607.06757#bib.bib29)\]\. In the configuration object of the vLLM, we uset​e​m​p​e​r​a​t​u​r​e=0\.2temperature=0\.2, a balanced setting for creative yet coherent outputs\. To increase efficiency in the hybrid system, we use batch processing for the LLM agents\.

Because the output of the LLM has been used as feedback for ABM, it is essential to generate structured output from the LLM\[[14](https://arxiv.org/html/2607.06757#bib.bib30)\]\. We use the recently developed Outlines framework to obtain structured output in a yes/no format\[[35](https://arxiv.org/html/2607.06757#bib.bib31)\]\. If the LLM output of an LLM agentlil\_\{i\}isyesfor a weekwwfor a certain activity, then the corresponding ABM agents contacts in the networkG​\(t\)G\(t\)will remain unchanged\. If the LLM output of an LLM agentlil\_\{i\}isnofor a week,ww, for a certain activity, then the corresponding ABM agents contacts in the networkG​\(t\)G\(t\)will be deactivated\. Therefore, an infected agentnin\_\{i\}cannot infect its susceptible neighbors becausenin\_\{i\}will no longer have connections with them, as it will stay home\. ForΔ​tL​L​M\\Delta t\_\{LLM\}, we take week \(seven days\) as the observation data from CDC is weekly for COVID\-19\[[9](https://arxiv.org/html/2607.06757#bib.bib45)\]\. However,Δ​tL​L​M\\Delta t\_\{LLM\}is flexible and can take any value depending on the research question\.

##### LLM\-ABM server\-client integration on an HPC platform:

We integrate the LLM as a GPU\-backed inference server and ABMs as clients\. The LLM runs as a single\-process, single\-GPU inference server\. The communication between the LLM and ABM is performed by FastAPI, Uvicorn, and python library requests\. The framework runs separate Slurm jobs for the LLM and ABM\. At each LLM time stepΔ​tL​L​M\\Delta t\_\{LLM\}, ABM sends an update and query to the LLM server for the inference\. The ABM pauses the simulation for the LLM inference\. After receiving a response from the LLM server, it updates G\(t\) and resumes the ABM simulation\. Notably, one server can respond to multiple ABM clients at a time\.

## 3Results and Discussion

In this section, we describe the HALE framework applied to a real\-world epidemic scenario: COVID\-19 in Salt Lake County, UT\. To perform our experiments, we use one compute node from the Frontier HPE Cray EX supercomputer hosted at Oak Ridge National Laboratory\. This compute node features AMD MI250X GPUs\.

Our experiments simulate the COVID\-19 outbreak from September 2021 to February 2022\. We use an individual\-based SIR epidemic model\. We adopt the average basic reproduction number and recovery rate that correspond to the Omicron variant\. Previous studies report that the average basic reproduction number for the Omicron variant is 9\.5 \(inter quartile range: 7\.25, 11\.88\)\[[20](https://arxiv.org/html/2607.06757#bib.bib32)\], with a recovery rate of15​d​a​y​s−1\\frac\{1\}\{5\}days^\{\-1\}\[[32](https://arxiv.org/html/2607.06757#bib.bib33)\]\. Because we use an individual\-based network framework, we adjust the transmissibility for each edge according to the average degree\.

### 3\.1HALE framework simulation

As the HALE framework is stochastic by nature, we simulate it for 30 times and compare the aggregated results with real\-world observations\. Figure[3](https://arxiv.org/html/2607.06757#S3.F3)shows new weekly case data\. The red\-dotted lines represent individual simulations, and the red\-solid line represents the mean of all simulations\. The green line shows weekly COVID\-19 case data in Salt Lake County, UT\[[9](https://arxiv.org/html/2607.06757#bib.bib45)\]\. The peak time of the mean over simulation line follows the observation\. However, the area under this red curve or the total number of infections in the experiments is 2\.21 times higher than the observation \(green line\)\. Because the observation data mostly included the symptomatic cases and missed many asymptomatic cases, it is lower than the actual infected population\. Wang et al\. found that asymptomatic population is almost 44% of the total infection\[[34](https://arxiv.org/html/2607.06757#bib.bib34)\]\. Another study estimated that presymptomatic and asymptomatic cases together comprise at least 50% of the force of infection at the outbreak’s peak\[[29](https://arxiv.org/html/2607.06757#bib.bib35)\], which is consistent with our simulation results\.

![Refer to caption](https://arxiv.org/html/2607.06757v1/figure/epi_curve_abm_llm.png)\(a\)HALE
![Refer to caption](https://arxiv.org/html/2607.06757v1/figure/epi_curve_only_abm.png)\(b\)ABM\-only

Figure 3:Weekly new COVID\-19 cases from September 2021 to February 2022 in Salt Lake County, UT\. The green line shows weekly case data from the Centers for Disease Control and Prevention\. The red lines represent simulation outputs\. Subfigure \(a\) shows simulation results from the HALE framework, and subfigure \(b\) shows simulation results from the ABM\. The red\-dotted lines represent individual simulations, and the red\-solid line represents the mean of all simulations\.Figure[3](https://arxiv.org/html/2607.06757#S3.F3)b presents the simulation output from ABM\-only with the same epidemic model, parameters, and contact networks\. Within the HALE framework, the average number of links decreases as people choose to stay at home rather than go out in response to increasing disease levels\. Therefore, for a fair comparison, we randomly deactivate 35% of contacts in the ABM\-only simulation\. The ABM\-only simulation still fails to follow the actual observation in terms of peak timing and epidemic size\.

### 3\.2LLM inference

We use the offline Llama model for behavioral prediction of the LLM agents\. Our observations find that LLM reasoning is sensitive to prompts and the temperature hyperparameter\. We use a temperature setting for the Llama model that balances randomness and determinism to avoid overly random or overly deterministic outputs\. Very Low temperature settings act deterministically and introduce very little randomness when generating the next token\. Att​e​m​p​e​r​a​t​u​r​e=0temperature=0, the model behaves very conservatively, and the probability of declining outdoor activities for all age groups is almost always close to 1\. Higher temperatures introduce more randomness in the inference and reduce meaningful patterns in the generation\. For example, att​e​m​p​e​r​a​t​u​r​e=0\.7temperature=0\.7, the probability for all groups at any time is approximately 0\.5\. Therefore, our experiment uses a temperature of 0\.2\.

A sample prompt for an LLM agent,lil\_\{i\}, in the experiment is provided below\.

Prompt:Behavior prediction taskLocation:Salt Lake County, UTContext:0\.02% people infected by COVID\-19Person:Asian female, age 25–34, lives in Bluffdale, Salt Lake County, UTQuestion:Considering the infected percentage and a person’s behavior based on demographic factors, will this person go to a public place?

We calculate the probability of declining an outdoor activity by LLM agents from the HALE simulations\.

#### 3\.2\.1Analysis of LLM agents across different age groups over time

Figure[4](https://arxiv.org/html/2607.06757#S3.F4)shows the probability of declining an outdoor activity with time for different age groups\. The probability of declining outdoor activities increases for all groups over time as the number of infected people increases \(Figure[3](https://arxiv.org/html/2607.06757#S3.F3)\)\. The over 60 and 0–18 age groups consistently exhibit higher probabilities, whereas the 18–24 and 25–34 age groups show lower probabilities among all groups\. The LLM inference behaves as expected for age group analysis\.

![Refer to caption](https://arxiv.org/html/2607.06757v1/figure/llm_response_age.png)Figure 4:LLM inference for different age groups in the HALE framework\. Theyy\-axis shows the probability to decline participation in an outdoor activity\.Experiments show that the probability of declining an outdoor activity for all age groups increases slowly and does not exhibit any sharp peak, unlike the weekly infected cases shown in Figure[3](https://arxiv.org/html/2607.06757#S3.F3)\. This indicates that human behavior changes dynamically in response to surroundings, but the changes occur gradually\.

#### 3\.2\.2Analysis of LLM agents across different municipalities

In the analysis of LLM inference \(Figure[5](https://arxiv.org/html/2607.06757#S3.F5)\), residents of different municipalities or cities show different levels of response to an epidemic\. Figure[5](https://arxiv.org/html/2607.06757#S3.F5)also shows that information referring only to the county, without specifying a city, is associated with a higher probability of declining outdoor activity\. In the experiments, agents living in Salt Lake County, UT, without any city\-specific information have a probability of declining outdoor activity of around 0\.7 \(last bar in the Figure[5](https://arxiv.org/html/2607.06757#S3.F5)\)\.

![Refer to caption](https://arxiv.org/html/2607.06757v1/figure/llm_response_city_onebar.png)Figure 5:LLM inference in the HALE framework across different municipalities\. Theyy\-axis shows the probability to decline participation in an outdoor activity\.However, when the city of the LLM agent within Salt Lake County is specified, the probability is lower and takes different values\.

In the experiment, Bluffdale, Cottonwood Heights, Herriman town, and Taylorsville city have a higher probability of decline\. Bluffdale and Harriman are rapidly growing suburbs with many young families\. These cities also have a higher proportion of children \(under 18\) compared to the Utah state average\. On the other hand, Cottonwood Heights has a relatively older population on average\. Due to their population structure and specfic locations, these cities have a higher probability of declining an outdoor activity during an outbreak\.

Small unincorporated or census\-designated communities—Alta, Emigration Canyon, Brighton, Copperton, and White City—exhibit very small probabilities, mostly near zero\. Only the over 60 age group shows some activity \(Figure[6](https://arxiv.org/html/2607.06757#S3.F6)\)\. It indicates that the model doesn’t have enough information about these communities\.

![Refer to caption](https://arxiv.org/html/2607.06757v1/figure/llm_response_city_small.png)Figure 6:LLM inference within the HALE framework across small municipalities\. Theyy\-axis shows the probability to decline participation in an outdoor activity\.

## 4Conclusion

The individual\-based ABM network models are well known for their utility in scenario analysis and for informing policy decisions\. On the other hand, LLMs are trained on massive amounts of text data that can be used to infer human behavior\. Therefore, leveraging both methods is an effective approach for modeling, in near real–time, a real\-world problem\. To this end, we develop the HALE framework to strengthen individual\-based temporal network modeling by reducing a data gap through LLM reasoning\. Notably, HALE is scalable and can model millions of agents\.

LLMs contribute to modeling and simulation efforts across a range of domains and applications\. In this research, we use an LLM to predict human behavior\. Previous studies have assessed the performance of LLMs, particularly GPT\-4, in predicting human decisions\[[37](https://arxiv.org/html/2607.06757#bib.bib38),[26](https://arxiv.org/html/2607.06757#bib.bib39)\]\. These studies find that LLMs tend to reflect population\-level patterns rather than individual decision\-making, while also exhibiting biases that vary by context and group\. Moreover, LLMs have the potential to tailor their outputs so they can capture the nuances of specific countries or locations\[[27](https://arxiv.org/html/2607.06757#bib.bib37)\]\. In this work, we use an LLM to predict human behavior at the group level rather than at the individual level\. We also provide contextual information about location and demographic characteristics\. Here, we use a zero\-shot LLM inference to address previously unseen scenarios\. In future, we will explore RAG \(retrieval\-augmented generation\) to improve the inference\. Additionally, we plan to include weather data to improve predictions of outdoor activity\.

In the HALE framework, we address two data gaps in the individual\-based network model: partially available activity data with unknown specific activity locations, and static activity data\. To address the first gap, we propose a systematic method to develop a temporal network model from activity data to leverage the available activity information while making realistic, data\-dependent assumptions\. The temporal networks display degree patterns consistent with those observed in prior studies\[[10](https://arxiv.org/html/2607.06757#bib.bib18)\]\.

ABM uses activity or mobility data\. Unfortunately, collecting this activity data is a rigorous and time\-consuming process; therefore, ABMs often use outdated or static activity data that do not reflect dynamic situations\. To address this second data gap, we incorporate an LLM component in the HALE framework\. In this approach, the ABM simulates disease spreading, and the LLM predicts the group\-level human decision\-making and behavior\. Overall, our results indicate that the epidemic curves generated by the HALE framework align more closely with observed data for Salt Lake County, UT than the curves produced by the ABM\-only model\.

The HALE framework opens the door to utilizing the powerful tool of LLM reasoning in network simulation for policy making, without replacing well\-developed computational models, but instead assisting in addressing data gaps\.

## Ethics

This work did not require any ethical approval from a human subject or animal welfare committee\.

## Data access

The code and supplementary data are available in the public repository in Zenodo:\[[23](https://arxiv.org/html/2607.06757#bib.bib48)\]\.

## Competing

We declare we have no competing interest\.

## Funding

This material is based upon work supported by the US Department of Energy, Office of Science, Office of Advanced Scientific Computing under Award Number DE\-SC\-ERKJ422\.

## acknowledgment

This manuscript has been authored by UT\-Battelle LLC under contract DE\-AC05\-00OR22725 with the US Department of Energy\. The publisher acknowledges the US government license to provide public access under the DOE Public Access Plan \(http://energy\.gov/downloads/doe\-public\-access\-plan\)\.

## References

- \[1\]F\. H\. A\. \(FHWA\)\(2025\)National household travel survey\.Note:[https://nhts\.ornl\.gov/](https://nhts.ornl.gov/)Accessed: 2025\-08\-21Cited by:[§2\.1](https://arxiv.org/html/2607.06757#S2.SS1.p2.1)\.
- \[2\]S\. Afroj Moon, A\. Marathe, and A\. Vullikanti\(2023\)Are all underimmunized measles clusters equally critical?\.Royal Society Open Science10\(8\)\.Cited by:[§1](https://arxiv.org/html/2607.06757#S1.p1.1)\.
- \[3\]M\. Akinlotan, N\. Khodakarami, K\. Primm, J\. Bolin, and A\. O\. Ferdinand\(2023\)Travel for medical or dental care by race/ethnicity and rurality in the us: findings from the 2001, 2009 and 2017 national household travel surveys\.Preventive Medicine Reports35,pp\. 102297\.Cited by:[§2\.2\.2](https://arxiv.org/html/2607.06757#S2.SS2.SSS2.p2.10)\.
- \[4\]R\. Bachman\(2024\)How close do you need to be to your gym?\.Note:[https://www\.wsj\.com/articles/how\-close\-do\-you\-need\-to\-be\-to\-your\-gym\-1490111186](https://www.wsj.com/articles/how-close-do-you-need-to-be-to-your-gym-1490111186)Accessed: 2024\-01\-01Cited by:[§2\.2\.2](https://arxiv.org/html/2607.06757#S2.SS2.SSS2.p2.10)\.
- \[5\]E\. Blumenberg, M\. Wander, and Z\. Yao\(2024\)Decisions & distance: the relationship between child care access and child care travel\.Journal of transport geography114,pp\. 103756\.Cited by:[§2\.2\.2](https://arxiv.org/html/2607.06757#S2.SS2.SSS2.p2.10)\.
- \[6\]I\. Brodsky\(2025\)H3: uber’s hexagonal hierarchical spatial index\.Note:[https://www\.uber\.com/blog/h3/](https://www.uber.com/blog/h3/)Accessed: 2025\-08\-21Cited by:[§2\.1](https://arxiv.org/html/2607.06757#S2.SS1.p2.1)\.
- \[7\]A\. Chopra, J\. Subramanian, B\. Krishnamurthy, and R\. RaskarAgentTorch: agent\-based modeling with automatic differentiation\.InSecond Agent Learning in Open\-Endedness Workshop,Cited by:[§1](https://arxiv.org/html/2607.06757#S1.p1.1),[§1](https://arxiv.org/html/2607.06757#S1.p3.1)\.
- \[8\]G\. Civitarese, M\. Fiori, P\. Choudhary, and C\. Bettini\(2025\)Large language models are zero\-shot recognizers for activities of daily living\.ACM Transactions on Intelligent Systems and Technology16\(4\),pp\. 1–32\.Cited by:[§2\.3\.3](https://arxiv.org/html/2607.06757#S2.SS3.SSS3.Px2.p3.1)\.
- \[9\]E\. Y\. Cramer, Y\. Huang, Y\. Wang, E\. L\. Ray, M\. Cornell, J\. Bracher, A\. Brennen, A\. J\. Castro Rivadeneira, A\. Gerding, K\. House, D\. Jayawardena, A\. H\. Kanji, A\. Khandelwal, K\. Le, J\. Niemi, A\. Stark, A\. Shah, N\. Wattanachit, M\. W\. Zorn, N\. G\. Reich, and U\. C\. F\. H\. Consortium\(2021\)The united states covid\-19 forecast hub dataset\.medRxiv\.External Links:[Document](https://dx.doi.org/10.1101/2021.11.04.21265886),[Link](https://www.medrxiv.org/content/10.1101/2021.11.04.21265886v1)Cited by:[§2\.3\.3](https://arxiv.org/html/2607.06757#S2.SS3.SSS3.Px2.p4.10),[§3\.1](https://arxiv.org/html/2607.06757#S3.SS1.p1.1)\.
- \[10\]S\. Y\. Del Valle, J\. M\. Hyman, H\. W\. Hethcote, and S\. G\. Eubank\(2007\)Mixing patterns between age groups in social networks\.Social Networks29\(4\),pp\. 539–554\.Cited by:[§2\.2\.1](https://arxiv.org/html/2607.06757#S2.SS2.SSS1.Px2.p3.1),[§2\.2\.1](https://arxiv.org/html/2607.06757#S2.SS2.SSS1.Px3.p1.5),[§2\.2\.4](https://arxiv.org/html/2607.06757#S2.SS2.SSS4.p2.1),[§4](https://arxiv.org/html/2607.06757#S4.p3.1)\.
- \[11\]A\. Dubey, A\. Jauhri, A\. Pandey, A\. Kadian, A\. Al\-Dahle, A\. Letman, A\. Mathur, A\. Schelten, A\. Yang, A\. Fan,et al\.\(2024\)The llama 3 herd of models\.arXiv e\-prints,pp\. arXiv–2407\.Cited by:[§2\.3\.3](https://arxiv.org/html/2607.06757#S2.SS3.SSS3.Px2.p3.1)\.
- \[12\]A\. Ferraro, A\. Galli, V\. La Gatta, M\. Postiglione, G\. M\. Orlando, D\. Russo, G\. Riccio, A\. Romano, and V\. Moscato\(2024\)Agent\-based modelling meets generative ai in social network simulations\.InInternational Conference on Advances in Social Networks Analysis and Mining,pp\. 155–170\.Cited by:[§1](https://arxiv.org/html/2607.06757#S1.p3.1)\.
- \[13\]N\. C\. for Education Statistics\(2025\)Average public school class size: average class size in public k–12 schools, by school level, class type, and state: 2020–21\.Note:[https://nces\.ed\.gov/surveys/ntps/estable/table/ntps/ntps2021\_sflt07\_t1s](https://nces.ed.gov/surveys/ntps/estable/table/ntps/ntps2021_sflt07_t1s)Accessed: 2025\-06\-01Cited by:[§2\.2\.1](https://arxiv.org/html/2607.06757#S2.SS2.SSS1.Px2.p2.1)\.
- \[14\]S\. Geng, H\. Cooper, M\. Moskal, S\. Jenkins, J\. Berman, N\. Ranchin, R\. West, E\. Horvitz, and H\. Nori\(2025\)Jsonschemabench: a rigorous benchmark of structured outputs for language models\.arXiv preprint arXiv:2501\.10868\.Cited by:[§2\.3\.1](https://arxiv.org/html/2607.06757#S2.SS3.SSS1.p3.7),[§2\.3\.3](https://arxiv.org/html/2607.06757#S2.SS3.SSS3.Px2.p4.10)\.
- \[15\]N\. Ghaffarzadegan, A\. Majumdar, R\. Williams, and N\. Hosseinichimeh\(2023\)Generative agent\-based modeling: unveiling social system dynamics through coupling mechanistic models with generative artificial intelligence\.arXiv preprint arXiv:2309\.11456\.Cited by:[§1](https://arxiv.org/html/2607.06757#S1.p1.1)\.
- \[16\]K\. O\. Kwok, T\. Huynh, W\. I\. Wei, S\. Y\. Wong, S\. Riley, and A\. Tang\(2024\)Utilizing large language models in infectious disease transmission modelling for public health preparedness\.Computational and Structural Biotechnology Journal23,pp\. 3254–3257\.Cited by:[§1](https://arxiv.org/html/2607.06757#S1.p3.1)\.
- \[17\]W\. Kwon, Z\. Li, S\. Zhuang, Y\. Sheng, L\. Zheng, C\. H\. Yu, J\. Gonzalez, H\. Zhang, and I\. Stoica\(2023\)Efficient memory management for large language model serving with pagedattention\.InProceedings of the 29th symposium on operating systems principles,pp\. 611–626\.Cited by:[§2\.3\.3](https://arxiv.org/html/2607.06757#S2.SS3.SSS3.Px2.p3.1)\.
- \[18\]J\. Liu, C\. Chu, Y\. Zhao, G\. Aoki, and Z\. Xiao\(2025\)Agentic ai for sustainable development: leveraging large language model\-enhanced agent\-based modeling for complex policy strategies\.Emerging Media3\(3\),pp\. 401–413\.Cited by:[§1](https://arxiv.org/html/2607.06757#S1.p5.1)\.
- \[19\]J\. L\. Liu, B\. Han, and D\. A\. Cohen\(2015\)Beyond neighborhood food environments: distance traveled to food establishments in 5 us cities, 2009–2011\.Preventing chronic disease12,pp\. E126\.Cited by:[§2\.2\.2](https://arxiv.org/html/2607.06757#S2.SS2.SSS2.p2.10)\.
- \[20\]Y\. Liu and J\. Rocklöv\(2022\)The effective reproductive number of the omicron variant of sars\-cov\-2 is several times relative to delta\.Journal of travel medicine29\(3\),pp\. taac037\.Cited by:[§3](https://arxiv.org/html/2607.06757#S3.p2.1)\.
- \[21\]S\. A\. Moon, L\. W\. Cohnstaedt, D\. S\. McVey, and C\. M\. Scoglio\(2019\)A spatio\-temporal individual\-based network framework for west nile virus in the usa: spreading pattern of west nile virus\.PLoS computational biology15\(3\),pp\. e1006875\.Cited by:[§2\.2\.2](https://arxiv.org/html/2607.06757#S2.SS2.SSS2.p1.1)\.
- \[22\]S\. A\. Moon and C\. M\. Scoglio\(2021\)Contact tracing evaluation for covid\-19 transmission in the different movement levels of a rural college town in the usa\.Scientific reports11\(1\),pp\. 4891\.Cited by:[§1](https://arxiv.org/html/2607.06757#S1.p4.1),[§2\.2\.1](https://arxiv.org/html/2607.06757#S2.SS2.SSS1.Px3.p1.5)\.
- \[23\]S\. A\. Moon\(2026\-04\)HALE data and simulation\.Zenodo\.External Links:[Document](https://dx.doi.org/10.5281/zenodo.19829913),[Link](https://doi.org/10.5281/zenodo.19829913)Cited by:[Data access](https://arxiv.org/html/2607.06757#Sx2.p1.1)\.
- \[24\]C\. C\. Mundt, K\. E\. Sackett, L\. D\. Wallace, C\. Cowger, and J\. P\. Dudley\(2009\)Long\-distance dispersal and accelerating waves of disease: empirical relationships\.The American Naturalist173\(4\),pp\. 456–466\.Cited by:[§2\.2\.2](https://arxiv.org/html/2607.06757#S2.SS2.SSS2.p1.1)\.
- \[25\]J\. Ozik, J\. M\. Wozniak, N\. Collier, C\. M\. Macal, and M\. Binois\(2021\)A population data\-driven workflow for covid\-19 modeling and learning\.The International Journal of High Performance Computing Applications35\(5\),pp\. 483–499\.Cited by:[§2\.3\.3](https://arxiv.org/html/2607.06757#S2.SS3.SSS3.Px1.p2.1)\.
- \[26\]H\. Peters and S\. C\. Matz\(2024\)Large language models can infer psychological dispositions of social media users\.PNAS nexus3\(6\),pp\. pgae231\.Cited by:[§4](https://arxiv.org/html/2607.06757#S4.p2.1)\.
- \[27\]Y\. Qu and J\. Wang\(2024\)Performance and biases of large language models in public opinion simulation\.Humanities and Social Sciences Communications11\(1\),pp\. 1–13\.Cited by:[§4](https://arxiv.org/html/2607.06757#S4.p2.1)\.
- \[28\]A\. Spannaus, S\. A\. Moon, J\. Gounley, and H\. A\. Hanson\(2025\)Data assimilation for robust uq within agent\-based simulation on hpc systems\.InProceedings of the Platform for Advanced Scientific Computing Conference,pp\. 1–11\.Cited by:[§2\.3\.3](https://arxiv.org/html/2607.06757#S2.SS3.SSS3.Px1.p1.2),[§2\.3\.3](https://arxiv.org/html/2607.06757#S2.SS3.SSS3.Px1.p2.1)\.
- \[29\]R\. Subramanian, Q\. He, and M\. Pascual\(2021\)Quantifying asymptomatic infection and transmission of covid\-19 in new york city using observed cases, serology, and testing capacity\.Proceedings of the National Academy of Sciences118\(9\),pp\. e2019716118\.Cited by:[§3\.1](https://arxiv.org/html/2607.06757#S3.SS1.p1.1)\.
- \[30\]J\. Tuccillo\(2024\)Downscaling synthetic populations to realistic residential locations\.Technical reportOak Ridge National Laboratory \(ORNL\), Oak Ridge, TN \(United States\)\.Cited by:[§2\.1](https://arxiv.org/html/2607.06757#S2.SS1.p3.1)\.
- \[31\]J\. Tuccillo, R\. Stewart, A\. Rose, N\. Trombley, J\. Moehl, N\. Nagle, and B\. Bhaduri\(2023\)UrbanPop: a spatial microsimulation framework for exploring demographic influences on human dynamics\.Applied Geography151,pp\. 102844\.Cited by:[§2\.1](https://arxiv.org/html/2607.06757#S2.SS1.p1.1)\.
- \[32\]UKHSA\(2023\)COVID\-19 omicron variant infectious period and transmission from people with asymptomatic compared with symptomatic infection: a rapid review\.GOV\-14430\.\{\\\{UK Health Security Agency\}\\\}\.Cited by:[§3](https://arxiv.org/html/2607.06757#S3.p2.1)\.
- \[33\]A\. S\. Vezhnevets, J\. P\. Agapiou, A\. Aharon, R\. Ziv, J\. Matyas, E\. A\. Duéñez\-Guzmán, W\. A\. Cunningham, S\. Osindero, D\. Karmon, and J\. Z\. Leibo\(2023\)Generative agent\-based modeling with actions grounded in physical, social, or digital space using concordia\.arXiv preprint arXiv:2312\.03664\.Cited by:[§1](https://arxiv.org/html/2607.06757#S1.p1.1)\.
- \[34\]B\. Wang, P\. Andraweera, S\. Elliott, H\. Mohammed, Z\. Lassi, A\. Twigger, C\. Borgas, S\. Gunasekera, S\. Ladhani, and H\. S\. Marshall\(2023\)Asymptomatic sars\-cov\-2 infection by age: a global systematic review and meta\-analysis\.The Pediatric Infectious Disease Journal42\(3\),pp\. 232–239\.Cited by:[§3\.1](https://arxiv.org/html/2607.06757#S3.SS1.p1.1)\.
- \[35\]B\. T\. Willard and R\. Louf\(2023\)Efficient guided generation for large language models\.arXiv preprint arXiv:2307\.09702\.Cited by:[§2\.3\.3](https://arxiv.org/html/2607.06757#S2.SS3.SSS3.Px2.p4.10)\.
- \[36\]R\. Williams, N\. Hosseinichimeh, A\. Majumdar, and N\. Ghaffarzadegan\(2023\)Epidemic modeling with generative agents\.arXiv preprint arXiv:2307\.04986\.Cited by:[§1](https://arxiv.org/html/2607.06757#S1.p3.1),[§1](https://arxiv.org/html/2607.06757#S1.p5.1)\.
- \[37\]F\. Xiao and X\. X\. Wang\(2025\)Evaluating the ability of large language models to predict human social decisions\.Scientific Reports15\(1\),pp\. 32290\.Cited by:[§4](https://arxiv.org/html/2607.06757#S4.p2.1)\.
- \[38\]H\. L\. Yang, M\. Laverdiere, T\. Hauser, B\. Swan, E\. Schmidt, J\. Moehl, A\. Reith, D\. Adams, B\. Morris, J\. McKee,et al\.\(2024\)A baseline structure inventory with critical attribution for the us and its territories\.Scientific Data11\(1\),pp\. 502\.Cited by:[§2\.1](https://arxiv.org/html/2607.06757#S2.SS1.p3.1)\.
- \[39\]L\. Zhang, Y\. Hu, W\. Li, Q\. Bai, and P\. Nand\(2025\)LLM\-aidsim: llm\-enhanced agent\-based influence diffusion simulation in social networks\.Systems13\(1\),pp\. 29\.Cited by:[§1](https://arxiv.org/html/2607.06757#S1.p1.1)\.
- \[40\]Y\. Zhang, S\. A\. Khan, A\. Mahmud, H\. Yang, A\. Lavin, M\. Levin, J\. Frey, J\. Dunnmon, J\. Evans, A\. Bundy,et al\.\(2025\)Exploring the role of large language models in the scientific method: from hypothesis to discovery\.npj Artificial Intelligence1\(1\),pp\. 14\.Cited by:[§1](https://arxiv.org/html/2607.06757#S1.p2.1)\.

Similar Articles

Simulate, Reason, Decide: Scientific Reasoning with LLMs for Simulation-Driven Decision Making

arXiv cs.AI

Researchers from the University of Michigan introduce MechSim, a mechanism-grounded neuro-symbolic reasoning framework that enables LLM agents to reason about the internal assumptions, dependencies, and execution behavior of scientific simulators rather than treating them as black boxes. The framework improves explanation quality and decision-making reliability across high-stakes domains like healthcare, finance, and public policy.

HawkesLLM: Semantic Uncertainty Propagation in Agentic Text Simulation

arXiv cs.CL

This paper introduces HawkesLLM, a framework that models semantic uncertainty propagation in multi-step agentic text simulations by combining a multivariate Hawkes process for temporal influence and memory selection with a language model for text generation. Evaluation on a GDELT news-cascade case study shows improved late-stage semantic alignment under compact prompt-memory constraints.