Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic framework”

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 1,369 records · Page 76

Signatures of Thermoreversible Associations in X-ray and Neutron Scattering from Dilute Polyzwitterion Solutions

In aqueous solutions of polyzwitterions (PZs), an interplay between dipole–dipole interactions and hydration of zwitterionic groups can lead to thermoreversible associations, which have been difficult to detect in experiments. Here, in this study, we investigated dilute aqueous solutions of poly(1-(3-sulphonatopropyl)-2-vinylpyridinium) (P2VPPS) using small-angle X-ray and neutron scattering (SAXS and SANS) to probe the structure and neutron spin-echo (NSE) spectroscopy to probe dynamics. The SAXS and SANS data show that the correlation length increases with an increase in the concentration of P2VPPS for three different molecular weights. The addition of 0.1 M NaCl to one of the solutions led to almost no dependence of the correlation length on the concentration. Such a concentration dependence of the correlation length suggests the formation of clusters driven by thermoreversible associations in the solutions. The NSE measurements show that the solutions with a larger correlation length display slower relaxation, reflecting the reduced mobility of larger clusters. These results should be considered as signatures of thermoreversible associations in the solutions of P2VPPS. To establish a quantitative link between local structure (clusters) and dynamics in dilute solutions of P2VPPS, we combined a thermoreversible gelation theory for the structure of PZ solutions (Li, S.-F.; Muthukumar, M. Theory of Thermoreversible Gelation and Anomalous Concentration Fluctuations in Polyzwitterion Solutions. J. Chem. Phys. 2024, 161, 024903) with a model for the dynamics of the clusters (generalized Zimm model), developed in this work. Using such a theoretical framework, we have predicted the distribution of clusters in the solutions probed with SANS and SAXS. With the distributions, the generalized Zimm model has been used to extract the diffusion constant of the clusters and their characteristic size from the NSE data, where the latter agrees with the values estimated from the SAXS/SANS data. These findings confirm the presence of thermoreversible associations and establish a quantitative link between local structure (clusters) and dynamics in dilute solutions of P2VPPS. Furthermore, with growing interest in technological applications, this work can provide useful insights into the structural and dynamical properties of other PZ solutions.

field theory↗

Extended Lagrangian Born–Oppenheimer molecular dynamics for orbital-free density-functional theory and polarizable charge equilibration models

We report extended Lagrangian Born–Oppenheimer molecular dynamics (XL-BOMD) is formulated for orbital-free Hohenberg–Kohn density-functional theory and for charge equilibration and polarizable force-field models that can be derived from the same orbital-free framework. The purpose is to introduce the most recent features of orbital-based XL-BOMD to molecular dynamics simulations based on charge equilibration and polarizable force-field models. These features include a metric tensor generalization of the extended harmonic potential, preconditioners, and the ability to use only a single Coulomb summation to determine the fully equilibrated charges and the interatomic forces in each time step for the shadow Born–Oppenheimer potential energy surface. The orbital-free formulation has a charge-dependent, short-range energy term that is separate from long-range Coulomb interactions. This enables local parameterizations of the short-range energy term, while the long-range electrostatic interactions can be treated separately. The theory is illustrated for molecular dynamics simulations of an atomistic system described by a charge equilibration model with periodic boundary conditions. The system of linear equations that determines the equilibrated charges and the forces is diagonal, and only a single Ewald summation is needed in each time step. The simulations exhibit the same features in accuracy, convergence, and stability as are expected from orbital-based XL-BOMD.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Boom Minimization Framework for Supersonic Aircraft Using CFD Analysis

A new framework is presented for shape optimization using analytical shape functions and high-fidelity computational fluid dynamics (CFD) via Cart3D. The focus of the paper is the system-level integration of several key enabling analysis tools and automation methods to perform shape optimization and reduce sonic boom footprint. A boom mitigation case study subject to performance, stability and geometrical requirements is presented to demonstrate a subset of the capabilities of the framework. Lastly, a design space exploration is carried out to assess the key parameters and constraints driving the design.

Ordaz, Irian↗

Surfactant-Specific AI-Driven Molecular Design: Integrating Generative Models, Predictive Modeling, and Reinforcement Learning for Tailored Surfactant Synthesis

Molecular design is a critical aspect of various scientific and industrial fields, where the properties of molecules hold significant importance. In this study, a 3-fold methodology design is presented that leverages the power of generative artificial intelligence (AI), predictive modeling, and reinforcement learning to create tailored molecules with desired properties. This model synergistically combines deep learning techniques with Self-Referencing Embedded Strings (SELFIES) molecular representation to build a generative model that generates valid molecules and a graphical neural network model that accurately forecasts molecular properties. The Variational Autoencoder (VAE) coupled with reinforcement learning helps refine molecule generation based on targeted attributes. Data from an experimental study involving surfactants were used to test the framework. A validation of the structural integrity of the molecules generated was conducted, and Tanimoto similarities were used to quantify the similarity and diversity between the original and generated molecular structures. Also, saliency maps for the generated surfactants were produced to identify the features explaining the property values. Lastly, molecular dynamics simulations were used to validate the stability of the generated molecules. The results showed that the proposed framework can effectively produce valid molecules within the set property threshold value.

36 MATERIALS SCIENCE↗

Enhanced Sulfur Dioxide Adsorption in UiO-66 Through Crystal Engineering and Chalcogen Bonding

Adsorption of corrosive SO 2 gas occurs in metal–organic frameworks (MOFs) including UiO-66. Improvements in SO 2 capacity is obtained through the incorporation of residual modulators in the UiO-66 framework by introducing new binding sites in the material, through residual modulators. Four residual modulators were investigated (acetic acid, trifluoroacetic acid, 3-DMAP acid, cyanoacetic acid), and the UiO-66 framework modulated with cyanoacetic acid exhibited nearly twice the SO 2 uptake for the 18:1 modulator/linker synthesis ratio compared with other modulated UiO-66 structures. Density functional theory investigations confirmed that targeted host–guest interactions were maintained after the modulator was incorporated into the framework. The strongest binding energy was between SO 2 and cyanoacetic acid, consistent with dynamic SO 2 adsorption data, and identified contributions from both the SO 2 reacting with the residual modulator and the coordinating linkers. The successful increase in dynamic SO 2 capacity illustrates how often-overlooked non-covalent interactions can be used in targeted adsorption applications. Further investigation into weak electrostatic interactions for adsorption properties is also needed to advance the potential selectivity and capacity in the adsorption sphere.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Optimal placement of excitations and sensors for verification of large dynamical systems

The computationally difficult problem of the optimal placement of excitations and sensors to maximize the observed measurements is studied within the framework of combinatorial optimization, and is solved numerically using a variation of the simulated annealing heuristic algorithm. Results of numerical experiments including a square plate and a 960 degrees-of-freedom Control of Flexible Structure (COFS) truss structure, are presented. Though the algorithm produces suboptimal solutions, its generality and simplicity allow the treatment of complex dynamical systems which would otherwise be difficult to handle.

Salama, M.↗

Dynamically Unveiling Metal–Nitrogen Coordination during Thermal Activation to Design High-Efficient Atomically Dispersed CoN 4 Active Sites

We elucidate the structural evolution of CoN 4 sites during thermal activation by developing a zeolitic imidazolate framework (ZIF)-8-derived carbon host as an ideal model for Co 2+ ion adsorption. Subsequent in situ X-ray absorption spectroscopy analysis can dynamically track the conversion from inactive Co-OH and Co-O species into active CoN 4 sites. The critical transition occurs at 700 °C and becomes optimal at 900 °C, generating the highest intrinsic activity and four-electron selectivity for the oxygen reduction reaction (ORR). DFT calculations elucidate that the ORR is kinetically favored by the thermal-induced compressive strain of Co-N bonds in CoN 4 active sites formed at 900 °C. Further, we developed a two-step (i.e., Co ion doping and adsorption) Co-N-C catalyst with increased CoN 4 site density and optimized porosity for mass transport, and demonstrated its outstanding fuel cell performance and durability.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A computational modeling framework for pre-clinical evaluation of cardiac mapping systems

There are a variety of difficulties in evaluating clinical cardiac mapping systems, most notably the inability to record the transmembrane potential throughout the entire heart during patient procedures which prevents the comparison to a relevant “gold standard”. Cardiac mapping systems are comprised of hardware and software elements including sophisticated mathematical algorithms, both of which continue to undergo rapid innovation. The purpose of this study is to develop a computational modeling framework to evaluate the performance of cardiac mapping systems. The framework enables rigorous evaluation of a mapping system’s ability to localize and characterize (i.e., focal or reentrant) arrhythmogenic sources in the heart. The main component of our tool is a library of computer simulations of various dynamic patterns throughout the entire heart in which the type and location of the arrhythmogenic sources are known. Our framework allows for performance evaluation for various electrode configurations, heart geometries, arrhythmias, and electrogram noise levels and involves blind comparison of mapping systems against a “silver standard” comprised of computer simulations in which the precise transmembrane potential patterns throughout the heart are known. A feasibility study was performed using simulations of patterns in the human left atria and three hypothetical virtual catheter electrode arrays. Activation times (AcT) and patterns (AcP) were computed for three virtual electrode arrays: two basket arrays with good and poor contact and one high-resolution grid with uniform spacing. The average root mean squared difference of AcTs of electrograms and those of the nearest endocardial action potential was less than 1 ms and therefore appears to be a poor performance metric. In an effort to standardize performance evaluation of mapping systems a novel performance metric is introduced based on the number of AcPs identified correctly and those considered spurious as well as misclassifications of arrhythmia type; spatial and temporal localization accuracy of correctly identified patterns was also quantified. This approach provides a rigorous quantitative analysis of cardiac mapping system performance. Proof of concept of this computational evaluation framework suggests that it could help safeguard that mapping systems perform as expected as well as provide estimates of system accuracy.

59 BASIC BIOLOGICAL SCIENCES↗

Empirical STORM-E Model

Auroral nighttime infrared emission observed by the Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) instrument onboard the Thermosphere-Ionosphere-Mesosphere Energetics and Dynamics (TIMED) satellite is used to develop an empirical model of geomagnetic storm enhancements to E-region peak electron densities. The empirical model is called STORM-E and will be incorporated into the 2012 release of the International Reference Ionosphere (IRI). The proxy for characterizing the E-region response to geomagnetic forcing is NO+(v) volume emission rates (VER) derived from the TIMED/SABER 4.3 lm channel limb radiance measurements. The storm-time response of the NO+(v) 4.3 lm VER is sensitive to auroral particle precipitation. A statistical database of storm-time to climatological quiet-time ratios of SABER-observed NO+(v) 4.3 lm VER are fit to widely available geomagnetic indices using the theoretical framework of linear impulse-response theory. The STORM-E model provides a dynamic storm-time correction factor to adjust a known quiescent E-region electron density peak concentration for geomagnetic enhancements due to auroral particle precipitation. Part II of this series describes the explicit development of the empirical storm-time correction factor for E-region peak electron densities, and shows comparisons of E-region electron densities between STORM-E predictions and incoherent scatter radar measurements. In this paper, Part I of the series, the efficacy of using SABER-derived NO+(v) VER as a proxy for the E-region response to solar-geomagnetic disturbances is presented. Furthermore, a detailed description of the algorithms and methodologies used to derive NO+(v) VER from SABER 4.3 lm limb emission measurements is given. Finally, an assessment of key uncertainties in retrieving NO+(v) VER is presented

Mertens, Christopher J.↗

RAVEN User Manual

RAVEN is a generic software framework to perform parametric and probabilistic analysis based on the response of complex system codes. The initial development was aimed to provide dynamic risk analysis capabilities to the Thermo-Hydraulic code RELAP-7, currently under development at the Idaho National Laboratory (INL). Although the initial goal has been fully accomplished, RAVEN is now a multi-purpose probabilistic and uncertainty quantification platform, capable to agnostically communicate with any system code. This agnosticism includes providing Application Programming Interfaces (APIs). These APIs are used to allow RAVEN to interact with any code as long as all the parameters that need to be perturbed are accessible by inputs files or via python interfaces. RAVEN is capable of investigating the system response, and investigating the input space using Monte Carlo, Grid, or Latin Hyper Cube sampling schemes, but its strength is focused to- ward system feature discovery, such as limit surfaces, separating regions of the input space leading to system failure, using dynamic supervised learning techniques. The development of RAVEN has started in 2012, when, within the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, the need to provide a modern risk evaluation framework became stronger. RAVEN principal assignment is to provide the necessary software and algorithms in order to employ the concept developed by the Risk Informed Safety Margin Characterization (RISMC) program. RISMC is one of the pathways defined within the Light Water Reactor Sustainability (LWRS) program. In the RISMC approach, the goal is not just the individuation of the frequency of an event potentially leading to a system failure, but the closeness (or not) to key safety-related events. Hence, the approach is interested in identifying and increasing the safety margins related to those events. A safety margin is a numerical value quantifying the probability that a safety metric (e.g. for an important process such as peak pressure in a pipe) is exceeded under certain conditions. The initial development of RAVEN has been focused on providing dynamic risk assessment capability to RELAP-7, currently under development at the INL and, likely, future replacement of the RELAP5-3D code. Most the capabilities that have been implemented having RELAP-7 as principal focus are easily deployable for other system codes. For this reason, several side activates are currently ongoing for coupling RAVEN with soft- ware such as RELAP5-3D, etc. The aim of this document is the explanation of the input requirements, focalizing on the input structure.

97 MATHEMATICS AND COMPUTING↗

Final Technical Report

The capture of CO2 and its simultaneously conversion to useful chemical fuels driven by solar energy represents one of the best solutions to resolve our growing energy and environmental concerns. The most critical challenge to this endeavor is the rational design of a photocatalytic architecture that can effectively couple a given photosensitizer (PS) with an appropriate catalyst, thereby enabling efficient photosensitization of a multi-electron reduction catalysis. This research program aims to address this challenge using an interdisciplinary approach that combines innovative material design and synthesis, fundamental mechanistic studies, and photocatalytic performance evaluation. The strategies include 1) constructing and investigating a novel class of 2D COF hybrid photocatalysts with an effective photoactive organic building block as PS and a precisely incorporated CO2 reduction molecular catalyst (MC); and 2) mechanistic origins of CO2 photoreduction using a set of complementary time-resolved and in situ spectroscopic techniques. The novelty of the proposed hybrid system lies in the unprecedented combination of the unique advantage of porous crystalline COF PS with the precise catalytic function of MC for photocatalytic CO2 reduction. In the periods of the support (09/01/2019-12/31/2022), we have made research progress in four projects: 1) Exploring 2D COFs with incorporated Mn complex for light driven CO2 reduction; 2) The dependence of excited state and charge transfer dynamics on monomer structure of 2D COFs; and 3) Control over Charge Separation by Imine Structural Isomerization in Covalent Organic Frameworks with Implications on CO2 Photoreduction; and 4) The impact of monomer structure on the photoluminescence properties of COFs. We found that both monomer structure and linker chemistry can effectively impact the excited state dynamics, charge transfer properties, and photoluminescence quantum yields, the important properties that dictate their applications in photocatalysis. In addition, we found that the direction of imine linker determines charge transfer direction and thus controls the types of catalytic reactions (e.g. water oxidation or CO2 reduction reactions). The result from these fundamental studies provides important information for correlating the structure of the COF photocatalysts with their photophysical properties and catalytic functions, paving the way for their novel application in photocatalysis. We expect that our findings will contribute to addressing current shortcomings of semiconductor- and molecular-based photocatalytic systems that suffer from inefficient light harvesting and charge separation and poor adsorption and activation of reactants. In turn, this research will contribute towards the development of novel photocatalytic systems for CO2 reduction to generate renewable chemical fuels and simultaneously address the problem of mitigating climate change due to CO2 accumulation. In addition, the experimental approaches employed in this research can be easily transferred to other energy technologies and are expected to broadly impact fields involving photocatalysis, optoelectronic devices, and solar energy conversion. The proposed research has also been integrated with educational activities and serve as a basis to raise awareness around the critical issues of global energy production and consumption, and to develop the next generation of solar energy scientists.

14 SOLAR ENERGY↗

Coupled Monte Carlo Transport and Conjugate Heat Transfer for Wire-Wrapped Bundles Within the MOOSE Framework

Cardinal is an open-source application that couples OpenMC Monte Carlo transport and NekRS computational fluid dynamics (CFD) to the Multiphysics Object-Oriented Simulation Environment (MOOSE), closing neutronics and thermal-fluid gaps in conducting high-resolution multiscale and multiphysics analyses of nuclear systems. We first provide a brief introduction to Cardinal's software design, data mapping, and coupling strategy to highlight our approach to overcoming common challenges in high-fidelity multiphysics simulations. Here we then present two Cardinal simulations for hexagonal pin bundles. The first is a validation of Cardinal's conjugate heat transfer coupling of NekRS's Reynolds-Averaged Navier Stokes model with MOOSE's heat conduction physics for a bare seven-pin Freon-12 bundle flow experiment. Predictions for pin surface temperatures under three different heating modes agree reasonably well with experimental data and similar CFD modeling from the literature. The second simulation is a multiphysics coupling of OpenMC, NekRS, and BISON for a reduced-scale, seven-pin wire-wrapped version of an Advanced Burner Reactor bundle. Wire wraps are approximated using a momentum source model, and coupled predictions are provided for velocity, temperature, and power distribution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Understanding Air Transportation Market Dynamics Using a Search Algorithm for Calibrating Travel Demand and Price

This paper presents a search algorithm based framework to calibrate origin-destination (O-D) market specific airline ticket demands and prices for the Air Transportation System (ATS). This framework is used for calibrating an agent based model of the air ticket buy-sell process - Airline Evolutionary Simulation (Airline EVOS) -that has fidelity of detail that accounts for airline and consumer behaviors and the interdependencies they share between themselves and the NAS. More specificially, this algorithm simultaneous calibrates demand and airfares for each O-D market, to within specified threshold of a pre-specified target value. The proposed algorithm is illustrated with market data targets provided by the Transportation System Analysis Model (TSAM) and Airline Origin and Destination Survey (DB1B). Although we specify these models and datasources for this calibration exercise, the methods described in this paper are applicable to calibrating any low-level model of the ATS to some other demand forecast model-based data. We argue that using a calibration algorithm such as the one we present here to synchronize ATS models with specialized forecast demand models, is a powerful tool for establishing credible baseline conditions in experiments analyzing the effects of proposed policy changes to the ATS.

Kumar, Vivek↗

NeuroSEM: A hybrid framework for simulating multiphysics problems by coupling PINNs and spectral elements

Multiphysics problems that are characterized by complex interactions among fluid dynamics, heat transfer, structural mechanics, and electromagnetics, are inherently challenging due to their coupled nature. While experimental data on certain state variables may be available, integrating these data with numerical solvers remains a significant challenge. Physics-informed neural networks (PINNs) have shown promising results in various engineering disciplines, particularly in handling noisy data and solving inverse problems in partial differential equations (PDEs). However, their effectiveness in forecasting nonlinear phenomena in multiphysics regimes, particularly involving turbulence, is yet to be fully established. Here, this study introduces NeuroSEM, a hybrid framework integrating PINNs with the highfidelity Spectral Element Method (SEM) solver, Nektar++. NeuroSEM leverages the strengths of both PINNs and SEM, providing robust solutions for multiphysics problems. PINNs are trained to assimilate data and model physical phenomena in specific subdomains, which are then integrated into the Nektar++ solver. We demonstrate the efficiency and accuracy of NeuroSEM for thermal convection in cavity flow and flow past a cylinder. The framework effectively handles data assimilation by addressing those subdomains and state variables where the data is available. We applied NeuroSEM to the Rayleigh-B´enard convection system, including cases with missing thermal boundary conditions and noisy datasets. Finally, we applied the proposed NeuroSEM framework to real particle image velocimetry (PIV) data to capture flow patterns characterized by horseshoe vortical structures. Our results indicate that NeuroSEM accurately models the physical phenomena and assimilates the data within the specified subdomains. The framework’s plug-and-play nature facilitates its extension to other multiphysics or multiscale problems. Furthermore, NeuroSEM is optimized for efficient execution on emerging integrated GPU-CPU architectures. This hybrid approach enhances the accuracy and efficiency of simulations, making it a powerful tool for tackling complex engineering challenges in various scientific domains.

42 ENGINEERING↗

CompLaB v1.0: a scalable pore-scale model for flow, biogeochemistry, microbial metabolism, and biofilm dynamics

Abstract. Microbial activity and chemical reactions in porous media depend on the local conditions at the pore scale and can involve complex feedback with fluid flow and mass transport. We present a modeling framework that quantitatively accounts for the interactions between the bio(geo)chemical and physical processes and that can integrate genome-scale microbial metabolic information into a dynamically changing, spatially explicit representation of environmental conditions. The model couples a lattice Boltzmann implementation of Navier–Stokes (flow) and advection–diffusion-reaction (mass conservation) equations. Reaction formulations can include both kinetic rate expressions and flux balance analysis, thereby integrating reactive transport modeling and systems biology. We also show that the use of surrogate models such as neural network representations of in silico cell models can speed up computations significantly, facilitating applications to complex environmental systems. Parallelization enables simulations that resolve heterogeneity at multiple scales, and a cellular automaton module provides additional capabilities to simulate biofilm dynamics. The code thus constitutes a platform suitable for a range of environmental, engineering and – potentially – medical applications, in particular ones that involve the simulation of microbial dynamics.

58 GEOSCIENCES↗

ADEPT: A Pedagogical Framework for Integrating Agentic AI with Deterministic Scientific Workflows

The integration of Large Language Models (LLMs) into scientific research promises to accelerate discovery, yet a significant gap remains between the dynamic reasoning of Artificial Intelligence (AI) agents and the static, deterministic nature of canonical scientific workflows. This paper introduces ADEPT (Agentic Discovery and Exploration Platform for Tools), a reference architecture and pedagogical framework explicitly designed to bridge this gap. ADEPT's primary mission is to provide a transparent, "glass-box" environment where researchers and engineers can learn to effectively wrap established scientific software (e.g., BLAST, Nextflow pipelines) and compose it into reliable, agent-driven workflows. We describe its modular, multi-server architecture, which leverages the Model Context Protocol (MCP) for tool serving, LangGraph for robust agentic orchestration, and a secure nsjail-based sandbox for safe code execution. By prioritizing architectural clarity, safety, and modularity, ADEPT serves as an extensible blueprint for building trustworthy AI-augmented systems and fosters the collaborative development necessary to responsibly employ agentic AI for science. We provide practical examples of how to adapt and extend this framework, highlighting its utility in workforce development and AI-readiness capabilities across research and development projects.

97 MATHEMATICS AND COMPUTING↗

On the measurement of hardness at high strain rates by nanoindentation impact testing

Recent advances in electronics have enabled nanomechanical measurements with very low noise, fast time constants and high data acquisition rates. Furthermore, these capabilities open the door for a wide range of ultra-fast nanomechanical testing. Given the inherent dynamic nature of high-speed testing, a thorough understanding of the testing system's dynamics and electronics is extremely important for accurate measurements. In this work, an analytical framework that includes the mechanical and electronic contributions of the instrument and the material constitutive response is presented to provide guidelines for performing high strain rate measurements of hardness by nanoindentation testing. Simple closed-form solutions that provide insights on the choice of test methodology, test parameters and instrument design are presented along with the strain rate range over which accurate measurements can be performed with the commercially available nanoindenters.

36 MATERIALS SCIENCE↗

Bayesian inference of structured latent spaces from neural population activity with the orthogonal stochastic linear mixing model

The brain produces diverse functions, from perceiving sounds to producing arm reaches, through the collective activity of populations of many neurons. Determining if and how the features of these exogenous variables (e.g., sound frequency, reach angle) are reflected in population neural activity is important for understanding how the brain operates. Often, high-dimensional neural population activity is confined to low-dimensional latent spaces. However, many current methods fail to extract latent spaces that are clearly structured by exogenous variables. This has contributed to a debate about whether or not brains should be thought of as dynamical systems or representational systems. Here, we developed a new latent process Bayesian regression framework, the orthogonal stochastic linear mixing model (OSLMM) which introduces an orthogonality constraint amongst time-varying mixture coefficients, and provide Markov chain Monte Carlo inference procedures. We demonstrate superior performance of OSLMM on latent trajectory recovery in synthetic experiments and show superior computational efficiency and prediction performance on several real-world benchmark data sets. We primarily focus on demonstrating the utility of OSLMM in two neural data sets: μ ECoG recordings from rat auditory cortex during presentation of pure tones and multi-single unit recordings form monkey motor cortex during complex arm reaching. We show that OSLMM achieves superior or comparable predictive accuracy of neural data and decoding of external variables (e.g., reach velocity). Most importantly, in both experimental contexts, we demonstrate that OSLMM latent trajectories directly reflect features of the sounds and reaches, demonstrating that neural dynamics are structured by neural representations. Together, these results demonstrate that OSLMM will be useful for the analysis of diverse, large-scale biological time-series datasets.

59 BASIC BIOLOGICAL SCIENCES↗