Publications

The geometry and interactions between the constituents of a liquid crystal, which are responsible for inducing the partial order in the fluid, may locally favor an attempted phase that could not be realized in R-3. While states that are incompatible with the geometry of R-3 were identified more than 50 years ago, the collection of compatible states remained poorly understood and not well characterized. Recently, the compatibility conditions for three-dimensional director fields were derived using the method of moving frames. These compatibility conditions take the form of six differential relations in five scalar fields locally characterizing the director field. In this work, we rederive these equations using a more transparent approach employing vector calculus. We then use these equations to characterize a wide collection of compatible phases.

Geometric frustration results from a discrepancy between the locally favored arrangement of the constituents of a system and the geometry of the embedding space. Geometric frustration can be either noncumulative, which implies an extensive energy growth, or cumulative, which implies superextensive energy scaling and highly cooperative ground-state configurations which may depend on the dimensions of the system. Cumulative geometric frustration was identified in a variety of continuous systems including liquid crystals, filament bundles, and molecular crystals. However, a spin-lattice model which clearly demonstrates cumulative geometric frustration was lacking. In this paper we describe a nonlinear variation of the XY-spin model on a triangular lattice that displays cumulative geometric frustration. The model is studied numerically and analyzed in three distinct parameter regimes, which are associated with different energy minimizing configurations. We show that, despite the difference in the ground-state structure in the different regimes, in all cases the superextensive power-law growth of the frustration energy for small domains grows with the same universal exponent that is predicted from the structure of the underlying compatibility condition.

Minimal surfaces arise as energy minimizers for fluid membranes and are thus found in a variety of biological systems. The tight lamellar structures of the endoplasmic reticulum and plant thylakoids are comprised of such minimal surfaces in which right- and left-handed helical motifs are embedded in stoichiometry suggesting global pitch balance. So far, the analytical treatment of helical motifs in minimal surfaces was limited to the small-slope approximation where motifs are represented by the graph of harmonic functions. However, in most biologically and physically relevant regimes the inter-motif separation is comparable with its pitch, and thus this approximation fails. Here, we present a recipe for constructing exact minimal surfaces with an arbitrary distribution of helical motifs, showing that any harmonic graph can be deformed into a minimal surface by exploiting lateral displacements only. We analyse in detail pairs of motifs of the similar and of opposite handedness and also an infinite chain of identical motifs with similar or alternating handedness. Last, we study the second variation of the area functional for collections of helical motifs with asymptotic helicoidal structure and show that in this subclass of minimal surfaces stability requires that the collection of motifs is pitch balanced.

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.

In flat space, changing a system's velocity requires the presence of an external force. However, an isolated nonrigid system can freely change its orientation due to the nonholonomic nature of the angular momentum conservation law. Such nonrigid isolated systems may thus manifest their internal dynamics as rotations. In this work, we show that for such systems chaotic internal dynamics may lead to macroscopic rotational random walk resembling thermally induced motion. We do so by studying the classical harmonic three-mass system in the strongly nonlinear regime, the simplest physical model capable of zero angular momentum rotation as well as chaotic dynamics. At low energies, the dynamics are regular and the system rotates at a constant rate with zero angular momentum. For sufficiently high energies a rotational random walk is observed. For intermediate energies the system performs ballistic bouts of constant rotation rates interrupted by unpredictable orientation reversal events, and the system constitutes a simple physical model for Levy walks. The orientation reversal statistics in this regime lead to a fractional rotational diffusion that interpolates smoothly between the ballistic and regular diffusive regimes.

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.