The OpenKIM processing pipeline: A cloud-based automatic material property computation engine.

The Open Knowledgebase of Interatomic Models (OpenKIM) is a framework intended to facilitate access to standardized implementations of interatomic models for molecular simulations along with computational protocols to evaluate them. These protocols include tests to compute material properties predicted by models and verification checks to assess their coding integrity. While housing this content in a unified, publicly available environment constitutes a major step forward for the molecular modeling community, it further presents the opportunity to understand the range of validity of interatomic models and their suitability for specific target applications. To this end, OpenKIM includes a computational pipeline that runs tests and verification checks using all available interatomic models contained within the OpenKIM Repository at https://openkim.org. The OpenKIM Processing Pipeline is built on a set of Docker images hosted on distributed, heterogeneous hardware and utilizes open-source software to automatically run test-model and verification check-model pairs and resolve dependencies between them. The design philosophy and implementation choices made in the development of the pipeline are discussed as well as an example of its application to interatomic model selection.

Intra-unit-cell electronic nematicity of the high-T(c) copper-oxide pseudogap states.

In the high-transition-temperature (high-T(c)) superconductors the pseudogap phase becomes predominant when the density of doped holes is reduced. Within this phase it has been unclear which electronic symmetries (if any) are broken, what the identity of any associated order parameter might be, and which microscopic electronic degrees of freedom are active. Here we report the determination of a quantitative order parameter representing intra-unit-cell nematicity: the breaking of rotational symmetry by the electronic structure within each CuO(2) unit cell. We analyse spectroscopic-imaging scanning tunnelling microscope images of the intra-unit-cell states in underdoped Bi(2)Sr(2)CaCu(2)O(8 +) (delta) and, using two independent evaluation techniques, find evidence for electronic nematicity of the states close to the pseudogap energy. Moreover, we demonstrate directly that these phenomena arise from electronic differences at the two oxygen sites within each unit cell. If the characteristics of the pseudogap seen here and by other techniques all have the same microscopic origin, this phase involves weak magnetic states at the O sites that break 90 degrees -rotational symmetry within every CuO(2) unit cell.

Optimal experimental design in an epidermal growth factor receptor signalling and down-regulation model.

We apply the methods of optimal experimental design to a differential equation model for epidermal growth factor receptor signalling, trafficking and down-regulation. The model incorporates the role of a recently discovered protein complex made up of the E3 ubiquitin ligase, Cbl, the guanine exchange factor (GEF), Cool-1 (beta -Pix) and the Rho family G protein Cdc42. The complex has been suggested to be important in disrupting receptor down-regulation. We demonstrate that the model interactions can accurately reproduce the experimental observations, that they can be used to make predictions with accompanying uncertainties, and that we can apply ideas of optimal experimental design to suggest new experiments that reduce the uncertainty on unmeasurable components of the system.

Bayesian error estimation in density-functional theory.

We present a practical scheme for performing error estimates for density-functional theory calculations. The approach, which is based on ideas from Bayesian statistics, involves creating an ensemble of exchange-correlation functionals by comparing with an experimental database of binding energies for molecules and solids. Fluctuations within the ensemble can then be used to estimate errors relative to experiment on calculated quantities such as binding energies, bond lengths, and vibrational frequencies. It is demonstrated that the error bars on energy differences may vary by orders of magnitude for different systems in good agreement with existing experience.

Vibration-assisted electron tunneling in C140 transistors.

We measure electron tunneling in transistors made from C(140), a molecule with a mass-spring-mass geometry chosen as a model system to study electron-vibration coupling. We observe vibration-assisted tunneling at an energy corresponding to the stretching mode of C(140). Molecular modeling provides explanations for why this mode couples more strongly to electron tunneling than to the other internal modes of the molecule. We make comparisons between the observed tunneling rates and those expected from the Franck-Condon model.

The statistical mechanics of complex signaling networks: nerve growth factor signaling.

The inherent complexity of cellular signaling networks and their importance to a wide range of cellular functions necessitates the development of modeling methods that can be applied toward making predictions and highlighting the appropriate experiments to test our understanding of how these systems are designed and function. We use methods of statistical mechanics to extract useful predictions for complex cellular signaling networks. A key difficulty with signaling models is that, while significant effort is being made to experimentally measure the rate constants for individual steps in these networks, many of the parameters required to describe their behavior remain unknown or at best represent estimates. To establish the usefulness of our approach, we have applied our methods toward modeling the nerve growth factor (NGF)-induced differentiation of neuronal cells. In particular, we study the actions of NGF and mitogenic epidermal growth factor (EGF) in rat pheochromocytoma (PC12) cells. Through a network of intermediate signaling proteins, each of these growth factors stimulates extracellular regulated kinase (Erk) phosphorylation with distinct dynamical profiles. Using our modeling approach, we are able to predict the influence of specific signaling modules in determining the integrated cellular response to the two growth factors. Our methods also raise some interesting insights into the design and possible evolution of cellular systems, highlighting an inherent property of these systems that we call 'sloppiness.'

Fracture in mode I using a conserved phase-field model.

We present a continuum phase-field model of crack propagation. It includes a phase-field that is proportional to the mass density and a displacement field that is governed by linear elastic theory. Generic macroscopic crack growth laws emerge naturally from this model. In contrast to classical continuum fracture mechanics simulations, our model avoids numerical front tracking. The added phase-field smooths the sharp interface, enabling us to use equations of motion for the material (grounded in basic physical principles) rather than for the interface (which often are deduced from complicated theories or empirical observations). The interface dynamics thus emerges naturally. In this paper, we look at stationary solutions of the model, mode I fracture, and also discuss numerical issues. We find that the Griffith's threshold underestimates the critical value at which our system fractures due to long wavelength modes excited by the fracture process.

Calculation of quantum tunneling for a spatially extended defect: the dislocation kink in copper has a low effective mass.

Several experiments indicate that there are atomic tunneling defects in plastically deformed metals. How this is possible has not been clear, given the large mass of the metal atoms. Using a classical molecular-dynamics calculation, we determine the structures, energy barriers, effective masses, and quantum tunneling rates for dislocation kinks and jogs in copper screw dislocations. We find that jogs are unlikely to tunnel, but the kinks should have large quantum fluctuations. The kink motion involves hundreds of atoms each shifting a tiny amount, leading to a small effective mass and tunneling barrier.

Crackling noise.

Crackling noise arises when a system responds to changing external conditions through discrete, impulsive events spanning a broad range of sizes. A wide variety of physical systems exhibiting crackling noise have been studied, from earthquakes on faults to paper crumpling. Because these systems exhibit regular behaviour over a huge range of sizes, their behaviour is likely to be independent of microscopic and macroscopic details, and progress can be made by the use of simple models. The fact that these models and real systems can share the same behaviour on many scales is called universality. We illustrate these ideas by using results for our model of crackling noise in magnets, explaining the use of the renormalization group and scaling collapses, and we highlight some continuing challenges in this still-evolving field.