Publications

We study the dynamics of charge transfer insulators after photoexcitation using the three-band Emery model and a nonequilibrium extension of
Hartree
−
Fock
+
EDMFT
(extended dynamical mean field theory) and
GW
+
EDMFT
. While the equilibrium properties are accurately reproduced by the Hartree-Fock treatment of the ligand bands, dynamical correlations are essential for a proper description of the photodoped state. Photodoping leads to a renormalization of the charge transfer gap and to a substantial broadening of the bands. We calculate the time-resolved photoemission spectrum and optical conductivity and find qualitative agreement with experiments. Our formalism enables the realistic description of nonequilibrium phenomena in materials with ligand bands. It provides a tool to explore the optical manipulation of interaction and correlation effects, including insulator-metal and magnetic transitions.

We present a statistical learning framework for robust identification of partial differential equations from noisy spatiotemporal data. Extending previous sparse regression approaches for inferring PDE models from simulated data, we address key issues that have thus far limited the application of these methods to noisy experimental data, namely their robustness against noise and the need for manual parameter tuning. We address both points by proposing a stability-based model selection scheme to determine the level of regularization required for reproducible recovery of the underlying PDE. This avoids manual parameter tuning and provides a principled way to improve the method's robustness against noise in the data. Our stability selection approach, termed PDE-STRIDE, can be combined with any sparsity-promoting penalized regression model and provides an interpretable criterion for model component importance. We show that in particular the combination of stability selection with the iterative hard-thresholding algorithm from compressed sensing provides a fast, parameter-free, and robust computational framework for PDE inference that outperforms previous algorithmic approaches with respect to recovery accuracy, amount of data required, and robustness to noise. We illustrate the performance of our approach on a wide range of noise-corrupted simulated benchmark problems, including 1D Burgers, 2D vorticity-transport, and 3D reaction-diffusion problems. We demonstrate the practical applicability of our method on real-world data by considering a purely data-driven re-evaluation of the advective triggering hypothesis for an embryonic polarization system in C. elegans. Using fluorescence microscopy images of C. elegans zygotes as input data, our framework is able to recover the PDE model for the regulatory reaction-diffusion-flow network of the associated proteins.

The Simons Observatory (SO) is a ground-based cosmic microwave background (CMB) experiment sited on Cerro Toco in the Atacama Desert in Chile that promises to provide breakthrough discoveries in fundamental physics, cosmology, and astrophysics. Supported by the Simons Foundation, the Heising-Simons Foundation, and with contributions from collaborating institutions, SO will see first light in 2021 and start a five year survey in 2022. SO has 287 collaborators from 12 countries and 53 institutions, including 85 students and 90 postdocs.
The SO experiment in its currently funded form ('SO-Nominal') consists of three 0.4 m Small Aperture Telescopes (SATs) and one 6 m Large Aperture Telescope (LAT). Optimized for minimizing systematic errors in polarization measurements at large angular scales, the SATs will perform a deep, degree-scale survey of 10% of the sky to search for the signature of primordial gravitational waves. The LAT will survey 40% of the sky with arc-minute resolution. These observations will measure (or limit) the sum of neutrino masses, search for light relics, measure the early behavior of Dark Energy, and refine our understanding of the intergalactic medium, clusters and the role of feedback in galaxy formation.
With up to ten times the sensitivity and five times the angular resolution of the Planck satellite, and roughly an order of magnitude increase in mapping speed over currently operating ("Stage 3") experiments, SO will measure the CMB temperature and polarization fluctuations to exquisite precision in six frequency bands from 27 to 280 GHz. SO will rapidly advance CMB science while informing the design of future observatories such as CMB-S4.

As we enter a new era of quantum technology, it is increasingly important to develop methods to aid in the accurate preparation of quantum states for a variety of materials, matter, and devices. Computational techniques can be used to reconstruct a state from data, however the growing number of qubits demands ongoing algorithmic advances in order to keep pace with experiments. In this paper, we present an open-source software package called QuCumber that uses machine learning to reconstruct a quantum state consistent with a set of projective measurements. QuCumber uses a restricted Boltzmann machine to efficiently represent the quantum wavefunction for a large number of qubits. New measurements can be generated from the machine to obtain physical observables not easily accessible from the original data.

In the present work, we introduce a Self-Consistent Density-Functional Embedding technique, which leaves the realm of standard energy-functional approaches in Density Functional Theory and targets directly the density-to-potential mapping that lies at its heart. Inspired by the Density Matrix Embedding Theory, we project the full system onto a set of small interacting fragments that can be solved accurately. Based on the rigorous relation of density and potential in Density Functional Theory, we then invert the fragment densities to local potentials. Combining these results in a continuous manner provides an update for the Kohn-Sham potential of the full system, which is then used to update the projection. The scheme proposed here converges to an accurate approximation for the density and the Kohn-Sham potential of the full system. Convergence to exact results can be achieved by increasing the fragment size. We find, however, that already for small embedded fragments accurate results are obtained. We benchmark our approach for molecular bond stretching in one and two dimensions and demonstrate that it reproduces the known steps and peaks that are present in the exact exchange-correlation potential with remarkable accuracy.

In time series classification and regression, signals are typically mapped into some intermediate representation used for constructing models. Since the underlying task is often insensitive to time shifts, these representations are required to be time-shift invariant. We introduce the joint time-frequency scattering transform, a time-shift invariant representation which characterizes the multiscale energy distribution of a signal in time and frequency. It is computed through wavelet convolutions and modulus non-linearities and may therefore be implemented as a deep convolutional neural network whose filters are not learned but calculated from wavelets. We consider the progression from mel-spectrograms to time scattering and joint time-frequency scattering transforms, illustrating the relationship between increased discriminability and refinements of convolutional network architectures. The suitability of the joint time-frequency scattering transform for time-shift invariant characterization of time series is demonstrated through applications to chirp signals and audio synthesis experiments. The proposed transform also obtains state-of-the-art results on several audio classification tasks, outperforming time scattering transforms and achieving accuracies comparable to those of fully learned networks.

The Hall coefficient RH of Sr2RuO4 exhibits a non-monotonic temperature dependence with two sign reversals. We show that this puzzling behavior is the signature of two crossovers, which are key to the physics of this material. The increase of RH and the first sign change upon cooling are associated with a crossover into a regime of coherent quasiparticles with strong orbital differentiation of the inelastic scattering rates. The eventual decrease and the second sign change at lower temperature are driven by the crossover from inelastic to impurity-dominated scattering. This qualitative picture is supported by quantitative calculations of RH(T) using the Boltzmann transport theory in combination with dynamical mean-field theory, taking into account the effect of spin–orbit coupling. Our insights shed new light on the temperature dependence of the Hall coefficient in materials with strong orbital differentiation, as observed in Hund’s metals.

We present the results of a search for short-duration gravitational-wave transients in the data from the second observing run of Advanced LIGO and Advanced Virgo. We search for gravitational-wave transients with a duration of milliseconds to approximately one second in the 32-4096 Hz frequency band with minimal assumptions about the signal properties, thus targeting a wide variety of sources. We also perform a matched-filter search for gravitational-wave transients from cosmic string cusps for which the waveform is well-modeled. The unmodeled search detected gravitational waves from several binary black hole mergers which have been identified by previous analyses. No other significant events have been found by either the unmodeled search or the cosmic string search. We thus present search sensitivity for a variety of signal waveforms and report upper limits on the source rate-density as function of the characteristic frequency of the signal. These upper limits are a factor of three lower than the first observing run, with a 50% detection probability for gravitational-wave emissions with energies of ∼10−9M⊙c2 at 153 Hz. For the search dedicated to cosmic string cusps we consider several loop distribution models, and present updated constraints from the same search done in the first observing run.

The importance of the Decadal Survey in astrophysics is great; it deserves attention and revision. We make recommendations to increase the Survey's transparency and political legitimacy. The Astro2020 charge asks the Survey to "generate consensus recommendations". It is healthy to re-evaluate how to achieve consensus as the community and context evolve. Our recommendations are the following: (R1) Appoint the Decadal panel chairs and panel members through a transparent process, or even a democratic process. (R2) Don't make panel members sign any kinds of non-disclosure agreements, or strictly limit these. (R3) Educate the community about the Decadal's decision-making and consensus-building. (R4) Provide written documentation about how white papers will be read and used. (R5) Give the community an opportunity to comment on and vote to approve the final reports. (R6) Ask the AAAC to help the agencies make these changes.

We develop the first algorithm able to jointly compute the maximum {\it a posteriori} estimate of the Cosmic Microwave Background (CMB) temperature and polarization fields, the gravitational potential by which they are lensed, and cosmological parameters such as the tensor-to-scalar ratio, r. This is an important step towards sampling from the joint posterior probability function of these quantities, which, assuming Gaussianity of the CMB fields and lensing potential, contains all available cosmological information and would yield theoretically optimal constraints. Attaining such optimal constraints will be crucial for next-generation CMB surveys like CMB-S4, where limits on r could be improved by factors of a few over currently used sub-optimal quadratic estimators. The maximization procedure described here depends on a newly developed lensing algorithm, which we term \textsc{LenseFlow}, and which lenses a map by solving a system of ordinary differential equations. This description has conceptual advantages, such as allowing us to give a simple non-perturbative proof that the lensing determinant is equal to unity in the weak-lensing regime. The algorithm itself maintains this property even on pixelized maps, which is crucial for our purposes and unique to \textsc{LenseFlow} as compared to other lensing algorithms we have tested. It also has other useful properties such as that it can be trivially inverted (i.e. delensing) for the same computational cost as the forward operation, and can be used to compute lensing adjoint, Jacobian, and Hessian operators. We test and validate the maximization procedure on flat-sky simulations covering up to 600\,deg2 with non-uniform noise and masking.