-
How to Build the Virtual Cell with Artificial Intelligence: Priorities and Opportunities
Authors:
Charlotte Bunne,
Yusuf Roohani,
Yanay Rosen,
Ankit Gupta,
Xikun Zhang,
Marcel Roed,
Theo Alexandrov,
Mohammed AlQuraishi,
Patricia Brennan,
Daniel B. Burkhardt,
Andrea Califano,
Jonah Cool,
Abby F. Dernburg,
Kirsty Ewing,
Emily B. Fox,
Matthias Haury,
Amy E. Herr,
Eric Horvitz,
Patrick D. Hsu,
Viren Jain,
Gregory R. Johnson,
Thomas Kalil,
David R. Kelley,
Shana O. Kelley,
Anna Kreshuk
, et al. (17 additional authors not shown)
Abstract:
The cell is arguably the most fundamental unit of life and is central to understanding biology. Accurate modeling of cells is important for this understanding as well as for determining the root causes of disease. Recent advances in artificial intelligence (AI), combined with the ability to generate large-scale experimental data, present novel opportunities to model cells. Here we propose a vision…
▽ More
The cell is arguably the most fundamental unit of life and is central to understanding biology. Accurate modeling of cells is important for this understanding as well as for determining the root causes of disease. Recent advances in artificial intelligence (AI), combined with the ability to generate large-scale experimental data, present novel opportunities to model cells. Here we propose a vision of leveraging advances in AI to construct virtual cells, high-fidelity simulations of cells and cellular systems under different conditions that are directly learned from biological data across measurements and scales. We discuss desired capabilities of such AI Virtual Cells, including generating universal representations of biological entities across scales, and facilitating interpretable in silico experiments to predict and understand their behavior using virtual instruments. We further address the challenges, opportunities and requirements to realize this vision including data needs, evaluation strategies, and community standards and engagement to ensure biological accuracy and broad utility. We envision a future where AI Virtual Cells help identify new drug targets, predict cellular responses to perturbations, as well as scale hypothesis exploration. With open science collaborations across the biomedical ecosystem that includes academia, philanthropy, and the biopharma and AI industries, a comprehensive predictive understanding of cell mechanisms and interactions has come into reach.
△ Less
Submitted 14 October, 2024; v1 submitted 17 September, 2024;
originally announced September 2024.
-
The distribution of the length of the longest path in random acyclic orientations of a complete bipartite graph
Authors:
Jessica Khera,
Erik Lundberg
Abstract:
Randomly sampling an acyclic orientation on the complete bipartite graph $K_{n,k}$ with parts of size $n$ and $k$, we investigate the length of the longest path. We provide a probability generating function for the distribution of the longest path length, and we use Analytic Combinatorics to perform asymptotic analysis of the probability distribution in the case of equal part sizes $n = k$ tending…
▽ More
Randomly sampling an acyclic orientation on the complete bipartite graph $K_{n,k}$ with parts of size $n$ and $k$, we investigate the length of the longest path. We provide a probability generating function for the distribution of the longest path length, and we use Analytic Combinatorics to perform asymptotic analysis of the probability distribution in the case of equal part sizes $n = k$ tending toward infinity. We show that the distribution is asymptotically Gaussian, and we obtain precise asymptotics for the mean and variance. These results address a question asked by Peter J. Cameron.
Keywords: bipartite graph, directed graph, random graph, acyclic orientation, poly-Bernoulli numbers, lonesum matrices, generating function, analytic combinatorics, asymptotics.
△ Less
Submitted 22 August, 2024;
originally announced August 2024.
-
Enabling Global Image Data Sharing in the Life Sciences
Authors:
Peter Bajcsy,
Sreenivas Bhattiprolu,
Katy Boerner,
Beth A Cimini,
Lucy Collinson,
Jan Ellenberg,
Reto Fiolka,
Maryellen Giger,
Wojtek Goscinski,
Matthew Hartley,
Nathan Hotaling,
Rick Horwitz,
Florian Jug,
Anna Kreshuk,
Emma Lundberg,
Aastha Mathur,
Kedar Narayan,
Shuichi Onami,
Anne L. Plant,
Fred Prior,
Jason Swedlow,
Adam Taylor,
Antje Keppler
Abstract:
Coordinated collaboration is essential to realize the added value of and infrastructure requirements for global image data sharing in the life sciences. In this White Paper, we take a first step at presenting some of the most common use cases as well as critical/emerging use cases of (including the use of artificial intelligence for) biological and medical image data, which would benefit tremendou…
▽ More
Coordinated collaboration is essential to realize the added value of and infrastructure requirements for global image data sharing in the life sciences. In this White Paper, we take a first step at presenting some of the most common use cases as well as critical/emerging use cases of (including the use of artificial intelligence for) biological and medical image data, which would benefit tremendously from better frameworks for sharing (including technical, resourcing, legal, and ethical aspects). In the second half of this paper, we paint an ideal world scenario for how global image data sharing could work and benefit all life sciences and beyond. As this is still a long way off, we conclude by suggesting several concrete measures directed toward our institutions, existing imaging communities and data initiatives, and national funders, as well as publishers. Our vision is that within the next ten years, most researchers in the world will be able to make their datasets openly available and use quality image data of interest to them for their research and benefit. This paper is published in parallel with a companion White Paper entitled Harmonizing the Generation and Pre-publication Stewardship of FAIR Image Data, which addresses challenges and opportunities related to producing well-documented and high-quality image data that is ready to be shared. The driving goal is to address remaining challenges and democratize access to everyday practices and tools for a spectrum of biomedical researchers, regardless of their expertise, access to resources, and geographical location.
△ Less
Submitted 9 August, 2024; v1 submitted 23 January, 2024;
originally announced January 2024.
-
Harmonizing the Generation and Pre-publication Stewardship of FAIR Image Data
Authors:
Nikki Bialy,
Frank Alber,
Brenda Andrews,
Michael Angelo,
Brian Beliveau,
Lacramioara Bintu,
Alistair Boettiger,
Ulrike Boehm,
Claire M. Brown,
Mahmoud Bukar Maina,
James J. Chambers,
Beth A. Cimini,
Kevin Eliceiri,
Rachel Errington,
Orestis Faklaris,
Nathalie Gaudreault,
Ronald N. Germain,
Wojtek Goscinski,
David Grunwald,
Michael Halter,
Dorit Hanein,
John W. Hickey,
Judith Lacoste,
Alex Laude,
Emma Lundberg
, et al. (22 additional authors not shown)
Abstract:
Together with the molecular knowledge of genes and proteins, biological images promise to significantly enhance the scientific understanding of complex cellular systems and to advance predictive and personalized therapeutic products for human health. For this potential to be realized, quality-assured image data must be shared among labs at a global scale to be compared, pooled, and reanalyzed, thu…
▽ More
Together with the molecular knowledge of genes and proteins, biological images promise to significantly enhance the scientific understanding of complex cellular systems and to advance predictive and personalized therapeutic products for human health. For this potential to be realized, quality-assured image data must be shared among labs at a global scale to be compared, pooled, and reanalyzed, thus unleashing untold potential beyond the original purpose for which the data was generated. There are two broad sets of requirements to enable image data sharing in the life sciences. One set of requirements is articulated in the companion White Paper entitled Enabling Global Image Data Sharing in the Life Sciences, which is published in parallel and addresses the need to build the cyberinfrastructure for sharing the digital array data. In this White Paper, we detail a broad set of requirements, which involves collecting, managing, presenting, and propagating contextual information essential to assess the quality, understand the content, interpret the scientific implications, and reuse image data in the context of the experimental details. We start by providing an overview of the main lessons learned to date through international community activities, which have recently made considerable progress toward generating community standard practices for imaging Quality Control (QC) and metadata. We then provide a clear set of recommendations for amplifying this work. The driving goal is to address remaining challenges and democratize access to everyday practices and tools for a spectrum of biomedical researchers, regardless of their expertise, access to resources, and geographical location.
△ Less
Submitted 30 August, 2024; v1 submitted 23 January, 2024;
originally announced January 2024.
-
On the average number of zeros of random harmonic polynomials with i.i.d. coefficients: precise asymptotics
Authors:
Erik Lundberg,
Andrew Thomack
Abstract:
Addressing a problem posed by W. Li and A. Wei (2009), we investigate the average number of (complex) zeros of a random harmonic polynomial $p(z) + \overline{q(z)}$ sampled from the Kac ensemble, i.e., where the coefficients are independent identically distributed centered complex Gaussian random variables. We establish a precise asymptotic, showing that when $\text{deg} p = \text{deg} q = n$ tend…
▽ More
Addressing a problem posed by W. Li and A. Wei (2009), we investigate the average number of (complex) zeros of a random harmonic polynomial $p(z) + \overline{q(z)}$ sampled from the Kac ensemble, i.e., where the coefficients are independent identically distributed centered complex Gaussian random variables. We establish a precise asymptotic, showing that when $\text{deg} p = \text{deg} q = n$ tends to infinity the average number of zeros is asymptotic to $\frac{1}{2} n \log n$. We further consider the average number of zeros restricted to various regions in the complex plane leading to interesting comparisons with the classically studied case of analytic Kac polynomials. We also consider deterministic extremal problems for harmonic polynomials with coefficient constraints; using an indirect probabilistic method we show the existence of harmonic polynomials with unimodular coefficients having at least $\frac{2}π n \log n + O(n)$ zeros. We conclude with a list of open problems.
△ Less
Submitted 20 August, 2023;
originally announced August 2023.
-
Multiplane gravitational lenses with an abundance of images
Authors:
Charles R. Keeton,
Erik Lundberg,
Sean Perry
Abstract:
We consider gravitational lensing of a background source by a finite system of point-masses. The problem of determining the maximum possible number of lensed images has been completely resolved in the single-plane setting (where the point masses all reside in a single lens plane), but this problem remains open in the multiplane setting. We construct examples of $K$-plane point-mass gravitational l…
▽ More
We consider gravitational lensing of a background source by a finite system of point-masses. The problem of determining the maximum possible number of lensed images has been completely resolved in the single-plane setting (where the point masses all reside in a single lens plane), but this problem remains open in the multiplane setting. We construct examples of $K$-plane point-mass gravitational lens ensembles that produce $\prod_{i=1}^K (5g_i-5)$ images of a single background source, where $g_i$ is the number of point masses in the $i^\text{th}$ plane. This gives asymptotically (for large $g_i$ with $K$ fixed) $5^K$ times the minimal number of lensed images. Our construction uses Rhie's single-plane examples and a structured parameter-rescaling algorithm to produce preliminary systems of equations with the desired number of solutions. Utilizing the stability principle from differential topology, we then show that the preliminary (nonphysical) examples can be perturbed to produce physically meaningful examples while preserving the number of solutions. We provide numerical simulations illustrating the result of our construction, including the positions of lensed images as well as the structure of the critical curves and caustics. We observe an interesting ``caustic of multiplicity'' phenomenon that occurs in the nonphysical case and has a noticeable effect on the caustic structure in the physically meaningful perturbative case.
△ Less
Submitted 22 February, 2023;
originally announced February 2023.
-
On the valence of logharmonic polynomials
Authors:
Dmitry Khavinson,
Erik Lundberg,
Sean Perry
Abstract:
Investigating a problem posed by W. Hengartner (2000), we study the maximal valence (number of preimages of a prescribed point in the complex plane) of logharmonic polynomials, i.e., complex functions that take the form $f(z) = p(z) \overline{q(z)}$ of a product of an analytic polynomial $p(z)$ of degree $n$ and the complex conjugate of another analytic polynomial $q(z)$ of degree $m$. In the case…
▽ More
Investigating a problem posed by W. Hengartner (2000), we study the maximal valence (number of preimages of a prescribed point in the complex plane) of logharmonic polynomials, i.e., complex functions that take the form $f(z) = p(z) \overline{q(z)}$ of a product of an analytic polynomial $p(z)$ of degree $n$ and the complex conjugate of another analytic polynomial $q(z)$ of degree $m$. In the case $m=1$, we adapt an indirect technique utilizing anti-holomorphic dynamics to show that the valence is at most $3n-1$. This confirms a conjecture of Bshouty and Hengartner (2000). Using a purely algebraic method based on Sylvester resultants, we also prove a general upper bound for the valence showing that for each $n,m \geq 1$ the valence is at most $n^2+m^2$. This improves, for every choice of $n,m \geq 1$, the previously established upper bound $(n+m)^2$ based on Bezout's theorem. We also consider the more general setting of polyanalytic polynomials where we show that this latter result can be extended under a nondegeneracy assumption.
△ Less
Submitted 8 February, 2023;
originally announced February 2023.
-
Inradius of random lemniscates
Authors:
Manjunath Krishnapur,
Erik Lundberg,
Koushik Ramachandran
Abstract:
A classically studied geometric property associated to a complex polynomial $p$ is the inradius (the radius of the largest inscribed disk) of its (filled) lemniscate $Λ:= \{z \in \mathbb{C}:|p(z)| < 1\}$.
In this paper, we study the lemniscate inradius when the defining polynomial $p$ is random, namely, with the zeros of $p$ sampled independently from a compactly supported probability measure…
▽ More
A classically studied geometric property associated to a complex polynomial $p$ is the inradius (the radius of the largest inscribed disk) of its (filled) lemniscate $Λ:= \{z \in \mathbb{C}:|p(z)| < 1\}$.
In this paper, we study the lemniscate inradius when the defining polynomial $p$ is random, namely, with the zeros of $p$ sampled independently from a compactly supported probability measure $μ$. If the negative set of the logarithmic potential $U_μ$ generated by $μ$ is non-empty, then the inradius is bounded from below by a positive constant with overwhelming probability. Moreover, the inradius has a determinstic limit if the negative set of $U_μ$ additionally contains the support of $μ$.
On the other hand, when the zeros are sampled independently and uniformly from the unit circle, then the inradius converges in distribution to a random variable taking values in $(0,1/2)$.
We also consider the characteristic polynomial of a Ginibre random matrix whose lemniscate we show is close to the unit disk with overwhelming probability.
△ Less
Submitted 31 January, 2023;
originally announced January 2023.
-
The valence of harmonic polynomials viewed through the probabilistic lens
Authors:
Erik Lundberg
Abstract:
We prove the existence of complex polynomials $p(z)$ of degree $n$ and $q(z)$ of degree $m<n$ such that the harmonic polynomial $ p(z) + \overline{q(z)}$ has at least $\lceil n \sqrt{m} \rceil$ many zeros. This provides an array of new counterexamples to Wilmshurst's conjecture that the maximum valence of harmonic polynomials $p(z)+\overline{q(z)}$ taken over polynomials $p$ of degree $n$ and $q$…
▽ More
We prove the existence of complex polynomials $p(z)$ of degree $n$ and $q(z)$ of degree $m<n$ such that the harmonic polynomial $ p(z) + \overline{q(z)}$ has at least $\lceil n \sqrt{m} \rceil$ many zeros. This provides an array of new counterexamples to Wilmshurst's conjecture that the maximum valence of harmonic polynomials $p(z)+\overline{q(z)}$ taken over polynomials $p$ of degree $n$ and $q$ of degree $m$ is $m(m-1)+3n-2$. More broadly, these examples show that there does not exist a linear (in $n$) bound on the valence with a uniform (in $m$) growth rate. The proof of this result uses a probabilistic technique based on estimating the average number of zeros of a certain family of random harmonic polynomials.
△ Less
Submitted 31 August, 2023; v1 submitted 3 January, 2022;
originally announced January 2022.
-
The number of limit cycles bifurcating from a randomly perturbed center
Authors:
Manjunath Krishnapur,
Erik Lundberg,
Oanh Nguyen
Abstract:
We consider the average number of limit cycles that bifurcate from a randomly perturbed linear center where the perturbation consists of random (bivariate) polynomials with independent coefficients. This problem reduces, by way of classical perturbation theory of the Poincaré first return map, to a problem on the real zeros of a random \emph{univariate} polynomial…
▽ More
We consider the average number of limit cycles that bifurcate from a randomly perturbed linear center where the perturbation consists of random (bivariate) polynomials with independent coefficients. This problem reduces, by way of classical perturbation theory of the Poincaré first return map, to a problem on the real zeros of a random \emph{univariate} polynomial $\displaystyle f_n(x) = \sum_{m=0}^n c_m ξ_m x^m$ with independent coefficients $ξ_m$ having mean zero, variance 1 and $c_m \sim m^{-1/2}$. This polynomial belongs to the class of {\it generalized Kac polynomials} at the critical regime. We provide asymptotics for the average number of real zeros and answer the question on bifurcating limit cycles. Additionally, we provide the correct order of the mean number of real roots in the subcritical regime.
△ Less
Submitted 19 August, 2022; v1 submitted 10 December, 2021;
originally announced December 2021.
-
New views of old proteins: clarifying the enigmatic proteome
Authors:
Participants in a NIH Workshop on Functional,
Integrative Proteomics,
:,
Kristin E. Burnum Johnson,
Thomas P. Conrads,
Richard R. Drake,
Amy E. Herr,
Ravi Iyengar,
Ryan T. Kelly,
Emma Lundberg,
Michael J. MacCoss,
Alexandra Naba,
Garry P. Nolan,
Pavel A. Pevzner,
Karin D. Rodland,
Salvatore Sechi,
Nikolai Slavov,
Jeffrey M. Spraggins,
Jennifer E. Van Eyk,
Marc Vidal,
Christine Vogel,
David R. Walt,
Neil L. Kelleher
Abstract:
All human diseases involve proteins, yet our current tools to characterize and quantify them are limited. To better elucidate proteins across space, time, and molecular composition, we provide provocative projections for technologies to meet the challenges that protein biology presents. With a broad perspective, we discuss grand opportunities to transition the science of proteomics into a more pro…
▽ More
All human diseases involve proteins, yet our current tools to characterize and quantify them are limited. To better elucidate proteins across space, time, and molecular composition, we provide provocative projections for technologies to meet the challenges that protein biology presents. With a broad perspective, we discuss grand opportunities to transition the science of proteomics into a more propulsive enterprise. Extrapolating recent trends, we offer potential futures for a next generation of disruptive approaches to define, quantify and visualize the multiple dimensions of the proteome, thereby transforming our understanding and interactions with human disease in the coming decade.
△ Less
Submitted 17 August, 2021;
originally announced August 2021.
-
Spatial mapping of protein composition and tissue organization: a primer for multiplexed antibody-based imaging
Authors:
John W. Hickey,
Elizabeth K. Neumann,
Andrea J. Radtke,
Jeannie M. Camarillo,
Rebecca T. Beuschel,
Alexandre Albanese,
Elizabeth McDonough,
Julia Hatler,
Anne E. Wiblin,
Jeremy Fisher,
Josh Croteau,
Eliza C. Small,
Anup Sood,
Richard M. Caprioli,
R. Michael Angelo,
Garry P. Nolan,
Kwanghun Chung,
Stephen M. Hewitt,
Ronald N. Germain,
Jeffrey M. Spraggins,
Emma Lundberg,
Michael P. Snyder,
Neil L. Kelleher,
Sinem K. Saka
Abstract:
Tissues and organs are composed of distinct cell types that must operate in concert to perform physiological functions. Efforts to create high-dimensional biomarker catalogs of these cells are largely based on transcriptomic single-cell approaches that lack the spatial context required to understand critical cellular communication and correlated structural organization. To probe in situ biology wi…
▽ More
Tissues and organs are composed of distinct cell types that must operate in concert to perform physiological functions. Efforts to create high-dimensional biomarker catalogs of these cells are largely based on transcriptomic single-cell approaches that lack the spatial context required to understand critical cellular communication and correlated structural organization. To probe in situ biology with sufficient coverage depth, several multiplexed protein imaging methods have recently been developed. Though these antibody-based technologies differ in strategy and mode of immunolabeling and detection tags, they commonly utilize antibodies directed against protein biomarkers to provide detailed spatial and functional maps of complex tissues. As these promising antibody-based multiplexing approaches become more widely adopted, new frameworks and considerations are critical for training future users, generating molecular tools, validating antibody panels, and harmonizing datasets. In this perspective, we provide essential resources and key considerations for obtaining robust and reproducible multiplexed antibody-based imaging data compiling specialized knowledge from domain experts and technology developers.
△ Less
Submitted 16 July, 2021;
originally announced July 2021.
-
On the Number of Equilibria Balancing Newtonian Point Masses with a Central Force
Authors:
Nickolas Arustamyan,
Christopher Cox,
Erik Lundberg,
Sean Perry,
Zvi Rosen
Abstract:
We consider the critical points (equilibria) of a planar potential generated by $n$ Newtonian point masses augmented with a quadratic term (such as arises from a centrifugal effect). Particular cases of this problem have been considered previously in studies of the circular restricted $n$-body problem. We show that the number of equilibria is finite for a generic set of parameters, and we establis…
▽ More
We consider the critical points (equilibria) of a planar potential generated by $n$ Newtonian point masses augmented with a quadratic term (such as arises from a centrifugal effect). Particular cases of this problem have been considered previously in studies of the circular restricted $n$-body problem. We show that the number of equilibria is finite for a generic set of parameters, and we establish estimates for the number of equilibria. We prove that the number of equilibria is bounded below by $n+1$, and we provide examples to show that this lower bound is sharp. We prove an upper bound on the number of equilibria that grows exponentially in $n$. In order to establish a lower bound on the maximum number of equilibria, we analyze a class of examples, referred to as ``ring configurations'', consisting of $n-1$ equal masses positioned at vertices of a regular polygon with an additional mass located at the center. Previous numerical observations indicate that these configurations can produce as many as $5n-5$ equilibria. We verify analytically that the ring configuration has at least $5n-5$ equilibria when the central mass is sufficiently small. We conjecture that the maximum number of equilibria grows linearly with the number of point masses. We also discuss some mathematical similarities to other equilibrium problems in mathematical physics, namely, Maxwell's problem from electrostatics and the image counting problem from gravitational lensing.
△ Less
Submitted 21 June, 2021;
originally announced June 2021.
-
Knowledge-Based Strategies for Multi-Agent Teams Playing Against Nature
Authors:
Dilian Gurov,
Valentin Goranko,
Edvin Lundberg
Abstract:
We study teams of agents that play against Nature towards achieving a common objective. The agents are assumed to have imperfect information due to partial observability, and have no communication during the play of the game. We propose a natural notion of higher-order knowledge of agents. Based on this notion, we define a class of knowledge-based strategies, and consider the problem of synthesis…
▽ More
We study teams of agents that play against Nature towards achieving a common objective. The agents are assumed to have imperfect information due to partial observability, and have no communication during the play of the game. We propose a natural notion of higher-order knowledge of agents. Based on this notion, we define a class of knowledge-based strategies, and consider the problem of synthesis of strategies of this class. We introduce a multi-agent extension, MKBSC, of the well-known Knowledge-Based Subset Construction applied to such games. Its iterative applications turn out to compute higher-order knowledge of the agents. We show how the MKBSC can be used for the design of knowledge-based strategy profiles and investigate the transfer of existence of such strategies between the original game and in the iterated applications of the MKBSC, under some natural assumptions. We also relate and compare the "intensional" view on knowledge-based strategies based on explicit knowledge representation and update, with the "extensional" view on finite memory strategies based on finite transducers and show that, in a certain sense, these are equivalent.
△ Less
Submitted 28 December, 2021; v1 submitted 29 December, 2020;
originally announced December 2020.
-
Limit cycle enumeration in random vector fields
Authors:
Erik Lundberg
Abstract:
We study the number and distribution of the limit cycles of a planar vector field whose component functions are random polynomials. We prove a lower bound on the average number of limit cycles when the random polynomials are sampled from the Kostlan-Shub-Smale ensemble. Investigating a problem introduced by Brudnyi [Annals of Mathematics (2001)] we also consider a special local setting of counting…
▽ More
We study the number and distribution of the limit cycles of a planar vector field whose component functions are random polynomials. We prove a lower bound on the average number of limit cycles when the random polynomials are sampled from the Kostlan-Shub-Smale ensemble. Investigating a problem introduced by Brudnyi [Annals of Mathematics (2001)] we also consider a special local setting of counting limit cycles near a randomly perturbed center focus, and when the perturbation has i.i.d. coefficients, we prove a limit law showing that the number of limit cycles situated within a disk of radius less than unity converges almost surely to the number of real zeros of a logarithmically-correlated random univariate power series. We also consider infinitesimal perturbations where we obtain precise asymptotics on the global average count of limit cycles for a family of models. The proofs of these results use novel combinations of techniques from dynamical systems and random analytic functions.
△ Less
Submitted 8 June, 2023; v1 submitted 1 July, 2020;
originally announced July 2020.
-
A note on arclength null quadrature domains
Authors:
Dmitry Khavinson,
Erik Lundberg
Abstract:
We prove the existence of a roof function for arclength null quadrature domains having finitely many boundary components. This bridges a gap toward classification of arclength null quadrature domains by removing an a priori assumption from previous classification results.
We prove the existence of a roof function for arclength null quadrature domains having finitely many boundary components. This bridges a gap toward classification of arclength null quadrature domains by removing an a priori assumption from previous classification results.
△ Less
Submitted 4 September, 2021; v1 submitted 19 May, 2020;
originally announced May 2020.
-
Asymptotic enumeration of lonesum matrices
Authors:
Jessica Khera,
Erik Lundberg,
Stephen Melczer
Abstract:
We provide bivariate asymptotics for the poly-Bernoulli numbers, a combinatorial array that enumerates lonesum matrices, using the methods of Analytic Combinatorics in Several Variables (ACSV). For the diagonal asymptotic (i.e., for the special case of square lonesum matrices) we present an alternative proof based on Parseval's identity. In addition, we provide an application in Algebraic Statisti…
▽ More
We provide bivariate asymptotics for the poly-Bernoulli numbers, a combinatorial array that enumerates lonesum matrices, using the methods of Analytic Combinatorics in Several Variables (ACSV). For the diagonal asymptotic (i.e., for the special case of square lonesum matrices) we present an alternative proof based on Parseval's identity. In addition, we provide an application in Algebraic Statistics on the asymptotic ML-degree of the bivariate multinomial missing data problem, and we strengthen an existing result on asymptotic enumeration of permutations having a specified excedance set.
△ Less
Submitted 6 October, 2020; v1 submitted 18 December, 2019;
originally announced December 2019.
-
Homotopy Types of Random Cubical Complexes
Authors:
Kenneth Dowling,
Erik Lundberg
Abstract:
We study the topology of a random cubical complex associated to Bernoulli site percolation on a cubical grid. We begin by establishing a limit law for homotopy types. More precisely, looking within an expanding window, we define a sequence of normalized counting measures (counting connected components according to homotopy type), and we show that this sequence of random probability measures conver…
▽ More
We study the topology of a random cubical complex associated to Bernoulli site percolation on a cubical grid. We begin by establishing a limit law for homotopy types. More precisely, looking within an expanding window, we define a sequence of normalized counting measures (counting connected components according to homotopy type), and we show that this sequence of random probability measures converges in probability to a deterministic probability measure. We then investigate the dependence of the limiting homotopy measure on the coloring probability $p$, and our results show a qualitative change in the homotopy measure as $p$ crosses the percolation threshold $p=p_c$. Specializing to the case of $d=2$ dimensions, we also present empirical results that raise further questions on the $p$-dependence of the limiting homotopy measure.
△ Less
Submitted 12 September, 2021; v1 submitted 28 October, 2019;
originally announced October 2019.
-
A note on the critical points of the localization landscape
Authors:
Erik Lundberg,
Koushik Ramachandran
Abstract:
Let $Ω\subset\mathbb{C}$ be a bounded domain. In this note, we use complex variable methods to study the number of critical points of the function $v=v_Ω$ that solves the elliptic problem $Δv = -2$ in $Ω,$ with boundary values $v=0$ on $\partialΩ.$ This problem has a classical flavor but is especially motivated by recent studies on localization of eigenfunctions. We provide an upper bound on the n…
▽ More
Let $Ω\subset\mathbb{C}$ be a bounded domain. In this note, we use complex variable methods to study the number of critical points of the function $v=v_Ω$ that solves the elliptic problem $Δv = -2$ in $Ω,$ with boundary values $v=0$ on $\partialΩ.$ This problem has a classical flavor but is especially motivated by recent studies on localization of eigenfunctions. We provide an upper bound on the number of critical points of $v$ when $Ω$ belongs to a special class of domains in the plane, namely, domains for which the boundary $\partialΩ$ is contained in $\{z:|z|^2 = f(z) + \overline{f(z)}\},$ where $f'(z)$ is a rational function. We furnish examples of domains where this bound is attained. We also prove a bound on the number of critical points in the case when $Ω$ is a quadrature domain, and conclude the note by stating some open problems and conjectures.
△ Less
Submitted 29 April, 2021; v1 submitted 19 July, 2019;
originally announced July 2019.
-
ImJoy: an open-source computational platform for the deep learning era
Authors:
Wei Ouyang,
Florian Mueller,
Martin Hjelmare,
Emma Lundberg,
Christophe Zimmer
Abstract:
Deep learning methods have shown extraordinary potential for analyzing very diverse biomedical data, but their dissemination beyond developers is hindered by important computational hurdles. We introduce ImJoy (https://imjoy.io/), a flexible and open-source browser-based platform designed to facilitate widespread reuse of deep learning solutions in biomedical research. We highlight ImJoy's main fe…
▽ More
Deep learning methods have shown extraordinary potential for analyzing very diverse biomedical data, but their dissemination beyond developers is hindered by important computational hurdles. We introduce ImJoy (https://imjoy.io/), a flexible and open-source browser-based platform designed to facilitate widespread reuse of deep learning solutions in biomedical research. We highlight ImJoy's main features and illustrate its functionalities with deep learning plugins for mobile and interactive image analysis and genomics.
△ Less
Submitted 30 May, 2019;
originally announced May 2019.
-
Topologies of random geometric complexes on Riemannian manifolds in the thermodynamic limit
Authors:
Antonio Auffinger,
Antonio Lerario,
Erik Lundberg
Abstract:
We investigate the topologies of random geometric complexes built over random points sampled on Riemannian manifolds in the so-called "thermodynamic" regime. We prove the existence of universal limit laws for the topologies; namely, the random normalized counting measure of connected components (counted according to homotopy type) is shown to converge in probability to a deterministic probability…
▽ More
We investigate the topologies of random geometric complexes built over random points sampled on Riemannian manifolds in the so-called "thermodynamic" regime. We prove the existence of universal limit laws for the topologies; namely, the random normalized counting measure of connected components (counted according to homotopy type) is shown to converge in probability to a deterministic probability measure. Moreover, we show that the support of the deterministic limiting measure equals the set of all homotopy types for Euclidean geometric complexes of the same dimension as the manifold.
△ Less
Submitted 19 November, 2020; v1 submitted 21 December, 2018;
originally announced December 2018.
-
EGBTER: Capturing degree distribution, clustering coefficients, and community structure in a single random graph model
Authors:
Omar El-daghar,
Erik Lundberg,
Robert A. Bridges
Abstract:
Random graph models are important constructs for data analytic applications as well as pure mathematical developments, as they provide capabilities for network synthesis and principled analysis. Several models have been developed with the aim of faithfully preserving important graph metrics and substructures. With the goal of capturing degree distribution, clustering coefficient, and communities i…
▽ More
Random graph models are important constructs for data analytic applications as well as pure mathematical developments, as they provide capabilities for network synthesis and principled analysis. Several models have been developed with the aim of faithfully preserving important graph metrics and substructures. With the goal of capturing degree distribution, clustering coefficient, and communities in a single random graph model, we propose a new model to address shortcomings in a progression of network modeling capabilities. The Block Two-Level Erd{ő}s-R{é}nyi (BTER) model of Seshadhri et al., designed to allow prescription of expected degree and clustering coefficient distributions, neglects community modeling, while the Generalized BTER (GBTER) model of Bridges et al., designed to add community modeling capabilities to BTER, struggles to faithfully represent all three characteristics simultaneously. In this work, we fit BTER and two GBTER configurations to several real-world networks and compare the results with that of our new model, the Extended GBTER (EGBTER) model. Our results support that EBGTER adds a community-modeling flexibility to BTER, while retaining a satisfactory level of accuracy in terms of degree and clustering coefficient. Our insights and empirical testing of previous models as well as the new model are novel contributions to the literature.
△ Less
Submitted 3 August, 2018;
originally announced August 2018.
-
The lemniscate tree of a random polynomial
Authors:
Michael Epstein,
Boris Hanin,
Erik Lundberg
Abstract:
To each generic complex polynomial $p(z)$ there is associated a labeled binary tree (here referred to as a "lemniscate tree") that encodes the topological type of the graph of $|p(z)|$. The branching structure of the lemniscate tree is determined by the configuration (i.e., arrangement in the plane) of the singular components of those level sets $|p(z)|=t$ passing through a critical point.
In th…
▽ More
To each generic complex polynomial $p(z)$ there is associated a labeled binary tree (here referred to as a "lemniscate tree") that encodes the topological type of the graph of $|p(z)|$. The branching structure of the lemniscate tree is determined by the configuration (i.e., arrangement in the plane) of the singular components of those level sets $|p(z)|=t$ passing through a critical point.
In this paper, we address the question "How many branches appear in a typical lemniscate tree?" We answer this question first for a lemniscate tree sampled uniformly from the combinatorial class and second for the lemniscate tree arising from a random polynomial generated by i.i.d. zeros. From a more general perspective, these results take a first step toward a probabilistic treatment (within a specialized setting) of Arnold's program of enumerating algebraic Morse functions.
△ Less
Submitted 1 June, 2018;
originally announced June 2018.
-
The arc length of a random lemniscate
Authors:
Erik Lundberg,
Koushik Ramachandran
Abstract:
A polynomial lemniscate is a curve in the complex plane defined by $\{z \in \mathbb{C}:|p(z)|=t\}$. Erdös, Herzog, and Piranian posed the extremal problem of determining the maximum length of a lemniscate $Λ=\{ z \in \mathbb{C}:|p(z)|=1\}$ when $p$ is a monic polynomial of degree $n$. In this paper, we study the length and topology of a random lemniscate whose defining polynomial has independent G…
▽ More
A polynomial lemniscate is a curve in the complex plane defined by $\{z \in \mathbb{C}:|p(z)|=t\}$. Erdös, Herzog, and Piranian posed the extremal problem of determining the maximum length of a lemniscate $Λ=\{ z \in \mathbb{C}:|p(z)|=1\}$ when $p$ is a monic polynomial of degree $n$. In this paper, we study the length and topology of a random lemniscate whose defining polynomial has independent Gaussian coefficients. In the special case of the Kac ensemble we show that the length approaches a nonzero constant as $n \rightarrow \infty$. We also show that the average number of connected components is asymptotically $n$, and we observe a positive probability (independent of $n$) of a giant component occurring.
△ Less
Submitted 1 February, 2017; v1 submitted 31 October, 2016;
originally announced October 2016.
-
Random fields and the enumerative geometry of lines on real and complex hypersurfaces
Authors:
Saugata Basu,
Antonio Lerario,
Erik Lundberg,
Chris Peterson
Abstract:
We derive a formula expressing the average number $E_n$ of real lines on a random hypersurface of degree $2n-3$ in $\mathbb{R}\textrm{P}^n$ in terms of the expected modulus of the determinant of a special random matrix. In the case $n=3$ we prove that the average number of real lines on a random cubic surface in $\mathbb{R}\textrm{P}^3$ equals: $$E_3=6\sqrt{2}-3.$$ Our technique can also be used t…
▽ More
We derive a formula expressing the average number $E_n$ of real lines on a random hypersurface of degree $2n-3$ in $\mathbb{R}\textrm{P}^n$ in terms of the expected modulus of the determinant of a special random matrix. In the case $n=3$ we prove that the average number of real lines on a random cubic surface in $\mathbb{R}\textrm{P}^3$ equals: $$E_3=6\sqrt{2}-3.$$ Our technique can also be used to express the number $C_n$ of complex lines on a generic hypersurface of degree $2n-3$ in $\mathbb{C}\textrm{P}^n$ in terms of the determinant of a random Hermitian matrix. As a special case we obtain a new proof of the classical statement $C_3=27.$
We determine, at the logarithmic scale, the asymptotic of the quantity $E_n$, by relating it to $C_n$ (whose asymptotic has been recently computed D. Zagier). Specifically we prove that: $$\lim_{n\to \infty}\frac{\log E_n}{\log C_n}=\frac{1}{2}.$$
Finally we show that this approach can be used to compute the number $R_n=(2n-3)!!$ of real lines, counted with their intrinsic signs, on a generic real hypersurface of degree $2n-3$ in $\mathbb{R}\textrm{P}^n$.
△ Less
Submitted 8 November, 2016; v1 submitted 4 October, 2016;
originally announced October 2016.
-
A four--dimensional Neumann ovaloid
Authors:
Lavi Karp,
Erik Lundberg
Abstract:
What is the shape of a uniformly massive object that generates a gravitational potential equivalent to that of two equal point-masses? If the weight of each point-mass is sufficiently small compared to the distance between the points then the answer is a pair of balls of equal radius, one centered at each of the two points, but otherwise it is a certain domain of revolution about the axis passing…
▽ More
What is the shape of a uniformly massive object that generates a gravitational potential equivalent to that of two equal point-masses? If the weight of each point-mass is sufficiently small compared to the distance between the points then the answer is a pair of balls of equal radius, one centered at each of the two points, but otherwise it is a certain domain of revolution about the axis passing through the two points. The existence and uniqueness of such a domain is known, but an explicit parameterization is known only in the plane where the region is referred to as a Neumann oval. We construct a four-dimensional "Neumann ovaloid", solving explicitly this inverse potential problem.
△ Less
Submitted 25 September, 2016;
originally announced September 2016.
-
The Bergman analytic content of planar domains
Authors:
Matthew Fleeman,
Erik Lundberg
Abstract:
Given a planar domain $Ω$, the Bergman analytic content measures the $L^{2}(Ω)$-distance between $\bar{z}$ and the Bergman space $A^{2}(Ω)$. We compute the Bergman analytic content of simply-connected quadrature domains with quadrature formula supported at one point, and we also determine the function $f \in A^2(Ω)$ that best approximates $\bar{z}$. We show that, for simply-connected domains, the…
▽ More
Given a planar domain $Ω$, the Bergman analytic content measures the $L^{2}(Ω)$-distance between $\bar{z}$ and the Bergman space $A^{2}(Ω)$. We compute the Bergman analytic content of simply-connected quadrature domains with quadrature formula supported at one point, and we also determine the function $f \in A^2(Ω)$ that best approximates $\bar{z}$. We show that, for simply-connected domains, the square of Bergman analytic content is equivalent to torsional rigidity from classical elasticity theory, while for multiply-connected domains these two domain constants are not equivalent in general.
△ Less
Submitted 11 February, 2016;
originally announced February 2016.
-
Dirichlet's problem with entire data posed on an ellipsoidal cylinder
Authors:
Dmitry Khavinson,
Erik Lundberg,
Hermann Render
Abstract:
We consider the Dirichlet problem in an ellipsoidal cylinder when the data function is entire. Under an additional assumption that the order of the data function is less than one, we show that there is a solution that extends as an entire harmonic function.
We consider the Dirichlet problem in an ellipsoidal cylinder when the data function is entire. Under an additional assumption that the order of the data function is less than one, we show that there is a solution that extends as an entire harmonic function.
△ Less
Submitted 8 February, 2016;
originally announced February 2016.
-
The Dirichlet problem for the slab with entire data and a difference equation for harmonic functions
Authors:
Dmitry Khavinson,
Erik Lundberg,
Hermann Render
Abstract:
It is shown that the Dirichlet problem for the slab $(a,b) \times \mathbb{R}^{d}$ with entire boundary data has an entire solution. The proof is based on a generalized Schwarz reflection principle. Moreover, it is shown that for a given entire harmonic function $g$ the inhomogeneous difference equation $h\left( t+1,y\right) -h\left(t,y\right) =g\left(t,y\right)$ has an entire harmonic solution…
▽ More
It is shown that the Dirichlet problem for the slab $(a,b) \times \mathbb{R}^{d}$ with entire boundary data has an entire solution. The proof is based on a generalized Schwarz reflection principle. Moreover, it is shown that for a given entire harmonic function $g$ the inhomogeneous difference equation $h\left( t+1,y\right) -h\left(t,y\right) =g\left(t,y\right)$ has an entire harmonic solution $h$.
△ Less
Submitted 4 February, 2016;
originally announced February 2016.
-
On the geometry of random lemniscates
Authors:
Antonio Lerario,
Erik Lundberg
Abstract:
We investigate the geometry of a random rational lemniscate $Γ$, the level set $\{|r(z)|=1\}$ on the Riemann sphere of the modulus of a random rational function $r$. We assign a probability distribution to the space of rational functions $r=p/q$ of degree $n$ by sampling $p$ and $q$ independently from the complex Kostlan ensemble of random polynomials of degree $n$.
We prove that the average \em…
▽ More
We investigate the geometry of a random rational lemniscate $Γ$, the level set $\{|r(z)|=1\}$ on the Riemann sphere of the modulus of a random rational function $r$. We assign a probability distribution to the space of rational functions $r=p/q$ of degree $n$ by sampling $p$ and $q$ independently from the complex Kostlan ensemble of random polynomials of degree $n$.
We prove that the average \emph{spherical length} of $Γ$ is $\frac{π^2}{2} \sqrt{n},$ which is proportional to the square root of the maximal spherical length. We also provide an asymptotic for the average number of points on the curve that are tangent to one of the meridians on the Riemann sphere (i.e. tangent to one of the radial directions in the plane).
Concerning the topology of $Γ$, on a local scale, we prove that for every disk $D$ of radius $O(n^{-1/2})$ in the Riemann sphere and any \emph{arrangement} (i.e. embedding) of finitely many circles $A\subset D$ there is a positive probability (independent of $n$) that $(D,Γ\cap D)$ is isotopic to $( D,A)$. (A local random version of Hilbert's Sixteenth Problem restricted to lemniscates.) Corollary: the average number of connected components of $Γ$ increases linearly (the maximum rate possible according to a deterministic upper bound).
△ Less
Submitted 10 January, 2016;
originally announced January 2016.
-
On the zeros of random harmonic polynomials: the truncated model
Authors:
Antonio Lerario,
Erik Lundberg
Abstract:
A probabilistic approach to the study of the number of zeros of complex harmonic polynomials was initiated by W. Li and A. Wei (2009), who derived a Kac-Rice type formula for the expected number of zeros of random harmonic polynomials with independent Gaussian coefficients. They also provided asymptotics for a complex version of the Kostlan ensemble. Here we determine asymptotics for the alternati…
▽ More
A probabilistic approach to the study of the number of zeros of complex harmonic polynomials was initiated by W. Li and A. Wei (2009), who derived a Kac-Rice type formula for the expected number of zeros of random harmonic polynomials with independent Gaussian coefficients. They also provided asymptotics for a complex version of the Kostlan ensemble. Here we determine asymptotics for the alternative truncated model that was recently proposed by J. Hauenstein, D. Mehta, and the authors. Our results confirm (and sharpen) their (3/2)-powerlaw conjecture that had been formulated on the basis of computer experiments.
△ Less
Submitted 12 May, 2016; v1 submitted 3 July, 2015;
originally announced July 2015.
-
Experiments on the zeros of harmonic polynomials using certified counting
Authors:
Jonathan D. Hauenstein,
Antonio Lerario,
Erik Lundberg,
Dhagash Mehta
Abstract:
Motivated by Wilmshurst's conjecture, we investigate the zeros of harmonic polynomials. We utilize a certified counting approach which is a combination of two methods from numerical algebraic geometry: numerical polynomial homotopy continuation to compute a numerical approximation of each zero and Smale's alpha-theory to certify the results. Using this approach, we provide new examples of harmonic…
▽ More
Motivated by Wilmshurst's conjecture, we investigate the zeros of harmonic polynomials. We utilize a certified counting approach which is a combination of two methods from numerical algebraic geometry: numerical polynomial homotopy continuation to compute a numerical approximation of each zero and Smale's alpha-theory to certify the results. Using this approach, we provide new examples of harmonic polynomials having the most extreme number of zeros known so far; we also study the mean and variance of the number of zeros of random harmonic polynomials.
△ Less
Submitted 20 June, 2014;
originally announced June 2014.
-
Quasi-exceptional domains
Authors:
Alexandre Eremenko,
Erik Lundberg
Abstract:
Exceptional domains are domains on which there exists a positive harmonic function, zero on the boundary and such that the normal derivative on the boundary is constant. Recent results classify exceptional domains as belonging to either a certain one-parameter family of simply periodic domains or one of its scaling limits.
We introduce quasi-exceptional domains by allowing the boundary values to…
▽ More
Exceptional domains are domains on which there exists a positive harmonic function, zero on the boundary and such that the normal derivative on the boundary is constant. Recent results classify exceptional domains as belonging to either a certain one-parameter family of simply periodic domains or one of its scaling limits.
We introduce quasi-exceptional domains by allowing the boundary values to be different constants on each boundary component. This relaxed definition retains the interesting property of being an \emph{arclength quadrature domain}, and also preserves the connection to the hollow vortex problem in fluid dynamics. We give a partial classification of such domains in terms of certain Abelian differentials. We also provide a new two-parameter family of periodic quasi-exceptional domains. These examples generalize the hollow vortex array found by Baker, Saffman, and Sheffield (1976). A degeneration of regions of this family provide doubly-connected examples.
△ Less
Submitted 29 May, 2014;
originally announced May 2014.
-
On the number of connected components of random algebraic hypersurfaces
Authors:
Yan Fyodorov,
Antonio Lerario,
Erik Lundberg
Abstract:
We study the expectation of the number of components $b_0(X)$ of a random algebraic hypersurface $X$ defined by the zero set in projective space $\mathbb{R}P^n$ of a random homogeneous polynomial $f$ of degree $d$. Specifically, we consider "invariant ensembles", that is Gaussian ensembles of polynomials that are invariant under an orthogonal change of variables.
The classification due to E. Kos…
▽ More
We study the expectation of the number of components $b_0(X)$ of a random algebraic hypersurface $X$ defined by the zero set in projective space $\mathbb{R}P^n$ of a random homogeneous polynomial $f$ of degree $d$. Specifically, we consider "invariant ensembles", that is Gaussian ensembles of polynomials that are invariant under an orthogonal change of variables.
The classification due to E. Kostlan shows that specifying an invariant ensemble is equivalent to assigning a weight to each eigenspace of the spherical Laplacian. Fixing $n$, we consider a family of invariant ensembles (choice of eigenspace weights) depending on the degree $d$. Under a rescaling assumption on the eigenspace weights (as $d \rightarrow \infty$), we prove that the order of growth of $\mathbb{E} b_0(X)$ satisfies: $$\mathbb{E} b_{0}(X)=Θ\left(\left[ \mathbb{E} b_0(X\cap \mathbb{R}P^1) \right]^{n} \right). $$ This relates the average number of components of $X$ to the classical problem of M. Kac (1943) on the number of zeros of the random univariate polynomial $f|_{\mathbb{R}P^1}.$
The proof requires an upper bound for $\mathbb{E} b_0(X)$, which we obtain by counting extrema using Random Matrix Theory methods from recent work of the first author, and it also requires a lower bound, which we obtain by a modification of the barrier method. We also provide a quantitative upper bound for the implied constant in the above asymptotic; for the real Fubini-Study model these estimates reveal super-exponential decay of the leading coefficient (in $d$) of $\mathbb{E} b_0(X)$ (as $n \rightarrow \infty$).
△ Less
Submitted 27 June, 2015; v1 submitted 21 April, 2014;
originally announced April 2014.
-
Asymptotics of the Extremal Excedance Set Statistic
Authors:
Rodrigo Ferraz de Andrade,
Erik Lundberg,
Brendan Nagle
Abstract:
Answering a question of Clark and Ehrenborg (2010), we determine asymptotics for the number of permutations of size n that admit the most common excedance set. In fact, we provide a more general bivariate asymptotic using the multivariate asymptotic methods of R. Pemantle and M. C. Wilson. We also consider two applications of our main result. First, we determine asymptotics on the number of permut…
▽ More
Answering a question of Clark and Ehrenborg (2010), we determine asymptotics for the number of permutations of size n that admit the most common excedance set. In fact, we provide a more general bivariate asymptotic using the multivariate asymptotic methods of R. Pemantle and M. C. Wilson. We also consider two applications of our main result. First, we determine asymptotics on the number of permutations of size n which simultaneously avoid the generalized patterns 21-34 and 34-21. Second, we determine asymptotics on the number of n-cycles that admit no stretching pairs.
△ Less
Submitted 26 May, 2016; v1 submitted 4 March, 2014;
originally announced March 2014.
-
On the growth of solutions to the minimal surface equation over domains containing a halfplane
Authors:
Erik Lundberg,
Allen Weitsman
Abstract:
We consider minimal graphs u(x,y)>0 over unbounded domains D (with u vanishing on the boundary of D). Assuming D contains a sector properly containing a halfplane, we obtain estimates on growth and provide examples illustrating a range of growth.
We consider minimal graphs u(x,y)>0 over unbounded domains D (with u vanishing on the boundary of D). Assuming D contains a sector properly containing a halfplane, we obtain estimates on growth and provide examples illustrating a range of growth.
△ Less
Submitted 13 November, 2013;
originally announced November 2013.
-
Gap probabilities and applications to geometry and random topology
Authors:
Antonio Lerario,
Erik Lundberg
Abstract:
We give an exact formula for the value of the derivative at zero of the gap probability in finite n x n Gaussian ensembles. As n goes to infinity our computation provides an asymptotic (with an explicit constant) of the order n^(1/2).
As a first application, we consider the set of n x n (Real, Complex or Quaternionic) Hermitian matrices with Frobenius norm one and determinant zero. We give an ex…
▽ More
We give an exact formula for the value of the derivative at zero of the gap probability in finite n x n Gaussian ensembles. As n goes to infinity our computation provides an asymptotic (with an explicit constant) of the order n^(1/2).
As a first application, we consider the set of n x n (Real, Complex or Quaternionic) Hermitian matrices with Frobenius norm one and determinant zero. We give an exact formula for the intrinsic volume of this set and as n goes to infinity its asymptotic (with an explicit constant) is of the order n^(1/2).
As a second application we consider the problem of computing Betti numbers of an intersection of k random Kostlan quadrics in RP^n. We show that the i-th Betti number is asymptotically expected to be one (for i sufficiently away from n/2).
In the case k=2 the the sum of all Betti numbers was recently shown by the first author to equal n+o(n). Here we sharpen this asymptotic proving an asymptotic with two orders of precision and explicit constants.
△ Less
Submitted 22 September, 2013;
originally announced September 2013.
-
Electrostatic skeletons
Authors:
Alexandre Eremenko,
Erik Lundberg,
Koushik Ramachandran
Abstract:
Let u be the equilibrium potential of a compact set K. An electrostatic skeleton of K is a positive measure whose closed support has connected complement and no interior, and whose potential is equal to u near infinity. We prove the existence of an electrostatic skeleton for any simplex.
Let u be the equilibrium potential of a compact set K. An electrostatic skeleton of K is a positive measure whose closed support has connected complement and no interior, and whose potential is equal to u near infinity. We prove the existence of an electrostatic skeleton for any simplex.
△ Less
Submitted 17 November, 2013; v1 submitted 21 September, 2013;
originally announced September 2013.
-
A tale of ellipsoids in potential theory
Authors:
Dmitry Khavinson,
Erik Lundberg
Abstract:
Ellipsoids possess several beautiful properties associated with classical potential theory. Some of them are well known, and some have been forgotten. In this article we hope to bring a few of the "lost" pieces of classical mathematics back to the limelight.
Ellipsoids possess several beautiful properties associated with classical potential theory. Some of them are well known, and some have been forgotten. In this article we hope to bring a few of the "lost" pieces of classical mathematics back to the limelight.
△ Less
Submitted 17 March, 2014; v1 submitted 9 September, 2013;
originally announced September 2013.
-
Remarks on Wilmshurst's theorem
Authors:
Seung-Yeop Lee,
Antonio Lerario,
Erik Lundberg
Abstract:
We demonstrate counterexamples to Wilmshurst's conjecture on the valence of harmonic polynomials in the plane, and we conjecture a bound that is linear in the analytic degree for each fixed anti-analytic degree. Then we initiate a discussion of Wilmshurt's theorem in more than two dimensions, showing that if the zero set of a polynomial harmonic field is bounded then it must have codimension at le…
▽ More
We demonstrate counterexamples to Wilmshurst's conjecture on the valence of harmonic polynomials in the plane, and we conjecture a bound that is linear in the analytic degree for each fixed anti-analytic degree. Then we initiate a discussion of Wilmshurt's theorem in more than two dimensions, showing that if the zero set of a polynomial harmonic field is bounded then it must have codimension at least two. Examples are provided to show that this conclusion cannot be improved.
△ Less
Submitted 29 August, 2013;
originally announced August 2013.
-
A solution to Sheil-Small's harmonic mapping problem for Jordan polygons
Authors:
Daoud Bshouty,
Erik Lundberg,
Allen Weitsman
Abstract:
The problem of mapping the interior of a Jordan polygon univalently by the Poisson integral of a step function was posed by T. Sheil-Small (1989). We describe a simple solution using "ear clipping" from computational geometry.
The problem of mapping the interior of a Jordan polygon univalently by the Poisson integral of a step function was posed by T. Sheil-Small (1989). We describe a simple solution using "ear clipping" from computational geometry.
△ Less
Submitted 22 February, 2013;
originally announced February 2013.
-
Statistics on Hilbert's Sixteenth Problem
Authors:
Antonio Lerario,
Erik Lundberg
Abstract:
We study the statistics of the number of connected components and the volume of a random real algebraic hypersurface in RP^n defined by a Real Bombieri-Weyl distributed homogeneous polynomial of degree d. We prove that the expectation of the number of connected components of such hypersurface has order d^n, the asymptotic being in d for n fixed. We do not restrict ourselves to the random homogeneo…
▽ More
We study the statistics of the number of connected components and the volume of a random real algebraic hypersurface in RP^n defined by a Real Bombieri-Weyl distributed homogeneous polynomial of degree d. We prove that the expectation of the number of connected components of such hypersurface has order d^n, the asymptotic being in d for n fixed. We do not restrict ourselves to the random homogeneous case and we consider more generally random polynomials belonging to a window of eigenspaces of the laplacian on the sphere S^n, proving that the same asymptotic holds. As for the volume, we prove its expectation is of order d. Both these behaviors exhibit expectation of maximal order in light of Milnor's bound and the a priori bound for the volume.
△ Less
Submitted 21 January, 2013; v1 submitted 16 December, 2012;
originally announced December 2012.
-
Self-commutators of Toeplitz operators and isoperimetric inequalities
Authors:
Steven R. Bell,
Timothy Ferguson,
Erik Lundberg
Abstract:
For a hyponormal operator, C. R. Putnam's inequality gives an upper bound on the norm of its self-commutator. In the special case of a Toeplitz operator with analytic symbol in the Smirnov space of a domain, there is also a geometric lower bound shown by D. Khavinson (1985) that when combined with Putnam's inequality implies the classical isoperimetric inequality. For a nontrivial domain, we compa…
▽ More
For a hyponormal operator, C. R. Putnam's inequality gives an upper bound on the norm of its self-commutator. In the special case of a Toeplitz operator with analytic symbol in the Smirnov space of a domain, there is also a geometric lower bound shown by D. Khavinson (1985) that when combined with Putnam's inequality implies the classical isoperimetric inequality. For a nontrivial domain, we compare these estimates to exact results. Then we consider such operators acting on the Bergman space of a domain, and we obtain lower bounds that also reflect the geometry of the domain. When combined with Putnam's inequality they give rise to the Faber-Krahn inequality for the fundamental frequency of a domain and the Saint-Venant inequality for the torsional rigidity (but with non-sharp constants). We conjecture an improved version of Putnam's inequality within this restricted setting.
△ Less
Submitted 11 November, 2014; v1 submitted 29 November, 2012;
originally announced November 2012.
-
An Overdetermined Problem in Potential Theory
Authors:
Dmitry Khavinson,
Erik Lundberg,
Razvan Teodorescu
Abstract:
We investigate a problem posed by L. Hauswirth, F. Hélein, and F. Pacard, namely, to characterize all the domains in the plane that admit a "roof function", i.e., a positive harmonic function which solves simultaneously a Dirichlet problem with null boundary data, and a Neumann problem with constant boundary data. Under some a priori assumptions, we show that the only three examples are the exteri…
▽ More
We investigate a problem posed by L. Hauswirth, F. Hélein, and F. Pacard, namely, to characterize all the domains in the plane that admit a "roof function", i.e., a positive harmonic function which solves simultaneously a Dirichlet problem with null boundary data, and a Neumann problem with constant boundary data. Under some a priori assumptions, we show that the only three examples are the exterior of a disk, a halfplane, and a nontrivial example. We show that in four dimensions the nontrivial simply connected example does not have any axially symmetric analog containing its own axis of symmetry.
△ Less
Submitted 27 March, 2013; v1 submitted 23 May, 2012;
originally announced May 2012.
-
Non-algebraic quadrature domains
Authors:
Alexandre Eremenko,
Erik Lundberg
Abstract:
It is well known that, in the plane, the boundary of any quadrature domain (in the classical sense) coincides with the zero set of a polynomial. We show, by explicitly constructing some four-dimensional examples, that this is not always the case. This confirms, in dimension 4, a conjecture of the second author. Our method is based on the Schwarz potential and involves elliptic integrals of the thi…
▽ More
It is well known that, in the plane, the boundary of any quadrature domain (in the classical sense) coincides with the zero set of a polynomial. We show, by explicitly constructing some four-dimensional examples, that this is not always the case. This confirms, in dimension 4, a conjecture of the second author. Our method is based on the Schwarz potential and involves elliptic integrals of the third kind.
△ Less
Submitted 22 February, 2012;
originally announced February 2012.
-
Laplacian Growth, Elliptic Growth, and Singularities of the Schwarz Potential
Authors:
Erik Lundberg
Abstract:
The Schwarz function has played an elegant role in understanding and in generating new examples of exact solutions to the Laplacian growth (or "Hele- Shaw") problem in the plane. The guiding principle in this connection is the fact that "non-physical" singularities in the "oil domain" of the Schwarz function are stationary, and the "physical" singularities obey simple dynamics. We give an elementa…
▽ More
The Schwarz function has played an elegant role in understanding and in generating new examples of exact solutions to the Laplacian growth (or "Hele- Shaw") problem in the plane. The guiding principle in this connection is the fact that "non-physical" singularities in the "oil domain" of the Schwarz function are stationary, and the "physical" singularities obey simple dynamics. We give an elementary proof that the same holds in any number of dimensions for the Schwarz potential, introduced by D. Khavinson and H. S. Shapiro [17] (1989). A generalization is also given for the so-called "elliptic growth" problem by defining a generalized Schwarz potential. New exact solutions are constructed, and we solve inverse problems of describing the driving singularities of a given flow. We demonstrate, by example, how \mathbb{C}^n - techniques can be used to locate the singularity set of the Schwarz potential. One of our methods is to prolong available local extension theorems by constructing "globalizing families". We make three conjectures in potential theory relating to our investigation.
△ Less
Submitted 27 September, 2010;
originally announced September 2010.
-
Transcendental Harmonic Mappings and Gravitational Lensing by Isothermal Galaxies
Authors:
Dmitry Khavinson,
Erik Lundberg
Abstract:
Using the Schwarz function of an ellipse, it was recently shown that galaxies with density constant on confocal ellipses can produce at most four ``bright'' images of a single source. The more physically interesting example of an isothermal galaxy has density that is constant on \emph{homothetic} ellipses. In that case bright images can be seen to correspond to zeros of a certain transcendental…
▽ More
Using the Schwarz function of an ellipse, it was recently shown that galaxies with density constant on confocal ellipses can produce at most four ``bright'' images of a single source. The more physically interesting example of an isothermal galaxy has density that is constant on \emph{homothetic} ellipses. In that case bright images can be seen to correspond to zeros of a certain transcendental harmonic mapping. We use complex dynamics to give an upper bound on the total number of such zeros.
△ Less
Submitted 23 August, 2009;
originally announced August 2009.