Publications

Observations indicate that nearly all galaxies contain supermassive black holes (SMBHs) at their centers. When galaxies merge, their component black holes form SMBH binaries (SMBHBs), which emit low-frequency gravitational waves (GWs) that can be detected by pulsar timing arrays (PTAs). We have searched the recently-released North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 11-year data set for GWs from individual SMBHBs in circular orbits. As we did not find strong evidence for GWs in our data, we placed 95\% upper limits on the strength of GWs from such sources as a function of GW frequency and sky location. We placed a sky-averaged upper limit on the GW strain of $\mathith_0<7.3(3)×10^−15 at \mathitf_gw=8 nHz$. We also developed a technique to determine the significance of a particular signal in each pulsar using ``dropout' parameters as a way of identifying spurious signals in measurements from individual pulsars. We used our upper limits on the GW strain to place lower limits on the distances to individual SMBHBs. At the most-sensitive sky location, we ruled out SMBHBs emitting GWs with fgw=8 nHz within 120 Mpc for =109M⊙, and within 5.5 Gpc for =1010M⊙. We also determined that there are no SMBHBs with >1.6×109M⊙ emitting GWs in the Virgo Cluster. Finally, we estimated the number of potentially detectable sources given our current strain upper limits based on galaxies in Two Micron All-Sky Survey (2MASS) and merger rates from the Illustris cosmological simulation project. Only 34 out of 75,000 realizations of the local Universe contained a detectable source, from which we concluded it was unsurprising that we did not detect any individual sources given our current sensitivity to GWs.

Improved distance measurements to millisecond pulsars can enhance pulsar timing array (PTA) sensitivity to gravitational waves, improve tests of general relativity with binary pulsars, improve constraints on fuzzy dark matter, and more. Here we report the parallax distance measurements to six Gaia DR2 objects associated with International PTA pulsars J0437-4715, J1012+5307, J1024-0719, J1732-5049, J1910+1256, and J1843-1113. By multiplying the posteriors of the PTA distances with the \gaia distance to the companion, we improve the distance measurements, and provide a tentative detection of a previously unknown binary companion to J1843-1113. Finally, we recommend that future Gaia data releases use J0437-4715 as a calibration point, since its distance estimate in Gaia DR2 is relatively poor compared to pulsar timing measurements.

Mechanical strain is a powerful technique for tuning electronic structure and interactions in quantum materials. In a system with tetragonal symmetry, a tunable uniaxial in-plane strain can be used to probe nematic correlations in the same way that a tunable magnetic field is used to probe magnetic correlations. Here, we present a new spectroscopic scanned probe technique that provides atomic-resolution insight into the effect of anisotropic strain on the electronic structure. We use this technique to study nematic fluctuations and nematic order across the phase diagram of a prototypical iron-based superconductor. By extracting quantitatively the electronic anisotropy as function of applied strain, we show that while true long range nematic order is established at the tetragonal to orthorhombic structural transition temperature, sizable nematic fluctuations persist to high temperatures and also to the overdoped end of the superconducting dome. Remarkably, we find that uniaxial strain in the pnictides significantly enhances the amplitude of the nematic fluctuations, indicating a strong nonlinear coupling between structure and electronic nematicity.

We present a novel representation of coupled matter-photon systems that allows the application of many-body methods developed for purely fermionic systems. We do so by rewriting the original coupled light-matter problem in a higher-dimensional configuration space and then use photon-dressed orbitals as a basis to expand the thus "fermionized" coupled system. As an application we present a dressed time-dependent density-functional theory approach. The resulting dressed Kohn-Sham scheme allows for straightforward non-adiabatic approximations to the unknown exchange-correlation potential that explicitly includes correlations. We illustrate this for simple model systems placed inside a high-Q optical cavity, and show also results for observables such as the photon-field fluctuations that are hard to capture in standard matter-photon Kohn-Sham. We finally highlight that the dressed-orbital approach extends beyond the context of cavity quantum electrodynamics and can be applied to, e.g., van-der-Waals problems.

We present a fluctuating boundary integral method (FBIM) for overdamped Brownian Dynamics (BD) of two-dimensional periodic suspensions of rigid particles of complex shape immersed in a Stokes fluid. We develop a novel approach for generating Brownian displacements that arise in response to the thermal fluctuations in the fluid. Our approach relies on a first-kind boundary integral formulation of a mobility problem in which a random surface velocity is prescribed on the particle surface, with zero mean and covariance proportional to the Green's function for Stokes flow (Stokeslet). This approach yields an algorithm that scales linearly in the number of particles for both deterministic and stochastic dynamics, handles particles of complex shape, achieves high order of accuracy, and can be generalized to three dimensions and other boundary conditions. We show that Brownian displacements generated by our method obey the discrete fluctuation-dissipation balance relation (DFDB). Based on a recently-developed Positively Split Ewald method [A. M. Fiore, F. Balboa Usabiaga, A. Donev and J. W. Swan, J. Chem. Phys., 146, 124116, 2017], near-field contributions to the Brownian displacements are efficiently approximated by iterative methods in real space, while far-field contributions are rapidly generated by fast Fourier-space methods based on fluctuating hydrodynamics. FBIM provides the key ingredient for time integration of the overdamped Langevin equations for Brownian suspensions of rigid particles. We demonstrate that FBIM obeys DFDB by performing equilibrium BD simulations of suspensions of starfish-shaped bodies using a random finite difference temporal integrator.

The brain is a massive neuronal network, organized into anatomically distributed sub-circuits, with functionally relevant activity occurring at timescales ranging from milliseconds to years. Current methods to monitor neural activity, however, lack the necessary conjunction of anatomical spatial coverage, temporal resolution, and long-term stability to measure this distributed activity. Here we introduce a large-scale, multi-site, extracellular recording platform that integrates polymer electrodes with a modular stacking headstage design supporting up to 1,024 recording channels in freely behaving rats. This system can support months-long recordings from hundreds of well-isolated units across multiple brain regions. Moreover, these recordings are stable enough to track large numbers of single units for over a week. This platform enables large-scale electrophysiological interrogation of the fast dynamics and long-timescale evolution of anatomically distributed circuits, and thereby provides a new tool for understanding brain activity.

In suspensions of microorganisms, pattern formation can arise from the interplay of chemotaxis and the fluid flows collectively generated by the organisms themselves. Here we investigate the resulting pattern formation in square and elongated domains in the context of two distinct models of locomotion in which the chemoattractant dynamics is fully coupled to the fluid flows and swimmer motion. Analyses for both models reveal an aggregative instability due to chemotaxis, independent of swimmer shape and type, and a hydrodynamic instability for “pusher” swimmers. We discuss the similarities and differences between the models. Simulations reveal a critical length scale of the swimmer aggregates and this feature can be utilized to stabilize swimmer concentration patterns into quasi-one-dimensional bands by varying the domain size. These concentration bands transition to traveling pulses under an external chemoattractant gradient, as observed in experiments with chemotactic bacteria.

A key unmet challenge in interpreting omics experiments is inferring biological meaning in the context of public functional genomics data. We developed a computational framework, Your Evidence Tailored Integration (YETI; http://yeti.princeton.edu/ ), which creates specialized functional interaction maps from large public datasets relevant to an individual omics experiment. Using this tailored integration, we predicted and experimentally confirmed an unexpected divergence in viral replication after seasonal or pandemic human influenza virus infection.

For smooth bounded domains in ℝn, we prove upper and lower L2 bounds on the boundary data of Neumann eigenfunctions, and we prove quasiorthogonality of this boundary data in a spectral window. The bounds are tight in the sense that both are independent of the eigenvalues; this is achieved by working with an appropriate norm for boundary functions, which includes a spectral weight, that is, a function of the boundary Laplacian. This spectral weight is chosen to cancel concentration at the boundary that can happen for whispering gallery-type eigenfunctions. These bounds are closely related to wave equation estimates due to Tataru. Using this, we bound the distance from an arbitrary Helmholtz parameter E>0 to the nearest Neumann eigenvalue in terms of boundary normal derivative data of a trial function u solving the Helmholtz equation (Δ−E)u=0. This inclusion bound improves over previously known bounds by a factor of E5/6, analogously to a recently improved inclusion bound in the Dirichlet case due to the first two authors. Finally, we apply our theory to present an improved numerical implementation of the method of particular solutions for computation of Neumann eigenpairs on smooth planar domains. We show that the new inclusion bound improves the relative accuracy in a computed Neumann eigenvalue (around the 42000th) from nine to fourteen digits, with negligible extra numerical effort.

Pulsar timing array (PTA) collaborations in North America, Australia, and Europe, have been exploiting the exquisite timing precision of millisecond pulsars over decades of observations to search for correlated timing deviations induced by gravitational waves (GWs). The discovery space of this nanohertz band is potentially rich with populations of inspiraling supermassive black-holes binaries, decaying cosmic string networks, relic post-inflation GWs, and even non-GW imprints of axionic dark matter.
This article aims to provide an understanding of the exciting open science questions in cosmology, galaxy evolution, and fundamental physics that will be addressed by the detection and study of GWs at nanohertz frequencies. The focus of the article is on providing an understanding of the mechanisms by which PTAs can address specific questions in these fields, and to outline some of the subtleties and difficulties in each case. The material included is weighted most heavily towards the questions which we expect will be answered in the coming months to decades with PTAs; however, we have made efforts to include most currently anticipated applications of nanohertz GWs.