Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA hybrid analytical-numerical algorithm for determining the neuronal current via EEG
362
0
0.0
(
0
)
Added by
Parham Hashemzadeh Dr
Publication date
2019
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
In this study, the neuronal current in the brain is represented using Helmholtz decomposition. It was shown in earlier work that data obtained via electroencephalography (EEG) are affected only by the irrotational component of the current. The irrotational component is denoted by $Psi$ and has support in the cerebrum. This inverse problem is severely ill-posed and requires that additional constraints are imposed. Here, we impose the requirement of the minimization of the $L_2$ norm of the current (energy). The function $Psi$ is expanded in terms of inverse multiquadric radial basis functions (RBF) on a uniform Cartesian grid inside the cerebrum. The minimal energy constraint in conjunction with the RBF parametrization of $Psi$ results in a Tikhonov regularized solution of $Psi$. The RBF shape parameter (regularization parameter), is computed by solving a 1-D non-linear maximization problem. Reconstructions are presented using synthetic data with a signal to noise ratio (SNR) of $20$ dB. The root mean square error (RMSE) between the exact and the reconstructed $Psi$ is RMSE=$0.1122$. The proposed reconstruction algorithm is computationally efficient and can be vectorized in MATLAB.
rate research
Read More
The pseudo-spectral analytical time-domain (PSATD) particle-in-cell (PIC) algorithm solves the vacuum Maxwells equations exactly, has no Courant time-step limit (as conventionally defined), and offers substantial flexibility in plasma and particle beam simulations. It is, however, not free of the usual numerical instabilities, including the numerical Cherenkov instability, when applied to relativistic beam simulations. This paper derives and solves the numerical dispersion relation for the PSATD algorithm and compares the results with corresponding behavior of the more conventional pseudo-spectral time-domain (PSTD) and finite difference time-domain (FDTD) algorithms. In general, PSATD offers superior stability properties over a reasonable range of time steps. More importantly, one version of the PSATD algorithm, when combined with digital filtering, is almost completely free of the numerical Cherenkov instability for time steps (scaled to the speed of light) comparable to or smaller than the axial cell size.
Tomography has made a radical impact on diverse fields ranging from the study of 3D atomic arrangements in matter to the study of human health in medicine. Despite its very diverse applications, the core of tomography remains the same, that is, a mathematical method must be implemented to reconstruct the 3D structure of an object from a number of 2D projections. In many scientific applications, however, the number of projections that can be measured is limited due to geometric constraints, tolerable radiation dose and/or acquisition speed. Thus it becomes an important problem to obtain the best-possible reconstruction from a limited number of projections. Here, we present the mathematical implementation of a tomographic algorithm, termed GENeralized Fourier Iterative REconstruction (GENFIRE). By iterating between real and reciprocal space, GENFIRE searches for a global solution that is concurrently consistent with the measured data and general physical constraints. The algorithm requires minimal human intervention and also incorporates angular refinement to reduce the tilt angle error. We demonstrate that GENFIRE can produce superior results relative to several other popular tomographic reconstruction techniques by numerical simulations, and by experimentally by reconstructing the 3D structure of a porous material and a frozen-hydrated marine cyanobacterium. Equipped with a graphical user interface, GENFIRE is freely available from our website and is expected to find broad applications across different disciplines.
A hybrid Maxwell solver for fully relativistic and electromagnetic (EM) particle-in-cell (PIC) codes is described. In this solver, the EM fields are solved in $k$ space by performing an FFT in one direction, while using finite difference operators in the other direction(s). This solver eliminates the numerical Cerenkov radiation for particles moving in the preferred direction. Moreover, the numerical Cerenkov instability (NCI) induced by the relativistically drifting plasma and beam can be eliminated using this hybrid solver by applying strategies that are similar to those recently developed for pure FFT solvers. A current correction is applied for the charge conserving current deposit to correctly account for the EM calculation in hybrid Yee-FFT solver. A theoretical analysis of the dispersion properties in vacuum and in a drifting plasma for the hybrid solver is presented, and compared with PIC simulations with good agreement obtained. This hybrid solver is applied to both 2D and 3D Cartesian and quasi-3D (in which the fields and current are decomposed into azimuthal harmonics) geometries. Illustrative results for laser wakefield accelerator simulation in a Lorentz boosted frame using the hybrid solver in the 2D Cartesian geometry are presented, and compared against results from 2D UPIC-EMMA simulation which uses a pure spectral Maxwell solver, and from OSIRIS 2D lab frame simulation using the standard Yee solver. Very good agreement is obtained which demonstrates the feasibility of using the hybrid solver for high fidelity simulation of relativistically drifting plasma with no evidence of the numerical Cerenkov instability.
This work presents a dynamic parallel distribution scheme for the Hartree-Fock exchange~(HFX) calculations based on the real-space NAO2GTO framework. The most time-consuming electron repulsion integrals~(ERIs) calculation is perfectly load-balanced with 2-level master-worker dynamic parallel scheme, the density matrix and the HFX matrix are both stored in the sparse format, the network communication time is minimized via only communicating the index of the batched ERIs and the final sparse matrix form of the HFX matrix. The performance of this dynamic scalable distributed algorithm has been demonstrated by several examples of large scale hybrid density-functional calculations on Tianhe-2 supercomputers, including both molecular and solid states systems with multiple dimensions, and illustrates good scalability.
This paper presents an extension of the hybrid scheme proposed by Wang et al. (J. Comput. Phys. 229 (2010) 169-180) for numerical simulation of compressible isotropic turbulence to flows with higher turbulent Mach numbers. The scheme still utilizes an 8th-order compact scheme with built-in hyperviscosity for smooth regions and a 7th-order WENO scheme for highly compressive regions, but now both in their conservation formulations and for the latter with the Roe type characteristic-wise reconstruction. To enhance the robustness of the WENO scheme without compromising its high-resolution and accuracy, the recursive-order-reduction procedure is adopted, where a new type of reconstruction-failure-detection criterion is constructed. To capture the upwind direction properly in extreme conditions, the global Lax-Friedrichs numerical flux is used. In addition, a new form of cooling function is proposed, which is proved to be positivity-preserving. With these techniques, the new scheme not only inherits the good properties of the original one but also extends largely the computable range of turbulent Mach number, which has been further confirmed by numerical results.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
