Publications

We consider the inverse acoustic obstacle problem for sound-soft star-shaped obstacles in two dimensions wherein the boundary of the obstacle is determined from measurements of the scattered field at a collection of receivers outside the object. One of the standard approaches for solving this problem is to reformulate it as an optimization problem: finding the boundary of the domain that minimizes the $L^2$ distance between computed values of the scattered field and the given measurement data. The optimization problem is computationally challenging since the local set of convexity shrinks with increasing frequency and results in an increasing number of local minima in the vicinity of the true solution. In many practical experimental settings, low frequency measurements are unavailable due to limitations of the experimental setup or the sensors used for measurement. Thus, obtaining a good initial guess for the optimization problem plays a vital role in this environment.
We present a neural network warm-start approach for solving the inverse scattering problem, where an initial guess for the optimization problem is obtained using a trained neural network. We demonstrate the effectiveness of our method with several numerical examples. For high frequency problems, this approach outperforms traditional iterative methods such as Gauss-Newton initialized without any prior (i.e., initialized using a unit circle), or initialized using the solution of a direct method such as the linear sampling method. The algorithm remains robust to noise in the scattered field measurements and also converges to the true solution for limited aperture data. However, the number of training samples required to train the neural network scales exponentially in frequency and the complexity of the obstacles considered. We conclude with a discussion of this phenomenon and potential directions for future research.

We introduce an efficient numerical method for second order linear ODEs whose solution may vary between highly oscillatory and slowly changing over the solution interval. In oscillatory regions the solution is generated via a nonoscillatory phase function that obeys the nonlinear Riccati equation. We propose a defect-correction iteration that gives an asymptotic series for such a phase function; this is numerically approximated on a Chebyshev grid with a small number of nodes. For analytic coefficients we prove that each iteration, up to a certain maximum number, reduces the residual by a factor of order of the local frequency. The algorithm adapts both the step size and the choice of method, switching to a conventional spectral collocation method away from oscillatory regions. In numerical experiments we find that our proposal outperforms other state-of-the-art oscillatory solvers, most significantly at low-to-intermediate frequencies and at low tolerances, where it may use up to 106 times fewer function evaluations. Even in high frequency regimes, our implementation is on average 10 times faster than other specialized solvers.

Cryo-electron microscopy (cryo-EM) has recently become a premier method for obtaining high-resolution structures of biological macromolecules. However, it is limited to biomolecular samples with low conformational heterogeneity, where all the conformations can be well-sampled at many projection angles. While cryo-EM technically provides single-molecule data for heterogeneous molecules, most existing reconstruction tools cannot extract the full distribution of possible molecular configurations. To overcome these limitations, we build on a prior Bayesian approach and develop an ensemble refinement framework that estimates the ensemble density from a set of cryo-EM particles by reweighting a prior ensemble of conformations, e.g., from molecular dynamics simulations or structure prediction tools. Our work is a general approach to recovering the equilibrium probability density of the biomolecule directly in conformational space from single-molecule data. To validate the framework, we study the extraction of state populations and free energies for a simple toy model and from synthetic cryo-EM images of a simulated protein that explores multiple folded and unfolded conformations.

Cryo-electron microscopy (cryo-EM) has recently become a premier method for obtaining high-resolution structures of biological macromolecules. However, it is limited to biomolecular samples with low conformational heterogeneity, where all the conformations can be well-sampled at many projection angles. While cryo-EM technically provides single-molecule data for heterogeneous molecules, most existing reconstruction tools cannot extract the full distribution of possible molecular configurations. To overcome these limitations, we build on a prior Bayesian approach and develop an ensemble refinement framework that estimates the ensemble density from a set of cryo-EM particles by reweighting a prior ensemble of conformations, e.g., from molecular dynamics simulations or structure prediction tools. Our work is a general approach to recovering the equilibrium probability density of the biomolecule directly in conformational space from single-molecule data. To validate the framework, we study the extraction of state populations and free energies for a simple toy model and from synthetic cryo-EM images of a simulated protein that explores multiple folded and unfolded conformations.

Recent experiments of chemical reactions in optical cavities have shown great promise to alter and steer chemical reactions but still remain poorly understood theoretically. In particular the origin of resonant effects between the cavity and certain vibrational modes in the collective limit is still subject to active research. In this paper, we study unimolecular dissociation reactions of many molecules collectively interacting with an infrared cavity mode through their vibrational dipole moment. We find that the reaction rate can slow down by increasing the number of aligned molecules if the cavity mode is resonant with a vibrational frequency of the molecules. We also discover a simple scaling relation that scales with the collective Rabi splitting to estimate the onset of reaction rate modification by collective vibrational strong coupling and numerically demonstrate these effects for up to 10,000 molecules.

We examine the ground-state phase diagram and thermal phase transitions in a plaquettized fully frustrated bilayer spin-1/2 Heisenberg model. Based on a combined analysis from sign-problem free quantum Monte Carlo simulations, perturbation theory and free-energy arguments, we identify a first-order quantum phase transition line that separates two competing quantum-disordered ground states with dominant singlet formations on inter-layer dimers and plaquettes, respectively. At finite temperatures, this line extends to form a wall of first-order thermal transitions, which terminates in a line of thermal critical points. From a perturbative approach in terms of an effective Ising model description, we identify a quadratic suppression of the critical temperature scale in the strongly plaquettized region. Based on free-energy arguments we furthermore obtain the full phase boundary of the low-temperature dimer-singlet regime, which agrees well with the quantum Monte Carlo data.

We examine the ground-state phase diagram and thermal phase transitions in a plaquettized fully frustrated bilayer spin-1/2 Heisenberg model. Based on a combined analysis from sign-problem free quantum Monte Carlo simulations, perturbation theory and free-energy arguments, we identify a first-order quantum phase transition line that separates two competing quantum-disordered ground states with dominant singlet formations on inter-layer dimers and plaquettes, respectively. At finite temperatures, this line extends to form a wall of first-order thermal transitions, which terminates in a line of thermal critical points. From a perturbative approach in terms of an effective Ising model description, we identify a quadratic suppression of the critical temperature scale in the strongly plaquettized region. Based on free-energy arguments we furthermore obtain the full phase boundary of the low-temperature dimer-singlet regime, which agrees well with the quantum Monte Carlo data.

Recent experiments of chemical reactions in optical cavities have shown great promise to alter and steer chemical reactions but still remain poorly understood theoretically. In particular the origin of resonant effects between the cavity and certain vibrational modes in the collective limit is still subject to active research. In this paper, we study unimolecular dissociation reactions of many molecules collectively interacting with an infrared cavity mode through their vibrational dipole moment. We find that the reaction rate can slow down by increasing the number of aligned molecules if the cavity mode is resonant with a vibrational frequency of the molecules. We also discover a simple scaling relation that scales with the collective Rabi splitting to estimate the onset of reaction rate modification by collective vibrational strong coupling and numerically demonstrate these effects for up to 10,000 molecules.

Theoretical research into many-body quantum systems has mostly focused on regular structures which have a small, simple unit cell and where a vanishingly small number of pairs of the constituents directly interact. Motivated by advances in control over the pairwise interactions in many-body simulators, we determine the fate of spin systems on more general, arbitrary graphs. Placing the minimum possible constraints on the underlying graph, we prove how, with certainty in the thermodynamic limit, such systems behave like a single collective spin. We thus understand the emergence of complex many-body physics as dependent on `exceptional', geometrically constrained structures such as the low-dimensional, regular ones found in nature. Within the space of dense graphs we identify hitherto unknown exceptions via their inhomogeneity and observe how complexity is heralded in these systems by entanglement and highly non-uniform correlation functions. Our work paves the way for the discovery and exploitation of a whole class of geometries which can host uniquely complex phases of matter.