An official website of the United States government
Neurons exhibit complex geometry in their branched networks of neurites which is essential to the function of individual neuron but also brings challenges to transport a wide variety of essential materials throughout their neurite networks for their survival and function. While numerical methods like isogeometric analysis (IGA) have been used for modeling the material transport process via solving partial differential equations (PDEs), they require long computation time and huge computation resources to ensure accurate geometry representation and solution, thus limit their biomedical application. Here we present a graph neural network (GNN)-based deep learning model to learn the IGA-based material transport simulation and provide fast material concentration prediction within neurite networks of any topology. Given input boundary conditions and geometry configurations, the well-trained model can predict the dynamical concentration change during the transport process with an average error less than 10% and 120∼330 times faster compared to IGA simulations. The effectiveness of the proposed model is demonstrated within several complex neurite networks.
more »
« less
Reaction diffusion system prediction based on convolutional neural network
Abstract The reaction-diffusion system is naturally used in chemistry to represent substances reacting and diffusing over the spatial domain. Its solution illustrates the underlying process of a chemical reaction and displays diverse spatial patterns of the substances. Numerical methods like finite element method (FEM) are widely used to derive the approximate solution for the reaction-diffusion system. However, these methods require long computation time and huge computation resources when the system becomes complex. In this paper, we study the physics of a two-dimensional one-component reaction-diffusion system by using machine learning. An encoder-decoder based convolutional neural network (CNN) is designed and trained to directly predict the concentration distribution, bypassing the expensive FEM calculation process. Different simulation parameters, boundary conditions, geometry configurations and time are considered as the input features of the proposed learning model. In particular, the trained CNN model manages to learn the time-dependent behaviour of the reaction-diffusion system through the input time feature. Thus, the model is capable of providing concentration prediction at certain time directly with high test accuracy (mean relative error <3.04%) and 300 times faster than the traditional FEM. Our CNN-based learning model provides a rapid and accurate tool for predicting the concentration distribution of the reaction-diffusion system.
more »
« less
- Award ID(s):
- 1804929
- PAR ID:
- 10154064
- Publisher / Repository:
- Nature Publishing Group
- Date Published:
- Journal Name:
- Scientific Reports
- Volume:
- 10
- Issue:
- 1
- ISSN:
- 2045-2322
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)Abstract Neurons exhibit complex geometry in their branched networks of neurites which is essential to the function of individual neuron but also brings challenges to transport a wide variety of essential materials throughout their neurite networks for their survival and function. While numerical methods like isogeometric analysis (IGA) have been used for modeling the material transport process via solving partial differential equations (PDEs), they require long computation time and huge computation resources to ensure accurate geometry representation and solution, thus limit their biomedical application. Here we present a graph neural network (GNN)-based deep learning model to learn the IGA-based material transport simulation and provide fast material concentration prediction within neurite networks of any topology. Given input boundary conditions and geometry configurations, the well-trained model can predict the dynamical concentration change during the transport process with an average error less than 10% and $$120 \sim 330$$ 120 ∼ 330 times faster compared to IGA simulations. The effectiveness of the proposed model is demonstrated within several complex neurite networks.more » « less
-
Purpose To develop a method of biologically guided deep learning for post-radiation 18 FDG-PET image outcome prediction based on pre-radiation images and radiotherapy dose information. Methods Based on the classic reaction–diffusion mechanism, a novel biological model was proposed using a partial differential equation that incorporates spatial radiation dose distribution as a patient-specific treatment information variable. A 7-layer encoder–decoder-based convolutional neural network (CNN) was designed and trained to learn the proposed biological model. As such, the model could generate post-radiation 18 FDG-PET image outcome predictions with breakdown biological components for enhanced explainability. The proposed method was developed using 64 oropharyngeal patients with paired 18 FDG-PET studies before and after 20-Gy delivery (2 Gy/day fraction) by intensity-modulated radiotherapy (IMRT). In a two-branch deep learning execution, the proposed CNN learns specific terms in the biological model from paired 18 FDG-PET images and spatial dose distribution in one branch, and the biological model generates post-20-Gy 18 FDG-PET image prediction in the other branch. As in 2D execution, 718/233/230 axial slices from 38/13/13 patients were used for training/validation/independent test. The prediction image results in test cases were compared with the ground-truth results quantitatively. Results The proposed method successfully generated post-20-Gy 18 FDG-PET image outcome prediction with breakdown illustrations of biological model components. Standardized uptake value (SUV) mean values in 18 FDG high-uptake regions of predicted images (2.45 ± 0.25) were similar to ground-truth results (2.51 ± 0.33). In 2D-based Gamma analysis, the median/mean Gamma Index (<1) passing rate of test images was 96.5%/92.8% using the 5%/5 mm criterion; such result was improved to 99.9%/99.6% when 10%/10 mm was adopted. Conclusion The developed biologically guided deep learning method achieved post-20-Gy 18 FDG-PET image outcome predictions in good agreement with ground-truth results. With the breakdown biological modeling components, the outcome image predictions could be used in adaptive radiotherapy decision-making to optimize personalized plans for the best outcome in the future.more » « less
-
Inference, optimization, and inverse problems are but three examples of mathematical operations that require the repeated solution of a complex system of mathematical equations. To this end, surrogates are often used to approximate the output of these large computer simulations, providing fast and cheap approximation solutions. Statistical emulators are surrogates that, in addition to predicting the mean behavior of the system, provide an estimate of the error in that prediction. Classical Gaussian stochastic process emulators predict scalar outputs based on a modest number of input parameters. Making predictions across a space-time field of input variables is not feasible using classical Gaussian process methods. Parallel partial emulation is a new statistical emulator methodology that predicts a field of outputs based on the input parameters. Parallel partial emulation is constructed as a Gaussian process in parameter space, but no correlation among space or time points is assumed. Thus the computational work of parallel partial emulation scales as the cube of the number of input parameters (as traditional Gaussian Process emulation) and linearly with a space-time grid. The numerical methods used in numerical simulations are often designed to exploit properties of the equations tobe solved. For example, modern solvers for hyperbolic conservation laws satisfy conservation at each time step, insuring overall conservation of the physical variables. Similarly, symplectic methods are used to solve Hamiltonian problems in physics. It is of interest, then, to study whether parallel partial emulation predictions inherit properties possessed by the simulation outputs. Does an emulated solution of a conservation law preserve the conserved quantities? Does an emulator of a Hamiltonian system preserve the energy? This paper investigates the properties of emulator predictions, in the context of systems of partial differential equations. We study conservation properties for three different kinds ofequations-conservation laws, reaction-diffusion systems, and a Hamiltonian system.We also investigate the effective convergence, in parameter space, of the predicted solution of a highly nonlinear system modeling shape memory alloys.more » « less
-
Abstract While critical to biological processes, molecular diffusion is difficult to quantify, and spatial mapping of local diffusivity is even more challenging. Here we report a machine-learning-enabled approach, pixels-to-diffusivity (Pix2D), to directly extract the diffusion coefficientDfrom single-molecule images, and consequently enable super-resolvedDspatial mapping. Working with single-molecule images recorded at a fixed framerate under typical single-molecule localization microscopy (SMLM) conditions, Pix2D exploits the often undesired yet evident motion blur,i.e., the convolution of single-molecule motion trajectory during the frame recording time with the diffraction-limited point spread function (PSF) of the microscope. Whereas the stochastic nature of diffusion imprints diverse diffusion trajectories to different molecules diffusing at the same givenD, we construct a convolutional neural network (CNN) model that takes a stack of single-molecule images as the input and evaluates aD-value as the output. We thus validate robustDevaluation and spatial mapping with simulated data, and with experimental data successfully characterizeDdifferences for supported lipid bilayers of different compositions and resolve gel and fluidic phases at the nanoscale.more » « less
