Publications

The acoustic inverse obstacle scattering problem consists of determining the shape of a domain from measurements of the scattered far field due to some set of incident fields (probes). For a penetrable object with known sound speed, this can be accomplished by treating the boundary alone as an unknown curve. Alternatively, one can treat the entire object as unknown and use a more general volumetric representation, without making use of the known sound speed. Both lead to strongly nonlinear and nonconvex optimization problems for which recursive linearization provides a useful framework for numerical analysis. After extending our shape optimization approach developed earlier for impenetrable bodies, we carry out a systematic study of both methods and compare their performance on a variety of examples. Our findings indicate that the volumetric approach is more robust, even though the number of degrees of freedom is significantly larger. We conclude with a discussion of this phenomenon and potential directions for further research.

A new scheme is presented for imposing periodic boundary conditions on unit cells with arbitrary source distributions. We restrict our attention here to the Poisson, modified Helmholtz, Stokes and modified Stokes equations. The approach extends to the oscillatory equations of mathematical physics, including the Helmholtz and Maxwell equations, but we will address these in a companion paper, since the nature of the problem is somewhat different and includes the consideration of quasiperiodic boundary conditions and resonances. Unlike lattice sum-based methods, the scheme is insensitive to the unit cell's aspect ratio and is easily coupled to adaptive fast multipole methods (FMMs). Our analysis relies on classical “plane-wave” representations of the fundamental solution, and yields an explicit low-rank representation of the field due to all image sources beyond the first layer of neighboring unit cells. When the aspect ratio of the unit cell is large, our scheme can be coupled with the nonuniform fast Fourier transform (NUFFT) to accelerate the evaluation of the induced field. Its performance is illustrated with several numerical examples.

First principle theoretical modeling of out-of-equilibrium processes observed in attosecond pump-probe transient absorption spectroscopy (TAS) triggering pure electron dynamics remains a challenging task, specially for heavy elements and/or core excitations containing fingerprints of scalar and spin-orbit relativistic effects. To address this, we formulate a methodology for simulating TAS within the relativistic real-time time-dependent density functional theory (RT-TDDFT) framework, for both the valence and core energy regime. Especially for TAS, full four-component (4c) RT simulations are feasible but computationally demanding. Therefore, in addition to the 4c approach, we also introduce the atomic mean-field exact two-component (amfX2C) Hamiltonian accounting for one- and two-electron picture-change corrections within RT-TDDFT. amfX2C preserves the accuracy of the parent 4c method at a fraction of its computational cost. Finally, we apply the methodology to study valence and near L

LaRh

To define the topology of driven systems, recent works have proposed synthetic dimensions as a way to uncover the underlying parameter space of topological invariants. Using time as a synthetic dimension, together with a momentum dimension, gives access to a synthetic 2D Chern number. It is, however, still unclear how the synthetic 2D Chern number is related to the Chern number that is defined from a parametric variable that evolves with time. Here we show that in periodically driven systems in the adiabatic limit, the synthetic 2D Chern number is a multiple of the Chern number defined from the parametric variable. The synthetic 2D Chern number can thus be engineered via how the parametric variable evolves in its own space. We justify our claims by investigating Thouless pumping in two 1D tight-binding models, a three-site chain model and a two-1D-sliding-chains model. The present findings could be extended to higher dimensions and other periodically driven configurations.

Using the truncated Wigner approximation (TWA) we study quench dynamics of two-dimensional lattice systems consisting of interacting spinless fermions with potential disorder. First, we demonstrate that the semiclassical dynamics generally relaxes faster than the full quantum dynamics. We obtain this result by comparing the semiclassical dynamics with exact diagonalization and Lanczos propagation of one-dimensional chains. Next, exploiting the TWA capabilities of simulating large lattices, we investigate how the relaxation rates depend on the dimensionality of the studied system. We show that strongly disordered one-dimensional and two-dimensional systems exhibit a transient, logarithmic-in-time relaxation, which was recently established for one-dimensional chains. Such relaxation corresponds to the infamous 1/f-noise at strong disorder.

We report on the strong field ionization of H2 by a corotating two-color laser field. We measure the electron momentum distribution in coincidence with the kinetic energy release (KER) of dissociating hydrogen molecules. In addition to a characteristic half-moon structure, we observe a low-energy structure in the electron momentum distribution at a KER of about 3.5 eV. We speculate that the outgoing electron interacts with the molecular ion, despite the absence of classical recollisions under these conditions. Time-dependent density functional theory simulations support our conclusions.

Recently the SP (Stochastic Polyak step size) method has emerged as a competitive adaptive method for setting the step sizes of SGD. SP can be interpreted as a method specialized to interpolated models, since it solves the interpolation equations. SP solves these equation by using local linearizations of the model. We take a step further and develop a method for solving the interpolation equations that uses the local second-order approximation of the model. Our resulting method SP2 uses Hessian-vector products to speed-up the convergence of SP. Furthermore, and rather uniquely among second-order methods, the design of SP2 in no way relies on positive definite Hessian matrices or convexity of the objective function. We show SP2 is competitive both in experiments and in theory.
We show SP2 is very competitive on matrix completion, non-convex test problems and logistic regression. We also provide a convergence theory on sums-of-quadratics.

Drosophila oogenesis is a powerful and tractable model for studies of cell and developmental biology due to the multitude of well-characterized events in both germline and somatic cells, the ease of genetic manipulation in fruit flies, and the large number of egg chambers produced by each fly. Recent improvements in live imaging and ex vivo culturing protocols have enabled researchers to conduct more detailed, longer-term studies of egg chamber development, enabling insights into fundamental biological processes. Here, we present a protocol for dissection, culturing, and imaging of late-stage egg chambers to study intercellular and directional cytoplasmic flow during “nurse cell dumping.” This critical developmental process towards the latter stages of oogenesis (stages 10b/11) results in rapid growth of the oocyte and shrinkage of the nurse cells and is accompanied by dynamic changes in cell shape.

This is a handbook of simple proofs of the convergence of gradient and stochastic gradient descent type methods. We consider functions that are Lipschitz, smooth, convex, strongly convex, and/or Polyak-Łojasiewicz functions. Our focus is on