Publications

We have constructed a new timescale, TT(IPTA16), based on observations of radio pulsars presented in the first data release from the International Pulsar Timing Array (IPTA). We used two analysis techniques with independent estimates of the noise models for the pulsar observations and different algorithms for obtaining the pulsar timescale. The two analyses agree within the estimated uncertainties and both agree with TT(BIPM17), a post-corrected timescale produced by the Bureau International des Poids et Mesures (BIPM). We show that both methods could detect significant errors in TT(BIPM17) if they were present. We estimate the stability of the atomic clocks from which TT(BIPM17) is derived using observations of four rubidium fountain clocks at the US Naval Observatory. Comparing the power spectrum of TT(IPTA16) with that of these fountain clocks suggests that pulsar-based timescales are unlikely to contribute to the stability of the best timescales over the next decade, but they will remain a valuable independent check on atomic timescales. We also find that the stability of the pulsar-based timescale is likely to be limited by our knowledge of solar-system dynamics, and that errors in TT(BIPM17) will not be a limiting factor for the primary goal of the IPTA, which is to search for the signatures of nano-Hertz gravitational waves.

Cosmological studies of large-scale structure have relied on two-point statistics, not fully exploiting the rich structure of the cosmic web. In this paper we show how to capture some of this cosmic web information by using the minimum spanning tree (MST), for the first time using it to estimate cosmological parameters in simulations. Discrete tracers of dark matter such as galaxies, N-body particles or haloes are used as nodes to construct a unique graph, the MST, that traces skeletal structure. We study the dependence of the MST on cosmological parameters using haloes from a suite of COLA simulations with a box size of 250 h−1Mpc, varying the amplitude of scalar fluctuations (As), matter density (Ωm), and neutrino mass (∑mν). The power spectrum P and bispectrum B are measured for wavenumbers between 0.125 and 0.5 hMpc−1, while a corresponding lower cut of ∼12.6 h−1Mpc is applied to the MST. The constraints from the individual methods are fairly similar but when combined we see improved 1σ constraints of ∼17% (∼12%) on Ωm and ∼12% (∼10%) on As with respect to P (P+B) thus showing the MST is providing additional information. The MST can be applied to current and future spectroscopic surveys (BOSS, DESI, Euclid, PSF, WFIRST, and 4MOST) in 3D and photometric surveys (DES and LSST) in tomographic shells to constrain parameters and/or test systematics.

Fish schools and bird flocks exhibit complex collective dynamics whose self-organization principles are largely unknown. The influence of hydrodynamics on such collectives has been relatively unexplored theoretically, in part due to the difficulty in modeling the temporally long-lived hydrodynamic interactions between many dynamic bodies. We address this through a novel discrete-time dynamical system (iterated map) that describes the hydrodynamic interactions between flapping swimmers arranged in one- and two-dimensional lattice formations. Our 1D results exhibit good agreement with previously published experimental data, in particular predicting the bistability of schooling states and new instabilities that can be probed in experimental settings. For 2D lattices, we determine the formations for which swimmers optimally benefit from hydrodynamic interactions. We thus obtain the following hierarchy: while a side-by-side single-row “phalanx” formation offers a small improvement over a solitary swimmer, 1D in-line and 2D rectangular lattice formations exhibit substantial improvements, with the 2D diamond lattice offering the largest hydrodynamic benefit. Generally, our self-consistent modeling framework may be broadly applicable to active systems in which the collective dynamics is primarily driven by a fluid-mediated memory.

Inventors of the original artificial neural networks (ANNs) derived their inspiration from biology [1]. However, today, most ANNs, such as backpropagation-based convolutional deeplearning networks, resemble natural NNs only superficially. Given that, on some tasks, such ANNs achieve human or even superhuman performance, why should one care about such dissimilarity with natural NNs? The algorithms of natural NNs are relevant if one's goal is not just to outperform humans on certain tasks but to develop general-purpose artificial intelligence rivaling that of a human. As contemporary ANNs are far from achieving this goal and natural NNs, by definition, achieve it, natural NNs must contain some "secret sauce" that ANNs lack. This is why we need to understand the algorithms implemented by natural NNs.

The nonequilibrium dynamics of correlated many-particle systems is of interest in connection with pump-probe experiments on molecular systems and solids, as well as theoretical investigations of transport properties and relaxation processes. Nonequilibrium Green's functions are a powerful tool to study interaction effects in quantum many-particle systems out of equilibrium, and to extract physically relevant information for the interpretation of experiments. We present the open-source software package NESSi (The Non-Equilibrium Systems Simulation package) which allows to perform many-body dynamics simulations based on Green's functions on the L-shaped Kadanoff-Baym contour. NESSi contains the library libcntr which implements tools for basic operations on these nonequilibrium Green's functions, for constructing Feynman diagrams, and for the solution of integral and integro-differential equations involving contour Green's functions. The library employs a discretization of the Kadanoff-Baym contour into time N points and a high-order implementation of integration routines. The total integrated error scales up to (N−7), which is important since the numerical effort increases at least cubically with the simulation time. A distributed-memory parallelization over reciprocal space allows large-scale simulations of lattice systems. We provide a collection of example programs ranging from dynamics in simple two-level systems to problems relevant in contemporary condensed matter physics, including Hubbard clusters and Hubbard or Holstein lattice models. The libcntr library is the basis of a follow-up software package for nonequilibrium dynamical mean-field theory calculations based on strong-coupling perturbative impurity solvers.

The Wide-Field Infrared Space Telescope (WFIRST) will be capable of delivering precise astrometry for faint sources over the enormous field of view of its main camera, the Wide-Field Imager (WFI). This unprecedented combination will be transformative for the many scientific questions that require precise positions, distances, and velocities of stars. We describe the expectations for the astrometric precision of the WFIRST WFI in different scenarios, illustrate how a broad range of science cases will see significant advances with such data, and identify aspects of WFIRST’s design where small adjustments could greatly improve its power as an astrometric instrument.

We describe the symmetry protected nodal points that can exist in magnetic space groups and show that only three-, six-, and eightfold degeneracies are possible (in addition to the two- and fourfold degeneracies that have already been studied). The three- and sixfold degeneracies are derived from “spin-1” Weyl fermions. The eightfold degeneracies come in different flavors. In particular, we distinguish between eightfold fermions that realize nonchiral “Rarita-Schwinger fermions” and those that can be described as four degenerate Weyl fermions. We list the (magnetic and nonmagnetic) space groups where these exotic fermions can be found. We further show that in several cases, a magnetic translation symmetry pins the Hamiltonian of the multifold fermion to an idealized exactly solvable point that is not achievable in nonmagnetic crystals without fine-tuning. Finally, we present known compounds that may host these fermions and methods for systematically finding more candidate materials.

Emergent quantum phases driven by electronic interactions can manifest in materials with narrowly dispersing, i.e. "flat", energy bands. Recently, flat bands have been realized in a variety of graphene-based heterostructures using the tuning parameters of twist angle, layer stacking and pressure, and resulting in correlated insulator and superconducting states. Here we report the experimental observation of similar correlated phenomena in twisted bilayer tungsten diselenide (tWSe2), a semiconducting transition metal dichalcogenide (TMD). Unlike twisted bilayer graphene where the flat band appears only within a narrow range around a "magic angle", we observe correlated states over a continuum of angles, spanning 4 degree to 5.1 degree. A Mott-like insulator appears at half band filling that can be sensitively tuned with displacement field. Hall measurements supported by ab initio calculations suggest that the strength of the insulator is driven by the density of states at half filling, consistent with a 2D Hubbard model in a regime of moderate interactions. At 5.1 degree twist, we observe evidence of superconductivity upon doping away from half filling, reaching zero resistivity around 3 K. Our results establish twisted bilayer TMDs as a model system to study interaction-driven phenomena in flat bands with dynamically tunable interactions.

Lei Wang, En-Min Shih, Augusto Ghiotto, Lede Xian, Daniel A. Rhodes, Cheng Tan, Martin Claassen, Dante M. Kennes, Yusong Bai, Bumho Kim, Kenji Watanabe, Takashi Taniguchi, Xiaoyang Zhu, James Hone, Angel Rubio, Abhay Pasupathy, Cory R. Dean

We explore the behavior of micron-scale autophoretic Janus (Au / Pt) rods, having various Au / Pt length ratios, swimming near a wall in an imposed background flow. We find that their ability to robustly orient and move upstream, i.e., to rheotax, depends strongly on the Au / Pt ratio, which is easily tunable in synthesis. Numerical simulations of swimming rods actuated by a surface slip show a similar rheotactic tunability when varying the location of the surface slip versus surface drag. The slip location determines whether swimmers are pushers (rear actuated), pullers (front actuated), or in between. Our simulations and modeling show that pullers rheotax most robustly due to their larger tilt angle to the wall, which makes them responsive to flow gradients. Thus, rheotactic response infers the nature of difficult to measure flow fields of an active particle, establishes its dependence on swimmer type, and shows how Janus rods can be tuned for flow responsiveness.

We show that optical excitation of the Mott insulating phase of the one-dimensional Hubbard model can create a state possessing two of the hallmarks of superconductivity: a nonvanishing charge stiffness and long-ranged pairing correlation. By employing the exact diagonalization method, we find that the superposition of the η-pairing eigenstates preferentially induced by the optical pump exhibits a nonvanishing charge stiffness and a pairing correlation that decays very slowly with system size. We show that the charge stiffness is indeed directly associated with the η-pairing correlation in the Hubbard model. Our finding demonstrates that optical pumping can actually lead to superconducting-like properties on the basis of the η-pairing states.