-
On odd covers of cliques and disjoint unions
Authors:
Calum Buchanan,
Alexander Clifton,
Eric Culver,
Péter Frankl,
Jiaxi Nie,
Kenta Ozeki,
Puck Rombach,
Mei Yin
Abstract:
Babai and Frankl posed the ``odd cover problem" of finding the minimum cardinality of a collection of complete bipartite graphs such that every edge of the complete graph of order $n$ is covered an odd number of times. In a previous paper with O'Neill, some of the authors proved that this value is always $\lceil n / 2 \rceil$ or $\lceil n / 2 \rceil + 1$ and that it is the former whenever $n$ is a…
▽ More
Babai and Frankl posed the ``odd cover problem" of finding the minimum cardinality of a collection of complete bipartite graphs such that every edge of the complete graph of order $n$ is covered an odd number of times. In a previous paper with O'Neill, some of the authors proved that this value is always $\lceil n / 2 \rceil$ or $\lceil n / 2 \rceil + 1$ and that it is the former whenever $n$ is a multiple of $8$. In this paper, we determine this value to be $\lceil n / 2 \rceil$ whenever $n$ is odd or equivalent to $18$ modulo $24$. We also further the study of odd covers of graphs which are not complete, wherein edges are covered an odd number of times and nonedges an even number of times by the complete bipartite graphs in the collection. Among various results on disjoint unions, we find the minimum cardinality of an odd cover of a union of odd cliques and of a union of cycles.
△ Less
Submitted 16 August, 2024;
originally announced August 2024.
-
Monomial identities in the Weyl algebra
Authors:
Darij Grinberg,
Tom Roby,
Stephan Wagner,
Mei Yin
Abstract:
Motivated by a question and some enumerative conjectures of Richard Stanley, we explore the equivalence classes of words in the Weyl algebra, $\mathbf{k} \langle D,U\rangle/(DU-UD=1)$. We show that each class is generated by the swapping of adjacent *balanced subwords*, i.e., those which have the same number of $D$'s as $U$'s, and give several other characterizations.
Armed with this we deduce a…
▽ More
Motivated by a question and some enumerative conjectures of Richard Stanley, we explore the equivalence classes of words in the Weyl algebra, $\mathbf{k} \langle D,U\rangle/(DU-UD=1)$. We show that each class is generated by the swapping of adjacent *balanced subwords*, i.e., those which have the same number of $D$'s as $U$'s, and give several other characterizations.
Armed with this we deduce a number of enumerative results about the number of such equivalence classes and their sizes. We extend these results to the class of $c$-Dyck words, where every prefix has at least $c$ times as many $U$'s as $D$'s. We also connect these results to previous work on bond percolation and rook theory, and generalize them to some other algebras.
△ Less
Submitted 30 May, 2024;
originally announced May 2024.
-
Fixed points and cycles of parking functions
Authors:
Martin Rubey,
Mei Yin
Abstract:
A parking function of length $n$ is a sequence $π=(π_1,\dots, π_n)$ of positive integers such that if $λ_1\leq\cdots\leq λ_n$ is the increasing rearrangement of $π_1,\dots,π_n$, then $λ_i\leq i$ for $1\leq i\leq n$. In this paper we obtain some exact results on the number of fixed points and cycles of parking functions. Our proofs will be based on generalizations of Pollak's argument. Extensions o…
▽ More
A parking function of length $n$ is a sequence $π=(π_1,\dots, π_n)$ of positive integers such that if $λ_1\leq\cdots\leq λ_n$ is the increasing rearrangement of $π_1,\dots,π_n$, then $λ_i\leq i$ for $1\leq i\leq n$. In this paper we obtain some exact results on the number of fixed points and cycles of parking functions. Our proofs will be based on generalizations of Pollak's argument. Extensions of our techniques are discussed.
△ Less
Submitted 25 March, 2024;
originally announced March 2024.
-
On the anticyclotomic Iwasawa theory of newforms at Eisenstein primes of semistable reduction
Authors:
Timo Keller,
Mulun Yin
Abstract:
Let $f$ be a newform of weight $k$ and level $N$ with trivial nebentypus. Let $\mathfrak{p}\nmid 2N$ be a maximal prime ideal of the coefficient ring of $f$ such that the self-dual twist of the mod-$\mathfrak{p}$ Galois representation of $f$ is reducible with constituents $φ,ψ$. Denote a decomposition group over the rational prime $p$ below $\mathfrak{p}$ by $G_p$. We remove the condition…
▽ More
Let $f$ be a newform of weight $k$ and level $N$ with trivial nebentypus. Let $\mathfrak{p}\nmid 2N$ be a maximal prime ideal of the coefficient ring of $f$ such that the self-dual twist of the mod-$\mathfrak{p}$ Galois representation of $f$ is reducible with constituents $φ,ψ$. Denote a decomposition group over the rational prime $p$ below $\mathfrak{p}$ by $G_p$. We remove the condition $φ|_{G_p} \neq \mathbf{1}, ω$ from [CGLS22], and generalize their results to newforms of arbitrary weights. As a consequence, we prove some Iwasawa main conjectures and get the $p$-part of the strong BSD conjecture for elliptic curves of analytic rank $0$ or $1$ over $\mathbf{Q}$ in this setting. In particular, non-trivial $p$-torsion is allowed in the Mordell--Weil group. Using Hida families, we prove a Iwasawa main conjecture for newforms of weight $2$ of multiplicative reduction at Eisenstein primes. In the above situations, we also get $p$-converse theorems to the theorems of Gross--Zagier--Kolyvagin. The $p$-converse theorems have applications to Goldfeld's conjecture in certain quadratic twist families of elliptic curves having a $3$-isogeny.
△ Less
Submitted 20 February, 2024;
originally announced February 2024.
-
Frequency-Domain Identification of Discrete-Time Systems using Sum-of-Rational Optimization
Authors:
Mohamed Abdalmoaty,
Jared Miller,
Mingzhou Yin,
Roy S. Smith
Abstract:
We propose a computationally tractable method for the identification of stable canonical discrete-time rational transfer function models, using frequency domain data. The problem is formulated as a global non-convex optimization problem whose objective function is the sum of weighted squared residuals at each observed frequency datapoint. Stability is enforced using a polynomial matrix inequality…
▽ More
We propose a computationally tractable method for the identification of stable canonical discrete-time rational transfer function models, using frequency domain data. The problem is formulated as a global non-convex optimization problem whose objective function is the sum of weighted squared residuals at each observed frequency datapoint. Stability is enforced using a polynomial matrix inequality constraint. The problem is solved globally by a moment-sum-of-squares hierarchy of semidefinite programs through a framework for sum-of-rational-functions optimization. Convergence of the moment-sum-of-squares program is guaranteed as the bound on the degree of the sum-of-squares polynomials approaches infinity. The performance of the proposed method is demonstrated using numerical simulation examples.
△ Less
Submitted 25 December, 2023;
originally announced December 2023.
-
Closed-Loop Identification of Stabilized Models Using Dual Input-Output Parameterization
Authors:
Ran Chen,
Amber Srivastava,
Mingzhou Yin,
Roy S. Smith
Abstract:
This paper introduces a dual input-output parameterization (dual IOP) for the identification of linear time-invariant systems from closed-loop data. It draws inspiration from the recent input-output parameterization developed to synthesize a stabilizing controller. The controller is parameterized in terms of closed-loop transfer functions, from the external disturbances to the input and output of…
▽ More
This paper introduces a dual input-output parameterization (dual IOP) for the identification of linear time-invariant systems from closed-loop data. It draws inspiration from the recent input-output parameterization developed to synthesize a stabilizing controller. The controller is parameterized in terms of closed-loop transfer functions, from the external disturbances to the input and output of the system, constrained to lie in a given subspace. Analogously, the dual IOP method parameterizes the unknown plant with analogous closed-loop transfer functions, also referred to as dual parameters. In this case, these closed-loop transfer functions are constrained to lie in an affine subspace guaranteeing that the identified plant is \emph{stabilized} by the known controller. Compared with existing closed-loop identification techniques guaranteeing closed-loop stability, such as the dual Youla parameterization, the dual IOP neither requires a doubly-coprime factorization of the controller nor a nominal plant that is stabilized by the controller. The dual IOP does not depend on the order and the state-space realization of the controller either, as in the dual system-level parameterization. Simulation shows that the dual IOP outperforms the existing benchmark methods.
△ Less
Submitted 15 November, 2023;
originally announced November 2023.
-
Enumerating pattern-avoiding permutations by leading terms
Authors:
Ömer Eğecioğlu,
Collier Gaiser,
Mei Yin
Abstract:
The number of 123-avoiding permutation on $\{1,2,\ldots,n\}$ with a fixed leading terms is counted by the ballot numbers. The same holds for $132$-avoiding permutations. These results were proved by Miner and Pak using the Robinson-Schensted-Knuth (RSK) correspondence to connect permutations with Dyck paths. In this paper, we first provide an alternate proof of these enumeration results via a dire…
▽ More
The number of 123-avoiding permutation on $\{1,2,\ldots,n\}$ with a fixed leading terms is counted by the ballot numbers. The same holds for $132$-avoiding permutations. These results were proved by Miner and Pak using the Robinson-Schensted-Knuth (RSK) correspondence to connect permutations with Dyck paths. In this paper, we first provide an alternate proof of these enumeration results via a direct counting argument. We then study the number of pattern-avoiding permutations with a fixed prefix of length $t\geq1$, generalizing the $t=1$ case. We find exact expressions for single and pairs of patterns of length three as well as the pair $3412$ and $3421$. These expressions depend on $t$, the extrema, and the order statistics. We also define $r$-Wilf equivalence for permutations with a single fixed leading term $r$, and classify the $r$-Wilf-equivalence classes for both classical and vincular patterns of length three.
△ Less
Submitted 24 June, 2024; v1 submitted 27 September, 2023;
originally announced September 2023.
-
Some enumerative properties of parking functions
Authors:
Richard P. Stanley,
Mei Yin
Abstract:
A parking function is a sequence $(a_1,\dots, a_n)$ of positive integers such that if $b_1\leq\cdots\leq b_n$ is the increasing rearrangement of $a_1,\dots,a_n$, then $b_i\leq i$ for $1\leq i\leq n$. In this paper we obtain some new results on the enumeration of parking functions. We will consider the joint distribution of several sets of statistics on parking functions. The distribution of most o…
▽ More
A parking function is a sequence $(a_1,\dots, a_n)$ of positive integers such that if $b_1\leq\cdots\leq b_n$ is the increasing rearrangement of $a_1,\dots,a_n$, then $b_i\leq i$ for $1\leq i\leq n$. In this paper we obtain some new results on the enumeration of parking functions. We will consider the joint distribution of several sets of statistics on parking functions. The distribution of most of these individual statistics is known, but the joint distributions are new. Parking functions of length $n$ are in bijection with labelled forests on the vertex set $[n]=\{1,2,\dots,n\}$ (or rooted trees on $[n]_0=\{0,1,\dots,n\}$ with root $0$), so our results can also be applied to labelled forests. Extensions of our techniques are discussed.
△ Less
Submitted 14 June, 2023;
originally announced June 2023.
-
Moments of Colored Permutation Statistics on Conjugacy Classes
Authors:
Jesse Campion Loth,
Michael Levet,
Kevin Liu,
Sheila Sundaram,
Mei Yin
Abstract:
In this paper, we consider the moments of statistics on conjugacy classes of the colored permutation groups $\mathfrak{S}_{n,r}=\mathbb{Z}_r\wr \mathfrak{S}_n$. We first show that any fixed moment coincides on all conjugacy classes where all cycles have sufficiently long length. Additionally, for permutation statistics that can be realized via a process we call symmetric extensions, these moments…
▽ More
In this paper, we consider the moments of statistics on conjugacy classes of the colored permutation groups $\mathfrak{S}_{n,r}=\mathbb{Z}_r\wr \mathfrak{S}_n$. We first show that any fixed moment coincides on all conjugacy classes where all cycles have sufficiently long length. Additionally, for permutation statistics that can be realized via a process we call symmetric extensions, these moments are polynomials in $n$. Finally, for the descent statistic on the hyperoctahedral group $B_n\cong \mathfrak{S}_{n,2}$, we show that its distribution on conjugacy classes without short cycles satisfies a central limit theorem. Our results build on and generalize previous work of Fulman (\textit{J. Comb. Theory Ser. A.}, 1998), Hamaker and Rhoades (arXiv, 2022), and Campion Loth, Levet, Liu, Stucky, Sundaram, and Yin (arXiv, 2023). In particular, our techniques utilize the combinatorial framework introduced by Campion Loth, Levet, Liu, Stucky, Sundaram, and Yin.
△ Less
Submitted 27 December, 2023; v1 submitted 19 May, 2023;
originally announced May 2023.
-
Blow-up phenomena for a class of extensible beam equations
Authors:
Gongwei Liu,
Mengyun Yin,
Suxia Xia
Abstract:
In this paper, we investigate the initial boundary value problem of the following nonlinear extensible beam equation with nonlinear damping term $$u_{t t}+Δ^2 u-M\left(\|\nabla u\|^2\right) Δu-Δu_t+\left|u_t\right|^{r-1} u_t=|u|^{p-1} u$$ which was considered by Yang et al. (Advanced Nonlinear Studies 2022; 22:436-468). We consider the problem with the nonlinear damping and establish the finite ti…
▽ More
In this paper, we investigate the initial boundary value problem of the following nonlinear extensible beam equation with nonlinear damping term $$u_{t t}+Δ^2 u-M\left(\|\nabla u\|^2\right) Δu-Δu_t+\left|u_t\right|^{r-1} u_t=|u|^{p-1} u$$ which was considered by Yang et al. (Advanced Nonlinear Studies 2022; 22:436-468). We consider the problem with the nonlinear damping and establish the finite time blow-up of the solution for the initial data at arbitrary high energy level, including the estimate lower and upper bounds of the blowup time. The result provides some affirmative answer to the open problems given in (Advanced Nonlinear Studies 2022; 22:436-468).
△ Less
Submitted 15 May, 2023;
originally announced May 2023.
-
Sampling planar tanglegrams and pairs of disjoint triangulations
Authors:
Alexander E. Black,
Kevin Liu,
Alex Mcdonough,
Garrett Nelson,
Michael C. Wigal,
Mei Yin,
Youngho Yoo
Abstract:
A tanglegram consists of two rooted binary trees and a perfect matching between their leaves, and a planar tanglegram is one that admits a layout with no crossings. We show that the problem of generating planar tanglegrams uniformly at random reduces to the corresponding problem for irreducible planar tanglegram layouts, which are known to be in bijection with pairs of disjoint triangulations of a…
▽ More
A tanglegram consists of two rooted binary trees and a perfect matching between their leaves, and a planar tanglegram is one that admits a layout with no crossings. We show that the problem of generating planar tanglegrams uniformly at random reduces to the corresponding problem for irreducible planar tanglegram layouts, which are known to be in bijection with pairs of disjoint triangulations of a convex polygon. We extend the flip operation on a single triangulation to a flip operation on pairs of disjoint triangulations. Interestingly, the resulting flip graph is both connected and regular, and hence a random walk on this graph converges to the uniform distribution. We also show that the restriction of the flip graph to the pairs with a fixed triangulation in either coordinate is connected, and give diameter bounds that are near optimal. Our results furthermore yield new insight into the flip graph of triangulations of a convex $n$-gon with a geometric interpretation on the associahedron.
△ Less
Submitted 11 April, 2023;
originally announced April 2023.
-
A Dual System-Level Parameterization for Identification from Closed-Loop Data
Authors:
Amber Srivastava,
Mingzhou Yin,
Andrea Iannelli,
Roy S. Smith
Abstract:
This work presents a dual system-level parameterization (D-SLP) method for closed-loop system identification. The recent system-level synthesis framework parameterizes all stabilizing controllers via linear constraints on closed-loop response functions, known as system-level parameters. It was demonstrated that several structural, locality, and communication constraints on the controller can be po…
▽ More
This work presents a dual system-level parameterization (D-SLP) method for closed-loop system identification. The recent system-level synthesis framework parameterizes all stabilizing controllers via linear constraints on closed-loop response functions, known as system-level parameters. It was demonstrated that several structural, locality, and communication constraints on the controller can be posed as convex constraints on these system-level parameters. In the current work, the identification problem is treated as a {\em dual} of the system-level synthesis problem. The plant model is identified from the dual system-level parameters associated to the plant. In comparison to existing closed-loop identification approaches (such as the dual-Youla parameterization), the D-SLP framework neither requires the knowledge of a nominal plant that is stabilized by the known controller, nor depends upon the choice of factorization of the nominal plant and the stabilizing controller. Numerical simulations demonstrate the efficacy of the proposed D-SLP method in terms of identification errors, compared to existing closed-loop identification techniques.
△ Less
Submitted 5 April, 2023;
originally announced April 2023.
-
Permutation Statistics in Conjugacy Classes of the Symmetric Group
Authors:
Jesse Campion Loth,
Michael Levet,
Kevin Liu,
Eric Nathan Stucky,
Sheila Sundaram,
Mei Yin
Abstract:
We introduce the notion of a weighted inversion statistic on the symmetric group, and examine its distribution on each conjugacy class. Our work generalizes the study of several common permutation statistics, including the number of inversions, the number of descents, the major index, and the number of excedances. As a consequence, we obtain explicit formulas for the first moments of several stati…
▽ More
We introduce the notion of a weighted inversion statistic on the symmetric group, and examine its distribution on each conjugacy class. Our work generalizes the study of several common permutation statistics, including the number of inversions, the number of descents, the major index, and the number of excedances. As a consequence, we obtain explicit formulas for the first moments of several statistics by conjugacy class. We also show that when the cycle lengths are sufficiently large, the higher moments of arbitrary permutation statistics are independent of the conjugacy class. Fulman (J. Comb. Theory Ser. A., 1998) previously established this result for major index and descents. We obtain these results, in part, by generalizing the techniques of Fulman (ibid.), and introducing the notion of permutation constraints. For permutation statistics that can be realized via symmetric constraints, we show that each moment is a polynomial in the degree of the symmetric group.
△ Less
Submitted 17 May, 2023; v1 submitted 2 January, 2023;
originally announced January 2023.
-
Probabilistic Parking Functions
Authors:
Irfan Durmić,
Alex Han,
Pamela E. Harris,
Rodrigo Ribeiro,
Mei Yin
Abstract:
We consider the notion of classical parking functions by introducing randomness and a new parking protocol, as inspired by the work presented in the paper ``Parking Functions: Choose your own adventure,'' (arXiv:2001.04817) by Carlson, Christensen, Harris, Jones, and Rodríguez. Among our results, we prove that the probability of obtaining a parking function, from a length $n$ preference vector, is…
▽ More
We consider the notion of classical parking functions by introducing randomness and a new parking protocol, as inspired by the work presented in the paper ``Parking Functions: Choose your own adventure,'' (arXiv:2001.04817) by Carlson, Christensen, Harris, Jones, and Rodríguez. Among our results, we prove that the probability of obtaining a parking function, from a length $n$ preference vector, is independent of the probabilistic parameter $p$. We also explore the properties of a preference vector given that it is a parking function and discuss the effect of the probabilistic parameter $p$. Of special interest is when $p=1/2$, where we demonstrate a sharp transition in some parking statistics. We also present several interesting combinatorial consequences of the parking protocol. In particular, we provide a combinatorial interpretation for the array described in OEIS A220884 as the expected number of preference sequences with a particular property related to occupied parking spots, which solves an open problem of Novelli and Thibon posed in 2020 (arXiv:1209.5959). Lastly, we connect our results to other weighted phenomena in combinatorics and provide further directions for research.
△ Less
Submitted 1 November, 2022;
originally announced November 2022.
-
Odd Covers of Graphs
Authors:
Calum Buchanan,
Alexander Clifton,
Eric Culver,
Jiaxi Nie,
Jason O'Neill,
Puck Rombach,
Mei Yin
Abstract:
Given a finite simple graph $G$, an odd cover of $G$ is a collection of complete bipartite graphs, or bicliques, in which each edge of $G$ appears in an odd number of bicliques and each non-edge of $G$ appears in an even number of bicliques. We denote the minimum cardinality of an odd cover of $G$ by $b_2(G)$ and prove that $b_2(G)$ is bounded below by half of the rank over $\mathbb{F}_2$ of the a…
▽ More
Given a finite simple graph $G$, an odd cover of $G$ is a collection of complete bipartite graphs, or bicliques, in which each edge of $G$ appears in an odd number of bicliques and each non-edge of $G$ appears in an even number of bicliques. We denote the minimum cardinality of an odd cover of $G$ by $b_2(G)$ and prove that $b_2(G)$ is bounded below by half of the rank over $\mathbb{F}_2$ of the adjacency matrix of $G$. We show that this lower bound is tight in the case when $G$ is a bipartite graph and almost tight when $G$ is an odd cycle. However, we also present an infinite family of graphs which shows that this lower bound can be arbitrarily far away from $b_2(G)$.
Babai and Frankl (1992) proposed the "odd cover problem," which in our language is equivalent to determining $b_2(K_n)$. Radhakrishnan, Sen, and Vishwanathan (2000) determined $b_2(K_n)$ for an infinite but density zero subset of positive integers $n$. In this paper, we determine $b_2(K_n)$ for a density $3/8$ subset of the positive integers.
△ Less
Submitted 20 February, 2022;
originally announced February 2022.
-
Ehrhart Theory of Paving and Panhandle Matroids
Authors:
Derek Hanely,
Jeremy L. Martin,
Daniel McGinnis,
Dane Miyata,
George D. Nasr,
Andrés R. Vindas-Meléndez,
Mei Yin
Abstract:
We show that the base polytope $P_M$ of any paving matroid $M$ can be systematically obtained from a hypersimplex by slicing off certain subpolytopes, namely base polytopes of lattice path matroids corresponding to panhandle-shaped Ferrers diagrams. We calculate the Ehrhart polynomials of these matroids and consequently write down the Ehrhart polynomial of $P_M$, starting with Katzman's formula fo…
▽ More
We show that the base polytope $P_M$ of any paving matroid $M$ can be systematically obtained from a hypersimplex by slicing off certain subpolytopes, namely base polytopes of lattice path matroids corresponding to panhandle-shaped Ferrers diagrams. We calculate the Ehrhart polynomials of these matroids and consequently write down the Ehrhart polynomial of $P_M$, starting with Katzman's formula for the Ehrhart polynomial of a hypersimplex. The method builds on and generalizes Ferroni's work on sparse paving matroids. Combinatorially, our construction corresponds to constructing a uniform matroid from a paving matroid by iterating the operation of stressed-hyperplane relaxation introduced by Ferroni, Nasr, and Vecchi, which generalizes the standard matroid-theoretic notion of circuit-hyperplane relaxation. We present evidence that panhandle matroids are Ehrhart positive and describe a conjectured combinatorial formula involving chain forests and Eulerian numbers from which Ehrhart positivity of panhandle matroids will follow. As an application of the main result, we calculate the Ehrhart polynomials of matroids associated with Steiner systems and finite projective planes, and show that they depend only on their design-theoretic parameters: for example, while projective planes of the same order need not have isomorphic matroids, their base polytopes must be Ehrhart equivalent.
△ Less
Submitted 25 July, 2023; v1 submitted 28 January, 2022;
originally announced January 2022.
-
Parking functions, multi-shuffle, and asymptotic phenomena
Authors:
Mei Yin
Abstract:
Given a positive-integer-valued vector $u=(u_1, \dots, u_m)$ with $u_1<\cdots<u_m$. A $u$-parking function of length $m$ is a sequence $π=(π_1, \dots, π_m)$ of positive integers whose non-decreasing rearrangement $(λ_1, \dots, λ_m)$ satisfies $λ_i\leq u_i$ for all $1\leq i\leq m$. We introduce a combinatorial construction termed a parking function multi-shuffle to generic $u$-parking functions and…
▽ More
Given a positive-integer-valued vector $u=(u_1, \dots, u_m)$ with $u_1<\cdots<u_m$. A $u$-parking function of length $m$ is a sequence $π=(π_1, \dots, π_m)$ of positive integers whose non-decreasing rearrangement $(λ_1, \dots, λ_m)$ satisfies $λ_i\leq u_i$ for all $1\leq i\leq m$. We introduce a combinatorial construction termed a parking function multi-shuffle to generic $u$-parking functions and obtain an explicit characterization of multiple parking coordinates. As an application, we derive various asymptotic probabilistic properties of a uniform $u$-parking function when $u_i=cm+ib$. The asymptotic scenario in the generic situation $c>0$ is in sharp contrast with that of the special situation $c=0$.
△ Less
Submitted 4 December, 2021;
originally announced December 2021.
-
A physics-informed variational DeepONet for predicting the crack path in brittle materials
Authors:
Somdatta Goswami,
Minglang Yin,
Yue Yu,
George Karniadakis
Abstract:
Failure trajectories, identifying the probable failure zones, and damage statistics are some of the key quantities of relevance in brittle fracture applications. High-fidelity numerical solvers that reliably estimate these relevant quantities exist but they are computationally demanding requiring a high resolution of the crack. Moreover, independent intensive simulations need to be carried out eve…
▽ More
Failure trajectories, identifying the probable failure zones, and damage statistics are some of the key quantities of relevance in brittle fracture applications. High-fidelity numerical solvers that reliably estimate these relevant quantities exist but they are computationally demanding requiring a high resolution of the crack. Moreover, independent intensive simulations need to be carried out even for a small change in domain parameters and/or material properties. Therefore, fast and generalizable surrogate models are needed to alleviate the computational burden but the discontinuous nature of fracture mechanics presents a major challenge to developing such models. We propose a physics-informed variational formulation of DeepONet (V-DeepONet) for brittle fracture analysis. V-DeepONet is trained to map the initial configuration of the defect to the relevant fields of interests (e.g., damage and displacement fields). Once the network is trained, the entire global solution can be rapidly obtained for any initial crack configuration and loading steps on that domain. While the original DeepONet is solely data-driven, we take a different path to train the V-DeepONet by imposing the governing equations in variational form and we also use some labelled data. We demonstrate the effectiveness of V-DeepOnet through two benchmarks of brittle fracture, and we verify its accuracy using results from high-fidelity solvers. Encoding the physical laws and also some data to train the network renders the surrogate model capable of accurately performing both interpolation and extrapolation tasks, considering that fracture modeling is very sensitive to fluctuations. The proposed hybrid training of V-DeepONet is superior to state-of-the-art methods and can be applied to a wide array of dynamical systems with complex responses.
△ Less
Submitted 13 September, 2021; v1 submitted 16 August, 2021;
originally announced August 2021.
-
Parking functions: Interdisciplinary connections
Authors:
Mei Yin
Abstract:
Suppose that $m$ drivers each choose a preferred parking space in a linear car park with $n$ spots. In order, each driver goes to their chosen spot and parks there if possible, and otherwise takes the next available spot if it exists. If all drivers park successfully, the sequence of choices is called a parking function. Classical parking functions correspond to the case $m=n$.
We investigate va…
▽ More
Suppose that $m$ drivers each choose a preferred parking space in a linear car park with $n$ spots. In order, each driver goes to their chosen spot and parks there if possible, and otherwise takes the next available spot if it exists. If all drivers park successfully, the sequence of choices is called a parking function. Classical parking functions correspond to the case $m=n$.
We investigate various probabilistic properties of a uniform parking function. Through a combinatorial construction termed a parking function multi-shuffle, we give a formula for the law of multiple coordinates in the generic situation $m \lesssim n$. We further deduce all possible covariances, between two coordinates, between a coordinate and an unattempted spot, and between two unattempted spots. This asymptotic scenario in the generic situation $m \lesssim n$ is in sharp contrast with that of the special situation $m=n$.
A generalization of parking functions called interval parking functions is also studied, in which each driver is willing to park only in a fixed interval of spots. We construct a family of bijections between interval parking functions with $n$ cars and $n$ spots and edge-labeled spanning trees with $n+1$ vertices and a specified root.
△ Less
Submitted 4 October, 2021; v1 submitted 4 July, 2021;
originally announced July 2021.
-
Parking functions: From combinatorics to probability
Authors:
Richard Kenyon,
Mei Yin
Abstract:
Suppose that $m$ drivers each choose a preferred parking space in a linear car park with $n$ spots. In order, each driver goes to their chosen spot and parks there if possible, and otherwise takes the next available spot if it exists. If all drivers park successfully, the sequence of choices is called a parking function. Classical parking functions correspond to the case $m=n$; we study here combi…
▽ More
Suppose that $m$ drivers each choose a preferred parking space in a linear car park with $n$ spots. In order, each driver goes to their chosen spot and parks there if possible, and otherwise takes the next available spot if it exists. If all drivers park successfully, the sequence of choices is called a parking function. Classical parking functions correspond to the case $m=n$; we study here combinatorial and probabilistic aspects of this generalized case.
We construct a family of bijections between parking functions $\text{PF}(m, n)$ with $m$ cars and $n$ spots and spanning forests $\mathscr{F}(n+1, n+1-m)$ with $n+1$ vertices and $n+1-m$ distinct trees having specified roots. This leads to a bijective correspondence between $\text{PF}(m, n)$ and monomial terms in the associated Tutte polynomial of a disjoint union of $n-m+1$ complete graphs. We present an identity between the "inversion enumerator" of spanning forests with fixed roots and the "displacement enumerator" of parking functions. The displacement is then related to the number of graphs on $n+1$ labeled vertices with a fixed number of edges, where the graph has $n+1-m$ disjoint rooted components with specified roots.
We investigate various probabilistic properties of a uniform parking function, giving a formula for the law of a single coordinate. As a side result we obtain a recurrence relation for the displacement enumerator. Adapting known results on random linear probes, we further deduce the covariance between two coordinates when $m=n$.
△ Less
Submitted 31 March, 2021;
originally announced March 2021.
-
Remarks on power-law random graphs
Authors:
Mei Yin
Abstract:
The theory of graphons is an important tool in understanding properties of large networks. We investigate a power-law random graph model and cast it in the graphon framework. The distinctively different structures of the limit graph are explored in detail in the sub-critical and super-critical regimes. In the sub-critical regime, the graph is empty with high probability, and in the rare event that…
▽ More
The theory of graphons is an important tool in understanding properties of large networks. We investigate a power-law random graph model and cast it in the graphon framework. The distinctively different structures of the limit graph are explored in detail in the sub-critical and super-critical regimes. In the sub-critical regime, the graph is empty with high probability, and in the rare event that it is non-empty, it consists of a single edge. Contrarily, in the super-critical regime, a non-trivial random graph exists in the limit, and it serves as an uncovered boundary case between different types of graph convergence.
△ Less
Submitted 18 August, 2022; v1 submitted 12 August, 2020;
originally announced August 2020.
-
Interval parking functions
Authors:
Emma Colaric,
Ryan DeMuse,
Jeremy L. Martin,
Mei Yin
Abstract:
Interval parking functions (IPFs) are a generalization of ordinary parking functions in which each car is willing to park only in a fixed interval of spaces. Each interval parking function can be expressed as a pair $(a,b)$, where $a$ is a parking function and $b$ is a dual parking function. We say that a pair of permutations $(x,y)$ is \emph{reachable} if there is an IPF $(a,b)$ such that $x,y$ a…
▽ More
Interval parking functions (IPFs) are a generalization of ordinary parking functions in which each car is willing to park only in a fixed interval of spaces. Each interval parking function can be expressed as a pair $(a,b)$, where $a$ is a parking function and $b$ is a dual parking function. We say that a pair of permutations $(x,y)$ is \emph{reachable} if there is an IPF $(a,b)$ such that $x,y$ are the outcomes of $a,b$, respectively, as parking functions. Reachability is reflexive and antisymmetric, but not in general transitive. We prove that its transitive closure, the \emph{pseudoreachability order}, is precisely the bubble-sort order on the symmetric group $\Sym_n$, which can be expressed in terms of the normal form of a permutation in the sense of du~Cloux; in particular, it is isomorphic to the product of chains of lengths $2,\dots,n$. It is thus seen to be a special case of Armstrong's sorting order, which lies between the Bruhat and (left) weak orders.
△ Less
Submitted 28 October, 2020; v1 submitted 16 June, 2020;
originally announced June 2020.
-
Robust Adaptive Model Predictive Control with Worst-Case Cost
Authors:
Anilkumar Parsi,
Andrea Iannelli,
Mingzhou Yin,
Mohammad Khosravi,
Roy S. Smith
Abstract:
A robust adaptive model predictive control (MPC) algorithm is presented for linear, time invariant systems with unknown dynamics and subject to bounded measurement noise. The system is characterized by an impulse response model, which is assumed to lie within a bounded set called the feasible system set. Online set-membership identification is used to reduce uncertainty in the impulse response. In…
▽ More
A robust adaptive model predictive control (MPC) algorithm is presented for linear, time invariant systems with unknown dynamics and subject to bounded measurement noise. The system is characterized by an impulse response model, which is assumed to lie within a bounded set called the feasible system set. Online set-membership identification is used to reduce uncertainty in the impulse response. In the MPC scheme, robust constraints are enforced to ensure constraint satisfaction for all the models in the feasible set. The performance objective is formulated as a worst-case cost with respect to the modeling uncertainties. That is, at each time step an optimization problem is solved in which the control input is optimized for the worst-case plant in the uncertainty set. The performance of the proposed algorithm is compared to an adaptive MPC algorithm from the literature using Monte-Carlo simulations.
△ Less
Submitted 19 November, 2019;
originally announced November 2019.
-
Dimension reduction in vertex-weighted exponential random graphs
Authors:
Ryan DeMuse,
Mei Yin
Abstract:
We investigate the behavior of vertex-weighted exponential random graphs. We show that vertex-weighted exponential random graphs with edge weights induced by products of independent vertex weights are approximate mixtures of graphs whose vertex weight vector is a near fixed point of a certain vector equation. For graphs with Hamiltonians counting cliques, it is demonstrated that, under appropriate…
▽ More
We investigate the behavior of vertex-weighted exponential random graphs. We show that vertex-weighted exponential random graphs with edge weights induced by products of independent vertex weights are approximate mixtures of graphs whose vertex weight vector is a near fixed point of a certain vector equation. For graphs with Hamiltonians counting cliques, it is demonstrated that, under appropriate conditions, every solution to this equation is close to a block vector with a small number of communities. We prove that for the cases of positive weights and small weights in the Hamiltonian in particular, the vector equation has a unique solution. Lastly, the behavior of vertex-weighted exponential random graphs counting triangles is studied in detail and the solution to the vector equation is shown to approach the zero vector as the weight diverges to negative infinity for sufficiently large networks.
△ Less
Submitted 22 August, 2019;
originally announced August 2019.
-
A Theoretical Case Study of Structured Variational Inference for Community Detection
Authors:
Mingzhang Yin,
Y. X. Rachel Wang,
Purnamrita Sarkar
Abstract:
Mean-field variational inference (MFVI) has been widely applied in large scale Bayesian inference. However MFVI, which assumes a product distribution on the latent variables, often leads to objective functions with many local optima, making optimization algorithms sensitive to initialization. In this paper, we study the advantage of structured variational inference for the two class Stochastic Blo…
▽ More
Mean-field variational inference (MFVI) has been widely applied in large scale Bayesian inference. However MFVI, which assumes a product distribution on the latent variables, often leads to objective functions with many local optima, making optimization algorithms sensitive to initialization. In this paper, we study the advantage of structured variational inference for the two class Stochastic Blockmodel. The variational distribution is constructed to have pairwise dependency structure on the nodes of the network. We prove that, in a broad density regime and for general random initializations, unlike MFVI, the class labels estimated from our method converge to the ground truth with high probability, when the model parameters are known, estimated within a reasonable range or jointly optimized with the variational parameters. In addition, empirically we demonstrate structured VI is more robust compared with MFVI when the graph is sparse and the signal to noise ratio is low. The paper takes a first step towards understanding the importance of dependency structure in variational inference for community detection.
△ Less
Submitted 29 February, 2020; v1 submitted 28 July, 2019;
originally announced July 2019.
-
Mixing Time of Vertex-Weighted Exponential Random Graphs
Authors:
Ryan DeMuse,
Terry Easlick,
Mei Yin
Abstract:
Exponential random graph models have become increasingly important in the study of modern networks ranging from social networks, economic networks, to biological networks. They seek to capture a wide variety of common network tendencies such as connectivity and reciprocity through local graph properties. Sampling from these exponential distributions is crucial for parameter estimation, hypothesis…
▽ More
Exponential random graph models have become increasingly important in the study of modern networks ranging from social networks, economic networks, to biological networks. They seek to capture a wide variety of common network tendencies such as connectivity and reciprocity through local graph properties. Sampling from these exponential distributions is crucial for parameter estimation, hypothesis testing, as well as understanding the features of the network in question. We inspect the efficiency of a popular sampling technique, the Glauber dynamics, for vertex-weighted exponential random graphs. Letting $n$ be the number of vertices in the graph, we identify a region in the parameter space where the mixing time for the Glauber dynamics is $Θ(n \log n)$ (the high temperature phase) and a complement region where the mixing time is exponentially slow on the order of $e^{Ω(n)}$ (the low temperature phase). Lastly, we give evidence that along a critical curve in the parameter space the mixing time is $O(n^{2/3})$.
△ Less
Submitted 20 December, 2017;
originally announced December 2017.
-
Ground States for Exponential Random Graphs
Authors:
Rajinder Mavi,
Mei Yin
Abstract:
We propose a perturbative method to estimate the normalization constant in exponential random graph models as the weighting parameters approach infinity. As an application, we give evidence of discontinuity in natural parametrization along the critical directions of the edge-triangle model.
We propose a perturbative method to estimate the normalization constant in exponential random graph models as the weighting parameters approach infinity. As an application, we give evidence of discontinuity in natural parametrization along the critical directions of the edge-triangle model.
△ Less
Submitted 4 January, 2018; v1 submitted 27 June, 2017;
originally announced June 2017.
-
Phase Transitions in Edge-Weighted Exponential Random Graphs: Near-Degeneracy and Universality
Authors:
Ryan DeMuse,
Danielle Larcomb,
Mei Yin
Abstract:
Conventionally used exponential random graphs cannot directly model weighted networks as the underlying probability space consists of simple graphs only. Since many substantively important networks are weighted, this limitation is especially problematic. We extend the existing exponential framework by proposing a generic common distribution for the edge weights. Minimal assumptions are placed on t…
▽ More
Conventionally used exponential random graphs cannot directly model weighted networks as the underlying probability space consists of simple graphs only. Since many substantively important networks are weighted, this limitation is especially problematic. We extend the existing exponential framework by proposing a generic common distribution for the edge weights. Minimal assumptions are placed on the distribution, that is, it is non-degenerate and supported on the unit interval. By doing so, we recognize the essential properties associated with near-degeneracy and universality in edge-weighted exponential random graphs.
△ Less
Submitted 23 February, 2018; v1 submitted 7 June, 2017;
originally announced June 2017.
-
Convergence Analysis of Gradient EM for Multi-component Gaussian Mixture
Authors:
Bowei Yan,
Mingzhang Yin,
Purnamrita Sarkar
Abstract:
In this paper, we study convergence properties of the gradient Expectation-Maximization algorithm \cite{lange1995gradient} for Gaussian Mixture Models for general number of clusters and mixing coefficients. We derive the convergence rate depending on the mixing coefficients, minimum and maximum pairwise distances between the true centers and dimensionality and number of components; and obtain a ne…
▽ More
In this paper, we study convergence properties of the gradient Expectation-Maximization algorithm \cite{lange1995gradient} for Gaussian Mixture Models for general number of clusters and mixing coefficients. We derive the convergence rate depending on the mixing coefficients, minimum and maximum pairwise distances between the true centers and dimensionality and number of components; and obtain a near-optimal local contraction radius. While there have been some recent notable works that derive local convergence rates for EM in the two equal mixture symmetric GMM, in the more general case, the derivations need structurally different and non-trivial arguments. We use recent tools from learning theory and empirical processes to achieve our theoretical results.
△ Less
Submitted 4 December, 2017; v1 submitted 23 May, 2017;
originally announced May 2017.
-
Phase transitions in edge-weighted exponential random graphs
Authors:
Mei Yin
Abstract:
The exponential family of random graphs represents an important and challenging class of network models. Despite their flexibility, conventionally used exponential random graphs have one shortcoming. They cannot directly model weighted networks as the underlying probability space consists of simple graphs only. Since many substantively important networks are weighted, this limitation is especially…
▽ More
The exponential family of random graphs represents an important and challenging class of network models. Despite their flexibility, conventionally used exponential random graphs have one shortcoming. They cannot directly model weighted networks as the underlying probability space consists of simple graphs only. Since many substantively important networks are weighted, this limitation is especially problematic. We extend the existing exponential framework by proposing a generic common distribution for the edge weights and rigorously analyze the associated phase transitions and critical phenomena. We then apply these general results to get concrete answers in exponential random graph models where the edge weights are uniformly distributed.
△ Less
Submitted 14 July, 2016;
originally announced July 2016.
-
Large deviations and exact asymptotics for constrained exponential random graphs
Authors:
Mei Yin
Abstract:
We present a technique for approximating generic normalization constants subject to constraints. The method is then applied to derive the exact asymptotics for the conditional normalization constant of constrained exponential random graphs.
We present a technique for approximating generic normalization constants subject to constraints. The method is then applied to derive the exact asymptotics for the conditional normalization constant of constrained exponential random graphs.
△ Less
Submitted 31 July, 2015; v1 submitted 18 December, 2014;
originally announced December 2014.
-
Reciprocity in directed networks
Authors:
Mei Yin,
Lingjiong Zhu
Abstract:
Reciprocity is an important characteristic of directed networks and has been widely used in the modeling of World Wide Web, email, social, and other complex networks. In this paper, we take a statistical physics point of view and study the limiting entropy and free energy densities from the microcanonical ensemble, the canonical ensemble, and the grand canonical ensemble whose sufficient statistic…
▽ More
Reciprocity is an important characteristic of directed networks and has been widely used in the modeling of World Wide Web, email, social, and other complex networks. In this paper, we take a statistical physics point of view and study the limiting entropy and free energy densities from the microcanonical ensemble, the canonical ensemble, and the grand canonical ensemble whose sufficient statistics are given by edge and reciprocal densities. The sparse case is also studied for the grand canonical ensemble. Extensions to more general reciprocal models including reciprocal triangle and star densities will likewise be discussed.
△ Less
Submitted 21 December, 2015; v1 submitted 5 December, 2014;
originally announced December 2014.
-
Asymptotics for Sparse Exponential Random Graph Models
Authors:
Mei Yin,
Lingjiong Zhu
Abstract:
We study the asymptotics for sparse exponential random graph models where the parameters may depend on the number of vertices of the graph. We obtain exact estimates for the mean and variance of the limiting probability distribution and the limiting log partition function of the edge-(single)-star model. They are in sharp contrast to the corresponding asymptotics in dense exponential random graph…
▽ More
We study the asymptotics for sparse exponential random graph models where the parameters may depend on the number of vertices of the graph. We obtain exact estimates for the mean and variance of the limiting probability distribution and the limiting log partition function of the edge-(single)-star model. They are in sharp contrast to the corresponding asymptotics in dense exponential random graph models. Similar analysis is done for directed sparse exponential random graph models parametrized by edges and multiple outward stars.
△ Less
Submitted 19 April, 2016; v1 submitted 17 November, 2014;
originally announced November 2014.
-
On the asymptotics of constrained exponential random graphs
Authors:
Richard Kenyon,
Mei Yin
Abstract:
The unconstrained exponential family of random graphs assumes no prior knowledge of the graph before sampling, but it is natural to consider situations where partial information about the graph is known, for example the total number of edges. What does a typical random graph look like, if drawn from an exponential model subject to such constraints? Will there be a similar phase transition phenomen…
▽ More
The unconstrained exponential family of random graphs assumes no prior knowledge of the graph before sampling, but it is natural to consider situations where partial information about the graph is known, for example the total number of edges. What does a typical random graph look like, if drawn from an exponential model subject to such constraints? Will there be a similar phase transition phenomenon (as one varies the parameters) as that which occurs in the unconstrained exponential model? We present some general results for this constrained model and then apply them to get concrete answers in the edge-triangle model with fixed density of edges.
△ Less
Submitted 21 June, 2016; v1 submitted 13 June, 2014;
originally announced June 2014.
-
Asymptotic quantization of exponential random graphs
Authors:
Mei Yin,
Alessandro Rinaldo,
Sukhada Fadnavis
Abstract:
We describe the asymptotic properties of the edge-triangle exponential random graph model as the natural parameters diverge along straight lines. We show that as we continuously vary the slopes of these lines, a typical graph drawn from this model exhibits quantized behavior, jumping from one complete multipartite graph to another, and the jumps happen precisely at the normal lines of a polyhedral…
▽ More
We describe the asymptotic properties of the edge-triangle exponential random graph model as the natural parameters diverge along straight lines. We show that as we continuously vary the slopes of these lines, a typical graph drawn from this model exhibits quantized behavior, jumping from one complete multipartite graph to another, and the jumps happen precisely at the normal lines of a polyhedral set with infinitely many facets. As a result, we provide a complete description of all asymptotic extremal behaviors of the model.
△ Less
Submitted 11 January, 2016; v1 submitted 7 November, 2013;
originally announced November 2013.
-
Critical phenomena in exponential random graphs
Authors:
Mei Yin
Abstract:
The exponential family of random graphs is one of the most promising class of network models. Dependence between the random edges is defined through certain finite subgraphs, analogous to the use of potential energy to provide dependence between particle states in a grand canonical ensemble of statistical physics. By adjusting the specific values of these subgraph densities, one can analyze the in…
▽ More
The exponential family of random graphs is one of the most promising class of network models. Dependence between the random edges is defined through certain finite subgraphs, analogous to the use of potential energy to provide dependence between particle states in a grand canonical ensemble of statistical physics. By adjusting the specific values of these subgraph densities, one can analyze the influence of various local features on the global structure of the network. Loosely put, a phase transition occurs when a singularity arises in the limiting free energy density, as it is the generating function for the limiting expectations of all thermodynamic observables. We derive the full phase diagram for a large family of 3-parameter exponential random graph models with attraction and show that they all consist of a first order surface phase transition bordered by a second order critical curve.
△ Less
Submitted 19 September, 2013; v1 submitted 14 August, 2012;
originally announced August 2012.
-
A Markov chain approach to renormalization group transformations
Authors:
Mei Yin
Abstract:
We aim at an explicit characterization of the renormalized Hamiltonian after decimation transformation of a one-dimensional Ising-type Hamiltonian with a nearest-neighbor interaction and a magnetic field term. To facilitate a deeper understanding of the decimation effect, we translate the renormalization flow on the Ising Hamiltonian into a flow on the associated Markov chains through the Markov-G…
▽ More
We aim at an explicit characterization of the renormalized Hamiltonian after decimation transformation of a one-dimensional Ising-type Hamiltonian with a nearest-neighbor interaction and a magnetic field term. To facilitate a deeper understanding of the decimation effect, we translate the renormalization flow on the Ising Hamiltonian into a flow on the associated Markov chains through the Markov-Gibbs equivalence. Two different methods are used to verify the well-known conjecture that the eigenvalues of the linearization of this renormalization transformation about the fixed point bear important information about all six of the critical exponents. This illustrates the universality property of the renormalization group map in this case.
△ Less
Submitted 6 December, 2012; v1 submitted 8 May, 2012;
originally announced May 2012.
-
A cluster expansion approach to exponential random graph models
Authors:
Mei Yin
Abstract:
The exponential family of random graphs is among the most widely-studied network models. We show that any exponential random graph model may alternatively be viewed as a lattice gas model with a finite Banach space norm. The system may then be treated by cluster expansion methods from statistical mechanics. In particular, we derive a convergent power series expansion for the limiting free energy i…
▽ More
The exponential family of random graphs is among the most widely-studied network models. We show that any exponential random graph model may alternatively be viewed as a lattice gas model with a finite Banach space norm. The system may then be treated by cluster expansion methods from statistical mechanics. In particular, we derive a convergent power series expansion for the limiting free energy in the case of small parameters. Since the free energy is the generating function for the expectations of other random variables, this characterizes the structure and behavior of the limiting network in this parameter region.
△ Less
Submitted 26 April, 2012; v1 submitted 24 February, 2012;
originally announced February 2012.
-
Phase transitions in exponential random graphs
Authors:
Charles Radin,
Mei Yin
Abstract:
We derive the full phase diagram for a large family of two-parameter exponential random graph models, each containing a first order transition curve ending in a critical point.
We derive the full phase diagram for a large family of two-parameter exponential random graph models, each containing a first order transition curve ending in a critical point.
△ Less
Submitted 5 December, 2013; v1 submitted 2 August, 2011;
originally announced August 2011.