Engineering PapersSearch

SEARCH · Engineering Papers

Results for “series”

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 37 records · Page 2

Systematic characterization of unknown compounds via dimensionality reduction of time series

Analysis of ambient aerosols provides valuable insight into particle sources and formation chemistry. However, due to the complexity of atmospheric data and the dynamic nature of aerosol composition, a substantial fraction of data often become discarded by conventional analysis methods. Furthermore, a large fraction of chemical species within those data are unidentifiable due to a lack of matching spectral information, resulting in suboptimal characterization of chemical composition. Previous work has demonstrated techniques for cataloging analytes in a chromatographic dataset by deconvolution of mass spectra, but integration of these analytes throughout a large dataset remains time consuming. Here, we present a method to automatically identify an ion for quantitation for single-ion chromatogram based peak fitting and integration, enabling comprehensive integration of analytes with minimal user interaction. The resulting time series are clustered with a machine-learning based dimensionality reduction technique to systematically investigate the underlying characteristics of the categorized analytes and gain new insights into the chemical composition and physicochemical properties of the unidentifiable analytes. We apply these methods to existing atmospheric datasets collected in Manacapuru, Brazil during the GoAmazon2014/5 campaign to identify new analytes and interpret their variability and transformations in the atmosphere. The analysis results generate 408 time series from cataloged analytes of interest, and the clustering of those time series with spherical k-means results in 8 distinct clusters. We find the analytes form clusters based on their distinct physicochemical properties, demonstrating the method’s ability to systematically identify and selectively filter contaminants and instrumental analytes and characterize the unidentifiable analytes.

54 ENVIRONMENTAL SCIENCES

NEUWAVE-12 at ESS: Workshop Series on Wavelength Dependent Neutron Imaging Back in Full Swing

The NEUWAVE workshop series, a cornerstone for advancements in energy- or wavelength-resolved neutron imaging, held its 12th meeting at the European Spallation Source (ESS) in Lund from September 1 to 5, 2024. The participants can be seen in Figure 1 during the site visit. The series focuses on developments in techniques using phenomena such as Bragg-edges, incoherent and inelastic scattering, or neutron absorption resonances as contrast mechanisms and is therefore far more specialized than more general conferences for neutron imaging such as the WCNR or ITMNR series [Citation1]. After a break during the pandemic and record attendance at NEUWAVE-11 in Tokyo in 2023, which underscored the immense interest in wavelength dependent neutron imaging, NEUWAVE returned to its roots with a focus on a smaller group of specialists meeting in a workshop format, allowing for in-depth discussions and the exchange of new ideas.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS

Circularity Futures Workshop Series: Summary Report

The aim of this report is to synthesize key feedback received from the three-part Circularity Futures workshop series held in Spring 2024. The workshop series was conducted by the National Renewable Energy Laboratory (NREL) on behalf of U.S. Department of Energy, Office Energy Efficiency and Renewable Energy (EERE), and was broken into three workshops: Workshop 1 - Circularity Analysis Needs and Priorities; Workshop 2 - Circularity Metrics and Indicators; and Workshop 3 - Circularity Data. Together, the workshops focused on identifying the existing priorities and gaps in the circularity modeling space, understanding different stakeholders' use and interpretation of circularity metrics and indicators, identifying common data gaps and data quality challenges, and assessing the robustness of available solutions. The workshop series brought a diverse group of stakeholders - including representatives from U.S. government offices, national labs, nonprofit organizations, industry, and academia - to collect first-hand feedback on needs, priorities, challenges and opportunities in the circularity modeling and analysis space. The workshop discussions highlighted numerous common needs, priorities and challenges among the interviewed groups. Several topics were frequently discussed, including: 1) Circularity as a pathway for sustainable economic growth: While circularity is generally defined in terms of resource conservation and reducing wasteful disposal of materials, participants agreed that circular strategies should serve broader economic, environmental, and social goals. It is therefore crucial for circularity analysis to look beyond waste reduction and instead evaluate a variety of impact metrics such as cost savings, job creation, air quality, and pollutant emissions. Mutli-criteria decision-making frameworks may be useful for making sense of disparate metrics and evaluating tradeoffs between impact categories.; 2) Economic and social factors are not well understood: Underdevelopment of existing end-of-life (EOL) management infrastructure, inconsistent standardization codes and policy space in reusing recycled content, and suboptimal collection and sorting strategies collectively contribute to uncertainty about the economic potential of circular pathways. The latter observation is consistent among all technologies but more emphasized for renewable energy systems. Social impacts of circularity practices are less understood and less researched than other sustainability aspects.; 3) Inconsistent methods for assessing emerging technologies: LCA and TEA results vary widely depending on the assumptions made with regards to market adoption of new technologies. Emerging technologies suffer limited availability of data needed to conduct a robust circularity analysis. Yet, understanding projected impacts of proposed nascent technology is a key need for different stakeholder groups.; and 4) Lack of temporally and geospatially explicit data: There is a need for open data that represents variations in circularity technologies over time and location. The lack thereof leads to aggregated and potentially misrepresented results in circularity analysis. Sensitivity analyses should be included to verify whether options perceived as more sustainable align with real-world practices.

29 ENERGY PLANNING, POLICY, AND ECONOMY

An interregional optimization approach for time series aggregation in continent-scale electricity system models

Modeling electric power systems with high shares of weather-dependent resources requires tradeoffs between temporal, spatial, and operational resolution. Many studies perform time series aggregation using clustering algorithms to reduce the temporal dimension, but when modeling continent-scale electricity systems that are large enough to contain multiple independent weather systems, this approach requires large numbers of representative periods to minimize errors in regional wind and solar capacity factors. Here, a new optimization-based approach for representative period selection and weighting is introduced that minimizes regional errors in average renewable capacity factors and electricity demand. The method delivers higher regional fidelity with fewer representative periods than alternative clustering methods when applied to wind, solar, and demand profiles for the contiguous United States. When representative periods are selected from multiple weather years, the optimized method reproduces regional averages with lower error than a complete 365-day time series from any single weather year. The method identifies only representative (as opposed to outlying) periods but can be combined with an iterative "stress period" identification approach to guide efficient decision-making considering both average and high-risk weather conditions.

24 POWER TRANSMISSION AND DISTRIBUTION

Is a Homologous Series Member Equal to or Greater than the Sum of its Subunits? A Focus on a Few Fe-based BaNiSn3 subunits: LnFeGe3 (Ln = La, Pr)

The homologous series Lnn+1MnGe3n+1 (Ln = Ce, Pr, M = Fe, Co) has proven to be a fruitful platform for tuning physical properties as a function of subunit stacking. In this study, we investigate the crystal growth, structure, and physical properties of PrFeGe3, aiming to deconstruct Pr4Fe3Ge10 (n = 3) into its subunits, including the BaNiSn3 and CeNiSi2 structure types. The electrical resistivity and magnetic susceptibility of PrFeGe3 were measured. The magnetic properties reveal PrFeGe3 to be antiferromagnetic (5 K) with an effective moment of eff = 4.38 B/F.U., larger than the spin-only moment from Pr3+ (3.58 B). Due to the large magnetic moment observed, the oxidation state of the elements was determined using spectroscopic methods including X-ray absorption spectroscopy and X-ray photoelectron spectroscopy. The neutral oxidation state of iron determined from X-ray photoelectron spectroscopy supports itinerant magnetic behavior of iron. Annealing, in tandem with differential scanning calorimetry, revealed PrFeGe3 to be a precursor for the formation of new homologous series members.

36 MATERIALS SCIENCE

Estimating an executive summary of a time series: the tendency

In this paper, we revisit the problem of decomposing a signal into a tendency and a residual. The tendency describes an executive summary of a signal that encapsulates its notable characteristics while disregarding seemingly random, less interesting aspects. Building upon the Intrinsic Time Decomposition (ITD) and information-theoretical analysis, we introduce two alternative procedures for selecting the tendency from the ITD baselines. The first is based on the maximum extrema prominence, namely the maximum difference between extrema within each baseline. Specifically this method selects the tendency as the baseline from which an ITD step would produce the largest decline of the maximum prominence. The second method uses the rotations from the ITD and selects the tendency as the last baseline for which the associated rotation is statistically stationary. We delve into a comparative analysis of the information content and interpretability of the tendencies obtained by our proposed methods and those obtained through conventional low-pass filtering schemes, particularly the Hodrik–Prescott (HP) filter. Our findings underscore a fundamental distinction in the nature and interpretability of these tendencies, highlighting their context-dependent utility with emphasis in multi-scale signals. Through a series of real-world applications, we demonstrate the computational robustness and practical utility of our proposed tendencies, emphasizing their adaptability and relevance in diverse time series contexts.

Time series analysis

Different Damp-Heat-Induced Series Resistance Degradation Behaviors in Fielded Crystalline Silicon Photovoltaic Modules Due to Difference in Bill of Materials

This case study investigates mono-crystalline silicon modules from underperforming portions of a utility-scale photovoltaic power plant. Field-collected I-V curves and electroluminescence imaging suggested that increased series resistance was a primary factor driving module degradation. Selected modules were removed from the field for further analysis, including incremental damp heat accelerated testing, which confirmed a progression in series resistance degradation. Two distinct cell degradation behaviors became apparent during the investigation. Cross-sectional scanning electron microscopy (with elemental analysis) and scanning spreading resistance microscopy identified key differences between the two degradation mechanisms, primarily grid finger width and contact resistance. Additionally, the study highlights the reliability implications of retest requirements in International Electrotechnical Commission 61215 for material changes and how they may have mitigated the degradation observed at this site.

14 SOLAR ENERGY

Refining Control, Charging, and Battery Chemistry for CO 2 e Savings in Heavy-Duty Off-Road Plug-In Series Hybrid

With current and future regulations continuing to drive reductions in carbon dioxide equivalent (CO 2 e) emissions in the on-road industry, the off-road industry is also likely to be regulated for fuel and CO 2 e savings. This work focuses on converting a heavy-duty off-road material handler from a conventional diesel powertrain to a plug-in series hybrid, achieving a 49% fuel reduction and 29% CO 2 e reduction via simulation. Control strategies were refined for energy savings, including a regenerative braking strategy to increase regenerative braking and a load-following hydraulic strategy to decrease electrical energy consumption. The load-following hydraulic control shuts off the hydraulic electric machine when it is not needed—an approach not previously seen in a load-sensing, pressure-compensated system. Furthermore, these strategies achieved a 24.1% fuel savings, resulting in total savings of 61% in fuel and 41% in CO 2 e in the plug-in series compared to the conventional machine. Beyond control strategies, this study evaluated battery chemistry and charging strategy refinements for total cost of ownership (TCO) and lifetime CO 2 e. LFP batteries emerged as the most cost-effective and least emitting due to their longer lifespan, which reduced replacement frequency. Charging comparisons showed that Level 2 charging (L2C) typically resulted in lower TCO but higher lifetime CO 2 e than DC fast charging (DCFC). DCFC costs were heavily influenced by local demand charges, and DCFC emissions were heavily influenced by local grid emissions.

33 ADVANCED PROPULSION SYSTEMS

Accidental symmetries, Hilbert series, and friends

Accidental symmetries in effective field theories can be established by computing and comparing Hilbert series. This invites us to study them with the tools of invariant theory. Applying this technology, we spotlight three classes of accidental symmetries that hold to all orders for non-derivative interactions. They are broken by derivative interactions and become ordinary finite-order accidental symmetries. To systematically understand the origin and the patterns of accidental symmetries, we introduce a novel mathematical construct — a (non-transitive) binary relation between subgroups that we call friendship. Equipped with this, we derive new criteria for all-order accidental symmetries in terms of friends, and criteria for finite-order accidental symmetries in terms of friends ma non troppo. They allow us to verify and identify accidental symmetries more efficiently without computing the Hilbert series. We demonstrate the success of our new criteria by applying them to a variety of sample accidental symmetries, including the custodial symmetry in the Higgs sector of the Standard Model effective field theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS

Stimulating β -Series Precipitation in Mg–Nd Alloys Via Microalloying: A Comparison of Electron Microscopy and Small-Angle Scattering Techniques

The Mg–Nd alloy system is of particular interest in the development of high-strength, lightweight structural alloys due to the formation of strengthening metastable Mg–Nd β-series precipitates during heat treatment. The strengthening is primarily provided by a combination of the β''' and β 1 precipitation. It is therefore important to understand how the precipitation behavior can be enhanced by other common alloying elements. In this work, the effects of 0.2 wt pct Zn and Ca on β-series precipitation were studied. Small-angle/ultra-small-angle X-ray scattering (SAXS/USAXS) techniques in combination with scanning transmission electron microscopy (STEM) were used to study the evolution of precipitation microstructure. Here, it is found that the Zn additions refine the precipitates, leading to an increase in hardness with aging at 200 °C. On the other hand, the Ca additions result in an acceleration in the formation of larger β 1 precipitates and chains which provides less strengthening. The β 1 chains are surrounded by precipitate-free zones (PFZs) that further contribute to the decreases in hardness observed in the over-aged condition. This paper demonstrated that SAXS/USAXS provides a powerful tool for an in situ study of the early stages of precipitation in the Mg–Nd-based alloys.

36 MATERIALS SCIENCE

Scaling of energy delivered through an electrostatic discharge to a small series load

We study the energy delivered through a small-resistance series “victim” load during electrostatic discharge events in air. For gap lengths over 1 mm, the fraction of the stored energy delivered is mostly gap-length independent, with a slight decrease at larger gaps due to electrode geometry. The energy to the victim scales linearly with circuit capacitance and victim load resistance but is not strongly dependent on circuit inductance. This scaling leads to a simple approach to predicting the maximum energy that will be delivered to a series resistance for the case where the victim load resistance is lower than the spark resistance.

42 ENGINEERING

Consistent thermodynamic properties for alicyclic components of jet fuels: Experimental data, estimation methods, and homologous series trends

Alkylcycloalkanes represent a significant fraction of jet fuel components. An evaluation of their thermodynamic properties, enthalpies of formation in liquid and gas phases and enthalpies of vaporization, was conducted. A combination of available experimental data, up-to-date group-contribution methods, high-level quantum-chemical calculations, and homologous series trends was used to identify outliers and to recommend the most reliable values. The group-contribution approach was found to work well for the enthalpies of vaporization. Its performance for the enthalpies of formation in the liquid and gas phases was found to be substantially less effective, especially considering notable differences in this property among stereoisomers. Computationally affordable high-level ab initio results and homologous series trend analysis appeared more reliable. In conclusion, the recommended property values for 212 individual compounds and their isomeric mixtures were provided.

09 BIOMASS FUELS

Study on the magnetothermal properties of the D y 1 - x T b x A l 2 series of compounds

Here, in this work, we developed a theoretical model Hamiltonian, in the mean field approximation, to describe the magnetic and magnetocaloric behavior of the series of compounds Dy 1-x Tb x Al 2 (x = 0.00, 0.15, 0.25, 0.30, 0.35, 0.40 and 0.75). We adjusted the exchange parameters λ DyDy , λ TbTb , and λ DyTb to obtain the spin reorientation temperatures (T SR ) and the critical temperature (T C ) for each compound in the series. The results obtained by the Hamiltonian model agree satisfactorily with the experimental results. The heat capacity curves with and without an applied magnetic field, adiabatic temperature variation and isothermal entropy variation were modeled and compared with the experimental data. As the experimental results show, our model was also able to reproduce the change in the spin reorientation process: a first order spin reorientation transition appears for concentrations x = 0.15, 0.25, 0.30, 0.35, and no spin reorientation transitions after x = 0.40.

36 MATERIALS SCIENCE

Integrating very-high-resolution imagery, Sentinel-2 time-series data, and machine learning to map shrub fractional abundance across arid and semi-arid ecosystems in China

Shrub fractional abundance (SFA), the proportion of shrub cover per unit area, serves as a critical indicator of environmental aridity and ecosystem health in arid and semi-arid regions, particularly across the Mongolian steppe. However, large-scale SFA mapping in Mongolian steppe ecosystems remains challenging due to the small crown size of shrubs, their sparse distribution, and spectral overlap with coexisting low vegetation (e.g., grasses and herbs), which hinders accurate detection using coarser-resolution satellite data or traditional field surveys. To address these challenges, we developed a two-step approach that integrates very-high-resolution (VHR) imagery, time-series Sentinel-2 data, and deep learning techniques. First, we generated high-accuracy benchmark maps of individual shrub crowns from 0.5 m VHR imagery by combining manual segmentation with a hybrid deep learning framework (Dino V2 and convolutional neural networks). Second, we used these shrub crown maps as training data to build an XGBoost model for predicting SFA from 20 m Sentinel-2 time-series data, leveraging phenological information to improve estimation. We validated our approach across 70 sites (1km 2 each) in the Inner Mongolia Autonomous Region, which is representative of Mongolian steppe ecosystems. From VHR imagery, we mapped 1.31 million shrub crowns with an accuracy of R 2 = 0.92. Scaling up with Sentinel-2 data yielded regional SFA maps with an R 2 = 0.60. Further SHAP (SHapley Additive exPlanations) analysis on the developed XGBoost model revealed that phenological metrics (particularly observations in early-May, mid-July, and late-September), which distinguish shrub phenology from that of other land cover types (e.g., grasses and bare soil), were the most influential predictors of SFA. Finally, our regional SFA maps uncovered unimodal relationships between shrub distribution and climate variables, peaking at mean annual minimum temperatures near 0 °C and annual precipitation around 200 mm. Collectively, these findings demonstrate how the integration of multi-source remote sensing and machine learning can overcome historical limitations in SFA mapping, enabling accurate, spatially continuous assessments across vast Inner-Mongolian steppe ecosystems. Our framework has the potential to be applied to other steppe ecosystems and dryland ecosystems across the Mongolian steppe and beyond, offering a foundation for improved monitoring and ecological impact assessments in the face of global climate changes.

Arid and semi-arid landscapes

Near-IR Luminescence Tuning in a Series of Chalcogenophene Carboxylate-Decorated Neodymium Dimers

Here, the solvothermal synthesis of a series of Nd dimers decorated with various chalcogenophene carboxylates and 2,2':6',2"-terpyridine of the general formula, [Nd 2 (µ-XC 5 H 3 O 2 ) 2 (XC 5 H 3 O 2 ) 4 (N 3 C 15 H 11 ) 2 (H 2 O) 2 ] where X = O, S, Se, and Te, is reported. The solid-state structures were characterized using single-crystal X-ray diffraction (scXRD) and all the complexes are isomorphous, despite substitution of the heterocyclic chalcogen; phase purity was confirmed via powder X-ray diffraction (pXRD). Vibrational spectroscopy was collected and correlations between chalcogen identity and the binding strength of the carboxylate groups of the chalcogenophene ligands with each metal center were shown to be independent of chalcogen identity. All four complexes displayed Nd(III)-based near-IR luminescence and exhibited ligand-sensitized emission. Varying the chalcogenophene chromophore enabled tuning of the sensitizing triplet state energy level, as evidenced by an 8-fold increase in the sensitization of the TeCA-decorated dimer relative to the other chalcogenophene congeners. This behavior was rationalized by comparing the Nd(III) acceptor and ligand donor states across the series. The donor triplet state of each ligand was estimated via low-temperature (77 K) phosphorescence measurements from 1:1 mixtures with Gd(III); these were found to be 24,631 cm –1 for furan-2-carboxylic acid (FCA), 23,764 cm –1 for thiophene-2-carboxylic acid (TCA), 22,548 cm –1 for selenophene-2-carboxylic acid (SeCA), and 21,186 cm –1 for tellurophene-2-carboxylic acid (TeCA). The greater sensitization efficiency of TeCA is the result of well-matched ligand donor and metal acceptor levels and thus suppression of nonradiative back-energy transfer. More broadly, triplet energy level information for these ligands serves as a guide for future application to other target metals based on the electronic properties necessary to effect efficient sensitization.

Coordination Chemistry

Dimensional Evolution Guides Property Control in the A n Cu 4– n TiS 4 Semiconductor Series

Through progressive reduction of the three-dimensional (3D) covalent network of Cu 4 TiS 4 , we isolate seven new members of the A n Cu 4–n TiS 4 family (A = alkali metal; n = 0–4), spanning 3D, 2D, 1D, and 0D structural fragments. The dimensional reduction is rational, as it preserves the edge-sharing connectivity between [CuS 4 ] 7– and [TiS 4 ] 4– tetrahedra across the series. This structural evolution is driven by the stepwise substitution of Cu with alkali metals, guiding the formation of fragments with reduced dimensionality. The effects of “n” and “A” on the crystal structures, stabilities, electronic structures, and optoelectronic properties are profound, demonstrating that the manipulation of alkali metal size and A n Cu 4–n TiS 4 stoichiometry enables predictable variations in structure and properties. For example, the n = 0 and n = 4 end members of the A n Cu 4–n TiS 4 family set the range of achievable band gaps with 2.00 eV for Cu 4 TiS 4 , 2.60 eV for Na 4 TiS 4 , and intermediate values for the n = 1–3 members. Notably, CsCu 3 TiS 4 exhibits exceptional air stability and congruent melting, with density functional theory (DFT) calculating moderate hole and electron effective masses in specific crystallographic directions (mh = 1.24m 0 , me = 0.87m 0 ). Additionally, A 3 CuTiS 4 (A = Na, K, Rb) displays direct band gap behavior and long photoluminescence lifetimes of 2.3–8.6 μs, and K 3 CuTiS 4 has a PLQY of 5.19%. These findings underscore the potential of the A n Cu 4–n TiS 4 family for applications in optoelectronics and demonstrate widely applicable design concepts that unveil rational stoichiometries within a given composition space to generate a series of crystal structures related through an evolving covalent dimensionality that corresponds to a predictable electronic structure and property progression.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

ATAT: Astronomical Transformer for time series and Tabular data

Context. The advent of next-generation survey instruments, such as theVera C. RubinObservatory and its Legacy Survey of Space and Time (LSST), is opening a window for new research in time-domain astronomy. The Extended LSST Astronomical Time-Series Classification Challenge (ELAsTiCC) was created to test the capacity of brokers to deal with a simulated LSST stream. Aims. Our aim is to develop a next-generation model for the classification of variable astronomical objects. We describe ATAT, the Astronomical Transformer for time series And Tabular data, a classification model conceived by the ALeRCE alert broker to classify light curves from next-generation alert streams. ATAT was tested in production during the first round of the ELAsTiCC campaigns. Methods. ATAT consists of two transformer models that encode light curves and features using novel time modulation and quantile feature tokenizer mechanisms, respectively. ATAT was trained on different combinations of light curves, metadata, and features calculated over the light curves. We compare ATAT against the current ALeRCE classifier, a balanced hierarchical random forest (BHRF) trained on human-engineered features derived from light curves and metadata. Results. When trained on light curves and metadata, ATAT achieves a macro F1 score of 82.9 ± 0.4 in 20 classes, outperforming the BHRF model trained on 429 features, which achieves a macro F1 score of 79.4 ± 0.1. Conclusions. The use of transformer multimodal architectures, combining light curves and tabular data, opens new possibilities for classifying alerts from a new generation of large etendue telescopes, such as theVera C. RubinObservatory, in real-world brokering scenarios.

Astronomy & Astrophysics

Uncovering heterogeneous intercommunity disease transmission from neutral allele frequency time series

The COVID-19 pandemic has underscored the need for accurate epidemic forecasting to predict pathogen spread, evolution, and evaluate intervention strategies. Forecast reliability hinges on detailed knowledge of disease transmission across population segments, which may be inferred from contact surveys or mobility data. However, these indirect approaches make it difficult to estimate rare transmissions between socially or geographically distant communities. We show that the steep ramp-up of genome sequencing surveillance during the pandemic can be leveraged to directly identify transmission patterns between geographically defined communities. Our approach uses a hidden Markov model to infer the fraction of infections a community imports from others based on how rapidly allele frequencies in the focal community converge to those in the donor communities. Applying this method to SARS-CoV-2 sequencing data from England and the United States, we uncover networks of intercommunity transmission that reflect geographical relationships while exposing significant long-range interactions. The scaling of importation rate with distance is consistent across both countries, yet weaker than expected based on mobility data, highlighting limitations of indirect inference. We show that transmission patterns can change between waves of variants of concern and analyze how the inferred heterogeneity in intercommunity transmission impacts evolutionary forecasts. While applied here to geographically defined communities, our approach could be applied to those defined by other traits (e.g., age, socioeconomic status), provided time-series data can be stratified accordingly. Overall, our study highlights population genomic time series data as a crucial record of epidemiological interactions, which can be deciphered using tree-free inference methods.

Okada, Takashi [Department of Physics; University