Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “visualization techniques”

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 631 records · Page 35

A Processing and Analytics System for Microscopy Data Workflows: The Pycroscopy Ecosystem of Packages

Major advancements in fields as diverse as biology and quantum computing have relied on a multitude of microscopy techniques. Despite the considerable proliferation of these instruments, significant bottlenecks remain in terms of processing, analysis, storage, and retrieval of the acquired datasets. Aside from lack of file standards, individual domain-specific analysis packages are often disjoint from the underlying datasets, and thus keeping track of analysis and processing steps remains tedious for the end-user, hampering reproducibility. Here, in this study, the pycroscopy ecosystem of packages is introduced, an open-source python-based ecosystem underpinned by a common data model. The data model, termed the N-dimensional spectral imaging data format, is realized in pycroscopy's sidpy package. This package is built on top of dask arrays, thus leveraging dask array attributes, but expanding them to accelerate microscopy relevant analysis and visualization. Several examples of the use of the pycroscopy ecosystem to create workflows for data ingestion and analysis of scanning transmission electron microscopy (STEM) and scanning probe microscopy data are shown. Adoption of such standardized routines will be critical to usher in the next generation of autonomous instruments where processing, computation, and meta-data storage will be critical to overall experimental operations.

97 MATHEMATICS AND COMPUTING↗

Framework Short-Range Order Observed in a Spinel-Type Li Superionic Conductor

Solid-state superionic conductors are characterized by rich structural disorders. Though structural complexities are central to their functionalities, they often give rise to short-range order that eludes detection by conventional diffraction-based techniques and is thus overlooked in establishing precise structure-property relationships. In this work, we synthesized single crystals of a recently discovered lithium (Li) superionic conductor Li16.2(1)In9.00(2)Sn1.10(1)O23.8 (LISO) for in-depth characterizations of structural subtleties. LISO exhibits an unusual spinel-like phase with significant Li overstoichiometry and a face-sharing Li network. Single-crystal neutron diffraction confirms significant Li disorder, as manifested in Li site splitting and partial occupancy. More importantly, synchrotron diffuse scattering combined with 3D-ΔPDF analysis and Monte Carlo simulations reveal short-range order in the nonalkali framework that might contribute to the phase stability and ionic conductivity. This work showcases an example in which subtle local energetics can be directly visualized in structurally disordered ionic conductors.

Chen, Yu↗

Evaluating Symmetry in 3-dimensional Simulation Data

The increasing use of three-dimensional simulations has led to an increased need for effective verification and validation, or V&V. Unfortunately, the high computational cost associated with three-dimensional simulations makes many standard V&V techniques cost-prohibitive. Given these circumstances, we must pay close attention to V&V metrics that can be evaluated using only a single simulation. One of the simplest of these metrics is the evaluation of physical symmetries. This report details a method of computing the deviation from symmetry in a three-dimensional dataset with an arbitrarily distributed mesh. Examples are presented for spherical and axisymmetric systems using the Paraview post-processing analysis workflow.

97 MATHEMATICS AND COMPUTING↗

Imaging solid–electrolyte interphase dynamics using operando reflection interference microscopy

The quality of the solid–electrolyte interphase is crucial for the performance of most battery chemistries, but its formation dynamics during operation are not well understood due to a lack of reliable operando characterization techniques. Herein, we report a dynamic, non-invasive, operando reflection interference microscope to enable the real-time imaging of the solid–electrolyte interphase during its formation and evolution processes with high sensitivity. The stratified structure of the solid–electrolyte interphase formed during four distinct steps includes the emergence of a permanent inner inorganic layer enriched in LiF, a transient assembly of an interfacial electrified double layer and a consequent emergence of a temporary outer organic-rich layer whose presence is reversible with electrochemical cycling. Reflection interference microscope imaging reveals an inverse correlation between the thicknesses of two interphasial subcomponents, implying that the permanent inorganic-rich inner layer dictates the organic-rich outer layer formation and lithium nucleation. In conclusion, the real-time visualization of solid–electrolyte interphase dynamics provides a powerful tool for the rational design of battery interphases.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Photoporomechanics: An Experimental Method to Visualize the Effective Stress Field in Fluid-Filled Granular Media

Effective stress governs the mechanical behavior of porous media. In this study, we use photoelasticimetry to visualize the evolving effective stress field in fluid-filled granular media in processes that couple fluid flow and mechanical deformation. We refer to this experimental method as photoporomechanics. We develop a fabrication process to produce millimeter-scale residual-stress-free photoelastic spheres with high geometric accuracy. We use color to quantify the forces acting on the particles over a wide range of forces, while using light intensity for a small range of forces. We then provide an application of photoporomechanics to illustrate the evolution of effective stress during one-dimensional consolidation: a process by which the stresses caused by a sudden load are gradually transmitted through a fluid-filled granular pack as the fluid drains and excess pore pressures dissipate. We report our technique provides a powerful experimental model system to study the grain-scale underpinning of coupled solid-fluid processes in granular media.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Nanoscale Antiferromagnetic Domain Imaging using Full-Field Resonant X-ray Magnetic Diffraction Microscopy

The physical properties of magnetic materials frequently depend not only on the microscopic spin and electronic structures, but also on the structures of mesoscopic length scales that emerge, for instance, from domain formations, or chemical and/or electronic phase separations. However, experimental access to such mesoscopic structures is currently limited, especially for antiferromagnets with net zero magnetization. Here, full-field microscopy and resonant magnetic X-ray diffraction are combined to visualize antiferromagnetic (AF) domains of the spin–orbit Mott insulator Sr 2 IrO 4 with area over ≈0.1 mm 2 and with spatial resolution as high as ≈150 nm. With the unprecedented wide field of views and high spatial resolution, an intertwining of two AF domains on a length comparable to the measured average AF domain wall width of 545 nm is revealed. This mesoscopic structure comprises a substantial portion of the sample surface, and thus can result in a macroscopic response unexpected from its microscopic magnetic structure. In particular, the symmetry analysis presented in this work shows that the inversion symmetry, which is preserved by the microscopic AF order, becomes ill-defined at the mesoscopic length scale. We report this result underscores the importance of this novel technique for a thorough understanding of the physical properties of antiferromagnets.

36 MATERIALS SCIENCE↗

Physics-Driven Convolutional Autoencoder Approach for CFD Data Compressions: Preprint

With the growing size and complexity of turbulent flow models, data compression approaches are of the utmost importance to analyze, visualize, or restart the simulations. Recently, in-situ autoencoder-based compression approaches have been proposed and shown to be effective at producing reduced representations of turbulent flow data. However, these approaches focus solely on training the model using point-wise sample reconstruction losses that do not take advantage of the physical properties of turbulent flows. In this paper, we show that training autoencoders with additional physics-informed regularizations, e.g., enforcing incompressibility and preserving enstrophy, improves the compression model in three ways: (i) the compressed data better conform to known physics for homogeneous isotropic turbulence without negatively impacting point-wise reconstruction quality, (ii) inspection of the gradients of the trained model uncovers changes to the learned compression mapping that can facilitate the use of explainability techniques, and (iii) as a performance byproduct, training losses are shown to converge up to 12x faster than the baseline model.

auto-encoders↗

Visualization of the surface distribution of photocatalytic activity at the anatase/rutile TiO2 interface using X-ray photoelectron emission microscopy

We applied X-ray photoelectron emission microscopy (XPEEM) to visualize the surface distribution of photocatalytic activity at the anatase/rutile (A/R) interface of titanium dioxide. Acetic acid was adsorbed on the surface, and its decomposition and desorption reactions under ultraviolet light irradiation were monitored by observing C 1s spectra, enabling direct observation of the spatial distribution of photocatalytic activity. Complementary crystal structure information was obtained using spatially resolved low-energy electron diffraction, low-energy electron microscopy, X-ray absorption spectroscopy, and micro-focused Raman spectroscopy. Comparison with these structural data allowed for the evaluation of the photocatalytic activity as a function of the distance from the A/R interface. The activity is enhanced in the vicinity of the interface and gradually decreases toward both the anatase and rutile regions. These results demonstrate that XPEEM is an effective technique for spatially resolved mapping of photocatalytic activity, which has\\r\\npreviously been inaccessible using conventional characterization techniques.

25 ENERGY STORAGE↗

Speeding-up fuzzing through directional seeds

Abstract Fuzzing is an automated process for discovering inputs in a program that may trigger unexpected behavior. Today, fuzzing has become a standard practice for the discovery of bugs and security vulnerabilities. However, the main issue with such practices is that the exploration of the input space of programs can often be prohibitively expensive. Therefore, several alternative fuzzing strategies have been introduced during the last few years. Some fuzzing techniques rely on human expertise to provide a plausible set of initial input examples, namely, seeds. However, the process of handcrafting seeds for fuzzing purposes often becomes strenuous for humans as it requires a deeper understanding of the Program-Under-Test (PUT). Also, the use of known inputs to programs often does not trigger vulnerable program behavior or may not reach potentially vulnerable code locations. To address those issues, we propose a seed generation framework that enables Human-In-The-Loop (HITL) directed fuzzing where the human assumes a more active role in the creation of seeds that can penetrate and assess desired locations of the PUT. Our proposed framework uses Symbolic Execution (SE) to generate seeds that exercise paths to target program locations. Moreover, our framework enables the visualization of the explored execution paths in the binary of the PUT for the generated seeds. We evaluated our approach on a set of 12 carefully designed C programs with diverse characteristics that mimic real-world programs. The experimental results show the effectiveness of the proposed approach in improving the performance of standard fuzzing tools such as the American Fuzzy Lop ("Image missing" <#comment/> ). Specifically, our solution can generate seeds that substantially enhance the performance of the fuzzer, achieving speedups ranging from $$1.46\times $$ 1.46 × to $$68.53\times $$ 68.53 × for branch conditions, $$1.39\times $$ 1.39 × to $$254.62\times $$ 254.62 × for branch depths, $$14,879.59\times $$ 14 , 879.59 × to $$30,295.88\times $$ 30 , 295.88 × for branch widths over traditional seeds. Additionally, the speedup increases with the number of target function ranging from $$12,260\times $$ 12 , 260 × to $$22,856.07\times $$ 22 , 856.07 × over traditional seeds while only requiring less than 15 seconds on average for the seed generation step.

97 MATHEMATICS AND COMPUTING↗

Quantum dynamics of cold atomic gas with SU(1,1) symmetry

Motivated by recent advances in quantum dynamics, we investigate the dynamics of the system with SU(1,1) symmetry. Instead of performing the time-ordered integral for the evolution operator of the time-dependent Hamiltonian, we show that the time evolution operator can be expressed as an SU(1,1) group element. Since the SU(1,1) group describes the “rotation” on a hyperbolic surface, the dynamics can be visualized on a Poincaré disk, a stereographic projection of the upper hyperboloid. As an example, we present the trajectory of the revival of Bose-Einstein condensation and that of the scale-invariant Fermi gas on the Poincaré disk. Further considering quantum gases in an oscillating lattice, we also study the dynamics of the system with time-dependent single-particle dispersion.

74 ATOMIC AND MOLECULAR PHYSICS↗

Lithium Metal Battery Quality Control via Transformer–CNN Segmentation

Lithium metal battery (LMB) has the potential to be the next-generation battery system because of its high theoretical energy density. However, defects known as dendrites are formed by heterogeneous lithium (Li) plating, which hinders the development and utilization of LMBs. Non-destructive techniques to observe the dendrite morphology often use X-ray computed tomography (XCT) to provide cross-sectional views. To retrieve three-dimensional structures inside a battery, image segmentation becomes essential to quantitatively analyze XCT images. This work proposes a new semantic segmentation approach using a transformer-based neural network called TransforCNN that is capable of segmenting out dendrites from XCT data. In addition, we compare the performance of the proposed TransforCNN with three other algorithms, U-Net, Y-Net, and E-Net, consisting of an ensemble network model for XCT analysis. Our results show the advantages of using TransforCNN when evaluating over-segmentation metrics, such as mean intersection over union (mIoU) and mean Dice similarity coefficient (mDSC), as well as through several qualitatively comparative visualizations.

Quenum, Jerome (ORCID:0000000271265853)↗

Nanoscale mapping of point defect concentrations with $\mathrm{4D-STEM}$

Vacancies are missing atoms in a crystalline material, and occur both at equilibrium (varying with temperature) and out of equilibrium such as when crystalline materials are damaged with radiation or corrosion. While we know of their importance, particularly regarding diffusive mechanisms, it is not straightforward to experimentally measure their concentration or to visualize them directly. Traditionally, measurements such as positron annihilation spectroscopy and to some degree X-ray diffraction can measure average concentrations, but generally lack the ability to visualize or quantify a heterogenous concentration of vacancies that can occur at the level of individual defects and microstructural features. Here, we present a method to map vacancy concentrations and their distribution using local lattice parameter measurements with a high-resolution electron microscope. Our method utilizes a Au thin film as a model to demonstrate the method via four-dimensional scanning transmission electron microscopy (4D-STEM) by correlating the differences between changes in lattice parameter and the volumetric thermal expansion during in situ heating experiments. The vacancy mapping methodology is also applied to non-equilibrium defects accumulated in pure Al via knock-on electron beam irradiation. Our method demonstrates the ability to map point defect concentrations in heterogeneous systems in situ with nanometer spatial resolution. The result is a technique that can provide direct measurements of vacancy concentrations at the level of individual defects in studies of materials in and out of equilibrium.

36 MATERIALS SCIENCE↗

A remarkable CHF of 345W/cm 2 is achieved in a wicked-microchannel using HFE-7100

Boiling heat transfer of dielectric fluids is a promising cooling technique for thermal management of microelectronic systems. However, the critical heat flux achievable is generally low, because of the poor thermophysical properties of these fluids. Here, to address this dilemma, we propose a new cooling concept to substantially enhance global liquid supply during phase-change process using enhanced capillary-driven force. Additionally, dedicated vapor pathways are designed among a bank of micro-pillars to facilitate vapor removal. In this work, new wicks comprised of silicon micro-pinfin arrays are explored to significantly enhance the flow boiling heat transfer performance. To examine the functionalities of this wick, experiments on HFE-7100 were carried out with mass velocities varying from 247 to 3,465 kg/m 2 s. To explore the enhancement mechanisms and to analyze the capillary-assisted flow boiling process, visualization studies were conducted. The results indicate that sustainable evaporation induced by wick microstructures and efficient liquid supply are the enhancement mechanisms compared to parallel microchannels with solid walls. It is found that the overall heat transfer coefficient is substantially increased up to 75%. Remarkably, a high critical heat flux (CHF) of approximately 345 W/cm 2 is recorded at G = 3,465 kg/m 2 s at coolant inlet temperature of ~20 °C. Equally importantly, this noticeable enhancement of CHF value is associated with drastically decreased pressure drops compared to microchannels decorated with μ-pinfin fences.

36 MATERIALS SCIENCE↗

Phase-contrast 3D tomography of HeLa cells grown in PLLA polymer electrospun scaffolds using synchrotron X-rays

Advanced imaging is useful for understanding the three-dimensional (3D) growth of cells. X-ray tomography serves as a powerful noninvasive, nondestructive technique that can fulfill these purposes by providing information about cell growth within 3D platforms. There are a limited number of studies taking advantage of synchrotron X-rays, which provides a large field of view and suitable resolution to image cells within specific biomaterials. In this study, X-ray synchrotron radiation microtomography at Diamond Light Source and advanced image processing were used to investigate cellular infiltration of HeLa cells within poly L-lactide (PLLA) scaffolds. This study demonstrates that synchrotron X-rays using phase contrast is a useful method to understand the 3D growth of cells in PLLA electrospun scaffolds. Two different fiber diameter (2 and 4 µm) scaffolds with different pore sizes, grown over 2, 5 and 8 days in vitro , were examined for infiltration and cell connectivity. After performing visualization by segmentation of the cells from the fibers, the results clearly show deeper cell growth and higher cellular interconnectivity in the 4 µm fiber diameter scaffold. Finally, this indicates the potential for using such 3D technology to study cell–scaffold interactions for future medical use.

3D culture↗

Nondestructive Evaluation of Stress Corrosion Cracking in a Welded Steel Plate Using Guided Ultrasonic Waves

Stress corrosion cracking (SCC) had occurred in early-generation high-level nuclear waste tanks constructed by welding carbon steel. This paper describes an ultrasonic inspection system and its fundamental ability to detect and quantify the length of SCC on thick welded steel plates. The finite element method (FEM) was applied to simulate the welding process to estimate the welding residual stress field. Growth of stress corrosion cracks is driven by crack stress intensities exceeding the subcritical cracking threshold intensity. The subject plate was experimentally inspected with ultrasonic nondestructive evaluation (NDE) techniques to characterize the extent of SCC. The NDE system uses a piezoelectric transducer to generate guided waves in the thick steel plate, and a scanning laser Doppler vibrometer (SLDV) to measure multidimensional time–space wavefield data over a user-defined scanning area in the plate surface. The measured wavefield data can show wave interactions in a localized area in the plate due to the presence of the discontinuities of the SCC. To generate an inspection image that can precisely show the crack’s location and/or the dimension, the wavefield data are further processed to generate inspection image that maps the entire sample plate so the crack can be clearly identified in the plate while its length can be readily estimated. In conclusion, the ultrasonic test results for crack length agree well with the visually estimated length and are close to that predicted by the FEM for cracks in the weld residual stress field.

36 MATERIALS SCIENCE↗

Sub-pilot-scale Production of High-Value Products from U.S. Coals

Investigators from the University of Utah, University of Wyoming and Marshall University pursued a program to study the conversion of raw coal to high-value products of carbon fiber and silicon carbide. Team members also developed an initial framework for a data portal that can incorporate laboratory data on coal processing and product quality, and also work with tools for machine learning for data analysis, data visualization and economic assessment. Experimental R&D efforts focused on the conversion of raw coal to coal tar and other byproducts, and the resulting tar intermediates were upgraded to form anisotropic and isotropic pitch materials. These pitch materials were produced from coal using both thermal (pyrolysis) and chemical (mild solvolysis liquefaction) decomposition of raw coal. Four different coals were studied: Utah bituminous coal (Sufco), Wyoming PRB coal (Black Thunder), Illinois bituminous coal (Illinois #6), and West Virginia bituminous coal (Flying Eagle). Both metallurgical-grade coking coals and lower-grade steam coals were investigated, and controlled secondary gas-phase reactions were used during a two-stage pyrolysis process to induce cracking and condensation reactions among the pyrolytic tar species. This approach successfully improved the performance of the lower grade coals for yielding pitch materials, with properties more consistent with a commercial-grade pitch that had previously demonstrated success for quality carbon fiber production. The use of waste plastic materials was also studied, to help improve physical and chemical characteristics of the intermediate tars and final pitch product; in particular, for lowering the pitch softening point to an acceptable level for melt spinning carbon fiber. Mild solvolysis liquefaction was also used as a method for producing pitch for carbon fiber production. As expected, significantly higher pitch yields were obtained using this approach, and waste plastic materials were also successfully used to reduce pitch softening point to an acceptable level. The plastic materials were also utilized to create a solvent for the mild solvolysis process, and this plastic-derived solvent was shown to provide results consistent with more expensive commercial chemical solvents, and could thus avoid the need for costly recovery and recycle of a liquefaction solvent. Additional experimental R&D focused on the production of silicon carbide (β-SiC) from the residual char byproduct from pitch production, and also on the production of carbon fiber from the anisotropic pitch. SiC was successfully synthesized using a mixture of residual char and sandstone at a ratio of 1:1. Reaction temperature and residence time were optimized and yielded a product purity of 81%. For carbon fiber production, the most successful pitch samples were obtained from the mild solvolysis liquefaction approach, combined with the use of a plastic (HDPE)-derived solvent. Fiber properties improved over time as laboratory fiber production methodologies improved, and final yields of carbon fiber were obtained with a diameter of 12.14 ± 1.10 um, Modulus of 173.73 ± 15.25 GPa, and Tensile Strength of 1.04 ± 0.10 GPa. A proof-of-concept Modern Community Research Data Portal (MCRDP) was developed and deployed for coal and coal-derived pitch characterization, with the full support of (i) remote web-based access, (ii) distributed analysis, (iii) interactive visualization and exploration, (iv) shared and long-term data access, (v) advanced query capabilities and (vi) real-time collaboration. The Coal to Products Data Portal “coaltoproducts.org” provides researchers with space to store and share data within a project, tools for analyzing and understanding data for scientific investigation, and the ability to publish data to the broader community for reproducibility. The portal leverages the Material Commons 2.0 (MC) platform developed by the Center for PRedictive Integrated Structural Materials Science (PRISMS) of the University of Michigan, to achieve long-term longevity of data collections and, more importantly, collaborative science. A number of data visualization tools were also assessed and implemented for interrogating the experimental and modeling data. The machine learning portion of this project analyzed datasets from two different coal conversion processes performed on a diverse set of coal samples from both the coal pyrolysis experiments and the solvent liquefaction experiments. The work was initiated by exploring standard regression models on the pyrolysis data, aiming to understand the impact of sample characteristics and processing conditions on key product metrics. Over the course of the project, the focus expanded to include a variety of machine learning tools, delving into both supervised and unsupervised learning methods. Models tested on the pyrolysis data included linear, ridge, lasso, elastic-net, Gaussian process, random forest regression, and AutoSklearn, and the approach was continually refined to enhance predictive accuracy and model interpretability. Similar techniques were applied to the liquefaction data with an additional focus on feature engineering. Along with mesophase content, additional outputs of interest were the pitch yield, softening point, and QI content. Insights derived from these analyses are crucial in determining the factors influencing the quality and yield of coal-derived products. As the work progressed, the research evolved from foundational model comparisons to analyses of random forests, decision paths, and feature importance scores. A thorough market analysis was performed to examine the prospects of coal-based carbon fibers. The best opportunities for coal come from its lower and more stable price relative to petroleum, particularly for subbituminous coals, which is the primary advantage that a coal refinery may have over a petroleum refinery. Before a commercial CTP production facility can be modeled, however, several things need to be understood regarding the nature of the would-be coal refinery. These include the technology to be deployed, the size of facility, the volume(s) of co-product(s), and the waste and emissions profile of the plant. The volume of co-products and waste may be substantial and will require separate market analysis to ensure viability. In the near-term, the importance of coal tar pitch, in the form of carbon pitch, to the aluminum and steel industries is likely to overshadow the alternative use of this material as an input for carbon fiber. The importance of steel and aluminum in building materials, and the need for carbon materials in their manufacturing, will ensure that demand for these products remains for the long run. In addition, carbon fiber may also be the best substitute for steel and aluminum well into the future. While society will eventually be able to shift production of much of its electricity needs to renewables, it will not be able to shift away from fossil fuels for production of high-strength construction and vehicular materials. Demand for carbon fiber is expected to increase quickly, but the volume of carbon fiber and the amount of coal that would be needed to produce even a sizeable share of this market may still be relatively small compared to current coal production. Thus, other coal-based products like graphene, graphite, carbon foams, resins, and carbon-based building products will play important roles in sustaining coal production as coal-fired power generation continues to decline.

01 COAL, LIGNITE, AND PEAT↗

A Method for Producing Hierarchical and Statistically Calibrated Predictions of Nuclear Material Properties from Existing Models

Computer vision-based analysis of micrographs of nuclear materials is an emerging technique for property prediction, synthetic route identification, and other material analysis tasks. These analysis tasks play a pivotal role in many material characterization applications such as signature development for treaty verification, process optimization, etc. The backbone in many of the recent computer vision-based techniques is a deep learning model, which takes a fixed-size set of pixels and provides a class prediction for that set of pixels. For example, previous work developed a deep convolutional neural network (CNN) to predict the synthetic route from a 256 px x 256 px patch taken from a larger image of uranium ore concentrates. In this work, we present several methods for first calibrating these models in a manner that they can provide accurate probabilities of their predictions’ veracity, and several methods of combining these probabilities. Overall, the combination of these two steps into a pipeline allows for full-image and even full-sample (where a sample has many images) predictions with associated confidence values. Finally, we show that one can also use the patch predictions and confidence to produce a visualization to map predicted constituents through the image. Results and examples for predicting and mapping uranium ore concentrates’ synthetic process from imagery will be presented.

artificial intelligence↗

Multimodal parameter spaces of a complex multi-channel neuron model

One of the most common types of models that helps us to understand neuron behavior is based on the Hodgkin–Huxley ion channel formulation (HH model). A major challenge with inferring parameters in HH models is non-uniqueness: many different sets of ion channel parameter values produce similar outputs for the same input stimulus. Such phenomena result in an objective function that exhibits multiple modes (i.e., multiple local minima). This non-uniqueness of local optimality poses challenges for parameter estimation with many algorithmic optimization techniques. HH models additionally have severe non-linearities resulting in further challenges for inferring parameters in an algorithmic fashion. To address these challenges with a tractable method in high-dimensional parameter spaces, we propose using a particular Markov chain Monte Carlo (MCMC) algorithm, which has the advantage of inferring parameters in a Bayesian framework. The Bayesian approach is designed to be suitable for multimodal solutions to inverse problems. We introduce and demonstrate the method using a three-channel HH model. We then focus on the inference of nine parameters in an eight-channel HH model, which we analyze in detail. We explore how the MCMC algorithm can uncover complex relationships between inferred parameters using five injected current levels. The MCMC method provides as a result a nine-dimensional posterior distribution, which we analyze visually with solution maps or landscapes of the possible parameter sets. The visualized solution maps show new complex structures of the multimodal posteriors, and they allow for selection of locally and globally optimal value sets, and they visually expose parameter sensitivities and regions of higher model robustness. We envision these solution maps as enabling experimentalists to improve the design of future experiments, increase scientific productivity and improve on model structure and ideation when the MCMC algorithm is applied to experimental data.

97 MATHEMATICS AND COMPUTING↗