Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Software behavior computation”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Scientific Core Library Stack (SCLS) v2026

SCLS (Scientific Core Library Stack) is an opinionated build and packaging system for scientific computing libraries developed at Lawrence Berkeley National Laboratory. It produces a coherent, reproducible stack of numerical libraries — including BLAS/LAPACK, MPI, sparse direct and iterative solvers, graph partitioners, and parallel I/O libraries (e.g., PETSc, SLEPc, HDF5, NetCDF, MUMPS, OpenBLAS) — that work together without manual repair by downstream scientific software. From a single recipe-and-flavor model, SCLS produces native RPM packages for RHEL-family Linux, DEB packages for Debian/Ubuntu, direct Unix-style prefix installs for HPC and locked-down environments, and native macOS builds. Multiple build "flavors" (e.g., GCC+OpenBLAS, GCC+MKL, Intel+MKL, debug) coexist in distinct prefixes on the same host. Compared to general-purpose meta-build frameworks, SCLS is deliberately curated rather than infinitely configurable. It enforces deterministic, audit-friendly behavior: explicit build dependencies, no silent feature autodetection, a clear open-source license policy, and rpath-based runtime linkage so installs integrate cleanly with standard package-manager workflows.

Messe, Christian [Lawrence Berkeley National Labor↗
DOE OSTI2025

Towards exact finite temperature electronic structure in solids and molecules (Final Technical Report)

This report describes the University of Iowa portion of a project that is now continuing at Michigan State University. We are developing novel methods, algorithms, and software to enable simulations of molecules and materials at high temperature. This is a key challenge in chemistry and materials science. By refining an approach called Density Matrix Quantum Monte Carlo (DMQMC), we developed faster and more accurate ways to conduct these simulations. These advances will help us understand how temperature affects the behavior of electrons, chemical bonds, and phase transitions in solids and molecules. These breakthroughs are especially important for applications where light and heat drive chemical reactions, superconductivity, and materials used in energy and sensing. In addition, this project involved the development of the open-source HANDE-QMC software package, supporting the broader community in benchmarking and developing finite-temperature electronic structure methods.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Reaction Mechanism Generator v3.0: Advances in Automatic Mechanism Generation

In chemical kinetics research, kinetic models containing hundreds of species and tens of thousands of elementary reactions are commonly used to understand and predict the behavior of reactive chemical systems. Reaction Mechanism Generator (RMG) is a software suite developed to automatically generate such models by incorporating and extrapolating from a database of known thermochemical and kinetic parameters. Here, we present the recent version 3 release of RMG and highlight improvements since the previously published description of RMG v1.0. Most notably, RMG can now generate heterogeneous catalysis models in addition to the previously available gas- and liquid-phase capabilities. For model analysis, new methods for local and global uncertainty analysis have been implemented to supplement first-order sensitivity analysis. The RMG database of thermochemical and kinetic parameters has been significantly expanded to cover more types of chemistry. The present release includes parallelization for faster model generation and a new molecule isomorphism approach to improve computational performance. RMG has also been updated to use Python 3, ensuring compatibility with the latest cheminformatics and machine learning packages. Overall, RMG v3.0 includes many changes which improve the accuracy of the generated chemical mechanisms and allow for exploration of a wider range of chemical systems.

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

Yearly report for the period Jan. 2018 – Feb. 2019 IC Project: w17_faultprediction “Critical Stress in Earth”

The consequence of earthquake is usually catastrophic. Currently, we cannot determine when a fault is approaching failure due to either anthropogenic activities or tectonic movement. In an anthropogenically perturbed setting, the rate of physical changes on faults far exceeds most natural processes with currently unpredictable and unmeasurable results. In a tectonic setting we generally do not know where a fault is in the earthquake cycle beyond noting how long has passed since the last rupture, which is itself sometimes unknown. The Laboratory has therefore funded a multi million multi-year LDRD-DR program that comprises experiments, numerical modeling, and machine learning in order to develop and test a highly novel approach to probe the earthquake cycle and in particular, to detect and locate stressed faults that are approaching failure. For this program, modeling of a series of sheared granular layers undergoing stick-slip behavior using the EES-17’s Combined Finite-Discrete Element Method (FDEM) based code HOSS (Hybrid Optimization Software Suite) plays a pivotal role since numerical simulations of granular fault gouge allow for analysis of the mechanical behavior of the system at a level of spatial and temporal resolution not accessible experimentally, and also offer the possibility of identifying the most important parameters governing earthquake dynamics. The simulated data also provide very detailed stress, strain and motion information for machine learning purpose to predict earthquake. The code – HOSS – has been developed in Geophysics Group in LANL during the past ten years and is used to simulate the stick-slip behavior of granular fault gouge. One of the key feature of HOSS is that it is superior to pure finite element method and discrete element method and can provide very detailed evolvement information of granular fault gouge. Additionally, its parallel computing capability makes large scale earthquake gouge simulation possible. Taking these facts into account, and also because few other fully parallel FDEM code exists, this IC project allows LANL to stay at the forefront of earthquake gouge and granular material simulation in the world.

42 ENGINEERING↗
DOE OSTI2011

BitPredator: A Discovery Algorithm for BitTorrent Initial Seeders and Peers

There is a large amount of illegal content being replicated through peer-to-peer (P2P) networks where BitTorrent is dominant; therefore, a framework to profile and police it is needed. The goal of this work is to explore the behavior of initial seeds and highly active peers to develop techniques to correctly identify them. We intend to establish a new methodology and software framework for profiling BitTorrent peers. This involves three steps: crawling torrent indexers for keywords in recently added torrents using Really Simple Syndication protocol (RSS), querying torrent trackers for peer list data and verifying Internet Protocol (IP) addresses from peer lists. We verify IPs using active monitoring methods. Peer behavior is evaluated and modeled using bitfield message responses. We also design a tool to profile worldwide file distribution by mapping IP-to-geolocation and linking to WHOIS server information in Google Earth.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2016

Recent Advances in Analytical Electron Microscopy for Irradiated Materials

Analytical electron microscopy (AEM) in the scanning transmission electron microscope (STEM) has long been a staple of post-irradiation examination (PIE) and materials development campaigns. Recently, new advances in both STEM hardware and post-processing data-analysis software, along with ancillary advances such as in focused ion beam (FIB) sample preparation, has led to exciting advances in the quantity and quality of data that can be obtained. This talk will summarize recent experiments performed at Oak Ridge National Laboratory (ORNL), especially work on neutron-irradiated materials and fuels taking advantage of a newly-installed latest-generation AEM tool in the ORNL irradiated materials laboratories. ORNL installed a new FEI Talos F200X microscope, a 200 keV, high-brightness Schottky emission AEM system optimized for high-speed, high-efficiency X-ray mapping experiments. The F200X features the FEI 'SuperX' four-sector silicon drift detector (SDD) X-ray detector, which provides 0.9 sr solid angle of collection. By combining the high-brightness electron source, latest-generation electron optics, and fast computers and pulse processors, X-ray mapping speed is increased by a factor of 50-100 x compared to the previous-generation instruments. In addition to improvements in hardware, data mining via multivariate statistical analysis (MVSA) allows for maximum insight to the data obtained from the new hardware. MVSA data mining of the X-ray maps allows, for instance, unbiased computational identification of the phases present, and discovery of small phases or precipitates not expected by the human analyst. This talk will present examples of experiments performed using the most modern AEM tools and software. Silicon carbide is a major candidate for next-generation nuclear applications. Nano-infiltration by transient eutectic (NITE) SiC requires low sintering temperatures. X-ray mapping of the NITE structure find pockets of different oxides, such as Zr-Y-oxide, Y-Al-oxide, and yttria. MVSA data mining finds the regions where different oxides are present and indicates that Y additions preferentially wet grain boundaries, where Zr additions are preferentially concentrated at grain triple points. SiC coatings from AGR-I experiment TRISO particles were found to contain significant quantities of fission products (e.g., Pd) as well as actinides (U, Pu). MVSA data mining found complex precipitates and grain boundary films of the fission products and actinides deep (∼several microns) into the SiC layer. Advanced AEM methods, combined with MVSA data mining, makes possible both rapid and detailed PIE. Ferritic alloys, such as accident-tolerant FeCrAl grades, benefit greatly from advanced AEM examination. Neutron-irradiated FeCrAl precipitates the detrimental Cr-rich α' phase. Because the SuperX SDD system has four detectors, it is possible to obtain high X-ray count rates even on a highly tilted specimen, as might be needed to obtain STEM images of dislocation loops. Results on neutron-irradiated FeCrAl indicates that α' size and number density can be obtained by AEM, and MVSA results determine that minor elements (such as Si) has complex radiation-induced segregation behavior at grain boundaries. We have also developed and validated computer models of the AEM experiment, allowing different microstructures' X-ray map results to be simulated for different STEM instruments. This allows X-ray mapping to be more quantitative, such as by estimating what size nanoclusters in an oxide-dispersion strengthened steel would be visible or invisible in STEM. In summary, new advanced AEM technology provides a major leap forward in experiments possible in PIE of irradiated materials. (authors)

36 MATERIALS SCIENCE↗
DOE OSTI2023

Thermal Modeling of Advanced Test Reactor Fuel in a Generalized Dry Storage System Under Hypothetical Accident Conditions

Work was performed to investigate the thermal behavior of a sealed aluminum-clad spent nuclear fuel (ASNF) dry storage configuration under hypothetical accident conditions during long-term storage. The considered system consists of a concrete dry storage overpack; a welded over canister, backfilled with argon; and nine Department of Energy (DOE) standardized canisters (DOESCs) loaded with ASNF, backfilled with helium gas. One technical concern associated with the long-term storage of ASNF includes radiolytic hydrogen generation. The yield of large quantities of hydrogen could lead to DOESC over pressurization or a flammable internal atmosphere. While flammability concerns can be resolved by preventing the ingress of oxygen, the canister pressure is partially controlled by the atmosphere temperature. This memorandum summarizes the results of modeling and simulation (M&S) work completed with Star CCM+ (a computational fluid dynamics [CFD] software) considering the thermal effects of a fully engulfing flame and another scenario with the vent ports of the concrete dry storage overpack completely blocked. The goal was to determine the bounding maximum temperature of the DOESC internal atmosphere to evaluate a worst-case pressure scenario, as well as the critical vent port blockage durations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Investigation of Best-Practices and Computationally Inexpensive Radiative Exchange Models for Discrete Element Method Modeling of Aluminosilicate Particles in Concentrating Solar Power Environments

Chemically inert, aluminosilicate based particles have been investigated as both a thermal transport and sensible energy storage medium for concentrating solar power facilities. These particles will experience a wide range of operating temperatures (300-1000 K) and handling conditions (dense to dilute falling particle curtains, dense granular flows, or dense structures), requiring specially-designed and optimized infrastructures. The relative influence of collisional and frictional interactions between particles varies based on temperature-dependent particulate properties and greatly impacts the bulk, granular flow behavior. These underlying physics are captured using discrete element method modeling tools. However, this modeling method is computationally expensive as each particle position and interaction is tracked during the simulation. These modeling methods are further complicated by introducing temperature-dependent particle properties, high-temperature radiative exchange, and directional irradiation sources experienced by granular flows in concentrating solar power environments. In this study, coupled experimental and numerical slump testing of aluminosilicate particles was performed and computationally efficient radiative exchange models were evaluated to establish best-practices for discrete element method models for concentrating solar power environments. The three particle types investigated included Carbobead HSP 30 /60, Carbobead CP 30/60, and Granusil 4030. Existing modeling limitations and computationally-efficient multi-modal heat transfer models were evaluated using Aspherix®, a commercial discrete element method software. High-temperature (< 1073 K) slump testing of aluminosilicate particles was performed to investigate the deviation between experimentally-observed and numerically-predicted angles of repose introduced by computation-time reduction practices including the relaxation of the particle elastic modulus and coarse-graining. Coarse-graining is used to use a single modeled particle that is representative of a collection of smaller particles, decreasing the computational cost at the expense of geometric accuracy. Additionally, relaxation of the elastic modulus is used to reduce computational time at the expense of an increased, modeled particle overlap. Prior studies have determined that aluminosilicate particles retain a high elastic modulus at high temperatures (< 1073 K), requiring small simulation timesteps to ensure resolved contact forces resemble appropriate solid mechanics. A parametric study was performed to evaluate the influence of computation time improvements on the deviation between experimental and modeled angle of repose across high temperatures < 1073 K. Additionally, numerical case studies were performed on candidate particle systems at varying porosities and temperatures. These studies were performed to investigate the influence of computationally-efficient radiative-exchange modeling methods coupled to Aspherix® on modeled accuracy and computation time. The recently-developed distance-based approximation was evaluated in estimating radiative exchange between particles and participating surfaces located in close proximity. The distance based approximation was developed to use tabulated estimates of the radiative distribution factor between individual particles and surfaces in close proximity (< 40 particle radii). These methods were expanded to the aluminosilicate particles of interest, including the influence of particle size distributions. To capture radiative exchange between particles and surfaces not in close proximity (> 40 particle radii) and to capture the absorption of directional irradiation from concentrating solar resources, a volumetrically-averaged radiative distribution factor was calculated between the modeled granular flow and surfaces using Monte Carlo ray-tracing for participating media. Volume-averaged absorption and scattering coefficients were predicted using a volumetric discretization of the modeled domain with monodisperse approximations based on geometric optics and experimentally-determined scattering phase functions for aluminosilicate particles.

14 SOLAR ENERGY↗
DOE OSTI2020

Battery Electric Vehicle Eco-Cooperative Adaptive Cruise Control in the Vicinity of Signalized Intersections

This study develops a connected eco-driving controller for battery electric vehicles (BEVs), the BEV Eco-Cooperative Adaptive Cruise Control at Intersections (Eco-CACC-I). The developed controller can assist BEVs while traversing signalized intersections with minimal energy consumption. The calculation of the optimal vehicle trajectory is formulated as an optimization problem under the constraints of (1) vehicle acceleration/deceleration behavior, defined by a vehicle dynamics model; (2) vehicle energy consumption behavior, defined by a BEV energy consumption model; and (3) the relationship between vehicle speed, location, and signal timing, defined by vehicle characteristics and signal phase and timing (SPaT) data shared under a connected vehicle environment. The optimal speed trajectory is computed in real-time by the proposed BEV eco-CACC-I controller, so that a BEV can follow the optimal speed while negotiating a signalized intersection. The proposed BEV controller was tested in a case study to investigate its performance under various speed limits, roadway grades, and signal timings. In addition, a comparison of the optimal speed trajectories for BEVs and internal combustion engine vehicles (ICEVs) was conducted to investigate the impact of vehicle engine types on eco-driving solutions. Lastly, the proposed controller was implemented in microscopic traffic simulation software to test its networkwide performance. The test results from an arterial corridor with three signalized intersections demonstrate that the proposed controller can effectively reduce stop-and-go traffic in the vicinity of signalized intersections and that the BEV Eco-CACC-I controller produces average savings of 9.3% in energy consumption and 3.9% in vehicle delays.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2014

Overview of New Tools to Perform Safety Analysis: BWR Station Black Out Test Case

Dynamic Probabilistic Risk Assessment (DPRA) methodologies couple system simulator codes (e.g., RELAP, MELCOR) with simulation controller codes (e.g., RAVEN, ADAPT). While system simulator codes accurately model system dynamics deterministically, simulation controller codes introduce both deterministic (e.g., system control logic, operating procedures) and stochastic (e.g., component failures, parameter uncertainties) elements into the simulation. Typically, a DPRA is performed by: 1) sampling values of a set of parameters from the uncertainty space of interest (using the simulation controller codes), and 2) simulating the system behavior for that specific set of parameter values (using the system simulator codes). For complex systems, one of the major challenges in using DPRA methodologies is to analyze the large amount of information (i.e., large number of scenarios ) generated, where clustering techniques are typically employed to allow users to better organize and interpret the data. In this paper, we focus on the analysis of a nuclear simulation dataset that is part of the Risk Informed Safety Margin Characterization (RISMC) Boiling Water Reactor (BWR) station blackout (SBO) case study. We apply a software tool that provides the domain experts with an interactive analysis and visualization environment for understanding the structures of such high-dimensional nuclear simulation datasets. Our tool encodes traditional and topology-based clustering techniques, where the latter partitions the data points into clusters based on their uniform gradient flow behavior. We demonstrate through our case study that both types of clustering techniques complement each other in bringing enhanced structural understanding of the data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

DXT: Darshan eXtended Tracing

As modern supercomputers evolve to exascale, their I/O subsystems are becoming increasingly complex, making optimization of I/O for scientific applications a daunting task. Although I/O profiling tools facilitate the process of optimizing application I/O performance, legacy profiling tools lack flexibility in their level of detail and ability to correlate traces with other sources of data. Additionally, a lack of robust trace analysis tools makes it difficult to derive actionable insights from large-scale I/O traces. Darshan is an HPC I/O characterization tool that records statistics in a lightweight manner that makes it appropriate for full-time production deployment. However, Darshan’s default characterization mechanism records information at a fixed granularity. We augment Darshan by proposing Darshan eXtended Tracing (DXT) for more detailed profiling of I/O software stacks. DXT enables users and administrators to vary the level of fidelity captured by Darshan at run time without modifying or recompiling applications. This capability facilitates systematic analysis on the I/O behavior of applications and can provide useful application kernel I/O traces to help advance parallel I/O research. We have demonstrated the power of DXT by obtaining a wide range of useful statistics for multiple case studies, and we further show that DXT is able to do so same with negligible overhead.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Studying performance portability of LAMMPS across diverse GPU-based platforms

The molecular dynamics simulation software, LAMMPS, utilizes the Kokkos acceleration library to port computation to a diverse set of architectures including those based on GPU accelerators. In addition to Kokkos, LAMMPS contains a vast code base that leverages the CUDA application programming interface using library functions such as cuFFT, CUDA's fast-fourier transform (FFT) library, and, more recently, also support for AMD's Heterogeneous Interface for Portability (HIP) that is rapidly growing. While preparing LAMMPS tests for the AMD GPU-based test system precursors to Frontier, we investigated several strategies for accelerating LAMMPS on AMD GPUs, using the AMD Instinct MI100 and MI250X. In this work, we integrated the HIP FFT library, hipFFT, into the particle-particle particle-mesh (PPPM) long-range solver, which allowed the porting of PPPM calculations to the GPUs. Kokkos behavior on the MI100 and MI250X was also investigated through the package kokkos command of LAMMPS, targeting communication, memory usage, and particle grid decomposition. The Tersoff, Reax, Lennard-Jones (LJ), EAM, Granular, and PPPM potentials were investigated in this effort, and results from these experiments are provided. In conclusion, the selected potentials were run on Spock (AMD Instinct MI100), Crusher (AMD Instinct MI250X), AFW HPC11 (NVIDIA A100) and Summit (NVIDIA V100), for comparison. Operational roofline models were constructed and analyzed for the Tersoff, Reax, and Lennard–Jones potentials on Crusher and Summit.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Studying Performance Portability of LAMMPS across Diverse GPU-based Platforms

The molecular dynamics simulation software, LAMMPS, utilizes the Kokkos acceleration library to port computation to a diverse set of architectures including those based on GPU accelerators. In addition to Kokkos, LAMMPS contains a vast code base that leverages the CUDA application programming interface using library functions such as cuFFT, CUDA’s fast-fourier transform (FFT) library, and, more recently, also support for AMD’s Heterogeneous Interface for Portability (HIP) that is rapidly growing. While preparing LAMMPS tests for the AMD GPU-based test system precursors to Frontier, we investigated several strategies for accelerating LAMMPS on AMD GPUs, using the AMD Instinct MI100 and MI250X. In this work, we integrated the HIP FFT library, hipFFT, into the particle-particle particle-mesh (PPPM) long-range solver, which allowed the porting of PPPM calculations to the GPUs. Kokkos behavior on the MI100 and MI250X was also investigated through the package kokkos command of LAMMPS, targeting com- munication, memory usage, and particle grid decomposition. The Tersoff, Reax, Lennard-Jones (LJ), EAM, Granular, and PPPM potentials were investigated in this effort, and results from these experiments are provided. The selected potentials were run on Spock (AMD Instinct MI100), Crusher (AMD Instinct MI250X), AFW HPC11 (NVIDIA A100) and Summit (NVIDIA V100), for comparison. Operational roofline models were constructed and analyzed for the Tersoff, Reax, and Lennard-Jones potentials on Crusher and Summit.

Hagerty, Nick↗
DOE OSTI2012

Building Agent™ Software (BA) [SWRs 12-05 and 17-21]

Energy interfaces provide a new level of motivational information. Many studies have reported energy savings after dashboards for operators are installed, but these provide little actionable information to the occupant and the occupant cannot report individual comfort. The Building Agent™ (BA), in development at the National Renewable Energy Laboratory (NREL), provides an interface for the occupant to converse with the building control system and aims to actively engage and affect occupant behavior to achieve positive energy outcomes. Energy savings achieved via occupant behavioral changes could be as high as 25%, and could rival savings achieved by technological measures. For example, energy audits often show that more than 50% of the energy used in a building is consumed at night. Occupants leaving equipment and lighting on after they leave work is a large contributing factor. Similar savings opportunities have been identified when occupants were away from their desks 50% of the time, while most of their computers, lights, and other equipment remained powered on. Building Agent aims to create cohesion between building automation system (BAS) data, local measurements, and occupant feedback to provide visualizations, empowering occupants and building engineers to take diagnostic, proactive, energy-saving, and comfort-improving actions. The multilayer BA architecture consists of central server, measurement, and distributed application layers. The server layer aggregates global and localized measurements with occupant feedback to correlate perceived comfort with actual data throughout the building. The measurement layer complements BAS data with localized, environmental condition measurements at workstations. The application layer provides a bidirectional communications mechanism between the building and the occupants that accepts comfort feedback and provides energy performance information. The interface layer also uses three displays to influence building energy use: empathetic and personalized gauges, graphical data, and competitive rankings. Building Agent was implemented at NREL in the Research Support Facility (RSF), a LEED Platinum and Net Zero office building, which has received many awards. The constant feedback from occupants has enabled BA to better enhance the BAS data controls. Continued work is being done on BA to increase visualization of data and enhanced the user experience.

Long, Nicholas↗
DOE OSTI2018

HyRAM V1.1 User Guide

Hydrogen Risk Assessment Models (HyRAM) is a software toolkit that provides a basis for quantitative risk assessment and consequence modeling for hydrogen infrastructure and transportation systems. HyRAM integrates validated, analytical models of hydrogen behavior, statistics, and a standardized QRA approach to generate useful, repeatable data for the safety analysis of various hydrogen systems. HyRAM is a software developed by Sandia National Laboratories for the U.S. Department of Energy. This document demonstrates how to use HyRAM to recreate a hydrogen system and obtain relevant data regarding potential risk. Specific examples are utilized throughout this document, providing detailed tutorials of HyRAM features with respect to hydrogen system safety analysis and risk assessment.

08 HYDROGEN↗
DOE OSTI2019

Quantifying Uncertainty in Emulations (LDRD Report)

This report summarizes the work performed under the project “Quantifying Uncertainty in Emulations.” Emulation can be used to model real-world systems, typically using virtualization to run the real software on virtualized hardware. Emulations are increasingly used to answer mission-oriented questions, but how well they represent the real-world systems is still an open area of research. The goal of the project was to quantify where and how emulations differ from the real world. To do so, we ran a representative workload on both, and collected and compared metrics to identify differences. We aimed to capture behaviors, rather than performance, differences as the latter is more well-understood in the literature.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Algorithm and program for precise determination of unit-cell parameters of single crystal taking into account the sample eccentricity

A technique has been developed to refine the unit-cell parameters of single crystals with minimization of the influence of instrumental errors on the result. The corresponding computational procedure HuberUB is added to the software package of Huber-5042 diffractometer with a point detector and closedcycle helium cryostat Displex DE-202. The parameters of unit cell, its orientation, the goniometer zero angles, the sample eccentricity, the distances in the goniometer, and the radiation wavelength were refined by the nonlinear least-squares method, which allows imposition of constraints on the unit-cell parameters, depending on the crystal symmetry. The technique is approved on a LuB{sub 12} single crystal. The unit-cell parameters are determined in a temperature range of 20–295 K, with an absolute error not larger than 0.0004 Å (the relative error is of ~5 × 10{sup –5}). The estimates of the unit-cell parameters obtained by the proposed method are evidenced to be unbiased. Some specific features of the behavior of parameters in the ranges of 120–140 and 20–50 K are revealed, which correlate with the anomalies of the physical properties of the crystal.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Causal network inference using biochemical kinetics

Motivation: Networks are widely used as structural summaries of biochemical systems. Statistical estimation of networks is usually based on linear or discrete models. However, the dynamics of biochemical systems are generally non-linear, suggesting that suitable non-linear formulations may offer gains with respect to causal network inference and aid in associated prediction problems. Results: We present a general framework for network inference and dynamical prediction using time course data that is rooted in nonlinear biochemical kinetics. This is achieved by considering a dynamical system based on a chemical reaction graph with associated kinetic parameters. Both the graph and kinetic parameters are treated as unknown; inference is carried out within a Bayesian framework. This allows prediction of dynamical behavior even when the underlying reaction graph itself is unknown or uncertain. Results, based on (i) data simulated from a mechanistic model of mitogen-activated protein kinase signaling and (ii) phosphoproteomic data from cancer cell lines, demonstrate that non-linear formulations can yield gains in causal network inference and permit dynamical prediction and uncertainty quantification in the challenging setting where the reaction graph is unknown. Availability and implementation: MATLAB R2014a software is available to download from warwick.ac.uk/chrisoates.

60 APPLIED LIFE SCIENCES↗