Publications

The rotation periods of planet-hosting stars can be used for modeling and mitigating the impact of magnetic activity in radial velocity measurements and can help constrain the high-energy flux environment and space weather of planetary systems. Millions of stars and thousands of planet hosts are observed with the Transiting Exoplanet Survey Satellite (TESS). However, most will only be observed for 27 contiguous days in a year, making it difficult to measure rotation periods with traditional methods. This is especially problematic for field M dwarfs, which are ideal candidates for exoplanet searches, but which tend to have periods in excess of the 27-day observing baseline. We present a new tool, Astraea, for predicting long rotation periods from short-duration light curves combined with stellar parameters from Gaia DR2. Using Astraea, we can predict the rotation periods from Kepler 4-year light curves with 13% uncertainty overall (and a 9% uncertainty for periods > 30 days). By training on 27-day Kepler light curve segments, Astraea can predict rotation periods up to 150 days with 9% uncertainty (5% for periods > 30 days). After training this tool on these 27-day Kepler light curve segments, we applied \texttt{Astraea} to real TESS data. For the 195 stars that were observed by both Kepler and TESS, we were able to predict the rotation periods with 55% uncertainty despite the wild differences in systematics.

Gravitational waves provide a unique tool for observational astronomy. While the first LIGO--Virgo catalogue of gravitational-wave transients (GWTC-1) contains eleven signals from black hole and neutron star binaries, the number of observations is increasing rapidly as detector sensitivity improves. To extract information from the observed signals, it is imperative to have fast, flexible, and scalable inference techniques. In a previous paper, we introduced BILBY: a modular and user-friendly Bayesian inference library adapted to address the needs of gravitational-wave inference. In this work, we demonstrate that BILBY produces reliable results for simulated gravitational-wave signals from compact binary mergers, and verify that it accurately reproduces results reported for the eleven GWTC-1 signals. Additionally, we provide configuration and output files for all analyses to allow for easy reproduction, modification, and future use. This work establishes that BILBY is primed and ready to analyse the rapidly growing population of compact binary coalescence gravitational-wave signals.

Registration of histology images from multiple sources is a pressing problem in large-scale studies of spatial -omics data. Researchers often perform “common coordinate registration,” akin to segmentation, in which samples are partitioned based on tissue type to allow for quantitative comparison of similar regions across samples. Accuracy in such registration requires both high image resolution and global awareness, which mark a difficult balancing act for contemporary deep learning architectures. We present a novel convolutional neural network (CNN) architecture that combines (1) a local classification CNN that extracts features from image patches sampled sparsely across the tissue surface, and (2) a global segmentation CNN that operates on these extracted features. This hybrid network can be trained in an end-to-end manner, and we demonstrate its relative merits over competing approaches on a reference histology dataset as well as two published spatial transcriptomics datasets. We believe that this paradigm will greatly enhance our ability to process spatial -omics data, and has general purpose applications for the processing of high-resolution histology images on commercially available GPUs.

Most cognitive processes are studied using abstract or synthetic stimuli with specific features to fully control what is presented to subjects. However, recent studies have revealed enhancements of cognitive capacities (such as working memory) when processing naturalistic versus abstract stimuli. Using abstract stimuli constructed from distinct visual features (e.g., color and shape), we have recently shown that human subjects can learn multidimensional stimulus-reward associations via initially estimating reward value of individual features (feature-based learning) before gradually switching to learning about reward value of individual stimuli (object-based learning). Here, we examined whether similar strategies are adopted during learning about naturalistic stimuli that are clearly perceived as objects (instead of a combination of features) and contain both task-relevant and irrelevant features. We found that similar to learning about abstract stimuli, subjects initially adopted feature-based learning more strongly before transitioning to object-based learning. However, there were three key differences between learning about naturalistic and abstract stimuli. First, compared with abstract stimuli, the initial learning strategy was less feature-based for naturalistic stimuli. Second, subjects transitioned to object-based learning faster for naturalistic stimuli. Third, unexpectedly, subjects were more likely to adopt feature-based learning for naturalistic stimuli, both at the steady state and overall. These results suggest that despite the stronger tendency to perceive naturalistic stimuli as objects, which leads to greater likelihood of using object-based learning as the initial strategy and a faster transition to object-based learning, the influence of individual features on learning is stronger for these stimuli such that ultimately the object-based strategy is adopted less. Overall, our findings suggest that feature-based learning is a general initial strategy for learning about reward value of all types of multi-dimensional stimuli.

Living matter moves, deforms, and organizes itself. In cells this is made possible by networks of polymer filaments and crosslinking molecules that connect filaments to each other and that act as motors to do mechanical work on the network. For the case of highly cross-linked filament networks, we discuss how the material properties of assemblies emerge from the forces exerted by microscopic agents. First, we introduce a phenomenological model that characterizes the forces that crosslink populations exert between filaments. Second, we derive a theory that predicts the material properties of highly crosslinked filament networks, given the crosslinks present. Third, we discuss which properties of crosslinks set the material properties and behavior of highly crosslinked cytoskeletal networks. The work presented here, will enable the better understanding of cytoskeletal mechanics and its molecular underpinnings. This theory is also a first step towards a theory of how molecular perturbations impact cytoskeletal organization, and provides a framework for designing cytoskeletal networks with desirable properties in the lab.

Multitaper estimators have enjoyed significant success in estimating spectral densities from finite samples using as tapers Slepian functions defined on the acquisition domain. Unfortunately, the numerical calculation of these Slepian tapers is only tractable for certain symmetric domains, such as rectangles or disks. In addition, no performance bounds are currently available for the mean squared error of the spectral density estimate. This situation is inadequate for applications such as cryo-electron microscopy, where noise models must be estimated from irregular domains with small sample sizes. We show that the multitaper estimator only depends on the linear space spanned by the tapers. As a result, Slepian tapers may be replaced by proxy tapers spanning the same subspace (validating the common practice of using partially converged solutions to the Slepian eigenproblem as tapers). These proxies may consequently be calculated using standard numerical algorithms for block diagonalization. We also prove a set of performance bounds for multitaper estimators on arbitrary domains. The method is demonstrated on synthetic and experimental datasets from cryo-electron microscopy, where it reduces mean squared error by a factor of two or more compared to traditional methods.

A generic flow distribution network typically does not deliver its load at a uniform rate across a service area, instead oversupplying regions near the nutrient source while leaving downstream regions undersupplied. In this work we demonstrate how a local adaptive rule coupling tissue growth with nutrient density results in a flow network that self-organizes to deliver nutrients uniformly. This geometric adaptive rule can be generalized and imported to mechanics-based adaptive models to address the effects spatial gradients in nutrients or growth factors in tissues.

A major obstacle to treating Alzheimer’s disease (AD) is our lack of understanding of the molecular mechanisms underlying selective neuronal vulnerability, a key characteristic of the disease. Here, we present a framework integrating high-quality neuron-type-specific molecular profiles across the lifetime of the healthy mouse, which we generated using bacTRAP, with postmortem human functional genomics and quantitative genetics data. We demonstrate human-mouse conservation of cellular taxonomy at the molecular level for neurons vulnerable and resistant in AD, identify specific genes and pathways associated with AD neuropathology, and pinpoint a specific functional gene module underlying selective vulnerability, enriched in processes associated with axonal remodeling, and affected by amyloid accumulation and aging. We have made all cell-type-specific profiles and functional networks available at http://alz.princeton.edu. Overall, our study provides a molecular framework for understanding the complex interplay between Aβ, aging, and neurodegeneration within the most vulnerable neurons in AD.

Hot spots, or plasmoids, which form due to magnetic reconnection in current sheets, are conjectured to power frequent X-ray and near-infrared flares from Sgr A*, the black hole in the center of our Galaxy. It is unclear how, where, and when current sheets form in black hole accretion disks. We perform axisymmetric general-relativistic resistive magnetohydrodynamics simulations to model reconnection and plasmoid formation in a range of accretion flows. Current sheets and plasmoids are ubiquitous features that form regardless of the initial magnetic field in the disk, the magnetization in the quasisteady-state phase of accretion, and the spin of the black hole. Within 10 Schwarzschild radii from the event horizon, we observe plasmoids forming, after which they can merge, grow to macroscopic scales of the order of a few Schwarzschild radii, and are ultimately advected along the jet’s sheath or into the disk. Large plasmoids are energized to relativistic temperatures via reconnection and contribute to the jet’s limb brightening. We find that only hot spots forming in magnetically arrested disks can potentially explain the energetics of Sgr A* flares. The flare period is determined by the reconnection rate, which we find to be between $0.01c$ and $0.03c$ in all cases, consistent with studies of reconnection in isolated Harris-type current sheets. We quantify magnetic dissipation and nonideal electric fields, which can efficiently inject nonthermal particles. We show that explicit resistivity allows for converged numerical solutions, such that the electromagnetic energy evolution and dissipation become independent of the grid scale for the extreme resolutions considered here.

With rheology applications in mind, we present a fast solver for the time-dependent effective viscosity of an infinite lattice containing one or more neutrally buoyant smooth rigid particles per unit cell, in a two-dimensional Stokes fluid with given shear rate. At each time, the mobility problem is reformulated as a 2nd-kind boundary integral equation, then discretized to spectral accuracy by the Nystrom method and solved iteratively, giving typically 10 digits of accuracy. Its solution controls the evolution of particle locations and angles in a first-order system of ordinary differential equations. The formulation is placed on a rigorous footing by defining a generalized periodic Green's function for the skew lattice. Numerically, the periodized integral operator is split into a near image sum|applied in linear time via the fast multipole method|plus a correction field solved cheaply via proxy Stokeslets. We use barycentric quadratures to evaluate particle interactions and velocity fields accurately, even at distances much closer than the node spacing. Using first-order time-stepping we simulate, eg, 25 ellipses per unit cell to 3-digit accuracy on a desktop in 1 hour per shear time. Our examples show equilibration at long times, force chains, and two types of blow-ups (jamming) whose power laws match lubrication theory asymptotics.