Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Population Dynamics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Intra- and inter-subtype HIV diversity between 1994 and 2018 in southern Uganda: a longitudinal population-based study

There is limited data on human immunodeficiency virus (HIV) evolutionary trends in African populations. We evaluated changes in HIV viral diversity and genetic divergence in southern Uganda over a 24-year period spanning the introduction and scale-up of HIV prevention and treatment programs using HIV sequence and survey data from the Rakai Community Cohort Study, an open longitudinal population-based HIV surveillance cohort. Gag (p24) and env (gp41) HIV data were generated from people living with HIV (PLHIV) in 31 inland semi-urban trading and agrarian communities (1994–2018) and four hyperendemic Lake Victoria fishing communities (2011–2018) under continuous surveillance. HIV subtype was assigned using the Recombination Identification Program with phylogenetic confirmation. Inter-subtype diversity was evaluated using the Shannon diversity index, and intra-subtype diversity with the nucleotide diversity and pairwise TN93 genetic distance. Genetic divergence was measured using root-to-tip distance and pairwise TN93 genetic distance analyses. Demographic history of HIV was inferred using a coalescent-based Bayesian Skygrid model. Evolutionary dynamics were assessed among demographic and behavioral population subgroups, including by migration status. 9931 HIV sequences were available from 4999 PLHIV, including 3060 and 1939 persons residing in inland and fishing communities, respectively. In inland communities, subtype A1 viruses proportionately increased from 14.3% in 1995 to 25.9% in 2017 (P < .001), while those of subtype D declined from 73.2% in 1995 to 28.2% in 2017 (P < .001). The proportion of viruses classified as recombinants significantly increased by nearly four-fold from 12.2% in 1995 to 44.8% in 2017. Inter-subtype HIV diversity has generally increased. While intra-subtype p24 genetic diversity and divergence leveled off after 2014, intra-subtype gp41 diversity, effective population size, and divergence increased through 2017. Intra- and inter-subtype viral diversity increased across all demographic and behavioral population subgroups, including among individuals with no recent migration history or extra-community sexual partners. This study provides insights into population-level HIV evolutionary dynamics following the scale-up of HIV prevention and treatment programs. Continued molecular surveillance may provide a better understanding of the dynamics driving population HIV evolution and yield important insights for epidemic control and vaccine development.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Quantitative multiplexed analysis of gene and protein expression patterns in Yarrowia lipolytica

In this report, we present coordinated observations of protein and mRNA transcript counts at the single-cell level in the oleaginous yeast model Yarrowia lipolytica. The transcription factor Xbp1p regulates entry into a quiescent state, representing a shift of resources to sequestration of nutrients rather than cell division. We observed the responses of wild-type and Δxbp1 cells to protein (by fluorescence) and transcript quantification and localization at both single-cell and population-averaged levels. Data were collected via single-molecule fluorescence in situ hybridization (smFISH) and qPCR under nitrogen depletion, a condition that drives lipid accumulation. These techniques reveal a complex and heterogeneous population of Xbp1p dynamics and downstream regulation. Our findings highlight the need for single-cell resolution analyses to describe cellular dynamics and regulatory processes.

Yarrowia lipolytica↗
DOE OSTI2025

Morphotype-resolved characterization of microalgal communities in a nutrient recovery process with ARTiMiS flow imaging microscopy

Microalgae-driven nutrient recovery represents a promising technology for phosphorus removal from wastewater while simultaneously generating biomass that can be valorized to offset treatment costs. As full-scale processes come online, system parameters including biomass composition must be carefully monitored to optimize performance and prevent culture crashes. In this study, flow imaging microscopy (FIM) was leveraged to characterize microalgal community composition in near real-time at a full-scale municipal wastewater treatment plant (WWTP) in Wisconsin, USA, and population and morphotype dynamics were examined to identify relationships between water chemistry, biomass composition, and system performance. Two FIM technologies, FlowCam and ARTiMiS, were evaluated as monitoring tools. ARTiMiS provided a more accurate estimate of total system biomass, and estimates derived from particle area as a proxy for biovolume yielded better approximations than particle counts. Deep learning classification models trained on annotated image libraries demonstrated equivalent performance between FlowCam and ARTiMiS, and convolutional neural network (CNN) classifiers proved significantly more accurate when compared to feature table-based dense neural network (DNN) models. Across a two-year study period, Scenedesmus spp. appeared most important for phosphorus removal, and were negatively impacted by elevated temperatures and increase in nitrite/nitrate concentrations. Chlorella and Monoraphidium also played an important role in phosphorus removal. For both Scenedesmus and Chlorella, smaller morphological types were more often associated with better system performance, whereas larger morphotypes likely associated with stress response(s) correlated with poor phosphorus recovery rates. Furthermore, these results demonstrate the potential of FIM as a critical technology for high-resolution characterization of industrial microalgal processes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Evidence for Photoinduced Polaron Generation in a High Persistence Length Low Bandgap Conjugated Polymer in Solution

Using ultrafast time-resolved infrared (TRIR) spectroscopy, we studied the solution-phase excited-state structural evolution of an indacenodithiophene-co-benzothiadiazole polymer (C8-IDTBT). Following band gap excitation, the TRIR spectra reveal vibrational features that develop within 10 ps and decay over 4 ns. Using pulse radiolysis measurements, charge-modulation spectroscopy, and quantum-chemical calculations, the IR features are assigned to polaron pairs. Interestingly, these features appear on an evolving broad mid-IR electronic absorption background, with kinetics correlating with the formation and decay of the cationradical vibrational bands. A three-state kinetic model successfully reproduces the spectral evolution, revealing that the polaron and exciton populations exist in dynamic equilibrium on picosecond time scales, with time constants for exciton dissociation in the range of 3− 5 ps and polaron-to-exciton reformation between 20 and 100 ps, while both species decay to the ground state on much slower nanosecond time scales (∼1 ns), yielding a remarkably high polaron-generation efficiency, higher than 50%. These findings provide fundamental insights into intramolecular charge photogeneration mechanisms in conjugated polymers, demonstrating efficient bound-polaron formation in isolated polymer chains.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Surface chlorination of IrO2(110) by HCl

The ability to controllably chlorinate metal-oxide surfaces can provide opportunities for designing selective oxidation catalysts. In the present study, we investigated the surface chlorination of IrO2(110) by HCl using temperature programmed reaction spectroscopy (TPRS), x-ray photoelectron spectroscopy (XPS), and density functional theory (DFT) calculations. We find that exposing IrO2(110) to HCl, followed by heating to 650 K in ultrahigh vacuum, produces nearly equal quantities of on-top and bridging Cl atoms on the surface, Clt and Clbr, where the Clbr atoms replace O-atoms that are removed from the surface by H2O formation. After HCl adsorption at 85 K, only H2O desorbs at low Cl coverages during TPRS, but HCl begins to desorb in increasing yields as the Cl coverage is increased above about 0.5 monolayer (ML). The desorption of Cl2 was not observed under any conditions, in good agreement with the high barrier for this reaction predicted by DFT. A maximum Cl coverage of 1 ML, with nearly equal coverages of Clt and Clbr atoms, could be generated by reacting HCl with IrO2(110) in UHV. Our results suggest that a kinetic competition between recombinative HCl and H2O desorption under the conditions studied limits the saturation Cl coverage to a value less than the 2 ML maximum predicted by thermodynamics. XPS further shows that the partitioning of Cl between the Clt and Clbr states can be altered by subjecting partially chlorinated IrO2(110) to reductive or oxidative treatments, demonstrating that the Cl site population can change dynamically in response to the gas environment. Our results provide insights for understanding the chlorination of IrO2(110) by HCl and can enable future experimental studies to determine how Cl-modification alters the surface chemical reactivity of IrO2(110) and potentially enhances selectivity toward partial oxidation chemistry.

Chemistry↗
DOE OSTI2026

Semi-Analytical Hierarchical Bayesian Inference of Nonlinear Model Structure in Stochastic Dynamics: Applied to Compartmental Models of Infectious Diseases

A Bayesian computational framework for parsimonious inference in stochastic nonlinear dynamical systems is presented. This framework enables the concurrent estimation of system states, time-varying parameters, time-invariant parameters, and the optimal sparsity structure of the model parameters. Because differential equation-based models are often simplified mechanistic or phenomenological representations, robust inference from noisy measurement data requires explicit treatment of model error and uncertainty. Model error and time-varying parameters can be represented as random processes, enabling inference while making minimal assumptions about the underlying sources of discrepancy and variability. Adopting stochastic differential equation representations affords the model significant flexibility, but can also render it susceptible to overfitting during statistical inversion, where the inferred model may track noise rather than the underlying signal. To alleviate the effects of overfitting and to enable the discovery of the optimal sparse representation of the time-invariant parameters, a Bayesian sparse learning algorithm is embedded within the framework. This sparse learning framework adopts an approximate hierarchical Bayesian setting defined by a series of semi-analytical expressions. The model structure inference framework is validated using a stochastic compartmental model for tracking and forecasting active cases of an infectious disease. Compartmental models describe population-level infectious disease dynamics through interactions among population fractions grouped by disease state. Mathematically, such models consist of a system of coupled ordinary differential equations. This example adopts an expressive compartmental model that includes multiple possible interactions between disease states, motivated by early uncertainty surrounding COVID-19 reinfection dynamics and their implications for long-term epidemic forecasting. The sparse learning exercise permits the inference of a priori unknown epidemiological dynamics from simulated public health data, discovering the nested compartmental model that optimizes the trade-off between average data-fit and model complexity. It is shown that inducing sparsity among the model parameters eliminates redundant interactions between compartments, equivalently revealing the optimal coupling structure between differential equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Diabatization with Electrostatic Embedding for Studying Photophysics in Organic Molecular Crystals

Highly emissive organic molecular crystals find applications in several areas, such as organic electronics, solar cells, and sensors. Understanding the excited-state mechanisms underlying these applications is essential for optimizing and controlling them effectively. Exciton models coupled with nonadiabatic dynamics, particularly quantum dynamics, provide crucial insights into photochemical and photophysical processes in molecular crystals. Nevertheless, there remains a lack of general tools and automated workflows to facilitate such simulations. In this paper, we present a computational strategy to investigate the photoactivated dynamics of organic molecular crystals, bridging methodologies traditionally used for molecular systems and materials science, with a particular focus on the interplay between local excitations and charge transfer (CT) processes. We have implemented an interface between the fromage and Overdia programs, enabling the construction of vibronic Hamiltonians for molecular crystals within an excited-state ONIOM(QM:QM′) framework, incorporating long-range electrostatics through a RESP-based Ewald summation. Fragment-based diabatization provides a route to quantum dynamics simulations in weak-to-intermediate coupling regimes. The method was applied to the photophysics of dibenzo[g,p]chrysene (DBC) crystals using time-dependent DFT. The fromage/ Overdia interface was employed to compute the couplings of local excitations and CT states for 18 unique DBC dimers in the crystal and to quantify the influence of electrostatic embedding, which was found to be modest (10−20%). Simulations on π-stacked dimers reproduced the small red shift observed experimentally from solution to crystal, attributed to electronic interactions among fixed monomers rather than crystal electrostatics. Quantum dynamics simulations revealed ultrafast population transfer from bright local excitations to CT states. This approach establishes a robust framework linking molecular and solid-state excited-state dynamics, with potential applications for studying excitations, defects, and impurities in molecular crystals.

Crystals↗
DOE OSTI2025

Mapping structures and dynamics with frequency-correlated diffusion exchange

Understanding molecular motion in diffusion-driven complex environments is critical for designing sustainable materials and improving chemical processes. Here, we introduce a multidimensional nuclear magnetic resonance (NMR) method that captures how molecular populations exchange across different dynamic regimes. By extending the modulated gradient spin-echo technique to include frequency-frequency correlations, our approach reveals diffusion pathways that are otherwise obscured in heterogeneous systems. Implemented on a unilateral NMR magnet, the method eliminates gradient pulsing constraints and accesses dynamics in the kilohertz regime. We apply this technique to swelling and acid-catalyzed deconstruction of cross-linked and linear polymers to observe how structural heterogeneity evolves over time. By linking molecular motion to topology and chemical state, we extract physical metrics such as fractal surface dimensionality and reaction wavefront velocity, properties inaccessible with standard diffusion measurements. This work expands the capabilities of NMR for probing soft matter, with implications for polymer recycling and materials design.

Fricke, Sophia N. [University of California, Berke↗
DOE OSTI2025

Using a recirculating anaerobic dynamic membrane bioreactor to treat hydrothermal liquefaction aqueous by-product

Hydrothermal liquefaction (HTL) has the potential to improve resource recovery at water resource recovery facilities (WRRF), but the production of a high-strength aqueous by-product (HTL-aq) is hampering HTL implementation. The formation of biofilms in anaerobic digestion have been shown to be useful when degrading recalcitrant compounds present in HTL-aq due to the promotion of direct interspecies electron transfer (DIET) and increase in the microbial activity of syntrophic and methanogenic populations. The Recirculating Anaerobic Dynamic Membrane Bioreactor (RAnDMBr) was able to degrade 65% of the chemical oxygen demand (COD) at 1.5 ± 0.2 g COD LR −1 day −1 and 5.6 ± 2.3 days producing 0.19 ± 0.02 LCH4 gCOD fed −1 . However, adding a solution rich in nutrients on a daily basis was necessary. The system presented microbial populations able to degrade aromatic compounds (i.e., Anaerolinaceae) to perform DIET and syntrophy (i.e., Syntrophus) and methanogens (i.e., Methanobacterium and Methanosarcina) with the biofilm having a higher relative abundance of methanogens than the suspended biomass. Increasing the organic loading rate to 2 g COD LR −1 day −1 caused inhibition in the system by accumulation of volatile fatty acids, probably due to an increase in phenol, N-heterocyclic and aromatic compounds. Overall, this research shows that the RAnDMBr can be used to treat HTL-aq in WRRF without inhibition at OLRs of 1.5 ± 0.2 g COD LR −1 day −1 or lower, making HTL-aq treatment more feasible. Future research should focus co-digestion of HTL-aq with a co-substrate rich in nutrients and on fouling mitigation strategies that will allow to increase the recirculation ratio to promote advective substrate transport through the biofilm.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

The Dark Matter Content of Milky Way Dwarf Spheroidal Galaxies: Draco, Sextans, and Ursa Minor

The Milky Way Survey of the Dark Energy Spectroscopic Instrument (DESI) has so far observed three classical dwarf spheroidal galaxies (dSphs): Draco, Sextans, and Ursa Minor. Based on the observed line-of-sight velocities and metallicities of their member stars, we apply the axisymmetric Jeans Anisotropic Multi-Gaussian Expansion modeling (JAM) approach to recover their inner dark matter distributions. In particular, both the traditional single-population Jeans model and the multiple population chemodynamical model are adopted. With the chemodynamical model, we divide member stars of each dSph into metal-rich and metal-poor populations. The metal-rich populations are more centrally concentrated and dynamically colder, featuring lower velocity dispersion profiles than the metal-poor populations. We find a diversity of the inner density slopes γ of dark matter halos, with the best constraints by the single-population or chemodynamical models consistent with each other. The inner density slopes are $0.7{1}_{-0.35}^{+0.34}$, $0.2{6}_{-0.12}^{+0.22}$, and $0.3{3}_{-0.16}^{+0.20}$ for Draco, Sextans, and Ursa Minor, respectively. We also present the measured astrophysical J and D factors of the three dSphs. Our results indicate that the study of the dark matter content of dSphs through stellar kinematics is still subject to uncertainties behind both the methodology and the observed data, through comparisons with previous measurements and datasets.

Dark matter↗
DOE OSTI2026

Noncanonical folding of peptoid oligomers: Formation of a closed conformation in nonpolar solvent

Peptoids provide a versatile platform for foldamer design, yet their conformational behavior in low-dielectric media remains poorly understood. The structural characteristics of N-1-phenylethylglycine (Nspe) homo-oligomers were inves-tigated in chloroform, a solvent that mimics the interior of lipid bilayers, using nuclear magnetic resonance (NMR) spec-troscopy and molecular dynamics (MD) simulations. Nspe7 populated two distinct compact conformations, while Nspe10 adopted a single homogeneous conformation related to the previously reported Nspe9 threaded-loop structure. Integrated experimental and computational analysis reveals that these structures are stabilized by cooperative end-to-end intramo-lecular hydrogen bonding, cis-trans backbone isomerism, and hydrophobic side-chain shielding. The resulting structures minimize exposed polar surface area, demonstrating a closed conformation in the low-dielectric environment. These find-ings establish specific chain-length requirements for achieving well-defined closed conformations. This work provides insights into peptoid folding in nonpolar media, enabling rational design strategies for solvent-directed conformational switching systems.

Oh, Jinyoung↗
DOE OSTI2025

Nonlinearities in Magnetic Confinement, Ionospheric Physics, and Population Explosion Leading to Profile Resilience Нелінійності в магнетному утриманні, фізиці іоносфери та процесі демографічного вибуху, які приводять до стійкості профілю

Nonlinearities play an important role in many fields. In the field of thermonuclear fusion, they are involved in questions such as profile resilience and fluid closure. A nonlinear phenomenon common to both fusion and astrophysical planets is the generation of zonal flows. These flows play a significant role in determining the level of turbulence and fluid closure in fusion. The effects of resonance broadening and nonlinearities are investigated, specifically focusing on the case of nonlinear instability that has appeared in drift waves. Similarities and differences between our systems are discussed, with population explosion and the dynamics of nonlinear systems for drift waves by different states in profile resilience described with great precision. The aim of our study is to put our fluid model for drift waves in tokamaks within the wider framework of statistical physics principles. This reinforces our belief in the broad application of our drift wave model, which encompasses current tokamaks, ITER, and the fusion pilot plant.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Design of facilitated dissociation enables timing of cytokine signalling

Protein design has focused on the design of ground states, ensuring that they are sufficiently low energy to be highly populated. Designing the kinetics and dynamics of a system requires, in addition, the design of excited states that are traversed in transitions from one low-lying state to another. This is a challenging task because such states must be sufficiently strained to be poorly populated, but not so strained that they are not populated at all, and because protein design methods have focused on generating near-ideal structures. Here we describe a general approach for designing systems that use an induced-fit power stroke to generate a structurally frustrated and strained excited state, allosterically driving protein complex dissociation. X-ray crystallography, double electron–electron resonance spectroscopy and kinetic binding measurements show that incorporating excited states enables the design of effector-induced increases in dissociation rates as high as 5,700-fold. We highlight the power of this approach by designing rapid biosensors, kinetically controlled circuits and cytokine mimics that can be dissociated from their receptors within seconds, enabling dissection of the temporal dynamics of interleukin-2 signalling.

deformation dynamics↗
DOE OSTI2025

Shake-Down Spectroscopy as State- and Site-Specific Probe of Ultrafast Chemical Dynamics

Tracking the multifarious ultrafast electronic and structural changes occurring in a molecule during a photochemical transformation is a challenging endeavor that benefits from recent experimental and computational progress in time-resolved techniques. Measurements of valence electronic states, which provide a global picture of the bonding structure of the molecule, and core electronic states, which provide insight into the local environment, traditionally require different approaches and are often studied separately. Here, we demonstrate that X-ray pulses from a seeded free-electron laser (FEL) enable the measurement of high-resolution, time-resolved X-ray photoelectron spectra (XPS) that capture weak satellite states resulting from shake-down processes in a valence-excited molecule. This approach effectively combines the advantages of both valence- and core-state investigations. We applied this method to investigate photoexcited CS2 molecules, where the role of internal conversion (IC) and intersystem crossing (ISC) in determining the predissociation dynamics is controversial. We present XPS spectra from photoexcited CS2, obtained at the FERMI FEL. High-resolution measurements, compared to the corresponding spectra obtained from accurate multireference quantum chemical calculations, reveal that shake-down satellite channels are highly sensitive to both valence electronic and geometric changes. Previous studies of the predissociation dynamics have led to uncertain assignments of the branching between singlet and triplet excited states. We derive a propensity rule that demonstrates the spin-selectivity of the shake-downs. This selectivity allows us to unequivocally assign contributions from the bright and dark singlet excited states, with populations tracked along the predissociation dynamic pathway.

Thompson, Henry J. [Univ. of Southampton (United K↗
DOE OSTI2022

Conformer-Specific Dissociation Dynamics in Dimethyl Methylphosphonate Radical Cation

The dynamics of the dimethyl methylphosphonate (DMMP) radical cation after production by strong field adiabatic ionization have been investigated. Pump-probe experiments using strong field 1300 nm pulses to adiabatically ionize DMMP and a 800 nm non-ionizing probe induce coherent oscillations of the parent ion yield with a period of about 45 fs. The yields of two fragments, PO 2 C 2 H 7 + and PO 2 CH 4 + , oscillate approximately out of phase with the parent ion, but with a slight phase shift relative to each other. We use electronic structure theory and nonadiabatic surface hopping dynamics to understand the underlying dynamics. The results show that while the cation oscillates on the ground state along the P=O bond stretch coordinate, the probe excites population to higher electronic states that can lead to fragments PO 2 C 2 H 7 + and PO 2 CH 4 + . The computational results combined with the experimental observations indicate that the two conformers of DMMP that are populated under experimental conditions exhibit different dynamics after being excited to the higher electronic states of the cation leading to different dissociation products. These results highlight the potential usefulness of these pump-probe measurements as a tool to study conformer-specific dynamics in molecules of biological interest.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Evaluating disease surveillance strategies for early outbreak detection in contact networks with varying community structure

Disease surveillance systems allow public health agencies to respond to emerging diseases before they become widespread. Developing such systems requires identifying optimal ways to monitor in the context of an epidemic outbreak; this problem is known as sensor selection. Contact networks represent the dynamics of interaction in a population and are used to model how a disease spreads in a population and to explore strategies of sensor selection. We evaluated five sensor selection strategies on their ability to provide an early warning of a COVID-like outbreak in synthetic contact networks encapsulated in four network scenarios. Three of these scenarios assessed different aspects of community structure. The fourth scenario employed a contact network representing the population and interactions of 6.8 million people in New York City, constructed from an agent-based simulation using census and transportation data. This scenario exemplifies how sensor selection strategies may perform in a real-world, urban context. Our findings suggest that the choice of the optimal strategy depends heavily on the community structure of the network. Strategies that select highly connected nodes or maximize network coverage are the optimal surveillance strategy for outbreak detection in many network community structures. However, a naive implementation of these strategies may fail to provide an early warning at all—including in the New York City scenario. Moreover, these methods are impractical for real-world use as they require knowledge of the underlying contact network. Instead, a selection strategy that starts with a set of random nodes and then performs a random walk through a chain of neighbors reliably provides early warnings without requiring prior knowledge of the network. We find this method, called “random chain”, to be the most pragmatic for implementation in a real-world disease surveillance context.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

The Stochastic Gause Predator-Prey model: Noise-induced extinctions and invariance

We consider the Gause predator-prey with general bounded or sub‑linear functional responses, – which includes those of Holling types Ⅰ–Ⅳ. – and multiplicative Gaussian noise. In contrast to previous studies, the prey in our model follows logistic dynamics while the predator's population is solely regulated by consumption of the prey. To ensure well-posedeness, we derive explicit Lyapunov‐type criteria ensuring global positivity and moment boundedness of solutions. We find conditions for noise‑induced extinctions, proving that stochasticity can drive either population to collapse even when the deterministic analogue predicts stable coexistence. In the case when the predator becomes extinct, we establish a limiting distribution for the predator's population. Last, for functional responses of Holling type Ⅰ, we provide sufficient conditions on the intensity of the noise for the existence and uniqueness of a stationary distribution.

Gause model↗
DOE OSTI2024

Efficient Scalable Contact Network Generation from Population Data

Modeling the contacts among a population is critical to understanding the dynamics of a disease outbreak. Contact networks, where nodes are individuals and edges are contacts among them, are used to represent these complex individual-level interactions. In this work, we are given the daily activity schedules of an urban population that represent the activity location and time of individuals in a population during a single twenty four hour period over multiple days. Using collocation to determine contact between individuals, our goal is to extract hourly contact networks from large-scale activity data. We improve upon the existing adjacency matrix-based method by implementing our custom sparse matrix multiplication algorithm. Starting with a Python implementation, we achieve a 1600x speed up in the computation with a fast custom designed sparse matrix multiplier algorithm implemented in the C++ language. This work is central to future parallel designs of the problem.

97 MATHEMATICS AND COMPUTING↗