Publications

We study the controlled manipulation of the Jahn-Teller metal state of fulleride compounds using nonequilibrium dynamical mean-field theory. This anomalous metallic state is a spontaneous orbital-selective Mott phase, which is characterized by one metallic and two insulating orbitals. Using protocols based on transiently reduced hopping amplitudes or periodic electric fields, we demonstrate the possibility to switch orbitals between Mott insulating and metallic on a subpicosecond time scale, and to rotate the order parameter between three equivalent states that can be distinguished by their anisotropic conductance. The Jahn-Teller metal phase of alkali-doped fullerides thus provides a platform for ultrafast persistent memory devices.

We present a proper motion measurement for the halo globular cluster Pyxis, using HST/ACS data as the first epoch, and GeMS/GSAOI Adaptive Optics data as the second, separated by a baseline of about 5 years. This is both the first measurement of the proper motion of Pyxis and the first calibration and use of Multi-Conjugate Adaptive Optics data to measure an absolute proper motion for a faint, distant halo object. Consequently, we present our analysis of the Adaptive Optics data in detail. We obtain a proper motion of mu_alpha cos(delta)=1.09+/-0.31 mas/yr and mu_delta=0.68+/-0.29 mas/yr. From the proper motion and the line-of-sight velocity we find the orbit of Pyxis is rather eccentric with its apocenter at more than 100 kpc and its pericenter at about 30 kpc. We also investigate two literature-proposed associations for Pyxis with the recently discovered ATLAS stream and the Magellanic system. Combining our measurements with dynamical modeling and cosmological numerical simulations we find it unlikely Pyxis is associated with either system. We examine other Milky Way satellites for possible association using the orbit, eccentricity, metallicity, and age as constraints and find no likely matches in satellites down to the mass of Leo II. We propose that Pyxis probably originated in an unknown galaxy, which today is fully disrupted. Assuming that Pyxis is bound and not on a first approach, we derive a 68% lower limit on the mass of the Milky Way of 0.95*10^12 M_sun.

The design of systems or models that work robustly under uncertainty and environmental fluctuations is a key challenge in both engineering and science. This is formalized in the design-centering problem, which is defined as finding a design that fulfills given specifications and has a high probability of still doing so if the system parameters or the specifications fluctuate randomly. Design centering is often accompanied by the problem of quantifying the robustness of a system. Here we present a novel adaptive statistical method to simultaneously address both problems. Our method, L p-Adaptation, is inspired by the evolution of robustness in biological systems and by randomized schemes for convex volume computation. It is able to address both problems in the general, non-convex case and at low computational cost. We describe the concept and the algorithm, test it on known benchmarks, and demonstrate its real-world applicability in electronic and biological systems. In all cases, the present method outperforms the previous state of the art. This enables re-formulating optimization problems in engineering and biology as design centering problems, taking global system robustness into account.

The swimming direction of biological or artificial microscale swimmers tends to be randomised over long time-scales by thermal fluctuations. Bacteria use various strategies to bias swimming behaviour and achieve directed motion against a flow, maintain alignment with gravity or travel up a chemical gradient. Herein, we explore a purely geometric means of biasing the motion of artificial nanorod swimmers. These artificial swimmers are bimetallic rods, powered by a chemical fuel, which swim on a substrate printed with teardrop-shaped posts. The artificial swimmers are hydrodynamically attracted to the posts, swimming alongside the post perimeter for long times before leaving. The rods experience a higher rate of departure from the higher curvature end of the teardrop shape, thereby introducing a bias into their motion. This bias increases with swimming speed and can be translated into a macroscopic directional motion over long times by using arrays of teardrop-shaped posts aligned along a single direction. This method provides a protocol for concentrating swimmers, sorting swimmers according to different speeds, and could enable artificial swimmers to transport cargo to desired locations.

Analysing computations in neural circuits often uses simplified models because the actual neuronal implementation is not known. For example, a problem in vision, how the eye detects image motion, has long been analysed using Hassenstein-Reichardt (HR) detector or Barlow-Levick (BL) models. These both simulate motion detection well, but the exact neuronal circuits undertaking these tasks remain elusive. We reconstructed a comprehensive connectome of the circuits of Drosophila‘s motion-sensing T4 cells using a novel EM technique. We uncover complex T4 inputs and reveal that putative excitatory inputs cluster at T4’s dendrite shafts, while inhibitory inputs localize to the bases. Consistent with our previous study, we reveal that Mi1 and Tm3 cells provide most synaptic contacts onto T4. We are, however, unable to reproduce the spatial offset between these cells reported previously. Our comprehensive connectome reveals complex circuits that include candidate anatomical substrates for both HR and BL types of motion detectors.

Analysing computations in neural circuits often uses simplified models because the actual neuronal implementation is not known. For example, a problem in vision, how the eye detects image motion, has long been analysed using Hassenstein-Reichardt (HR) detector or Barlow-Levick (BL) models. These both simulate motion detection well, but the exact neuronal circuits undertaking these tasks remain elusive. We reconstructed a comprehensive connectome of the circuits of Drosophila‘s motion-sensing T4 cells using a novel EM technique. We uncover complex T4 inputs and reveal that putative excitatory inputs cluster at T4’s dendrite shafts, while inhibitory inputs localize to the bases. Consistent with our previous study, we reveal that Mi1 and Tm3 cells provide most synaptic contacts onto T4. We are, however, unable to reproduce the spatial offset between these cells reported previously. Our comprehensive connectome reveals complex circuits that include candidate anatomical substrates for both HR and BL types of motion detectors.

Over the past decade, the Rosetta biomolecular modeling suite has informed diverse biological questions and engineering challenges ranging from interpretation of low-resolution structural data to design of nanomaterials, protein therapeutics, and vaccines. Central to Rosetta’s success is the energy function: a model parametrized from small-molecule and X-ray crystal structure data used to approximate the energy associated with each biomolecule conformation. This paper describes the mathematical models and physical concepts that underlie the latest Rosetta energy function, called the Rosetta Energy Function 2015 (REF15). Applying these concepts, we explain how to use Rosetta energies to identify and analyze the features of biomolecular models. Finally, we discuss the latest advances in the energy function that extend its capabilities from soluble proteins to also include membrane proteins, peptides containing noncanonical amino acids, small molecules, carbohydrates, nucleic acids, and other macromolecules.

We study the relaxation dynamics of photocarriers in the paramagnetic Mott insulating phase of the half-filled two-band Hubbard model. Using nonequilibrium dynamical mean-field theory, we excite charge carriers across the Mott gap by a short hopping modulation, and simulate the evolution of the photodoped population within the Hubbard bands. We observe an ultrafast charge-carrier relaxation driven by the emission of local spin excitations with an inverse relaxation time proportional to the Hund's coupling. The photodoping generates additional side-bands in the spectral function, and for strong Hund's coupling, the photodoped population also splits into several resonances. The dynamics of the local many-body states reveals two effects, thermal blocking and kinetic freezing, which manifest themselves when the Hund's coupling becomes of the order of the temperature or the bandwidth, respectively. These effects, which are absent in the single-band Hubbard model, should be relevant for the interpretation of experiments on correlated materials with multiple active orbitals. In particular, the features revealed in the nonequilibrium energy distribution of the photocarriers are experimentally accessible, and provide information on the role of the Hund's coupling in these materials.

Standard approaches to Bayesian parameter inference in large-scale structure assume a Gaussian functional form (chi-squared form) for the likelihood. This assumption, in detail, cannot be correct. Likelihood free inferences such as approximate Bayesian computation (ABC) relax these restrictions and make inference possible without making any assumptions on the likelihood. Instead ABC relies on a forward generative model of the data and a metric for measuring the distance between the model and data. In this work, we demonstrate that ABC is feasible for LSS parameter inference by using it to constrain parameters of the halo occupation distribution (HOD) model for populating dark matter haloes with galaxies. Using specific implementation of ABC supplemented with population Monte Carlo importance sampling, a generative forward model using HOD and a distance metric based on galaxy number density, two-point correlation function and galaxy group multiplicity function, we constrain the HOD parameters of mock observation generated from selected ‘true' HOD parameters. The parameter constraints we obtain from ABC are consistent with the ‘true' HOD parameters, demonstrating that ABC can be reliably used for parameter inference in LSS. Furthermore, we compare our ABC constraints to constraints we obtain using a pseudo-likelihood function of Gaussian form with MCMC and find consistent HOD parameter constraints. Ultimately, our results suggest that ABC can and should be applied in parameter inference for LSS analyses.

The instant invention provides methods and related compositions for identifying polypeptides with improved stability and/or enzymatic activity in comparison to native forms, wherein the identified polypeptides comprise one or more non-natural amino acids. In certain embodiments, the present invention relates to novel phosphotriesterase enzymes comprising one or more non-natural amino acids. In a particular embodiment, the instant invention provides novel phosphotriesterase enzymes with greater stability and/or enhanced activity in comparison to native forms of the enzyme. The present invention also relates to compositions comprising novel phophotriesterase enzymes, such as prophylactics, decontaminants, animal feedstocks, and assay kits.