Bifurcation into functional niches in adaptation

One of the central questions in evolutionary biology concerns the dynamics of adaptation and diversification. This issue can be addressed experimentally if replicate populations adapting to identical environments Call be investigated in detail. We have studied 501 such replicas Using digital organisms adapting to at least two fundamentally different functional niches (survival strategies) present in the same environment: one in which fast replication is the way to live, and another where exploitation of the environment's complexity leads to complex organisms with longer life spans and smaller replication rates. While these two modes of survival are closely analogous to those expected to emerge in so-called r and K selection scenarios respectively, the bifurcation of evolutionary histories according to these functional niches occurs in identical environments, under identical selective pressures. We find that the branching occurs early, and leads to drastic phenotypic differences (in fitness, sequence length, and gestation time) that are permanent and irreversible. This study confirms an earlier experimental effort using microorganisms, in that diversification can be understood at least in part in terms of bifurcations on saddle points leading to peak shifts, as in the picture drawn by Sewall Wright

Correct ordering in the Zipf-Poisson ensemble

We consider a Zipf--Poisson ensemble in which X_i\sim\poi(Ni^{-\alpha}) for $\alpha>1$ and $N>0$ and integers $i\ge 1$. As $N\to\infty$ the first $n'(N)$ random variables have their proper order $X_1>X_2>...>X_{n'}$ relative to each other, with probability tending to 1 for $n'$ up to $(AN/\log(N))^{1/(\alpha+2)}$ for an explicit constant $A(\alpha)\ge 3/4$. The rate $N^{1/(\alpha+2)}$ cannot be achieved. The ordering of the first $n'(N)$ entities does not preclude $X_m>X_{n'}$ for some interloping $m>n'$. The first $n"$ random variables are correctly ordered exclusive of any interlopers, with probability tending to 1 if $n"\le (BN/\log(N))^{1/(\alpha+2)}$ for $B<A$. For a Zipf--Poisson model of the British National Corpus, which has a total word count of $100{,}000{,}000$, our result estimates that the 72 words with the highest counts are properly ordered

Parametric inference in the large data limit using maximally informative models

Motivated by data-rich experiments in transcriptional regulation and sensory neuroscience, we consider the following general problem in statistical inference. When exposed to a high-dimensional signal S, a system of interest computes a representation R of that signal which is then observed through a noisy measurement M. From a large number of signals and measurements, we wish to infer the "filter" that maps S to R. However, the standard method for solving such problems, likelihood-based inference, requires perfect a priori knowledge of the "noise function" mapping R to M. In practice such noise functions are usually known only approximately, if at all, and using an incorrect noise function will typically bias the inferred filter. Here we show that, in the large data limit, this need for a pre-characterized noise function can be circumvented by searching for filters that instead maximize the mutual information I[M;R] between observed measurements and predicted representations. Moreover, if the correct filter lies within the space of filters being explored, maximizing mutual information becomes equivalent to simultaneously maximizing every dependence measure that satisfies the Data Processing Inequality. It is important to note that maximizing mutual information will typically leave a small number of directions in parameter space unconstrained. We term these directions "diffeomorphic modes" and present an equation that allows these modes to be derived systematically. The presence of diffeomorphic modes reflects a fundamental and nontrivial substructure within parameter space, one that is obscured by standard likelihood-based inference.Comment: To appear in Neural Computatio

Equitability, mutual information, and the maximal information coefficient

Reshef et al. recently proposed a new statistical measure, the "maximal information coefficient" (MIC), for quantifying arbitrary dependencies between pairs of stochastic quantities. MIC is based on mutual information, a fundamental quantity in information theory that is widely understood to serve this need. MIC, however, is not an estimate of mutual information. Indeed, it was claimed that MIC possesses a desirable mathematical property called "equitability" that mutual information lacks. This was not proven; instead it was argued solely through the analysis of simulated data. Here we show that this claim, in fact, is incorrect. First we offer mathematical proof that no (non-trivial) dependence measure satisfies the definition of equitability proposed by Reshef et al.. We then propose a self-consistent and more general definition of equitability that follows naturally from the Data Processing Inequality. Mutual information satisfies this new definition of equitability while MIC does not. Finally, we show that the simulation evidence offered by Reshef et al. was artifactual. We conclude that estimating mutual information is not only practical for many real-world applications, but also provides a natural solution to the problem of quantifying associations in large data sets

Flight Mechanics Modeling and Post-Flight Analysis of ADEPT SR-1

Sounding Rocket One (SR-1), the first flight test of the Adaptable Deployable Entry and Placement Technology (ADEPT), was performed on Sept. 12, 2018. ADEPT is a deployable aeroshell that can be stowed during launch and then opened after launch to increase the drag area of the spacecraft when entering into a planetary atmosphere. The main objectives of the SR-1 flight test were to demonstrate that the ADEPT vehicle can be opened exo-atmospherically and to characterize the stability of the vehicle during atmospheric flight. The SR-1 test vehicle was a 0.7 m diameter 70 degree half-angle, faceted, sphere-cone, which was the primary payload on an UP Aerospace Spaceloft (SL) launch vehicle from the White Sands Missile Range (WSMR). ADEPT successfully separated from the spent booster in its stowed configuration, opened above 100 km altitude, and then landed in the deployed configuration within WSMR. The flight mechanics of the vehicle was modeled pre-flight for performance and range safety predictions. This paper describes the pre-flight ADEPT trajectory simulation and how the flight data compared with the predictions from the simulations

Berryogenesis: self-induced Berry flux and spontaneous non-equilibrium magnetism

Spontaneous symmetry breaking is central to the description of interacting phases of matter. Here we reveal a new mechanism through which a driven interacting system subject to a time-reversal symmetric driving field can spontaneously magnetize. We show that the strong internal ac fields of a metal driven close to its plasmon resonance may enable Berryogenesis: the spontaneous generation of a self-induced Bloch band Berry flux. The self-induced Berry flux supports and is sustained by a circulating plasmonic motion, which may arise even for a linearly polarized driving field. This non-equilibrium phase transition occurs above a critical driving amplitude, and may be of either continuous or discontinuous type. Berryogenesis relies on feedback due to interband coherences induced by internal fields, and may readily occur in a wide variety of multiband systems. We anticipate that graphene devices, in particular, provide a natural platform to achieve Berryogenesis and plasmon-mediated spontaneous non-equilibrium magnetization in present-day devices

Energy-driven Drag at Charge Neutrality in Graphene

Coulomb coupling between proximal layers in graphene heterostructures results in efficient energy transfer between the layers. We predict that, in the presence of correlated density inhomogeneities in the layers, vertical energy transfer has a strong impact on lateral charge transport. In particular, for Coulomb drag it dominates over the conventional momentum drag near zero doping. The dependence on doping and temperature, which is different for the two drag mechanisms, can be used to separate these mechanisms in experiment. We predict distinct features such as a peak at zero doping and a multiple sign reversal, which provide diagnostics for this new drag mechanism.Comment: 6 pgs, 3 fg