BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Jacob K. White

Jacob K. White appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
Computational Engineering, Finance, and Scienceeess.SPmath.NAmath.OCphysics.comp-phcond-mat.mtrl-scimath-phmath.MPMolecular Networksphysics.med-phphysics.opticsquant-phQuantitative MethodsTissues and Organs

Catalog footprint

What is connected

12works
14topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

12 published item(s)

preprint2022arXiv

A hybrid volume-surface integral equation method for rapid electromagnetic simulations in MRI

Objective: We developed a hybrid volume surface integral equation (VSIE) method based on domain decomposition to perform fast and accurate magnetic resonance imaging (MRI) simulations that include both remote and local conductive elements. Methods: We separated the conductive surfaces present in MRI setups into two domains and optimized electromagnetic (EM) modeling for each case. Specifically, interactions between the body and EM waves originating from local radiofrequency (RF) coils were modeled with the precorrected fast Fourier transform, whereas the interactions with remote conductive surfaces (RF shield, scanner bore) were modeled with a novel cross tensor train-based algorithm. We compared the hybrid- VSIE with other VSIE methods for realistic MRI simulation setups. Results: The hybrid-VSIE was the only practical method for simulation using 1 mm voxel isotropic resolution (VIR). For 2 mm VIR, our method could be solved at least 23 times faster and required 760 times lower memory than traditional VSIE methods. Conclusion: The hybrid-VSIE demonstrated a marked improvement in terms of convergence times of the numerical EM simulation compared to traditional approaches in multiple realistic MRI scenarios. Significance: The efficiency of the novel hybrid-VSIE method could enable rapid simulations of complex and comprehensive MRI setups.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Compression of volume-surface integral equation matrices via Tucker decomposition for magnetic resonance applications

In this work, we propose a method for the compression of the coupling matrix in volume\hyp surface integral equation (VSIE) formulations. VSIE methods are used for electromagnetic analysis in magnetic resonance imaging (MRI) applications, for which the coupling matrix models the interactions between the coil and the body. We showed that these effects can be represented as independent interactions between remote elements in 3D tensor formats, and subsequently decomposed with the Tucker model. Our method can work in tandem with the adaptive cross approximation technique to provide fast solutions of VSIE problems. We demonstrated that our compression approaches can enable the use of VSIE matrices of prohibitive memory requirements, by allowing the effective use of modern graphical processing units (GPUs) to accelerate the arising matrix\hyp vector products. This is critical to enable numerical MRI simulations at clinical voxel resolutions in a feasible computation time. In this paper, we demonstrate that the VSIE matrix\hyp vector products needed to calculate the electromagnetic field produced by an MRI coil inside a numerical body model with $1$ mm$^3$ voxel resolution, could be performed in $\sim 33$ seconds in a GPU, after compressing the associated coupling matrix from $\sim 80$ TB to $\sim 43$ MB.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Magnetic-resonance-based electrical property mapping using Global Maxwell Tomography with an 8-channel head coil at 7 Tesla: a simulation study

Objective: Global Maxwell Tomography (GMT) is a recently introduced volumetric technique for noninvasive estimation of electrical properties (EP) from magnetic resonance measurements. Previous work evaluated GMT using ideal radiofrequency (RF) excitations. The aim of this simulation study was to assess GMT performance with a realistic RF coil. Methods: We designed a transmit-receive RF coil with $8$ decoupled channels for $7$T head imaging. We calculated the RF transmit field ($B_1^+$) inside heterogeneous head models for different RF shimming approaches, and used them as input for GMT to reconstruct EP for all voxels. Results: Coil tuning/decoupling remained relatively stable when the coil was loaded with different head models. Mean error in EP estimation changed from $7.5\%$ to $9.5\%$ and from $4.8\%$ to $7.2\%$ for relative permittivity and conductivity, respectively, when changing head model without re-tuning the coil. Results slightly improved when an SVD-based RF shimming algorithm was applied, in place of excitation with one coil at a time. Despite errors in EP, RF transmit field ($B_1^+$) and absorbed power could be predicted with less than $0.5\%$ error over the entire head. GMT could accurately detect a numerically inserted tumor. Conclusion: This work demonstrates that GMT can reliably reconstruct EP in realistic simulated scenarios using a tailored 8-channel RF coil design at $7$T. Future work will focus on construction of the coil and optimization of GMT's robustness to noise, to enable in-vivo GMT experiments. Significance: GMT could provide accurate estimations of tissue EP, which could be used as biomarkers and could enable patient-specific estimation of RF power deposition, which is an unsolved problem for ultra-high-field magnetic resonance imaging.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

VoxCap: FFT-Accelerated and Tucker-Enhanced Capacitance Extraction Simulator for Voxelized Structures

VoxCap, a fast Fourier transform (FFT)-accelerated and Tucker-enhanced integral equation simulator for capacitance extraction of voxelized structures, is proposed. The VoxCap solves the surface integral equations (SIEs) for conductor and dielectric surfaces with three key attributes that make the VoxCap highly CPU and memory efficient for the capacitance extraction of the voxelized structures: (i) VoxCap exploits the FFTs for accelerating the matrix-vector multiplications during the iterative solution of linear system of equations arising due to the discretization of SIEs. (ii) During the iterative solution, VoxCap uses a highly effective and memory-efficient preconditioner that reduces the number of iterations significantly. (iii) VoxCap employs Tucker decompositions to compress the block Toeplitz and circulant tensors, requiring the largest memory in the simulator. By doing so, it reduces the memory requirement of these tensors from hundreds of gigabytes to a few megabytes and the CPU time required to obtain Toeplitz tensors from tens of minutes (even hours) to a few seconds for very large scale problems. VoxCap is capable of accurately computing capacitance of arbitrarily shaped and large-scale voxelized structures on a desktop computer.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

GMRES-Accelerated ADMM for Quadratic Objectives

We consider the sequence acceleration problem for the alternating direction method-of-multipliers (ADMM) applied to a class of equality-constrained problems with strongly convex quadratic objectives, which frequently arise as the Newton subproblem of interior-point methods. Within this context, the ADMM update equations are linear, the iterates are confined within a Krylov subspace, and the General Minimum RESidual (GMRES) algorithm is optimal in its ability to accelerate convergence. The basic ADMM method solves a $κ$-conditioned problem in $O(\sqrtκ)$ iterations. We give theoretical justification and numerical evidence that the GMRES-accelerated variant consistently solves the same problem in $O(κ^{1/4})$ iterations for an order-of-magnitude reduction in iterations, despite a worst-case bound of $O(\sqrtκ)$ iterations. The method is shown to be competitive against standard preconditioned Krylov subspace methods for saddle-point problems. The method is embedded within SeDuMi, a popular open-source solver for conic optimization written in MATLAB, and used to solve many large-scale semidefinite programs with error that decreases like $O(1/k^{2})$, instead of $O(1/k)$, where $k$ is the iteration index.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Parameter Insensitivity in ADMM-Preconditioned Solution of Saddle-Point Problems

We consider the solution of linear saddle-point problems, using the alternating direction method-of-multipliers (ADMM) as a preconditioner for the generalized minimum residual method (GMRES). We show, using theoretical bounds and empirical results, that ADMM is made remarkably insensitive to the parameter choice with Krylov subspace acceleration. We prove that ADMM-GMRES can consistently converge, irrespective of the exact parameter choice, to an $ε$-accurate solution of a $κ$-conditioned problem in $O(κ^{2/3}\logε^{-1})$ iterations. The accelerated method is applied to randomly generated problems, as well as the Newton direction computation for the interior-point solution of semidefinite programs in the SDPLIB test suite. The empirical results confirm this parameter insensitivity, and suggest a slightly improved iteration bound of $O(\sqrtκ\logε^{-1})$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Fluctuating volume-current formulation of electromagnetic fluctuations in inhomogeneous media: incandecence and luminescence in arbitrary geometries

We describe a fluctuating volume--current formulation of electromagnetic fluctuations that extends our recent work on heat exchange and Casimir interactions between arbitrarily shaped homogeneous bodies [Phys. Rev. B. 88, 054305] to situations involving incandescence and luminescence problems, including thermal radiation, heat transfer, Casimir forces, spontaneous emission, fluorescence, and Raman scattering, in inhomogeneous media. Unlike previous scattering formulations based on field and/or surface unknowns, our work exploits powerful techniques from the volume--integral equation (VIE) method, in which electromagnetic scattering is described in terms of volumetric, current unknowns throughout the bodies. The resulting trace formulas (boxed equations) involve products of well-studied VIE matrices and describe power and momentum transfer between objects with spatially varying material properties and fluctuation characteristics. We demonstrate that thanks to the low-rank properties of the associatedmatrices, these formulas are susceptible to fast-trace computations based on iterative methods, making practical calculations tractable. We apply our techniques to study thermal radiation, heat transfer, and fluorescence in complicated geometries, checking our method against established techniques best suited for homogeneous bodies as well as applying it to obtain predictions of radiation from complex bodies with spatially varying permittivities and/or temperature profiles.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Generalized Taylor-Duffy Method for Efficient Evaluation of Galerkin Integrals in Boundary-Element Method Computations

We present a generic technique, automated by computer-algebra systems and available as open-source software \cite{scuff-em}, for efficient numerical evaluation of a large family of singular and nonsingular 4-dimensional integrals over triangle-product domains, such as those arising in the boundary-element method (BEM) of computational electromagnetism. To date, practical implementation of BEM solvers has often required the aggregation of multiple disparate integral-evaluation schemes to treat all of the distinct types of integrals needed for a given BEM formulation; in contrast, our technique allows many different types of integrals to be handled by the \emph{same} algorithm and the same code implementation. Our method is a significant generalization of the Taylor--Duffy approach \cite{Taylor2003,Duffy1982}, which was originally presented for just a single type of integrand; in addition to generalizing this technique to a broad class of integrands, we also achieve a significant improvement in its efficiency by showing how the \emph{dimension} of the final numerical integral may often be reduced by one. In particular, if $n$ is the number of common vertices between the two triangles, in many cases we can reduce the dimension of the integral from $4-n$ to $3-n$, obtaining a closed-form analytical result for $n=3$ (the common-triangle case).

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the use of rational-function fitting methods for the solution of 2D Laplace boundary-value problems

A computational scheme for solving 2D Laplace boundary-value problems using rational functions as the basis functions is described. The scheme belongs to the class of desingularized methods, for which the location of singularities and testing points is a major issue that is addressed by the proposed scheme, in the context of the 2D Laplace equation. Well-established rational-function fitting techniques are used to set the poles, while residues are determined by enforcing the boundary conditions in the least-squares sense at the nodes of rational Gauss-Chebyshev quadrature rules. Numerical results show that errors approaching the machine epsilon can be obtained for sharp and almost sharp corners, nearly-touching boundaries, and almost-singular boundary data. We show various examples of these cases in which the method yields compact solutions, requiring fewer basis functions than the Nyström method, for the same accuracy. A scheme for solving fairly large-scale problems is also presented.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Estimating Reaction Rate Parameters in Cell Signaling Pathways Using Extreme Decomposition and Belief Propagation Tailored for Data-Rich Cases

Motivation: Modeling biological signaling networks using ordinary differential equations (ODEs) has proven to be a powerful technique for generating insight into cellular dynamics, but it typically requires estimating rate parameters based on experimentally observed concentrations. New measurement methods can measure concentrations for all molecular species in a pathway, which creates a new opportunity to decompose the optimization of rate parameters. Results: In contrast with conventional methods that minimize the disagreement between simulated and observed concentrations, the BPPE method fits a spline curve through the observed concentration points, and then matches the derivatives of the spline-curve to the production and consumption of each species. Whereas traditional methods follow the ODEs exactly and then attempt to match the data, BPPE follows the data exactly and then attempts to match the ODEs. The new objective function is an extreme decomposition of the problem because each factor in the function is enforcing the equality of one ODE at one timeslice. A "loopy belief propagation" algorithm solves this factorized approximation of the parameter estimation problem providing systematic coverage of the search space and unique asymptotic behavior; the run time is polynomial in the number of molecules and timepoints, but exponential in the degree of the biochemical network. The implementation is a global-local hybrid optimization, and we compare with the performance of local, global, and hybrid methods. BPPE is demonstrated for a novel model of Akt activation dynamics including redox-mediated inactivation of PTEN. Availability: Software and supplementary information are available at http://webbppe.nus.edu.sg:8080/opal2/WebBPPE . Contact: LisaTK@nus.edu.sg . Keywords: probabilistic graphical models, physico-chemical modeling, systems biology, signal transduction.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A novel boundary element method using surface conductive absorbers for full-wave analysis of 3-D nanophotonics

Fast surface integral equation (SIE) solvers seem to be ideal approaches for simulating 3-D nanophotonic devices, as these devices generate fields both in an interior channel and in the infinite exterior domain. However, many devices of interest, such as optical couplers, have channels that can not be terminated without generating reflections. Generating absorbers for these channels is a new problem for SIE methods, as the methods were initially developed for problems with finite surfaces. In this paper we show that the obvious approach for eliminating reflections, making the channel mildly conductive outside the domain of interest, is inaccurate. We describe a new method, in which the absorber has a gradually increasing surface conductivity; such an absorber can be easily incorporated in fast integral equation solvers. Numerical experiments from a surface-conductivity modified FFT-accelerated PMCHW-based solver are correlated with analytic results, demonstrating that this new method is orders of magnitude more effective than a volume absorber, and that the smoothness of the surface conductivity function determines the performance of the absorber. In particular, we show that the magnitude of the transition reflection is proportional to 1/L^(2d+2), where L is the absorber length and d is the order of the differentiability of the surface conductivity function.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Calculation of nonzero-temperature Casimir forces in the time domain

We show how to compute Casimir forces at nonzero temperatures with time-domain electromagnetic simulations, for example using a finite-difference time-domain (FDTD) method. Compared to our previous zero-temperature time-domain method, only a small modification is required, but we explain that some care is required to properly capture the zero-frequency contribution. We validate the method against analytical and numerical frequency-domain calculations, and show a surprising high-temperature disappearance of a non-monotonic behavior previously demonstrated in a piston-like geometry.

0 citations0 reviewsNo code linkedOpen access