Showing 1–37 of 37 results for author: Daniels, C

Continuous Attractor Networks for Laplace Neural Manifolds

Authors: Bryan C. Daniels, Marc W. Howard

Abstract: Many cognitive models, including those for predicting the time of future events, can be mapped onto a particular form of neural representation in which activity across a population of neurons is restricted to manifolds that specify the Laplace transform of functions of continuous variables. These populations coding Laplace transform are associated with another population that inverts the transform… ▽ More Many cognitive models, including those for predicting the time of future events, can be mapped onto a particular form of neural representation in which activity across a population of neurons is restricted to manifolds that specify the Laplace transform of functions of continuous variables. These populations coding Laplace transform are associated with another population that inverts the transform, approximating the original function. This paper presents a neural circuit that uses continuous attractor dynamics to represent the Laplace transform of a delta function evolving in time. One population places an edge at any location along a 1-D array of neurons; another population places a bump at a location corresponding to the edge. Together these two populations can estimate a Laplace transform of delta functions in time along with an approximate inverse transform. Building the circuit so the edge moves at an appropriate speed enables the network to represent events as a function of log time. Choosing the connections appropriately within the edge network make the network states map onto Laplace transform with exponential change as a function of time. In this paper we model a learned temporal association in which one stimulus predicts another at some fixed delay $T$. Shortly after $t=0$ the first stimulus recedes into the past. The Laplace Neural Manifold representing the past maintains the Laplace transform $\exp(-st)$. Another Laplace Neural Manifold represents the predicted future. At $t=0$, the second stimulus is represented a time $T$ in the future. At each moment between 0 and $T$, firing over the Laplace transform predicting the future changes as $\exp[-s(T-t)]$. Despite exponential growth in firing, the circuit is robust to noise, making it a practical means to implement Laplace Neural Manifolds in populations of neurons for a variety of cognitive models. △ Less

Submitted 6 June, 2024; originally announced June 2024.

Comments: 31 pages, 9 figures

Difficult control is related to instability in biologically inspired Boolean networks

Authors: Bryan C. Daniels, Enrico Borriello

Abstract: Previous work in Boolean dynamical networks has suggested that the number of components that must be controlled to select an existing attractor is typically set by the number of attractors admitted by the dynamics, with no dependence on the size of the network. Here we study the rare cases of networks that defy this expectation, with attractors that require controlling most nodes. We find empirica… ▽ More Previous work in Boolean dynamical networks has suggested that the number of components that must be controlled to select an existing attractor is typically set by the number of attractors admitted by the dynamics, with no dependence on the size of the network. Here we study the rare cases of networks that defy this expectation, with attractors that require controlling most nodes. We find empirically that unstable fixed points are the primary recurring characteristic of networks that prove more difficult to control. We describe an efficient way to identify unstable fixed points and show that, in both existing biological models and ensembles of random dynamics, we can better explain the variance of control kernel sizes by incorporating the prevalence of unstable fixed points. In the end, the fact that these exceptions are associated with dynamics that are unstable to small perturbations hints that they are likely an artifact of using deterministic models. These exceptions are likely to be biologically irrelevant, supporting the generality of easy controllability in biological networks. △ Less

Submitted 28 February, 2024; originally announced February 2024.

Comments: 26 pages, 6 figures

Polynomial Chaos Surrogate Construction for Random Fields with Parametric Uncertainty

Authors: Joy N. Mueller, Khachik Sargsyan, Craig J. Daniels, Habib N. Najm

Abstract: Engineering and applied science rely on computational experiments to rigorously study physical systems. The mathematical models used to probe these systems are highly complex, and sampling-intensive studies often require prohibitively many simulations for acceptable accuracy. Surrogate models provide a means of circumventing the high computational expense of sampling such complex models. In partic… ▽ More Engineering and applied science rely on computational experiments to rigorously study physical systems. The mathematical models used to probe these systems are highly complex, and sampling-intensive studies often require prohibitively many simulations for acceptable accuracy. Surrogate models provide a means of circumventing the high computational expense of sampling such complex models. In particular, polynomial chaos expansions (PCEs) have been successfully used for uncertainty quantification studies of deterministic models where the dominant source of uncertainty is parametric. We discuss an extension to conventional PCE surrogate modeling to enable surrogate construction for stochastic computational models that have intrinsic noise in addition to parametric uncertainty. We develop a PCE surrogate on a joint space of intrinsic and parametric uncertainty, enabled by Rosenblatt transformations, and then extend the construction to random field data via the Karhunen-Loeve expansion. We then take advantage of closed-form solutions for computing PCE Sobol indices to perform a global sensitivity analysis of the model which quantifies the intrinsic noise contribution to the overall model output variance. Additionally, the resulting joint PCE is generative in the sense that it allows generating random realizations at any input parameter setting that are statistically approximately equivalent to realizations from the underlying stochastic model. The method is demonstrated on a chemical catalysis example model. △ Less

Submitted 17 June, 2024; v1 submitted 1 November, 2023; originally announced November 2023.

MSC Class: 60G99; 65C20; 33C45; 62G07; 62J02

Phase Transitions and Criticality in the Collective Behavior of Animals -- Self-organization and biological function

Authors: Pawel Romanczuk, Bryan C. Daniels

Abstract: Collective behaviors exhibited by animal groups, such as fish schools, bird flocks, or insect swarms are fascinating examples of self-organization in biology. Concepts and methods from statistical physics have been used to argue theoretically about the potential consequences of collective effects in such living systems. In particular, it has been proposed that such collective systems should operat… ▽ More Collective behaviors exhibited by animal groups, such as fish schools, bird flocks, or insect swarms are fascinating examples of self-organization in biology. Concepts and methods from statistical physics have been used to argue theoretically about the potential consequences of collective effects in such living systems. In particular, it has been proposed that such collective systems should operate close to a phase transition, specifically a (pseudo-)critical point, in order to optimize their capability for collective computation. In this chapter, we will first review relevant phase transitions exhibited by animal collectives, pointing out the difficulties of applying concepts from statistical physics to biological systems. Then we will discuss the current state of research on the "criticality hypothesis", including methods for how to measure distance from criticality and specific functional consequences for animal groups operating near a phase transition. We will highlight the emerging view that de-emphasizes the optimality of being exactly at a critical point and instead explores the potential benefits of living systems being able to tune to an optimal distance from criticality. We will close by laying out future challenges for studying collective behavior at the interface of physics and biology. △ Less

Submitted 7 November, 2022; originally announced November 2022.

Comments: to appear in "Order, disorder, and criticality", vol. VII, World Scientific Publishing

Growth optimization and device integration of narrow-bandgap graphene nanoribbons

Authors: Gabriela Borin Barin, Qiang Sun, Marco Di Giovannantonio, Cheng-Zhuo Du, Xiao-Ye Wang, Juan Pablo Llinas, Zafer Mutlu, Yuxuan Lin, Jan Wilhelm, Jan Overbeck, Colin Daniels, Michael Lamparski, Hafeesudeen Sahabudeen, Mickael L. Perrin, José I. Urgel, Shantanu Mishra, Amogh Kinikar, Roland Widmer, Samuel Stolz, Max Bommert, Carlo Pignedoli, Xinliang Feng, Michel Calame, Klaus Müllen, Akimitsu Narita , et al. (4 additional authors not shown)

Abstract: The electronic, optical and magnetic properties of graphene nanoribbons (GNRs) can be engineered by controlling their edge structure and width with atomic precision through bottom-up fabrication based on molecular precursors. This approach offers a unique platform for all-carbon electronic devices but requires careful optimization of the growth conditions to match structural requirements for succe… ▽ More The electronic, optical and magnetic properties of graphene nanoribbons (GNRs) can be engineered by controlling their edge structure and width with atomic precision through bottom-up fabrication based on molecular precursors. This approach offers a unique platform for all-carbon electronic devices but requires careful optimization of the growth conditions to match structural requirements for successful device integration, with GNR length being the most critical parameter. In this work, we study the growth, characterization, and device integration of 5-atom wide armchair GNRs (5-AGNRs), which are expected to have an optimal band gap as active material in switching devices. 5-AGNRs are obtained via on-surface synthesis under ultra-high vacuum conditions from Br- and I-substituted precursors. We show that the use of I-substituted precursors and the optimization of the initial precursor coverage quintupled the average 5-AGNR length. This significant length increase allowed us to integrate 5-AGNRs into devices and to realize the first field-effect transistor based on narrow bandgap AGNRs that shows switching behavior at room temperature. Our study highlights that optimized growth protocols can successfully bridge between the sub-nanometer scale, where atomic precision is needed to control the electronic properties, and the scale of tens of nanometers relevant for successful device integration of GNRs. △ Less

Submitted 2 February, 2022; originally announced February 2022.

Subcritical escape waves in schooling fish

Authors: Winnie Poel, Bryan C. Daniels, Matthew M. G. Sosna, Colin R. Twomey, Simon P. Leblanc, Iain D. Couzin, Pawel Romanczuk

Abstract: Living systems such as neuronal networks and animal groups process information about their environment via the dynamics of interacting units. These can transition between distinct macroscopic behaviors. Near such a transition (or critical point) collective computation is generally thought to be optimized, due to the associated maximal sensitivity to perturbations and fast dissemination of informat… ▽ More Living systems such as neuronal networks and animal groups process information about their environment via the dynamics of interacting units. These can transition between distinct macroscopic behaviors. Near such a transition (or critical point) collective computation is generally thought to be optimized, due to the associated maximal sensitivity to perturbations and fast dissemination of information. For biological systems, however, optimality depends on environmental context, making the flexible, context-dependent adoption of different distances to a critical point potentially more beneficial than its unique properties. Here, studying escape waves in schooling fish at two levels of perceived environmental risk, we investigate a) if and how distance to criticality is regulated in response to environmental changes and b) how the individual level benefits derived from special properties of the critical point compare to those achieved via regulation of the group's distance to it. We find that the observed fish schools are subcritical (not maximally responsive and sensitive to environmental cues), but decrease their distance to criticality with increased perceived risk. Considering an individual's hypothetical costs of two detection error types, we find that optimal distance to criticality depends on the riskiness and noisiness of the environment, which may explain the observed behavior. Our results highlight the benefit of evaluating biological consequences of different distances to criticality for individuals within animal collectives. This provides insights into the adaptive function of a collective system and motivates future questions about the evolutionary forces that brought the system to make this particular trade-off. △ Less

Submitted 12 August, 2021; originally announced August 2021.

Discovering sparse control strategies in C. elegans

Authors: Edward D. Lee, Xiaowen Chen, Bryan C. Daniels

Abstract: Biological circuits such as neural or gene regulation networks use internal states to map sensory input to an adaptive repertoire of behavior. Characterizing this mapping is a major challenge for systems biology, and though experiments that probe internal states are developing rapidly, organismal complexity presents a fundamental obstacle given the many possible ways internal states could map to b… ▽ More Biological circuits such as neural or gene regulation networks use internal states to map sensory input to an adaptive repertoire of behavior. Characterizing this mapping is a major challenge for systems biology, and though experiments that probe internal states are developing rapidly, organismal complexity presents a fundamental obstacle given the many possible ways internal states could map to behavior. Using C. elegans as an example, we propose a protocol for systematic perturbation of neural states that limits experimental complexity but still characterizes collective aspects of the neural-behavioral map. We consider experimentally motivated small perturbations -- ones that are most likely to preserve natural dynamics and are closer to internal control mechanisms -- to neural states and their impact on collective neural behavior. Then, we connect such perturbations to the local information geometry of collective statistics, which can be fully characterized using pairwise perturbations. Applying the protocol to a minimal model of C. elegans neural activity, we find that collective neural statistics are most sensitive to a few principal perturbative modes. Dominant eigenvalues decay initially as a power law, unveiling a hierarchy that arises from variation in individual neural activity and pairwise interactions. Highest-ranking modes tend to be dominated by a few, "pivotal" neurons that account for most of the system's sensitivity, suggesting a sparse mechanism for control of collective behavior. △ Less

Submitted 2 August, 2021; originally announced August 2021.

A Study of Policy Gradient on a Class of Exactly Solvable Models

Authors: Gavin McCracken, Colin Daniels, Rosie Zhao, Anna Brandenberger, Prakash Panangaden, Doina Precup

Abstract: Policy gradient methods are extensively used in reinforcement learning as a way to optimize expected return. In this paper, we explore the evolution of the policy parameters, for a special class of exactly solvable POMDPs, as a continuous-state Markov chain, whose transition probabilities are determined by the gradient of the distribution of the policy's value. Our approach relies heavily on rando… ▽ More Policy gradient methods are extensively used in reinforcement learning as a way to optimize expected return. In this paper, we explore the evolution of the policy parameters, for a special class of exactly solvable POMDPs, as a continuous-state Markov chain, whose transition probabilities are determined by the gradient of the distribution of the policy's value. Our approach relies heavily on random walk theory, specifically on affine Weyl groups. We construct a class of novel partially observable environments with controllable exploration difficulty, in which the value distribution, and hence the policy parameter evolution, can be derived analytically. Using these environments, we analyze the probabilistic convergence of policy gradient to different local maxima of the value function. To our knowledge, this is the first approach developed to analytically compute the landscape of policy gradient in POMDPs for a class of such environments, leading to interesting insights into the difficulty of this problem. △ Less

Submitted 3 November, 2020; originally announced November 2020.

The basis of easy controllability in Boolean networks

Authors: Enrico Borriello, Bryan C. Daniels

Abstract: Effective control of biological systems can often be achieved through the control of a surprisingly small number of distinct variables. We bring clarity to such results using the formalism of Boolean dynamical networks, analyzing the effectiveness of external control in selecting a desired final state when that state is among the original attractors of the dynamics. Analyzing 49 existing biologica… ▽ More Effective control of biological systems can often be achieved through the control of a surprisingly small number of distinct variables. We bring clarity to such results using the formalism of Boolean dynamical networks, analyzing the effectiveness of external control in selecting a desired final state when that state is among the original attractors of the dynamics. Analyzing 49 existing biological network models, we find strong numerical evidence that the average number of nodes that must be forced scales logarithmically with the number of original attractors. This suggests that biological networks may be typically easy to control even when the number of interacting components is large. We provide a theoretical explanation of the scaling by separating controlling nodes into three types: those that act as inputs, those that distinguish among attractors, and any remaining nodes. We further identify characteristics of dynamics that can invalidate this scaling, and speculate about how this relates more broadly to non-biological systems. △ Less

Submitted 9 September, 2021; v1 submitted 22 October, 2020; originally announced October 2020.

Comments: 44 pages, 17 figures, 2 table

Journal ref: Nat Commun 12, 5227 (2021)

Quantifying dynamical high-order interdependencies from the O-information: an application to neural spiking dynamics

Authors: Sebastiano Stramaglia, Tomas Scagliarini, Bryan C. Daniels, Daniele Marinazzo

Abstract: We address the problem of efficiently and informatively quantifying how multiplets of variables carry information about the future of the dynamical system they belong to. In particular we want to identify groups of variables carrying redundant or synergistic information, and track how the size and the composition of these multiplets changes as the collective behavior of the system evolves. In orde… ▽ More We address the problem of efficiently and informatively quantifying how multiplets of variables carry information about the future of the dynamical system they belong to. In particular we want to identify groups of variables carrying redundant or synergistic information, and track how the size and the composition of these multiplets changes as the collective behavior of the system evolves. In order to afford a parsimonious expansion of shared information, and at the same time control for lagged interactions and common effect, we develop a dynamical, conditioned version of the O-information, a framework recently proposed to quantify high-order interdependencies via multivariate extension of the mutual information. We thus obtain an expansion of the transfer entropy in which synergistic and redundant effects are separated. We apply this framework to a dataset of spiking neurons from a monkey performing a perceptual discrimination task. The method identifies synergistic multiplets that include neurons previously categorized as containing little relevant information individually. △ Less

Submitted 31 July, 2020; originally announced July 2020.

Machine-learning models for Raman spectra analysis of twisted bilayer graphene

Authors: Natalya Sheremetyeva, Michael Lamparski, Colin Daniels, Benoit Van Troeye, Vincent Meunier

Abstract: The vibrational properties of twisted bilayer graphene (tBLG) show complex features, due to the intricate energy landscape of its low-symmetry configurations. A machine learning-based approach is developed to provide a continuous model between the twist angle and the simulated Raman spectra of tBLGs. Extracting the structural information of the twist angle from Raman spectra corresponds to solving… ▽ More The vibrational properties of twisted bilayer graphene (tBLG) show complex features, due to the intricate energy landscape of its low-symmetry configurations. A machine learning-based approach is developed to provide a continuous model between the twist angle and the simulated Raman spectra of tBLGs. Extracting the structural information of the twist angle from Raman spectra corresponds to solving a complicated inverse problem. Once trained, the machine learning regressors (MLRs) quickly provide predictions without human bias and with an average of 98% of the data variance being explained by the model. The significant spectral features learned by MLRs are analyzed revealing the intensity profile near the calculated G-band to be the most important feature. The trained models are tested on noise-containing test data demonstrating their robustness. The transferability of the present models to experimental Raman spectra is discussed in the context of validation of the level of theory used for the construction of the analyzed database. This work serves as a proof of concept that machine-learning analysis is a potentially powerful tool for interpretation of Raman spectra of tBLG and other 2D materials. △ Less

Submitted 7 July, 2020; originally announced July 2020.

Comments: 11 pages, 5 figures

A scaling theory of armed conflict avalanches

Authors: Edward D. Lee, Bryan C. Daniels, Christopher R. Myers, David C. Krakauer, Jessica C. Flack

Abstract: Armed conflict data display scaling and universal dynamics in both social and physical properties like fatalities and geographic extent. We propose a randomly branching, armed-conflict model that relates multiple properties to one another in a way consistent with data. The model incorporates a fractal lattice on which conflict spreads, uniform dynamics driving conflict growth, and regional virulen… ▽ More Armed conflict data display scaling and universal dynamics in both social and physical properties like fatalities and geographic extent. We propose a randomly branching, armed-conflict model that relates multiple properties to one another in a way consistent with data. The model incorporates a fractal lattice on which conflict spreads, uniform dynamics driving conflict growth, and regional virulence that modulates local conflict intensity. The quantitative constraints on scaling and universal dynamics we use to develop our minimal model serve more generally as a set of constraints for other models for armed conflict dynamics. We show how this approach akin to thermodynamics imparts mechanistic intuition and unifies multiple conflict properties, giving insight into causation, prediction, and intervention timing. △ Less

Submitted 29 April, 2020; originally announced April 2020.

Journal ref: Phys. Rev. E 102, 042312 (2020)

Towards Digital Engineering -- The Advent of Digital Systems Engineering

Authors: Jingwei Huang, Adrian Gheorghe, Holly Handley, Pilar Pazos, Ariel Pinto, Samuel Kovacic, Andy Collins, Charles Keating, Andres Sousa-Poza, Ghaith Rabadi, Resit Unal, Teddy Cotter, Rafael Landaeta, Charles Daniels

Abstract: Digital Engineering, the digital transformation of engineering to leverage digital technologies, is coming globally. This paper explores digital systems engineering, which aims at developing theory, methods, models, and tools to support the emerging digital engineering. A critical task is to digitalize engineering artifacts, thus enabling information sharing across platform, across life cycle, and… ▽ More Digital Engineering, the digital transformation of engineering to leverage digital technologies, is coming globally. This paper explores digital systems engineering, which aims at developing theory, methods, models, and tools to support the emerging digital engineering. A critical task is to digitalize engineering artifacts, thus enabling information sharing across platform, across life cycle, and across domains. We identify significant challenges and enabling digital technologies; analyze the transition from traditional engineering to digital engineering; define core concepts, including "digitalization", "unique identification", "digitalized artifacts", "digital augmentation", and others; present a big picture of digital systems engineering in four levels: vision, strategy, action, and foundation; briefly discuss each of main areas of research issues. Digitalization enables fast infusing and leveraging novel digital technologies; unique identification enables information traceability and accountability in engineering lifecycle; provenance enables tracing dependency relations among engineering artifacts; supporting model reproducibility and replicability; helping with trustworthiness evaluation of digital engineering artifacts. △ Less

Submitted 30 August, 2020; v1 submitted 20 February, 2020; originally announced February 2020.

Comments: 23 pages, 6 figures International Journal of System of Systems Engineering, (in press)

Journal ref: International Journal of System of Systems Engineering, 2020

A Universal Length-Dependent Vibrational Mode in Graphene Nanoribbons

Authors: Jan Overbeck, Gabriela Borin Barin, Colin Daniels, Mickael L. Perrin, Oliver Braun, Qiang Sun, Rimah Darawish, Marta De Luca, Xiao-Ye Wang, Tim Dumslaff, Akimitsu Narita, Klaus Müllen, Pascal Ruffieux, Vincent Meunier, Roman Fasel, Michel Calame

Abstract: Graphene nanoribbons (GNRs) have attracted considerable interest as their atomically tunable structure makes them promising candidates for future electronic devices. However, obtaining detailed information about the length of GNRs has been challenging and typically relies on low-temperature scanning tunneling microscopy. Such methods are ill-suited for practical device application and characteriza… ▽ More Graphene nanoribbons (GNRs) have attracted considerable interest as their atomically tunable structure makes them promising candidates for future electronic devices. However, obtaining detailed information about the length of GNRs has been challenging and typically relies on low-temperature scanning tunneling microscopy. Such methods are ill-suited for practical device application and characterization. In contrast, Raman spectroscopy is a sensitive method for the characterization of GNRs, in particular for investigating their width and structure. Here, we report on a length-dependent, Raman active low-energy vibrational mode that is present in atomically precise, bottom-up synthesized armchair graphene nanoribbons (AGNRs). Our Raman study demonstrates that this mode is present in all families of AGNRs and provides information on their length. Our spectroscopic findings are corroborated by scanning tunneling microscopy images and supported by first-principles calculations that allow us to attribute this mode to a longitudinal acoustic phonon. Finally, we show that this mode is a sensitive probe for the overall structural integrity of the ribbons and their interaction with technologically relevant substrates. △ Less

Submitted 14 December, 2019; originally announced December 2019.

Comments: Unedited author's version

Journal ref: ACS Nano 2019 13 (11), 13083-13091

In-Plane Breathing and Shear Modes in Low-Dimensional Nanostructures

Authors: Dan Liu, Colin Daniels, Vincent Meunier, Arthur G. Every, David Tomanek

Abstract: We use continuum elasticity theory to revise scaling laws for radial breathing-like and shear-like vibration modes in quasi-2D nanostructures including finite-width nanoribbons and finite-size thin circular discs. Such modes can be observed spectroscopically in corresponding nanostructures of graphene and phosphorene and can be determined numerically by atomistic <i> ab initio</i> density function… ▽ More We use continuum elasticity theory to revise scaling laws for radial breathing-like and shear-like vibration modes in quasi-2D nanostructures including finite-width nanoribbons and finite-size thin circular discs. Such modes can be observed spectroscopically in corresponding nanostructures of graphene and phosphorene and can be determined numerically by atomistic <i> ab initio</i> density functional theory and classical force-field calculation. The revised scaling laws differ from previously used expressions, some of which display an unphysical asymptotic behavior. Apart from model assumptions describing the effect of edge termination, the continuum scaling laws have no adjustable parameters and display correct asymptotic behavior. These scaling laws yield excellent agreement with experimental and numerical results for vibration frequencies in both isotropic and anisotropic structures as well as useful expressions for the frequency dependence on structure size and edge termination. △ Less

Submitted 17 October, 2019; originally announced October 2019.

Comments: accepted by Carbon (2019)

Optimized substrates and measurement approaches for Raman spectroscopy of graphene nanoribbons

Authors: Jan Overbeck, Gabriela Borin Barin, Colin Daniels, Mickael Perrin, Liangbo Liang, Oliver Braun, Rimah Darawish, Bryanna Burkhardt, Tim Dumslaff, Xiao-Ye Wang, Akimitsu Narita, Klaus Müllen, Vincent Meunier, Roman Fasel, Michel Calame, Pascal Ruffieux

Abstract: The on-surface synthesis of graphene nanoribbons (GNRs) allows for the fabrication of atomically precise narrow GNRs. Despite their exceptional properties which can be tuned by ribbon width and edge structure, significant challenges remain for GNR processing and characterization. In this contribution, we use Raman spectroscopy to characterize different types of GNRs on their growth substrate and t… ▽ More The on-surface synthesis of graphene nanoribbons (GNRs) allows for the fabrication of atomically precise narrow GNRs. Despite their exceptional properties which can be tuned by ribbon width and edge structure, significant challenges remain for GNR processing and characterization. In this contribution, we use Raman spectroscopy to characterize different types of GNRs on their growth substrate and to track their quality upon substrate transfer. We present a Raman-optimized (RO) device substrate and an optimized mapping approach that allows for acquisition of high-resolution Raman spectra, achieving enhancement factors as high as 120 with respect to signals measured on standard SiO2/Si substrates. We show that this approach is well-suited to routinely monitor the geometry-dependent low-frequency modes of GNRs. In particular, we track the radial breathing-like mode (RBLM) and the shear-like mode (SLM) for 5-, 7- and 9-atom wide armchair GNRs (AGNRs) and compare their frequencies with first-principles calculations. △ Less

Submitted 22 July, 2019; v1 submitted 3 July, 2019; originally announced July 2019.

Comments: *contributed equally 190722 - corrected typos

Quantifying the impact of network structure on speed and accuracy in collective decision-making

Authors: Bryan C. Daniels, Pawel Romanczuk

Abstract: Found in varied contexts from neurons to ants to fish, binary decision-making is one of the simplest forms of collective computation. In this process, information collected by individuals about an uncertain environment is accumulated to guide behavior at the aggregate scale. We study binary decision-making dynamics in networks responding to inputs with small signal-to-noise ratios, looking for qua… ▽ More Found in varied contexts from neurons to ants to fish, binary decision-making is one of the simplest forms of collective computation. In this process, information collected by individuals about an uncertain environment is accumulated to guide behavior at the aggregate scale. We study binary decision-making dynamics in networks responding to inputs with small signal-to-noise ratios, looking for quantitative measures of collectivity that control decision-making performance. We find that decision accuracy is controlled largely by three factors: the leading eigenvalue of the network adjacency matrix, the corresponding eigenvector's participation ratio, and distance from the corresponding symmetry-breaking bifurcation. This allows us to predict how decision-making performance scales in large networks based on their spectral properties. Specifically, we explore the effects of localization caused by the hierarchical assortative structure of a "rich club" topology. This gives insight into the tradeoffs involved in the higher-order structure found in living networks performing collective computations. △ Less

Submitted 22 March, 2019; originally announced March 2019.

Comments: 16 pages, 4 figures

Emergent regularities and scaling in armed conflict data

Authors: Edward D. Lee, Bryan C. Daniels, Christopher R. Myers, David C. Krakauer, Jessica C. Flack

Abstract: Armed conflict exhibits regularities beyond known power law distributions of fatalities and duration over varying culture and geography. We systematically cluster conflict reports from a database of $10^5$ events from Africa spanning 20 years into conflict avalanches. Conflict profiles collapse over a range of scales. Duration, diameter, extent, fatalities, and report totals satisfy mutually consi… ▽ More Armed conflict exhibits regularities beyond known power law distributions of fatalities and duration over varying culture and geography. We systematically cluster conflict reports from a database of $10^5$ events from Africa spanning 20 years into conflict avalanches. Conflict profiles collapse over a range of scales. Duration, diameter, extent, fatalities, and report totals satisfy mutually consistent scaling relations captured with a model combining geographic spread and local conflict-site growth. The emergence of such social scaling laws hints at principles guiding conflict evolution. △ Less

Submitted 29 April, 2020; v1 submitted 18 March, 2019; originally announced March 2019.

Automated, predictive, and interpretable inference of C. elegans escape dynamics

Authors: Bryan C. Daniels, William S. Ryu, Ilya Nemenman

Abstract: The roundworm C. elegans exhibits robust escape behavior in response to rapidly rising temperature. The behavior lasts for a few seconds, shows history dependence, involves both sensory and motor systems, and is too complicated to model mechanistically using currently available knowledge. Instead we model the process phenomenologically, and we use the Sir Isaac dynamical inference platform to infe… ▽ More The roundworm C. elegans exhibits robust escape behavior in response to rapidly rising temperature. The behavior lasts for a few seconds, shows history dependence, involves both sensory and motor systems, and is too complicated to model mechanistically using currently available knowledge. Instead we model the process phenomenologically, and we use the Sir Isaac dynamical inference platform to infer the model in a fully automated fashion directly from experimental data. The inferred model requires incorporation of an unobserved dynamical variable, and is biologically interpretable. The model makes accurate predictions about the dynamics of the worm behavior, and it can be used to characterize the functional logic of the dynamical system underlying the escape response. This work illustrates the power of modern artificial intelligence to aid in discovery of accurate and interpretable models of complex natural systems. △ Less

Submitted 25 September, 2018; originally announced September 2018.

Comments: 19 pages, 5 figures

Engineering of robust topological quantum phases in graphene nanoribbons

Authors: Oliver Gröning, Shiyong Wang, Xuelin Yao, Carlo A. Pignedoli, Gabriela Borin Barin, Colin Daniels, Andrew Cupo, Vincent Meunier, Xinliang Feng, Akimitsu Narita, Klaus Müllen, Pascal Ruffieux, Roman Fasel

Abstract: Here we present a flexible strategy to realize robust nanomaterials exhibiting valence electronic structures whose fundamental physics is described by the SSH-Hamiltonian. These solid-state materials are realized using atomically precise graphene nanoribbons (GNR). We demonstrate the controlled periodic coupling of topological boundary states at junctions of armchair GNRs of different widths to cr… ▽ More Here we present a flexible strategy to realize robust nanomaterials exhibiting valence electronic structures whose fundamental physics is described by the SSH-Hamiltonian. These solid-state materials are realized using atomically precise graphene nanoribbons (GNR). We demonstrate the controlled periodic coupling of topological boundary states at junctions of armchair GNRs of different widths to create quasi-1D trivial and non-trivial electronic quantum phases. Their topological class is experimentally determined by drawing upon the bulk-boundary correspondence and measuring the presence (non-trivial) or absence (trivial) of localized end states by scanning tunneling spectroscopy (STS). The strategy we propose has the potential to tune the band width of the topological electronic bands close to the energy scale of proximity induced spin-orbit coupling or superconductivity, and may allow the realization of Kitaev like Hamiltonians and Majorana type end states. △ Less

Submitted 17 May, 2018; originally announced May 2018.

Logic and connectivity jointly determine criticality in biological gene regulatory networks

Authors: Bryan C. Daniels, Hyunju Kim, Douglas Moore, Siyu Zhou, Harrison Smith, Bradley Karas, Stuart A. Kauffman, Sara I. Walker

Abstract: The complex dynamics of gene expression in living cells can be well-approximated using Boolean networks. The average sensitivity is a natural measure of stability in these systems: values below one indicate typically stable dynamics associated with an ordered phase, whereas values above one indicate chaotic dynamics. This yields a theoretically motivated adaptive advantage to being near the critic… ▽ More The complex dynamics of gene expression in living cells can be well-approximated using Boolean networks. The average sensitivity is a natural measure of stability in these systems: values below one indicate typically stable dynamics associated with an ordered phase, whereas values above one indicate chaotic dynamics. This yields a theoretically motivated adaptive advantage to being near the critical value of one, at the boundary between order and chaos. Here, we measure average sensitivity for 66 publicly available Boolean network models describing the function of gene regulatory circuits across diverse living processes. We find the average sensitivity values for these networks are clustered around unity, indicating they are near critical. In many types of random networks, mean connectivity <K> and the average activity bias of the logic functions <p> have been found to be the most important network properties in determining average sensitivity, and by extension a network's criticality. Surprisingly, many of these gene regulatory networks achieve the near-critical state with <K> and <p> far from that predicted for critical systems: randomized networks sharing the local causal structure and local logic of biological networks better reproduce their critical behavior than controlling for macroscale properties such as <K> and <p> alone. This suggests the local properties of genes interacting within regulatory networks are selected to collectively be near-critical, and this non-local property of gene regulatory network dynamics cannot be predicted using the density of interactions alone. △ Less

Submitted 3 May, 2018; originally announced May 2018.

Comments: 10 pages, 7 figures

Journal ref: Phys. Rev. Lett. 121, 138102 (2018)

Convenient Interface to Inverse Ising (ConIII): A Python 3 Package for Solving Ising-Type Maximum Entropy Models

Authors: Edward D. Lee, Bryan C Daniels

Abstract: ConIII (pronounced CON-ee) is an open-source Python project providing a simple interface to solving the pairwise and higher order Ising model and a base for extension to other maximum entropy models. We describe the maximum entropy problem and give an overview of the algorithms that are implemented as part of ConIII (https://github.com/eltrompetero/coniii) including Monte Carlo histogram, pseudoli… ▽ More ConIII (pronounced CON-ee) is an open-source Python project providing a simple interface to solving the pairwise and higher order Ising model and a base for extension to other maximum entropy models. We describe the maximum entropy problem and give an overview of the algorithms that are implemented as part of ConIII (https://github.com/eltrompetero/coniii) including Monte Carlo histogram, pseudolikelihood, minimum probability flow, a regularized mean field method, and a cluster expansion method. Our goal is to make a variety of maximum entropy techniques accessible to those unfamiliar with the techniques and accelerate workflow for users. △ Less

Submitted 10 March, 2019; v1 submitted 24 January, 2018; originally announced January 2018.

Investigation of Deuterium Loaded Materials Subject to X-Ray Exposure

Authors: Theresa L. Benyo, Bruce M. Steinetz, Robert C. Hendricks, Richard E. Martin, Lawrence P. Forsley, Christopher C. Daniels, Arnon Chait, Vladimir Pines, Marianna Pines, Nicholas Penney, Tracy R. Kamm, Michael D. Becks

Abstract: Results are presented from an exploratory study involving x-ray irradiation of select deuterated materials. Titanium deuteride (TiD2) plus deuterated polyethylene ([-CD2-]n; DPE), DPE alone, and for control, hydrogen-based polyethylene ([-CH2-]n; HPE) samples and nondeuterated titanium samples were exposed to x-ray irradiation. These samples were exposed to various energy levels from 65 to 280 kV… ▽ More Results are presented from an exploratory study involving x-ray irradiation of select deuterated materials. Titanium deuteride (TiD2) plus deuterated polyethylene ([-CD2-]n; DPE), DPE alone, and for control, hydrogen-based polyethylene ([-CH2-]n; HPE) samples and nondeuterated titanium samples were exposed to x-ray irradiation. These samples were exposed to various energy levels from 65 to 280 kV with prescribed electron flux from 500 to 9000 micro-A impinging on a tungsten braking target, with total exposure times ranging from 55 to 280 min. Gamma activity was measured using a high-purity germanium (HPGe) detector, and for all samples no gamma activity above background was detected. Alpha and beta activities were measured using a gas proportional counter, and for select samples beta activity was measured with a liquid scintillator spectrometer. The majority of the deuterated materials subjected to the microfocus x-ray irradiation exhibited postexposure beta activity above background and several showed short-lived alpha activity. The HPE and nondeuterated titanium control samples exposed to the x-ray irradiation showed no postexposure alpha or beta activities above background. Several of the samples (SL10A, SL16, SL17A) showed beta activity above background with a greater than 4-sigma confidence level, months after exposure. Portions of SL10A, SL16, and SL17A samples were also scanned using a beta scintillator and found to have beta activity in the tritium energy band, continuing without noticeable decay for over 12 months. Beta scintillation investigation of as-received materials (before x-ray exposure) showed no beta activity in the tritium energy band, indicating the beta emitters were not in the starting materials. △ Less

Submitted 3 April, 2017; originally announced April 2017.

Comments: 30 pages, 10 figures, Created as work for NASA

Report number: NASA/TM--2015-218491/REV1

Forward Collision Vehicular Radar with IEEE 802.11: Feasibility Demonstration through Measurements

Authors: Enoch R. Yeh, Robert C. Daniels, Robert W. Heath, Jr

Abstract: Increasing safety and automation in transportation systems has led to the proliferation of radar and IEEE 802.11 dedicated short range communication (DSRC) in vehicles. Current implementations of vehicular radar devices, however, are expensive, use a substantial amount of bandwidth, and are susceptible to multiple security risks. Consider the feasibility of using an IEEE 802.11 orthogonal frequenc… ▽ More Increasing safety and automation in transportation systems has led to the proliferation of radar and IEEE 802.11 dedicated short range communication (DSRC) in vehicles. Current implementations of vehicular radar devices, however, are expensive, use a substantial amount of bandwidth, and are susceptible to multiple security risks. Consider the feasibility of using an IEEE 802.11 orthogonal frequency division multiplexing (OFDM) communications waveform to perform radar functions. In this paper, we present an approach that determines the mean-normalized channel energy from frequency domain channel estimates and models it as a direct sinusoidal function of target range, enabling closest target range estimation. In addition, we propose an alternative to vehicular forward collision detection by extending IEEE 802.11 dedicated short-range communications (DSRC) and WiFi technology to radar, providing a foundation for joint communications and radar framework. Furthermore, we perform an experimental demonstration using existing IEEE 802.11 devices with minimal modification through algorithm processing on frequency-domain channel estimates. The results of this paper show that our solution delivers similar accuracy and reliability to mmWave radar devices with as little as 20 MHz of spectrum (doubling DSRC's 10 MHz allocation), indicating significant potential for industrial devices with joint vehicular communications and radar capabilities. △ Less

Submitted 5 June, 2017; v1 submitted 10 February, 2017; originally announced February 2017.

Control of critical behavior in a small-scale social system

Authors: Bryan C. Daniels, David C. Krakauer, Jessica C. Flack

Abstract: Over the last decade new technologies for making large numbers of fine-grained measurements have led to the surprising discovery that many biological systems sit near a critical point. These systems are potentially more adaptive in that small changes to component behavior can induce large-scale changes in aggregate structure and function. Accounting for criticality remains a challenge as sensitivi… ▽ More Over the last decade new technologies for making large numbers of fine-grained measurements have led to the surprising discovery that many biological systems sit near a critical point. These systems are potentially more adaptive in that small changes to component behavior can induce large-scale changes in aggregate structure and function. Accounting for criticality remains a challenge as sensitivity to perturbation suggests a lack of robustness. Furthermore, change induced by perturbation may not be adaptive. Complicating matters further critical phenomena can result from history-dependent stochastic processes. A question central to distinguishing among these conflicting views of criticality is to what degree criticality can be controlled by the components of the system. We address the control of criticality using data on conflict dynamics and fight sizes from an animal society model system (Macaca nemestrina, n=48). The system is fundamentally finite so we operationalize criticality in information theoretic terms using Fisher information and a measure of instability. We analyze criticality using empirically-grounded equilibrium (maximum entropy) and dynamic (branching process) models of the monkeys' fight-joining behavior. We find that (1) this heterogeneous, socially organized system, like homogeneous, spatial systems (flocks and schools), sits near a critical point, (2) the contributions individuals make to how critical the system is can be quantified and vary, and (3) the distance from the critical point (DFC) can be controlled through biologically plausible mechanisms operating on this heterogeneity. These mechanisms include third-party policing, which dampens fight participation of the individuals with the largest effect on DFC. Control of DFC allows biological systems to balance the tradeoff between robustness and need for rapid change. △ Less

Submitted 13 January, 2016; originally announced January 2016.

Comments: 34 pages, 8 figures

Percolation of even sites for enhanced random sequential adsorption

Authors: Christopher J. E. Daniels, Mathew D. Penrose

Abstract: Consider random sequential adsorption on a chequerboard lattice with arrivals at rate $1$ on light squares and at rate $λ$ on dark squares. Ultimately, each square is either occupied, or blocked by an occupied neighbour. Colour the occupied dark squares and blocked light sites {\em black}, and the remaining squares {\em white}. Independently at each meeting-point of four squares, allow diagonal co… ▽ More Consider random sequential adsorption on a chequerboard lattice with arrivals at rate $1$ on light squares and at rate $λ$ on dark squares. Ultimately, each square is either occupied, or blocked by an occupied neighbour. Colour the occupied dark squares and blocked light sites {\em black}, and the remaining squares {\em white}. Independently at each meeting-point of four squares, allow diagonal connections between black squares with probability $p$; otherwise allow diagonal connections between white squares. We show that there is a critical surface of pairs $(λ, p)$, containing the pair $(1,0.5)$, such that for $(λ, p)$ lying above (respectively, below) the critical surface the black (resp. white) phase percolates, and on the critical surface neither phase percolates. △ Less

Submitted 22 June, 2016; v1 submitted 31 August, 2015; originally announced August 2015.

Comments: 29 Pages

MSC Class: 60K35; 82B43

Sloppiness and Emergent Theories in Physics, Biology, and Beyond

Authors: Mark K. Transtrum, Benjamin Machta, Kevin Brown, Bryan C. Daniels, Christopher R. Myers, James P. Sethna

Abstract: Large scale models of physical phenomena demand the development of new statistical and computational tools in order to be effective. Many such models are `sloppy', i.e., exhibit behavior controlled by a relatively small number of parameter combinations. We review an information theoretic framework for analyzing sloppy models. This formalism is based on the Fisher Information Matrix, which we inter… ▽ More Large scale models of physical phenomena demand the development of new statistical and computational tools in order to be effective. Many such models are `sloppy', i.e., exhibit behavior controlled by a relatively small number of parameter combinations. We review an information theoretic framework for analyzing sloppy models. This formalism is based on the Fisher Information Matrix, which we interpret as a Riemannian metric on a parameterized space of models. Distance in this space is a measure of how distinguishable two models are based on their predictions. Sloppy model manifolds are bounded with a hierarchy of widths and extrinsic curvatures. We show how the manifold boundary approximation can extract the simple, hidden theory from complicated sloppy models. We attribute the success of simple effective models in physics as likewise emerging from complicated processes exhibiting a low effective dimensionality. We discuss the ramifications and consequences of sloppy models for biochemistry and science more generally. We suggest that the reason our complex world is understandable is due to the same fundamental reason: simple theories of macroscopic behavior are hidden inside complicated microscopic processes. △ Less

Submitted 30 January, 2015; originally announced January 2015.

Efficient inference of parsimonious phenomenological models of cellular dynamics using S-systems and alternating regression

Authors: Bryan C. Daniels, Ilya Nemenman

Abstract: The nonlinearity of dynamics in systems biology makes it hard to infer them from experimental data. Simple linear models are computationally efficient, but cannot incorporate these important nonlinearities. An adaptive method based on the S-system formalism, which is a sensible representation of nonlinear mass-action kinetics typically found in cellular dynamics, maintains the efficiency of linear… ▽ More The nonlinearity of dynamics in systems biology makes it hard to infer them from experimental data. Simple linear models are computationally efficient, but cannot incorporate these important nonlinearities. An adaptive method based on the S-system formalism, which is a sensible representation of nonlinear mass-action kinetics typically found in cellular dynamics, maintains the efficiency of linear regression. We combine this approach with adaptive model selection to obtain efficient and parsimonious representations of cellular dynamics. The approach is tested by inferring the dynamics of yeast glycolysis from simulated data. With little computing time, it produces dynamical models with high predictive power and with structural complexity adapted to the difficulty of the inference problem. △ Less

Submitted 16 June, 2014; originally announced June 2014.

Comments: 14 pages, 2 figures

Automated adaptive inference of coarse-grained dynamical models in systems biology

Authors: Bryan C. Daniels, Ilya Nemenman

Abstract: Cellular regulatory dynamics is driven by large and intricate networks of interactions at the molecular scale, whose sheer size obfuscates understanding. In light of limited experimental data, many parameters of such dynamics are unknown, and thus models built on the detailed, mechanistic viewpoint overfit and are not predictive. At the other extreme, simple ad hoc models of complex processes ofte… ▽ More Cellular regulatory dynamics is driven by large and intricate networks of interactions at the molecular scale, whose sheer size obfuscates understanding. In light of limited experimental data, many parameters of such dynamics are unknown, and thus models built on the detailed, mechanistic viewpoint overfit and are not predictive. At the other extreme, simple ad hoc models of complex processes often miss defining features of the underlying systems. Here we propose an approach that instead constructs phenomenological, coarse-grained models of network dynamics that automatically adapt their complexity to the amount of available data. Such adaptive models lead to accurate predictions even when microscopic details of the studied systems are unknown due to insufficient data. The approach is computationally tractable, even for a relatively large number of dynamical variables, allowing its software realization, named Sir Isaac, to make successful predictions even when important dynamic variables are unobserved. For example, it matches the known phase space structure for simulated planetary motion data, avoids overfitting in a complex biological signaling system, and produces accurate predictions for a yeast glycolysis model with only tens of data points and over half of the interacting species unobserved. △ Less

Submitted 24 April, 2014; originally announced April 2014.

Comments: 38 pages, 12 figures

Nucleation at the DNA supercoiling transition

Authors: Bryan C. Daniels, James P. Sethna

Abstract: Twisting DNA under a constant applied force reveals a thermally activated transition into a state with a supercoiled structure known as a plectoneme. Using transition state theory, we predict the rate of this plectoneme nucleation to be of order 10^4 Hz. We reconcile this with experiments that have measured hopping rates of order 10 Hz by noting that the viscosity of the bead used to manipulate th… ▽ More Twisting DNA under a constant applied force reveals a thermally activated transition into a state with a supercoiled structure known as a plectoneme. Using transition state theory, we predict the rate of this plectoneme nucleation to be of order 10^4 Hz. We reconcile this with experiments that have measured hopping rates of order 10 Hz by noting that the viscosity of the bead used to manipulate the DNA limits the measured rate. We find that the intrinsic bending caused by disorder in the base-pair sequence is important for understanding the free energy barrier that governs the transition. Both analytic and numerical methods are used in the calculations. We provide extensive details on the numerical methods for simulating the elastic rod model with and without disorder. △ Less

Submitted 14 March, 2011; v1 submitted 12 October, 2010; originally announced October 2010.

Comments: 18 pages, 15 figures

Discontinuities at the DNA supercoiling transition

Authors: Bryan C. Daniels, Scott Forth, Maxim Y. Sheinin, Michelle D. Wang, James P. Sethna

Abstract: While slowly turning the ends of a single molecule of DNA at constant applied force, a discontinuity was recently observed at the supercoiling transition, when a small plectoneme is suddenly formed. This can be understood as an abrupt transition into a state in which stretched and plectonemic DNA coexist. We argue that there should be discontinuities in both the extension and the torque at the t… ▽ More While slowly turning the ends of a single molecule of DNA at constant applied force, a discontinuity was recently observed at the supercoiling transition, when a small plectoneme is suddenly formed. This can be understood as an abrupt transition into a state in which stretched and plectonemic DNA coexist. We argue that there should be discontinuities in both the extension and the torque at the transition, and provide experimental evidence for both. To predict the sizes of these discontinuities and how they change with the overall length of DNA, we organize a theory for the coexisting plectonemic state in terms of four length-independent parameters. We also test plectoneme theories, including our own elastic rod simulation, finding discrepancies with experiment that can be understood in terms of the four coexisting state parameters. △ Less

Submitted 21 July, 2009; v1 submitted 21 November, 2008; originally announced November 2008.

Comments: 11 pages, 5 figures; revised version, with added supplemental material

Statistical properties of multistep enzyme-mediated reactions

Authors: Wiet H. de Ronde, Bryan C. Daniels, Andrew Mugler, Nikolai A. Sinitsyn, Ilya Nemenman

Abstract: Enzyme-mediated reactions may proceed through multiple intermediate conformational states before creating a final product molecule, and one often wishes to identify such intermediate structures from observations of the product creation. In this paper, we address this problem by solving the chemical master equations for various enzymatic reactions. We devise a perturbation theory analogous to tha… ▽ More Enzyme-mediated reactions may proceed through multiple intermediate conformational states before creating a final product molecule, and one often wishes to identify such intermediate structures from observations of the product creation. In this paper, we address this problem by solving the chemical master equations for various enzymatic reactions. We devise a perturbation theory analogous to that used in quantum mechanics that allows us to determine the first (<n>) and the second (variance) cumulants of the distribution of created product molecules as a function of the substrate concentration and the kinetic rates of the intermediate processes. The mean product flux V=d<n>/dt (or "dose-response" curve) and the Fano factor F=variance/<n> are both realistically measurable quantities, and while the mean flux can often appear the same for different reaction types, the Fano factor can be quite different. This suggests both qualitative and quantitative ways to discriminate between different reaction schemes, and we explore this possibility in the context of four sample multistep enzymatic reactions. We argue that measuring both the mean flux and the Fano factor can not only discriminate between reaction types, but can also provide some detailed information about the internal, unobserved kinetic rates, and this can be done without measuring single-molecule transition events. △ Less

Submitted 20 November, 2008; originally announced November 2008.

Comments: 8 pages, 3 figures

Journal ref: IET Sys. Bio. 3 (2009) 429-437

Sloppiness, robustness, and evolvability in systems biology

Authors: Bryan C. Daniels, Yan-Jiun Chen, James P. Sethna, Ryan N. Gutenkunst, Christopher R. Myers

Abstract: The functioning of many biochemical networks is often robust -- remarkably stable under changes in external conditions and internal reaction parameters. Much recent work on robustness and evolvability has focused on the structure of neutral spaces, in which system behavior remains invariant to mutations. Recently we have shown that the collective behavior of multiparameter models is most often '… ▽ More The functioning of many biochemical networks is often robust -- remarkably stable under changes in external conditions and internal reaction parameters. Much recent work on robustness and evolvability has focused on the structure of neutral spaces, in which system behavior remains invariant to mutations. Recently we have shown that the collective behavior of multiparameter models is most often 'sloppy': insensitive to changes except along a few 'stiff' combinations of parameters, with an enormous sloppy neutral subspace. Robustness is often assumed to be an emergent evolved property, but the sloppiness natural to biochemical networks offers an alternative non-adaptive explanation. Conversely, ideas developed to study evolvability in robust systems can be usefully extended to characterize sloppy systems. △ Less

Submitted 19 May, 2008; originally announced May 2008.

Comments: 12 pages, 7 figures, submitted to Current Opinion in Biotechnology

Gamma-Rays Produced in Cosmic-Ray Interactions and the TeV-band Spectrum of RX J1713.7-3946

Authors: C. -Y. Huang, S. -E. Park, M. Pohl, C. D. Daniels

Abstract: We employ the Monte Carlo particle collision code DPMJET3.04 to determine the multiplicity spectra of various secondary particles (in addition to $π^0$'s) with $γ$'s as the final decay state, that are produced in cosmic-ray ($p$'s and $α$'s) interactions with the interstellar medium. We derive an easy-to-use $γ$-ray production matrix for cosmic rays with energies up to about 10 PeV. This $γ$-ray… ▽ More We employ the Monte Carlo particle collision code DPMJET3.04 to determine the multiplicity spectra of various secondary particles (in addition to $π^0$'s) with $γ$'s as the final decay state, that are produced in cosmic-ray ($p$'s and $α$'s) interactions with the interstellar medium. We derive an easy-to-use $γ$-ray production matrix for cosmic rays with energies up to about 10 PeV. This $γ$-ray production matrix is applied to the GeV excess in diffuse Galactic $γ$-rays observed by EGRET, and we conclude the non-$π^0$ decay components are insufficient to explain the GeV excess, although they have contributed a different spectrum from the $π^0$-decay component. We also test the hypothesis that the TeV-band $γ$-ray emission of the shell-type SNR RX J1713.7-3946 observed with HESS is caused by hadronic cosmic rays which are accelerated by a cosmic-ray modified shock. By the $χ^2$ statistics, we find a continuously softening spectrum is strongly preferred, in contrast to expectations. A hardening spectrum has about 1% probability to explain the HESS data, but then only if a hard cutoff at 50-100 TeV is imposed on the particle spectrum. △ Less

Submitted 4 June, 2007; originally announced June 2007.

Comments: 3 pages; 4 figures; Contribution to the First GLAST Symposium, Standord, 2007

Journal ref: AIPConf.Proc.921:169-171,2007

Diffuse Gamma-Rays Produced in Cosmic-Ray Interactions and the TeV-band Spectrum of RX J1713.7-3946

Authors: C. -Y. Huang, M. Pohl, S. -E. Park, C. D. Daniels

Abstract: We employ the Monte Carlo particle collision code DPMJET3.04 to determine the multiplicity spectra of various secondary particles (in addition to $π^0$'s) with $γ$'s as the final decay state, that are produced in cosmic-ray ($p$'s and $α$'s) interactions with the interstellar medium. We derive an easy-to-use $γ$-ray production matrix for cosmic rays with energies up to about 10 PeV. This $γ$-ray… ▽ More We employ the Monte Carlo particle collision code DPMJET3.04 to determine the multiplicity spectra of various secondary particles (in addition to $π^0$'s) with $γ$'s as the final decay state, that are produced in cosmic-ray ($p$'s and $α$'s) interactions with the interstellar medium. We derive an easy-to-use $γ$-ray production matrix for cosmic rays with energies up to about 10 PeV. This $γ$-ray production matrix is applied to the GeV excess in diffuse Galactic $γ$-rays observed by EGRET, and we conclude the non-$π^0$ decay components are insufficient to explain the GeV excess, although they have contributed a different spectrum from the $π^0$-decay component. We also test the hypothesis that the TeV-band $γ$-ray emission of the shell-type SNR RX J1713.7-3946 observed with HESS is caused by hadronic cosmic rays which are accelerated by a cosmic-ray modified shock. By the $χ^2$ statistics, we find a continuously softening spectrum is strongly preferred, in contrast to expectations. A hardening spectrum has about 1% probability to explain the HESS data, but then only if a hard cutoff at 50-100 TeV is imposed on the particle spectrum. △ Less

Submitted 5 June, 2007; v1 submitted 4 June, 2007; originally announced June 2007.

Comments: 4 pages; 5 figures; Contribution to the 30th ICRC, Merida, Mexico, 2007; References included

Gamma-Rays Produced in Cosmic-Ray Interactions and the TeV-band Spectrum of RX J1713-3946

Authors: C. -Y. Huang, S. -E. Park, M. Pohl, C. D. Daniels

Abstract: In this work we study the individual contribution to diffuse $γ$-ray emission from the secondary products in hadronic interactions generated by cosmic rays (CRs), in addition to the contribution of $π^0$ decay via the decay mode $π^0 \to 2γ$. For that purpose we employ the Monte Carlo particle collision code DPMJET3.04 to determine the multiplicity spectra of various secondary particles with… ▽ More In this work we study the individual contribution to diffuse $γ$-ray emission from the secondary products in hadronic interactions generated by cosmic rays (CRs), in addition to the contribution of $π^0$ decay via the decay mode $π^0 \to 2γ$. For that purpose we employ the Monte Carlo particle collision code DPMJET3.04 to determine the multiplicity spectra of various secondary particles with $γ$'s as the final decay state, that result from inelastic collisions between cosmic-ray protons and Helium nuclei and the interstellar medium with standard composition. We thus derive an easy-to-use $γ$-ray production matrix for cosmic ray up to about 10 PeV, that can be used to interpret the $γ$-ray spectra of diffuse galactic emission and supernova remnants (SNR). We apply the $γ$-ray production matrix to the GeV excess in diffuse galactic $γ$-rays that was seen with EGRET. Although the non-$π^0$ contributions to the total emission have a different spectrum than the $π^0$-decay component, they are insufficient to explain the GeV excess. We also test the hypothesis that the TeV-band $γ$-ray emission of the shell-type SNR RX J1713-3946, that was observed with HESS, is caused by shock-accelerated hadronic cosmic rays. This scenario implies a very high efficacy of particle acceleration, so the particle spectrum is expected to continuously harden toward high energies on account of cosmic-ray modification of the shock. Using the $χ^2$ statistic we find that a continuously softening spectrum is strongly preferred, in contrast to expectations. A hardening spectrum has about 1% probability to explain the HESS data, but then only if a hard cut-off at 50-100 TeV is imposed on the particle spectrum. △ Less

Submitted 11 January, 2007; v1 submitted 28 November, 2006; originally announced November 2006.

Comments: minor revision; Astropart. Phys. in press

Journal ref: Astropart.Phys.27:429-439,2007

Absence of lattice coherence effects in Ce0.6La0.4Pb3: A magnetic field study

Authors: R. Pietri, C. R. Rotundu, B. Andraka, B. C. Daniels, K. Ingersent

Abstract: The specific heat of polycrystalline Ce0.6La0.4Pb3 has been measured in magnetic fields ranging from 0 to 14 T. After subtraction of a lattice contribution, the specific heat between 1 K and 10 K is well described by the S=1/2 single-impurity Kondo model with just one adjustable parameter: the zero-field Kondo temperature. In particular, the variation in the temperature and the height of the pea… ▽ More The specific heat of polycrystalline Ce0.6La0.4Pb3 has been measured in magnetic fields ranging from 0 to 14 T. After subtraction of a lattice contribution, the specific heat between 1 K and 10 K is well described by the S=1/2 single-impurity Kondo model with just one adjustable parameter: the zero-field Kondo temperature. In particular, the variation in the temperature and the height of the peak in C vs T is captured with good accuracy. This fit suggests that lattice coherence effects play no significant role in the magnetic-field response of this concentrated Kondo system. △ Less

Submitted 19 May, 2005; originally announced May 2005.

Comments: 4 pages, 2 figures, Proceedings of MMM'04

Journal ref: J. Appl. Phys. 97, 10A510 (2005)