
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 KarhumenLo\'{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 nonbroken rays,
brokenonce rays, and others, respectively, 2) use the the nonbroken ray data
and an effective datadriven layer stripping strategyan optimization based
iterative imaging methodto recover the medium velocity outside the convex
hull of the scatterers, and 3) use selected brokenonce ray data to recover the
boundary of the scatterersa direct imaging method. By numerical tests, we
show that our hybrid method can recover both the unknown medium and the
nottooconcave scatterers efficiently and robustly.

A numerical method using implicit surface representations is proposed to
solve the linearized PoissonBoltzmann 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 scatteringfree 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 highresolution 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 [RenZhao, SIAM J. Imag. Sci., 2013] on the same problem
but in the diffusive regime.