Parameter estimation by nonlinear least-squares minimization is a common problem that has an elegant geometric interpretation: the possible parameter values of a model induce a manifold within the space of data predictions. The minimization problem is then to find the point on the manifold closest to the experimental data. We show that the model manifolds of a large class of models, known as sloppy models, have many universal features; they are characterized by a geometric series of widths, extrinsic curvatures, and parameter-effect curvatures, which we describe as a hyper-ribbon.

We provide a minimal continuum model for mesoscale plasticity, explaining the cellular dislocation structures observed in deformed crystals. Our dislocation density tensor evolves from random, smooth initial conditions to form self-similar structures strikingly similar to those seen experimentally-reproducing both the fractal morphologies and some features of the scaling of cell sizes and misorientations analyzed experimentally.

In the high-transition-temperature (high-Tc) 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 CuO2 unit cell.

Fitting model parameters to experimental data is a common yet often challenging task, especially if the model contains many parameters. Typically, algorithms get lost in regions of parameter space in which the model is unresponsive to changes in parameters, and one is left to make adjustments by hand. We explain this difficulty by interpreting the fitting process as a generalized interpolation procedure. By considering the manifold of all model predictions in data space, we find that cross sections have a hierarchy of widths and are typically very narrow.

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.

Intergranular fracture in polycrystals is often simulated by finite elements coupled to a cohesive-zone model for the interfaces, requiring cohesive laws for grain boundaries as a function of their geometry. We discuss three challenges in understanding intergranular fracture in polycrystals. First, 3D grain boundary geometries comprise a five dimensional space. Second, the energy and peak stress of grain boundaries have singularities for all commensurate grain boundaries, especially those with short repeat distances.

Magneto-optical methods allow us to observe the dynamics of domain wall motion, but this is intrinsically a very noisy process. We discuss a new method allowing us to reduce the measurement noise, taking advantage of the acquisition of a whole temporal sequence of images. The resulting avalanche distributions give interesting hints as to the magnetization dynamics, but are strongly dependent on the size of the observation windows chosen. We investigate the effects of window size by studying finite-size scaling, and use this to extract the fractal dimension critical exponent 1/σν.

We study the metastability of the superheated Meissner state in type II superconductors with κ≫1 beyond Ginzburg-Landau theory, which is applicable only in the vicinity of the critical temperature. Within Eilenberger's semiclassical approximation, we use the local electrodynamic response of the superconductor to derive a generalized thermodynamic potential valid at any temperature. The stability analysis of this functional yields the temperature dependence of the superheating field.

We demonstrate the use of a variational method to determine a quantitative lower bound on the rate of convergence of Markov chain Monte Carlo (MCMC) algorithms as a function of the target density and proposal density. The bound relies on approximating the second largest eigenvalue in the spectrum of the MCMC operator using a variational principle and the approach is applicable to problems with continuous state spaces.

Are the cohesive laws of interfaces sufficient for modelling fracture in polycrystals using the cohesive zone model? We examine this question by comparing a fully atomistic simulation of a silicon polycrystal with a finite element simulation with a similar overall geometry. The cohesive laws used in the finite element simulation are measured atomistically. We describe in detail how to convert the output of atomistic grain boundary fracture simulations into the piecewise linear form needed by a cohesive zone model.