Giant director fluctuations in liquid crystal drops.

We report the discovery and elucidation of giant spatiotemporal orientational fluctuations in nematic liquid crystal drops with radial orientation of the nematic anisotropy axis producing a central "hedgehog" defect. We study the spatial and temporal properties of the fluctuations experimentally using polarized optical microscopy, and theoretically, by calculating the eigenspectrum of the Frank elastic free energy of a nematic drop of radius R_{2}, containing a spherical central core of radius R_{1} and constrained by perpendicular boundary conditions on all surfaces. We find that the hedgehog defect with radial orientation has a complex excitation spectrum with a single critical mode whose energy vanishes at a critical value µ_{c} of the ratio µ=R_{2}/R_{1}. When µ<µ_{c}, the mode has positive energy, indicating that the radial hedgehog state is stable; when µ>µ_{c}, it has negative energy indicating that the radial state is unstable to the formation of a lower-energy state. This mode gives rise to the large-amplitude director fluctuations we observe near the core, for µ near µ_{c}. A collapse of the experimental data corroborates model predictions for µ<µ_{c} and provides an estimate of the defect core size.

Theory of director fluctuations about a hedgehog defect in a nematic drop.

We present calculations of eigenmode energies and wave functions of both azimuthal and polar distortions of the nematic director relative to a radial hedgehog trapped in a spherical drop with a smaller concentric spherical droplet at its core. All surfaces interior to the drop have perpendicular (homeotropic) boundary conditions. We also calculate director correlation functions and their relaxation times. Of particular interest is a critical mode whose energy, with fixed Frank constants, vanishes as the ratio µ=R_{2}/R_{1} increases toward a critical value µ_{c}, where R_{2} is the radius of the drop and R_{1} that of the inner droplet, and then becomes negative for µ>µ_{c}. Our calculations form a basis for interpreting experimental measurements of director fluctuations relative to a radial hedgehog state in a spherical drop. We compare results with those obtained by previous investigations, which use a calculational approach different from ours, and with our experimental observations.

Signatures of Topological Phonons in Superisostatic Lattices.

Soft topological surface phonons in idealized ball-and-spring lattices with coordination number z=2d in d dimensions become finite-frequency surface phonons in physically realizable superisostatic lattices with z>2d. We study these finite-frequency modes in model lattices with added next-nearest-neighbor springs or bending forces at nodes with an eye to signatures of the topological surface modes that are retained in the physical lattices. Our results apply to metamaterial lattices, prepared with modern printing techniques, that closely approach isostaticity.

Jamming as a Multicritical Point.

The discontinuous jump in the bulk modulus B at the jamming transition is a consequence of the formation of a critical contact network of spheres that resists compression. We introduce lattice models with underlying undercoordinated compression-resistant spring lattices to which next-nearest-neighbor springs can be added. In these models, the jamming transition emerges as a kind of multicritical point terminating a line of rigidity-percolation transitions. Replacing the undercoordinated lattices with the critical network at jamming yields a faithful description of jamming and its relation to rigidity percolation.

Directed percolation with a conserved field and the depinning transition.

Conserved directed percolation (C-DP) and the depinning transition of a disordered elastic interface belong to the same universality class, as has been proven very recently by Le Doussal and Wiese [Phys. Rev. Lett. 114, 110601 (2015)PRLTAO0031-900710.1103/PhysRevLett.114.110601] through a mapping of the field theory for C-DP onto that of the quenched Edwards-Wilkinson model. Here, we present an alternative derivation of the C-DP field theoretic functional, starting with the coherent-state path integral formulation of the C-DP and then applying the Grassberger transformation, which avoids the disadvantages of the so-called Doi shift. We revisit the aforementioned mapping with focus on a specific term in the field theoretic functional that has been problematic in the past when it came to assessing its relevance. We show that this term is redundant in the sense of the renormalization group.

Topological Phonons and Weyl Lines in Three Dimensions.

Topological mechanics and phononics have recently emerged as an exciting field of study. Here we introduce and study generalizations of the three-dimensional pyrochlore lattice that have topologically protected edge states and Weyl lines in their bulk phonon spectra, which lead to zero surface modes that flip from one edge to the opposite as a function of surface wave number.

Topological boundary modes in jammed matter.

Granular matter at the jamming transition is poised on the brink of mechanical stability, and hence it is possible that these random systems have topologically protected surface phonons. Studying two model systems for jammed matter, we find states that exhibit distinct mechanical topological classes, protected surface modes, and ubiquitous Weyl points. The detailed statistics of the boundary modes shed surprising light on the properties of the jamming critical point and help inform a common theoretical description of the detailed features of the transition.

Penrose tilings as jammed solids.

Penrose tilings form lattices, exhibiting fivefold symmetry and isotropic elasticity, with inhomogeneous coordination much like that of the force networks in jammed systems. Under periodic boundary conditions, their average coordination is exactly four. We study the elastic and vibrational properties of rational approximants to these lattices as a function of unit-cell size N(S) and find that they have of order sqrt[N(S)] zero modes and states of self-stress and yet all their elastic moduli vanish. In their generic form, obtained by randomizing site positions, their elastic and vibrational properties are similar to those of particulate systems at jamming with a nonzero bulk modulus, vanishing shear modulus, and a flat density of states.

Driven surface diffusion with detailed balance and elastic phase transitions.

Driven surface diffusion occurs, for example, in molecular beam epitaxy when particles are deposited under an oblique angle. Elastic phase transitions happen when normal modes in crystals become soft due to the vanishing of certain elastic constants. We show that these seemingly entirely disparate systems fall under appropriate conditions into the same universality class. We derive the field-theoretic Hamiltonian for this universality class, and we use renormalized field theory to calculate critical exponents and logarithmic corrections for several experimentally relevant quantities.

Effective-medium theory of a filamentous triangular lattice.

We present an effective-medium theory that includes bending as well as stretching forces, and we use it to calculate the mechanical response of a diluted filamentous triangular lattice. In this lattice, bonds are central-force springs, and there are bending forces between neighboring bonds on the same filament. We investigate the diluted lattice in which each bond is present with a probability p. We find a rigidity threshold p(b) which has the same value for all positive bending rigidity and a crossover characterizing bending, stretching, and bend-stretch coupled elastic regimes controlled by the central-force rigidity percolation point at p(CF)=/~2/3 of the lattice when fiber bending rigidity vanishes.

Elasticity of a filamentous kagome lattice.

The diluted kagome lattice, in which bonds are randomly removed with probability 1-p, consists of straight lines that intersect at points with a maximum coordination number of 4. If lines are treated as semiflexible polymers and crossing points are treated as cross-links, this lattice provides a simple model for two-dimensional filamentous networks. Lattice-based effective-medium theories and numerical simulations for filaments modeled as elastic rods, with stretching modulus µ and bending modulus κ, are used to study the elasticity of this lattice as functions of p and κ. At p=1, elastic response is purely affine, and the macroscopic elastic modulus G is independent of κ. When κ=0, the lattice undergoes a first-order rigidity-percolation transition at p=1. When κ>0, G decreases continuously as p decreases below one, reaching zero at a continuous rigidity-percolation transition at p=p(b)≈0.605 that is the same for all nonzero values of κ. The effective-medium theories predict scaling forms for G, which exhibit crossover from bending-dominated response at small κ/µ to stretching-dominated response at large κ/µ near both p=1 and p(b), that match simulations with no adjustable parameters near p=1. The affine response as p→1 is identified with the approach to a state with sample-crossing straight filaments treated as elastic rods.

Scaling exponents for a monkey on a tree: fractal dimensions of randomly branched polymers.

We study asymptotic properties of diffusion and other transport processes (including self-avoiding walks and electrical conduction) on large, randomly branched polymers using renormalized dynamical field theory. We focus on the swollen phase and the collapse transition, where loops in the polymers are irrelevant. Here the asymptotic statistics of the polymers is that of lattice trees, and diffusion on them is reminiscent of the climbing of a monkey on a tree. We calculate a set of universal scaling exponents including the diffusion exponent and the fractal dimension of the minimal path to two-loop order and, where available, compare them to numerical results.

Linear polymers in disordered media: the shortest, the longest, and the mean self-avoiding walk on percolation clusters.

Long linear polymers in strongly disordered media are well described by self-avoiding walks (SAWs) on percolation clusters and a lot can be learned about the statistics of these polymers by studying the length distribution of SAWs on percolation clusters. This distribution encompasses 2 distinct averages, viz., the average over the conformations of the underlying cluster and the SAW conformations. For the latter average, there are two basic options, one being static and one being kinetic. It is well known for static averaging that if the disorder of the underlying medium is weak, this disorder is redundant in the sense the renormalization group; i.e., differences to the ordered case appear merely in nonuniversal quantities. Using dynamical field theory, we show that the same holds true for kinetic averaging. Our main focus, however, lies on strong disorder, i.e., the medium being close to the percolation point, where disorder is relevant. Employing a field theory for the nonlinear random resistor network in conjunction with a real-world interpretation of the corresponding Feynman diagrams, we calculate the scaling exponents for the shortest, the longest, and the mean or average SAW to 2-loop order. In addition, we calculate to 2-loop order the entire family of multifractal exponents that governs the moments of the the statistical weights of the elementary constituents (bonds or sites of the underlying fractal cluster) contributing to the SAWs. Our RG analysis reveals that kinetic averaging leads to renormalizability whereas static averaging does not, and hence, we argue that the latter does not lead to a well-defined scaling limit. We discuss the possible implications of this finding for experiments and numerical simulations which have produced widespread results for the exponent of the average SAW. To corroborate our results, we also study the well-known Meir-Harris model for SAWs on percolation clusters. We demonstrate that the Meir-Harris model leads back up to 2-loop order to the renormalizable real-world formulation with kinetic averaging if the replica limit is consistently performed at the first possible instant in the course of the calculation.

Collapse transition of randomly branched polymers: renormalized field theory.

We present a minimal dynamical model for randomly branched isotropic polymers, and we study this model in the framework of renormalized field theory. For the swollen phase, we show that our model provides a route to understand the well-established dimensional-reduction results from a different angle. For the collapse θ transition, we uncover a hidden Becchi-Rouet-Stora supersymmetry, signaling the sole relevance of tree configurations. We correct the long-standing one-loop results for the critical exponents, and we push these results on to two-loop order. For the collapse θ' transition, we find a runaway of the renormalization group flow, which lends credence to the possibility that this transition is a fluctuation-induced first-order transition. Our dynamical model allows us to calculate for the first time the fractal dimension of the shortest path on randomly branched polymers in the swollen phase as well as at the collapse transition and related fractal dimensions.

Renormalized field theory of collapsing directed randomly branched polymers.

We present a dynamical field theory for directed randomly branched polymers and in particular their collapse transition. We develop a phenomenological model in the form of a stochastic response functional that allows us to address several interesting problems such as the scaling behavior of the swollen phase and the collapse transition. For the swollen phase, we find that by choosing model parameters appropriately, our stochastic functional reduces to the one describing the relaxation dynamics near the Yang-Lee singularity edge. This corroborates that the scaling behavior of swollen branched polymers is governed by the Yang-Lee universality class as has been known for a long time. The main focus of our paper lies on the collapse transition of directed branched polymers. We show to arbitrary order in renormalized perturbation theory with epsilon expansion that this transition belongs to the same universality class as directed percolation.

Distribution functions in percolation problems.

Percolation clusters are random fractals whose geometrical and transport properties can be characterized with the help of probability distribution functions. Using renormalized field theory, we determine the asymptotic form of various such distribution functions in the limits where certain scaling variables become small or large. Our study includes the pair-connection probability, the distributions of the fractal masses of the backbone, the red bonds, and the shortest, the longest, and the average self-avoiding walk between any two points on a cluster, as well as the distribution of the total resistance in the random resistor network. Our analysis draws solely on general, structural features of the underlying diagrammatic perturbation theory, and hence our main results are valid to arbitrary loop order.

Smectic-C tilt under shear in smectic-A elastomers.

Stenull and Lubensky [Phys. Rev. E 76, 011706 (2007)] have argued that shear strain and tilt of the director relative to the layer normal are coupled in smectic elastomers and that the imposition of one necessarily leads to the development of the other. This means, in particular, that a smectic-A elastomer subjected to a simple shear will develop smectic-C-like tilt of the director. Recently, Kramer and Finkelmann [e-print arXiv:0708.2024; Phys. Rev. E 78, 021704 (2008)], performed shear experiments on smectic-A elastomers using two different shear geometries. One of the experiments, which implements simple shear, produces clear evidence for the development of smectic-C-like tilt. Here, we generalize a model for smectic elastomers introduced by Adams and Warner [Phys. Rev. E 71, 021708 (2005)] and use it to study the magnitude of SmC-like tilt under shear for the two geometries investigated by Kramer and Finkelmann. Using reasonable estimates of model parameters, we estimate the tilt angle for both geometries, and we compare our estimates to the experimental results. The other shear geometry is problematic since it introduces additional in-plane compressions in a sheetlike sample, thus inducing instabilities that we discuss.

Anomalous elasticity in nematic and smectic elastomer tubule phases.

We study anomalous elasticity in the tubule phases of nematic and smectic elastomer membranes, which are flat in one direction and crumpled in another. These phases share the same macroscopic symmetry properties including spontaneously broken in-plane isotropy and hence belong to the same universality class. Below an upper critical value D_{c}=3 of the membranes' intrinsic dimension D , thermal fluctuations renormalize the elasticity with respect to elastic displacements along the tubule axis so that elastic moduli for compression along the tubule axis and for bending the tubule axis become length-scale dependent. This anomalous elasticity belongs to the same universality class as that of d -dimensional conventional smectic liquid crystals with D taking on the role of d . For physical tubule phases, D=2 , this anomaly is of power-law type and thus might by easier to detect experimentally than the logarithmic anomaly in conventional smectics.

Smectic-A elastomers with weak director anchoring.

Experimentally it is possible to manipulate the director in a (chiral) smectic- A elastomer using an electric field. This suggests that the director is not necessarily locked to the layer normal, as described in earlier papers that extended rubber elasticity theory to smectics. Here, we consider the case that the director is weakly anchored to the layer normal assuming that there is a free energy penalty associated with relative tilt between the two. We use a recently developed weak-anchoring generalization of rubber elastic approaches to smectic elastomers and study shearing in the plane of the layers, stretching in the plane of the layers, and compression and elongation parallel to the layer normal. We calculate, inter alia, the engineering stress and the tilt angle between director and layer normal as functions of the applied deformation. For the latter three deformations, our results predict the existence of an instability towards the development of shear accompanied by smectic- C-like order.