Publications

The resolution of geometric frustration in systems with continuous degrees of freedom often involves a cooperative inhomogeneous response and superextensive energy scaling. In contrast, the frustration in frustrated Ising-like spin systems is resolved uniformly. In this work we bridge between these two extremes by studying a frustrated model composed of N-state spins, and varying N. The expected cooperative response, observed for large N, is strongly attenuated as N is reduced, in a nontrivial way. Moderate N values show unique topological-like phases not observed before in frustrated models.

The geometry and topology of the region in which a director field is embedded impose limitations on the kind of supported orientational order. These limitations manifest as compatibility conditions that relate the quantities describing the director field to the geometry of the embedding space. For example, in two dimensions the splay and bend fields suffice to determine a director uniquely (up to rigid motions) and must comply with one relation linear in the Gaussian curvature of the embedding manifold. In 3D there are additional local fields describing the director, i.e. fields available to a local observer residing within the material, and a number of distinct ways to yield geometric frustration. So far it was unknown how many such local fields are required to uniquely describe a 3D director field, nor what are the compatibility relations they must satisfy. In this work, we address these questions directly. We employ the method of moving frames to show that a director field is fully determined by five local fields. These fields are shown to be related to each other and to the curvature of the embedding space through six differential relations. As an application of our method, we characterize all uniform distortion director fields, i.e., directors for which all the local characterizing fields are constant in space, in manifolds of constant curvature. The classification of such phases has been recently provided for directors in Euclidean space, where the textures correspond to foliations of space by parallel congruent helices. For non-vanishing curvature, we show that the pure twist phase is the only solution in positively curved space, while in the hyperbolic space uniform distortion fields correspond to foliations of space by (non-necessarily parallel) congruent helices. Further analysis of the obtained compatibility fields is expected to allow to also construct new non-uniform director fields.

Determining the stability of a viscoelastic structure is a difficult task. Seemingly stable conformations of viscoelastic structures may gradually creep until their stability is lost, while a discernible creeping in viscoelastic solids does not necessarily lead to instability. In lieu of theoretical predictive tools for viscoelastic instabilities, we are presently limited to numerical simulation to predict future stability. In this work, we describe viscoelastic solids through a temporally evolving instantaneous reference metric with respect to which elastic strains are measured. We show that for incompressible viscoelastic solids, this transparent and intuitive description allows to reduce the question of future stability to static calculations. We demonstrate the predictive power of the approach by elucidating the subtle mechanism of delayed instability in thin elastomeric shells, showing quantitative agreement with experiments.

More than one-quarter of molecular crystals that are able to be melted can be made to grow in the form of twisted lamellae or fibers. The mechanisms leading to such unusual crystal morphologies lacking long-range translational symmetry on the mesoscale are poorly understood. Benzil (C6H5C(O)-C(O)C6H5) is one such crystal. Here, we calculate the morphology of rod-shaped benzil nanocrystals and other related structures. The ground states of these ensembles were twisted by 0.05-0.75°/Å for rods with cross sections of 50-10 nm2, respectively; the degree of twisting decreased inversely proportional to the crystal cross-sectional area. In the aggregate, our computational studies, combined with earlier observations by light microscopy, suggest that in some cases very small crystals acquire 3D translational periodicity only after reaching a certain size. Twisting is accompanied by conformational changes of molecules on the {101¯ 0} surfaces of the six-sided rods, although it is not easily answered from our data whether such changes are causes of the twisting, consequences of surface stress where symmetry is broken, or consequences of intrinsic dissymmetry when two or more geometric tendencies are in conflict. Nevertheless, it has become clear that, in some cases, the development of a crystal with a lattice having long-range translational symmetry is not foretold in the thermodynamics of aggregates of molecules. Rather, a lattice is sometimes a device for allowing a growing crystal to take advantage of the thermodynamic driving force of growth, the best compromise for a large number of molecules, which on a smaller scale would be dissymmetric (have a point symmetry only). The relationship between these calculations and the ubiquity of crystal twisting on the mesoscale are discussed.

Plant photosynthetic (thylakoid) membranes are organized into complex networks that are differentiated into 2 distinct morphological and functional domains called grana and stroma lamellae. How the 2 domains join to form a continuous lamellar system has been the subject of numerous studies since the mid-1950s. Using different electron tomography techniques, we found that the grana and stroma lamellae are connected by an array of pitch-balanced right- and left-handed helical membrane surfaces of different radii and pitch. Consistent with theoretical predictions, this arrangement is shown to minimize the surface and bending energies of the membranes. Related configurations were proposed to be present in the rough endoplasmic reticulum and in dense nuclear matter phases theorized to exist in neutron star crusts, where the right- and left-handed helical elements differ only in their handedness. Pitch-balanced helical elements of alternating handedness may thus constitute a fundamental geometry for the efficient packing of connected layers or sheets.

Symmetry considerations preclude the possibility of twist or continuous helical symmetry in bulk crystalline structures. However, as has been shown nearly a century ago, twisted molecular crystals are ubiquitous and can be formed by about 1/4 of organic substances. Despite its ubiquity, this phenomenon has so far not been satisfactorily explained. In this work we study twisted molecular crystals as geometrically frustrated assemblies. We model the molecular constituents as uniaxially twisted cubes and examine their crystalline assembly. We exploit a renormalization group (RG) approach to follow the growth of the rod-like twisted crystals these constituents produce, inquiring in every step into the evolution of their morphology, response functions and residual energy. The gradual untwisting of the rod-like frustrated crystals predicted by the RG approach is verified experimentally using silicone rubber models of similar geometry. Our theory provides a mechanism for the conveyance of twist across length-scales observed experimentally and reconciles the apparent paradox of a twisted single crystal as a finite size effect.

Bent core (or banana shaped) liquid-crystal-forming-molecules locally favor an ordered state of zero splay and constant bend. Such a state, however, cannot be realized in the plane and the resulting liquid-crystalline phase is frustrated and must exhibit some compromise of these two mutually contradicting local intrinsic tendencies. This constitutes one of the most well-studied examples in which the intrinsic geometry of the constituents of a material gives rise to a geometrically frustrated assembly. Such geometric frustration is not only natural and ubiquitous but also leads to a striking variety of morphologies of ground states and exotic response properties. In this work we establish the necessary and sufficient conditions for two scalar functions, s and b to describe the splay and bend of a director field in the plane. We generalize these compatibility conditions for geometries with non-vanishing constant Gaussian curvature, and provide a reconstruction formula for the director field depending only on the splay and bend fields and their derivatives. Finally, we discuss optimal compromises for simple incompatible cases where the locally preferred values of the splay and bend cannot be simultaneously achieved.

A crumpled sheet of paper displays an intricate pattern of creases and point-like singular structures, termed d-cones. It is typically assumed that elongated creases form when ridges connecting two d-cones fold beyond the material yielding threshold, and scarring is thus a by-product of the folding dynamics that seek to minimize elastic energy. Here we show that rather than merely being the consequence of folding, plasticity can act as its instigator. We introduce and characterize a different type of crease that is inherently plastic and is formed by the propagation of a single point defect. When a pre-existing d-cone is strained beyond a certain threshold, the singular structure at its apex sharpens abruptly. The resulting focusing of strains yields the material just ahead of the singularity, allowing it to propagate, leaving a furrow-like scar in its wake. We suggest an intuitive fracture analogue to explain the creation of furrows.

We examine the shape change of a thin disk with an inserted wedge of material when it is pushed against a plane, using analytical, numerical, and experimental methods. Such sheets occur in packaging, surgery, and nanotechnology. We approximate the sheet as having vanishing strain, so that it takes a conical form in which straight generators converge to a disclination singularity. Then, its shape is that which minimizes elastic bending energy alone. Real sheets are expected to approach this limiting shape as their thickness approaches zero. The planar constraint forces a sector of the sheet to buckle into the third dimension. We find that the unbuckled sector is precisely semicircular, independent of the angle δ of the inserted wedge. We generalize the analysis to include conical as well as planar constraints and thereby establish a law of corresponding states for shallow cones of slope ε and thin wedges. In this regime, the single parameter δ/ε2 determines the shape. We discuss the singular limit in which the cone becomes a plane, and the unexpected slow convergence to the semicircular buckling observed in real sheets.

This review compares the conceptualization and practice of early real-space renormalization group methods with the conceptualization of more recent real-space transformations based on tensor networks. For specificity, it focuses upon two basic methods: the "potential-moving" approach most used in the period 1975 1980 and the "rewiring method" as it has been developed in the last five years. The newer method, part of a development called the tensor renormalization group, was originally based on principles of quantum entanglement. It is specialized for computing approximations for tensor products constituting, for example, the free energy or the ground state energy of a large system. It can attack a wide variety of problems, including quantum problems, which would otherwise be intractable. The older method is formulated in terms of spin variables and permits a straightforward construction and analysis of fixed points in rather transparent terms. However, in the form described here it is unsystematic, offers no path for improvement, and of unknown reliability. The new method is formulated in terms of index variables which may be considered as linear combinations of the statistical variables. Free energies emerge naturally, but fixed points are more subtle. Further, physical interpretations of the index variables are often elusive due to a gauge symmetry which allows only selected combinations of tensor entries to have physical significance. In applications, both methods employ analyses with varying degrees of complexity. The complexity is parametrized by an integer called chi (or D in the recent literature). Both methods are examined in action by using them to compute fixed points related to Ising models for small values of the complexity parameter. They behave quite differently. The old method gives a reasonably good picture of the fixed point, as measured, for example, by the accuracy of the measured critical indices. This happens at low values of chi, but there i

Background: The objective of this study was to assess the regional geometry of the Heineke-Mikulicz (HM) strictureplasty. The HM intestinal strictureplasty is commonly performed for the treatment of stricturing Crohn's disease of the small intestine. This procedure shifts relatively normal proximal and distal tissue to the point of narrowing and thus increases the luminal diameter. The overall effect on the regional geometry of the HM strictureplasty, however, has not been previously described in detail. Methods: HM strictureplasties were created in latex tubing and cast with an epoxy resin. The resultant casts of the lumens were then imaged using computed tomography. Using 3-dimensional vascular reconstruction software, the cross-sectional areas were determined and the surface geometry was examined. Results: The HM strictureplasty, while increasing the lumen at the point of the stricture, also results in a counterproductive luminal narrowing proximal and distal to the strictureplasty. Within the model used, cross-sectional area was diminished 25% to 50% below baseline. This effect is enhanced when 2 strictureplasties are placed in close proximity to each other. Conclusions: The HM strictureplasty results in alterations in the regional geometry that may result in a compromise of the lumen proximal and distal to the location of the strictureplasty. When 2 HM strictureplasties are created in close proximity to each other, care should be undertaken to assure that the lumen of the intervening segment is adequate.

Geometric frustration arises whenever the constituents of a physical assembly locally favor an arrangement that cannot be realized globally. Recently, such frustrated assemblies were shown to exhibit filamentation, size limitation, large morphological variations and other exotic response properties. While these unique characteristics can be shown to be a direct outcome of the geometric frustration, some geometrically frustrated systems do not exhibit any of the above phenomena. In this work we exploit the intrinsic approach to provide a framework for directly addressing the frustration in physical assemblies. The framework highlights the role of the compatibility conditions associated with the intrinsic fields describing the physical assembly. We show that the structure of the compatibility conditions determines the behavior of small assemblies and in particular predicts their superextensive energy growth exponent. We illustrate the use of this framework to several well-known frustrated assemblies.

We study equilibrium configurations of thin and elongated non-Euclidean elastic strips with hyperbolic two-dimensional reference metrics ? which are invariant along the strip. In the vanishing thickness limit energy minima are obtained by minimizing the integral of the mean curvature squared among all isometric embeddings of ?. For narrow strips these minima are very close to minimal surfaces regardless of the specific form of the metric. We study the properties of these "almost minimal" surfaces and find a rich range of three-dimensional stable configurations. We provide some explicit solutions as well as a framework for the incorporation of additional forces and constraints.

Non-Euclidean plates are thin elastic bodies having no stress-free configuration, hence exhibiting residual stresses in the absence of external constraints. These bodies are endowed with a three-dimensional reference metric, which may not necessarily be immersible in physical space. Here, based on a recently developed theory for such bodies, we characterize the transition from flat to buckled equilibrium configurations at a critical value of the plate thickness. Depending on the reference metric, the buckling transition may be either continuous or discontinuous. In the infinitely thin plate limit, under the assumption that a limiting configuration exists, we show that the limit is a configuration that minimizes the bending content, among all configurations with zero stretching content (isometric immersions of the midsurface). For small but finite plate thickness, we show the formation of a boundary layer, whose size scales with the square root of the plate thickness and whose shape is determined by a balance between stretching and bending energies.

The connection between a surface's metric and its Gaussian curvature (Gauss theorem) provides the base for a shaping principle of locally growing or shrinking elastic sheets. We constructed thin gel sheets that undergo laterally nonuniform shrinkage. This differential shrinkage prescribes non-Euclidean metrics on the sheets. To minimize their elastic energy, the free sheets form three-dimensional structures that follow the imposed metric. We show how both large-scale buckling and multiscale wrinkling structures appeared, depending on the nature of possible embeddings of the prescribed metrics. We further suggest guidelines for how to generate each type of feature.