Source author record

Meera Sitharam

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

Computational Geometry Biological Physics Machine Learning math.CO physics.chem-ph Artificial Intelligence Biomolecules Computer Vision cond-mat.soft Discrete Mathematics Graphics math.MG math.PR

Catalog footprint

Research overview

17works

13topics

4close collaborators

Actions

Is this your work?

Open graph Browse papers

preprint2020arXiv

Bounds on the Jensen Gap, and Implications for Mean-Concentrated Distributions

This paper gives upper and lower bounds on the gap in Jensen's inequality, i.e., the difference between the expected value of a function of a random variable and the value of the function at the expected value of the random variable. The bounds depend only on growth properties of the function and specific moments of the random variable. The bounds are particularly useful for distributions that are concentrated around the mean, a commonly occurring scenario such as the average of i.i.d. samples and in statistical mechanics.

preprint2020arXiv

Fast and Flexible Geometric Method For Enhancing MC Sampling of Compact Configurations For Protein Docking Problem

EASAL (Efficient Atlasing and Sampling of Assembly Landscapes) is a geometric method for sampling and computing integrals over the potential energy landscape of small molecular assemblies. EASAL's efficiency arises from the fact that small assembly landscapes permit the use of so-called Cayley (inter-atomic distance based) parameters for geometric representation and sampling of the assembly configuration space regions; this results in their isolation, convexification, customized sampling and systematic traversal using a comprehensive topological roadmap. We define custom-designed measurements to investigate and compare various sampling characteristics of EASAL and the traditional Monte Carlo (MC) sampling, including (i) sampling speed, (ii) efficiency and accuracy of uniform grid coverage, (iii) accuracy of weighted coverage at covering low energy regions, (iv) ability to localize sampling to macrostates, and (v) flexibility in sampling distributions. In particular, we compare the sampling characteristics of EASAL and MC in sampling the assembly landscape of 2 trans-membrane helices, with short-range pair-potentials. We demonstrate that EASAL provides a reasonable coverage of crucial but narrow regions of the energy landscape of low effective dimension, with much fewer samples and computational resources than MC sampling. Promising avenues for combining the complementary advantages of the two methods are discussed.

preprint2020arXiv

Plane-Activated Mapped Microstructure

Querying and interacting with models of massive material micro-structure requires localized on-demand generation of the micro-structure since the full-scale storing and retrieving is cost prohibitive. When the micro-structure is efficiently represented as the image of a canonical structure under a non-linear space deformation to allow it to conform to curved shape, the additional challenge is to relate the query of the mapped micro-structure back to its canonical structure. This paper presents an efficient algorithm to pull back a mapped micro-structure to a partition of the canonical domain structure into boxes and only activates boxes whose image is likely intersected by a plane. The active boxes are organized into a forest whose trees are traversed depth first to generate mapped micro-structure only of the active boxes. The traversal supports, for example, 3D print slice generation in additive manufacturing.

preprint2020arXiv

Rapid prediction of crucial hotspot interactions for icosahedral viral capsid self-assembly by energy landscape atlasing validated by mutagenesis

Icosahedral viruses have their infectious genome encapsulated by a shell assembled by a multiscale process, starting from an integer multiple of 60 viral capsid or coat protein (VP) monomers. We predict and validate inter-atomic hotspot interactions between VP monomers that are important for the assembly of 3 icosahedral viral capsids: Adeno Associated Virus serotype 2 (AAV2) and Minute Virus of Mice (MVM), both T=1 single stranded DNA viruses, and Bromo Mosaic Virus (BMV), a T=3 single stranded RNA virus. Experimental validation is by in-vitro, site-directed mutagenesis data found in literature. We combine ab-initio predictions at two scales: at the interface-scale, we predict the importance (cruciality) of an interaction for successful subassembly across each interface between VP monomers; and at the capsid-scale, we predict the cruciality of an interface for successful capsid assembly. At the interface-scale, we measure cruciality by changes in the capsid free-energy landscape partition function when an interaction is removed. The partition function computation uses atlases of interface subassembly landscapes, rapidly generated by a novel geometric method and curated opensource software EASAL (efficient atlasing and search of assembly landscapes). At the capsid-scale, cruciality of an interface for successful assembly of the capsid is based on combinatorial entropy. Our study goes from resource-light, multiscale computational predictions of crucial hotspot inter-atomic interactions to validation using data on site-directed mutagenesis' effect on capsid assembly. By reliably and rapidly narrowing down target interactions, (no more than 1.5 hours per interface on a laptop with Intel Core i5-2500K 3.2Ghz CPU and 8GB of RAM) our predictions can inform and reduce time-consuming in-vitro and in-vivo experiments, or more computationally intensive in-silico analyses.

preprint2018arXiv

Corner-Sharing Tetrahedra for Modeling Micro-Structure

State-of-the-art representations of volumetric multi-scale shape and structure can be classified into three broad categories: continuous, continuous-from-discrete, and discrete representations. We propose modeling micro-structure with a class of discrete Corner-Sharing Tetrahedra (CoSTs). CoSTs can represent bar-joint, tensegrity, line-incidence, and similar constraint systems that capture local physical constraints and global multi-scale properties for design and analysis. The paper develops a palette of simple geometry processing operations on CoSTs including graph manipulation, hierarchical refinement, randomization, and generating associated continuous representations.

preprint2016arXiv

Combinatorial rigidity of Incidence systems and Application to Dictionary learning

Given a hypergraph $H$ with $m$ hyperedges and a set $Q$ of $m$ \emph{pinning subspaces}, i.e.\ globally fixed subspaces in Euclidean space $\mathbb{R}^d$ , a \emph{pinned subspace-incidence system} is the pair $(H, Q)$ , with the constraint that each pinning subspace in $Q$ is contained in the subspace spanned by the point realizations in $\mathbb{R}^d$ of vertices of the corresponding hyperedge of $H$ . This paper provides a combinatorial characterization of pinned subspace-incidence systems that are \emph{minimally rigid}, i.e.\ those systems that are guaranteed to generically yield a locally unique realization. Pinned subspace-incidence systems have applications in the \emph{Dictionary Learning (aka sparse coding)} problem, i.e.\ the problem of obtaining a sparse representation of a given set of data vectors by learning \emph{dictionary vectors} upon which the data vectors can be written as sparse linear combinations. Viewing the dictionary vectors from a geometry perspective as the spanning set of a subspace arrangement, the result gives a tight bound on the number of dictionary vectors for sufficiently randomly chosen data vectors, and gives a way of constructing a dictionary that meets the bound. For less stringent restrictions on data, but a natural modification of the dictionary learning problem, a further dictionary learning algorithm is provided. Although there are recent rigidity based approaches for low rank matrix completion, we are unaware of prior application of combinatorial rigidity techniques in the setting of Dictionary Learning. We also provide a systematic classification of problems related to dictionary learning together with various algorithms, their assumptions and performance.

preprint2016arXiv

Symmetry in Sphere-based Assembly Configuration Spaces

Many remarkably robust, rapid and spontaneous self-assembly phenomena in nature can be modeled geometrically starting from a collection of rigid bunches of spheres. This paper highlights the role of symmetry in sphere-based assembly processes. Since spheres within bunches could be identical and bunches could be identical as well, the underlying symmetry groups could be of large order that grows with the number of participating spheres and bunches. Thus, understanding symmetries and associated isomorphism classes of microstates correspond to various types of macrostates can significantly reduce the complexity of computing entropy and free energy, as well as paths and kinetics, in high dimensional configuration spaces. In addition, a precise understanding of symmetries is crucial for giving provable guarantees of algorithmic accuracy and efficiency in such computations. In particular, this may aid in predicting crucial assembly-driving interactions. This is a primarily expository paper that develops a novel, original framework for dealing with symmetries in configuration spaces of assembling spheres with the following goals. (1) We give new, formal definitions of various concepts relevant to sphere-based assembly that occur in previous work, and in turn, formal definitions of their relevant symmetry groups leading to the main theorem concerning their symmetries. These previously developed concepts include, for example, (a) assembly configuration spaces, (b) stratification of assembly configuration space into regions defined by active constraint graphs, (c) paths through the configurational regions, and (d) coarse assembly pathways. (2) We demonstrate the new symmetry concepts to compute sizes and numbers of orbits in two example settings appearing in previous work. (3) We give formal statements of a variety of open problems and challenges using the new conceptual definitions.

preprint2015arXiv

An Incidence Geometry approach to Dictionary Learning

We study the Dictionary Learning (aka Sparse Coding) problem of obtaining a sparse representation of data points, by learning \emph{dictionary vectors} upon which the data points can be written as sparse linear combinations. We view this problem from a geometry perspective as the spanning set of a subspace arrangement, and focus on understanding the case when the underlying hypergraph of the subspace arrangement is specified. For this Fitted Dictionary Learning problem, we completely characterize the combinatorics of the associated subspace arrangements (i.e.\ their underlying hypergraphs). Specifically, a combinatorial rigidity-type theorem is proven for a type of geometric incidence system. The theorem characterizes the hypergraphs of subspace arrangements that generically yield (a) at least one dictionary (b) a locally unique dictionary (i.e.\ at most a finite number of isolated dictionaries) of the specified size. We are unaware of prior application of combinatorial rigidity techniques in the setting of Dictionary Learning, or even in machine learning. We also provide a systematic classification of problems related to Dictionary Learning together with various algorithms, their assumptions and performance.

preprint2015arXiv

Combinatorial rigidity and independence of generalized pinned subspace-incidence constraint systems

Given a hypergraph $H$ with $m$ hyperedges and a set $X$ of $m$ \emph{pins}, i.e.\ globally fixed subspaces in Euclidean space $\mathbb{R}^d$ , a \emph{pinned subspace-incidence system} is the pair $(H, X)$ , with the constraint that each pin in $X$ lies on the subspace spanned by the point realizations in $\mathbb{R}^d$ of vertices of the corresponding hyperedge of $H$ . We are interested in combinatorial characterization of pinned subspace-incidence systems that are \emph{minimally rigid}, i.e.\ those systems that are guaranteed to generically yield a locally unique realization. As is customary, this is accompanied by a characterization of generic independence as well as rigidity. In a previous paper \cite{sitharam2014incidence}, we used pinned subspace-incidence systems towards solving the \emph{fitted dictionary learning} problem, i.e.\ dictionary learning with specified underlying hypergraph, and gave a combinatorial characterization of minimal rigidity for a more restricted version of pinned subspace-incidence system, with $H$ being a uniform hypergraph and pins in $X$ being 1-dimension subspaces. Moreover in a recent paper \cite{Baker2015}, the special case of pinned line incidence systems was used to model biomaterials such as cellulose and collagen fibrils in cell walls. In this paper, we extend the combinatorial characterization to general pinned subspace-incidence systems, with $H$ being a non-uniform hypergraph and pins in $X$ being subspaces with arbitrary dimension. As there are generally many data points per subspace in a dictionary learning problem, which can only be modeled with pins of dimension larger than $1$ , such an extension enables application to a much larger class of fitted dictionary learning problems.

preprint2015arXiv

On Flattenability of Graphs

We consider a generalization of the concept of $d$ -flattenability of graphs - introduced for the $l_2$ norm by Belk and Connelly - to general $l_p$ norms, with integer $P$ , $1 \le p < \infty$ , though many of our results work for $l_\infty$ as well. The following results are shown for graphs $G$ , using notions of genericity, rigidity, and generic $d$ -dimensional rigidity matroid introduced by Kitson for frameworks in general $l_p$ norms, as well as the cones of vectors of pairwise $l_p^p$ distances of a finite point configuration in $d$ -dimensional, $l_p$ space: (i) $d$ -flattenability of a graph $G$ is equivalent to the convexity of $d$ -dimensional, inherent Cayley configurations spaces for $G$ , a concept introduced by the first author; (ii) $d$ -flattenability and convexity of Cayley configuration spaces over specified non-edges of a $d$ -dimensional framework are not generic properties of frameworks (in arbitrary dimension); (iii) $d$ -flattenability of $G$ is equivalent to all of $G$ 's generic frameworks being $d$ -flattenable; (iv) existence of one generic $d$ -flattenable framework for $G$ is equivalent to the independence of the edges of $G$ , a generic property of frameworks; (v) the rank of $G$ equals the dimension of the projection of the $d$ -dimensional stratum of the $l_p^p$ distance cone. We give stronger results for specific norms for $d = 2$ : we show that (vi) 2-flattenable graphs for the $l_1$ -norm (and $l_\infty$ -norm) are a larger class than 2-flattenable graphs for Euclidean $l_2$ -norm case and finally (vii) prove further results towards characterizing 2-flattenability in the $l_1$ -norm. A number of conjectures and open problems are posed.

preprint2015arXiv

Optimal Decomposition and Recombination of Isostatic Geometric Constraint Systems for Designing Layered Materials

Optimal recursive decomposition (or DR-planning) is crucial for analyzing, designing, solving or finding realizations of geometric constraint sytems. While the optimal DR-planning problem is NP-hard even for general 2D bar-joint constraint systems, we describe an O(n^3) algorithm for a broad class of constraint systems that are isostatic or underconstrained. The algorithm achieves optimality by using the new notion of a canonical DR-plan that also meets various desirable, previously studied criteria. In addition, we leverage recent results on Cayley configuration spaces to show that the indecomposable systems---that are solved at the nodes of the optimal DR-plan by recombining solutions to child systems---can be minimally modified to become decomposable and have a small DR-plan, leading to efficient realization algorithms. We show formal connections to well-known problems such as completion of underconstrained systems. Well suited to these methods are classes of constraint systems that can be used to efficiently model, design and analyze quasi-uniform (aperiodic) and self-similar, layered material structures. We formally illustrate by modeling silica bilayers as body-hyperpin systems and cross-linking microfibrils as pinned line-incidence systems. A software implementation of our algorithms and videos demonstrating the software are publicly available online (visit http://cise.ufl.edu/~tbaker/drp/index.html.)

preprint2014arXiv

Best of Both Worlds: Uniform sampling in Cartesian and Cayley Molecular Assembly Configuration Space

EASAL (efficient atlasing and sampling of assembly landscapes) is a recently reported geometric method for representing, visualizing, sampling and computing integrals over the potential energy landscape tailored for small molecular assemblies. EASAL's efficiency arises from the fact that small assembly landscapes permit the use of so-called Cayley parameters (inter-atomic distances) for geometric representation and sampling of the assembly configuration space regions; this results in their isolation, convexification, customized sampling and systematic traversal using a comprehensive topological roadmap, ensuring reasonable coverage of crucial but narrow regions of low effective dimension. However, this alone is inadequate for accurate computation of configurational entropy and other integrals, required for estimation of both free energy and kinetics - where it is essential to obtain uniform sampling in appropriate cartesian or moduli space parameterization. Standard adjustment of Cayley sampling via the Jacobian of the map between the two parameterizations is fraught with challenges stemming from an illconditioned Jacobian. This paper formalizes and analyzes these challenges to provide modifications to EASAL that secure the advantages of Cayley sampling while ensuring certain minimum distance and coverage relationships between sampled configurations - in Cartesian space. The modified EASAL's performance is compared with the basic EASAL and the data are presented for Human and Rat Islet Amylin Polypeptide (HiAPP, PDB-2KJ7 and RiAPP PDB-2KB8) dimerization (the two differ in only 6 out of 37 residues, but the former aggregates into fibrils, while the latter does not).

preprint2014arXiv

Cayley Analysis of Mechanism Configuration Spaces using CayMos: Software Functionalities and Architecture

For a common class of 2D mechanisms called 1-dof tree decomposable linkages, we present a software CayMos which uses new theoretical results to implement efficient algorithmic solutions for: (a) meaningfully representing and visualizing the connected components in the Euclidean realization space; (b) finding a path of continuous motion between two realizations in the same connected component, with or without restricting the realization type (sometimes called orientation type); (c) finding two ``closest'' realizations in different connected components.

preprint2013arXiv

Maxwell-independence: a new rank estimate for 3D rigidity matroids

The problem of combinatorially determining the rank of the 3-dimensional bar-joint {\em rigidity matroid} of a graph is an important open problem in combinatorial rigidity theory. Maxwell's condition states that the edges of a graph $G=(V, E)$ are {\em independent} in its $d$ -dimensional generic rigidity matroid only if $(a)$ the number of edges $|E|$ $\le$ $d|V| - {d+1\choose 2}$ , and $(b)$ this holds for every induced subgraph with at least $d$ vertices. We call such graphs {\em Maxwell-independent} in $d$ dimensions. Laman's theorem shows that the converse holds for $d=2$ and thus every maximal Maxwell-independent set of $G$ has size equal to the rank of the 2-dimensional generic rigidity matroid. While this is false for $d=3$ , we show that every maximal, Maxwell-independent set of a graph $G$ has size at least the rank of the 3-dimensional generic rigidity matroid of $G$ . This answers a question posed by Tibór Jordán at the 2008 rigidity workshop at BIRS \cite{bib:birs}. Along the way, we construct subgraphs (1) that yield alternative formulae for a rank upper bound for Maxwell-independent graphs and (2) that contain a maximal (true) independent set. We extend this bound to special classes of non-Maxwell-independent graphs. One further consequence is a simpler proof of correctness for existing algorithms that give rank bounds.

preprint2013arXiv

Nucleation-free $3D$ rigidity

When all non-edge distances of a graph realized in $\mathbb{R}^{d}$ as a {\em bar-and-joint framework} are generically {\em implied} by the bar (edge) lengths, the graph is said to be {\em rigid} in $\mathbb{R}^{d}$ . For $d=3$ , characterizing rigid graphs, determining implied non-edges and {\em dependent} edge sets remains an elusive, long-standing open problem. One obstacle is to determine when implied non-edges can exist without non-trivial rigid induced subgraphs, i.e., {\em nucleations}, and how to deal with them. In this paper, we give general inductive construction schemes and proof techniques to generate {\em nucleation-free graphs} (i.e., graphs without any nucleation) with implied non-edges. As a consequence, we obtain (a) dependent graphs in $3D$ that have no nucleation; and (b) $3D$ nucleation-free {\em rigidity circuits}, i.e., minimally dependent edge sets in $d=3$ . It additionally follows that true rigidity is strictly stronger than a tractable approximation to rigidity given by Sitharam and Zhou \cite{sitharam:zhou:tractableADG:2004}, based on an inductive combinatorial characterization. As an independently interesting byproduct, we obtain a new inductive construction for independent graphs in $3D$ . Currently, very few such inductive constructions are known, in contrast to $2D$ .

preprint2012arXiv

Cayley Configuration Spaces of 1-dof Tree-decomposable Linkages, Part II: Combinatorial Characterization of Complexity

We continue to study Cayley configuration spaces of 1-dof linkages in 2D begun in Part I of this paper, i.e. the set of attainable lengths for a non-edge. In Part II, we focus on the algebraic complexity of describing endpoints of the intervals in the set, i.e., the Cayley complexity. Specifically, We focus on Cayley configuration spaces of a natural class of 1-dof linkages, called 1-dof tree-decomposable linkages. The underlying graphs G satisfy the following: for some base non-edge f, G \cup f is quadratic-radically solvable (QRS), meaning that G \cup f is minimally rigid, and given lengths \bar{l} of all edges, the corresponding linkage (G \cup f, \bar{l}) can be simply realized by ruler and compass starting from f. It is clear that the Cayley complexity only depends on the graph G and possibly the non-edge f. Here we ask whether the Cayley complexity depends on the choice of a base non-edge f. We answer this question in the negative, thereby showing that low Cayley complexity is a property of the graph G (independent of the non-edge f). Then, we give a simple characterization of graphs with low Cayley complexity, leading to an efficient algorithmic characterization, i.e. an efficient algorithm for recognizing such graphs. Next, we show a surprising result that (graph) planarity is equivalent to low Cayley complexity for a natural subclass of 1-dof triangle-decomposable graphs. While this is a finite forbidden minor graph characterization of low Cayley complexity, we provide counterexamples showing impossibility of such finite forbidden minor characterizations when the above subclass is enlarged.

preprint2011arXiv

Learning Hierarchical Sparse Representations using Iterative Dictionary Learning and Dimension Reduction

This paper introduces an elemental building block which combines Dictionary Learning and Dimension Reduction (DRDL). We show how this foundational element can be used to iteratively construct a Hierarchical Sparse Representation (HSR) of a sensory stream. We compare our approach to existing models showing the generality of our simple prescription. We then perform preliminary experiments using this framework, illustrating with the example of an object recognition task using standard datasets. This work introduces the very first steps towards an integrated framework for designing and analyzing various computational tasks from learning to attention to action. The ultimate goal is building a mathematically rigorous, integrated theory of intelligence.

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.

Meera Sitharam

Research overview

Is this your work?

17 publications

Bounds on the Jensen Gap, and Implications for Mean-Concentrated Distributions

Fast and Flexible Geometric Method For Enhancing MC Sampling of Compact Configurations For Protein Docking Problem

Plane-Activated Mapped Microstructure

Rapid prediction of crucial hotspot interactions for icosahedral viral capsid self-assembly by energy landscape atlasing validated by mutagenesis

Corner-Sharing Tetrahedra for Modeling Micro-Structure

Combinatorial rigidity of Incidence systems and Application to Dictionary learning

Symmetry in Sphere-based Assembly Configuration Spaces

An Incidence Geometry approach to Dictionary Learning

Combinatorial rigidity and independence of generalized pinned subspace-incidence constraint systems

On Flattenability of Graphs

Optimal Decomposition and Recombination of Isostatic Geometric Constraint Systems for Designing Layered Materials

Best of Both Worlds: Uniform sampling in Cartesian and Cayley Molecular Assembly Configuration Space

Cayley Analysis of Mechanism Configuration Spaces using CayMos: Software Functionalities and Architecture

Maxwell-independence: a new rank estimate for 3D rigidity matroids

Nucleation-free $3D$ rigidity

Cayley Configuration Spaces of 1-dof Tree-decomposable Linkages, Part II: Combinatorial Characterization of Complexity

Learning Hierarchical Sparse Representations using Iterative Dictionary Learning and Dimension Reduction

See the researcher in context

Building this map preview

Meera Sitharam

Research overview

Is this your work?

17 publications

Bounds on the Jensen Gap, and Implications for Mean-Concentrated Distributions

Fast and Flexible Geometric Method For Enhancing MC Sampling of Compact Configurations For Protein Docking Problem

Plane-Activated Mapped Microstructure

Rapid prediction of crucial hotspot interactions for icosahedral viral capsid self-assembly by energy landscape atlasing validated by mutagenesis

Corner-Sharing Tetrahedra for Modeling Micro-Structure

Combinatorial rigidity of Incidence systems and Application to Dictionary learning

Symmetry in Sphere-based Assembly Configuration Spaces

An Incidence Geometry approach to Dictionary Learning

Combinatorial rigidity and independence of generalized pinned subspace-incidence constraint systems

On Flattenability of Graphs

Optimal Decomposition and Recombination of Isostatic Geometric Constraint Systems for Designing Layered Materials

Best of Both Worlds: Uniform sampling in Cartesian and Cayley Molecular Assembly Configuration Space

Cayley Analysis of Mechanism Configuration Spaces using CayMos: Software Functionalities and Architecture

Maxwell-independence: a new rank estimate for 3D rigidity matroids

Nucleation-free 3D3D3D rigidity

Cayley Configuration Spaces of 1-dof Tree-decomposable Linkages, Part II: Combinatorial Characterization of Complexity

Learning Hierarchical Sparse Representations using Iterative Dictionary Learning and Dimension Reduction

See the researcher in context

Building this map preview

Nucleation-free $3D$ rigidity