Source author record

Robert J. McCann

Robert J. McCann appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

18works
11topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

18 published item(s)

preprint2022arXiv

Classifying minimum energy states for interacting particles: Spherical Shells

Particles interacting through long-range attraction and short-range repulsion given by power-laws have been widely used to model physical and biological systems, and to predict or explain many of the patterns they display. Apart from rare values of the attractive and repulsive exponents $(α,β)$, the energy minimizing configurations of particles are not explicitly known, although simulations and local stability considerations have led to conjectures with strong evidence over a much wider region of parameters. For a segment $β=2<α<4$ on the mildly repulsive frontier we employ strict convexity to conclude that the energy is uniquely minimized (up to translation) by a spherical shell. In a companion work, we show that in the mildly repulsive range $α>β\ge2$, a unimodal threshold $2<α_{Δ^n}(β) \le \max\{β,4\}$ exists such that equidistribution of particles over a unit diameter regular $n$-simplex minimizes the energy if and only if $α\ge α_{Δ^n}(β)$ (and minimizes uniquely up to rigid motions if strict inequality holds). At the point $(α,β)=(2,4)$ separating these regimes, we show the minimizers all lie on a sphere and are precisely characterized by sharing all first and second moments with the spherical shell. Although the minimizers need not be asymptotically stable, our approach establishes $d_α$-Lyapunov nonlinear stability of the associated ($d_2$-gradient) aggregation dynamics near the minimizer in both of these adjacent regimes -- without reference to linearization. The $L^α$-Kantorovich-Rubinstein distance $d_α$ which quantifies stability is chosen to match the attraction exponent.

preprint2022arXiv

On the cardinality of sets in ${\bf R}^d$ obeying a slightly obtuse angle bound

In this paper we explicitly estimate the number of points in a subset $A \subset \R^{d}$ as a function of the maximum angle $\angle A$ that any three of these points form, provided $\angle A < θ_d := \arccos(-\frac 1 {d}) \in (π/2,π)$. We also show $\angle A < θ_d$ ensures that $A$ coincides with the vertex set of a convex polytope. This study is motivated by a question of Paul Erdős and indirectly by a conjecture of László Fejes Tóth.

preprint2021arXiv

Maximizing expected powers of the angle between pairs of points in projective space

Among probability measures on $d$-dimensional real projective space, one which maximizes the expected angle $\arccos(\frac{x}{|x|}\cdot \frac{y}{|y|})$ between independently drawn projective points $x$ and $y$ was conjectured to equidistribute its mass over the standard Euclidean basis $\{e_0,e_1,\ldots, e_d\}$ by Fejes Tóth \cite{FT59}. If true, this conjecture evidently implies the same measure maximizes the expectation of $\arccos^α(\frac{x}{|x|}\cdot \frac{y}{|y|})$ for any exponent $α> 1$. The kernel $\arccos^α(\frac{x}{|x|}\cdot \frac{y}{|y|})$ represents the objective of an infinite-dimensional quadratic program. We verify discrete and continuous versions of this {milder} conjecture in a non-empty range $α> α_{Δ^d} \ge 1$, and establish uniqueness of the resulting maximizer $\hat μ$ up to rotation. We show $\hat μ$ no longer maximizes when $α<α_{Δ^d}$. At the endpoint $α=α_{Δ^d}$ of this range, we show another maximizer $μ$ must also exist which is not a rotation of $\hat μ$. For the continuous version of the conjecture, an appendix provided by Bilyk et al in response to an earlier draft of this work combines with the present improvements to yield $α_{Δ^d}<2$. The original conjecture $\ald=1$ remains open (unless $d=1$). However, in the maximum possible range $α>1$, we show $\hat μ$ and its rotations maximize the aforementioned expectation uniquely on a sufficiently small ball in the $L^\infty$-Kantorovich-Rubinstein-Wasserstein metric $d_\infty$ from optimal transportation; the same is true for any measure $μ$ which is mutually absolutely continuous with respect to $\hat μ$, but the size of the ball depends on {$α,d$, and} $\|\frac{d \hat μ}{dμ}\|_{\infty}$.

preprint2020arXiv

Geometrical bounds for the variance and recentered moments

We bound the variance and other moments of a random vector based on the range of its realizations, thus generalizing inequalities of Popoviciu (1935) and Bhatia and Davis (2000) concerning measures on the line to several dimensions. This is done using convex duality and (infinite-dimensional) linear programming. The following consequence of our bounds exhibits symmetry breaking, provides a new proof of Jung's theorem (1901), and turns out to have applications to the aggregation dynamics modelling attractive-repulsive interactions: among probability measures on ${\mathbf R}^n$ whose support has diameter at most $\sqrt{2}$, we show that the variance around the mean is maximized precisely by those measures which assign mass $1/(n+1)$ to each vertex of a standard simplex. For $1 \le p <\infty$, the $p$-th moment --- optimally centered --- is maximized by the same measures among those satisfying the diameter constraint.

preprint2017arXiv

On concavity of the monopolist's problem facing consumers with nonlinear price preferences

A monopolist wishes to maximize her profits by finding an optimal price policy. After she announces a menu of products and prices, each agent $x$ will choose to buy that product $y(x)$ which maximizes his own utility, if positive. The principal's profits are the sum of the net earnings produced by each product sold. These are determined by the costs of production and the distribution of products sold, which in turn are based on the distribution of anonymous agents and the choices they make in response to the principal's price menu. In this paper, we provide a necessary and sufficient condition for the convexity or concavity of the principal's (bilevel) optimization problem, assuming each agent's disutility is a strictly increasing but not necessarily affine (i.e. quasilinear) function of the price paid. Concavity when present, makes the problem more amenable to computational and theoretical analysis; it is key to obtaining uniqueness and stability results for the principal's strategy in particular. Even in the quasilinear case, our analysis goes beyond previous work by addressing convexity as well as concavity, by establishing conditions which are not only sufficient but necessary, and by requiring fewer hypotheses on the agents' preferences.

preprint2015arXiv

Academic wages and pyramid schemes: a mathematical model

This paper analyzes a steady state matching model interrelating the education and labor sectors. In this model, a heterogeneous population of students match with teachers to enhance their cognitive skills. As adults, they then choose to become workers, managers, or teachers, who match in the labor or educational market to earn wages by producing output. We study the competitive equilibrium which results from the steady state requirement that the educational process replicate the same endogenous distribution of cognitive skills among adults in each generation (assuming the same distribution of student skills). We show such an equilibrium can be found by solving an infinite-dimensional linear program and its dual. We analyze the structure of our solutions, and give sufficient conditions for them to be unique. Whether or not the educational matching is positive assortative turns out to depend on convexity of the equilibrium wages as a function of ability, suitably parameterized; we identity conditions which imply this convexity. Moreover, due to the recursive nature of the education market, it is a priori conceivable that a pyramid scheme leads to greater and greater discrepancies in the wages of the most talented teachers at the top of the market. Assuming each teacher teaches $N$ students, and contributes a fraction $θ\in]0,1[$ to their cognitive skill, we show a phase transition occurs at $Nθ=1$, which determines whether or not the wage gradients of these teachers remain bounded as market size grows, and make a quantitative prediction for their asymptotic behaviour in both regimes: $Nθ\ge 1$ and $Nθ<1$.

preprint2014arXiv

Dual potentials for capacity constrained optimal transport

Optimal transportation with capacity constraints, a variant of the well-known optimal transportation problem, is concerned with transporting one probability density $f \in L^1(\mathbb{R}^m)$ onto another one $g \in L^1(\mathbb{R}^n)$ so as to optimize a cost function $c \in L^1(\mathbb{R}^{m+n})$ while respecting the capacity constraints $0\le h \le \bar h\in L^\infty(\mathbb{R}^{m+n})$. A linear programming duality theorem for this problem was first established by Levin. In this note, we prove under mild assumptions on the given data, the existence of a pair of $L^1$-functions optimizing the dual problem. Using these functions, which can be viewed as Lagrange multipliers to the marginal constraints $f$ and $g$, we characterize the solution $h$ of the primal problem. We expect these potentials to play a key role in any further analysis of $h$. Moreover, starting from Levin's duality, we derive the classical Kantorovich duality for unconstrained optimal transport. In tandem with results obtained in our companion paper (arXiv:1309.3022), this amounts to a new and elementary proof of Kantorovich's duality.

preprint2013arXiv

The spectrum of a family of fourth-order equations near the global attractor

The thin-film and quantum drift diffusion equations belong to a fourth-order family of evolution equations proposed in ref. 16 to be analogous to the (second-order) porous medium family. They are $2$-Wasserstein gradient ($W_2$) flows of the generalized Fisher information $I(u)$ just as the porous medium family was shown to be the $W_2$ gradient flow of the generalized entropy $E(u)$ by Otto. The identity $I(u) = |\nabla_{W_2} E(u)|^2/2$ (formally) becomes $\Hess_{W_2} I(u_*) =\Hess_{W_2}^2 E(u_*)$ when linearizing the equation around its self-similar solution $u_*$. We couple this relation with the diagonalization of $\Hess_{W_2} E(u_*)$ for the porous medium flow computed in ref. 38. This yields information about the leading- and higher-order asymptotics of the equation on $\R^N$ which --- outside of special cases --- was inaccessible previously.

preprint2012arXiv

Optimal Transportation with Capacity Constraints

The classical problem of optimal transportation can be formulated as a linear optimization problem on a convex domain: among all joint measures with fixed marginals find the optimal one, where optimality is measured against a cost function. Here we consider a natural but largely unexplored variant of this problem by imposing a pointwise constraint on the joint (absolutely continuous) measures: among all joint densities with fixed marginals and which are dominated by a given density, find the optimal one. For this variant, we show local non-degeneracy of the cost function implies every minimizer is extremal in the convex set of competitors, hence unique. An appendix develops rudiments of a duality theory for this problem, which allows us to compute several suggestive examples.

preprint2011arXiv

Hölder continuity and injectivity of optimal maps

Consider transportation of one distribution of mass onto another, chosen to optimize the total expected cost, where cost per unit mass transported from x to y is given by a smooth function c(x,y). If the source density f^+(x) is bounded away from zero and infinity in an open region U' \subset R^n, and the target density f^-(y) is bounded away from zero and infinity on its support \cl V \subset R^n, which is strongly c-convex with respect to U', and the transportation cost c satisfies the A3w condition of Trudinger and Wang \cite{TrudingerWang07p}, we deduce local Hölder continuity and injectivity of the optimal map inside U' (so that the associated potential u belongs to C^{1,α}_{loc}(\U')). Here the exponent α>0 depends only on the dimension and the bounds on the densities, but not on c. Our result provides a crucial step in the low/interior regularity setting: in a sequel \cite{FigalliKimMcCann09p}, we use it to establish regularity of optimal maps with respect to the Riemannian distance squared on arbitrary products of spheres. Three key tools are introduced in the present paper. Namely, we first find a transformation that under A3w makes c-convex functions level-set convex (as was also obtained independently from us by Liu \cite{Liu09}). We then derive new Alexandrov type estimates for the level-set convex c-convex functions, and a topological lemma showing optimal maps do not mix interior with boundary. This topological lemma, which does not require A3w, is needed by Figalli and Loeper \cite{FigalliLoeper08p} to conclude continuity of optimal maps in two dimensions. In higher dimensions, if the densities f^\pm are Hölder continuous, our result permits continuous differentiability of the map inside U' (in fact, C^{2,α}_{loc} regularity of the associated potential) to be deduced from the work of Liu, Trudinger and Wang \cite{LiuTrudingerWang09p}.

preprint2011arXiv

Unique equilibria and substitution effects in a stochastic model of the marriage market

Choo-Siow (2006) proposed a model for the marriage market which allows for random identically distributed noise in the preferences of each of the participants. The randomness is McFadden-type, which permits an explicit resolution of the equilibrium preference probabilities. The purpose of this note is to prove uniqueness of the resulting equilibrium marriage distribution, and find a representation of it in closed form. This allows us to derive smooth dependence of this distribution on exogenous preference and population parameters, and establish sign, symmetry, and size of the various substitution effects, facilitating comparative statics. For example, we show that an increase in the population of men of any given type in this model leads to an increase in single men of each type, and a decrease in single women of each type. We show that an increase in the number of men of a given type increases the equilibrium transfer paid by such men to their spouses, and also increases the percentage of men of that type who choose to remain unmarried. The verification of such properties helps to substantiate the validity of the model. Moreover, we make unexpected predictions which could be tested: the percentage change of type $i$ unmarrieds with respect to fluctuations in the total number of type $j$ men or women turns out to form a symmetric positive-definite matrix $r_{ij}=r_{ji}$ in this model, and thus to satisfy bounds like $|r_{ij}| \le (r_{ii}r_{jj})^{1/2}$. Along the way, we give a new proof for the existence of an equilibrium, based on a strictly convex variational principle and a simple estimate, rather than a fixed point theorem. Fixed point approaches to the existence part of our result have been explored by others \cite{CSS} \cite{Dag} \cite{Fox}, but are much more complicated and yield neither uniqueness, nor comparative statics, nor an explicit representation of the solution.

preprint2010arXiv

Extremal doubly stochastic measures and optimal transportation

This article connects the theory of extremal doubly stochastic measures to the geometry and topology of optimal transportation. We begin by reviewing an old question (# 111) of Birkhoff in probability and statistics [4], which is to give a necessary and sufficient condition on the support of a joint probability to guarantee extremality among all measures which share its marginals. Following work of Douglas, Lindenstrauss, and Beneš and Štěpán, Hestir and Williams [15] found a necessary condition which is nearly sufficient; we relax their subtle measurability hypotheses separating necessity from sufficiency slightly, yet demonstrate by example that to be sufficient certainly requires some measurability. Their condition amounts to the vanishing of γ outside a countable alternating sequence of graphs and antigraphs in which no two graphs (or two antigraphs) have domains that overlap, and where the domain of each graph / antigraph in the sequence contains the range of the succeeding antigraph (respectively, graph). Such sequences are called numbered limb systems. Surprisingly, this characterization can be used to resolve the uniqueness question for optimal transportation on manifolds with the topology of the sphere.

preprint2010arXiv

Optimal transportation, topology and uniqueness

The Monge-Kantorovich transportation problem involves optimizing with respect to a given a cost function. Uniqueness is a fundamental open question about which little is known when the cost function is smooth and the landscapes containing the goods to be transported possess (non-trivial) topology. This question turns out to be closely linked to a delicate problem (# 111) of Birkhoff [14]: give a necessary and sufficient condition on the support of a joint probability to guarantee extremality among all measures which share its marginals. Fifty years of progress on Birkhoff's question culminate in Hestir and Williams' necessary condition which is nearly sufficient for extremality; we relax their subtle measurability hypotheses separating necessity from sufficiency slightly, yet demonstrate by example that to be sufficient certainly requires some measurability. Their condition amounts to the vanishing of the measure γoutside a countable alternating sequence of graphs and antigraphs in which no two graphs (or two antigraphs) have domains that overlap, and where the domain of each graph / antigraph in the sequence contains the range of the succeeding antigraph (respectively, graph). Such sequences are called numbered limb systems. We then explain how this characterization can be used to resolve the uniqueness of Kantorovich solutions for optimal transportation on a manifold with the topology of the sphere.

preprint2010arXiv

Regularity of optimal transport maps on multiple products of spheres

This article addresses regularity of optimal transport maps for cost="squared distance" on Riemannian manifolds that are products of arbitrarily many round spheres with arbitrary sizes and dimensions. Such manifolds are known to be non-negatively cross-curved [KM2]. Under boundedness and non-vanishing assumptions on the transfered source and target densities we show that optimal maps stay away from the cut-locus (where the cost exhibits singularity), and obtain injectivity and continuity of optimal maps. This together with the result of Liu, Trudinger and Wang [LTW] also implies higher regularity (C^{1,α}/C^\infty) of optimal maps for more smooth (C^α/C^\infty)) densities. These are the first global regularity results which we are aware of concerning optimal maps on non-flat Riemannian manifolds which possess some vanishing sectional curvatures. Moreover, such product manifolds have potential relevance in statistics (see [S]) and in statistical mechanics (where the state of a system consisting of many spins is classically modeled by a point in the phase space obtained by taking many products of spheres). For the proof we apply and extend the method developed in [FKM1], where we showed injectivity and continuity of optimal maps on domains in R^n for smooth non-negatively cross-curved cost. The major obstacle in the present paper is to deal with the non-trivial cut-locus and the presence of flat directions.