Source author record

T. C. Lubensky

T. C. Lubensky appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

23works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

23 published item(s)

preprint2020arXiv

Multi-functional Twisted-Kagome lattices: Tuning by Pruning Mechanical Metamaterials

This article investigates phonons and elastic response in randomly diluted lattices constructed by combining (via the addition of next-nearest bonds) a twisted kagome lattice, with bulk modulus $B=0$ and shear modulus $G>0$, with either a generalized untwisted kagome lattice with $B>0$ and $G>0$ or with a honeycomb lattice with $B>0$ and $G=0$. These lattices exhibit jamming-like critical end-points at which $B$, $G$, or both $B$ and $G$ jump discontinuously from zero while the remaining moduli (if any) begin to grow continuously from zero. Pairs of these jamming points are joined by lines of continuous rigidity percolation transitions at which both $B$ and $G$ begin to grow continuously from zero. The Poisson ratio and $G/B$ can be continuously tuned throughout their physical range via random dilution in a manner analogous to "tuning by pruning" in random jammed lattices. These lattices can be produced with modern techniques, such as 3D printing, for constructing metamaterials.

preprint2016arXiv

Elasticity of randomly diluted honeycomb and diamond lattices with bending forces

We use numerical simulations and an effective-medium theory to study the rigidity percolation transition of the honeycomb and diamond lattices when weak bond-bending forces are included. We use a rotationally invariant bond-bending potential, which, in contrast to the Keating potential, does not involve any stretching. As a result, the bulk modulus does not depend on the bending stiffness $κ$. We obtain scaling functions for the behavior of some elastic moduli in the limits of small $Δ\mathcal{P} = 1 - \mathcal{P}$, and small $δ\mathcal{P} = \mathcal{P} - \mathcal{P}_c$, where $\mathcal{P}$ is an occupation probability of each bond, and $\mathcal{P}_c$ is the critical probability at which rigidity percolation occurs. We find good quantitative agreement between effective-medium theory and simulations for both lattices for $\mathcal{P}$ close to one.

preprint2016arXiv

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 enable tests of a standard understanding of the detailed features of the jamming transition, and show that parts of this argument are invalid.

preprint2015arXiv

Entropic forces stabilize diverse emergent structures in colloidal membranes

The depletion interaction mediated by non-adsorbing polymers promotes condensation and assembly of repulsive colloidal particles into diverse higher-order structures and materials. One example, with particularly rich emergent behaviors, is the formation of two-dimensional colloidal membranes from a suspension of filamentous $\it{fd}$ viruses, which act as rods with effective repulsive interactions, and dextran, which acts as a condensing, depletion-inducing agent. Colloidal membranes exhibit chiral twist even when the constituent virus mixture lacks macroscopic chirality, change from a circular shape to a striking starfish shape upon changing the chirality of constituent rods, and partially coalesce via domain walls through which the viruses twist by $180^\circ$. We formulate an entropically-motivated theory that can quantitatively explain these experimental structures and measurements, both previously published and newly performed, over a wide range of experimental conditions. Our results elucidate how entropy alone, manifested through the viruses as Frank elastic energy and through the depletants as an effective surface tension, drives the formation and behavior of these diverse structures. Our generalizable principles propose the existence of analogous effects in molecular membranes and can be exploited in the design of reconfigurable colloidal structures.

preprint2015arXiv

Mechanical Weyl Modes in Topological Maxwell Lattices

Topological mechanical structures exhibit robust properties protected by topological invariants. In this letter, we study a family of deformed square lattices that display topologically protected zero-energy bulk modes analogous to the massless fermion modes of Weyl semimetals. Our findings apply to sufficiently complex lattices satisfying the Maxwell criterion of equal numbers of constraints and degrees of freedom. We demonstrate that such systems exhibit pairs of oppositely charged Weyl points, corresponding to zero-frequency bulk modes, that can appear at the origin of the Brillouin zone and move away to the zone edge (or return to the origin) where they annihilate. We prove that the existence of these Weyl points leads to a wavenumber-dependent count of topological mechanical states at free surfaces and domain walls.

preprint2014arXiv

Mechanical instability at finite temperature

Many physical systems including lattices near structural phase transitions, glasses, jammed solids, and bio-polymer gels have coordination numbers that place them at the edge of mechanical instability. Their properties are determined by an interplay between soft mechanical modes and thermal fluctuations. In this paper we investigate a simple square-lattice model with a $ϕ^4$ potential between next-nearest-neighbor sites whose quadratic coefficient $κ$ can be tuned from positive negative. We show that its zero-temperature ground state for $κ<0$ is highly degenerate, and we use analytical techniques and simulation to explore its finite temperature properties. We show that a unique rhombic ground state is entropically favored at nonzero temperature at $κ<0$ and that the existence of a subextensive number of "floppy" modes whose frequencies vanish at $κ= 0$ leads to singular contributions to the free energy that render the square-to-rhombic transition first order and lead to power-law behavior of the shear modulus as a function of temperature. We expect our study to provide a general framework for the study of finite-temperature mechanical and phase behavior of other systems with a large number of floppy modes.

preprint2014arXiv

Penrose Tilings as Jammed Solids

Penrose tilings form lattices, exhibiting 5-fold 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.

preprint2013arXiv

Effective Medium Theory of Filamentous Triangular Lattice

We present an effective medium theory that includes bending as well as stretching forces, and we use it to calculate 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_{\textrm{CF}} \simeq 2/3$ of the lattice when fiber bending rigidity vanishes.

preprint2013arXiv

Elasticity of 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 four. If lines are treated as semi-flexible polymers and crossing points are treated as crosslinks, 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 \approx 0.605$ that is the same for all non-zero 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=p_b$, that match simulations with no adjustable parameters near $p=1$. The affine response as $p\rightarrow 1$ is identified with the approach to a state with sample-crossing straight filaments treated as elastic rods.

preprint2013arXiv

Topological Boundary Modes in Isostatic Lattices

Frames, or lattices consisting of mass points connected by rigid bonds or central force springs, are important model constructs that have applications in such diverse fields as structural engineering, architecture, and materials science. The difference between the number of bonds and the number of degrees of freedom of these lattices determines the number of their zero-frequency "floppy modes". When these are balanced, the system is on the verge of mechanical instability and is termed isostatic. It has recently been shown that certain extended isostatic lattices exhibit floppy modes localized at their boundary. These boundary modes are insensitive to local perturbations, and appear to have a topological origin, reminiscent of the protected electronic boundary modes that occur in the quantum Hall effect and in topological insulators. In this paper we establish the connection between the topological mechanical modes and the topological band theory of electronic systems, and we predict the existence of new topological bulk mechanical phases with distinct boundary modes. We introduce model systems in one and two dimensions that exemplify this phenomenon.

preprint2011arXiv

Isostaticity, auxetic response, surface modes, and conformal invariance in twisted kagome lattices

Model lattices consisting of balls connected by central-force springs provide much of our understanding of mechanical response and phonon structure of real materials. Their stability depends critically on their coordination number $z$. $d$-dimensional lattices with $z=2d$ are at the threshold of mechanical stability and are isostatic. Lattices with $z<2d$ exhibit zero-frequency "floppy" modes that provide avenues for lattice collapse. The physics of systems as diverse as architectural structures, network glasses, randomly packed spheres, and biopolymer networks is strongly influenced by a nearby isostatic lattice. We explore elasticity and phonons of a special class of two-dimensional isostatic lattices constructed by distorting the kagome lattice. We show that the phonon structure of these lattices, characterized by vanishing bulk moduli and thus negative Poisson ratios and auxetic elasticity, depends sensitively on boundary conditions and on the nature of the kagome distortions. We construct lattices that under free boundary conditions exhibit surface floppy modes only or a combination of both surface and bulk floppy modes; and we show that bulk floppy modes present under free boundary conditions are also present under periodic boundary conditions but that surface modes are not. In the the long-wavelength limit, the elastic theory of all these lattices is a conformally invariant field theory with holographic properties, and the surface waves are Rayleigh waves. We discuss our results in relation to recent work on jammed systems. Our results highlight the importance of network architecture in determining floppy-mode structure.

preprint2011arXiv

Nonlinear smectic elasticity of a helical state in cholesterics and helimagnets

General symmetry arguments, dating back to de Gennes dictate that at scales longer than the pitch, the low-energy elasticity of a chiral nematic liquid crystal (cholesteric) and of a Dzyaloshinskii-Morya (DM) spiral state in a helimagnet with negligible crystal symmetry fields (e.g., MnSi, FeGe) is identical to that of a smectic liquid crystal, thereby inheriting its rich phenomenology. Starting with a chiral Frank free-energy (exchange and DM interactions of a helimagnet) we present a transparent derivation of the fully nonlinear Goldstone mode elasticity, which involves an analog of the Anderson-Higgs mechanism that locks the spiral orthonormal (director/magnetic moment) frame to the cholesteric (helical) layers. This shows explicitly the reduction of three orientational modes of a cholesteric down to a single phonon Goldstone mode that emerges on scales longer than the pitch. At a harmonic level our result reduces to that derived many years ago by Lubensky and collaborators.

preprint2011arXiv

Simple lattice model for biological gels

We construct a three-dimensional lattice model for biological gels in which straight lines of bonds correspond to filamentous semi-flexible polymers and lattice sites, which are exactly four-fold coordinated, to crosslinks. With only stretching central forces between nearest neighbors, this lattice is sub-isostatic with an extensive number of zero modes; but all of its elastic constants are nonzero, and its elastic response is affine. Removal of bonds with probability $1-p$ leads to a lattice with average coordination number less than four and a distribution of polymer lengths. When bending forces are added, the diluted lattice exhibits a rigidity threshold at $p=p_b<1$ and crossover from bending-dominated nonaffine to stretching-dominated affine response between $p_b$ and $p=1$.

preprint2010arXiv

Centrifugal compression of soft particle packings - theory and experiment

An exact method is developed for computing the height of an elastic medium subjected to centrifugal compression, for arbitrary constitutive relation between stress and strain. Example solutions are obtained for power-law media and for cases where the stress diverges at a critical strain -- for example as required by packings composed of deformable but incompressible particles. Experimental data are presented for the centrifugal compression of thermo-responsive N-isopropylacrylamide (NIPA) microgel beads in water. For small radial acceleration, the results are consistent with Hertzian elasticity, and are analyzed in terms of the Young elastic modulus of the bead material. For large radial acceleration, the sample compression asymptotes to a value corresponding to a space-filling particle volume fraction of unity. Therefore we conclude that the gel beads are incompressible, and deform without deswelling. In addition, we find that the Young elastic modulus of the particulate gel material scales with cross-link density raised to the power 3.3+-0.8, somewhat larger than the Flory expectation.

preprint2010arXiv

Coherent potential approximation of random nearly isostatic kagome lattice

The kagome lattice has coordination number $4$, and it is mechanically isostatic when nearest neighbor ($NN$) sites are connected by central force springs. A lattice of $N$ sites has $O(\sqrt{N})$ zero-frequency floppy modes that convert to finite-frequency anomalous modes when next-nearest-neighbor ($NNN$) springs are added. We use the coherent potential approximation (CPA) to study the mode structure and mechanical properties of the kagome lattice in which $NNN$ springs with spring constant $κ$ are added with probability $\Prob= Δz/4$, where $Δz= z-4$ and $z$ is the average coordination number. The effective medium static $NNN$ spring constant $κ_m$ scales as $\Prob^2$ for $\Prob \ll κ$ and as $\Prob$ for $\Prob \gg κ$, yielding a frequency scale $ω^* \sim Δz$ and a length scale $l^*\sim (Δz)^{-1}$. To a very good approximation at at small nonzero frequency, $κ_m(\Prob,ω)/κ_m(\Prob,0)$ is a scaling function of $ω/ω^*$. The Ioffe-Regel limit beyond which plane-wave states becomes ill-define is reached at a frequency of order $ω^*$.

preprint2010arXiv

Criticality and isostaticity in fiber networks

The rigidity of elastic networks depends sensitively on their internal connectivity and the nature of the interactions between constituents. Particles interacting via central forces undergo a zero-temperature rigidity-percolation transition near the isostatic threshold, where the constraints and internal degrees of freedom are equal in number. Fibrous networks, such as those that form the cellular cytoskeleton, become rigid at a lower threshold due to additional bending constraints. However, the degree to which bending governs network mechanics remains a subject of considerable debate. We study disordered fibrous networks with variable coordination number, both above and below the central-force isostatic point. This point controls a broad crossover from stretching- to bending-dominated elasticity. Strikingly, this crossover exhibits an anomalous power-law dependence of the shear modulus on both stretching and bending rigidities. At the central-force isostatic point---well above the rigidity threshold---we find divergent strain fluctuations together with a divergent correlation length $ξ$, implying a breakdown of continuum elasticity in this simple mechanical system on length scales less than $ξ$.

preprint2010arXiv

Non-affine displacements in flexible polymer networks

The validity of the affine assumption in model flexible polymer networks is explored. To this end, the displacements of fluorescent tracer beads embedded in polyacrylamide gels are quantified by confocal microscopy under shear deformation, and the deviations of these displacements from affine responses are recorded. Non-affinity within the gels is quantified as a function of polymer chain density and cross-link concentration. Observations are in qualitative agreement with current theories of polymer network non-affinity. The measured degree of non-affinity in the polyacrylamide gels suggests the presence of structural inhomogeneities which likely result from heterogeneous reaction kinetics during gel preparation. In addition, the macroscopic elasticity of the polyacrylamide gels is confirmed to behave in accordance with standard models of flexible polymer network elasticity.

preprint2010arXiv

Soft modes and elasticity of nearly isostatic lattices: randomness and dissipation

The square lattice with central-force springs on nearest-neighbor bonds is isostatic. It has a zero mode for each row and column, and it does not support shear. Using the Coherent Potential Approximation (CPA), we study how the random addition, with probability $\mathcal{P}=(z-4)/4$ ($z$ = average number of nearest neighbors), of springs on next-nearest-neighbor ($NNN$) bonds restores rigidity and affects phonon structure. We find that the CPA effective $NNN$ spring constant $\tildeκ_m(ω)$, equivalent to the complex shear modulus $G(ω)$, obeys the scaling relation, $\tildeκ_m(ω) = κ_m h(ω/ω^*)$, at small $\mathcal{P}$, where $κ_m = \tildeκ'_m(0)\sim \mathcal{P}^2$ and $ω^* \sim \mathcal{P}$, implying that elastic response is nonaffine at small $\mathcal{P}$ and that plane-wave states are ill-defined beyond the Ioffe-Regel limit at $ω\approx ω^*$. We identify a divergent length $l^* \sim \mathcal{P}^{-1}$, and we relate these results to jamming.

preprint2008arXiv

Dynamics of Gas-Fluidized Granular Rods

We study a quasi-two-dimensional monolayer of granular rods fluidized by a spatially and temporally homogeneous upflow of air. By tracking the position and orientation of the particles, we characterize the dynamics of the system with sufficient resolution to observe ballistic motion at the shortest time scales. Particle anisotropy gives rise to dynamical anisotropy and superdiffusive dynamics parallel to the rod's long axis, causing the parallel and perpendicular mean squared displacements to become diffusive on different timescales. The distributions of free times and free paths between collisions deviate from exponential behavior, underscoring the non-thermal character of the particle motion. The dynamics show evidence of rotational-translational coupling similar to that of an anisotropic Brownian particle. We model rotational-translation coupling in the single-particle dynamics with a modified Langevin model using non-thermal noise sources. This suggests a phenomenological approach to thinking about collections of self-propelling particles in terms of enhanced memory effects.

preprint2006arXiv

Dynamics of smectic elastomers

We study the low-frequency, long-wavelength dynamics of liquid crystal elastomers, crosslinked in the smectic-$A$ phase, in their smectic-$A$, biaxial smectic and smectic-$C$ phases. Two different yet related formulations are employed. One formulation describes the pure hydrodynamics and does not explicitly involve the Frank director, which relaxes to its local equilibrium value in a non-hydrodynamic time. The other formulation explicitly treats the director and applies beyond the hydrodynamic limit. We compare the low-frequency, long-wavelength dynamics of smectic-$A$ elastomers to that of nematics and show that the two are closely related. For the biaxial smectic and the smectic-$C$ phases, we calculate sound velocities and the mode structure in certain symmetry directions. For the smectic-$C$ elastomers, in addition, we discuss in some detail their possible behavior in rheology experiments.

preprint1999arXiv

Sliding Phases in XY-Models, Crystals, and Cationic Lipid-DNA Complexes

We predict the existence of a totally new class of phases in weakly coupled, three-dimensional stacks of two-dimensional (2D) XY-models. These ``sliding phases'' behave essentially like decoupled, independent 2D XY-models with precisely zero free energy cost associated with rotating spins in one layer relative to those in neighboring layers. As a result, the two-point spin correlation function decays algebraically with in-plane separation. Our results, which contradict past studies because we include higher-gradient couplings between layers, also apply to crystals and may explain recently observed behavior in cationic lipid-DNA complexes.

preprint1998arXiv

Stability of Texture and Shape of Circular Domains of Langmuir Monolayers

Finite domains of a Langmuir monolayer in a phase with tilted molecules can be modeled by a simple elastic free energy of an XY order parameter with isotropic and anisotropic line tension terms. The domains can and often do contain nontrivial textures, which in turn influence the shape of the domains. Herein we investigate the properties of a simplified isotropic model with a single elastic constant. For circular domains a first-order phase transition is found between two distinct textures: an exterior defect (or ``virtual boojum'') texture, and an interior defect texture. Starting with a circular domain and either of these two textures as a ground state, we find shape instabilities develop that depend on the elastic constants and line tensions in the simplified model. In both cases a necessary but not sufficient condition for the onset of shape instabilities is the possibility for a local negative effective line tension to develop from the anisotropic line tension term.