Publications

Large stellar surveys are revealing the chemodynamical structure of the Galaxy across a vast spatial extent. However, the many millions of low-resolution spectra observed to date are yet to be fully exploited. We employ The Cannon, a data-driven approach to estimating abundances, to obtain detailed abundances from low-resolution (R = 1800) LAMOST spectra, using the GALAH survey as our reference. We deliver five (for dwarfs) or six (for giants) estimated abundances representing five different nucleosynthetic channels, for 3.9 million stars, to a precision of 0.05 - 0.23 dex. Using wide binary pairs, we demonstrate that our abundance estimates provide chemical discriminating power beyond metallicity alone. We show the coverage of our catalogue with radial, azimuthal and dynamical abundance maps, and examine the neutron capture abundances across the disk and halo, which indicate different origins for the in-situ and accreted halo populations. LAMOST has near-complete Gaia coverage and provides an unprecedented perspective on chemistry across the Milky Way.

Animal embryos are patterned by a handful of highly conserved inductive signals. Yet, in most cases, it is unknown which pattern features (i.e., spatial gradients or temporal dynamics) are required to support normal development. An ideal experiment to address this question would be to “paint” arbitrary synthetic signaling patterns on “blank canvas” embryos to dissect their requirements. Here, we demonstrate exactly this capability by combining optogenetic control of Ras/extracellular signal-related kinase (ERK) signaling with the genetic loss of the receptor tyrosine-kinase-driven terminal signaling patterning in early Drosophila embryos. Blue-light illumination at the embryonic termini for 90 min was sufficient to rescue normal development, generating viable larvae and fertile adults from an otherwise lethal terminal signaling mutant. Optogenetic rescue was possible even using a simple, all-or-none light input that reduced the gradient of Erk activity and eliminated spatiotemporal differences in terminal gap gene expression. Systematically varying illumination parameters further revealed that at least three distinct developmental programs are triggered at different signaling thresholds and that the morphogenetic movements of gastrulation are robust to a 3-fold variation in the posterior pattern width. These results open the door to controlling tissue organization with simple optical stimuli, providing new tools to probe natural developmental processes, create synthetic tissues with defined organization, or directly correct the patterning errors that underlie developmental defects.

High order perturbation theory has seen an unexpected recent revival for controlled calculations of quantum many-body systems, even at strong coupling. We adapt integration methods using low-discrepancy sequences to this problem. They greatly outperform state-of-the-art diagrammatic Monte Carlo. In practical applications, we show speed-ups of several orders of magnitude with scaling as fast as $1/N$ in sample number $N$; parametrically faster than $1/\sqrt{N}$ in Monte Carlo. We illustrate our technique with a solution of the Kondo ridge in quantum dots, where it allows large parameter sweeps.

We study the superconducting pairing correlations in the ground state of the doped Hubbard model -- in its original form without hopping beyond nearest neighbor or other perturbing parameters -- in two dimensions at intermediate to strong coupling and near optimal doping. The nature of such correlations has been a central question ever since the discovery of cuprate high-temperature superconductors. Despite unprecedented effort and tremendous progress in understanding the properties of this fundamental model, a definitive answer to whether the ground state is superconducting in the parameter regime most relevant to cuprates has proved exceedingly difficult to establish. In this work, we employ two complementary, state-of-the-art many-body computational methods, auxiliary-field quantum Monte Carlo (AFQMC) and density matrix renormalization group (DMRG) methods, deploying the most recent algorithmic advances in each. Systematic and detailed comparisons between the two methods are performed. The DMRG is extremely reliable on small width cylinders, where we use it to validate the AFQMC. The AFQMC is then used to study wide systems as well as fully periodic systems, to establish that we have reached the thermodynamic limit. The ground state is found to be non-superconducting in the moderate to strong coupling regime in the vicinity of optimal hole doping.

Accurate and computationally efficient modeling of systems of interacting electrons is an outstanding problem in theoretical and computational materials science. For materials where strong electronic interactions are primarily of a localized character and act within a subspace of localized quantum states on separate atomic sites (e.g., in transition metal and rare-earth compounds), their electronic behaviors are typically described by the Hubbard model and its extensions. In this work, we describe BoSS (Boson Slave Solver), a software implementation of the slave-boson method appropriate for describing a variety of extended Hubbard models, namely p-d models that include both the interacting atomic sites ("d" states) and non-interacting or ligand sites ("p" states). We provide a theoretical background, a description of the equations solved by BoSS, an overview of the algorithms used, the key input/output and control variables of the software program, and tutorial examples of its use featuring band renormalization in SrVO

Self-diffusion in dense rod suspensions are subject to strong geometric constraints because of steric interactions. This topological effect is essentially anisotropic when rods are nematically aligned with their neighbors, raising considerable challenges in understanding and analyzing their impacts on the bulk physical properties. Via a classical Doi-Onsager kinetic model with the Maier-Saupe potential, we characterize the long-time rotational Brownian diffusivity for dense suspensions of hard rods of finite aspect ratios, based on quadratic orientation autocorrelation functions. Furthermore, we show that the computed nonmonotonic scalings of the diffusivity as a function of volume fraction can be accurately predicted by an alternative tube model in the nematic phase.

The nonlinear optical response of an excitonic insulator coupled to lattice degrees of freedom is shown to depend in strong and characteristic ways on whether the insulating behavior originates primarily from electron-electron or electron-lattice interactions. Linear response optical signatures of the massive phase mode and the amplitude (Higgs) mode are identified. Upon nonlinear excitation resonant to the phase mode, a new in-gap mode at twice the phase mode frequency is induced, leading to a huge second harmonic response. Excitation of in-gap phonon modes leads to different and much smaller effects. A Landau-Ginzburg theory analysis explain these different behavior and reveals that a parametric resonance of the strongly excited phase mode is the origin of the photo-induced mode in the electron-dominant case. The difference in the nonlinear optical response serve as a measure of the dominant mechanism of the ordered phase.

The famous Arrhenius equation is well motivated to describe the temperature dependence of chemical reactions but has also been used for complicated biological processes. Here, we evaluate how well the simple Arrhenius equation predicts complex multistep biological processes, using frog and fruit fly embryogenesis as two canonical models. We find the Arrhenius equation provides a good approximation for the temperature dependence of embryogenesis, even though individual developmental stages scale differently with temperature. At low and high temperatures, however, we observed significant departures from idealized Arrhenius Law behavior. When we model multistep reactions of idealized chemical networks we are unable to generate comparable deviations from linearity. In contrast, we find the single enzyme GAPDH shows non-linearity in the Arrhenius plot similar to our observations of embryonic development. Thus, we find that complex embryonic development can be well approximated by the simple Arrhenius Law and propose that the observed departure from this law results primarily from non-idealized individual steps rather than the complexity of the system.

The remarkable null pulse coincident with the 2016 glitch in Vela rotation indicates a dynamical event involving the crust and the magnetosphere of the neutron star. We propose that a crustal quake associated with the glitch strongly disturbed the Vela magnetosphere and thus interrupted its radio emission. We present the first global numerical simulations of a neutron starquake. Our code resolves the elasto-dynamics of the entire crust and follows the evolution of Alfvén waves excited in the magnetosphere. We observe Rayleigh surface waves propagating away from the epicenter of the quake, around the circumference of the crust - an instance of the so-called whispering gallery modes. The Rayleigh waves set the initial spatial scale of the magnetospheric disturbance. Once launched, the Alfvén waves bounce in the closed magnetosphere, become de-phased, and generate strong electric currents, capable of igniting electric discharge. Most likely, the discharge floods the magnetosphere with electron-positron plasma, quenching the radio emission. We find that the observed ∼0.2 s disturbance is consistent with the damping time of the crustal waves if the crust is magnetically coupled to the superconducting core of the neutron star. The quake is expected to produce a weak X-ray burst of short duration.

Dorsal root ganglia (DRG) are promising sites for recording sensory activity. Current technologies for DRG recording are stiff and typically do not have sufficient site density for high-density neural data techniques. We demonstrate neural recordings in feline sacral DRG using a flexible polyimide microelectrode array with 30-40 μm site spacing. We delivered arrays into DRG with ultrananocrystalline diamond shuttles designed for high stiffness with small footprint. We recorded neural activity during sensory activation, including cutaneous brushing and bladder filling. We successfully delivered arrays in 5/6 experiments and recorded sensory activity in 4. Median signal amplitude was 55 μV and the maximum unique units recorded at one array position was 260, with 157 driven by sensory or electrical stimulation. We used specialized high-density neural signal analysis software to sort neural signals and, in one experiment, track 8 signals as the array was retracted 500 μm. This study is the first demonstration of ultrathin, flexible, high-density electronics delivered into DRG, with capabilities for recording and tracking sensory information that are a significant improvement over conventional DRG interfaces.Competing Interest StatementT.M.B. is a named inventors on a granted patent (US9622671B2; assigned to University of Pittsburgh) which is on the monitoring of physiological states via microelectrodes at DRG. The authors declare no other personal or institutional interest with regards to the authorship and/or publication of this manuscript.