• Random fields are commonly used for modeling of spatially (or timely) dependent stochastic processes. In this study, we provide a characterization of the intrinsic complexity of a random field in terms of its second order statistics, e.g., the covariance function, based on the Karhumen-Lo\'{e}ve expansion. We then show scaling laws for the intrinsic complexity of a random field in terms of the correlation length as it goes to 0. In the discrete setting, it becomes approximate embeddings of a set of random vectors. We provide a precise scaling law when the random vectors have independent and identically distributed entires using random matrix theory as well as when the random vectors has a specific covariance structure.
  • We investigate a hybrid inverse problem in fluorescence ultrasound modulated optical tomography (fUMOT) in the diffusive regime. We prove that the absorption coefficient of the fluorophores at the excitation frequency and the quantum efficiency coefficient can be uniquely and stably reconstructed from boundary measurement of the photon currents, provided that some background medium parameters are known. Reconstruction algorithms are proposed and numerically implemented as well.
  • We present a hybrid imaging method for a challenging travel time tomography problem which includes both unknown medium and unknown scatterers in a bounded domain. The goal is to recover both the medium and the boundary of the scatterers from the scattering relation data on the domain boundary. Our method is composed of three steps: 1) preprocess the data to classify them into three different categories of measurements corresponding to non-broken rays, broken-once rays, and others, respectively, 2) use the the non-broken ray data and an effective data-driven layer stripping strategy--an optimization based iterative imaging method--to recover the medium velocity outside the convex hull of the scatterers, and 3) use selected broken-once ray data to recover the boundary of the scatterers--a direct imaging method. By numerical tests, we show that our hybrid method can recover both the unknown medium and the not-too-concave scatterers efficiently and robustly.
  • A numerical method using implicit surface representations is proposed to solve the linearized Poisson-Boltzmann equations that arise in mathematical models for the electrostatics of molecules in solvent. The proposed method used an implicit boundary integral formulation to derived a linear system defined on Cartesian nodes in a narrowband surrounding the closed surface that separate the molecule and the solvent. The needed implicit surfaces is constructed from the given atomic description of the molecules, by a sequence of standard level set algorithms. A fast multipole method is applied to accelerate the solution of the linear system. A few numerical studies involving some standard test cases are presented and compared to other existing results.
  • We propose in this work a fast numerical algorithm for solving the equation of radiative transfer (ERT) in isotropic media. The algorithm has two steps. In the first step, we derive an integral equation for the angularly averaged ERT solution by taking advantage of the isotropy of the scattering kernel, and solve the integral equation with a fast multipole method (FMM). In the second step, we solve a scattering-free transport equation to recover the original ERT solution. Numerical simulations are presented to demonstrate the performance of the algorithm for both homogeneous and inhomogeneous media.
  • Fluorescence photoacoustic tomography (fPAT) is a molecular imaging modality that combines photoacoustic tomography (PAT) with fluorescence imaging to obtain high-resolution imaging of fluorescence distributions inside heterogeneous media. The objective of this work is to study inverse problems in the quantitative step of fPAT where we intend to reconstruct physical coefficients in a coupled system of radiative transport equations using internal data recovered from ultrasound measurements. We derive uniqueness and stability results on the inverse problems and develop some efficient algorithms for image reconstructions. Numerical simulations based on synthetic data are presented to validate the theoretical analysis. The results we present here complement these in [Ren-Zhao, SIAM J. Imag. Sci., 2013] on the same problem but in the diffusive regime.