Publications

In this work, we provide a solution to the problem of computing collision stress in particle-tracking simulations. First, a formulation for the collision stress between particles is derived as an extension of the virial stress formula to general-shaped particles with uniform or non-uniform mass density. Second, we describe a collision-resolution algorithm based on geometric constraint minimization which eliminates the stiff pairwise potentials in traditional methods. The method is validated with a comparison to the equation of state of Brownian spherocylinders. Then we demonstrate the application of this method in several emerging problems of soft active matter.

Although core helium-burning red clump (RC) stars are faint at ultraviolet wavelengths, their ultraviolet (UV)–optical color is a unique and accessible probe of their physical properties. Using data from the Galaxy Evolution Explorer All Sky Imaging Survey, Gaia Data Release 2, and the Sloan Digital Sky Survey Apache Point Observatory Galactic Evolution Experiment (APOGEE) DR14 survey, we find that spectroscopic metallicity is strongly correlated with the location of an RC star in the UV–optical color–magnitude diagram. The RC has a wide spread in (NUV–G)0 color of over 4 mag compared to a 0.7 mag range in (GBP–GRP)0. We propose a photometric, dust-corrected, UV–optical (NUV–G)0 color–metallicity [Fe/H] relation using a sample of 5,175 RC stars from APOGEE. We show that this relation has a scatter of 0.16 dex and is easier to obtain for large, wide-field samples than for spectroscopic metallicities. Importantly, the effect may be comparable to the spread in RC color attributed to extinction in other studies.

Rare-earth nickelates exhibit a remarkable metal-insulator transition accompanied by a symmetry-lowering structural distortion. Using model considerations and first-principles calculations, we present a theory of this phase transition which reveals the key role of the coupling between electronic and lattice instabilities. We show that the transition is driven by the proximity to an instability towards electronic disproportionation which couples to a specific structural distortion mode, cooperatively driving the system into the insulating state. This allows us to identify two key control parameters of the transition: the susceptibility to electronic disproportionation and the stiffness of the lattice mode. We show that our findings can be rationalized in terms of a Landau theory involving two coupled order parameters, with general implications for transition-metal oxides.

In the study of strongly-correlated, many-electron systems, the Hubbard Kanamori (HK) model has emerged as one of the prototypes for transition metal oxide physics. The model is multi-band in nature and contains Hund's coupling terms, which have pronounced effects on metal-insulator transitions, high-temperature superconductivity, and other physical properties. In the following, we present a complete theoretical framework for treating the HK model using the ground state Auxiliary Field Quantum Monte Carlo (AFQMC) method and analyze its performance on few-band models whose parameters approximate those observed in ruthenate, rhodates, and other materials exhibiting Hund's physics. Unlike previous studies, the constrained path and phaseless approximations are used to respectively control the sign and phase problems, which enables high accuracy modeling of the HK model's ground state properties within parameter regimes of experimental interest. We demonstrate that, after careful consideration of the Hubbard-Stratonovich transformations and trial wave functions employed, relative errors in the energy of less than 1% can routinely be achieved for moderate to large values of the Hund's coupling constant. Crucially, our methodology also accurately predicts magnetic ordering and phase transitions. The results presented open the door to more predictive modeling of Hund's physics within a wide range of strongly-correlated materials using AFQMC.

Many subcellular structures contain large numbers of cytoskeletal filaments. Such assemblies underlie much of cell division, motility, signaling, metabolism, and growth. Thus, understanding cell biology requires understanding the properties of networks of cytoskeletal filaments. While there are well established disciplines in biology dedicated to studying isolated proteins — their structure (Structural Biology) and behaviors (Biochemistry) — it is much less clear how to investigate, or even just describe, the structure and behaviors of collections of cytoskeletal filaments. One approach is to use methodologies from Mechanics and Soft Condensed Matter Physics, which have been phenomenally successful in the domains where they have been traditionally applied. From this perspective, collections of cytoskeletal filaments are viewed as materials, albeit very complex, ‘active’ materials, composed of molecules which use chemical energy to perform mechanical work. A major challenge is to relate these material level properties to the behaviors of the molecular constituents. Here we discuss this materials perspective and review recent work bridging molecular and network scale properties of the cytoskeleton, focusing on the organization of microtubules by dynein as an illustrative example.

An important class of fluid-structure problems involve the dynamics of ordered arrays of immersed, flexible fibers. While specialized numerical methods have been developed to study fluid-fiber systems, they become infeasible when there are many, rather than a few, fibers present, nor do these methods lend themselves to analytical calculation. Here, we introduce a coarse-grained continuum model, based on local-slender body theory, for elastic fibers immersed in a viscous Newtonian fluid. It takes the form of an anisotropic Brinkman equation whose skeletal drag is coupled to elastic forces. This model has two significant benefits: (1) the density effects of the fibers in a suspension become analytically manifest, and (2) it allows for the rapid simulation of dense suspensions of fibers in regimes inaccessible to standard methods. As a first validation, without fitting parameters, we achieve very reasonable agreement with 3D Immersed Boundary simulations of a bed of anchored fibers bent by a shear flow. Secondly, we characterize the effect of density on the relaxation time of fiber beds under oscillatory shear, and find close agreement to results from full numerical simulations. We then study buckling instabilities in beds of fibers, using our model both numerically and analytically to understand the role of fiber density and the structure of buckling transitions. We next apply our model to study the flow-induced bending of inclined fibers in a channel, as has been recently studied as a flow rectifier, examining the nature of the internal flows within the bed, and the emergence of inhomogeneous permeability. Finally, we extend the method to study a simple model of metachronal waves on beds of actuated fibers, as a model for ciliary beds. Our simulations reproduce qualitatively the pumping action of coordinated waves of compression through the bed.

Dissipationless charge transport is one of the defining properties of superconductors, but the interplay between dimensionality and disorder in determining the onset of dissipation remains an open theoretical and experimental problem. Here, we present measurements of the dissipation phase diagrams of superconductors in the two-dimensional limit, layer by layer, down to a monolayer in the presence of temperature (T), magnetic field (B) and current (I) in 2H-NbSe2. Our results show that the phase diagram strongly depends on the thickness even in the two-dimensional limit. At four layers we can define a finite region in the I–B phase diagram where dissipationless transport exists at T = 0. At even smaller thicknesses, this region shrinks in area until in a monolayer it approaches a single point defined by T = B = I = 0. In applied field, we show that time-dependent Ginzburg–Landau simulations that describe dissipation by vortex motion qualitatively reproduce our experimental I–B phase diagram. Last, we show that by using non-local transport and time-dependent Ginzburg–Landau calculations that we can engineer charge flow and create phase boundaries between dissipative and dissipationless transport regions in a single sample, demonstrating control over non-equilibrium states of matter.

In this work, we examine the topological phases that can arise in triangular lattices with disconnected elementary band representations. We show that, although these phases may be “fragile” with respect to the addition of extra bands, their topological properties are manifest in certain nontrivial holonomies (Wilson loops) in the space of nontrivial bands. We introduce an eigenvalue index for fragile topology, and we show how a nontrivial value of this index manifests as the winding of a hexagonal Wilson loop; this remains true even in the absence of time-reversal or sixfold rotational symmetry. Additionally, when time-reversal and twofold rotational symmetry are present, we show directly that there is a protected nontrivial winding in more conventional Wilson loops. Crucially, we emphasize that these Wilson loops cannot change without closing a gap to the nontrivial bands. By studying the entanglement spectrum for the fragile bands, we comment on the relationship between fragile topology and the “obstructed atomic limit” of Bradlyn et al. [Nature (London) 547, 298 (2017)]. We conclude with some perspectives on topological matter beyond the K-theory classification.

We present a search for prompt gamma-ray counterparts to compact binary coalescence gravitational wave (GW) candidates from Advanced LIGO's first observing run (O1). As demonstrated by the multimessenger observations of GW170817/GRB 170817A, electromagnetic and GW observations provide complementary information about the astrophysical source, and in the case of weaker candidates, may strengthen the case for an astrophysical origin. Here we investigate low-significance GW candidates from the O1 compact binary coalescence searches using the Fermi Gamma-Ray Burst Monitor (GBM), leveraging its all sky and broad energy coverage. Candidates are ranked and compared to background to measure the significance. Those with false alarm rates (FARs) of less than 10−5 Hz (about one per day, yielding a total of 81 candidates) are used as the search sample for gamma-ray follow-up. No GW candidates were found to be coincident with gamma-ray transients independently identified by blind searches of the GBM data. In addition, GW candidate event times were followed up by a separate targeted search of GBM data. Among the resulting GBM events, the two with the lowest FARs were the gamma-ray transient GW150914-GBM presented in Connaughton et al. and a solar flare in chance coincidence with a GW candidate.