-
A charged liquid drop model with Willmore energy
Authors:
Michael Goldman,
Matteo Novaga,
Berardo Ruffini
Abstract:
We consider a variational model of electrified liquid drops, involving competition between surface tension and charge repulsion. Since the natural model happens to be ill-posed, we show that by adding to the perimeter a Willmore-type energy, the problem turns back to be well-posed. We also prove that for small charge the droplets is spherical.
We consider a variational model of electrified liquid drops, involving competition between surface tension and charge repulsion. Since the natural model happens to be ill-posed, we show that by adding to the perimeter a Willmore-type energy, the problem turns back to be well-posed. We also prove that for small charge the droplets is spherical.
△ Less
Submitted 2 September, 2024;
originally announced September 2024.
-
Partial regularity for optimal transport with p-cost away from fixed points
Authors:
Michael Goldman,
Lukas Koch
Abstract:
We consider maps $T$ solving the optimal transport problem with a cost $c(x-y)$ modeled on the $p$-cost. For Hölder continuous marginals, we prove a $C^{1,α}$-partial regularity result for $T $in the set $\{|T(x)-x|>0\}$.
We consider maps $T$ solving the optimal transport problem with a cost $c(x-y)$ modeled on the $p$-cost. For Hölder continuous marginals, we prove a $C^{1,α}$-partial regularity result for $T $in the set $\{|T(x)-x|>0\}$.
△ Less
Submitted 11 July, 2024;
originally announced July 2024.
-
Subadditivity and optimal matching of unbounded samples
Authors:
Emanuele Caglioti,
Michael Goldman,
Francesca Pieroni,
Dario Trevisan
Abstract:
We obtain new bounds for the optimal matching cost for empirical measures with unbounded support. For a large class of radially symmetric and rapidly decaying probability laws, we prove for the first time the asymptotic rate of convergence for the whole range of power exponents $p$ and dimensions $d$. Moreover we identify the exact prefactor when $p\le d$. We cover in particular the Gaussian case,…
▽ More
We obtain new bounds for the optimal matching cost for empirical measures with unbounded support. For a large class of radially symmetric and rapidly decaying probability laws, we prove for the first time the asymptotic rate of convergence for the whole range of power exponents $p$ and dimensions $d$. Moreover we identify the exact prefactor when $p\le d$. We cover in particular the Gaussian case, going far beyond the currently known bounds. Our proof technique is based on approximate sub- and super-additivity bounds along a geometric decomposition adapted to some features the density, such as its radial symmetry and its decay at infinity.
△ Less
Submitted 8 July, 2024;
originally announced July 2024.
-
A Lyapunov theory demonstrating a fundamental limit on the speed of systems consolidation
Authors:
Alireza Alemi,
Emre R. F. Aksay,
Mark S. Goldman
Abstract:
The nervous system reorganizes memories from an early site to a late site, a commonly observed feature of learning and memory systems known as systems consolidation. Previous work has suggested learning rules by which consolidation may occur. Here, we provide conditions under which such rules are guaranteed to lead to stable convergence of learning and consolidation. We use the theory of Lyapunov…
▽ More
The nervous system reorganizes memories from an early site to a late site, a commonly observed feature of learning and memory systems known as systems consolidation. Previous work has suggested learning rules by which consolidation may occur. Here, we provide conditions under which such rules are guaranteed to lead to stable convergence of learning and consolidation. We use the theory of Lyapunov functions, which enforces stability by requiring learning rules to decrease an energy-like (Lyapunov) function. We present the theory in the context of a simple circuit architecture motivated by classic models of learning in systems consolidation mediated by the cerebellum. Stability is only guaranteed if the learning rate in the late stage is not faster than the learning rate in the early stage. Further, the slower the learning rate at the late stage, the larger the perturbation the system can tolerate with a guarantee of stability. We provide intuition for this result by mapping the consolidation model to a damped driven oscillator system, and showing that the ratio of early- to late-stage learning rates in the consolidation model can be directly identified with the (square of the) oscillator's damping ratio. This work suggests the power of the Lyapunov approach to provide constraints on nervous system function.
△ Less
Submitted 2 February, 2024;
originally announced February 2024.
-
Almost sharp rates of convergence for the average cost and displacement in the optimal matching problem
Authors:
Michael Goldman,
Martin Huesmann,
Felix Otto
Abstract:
In this note we prove estimates for the average cost in the quadratic optimal transport problem on the two-dimensional flat torus which are optimal up to a double logarithm. We also prove sharp estimates on the displacement. This is based on the combination of a post-processing of our quantitative linearization result together with a quasi-orthogonality property.
In this note we prove estimates for the average cost in the quadratic optimal transport problem on the two-dimensional flat torus which are optimal up to a double logarithm. We also prove sharp estimates on the displacement. This is based on the combination of a post-processing of our quantitative linearization result together with a quasi-orthogonality property.
△ Less
Submitted 13 December, 2023;
originally announced December 2023.
-
Non-convex functionals penalizing simultaneous oscillations along two independent directions: structure of the defect measure
Authors:
Michael Goldman,
Benoît Merlet
Abstract:
We continue the analysis of a family of energies penalizing oscillations in oblique directions: they apply to functions $u(x_1,x_2)$ with $x_l\in\mathbb{R}^{n_l}$ and vanish when $u(x)$ is of the form $u_1(x_1)$ or $u_2(x_2)$. We mainly study the rectifiability properties of the defect measure $\nabla_1\nabla_2u$ of functions with finite energy. The energies depend on a parameter $θ\in(0,1]$ and t…
▽ More
We continue the analysis of a family of energies penalizing oscillations in oblique directions: they apply to functions $u(x_1,x_2)$ with $x_l\in\mathbb{R}^{n_l}$ and vanish when $u(x)$ is of the form $u_1(x_1)$ or $u_2(x_2)$. We mainly study the rectifiability properties of the defect measure $\nabla_1\nabla_2u$ of functions with finite energy. The energies depend on a parameter $θ\in(0,1]$ and the set of functions with finite energy grows with $θ$. For $θ<1$ we prove that the defect measure is $(n_1-1,n_2-1)$-tensor rectifiable in $Ω_1\timesΩ_2$. We first get the result for $n_1=n_2=1$ and deduce the general case through slicing using White's rectifiability criterion. When $θ=1$ the situation is less clear as measures of arbitrary dimensions from zero to $n_1+n_2-1$ are possible. We show however, in the case $n_1=n_2=1$ and for Lipschitz continuous functions, that the defect measures are $1\,$-rectifiable. This case bears strong analogies with the study of entropic solutions of the eikonal equation.
△ Less
Submitted 29 September, 2023;
originally announced September 2023.
-
An exterior optimal transport problem
Authors:
Jules Candau-Tilh,
Michael Goldman,
Benoit Merlet
Abstract:
This paper deals with a variant of the optimal transportation problem. Given f $\in$ L 1 (R d , [0, 1]) and a cost function c $\in$ C(R d x R d) of the form c(x, y) = k(y -- x), we minimise $\int$ c d$γ$ among transport plans $γ$ whose first marginal is f and whose second marginal is not prescribed but constrained to be smaller than 1 -- f. Denoting by $Υ$(f) the infimum of this problem, we then c…
▽ More
This paper deals with a variant of the optimal transportation problem. Given f $\in$ L 1 (R d , [0, 1]) and a cost function c $\in$ C(R d x R d) of the form c(x, y) = k(y -- x), we minimise $\int$ c d$γ$ among transport plans $γ$ whose first marginal is f and whose second marginal is not prescribed but constrained to be smaller than 1 -- f. Denoting by $Υ$(f) the infimum of this problem, we then consider the maximisation problem sup{$Υ$(f) : $\int$ f = m} where m > 0 is given. We prove that maximisers exist under general assumptions on k, and that for k radial, increasing and coercive these maximisers are the characteristic functions of the balls of volume m.
△ Less
Submitted 6 September, 2023;
originally announced September 2023.
-
Evaluation of software impact designed for biomedical research: Are we measuring what's meaningful?
Authors:
Awan Afiaz,
Andrey Ivanov,
John Chamberlin,
David Hanauer,
Candace Savonen,
Mary J Goldman,
Martin Morgan,
Michael Reich,
Alexander Getka,
Aaron Holmes,
Sarthak Pati,
Dan Knight,
Paul C. Boutros,
Spyridon Bakas,
J. Gregory Caporaso,
Guilherme Del Fiol,
Harry Hochheiser,
Brian Haas,
Patrick D. Schloss,
James A. Eddy,
Jake Albrecht,
Andrey Fedorov,
Levi Waldron,
Ava M. Hoffman,
Richard L. Bradshaw
, et al. (2 additional authors not shown)
Abstract:
Software is vital for the advancement of biology and medicine. Analysis of usage and impact metrics can help developers determine user and community engagement, justify additional funding, encourage additional use, identify unanticipated use cases, and help define improvement areas. However, there are challenges associated with these analyses including distorted or misleading metrics, as well as e…
▽ More
Software is vital for the advancement of biology and medicine. Analysis of usage and impact metrics can help developers determine user and community engagement, justify additional funding, encourage additional use, identify unanticipated use cases, and help define improvement areas. However, there are challenges associated with these analyses including distorted or misleading metrics, as well as ethical and security concerns. More attention to the nuances involved in capturing impact across the spectrum of biological software is needed. Furthermore, some tools may be especially beneficial to a small audience, yet may not have compelling typical usage metrics. We propose more general guidelines, as well as strategies for more specific types of software. We highlight outstanding issues regarding how communities measure or evaluate software impact. To get a deeper understanding of current practices for software evaluations, we performed a survey of participants in the Informatics Technology for Cancer Research (ITCR) program funded by the National Cancer Institute (NCI). We also investigated software among this community and others to assess how often infrastructure that supports such evaluations is implemented and how this impacts rates of papers describing usage of the software. We find that developers recognize the utility of analyzing software usage, but struggle to find the time or funding for such analyses. We also find that infrastructure such as social media presence, more in-depth documentation, the presence of software health metrics, and clear information on how to contact developers seem to be associated with increased usage rates. Our findings can help scientific software developers make the most out of evaluations of their software.
△ Less
Submitted 5 June, 2023;
originally announced June 2023.
-
On the concave one-dimensional random assignment problem and Young integration theory
Authors:
Michael Goldman,
Dario Trevisan
Abstract:
We investigate the one-dimensional random assignment problem in the concave case, i.e., the assignment cost is a concave power function, with exponent $0<p<1$, of the distance between $n$ source and $n$ target points, that are i.i.d. random variables with a common law on an interval. We prove that the limit of a suitable renormalization of the costs exists if the exponent $p$ is different than…
▽ More
We investigate the one-dimensional random assignment problem in the concave case, i.e., the assignment cost is a concave power function, with exponent $0<p<1$, of the distance between $n$ source and $n$ target points, that are i.i.d. random variables with a common law on an interval. We prove that the limit of a suitable renormalization of the costs exists if the exponent $p$ is different than $1/2$. Our proof in the case $1/2<p<1$ makes use of a novel version of the Kantorovich optimal transport problem based on Young integration theory, where the difference between two measures is replaced by the weak derivative of a function with finite $q$-variation, which may be of independent interest. We also prove a similar result for the random bipartite Traveling Salesperson Problem.
△ Less
Submitted 16 May, 2023;
originally announced May 2023.
-
From energy bounds to dimensional estimates in a branched transport model for type-I superconductors
Authors:
Guido De Philippis,
Michael Goldman,
Berardo Ruffini
Abstract:
We consider a branched transport type problem which describes the magnetic flux through type-I superconductors in a regime of very weak applied fields. At the boundary of the sample, deviation of the magnetization from being uniform is penalized through a negative Sobolev norm. It was conjectured by S. Conti, F. Otto and S. Serfaty that as a result, the trace of the magnetization on the boundary s…
▽ More
We consider a branched transport type problem which describes the magnetic flux through type-I superconductors in a regime of very weak applied fields. At the boundary of the sample, deviation of the magnetization from being uniform is penalized through a negative Sobolev norm. It was conjectured by S. Conti, F. Otto and S. Serfaty that as a result, the trace of the magnetization on the boundary should be a measure of Hausdorff dimension $8/5$. We prove that this conjecture is equivalent to the proof of local energy bounds with an optimal exponent. We then obtain local bounds which are however not optimal. These yield improved lower bounds on the dimension of the irrigated measure but unfortunately does not improve on the trivial upper bound. In order to illustrate the dependence of this dimension on the choice of penalization, we consider in the last part of the paper a toy model where the boundary energy is given by a Wasserstein distance to Lebesgue. In this case minimizers are finite graphs and thus the trace is atomic.
△ Less
Submitted 25 April, 2023;
originally announced April 2023.
-
A deformation theorem for tensor flat chains and applications (complement to ''Tensor rectifiable G-flat chains'')
Authors:
Michael Goldman,
Benoît Merlet
Abstract:
In this note we extend White's deformation theorem for G-flat chains to the setting of G-flat tensor chains. As a corollary we obtain that the groups of normal tensor chains identify with some subgroups of normal chains. Moreover the corresponding natural group isomorphisms are isometric with respect to norms based on the coordinate slicing mass. The coordinate slicing mass of a k-chain is the int…
▽ More
In this note we extend White's deformation theorem for G-flat chains to the setting of G-flat tensor chains. As a corollary we obtain that the groups of normal tensor chains identify with some subgroups of normal chains. Moreover the corresponding natural group isomorphisms are isometric with respect to norms based on the coordinate slicing mass. The coordinate slicing mass of a k-chain is the integral of the mass of its 0-slices along all coordinate-planes of codimension k. The fact that this quantity is equivalent to the usual mass is not straightforward. To prove it, we use the deformation theorem and a partial extension of the restriction operator defined for all chains (not only of finite mass). On the contrary, except in some limit or degenerate cases, the whole groups of tensor chains and of finite mass tensor chains do not identify naturally with subgroups of chains.
△ Less
Submitted 13 December, 2022;
originally announced December 2022.
-
Tensor rectifiable G-flat chains
Authors:
Michael Goldman,
Benoît Merlet
Abstract:
A rigidity result for normal rectifiable $k$-chains in $\mathbb{R}^n$ with coefficients in an Abelian normed group is established. Given some decompositions $k=k_1+k_2$, $n=n_1+n_2$ and some rectifiable $k$-chain $A$ in $\mathbb{R}^n$, we consider the properties:(1) The tangent planes to $μ_A$ split as $T_xμ_A=L^1(x)\times L^2(x)$ for some $k_1$-plane $L^1(x)\subset\mathbb{R}^{n_1}$ and some…
▽ More
A rigidity result for normal rectifiable $k$-chains in $\mathbb{R}^n$ with coefficients in an Abelian normed group is established. Given some decompositions $k=k_1+k_2$, $n=n_1+n_2$ and some rectifiable $k$-chain $A$ in $\mathbb{R}^n$, we consider the properties:(1) The tangent planes to $μ_A$ split as $T_xμ_A=L^1(x)\times L^2(x)$ for some $k_1$-plane $L^1(x)\subset\mathbb{R}^{n_1}$ and some $k_2$-plane $L^2(x)\subset\mathbb{R}^{n_2}$.(2) $A=A_{\vertΣ^1\timesΣ^2}$ for some sets $Σ^1\subset\mathbb{R}^{n_1}$, $Σ^2\subset\mathbb{R}^{n_2}$ such that $Σ^1$ is $k_1$-rectifiable and $Σ^2$ is $k_2$-rectifiable (we say that $A$ is $(k_1,k_2)$-rectifiable).The main result is that for normal chains, (1) implies (2), the converse is immediate. In the proof we introduce the new groups of tensor flat chains (or $(k_1,k_2)$-chains) in $\mathbb{R}^{n_1}\times\mathbb{R}^{n_2}$ which generalize Fleming's $G$-flat chains. The other main tool is White's rectifiable slices theorem. We show that on the one hand any normal rectifiable chain satisfying~(1) identifies with a normal rectifiable $(k_1,k_2)$-chain and that on the other hand any normal rectifiable $(k_1,k_2)$-chain is $(k_1,k_2)$-rectifiable.
△ Less
Submitted 9 December, 2022;
originally announced December 2022.
-
Set-decomposition of normal rectifiable G-chains via an abstract decomposition principle
Authors:
Michael Goldman,
Benoît Merlet
Abstract:
We introduce the notion of set-decomposition of a normal G-flat chain. We show that any normal rectifiable $G$-flat chain admits a decomposition in set-indecomposable sub-chains. This generalizes the decomposition of sets of finite perimeter in their ``measure theoretic'' connected components due to Ambrosio, Caselles, Masnou and Morel. It can also be seen as a variant of the decomposition of int…
▽ More
We introduce the notion of set-decomposition of a normal G-flat chain. We show that any normal rectifiable $G$-flat chain admits a decomposition in set-indecomposable sub-chains. This generalizes the decomposition of sets of finite perimeter in their ``measure theoretic'' connected components due to Ambrosio, Caselles, Masnou and Morel. It can also be seen as a variant of the decomposition of integral currents in indecomposable components by Federer.As opposed to previous results, we do not assume that G is boundedly compact. Therefore we cannot rely on the compactness of sequences of chains with uniformly bounded N-norms. We deduce instead the result from a new abstract decomposition principle. As in earlier proofs a central ingredient is the validity of an isoperimetric inequality. We obtain it here using the finiteness of some h-mass to replace integrality.
△ Less
Submitted 19 June, 2023; v1 submitted 9 December, 2022;
originally announced December 2022.
-
Optimal transport methods for combinatorial optimization over two random point sets
Authors:
Michael Goldman,
Dario Trevisan
Abstract:
We investigate the minimum cost of a wide class of combinatorial optimization problems over random bipartite geometric graphs in $\mathbb{R}^d$ where the edge cost between two points is given by a $p$-th power of their Euclidean distance. This includes e.g.\ the travelling salesperson problem and the bounded degree minimum spanning tree. We establish in particular almost sure convergence, as $n$ g…
▽ More
We investigate the minimum cost of a wide class of combinatorial optimization problems over random bipartite geometric graphs in $\mathbb{R}^d$ where the edge cost between two points is given by a $p$-th power of their Euclidean distance. This includes e.g.\ the travelling salesperson problem and the bounded degree minimum spanning tree. We establish in particular almost sure convergence, as $n$ grows, of a suitable renormalization of the random minimum cost, if the points are uniformly distributed and $d \ge 3$, $1\le p<d$. Previous results were limited to the range $p<d/2$. Our proofs are based on subadditivity methods and build upon new bounds for random instances of the Euclidean bipartite matching problem, obtained through its optimal transport relaxation and functional analytic techniques.
△ Less
Submitted 19 July, 2023; v1 submitted 29 September, 2022;
originally announced September 2022.
-
Uniform $C^{1,α}$-regularity for almost-minimizers of some nonlocal perturbations of the perimeter
Authors:
Michael Goldman,
Benoît Merlet,
Marc Pegon
Abstract:
In this paper, we establish a $C^{1,α}$-regularity theorem for almost-minimizers of the functional $\mathcal{F}_{\varepsilon,γ}=P-γP_{\varepsilon}$, where $γ\in(0,1)$ and $P_{\varepsilon}$ is a nonlocal energy converging to the perimeter as $\varepsilon$ vanishes. Our theorem provides a criterion for $C^{1,α}$-regularity at a point of the boundary which is uniform as the parameter $\varepsilon$ go…
▽ More
In this paper, we establish a $C^{1,α}$-regularity theorem for almost-minimizers of the functional $\mathcal{F}_{\varepsilon,γ}=P-γP_{\varepsilon}$, where $γ\in(0,1)$ and $P_{\varepsilon}$ is a nonlocal energy converging to the perimeter as $\varepsilon$ vanishes. Our theorem provides a criterion for $C^{1,α}$-regularity at a point of the boundary which is uniform as the parameter $\varepsilon$ goes to $0$. As a consequence we obtain that volume-constrained minimizers of $\mathcal{F}_{\varepsilon,γ}$ are balls for any $\varepsilon$ small enough. For small $\varepsilon$, this minimization problem corresponds to the large mass regime for a Gamow-type problem where the nonlocal repulsive term is given by an integrable kernel $G$ with sufficiently fast decay at infinity.
△ Less
Submitted 22 September, 2022;
originally announced September 2022.
-
Rigidity of the ball for an isoperimetric problem with strong capacitary repulsion
Authors:
Michael Goldman,
Matteo Novaga,
Berardo Ruffini
Abstract:
We consider a variational problem involving competition between surface tension and charge repulsion. We show that, as opposed to the case of weak (short-range) interactions where we proved ill-posedness of the problem in a previous paper, when the repulsion is stronger the perimeter dominates the capacitary term at small scales. In particular we prove existence of minimizers for small charges as…
▽ More
We consider a variational problem involving competition between surface tension and charge repulsion. We show that, as opposed to the case of weak (short-range) interactions where we proved ill-posedness of the problem in a previous paper, when the repulsion is stronger the perimeter dominates the capacitary term at small scales. In particular we prove existence of minimizers for small charges as well as their regularity. Combining this with the stability of the ball under small $C^{1,γ}$ perturbations, this ultimately leads to the minimality of the ball for small charges. We cover in particular the borderline case of the $1-$capacity where both terms in the energy are of the same order.
△ Less
Submitted 12 January, 2022;
originally announced January 2022.
-
Compactness and structure of zero-states for unoriented Aviles-Giga functionals
Authors:
Michael Goldman,
Benoit Merlet,
Marc Pegon,
Sylvia Serfaty
Abstract:
Motivated by some models of pattern formation involving an unoriented director field in the plane, we study a family of unoriented counterparts to the Aviles-Giga functional. We introduce a nonlinear curl operator for such unoriented vector fields as well as a family of even entropies which we call "trigonometric entropies". Using these tools we show two main theorems which parallel some results i…
▽ More
Motivated by some models of pattern formation involving an unoriented director field in the plane, we study a family of unoriented counterparts to the Aviles-Giga functional. We introduce a nonlinear curl operator for such unoriented vector fields as well as a family of even entropies which we call "trigonometric entropies". Using these tools we show two main theorems which parallel some results in the literature on the classical Aviles-Giga energy. The first is a compactness result for sequences of configurations with uniformly bounded energies. The second is a complete characterization of zero-states, that is, the limit configurations when the energies go to 0. These are Lipschitz continuous away from a locally finite set of points, near which they form either a vortex pattern or a disclination with degree 1/2. The proof is based on a combination of regularity theory together with techniques coming from the study of the Ginzburg-Landau energy. Our methods provide alternative proofs in the classical Aviles-Giga context.
△ Less
Submitted 9 December, 2021;
originally announced December 2021.
-
On the quadratic random matching problem in two-dimensional domains
Authors:
Luigi Ambrosio,
Michael Goldman,
Dario Trevisan
Abstract:
We investigate the average minimum cost of a bipartite matching, with respect to the squared Euclidean distance, between two samples of n i.i.d. random points on a bounded Lipschitz domain in the Euclidean plane, whose common law is absolutely continuous with strictly positive H{ö}lder continuous density. We confirm in particular the validity of a conjecture by D. Benedetto and E. Caglioti stating…
▽ More
We investigate the average minimum cost of a bipartite matching, with respect to the squared Euclidean distance, between two samples of n i.i.d. random points on a bounded Lipschitz domain in the Euclidean plane, whose common law is absolutely continuous with strictly positive H{ö}lder continuous density. We confirm in particular the validity of a conjecture by D. Benedetto and E. Caglioti stating that the asymptotic cost as n grows is given by the logarithm of n multiplied by an explicit constant times the volume of the domain. Our proof relies on a reduction to the optimal transport problem between the associated empirical measures and a Whitney-type decomposition of the domain, together with suitable upper and lower bounds for local and global contributions, both ultimately based on PDE tools. We further show how to extend our results to more general settings, including Riemannian manifolds, and also give an application to the asymptotic cost of the random quadratic bipartite travelling salesperson problem.
△ Less
Submitted 27 October, 2021;
originally announced October 2021.
-
Existence and stability results for an isoperimetric problem with a non-local interaction of Wasserstein type
Authors:
Jules Candau-Tilh,
Michael Goldman
Abstract:
The aim of this paper is to prove the existence of minimizers for a variational problem involving the minimization under volume constraint of the sum of the perimeter and a non-local energy of Wasserstein type. This extends previous partial results to the full range of parameters. We also show that in the regime where the perimeter is dominant, the energy is uniquely minimized by balls.
The aim of this paper is to prove the existence of minimizers for a variational problem involving the minimization under volume constraint of the sum of the perimeter and a non-local energy of Wasserstein type. This extends previous partial results to the full range of parameters. We also show that in the regime where the perimeter is dominant, the energy is uniquely minimized by balls.
△ Less
Submitted 25 August, 2021;
originally announced August 2021.
-
Machine Learning for Variance Reduction in Online Experiments
Authors:
Yongyi Guo,
Dominic Coey,
Mikael Konutgan,
Wenting Li,
Chris Schoener,
Matt Goldman
Abstract:
We consider the problem of variance reduction in randomized controlled trials, through the use of covariates correlated with the outcome but independent of the treatment. We propose a machine learning regression-adjusted treatment effect estimator, which we call MLRATE. MLRATE uses machine learning predictors of the outcome to reduce estimator variance. It employs cross-fitting to avoid overfittin…
▽ More
We consider the problem of variance reduction in randomized controlled trials, through the use of covariates correlated with the outcome but independent of the treatment. We propose a machine learning regression-adjusted treatment effect estimator, which we call MLRATE. MLRATE uses machine learning predictors of the outcome to reduce estimator variance. It employs cross-fitting to avoid overfitting biases, and we prove consistency and asymptotic normality under general conditions. MLRATE is robust to poor predictions from the machine learning step: if the predictions are uncorrelated with the outcomes, the estimator performs asymptotically no worse than the standard difference-in-means estimator, while if predictions are highly correlated with outcomes, the efficiency gains are large. In A/A tests, for a set of 48 outcome metrics commonly monitored in Facebook experiments the estimator has over 70% lower variance than the simple difference-in-means estimator, and about 19% lower variance than the common univariate procedure which adjusts only for pre-experiment values of the outcome.
△ Less
Submitted 6 January, 2022; v1 submitted 14 June, 2021;
originally announced June 2021.
-
A fluctuation result for the displacement in the optimal matching problem
Authors:
Michael Goldman,
Martin Huesmann
Abstract:
The aim of this paper is to justify in dimensions two and three the ansatz of Caracciolo et al. stating that the displacement in the optimal matching problem is essentially given by the solution to the linearized equation i.e. the Poisson equation. Moreover, we prove that at all mesoscopic scales, this displacement is close in suitable negative Sobolev spaces to a curl-free Gaussian free field. Fo…
▽ More
The aim of this paper is to justify in dimensions two and three the ansatz of Caracciolo et al. stating that the displacement in the optimal matching problem is essentially given by the solution to the linearized equation i.e. the Poisson equation. Moreover, we prove that at all mesoscopic scales, this displacement is close in suitable negative Sobolev spaces to a curl-free Gaussian free field. For this we combine a quantitative estimate on the difference between the displacement and the linearized object, which is based on the large-scale regularity theory recently developed in collaboration with F. Otto, together with a qualitative convergence result for the linearized problem.
△ Less
Submitted 6 May, 2021;
originally announced May 2021.
-
Reifenberg flatness for almost-minimizers of the perimeter under minimal assumptions
Authors:
Michael Goldman,
Matteo Novaga,
Berardo Ruffini
Abstract:
The aim of this note is to prove that almost-minimizers of the perimeter are Reifenberg flat, for a very weak notion of minimality. The main observation is that smallness of the excess at some scale implies smallness of the excess at all smaller scales.
The aim of this note is to prove that almost-minimizers of the perimeter are Reifenberg flat, for a very weak notion of minimality. The main observation is that smallness of the excess at some scale implies smallness of the excess at all smaller scales.
△ Less
Submitted 17 June, 2021; v1 submitted 20 April, 2021;
originally announced April 2021.
-
Multi-Beam Energy Moments of Compound Measured Ion Velocity Distributions
Authors:
Martin V. Goldman,
David. L. Newman,
Jonathan P. Eastwood,
Giovanni Lapenta,
James L. Burch,
Barbara Giles
Abstract:
Compound ion distributions, fi(v), have been measured by NASA's Magnetospheric Multi-Scale Mission (MMS) and have been found in reconnection simulations. A complex distribution, fi(v), consisting, for example, of essentially disjoint pieces will be called a multi-beam distribution and modeled as a sum of "beams," fi(v) = f1(v) + ... +fN(v). Velocity moments of fi(v) are taken beam by beam and summ…
▽ More
Compound ion distributions, fi(v), have been measured by NASA's Magnetospheric Multi-Scale Mission (MMS) and have been found in reconnection simulations. A complex distribution, fi(v), consisting, for example, of essentially disjoint pieces will be called a multi-beam distribution and modeled as a sum of "beams," fi(v) = f1(v) + ... +fN(v). Velocity moments of fi(v) are taken beam by beam and summed. Such multi-beam moments of fi(v) have advantages over the customary standard velocity moments of fi(v), forwhich there is only one mean flow velocity. For example, the standard thermal energy momentof a pair of equal and opposite cold particle beams is non-zero even though each beam has zero thermal energy. We therefore call this thermal energy pseudo-thermal. By contrast, a multi-beam moment of two or more beams has no pseudo-thermal energy. We develop three different ways of decomposing into a sum and finding multi-beam moments for both a multi-beam fi(v) measured by MMS in the dayside magnetosphere during reconnection and a multi-beam fi(v) found in a PIC simulation of magnetotail reconnection. The three methods are: A visual method in which the velocity centroid of each beam is estimated and its density determined self-consistently; A k-means method in which particles in a particle-representation of fi(v) are sorted into a minimum energy configuration of N (= k) clusters; A nonlinear least squares method based on a fit to a sum of N kappa functions.
Multi-beam energy moments are calculated and compared with standard moments for the thermal energy density, pressure tensor, thermal energy flux (heat plus enthalpy fluxes), bulk kinetic energy density, RAM pressure and bulk kinetic energy flux. Applying this new formalism to real data demonstrates in detail how multi-beam techniques may provide significant insight into the properties of observed space plasmas.
△ Less
Submitted 11 July, 2021; v1 submitted 25 January, 2021;
originally announced January 2021.
-
Calculations of the norms for monotone operators on the cones of functions with monotonicity properties
Authors:
E. G. Bakhtigareeva,
M. L. Goldman
Abstract:
The paper is devoted to the problem of exact calculation of the norms in ideal spaces for monotone operators on the cones of functions with monotonicity properties. We implement a general approach to this problem that covers many concrete variants of monotone operators in ideal spaces and different monotonicity conditions for functions. As applications, we calculate the norms of some integral oper…
▽ More
The paper is devoted to the problem of exact calculation of the norms in ideal spaces for monotone operators on the cones of functions with monotonicity properties. We implement a general approach to this problem that covers many concrete variants of monotone operators in ideal spaces and different monotonicity conditions for functions. As applications, we calculate the norms of some integral operators on the cones, associate norms over some cones in Lebesgue spaces, the norms of the dilation operator and embedding operators on weighted Lorentz spaces with general weights. Under some more general conditions, we present order sharp estimates for Hardy-type operators on the cones.
△ Less
Submitted 18 January, 2021;
originally announced January 2021.
-
Network experimentation at scale
Authors:
Brian Karrer,
Liang Shi,
Monica Bhole,
Matt Goldman,
Tyrone Palmer,
Charlie Gelman,
Mikael Konutgan,
Feng Sun
Abstract:
We describe our framework, deployed at Facebook, that accounts for interference between experimental units through cluster-randomized experiments. We document this system, including the design and estimation procedures, and detail insights we have gained from the many experiments that have used this system at scale. We introduce a cluster-based regression adjustment that substantially improves pre…
▽ More
We describe our framework, deployed at Facebook, that accounts for interference between experimental units through cluster-randomized experiments. We document this system, including the design and estimation procedures, and detail insights we have gained from the many experiments that have used this system at scale. We introduce a cluster-based regression adjustment that substantially improves precision for estimating global treatment effects as well as testing for interference as part of our estimation procedure. With this regression adjustment, we find that imbalanced clusters can better account for interference than balanced clusters without sacrificing accuracy. In addition, we show how logging exposure to a treatment can be used for additional variance reduction. Interference is a widely acknowledged issue with online field experiments, yet there is less evidence from real-world experiments demonstrating interference in online settings. We fill this gap by describing two case studies that capture significant network effects and highlight the value of this experimentation framework.
△ Less
Submitted 15 December, 2020;
originally announced December 2020.
-
Beyond Weak Perspective for Monocular 3D Human Pose Estimation
Authors:
Imry Kissos,
Lior Fritz,
Matan Goldman,
Omer Meir,
Eduard Oks,
Mark Kliger
Abstract:
We consider the task of 3D joints location and orientation prediction from a monocular video with the skinned multi-person linear (SMPL) model. We first infer 2D joints locations with an off-the-shelf pose estimation algorithm. We use the SPIN algorithm and estimate initial predictions of body pose, shape and camera parameters from a deep regression neural network. We then adhere to the SMPLify al…
▽ More
We consider the task of 3D joints location and orientation prediction from a monocular video with the skinned multi-person linear (SMPL) model. We first infer 2D joints locations with an off-the-shelf pose estimation algorithm. We use the SPIN algorithm and estimate initial predictions of body pose, shape and camera parameters from a deep regression neural network. We then adhere to the SMPLify algorithm which receives those initial parameters, and optimizes them so that inferred 3D joints from the SMPL model would fit the 2D joints locations. This algorithm involves a projection step of 3D joints to the 2D image plane. The conventional approach is to follow weak perspective assumptions which use ad-hoc focal length. Through experimentation on the 3D Poses in the Wild (3DPW) dataset, we show that using full perspective projection, with the correct camera center and an approximated focal length, provides favorable results. Our algorithm has resulted in a winning entry for the 3DPW Challenge, reaching first place in joints orientation accuracy.
△ Less
Submitted 14 September, 2020;
originally announced September 2020.
-
Convergence of asymptotic costs for random Euclidean matching problems
Authors:
Michael Goldman,
Dario Trevisan
Abstract:
We investigate the average minimum cost of a bipartite matching between two samples of n independent random points uniformly distributed on a unit cube in d $\ge$ 3 dimensions, where the matching cost between two points is given by any power p $\ge$ 1 of their Euclidean distance. As n grows, we prove convergence, after a suitable renormalization, towards a finite and positive constant. We also con…
▽ More
We investigate the average minimum cost of a bipartite matching between two samples of n independent random points uniformly distributed on a unit cube in d $\ge$ 3 dimensions, where the matching cost between two points is given by any power p $\ge$ 1 of their Euclidean distance. As n grows, we prove convergence, after a suitable renormalization, towards a finite and positive constant. We also consider the analogous problem of optimal transport between n points and the uniform measure. The proofs combine sub-additivity inequalities with a PDE ansatz similar to the one proposed in the context of the matching problem in two dimensions and later extended to obtain upper bounds in higher dimensions.
△ Less
Submitted 17 December, 2020; v1 submitted 9 September, 2020;
originally announced September 2020.
-
Optimal Calderón Spaces for generalized Bessel potentials
Authors:
Elza Bakhtigareeva,
Mikhail L. Goldman,
Dorothee D. Haroske
Abstract:
In the paper we investigate the properties of spaces with generalized smoothness, such as Calderón spaces that include the classical Nikolskii-Besov spaces and many of their generalizations, and describe differential properties of generalized Bessel potentials that include classical Bessel potentials and Sobolev spaces. Kernels of potentials may have non-power singularity at the origin. With the h…
▽ More
In the paper we investigate the properties of spaces with generalized smoothness, such as Calderón spaces that include the classical Nikolskii-Besov spaces and many of their generalizations, and describe differential properties of generalized Bessel potentials that include classical Bessel potentials and Sobolev spaces. Kernels of potentials may have non-power singularity at the origin. With the help of order-sharp estimates for moduli of continuity of potentials, we establish the criteria of embeddings of potentials into Calderón spaces, and describe the optimal spaces for such embeddings.
△ Less
Submitted 16 July, 2020;
originally announced July 2020.
-
Quantum Gates on Individually-Addressed Atomic Qubits Subject to Noisy Transverse Motion
Authors:
M. Cetina,
L. N. Egan,
C. A. Noel,
M. L. Goldman,
A. R. Risinger,
D. Zhu,
D. Biswas,
C. Monroe
Abstract:
Individual trapped atomic qubits represent one of the most promising technologies to scale quantum computers, owing to their negligible idle errors and the ability to implement a full set of reconfigurable gate operations via focused optical fields. However, the fidelity of quantum gate operations can be limited by weak confinement of the atoms transverse to the laser. We present measurements of t…
▽ More
Individual trapped atomic qubits represent one of the most promising technologies to scale quantum computers, owing to their negligible idle errors and the ability to implement a full set of reconfigurable gate operations via focused optical fields. However, the fidelity of quantum gate operations can be limited by weak confinement of the atoms transverse to the laser. We present measurements of this effect by performing individually-addressed entangling gates in chains of up to 25 trapped atomic ions that are weakly confined along the chain axis. We present a model that accurately describes the observed decoherence from the residual heating of the ions caused by noisy electric fields. We propose to suppress these effects through the use of ancilla ions interspersed in the chain to sympathetically cool the qubit ions throughout a quantum circuit.
△ Less
Submitted 13 July, 2020;
originally announced July 2020.
-
Multi-beam Energy Moments of Multibeam Particle Velocity Distributions
Authors:
M. V. Goldman,
D. L. Newman,
J. P. Eastwood,
G. Lapenta
Abstract:
High resolution electron and ion velocity distributions, f(v), which consist of N effectively disjoint beams, have been measured by NASA's Magnetospheric Multi-Scale Mission (MMS) observatories and in reconnection simulations. Commonly used standard velocity moments generally assume a single mean-flow-velocity for the entire distribution, which can lead to counterintuitive results for a multibeam…
▽ More
High resolution electron and ion velocity distributions, f(v), which consist of N effectively disjoint beams, have been measured by NASA's Magnetospheric Multi-Scale Mission (MMS) observatories and in reconnection simulations. Commonly used standard velocity moments generally assume a single mean-flow-velocity for the entire distribution, which can lead to counterintuitive results for a multibeam f(v). An example is the (false) standard thermal energy moment of a pair of equal and opposite cold particle beams, which is nonzero even though each beam has zero thermal energy. By contrast, a multibeam moment of two or more beams has no false thermal energy. A multibeam moment is obtained by taking a standard moment of each beam and then summing over beams. In this paper we will generalize these notions, explore their consequences and apply them to an f(v) which is sum of tri-Maxwellians. Both standard and multibeam energy moments have coherent and incoherent forms. Examples of incoherent moments are the thermal energy density, the pressure and the thermal energy flux (enthalpy flux plus heat flux). Corresponding coherent moments are the bulk kinetic energy density, the RAM pressure and the bulk kinetic energy flux. The false part of an incoherent moment is defined as the difference between the standard incoherent moment and the corresponding multibeam moment. The sum of a pair of corresponding coherent and incoherent moments will be called the undecomposed moment. Undecomposed moments are independent of whether the sum is standard or multibeam and therefore have advantages when studying moments of measured f(v).
△ Less
Submitted 18 May, 2020;
originally announced May 2020.
-
An $ε$--regularity result for optimal transport maps between continuous densities
Authors:
Michael Goldman
Abstract:
The aim of this short note is to extend the recent variational proof of partial regularity for optimal transport maps to the case of continuous densities.
The aim of this short note is to extend the recent variational proof of partial regularity for optimal transport maps to the case of continuous densities.
△ Less
Submitted 20 November, 2020; v1 submitted 14 May, 2020;
originally announced May 2020.
-
Proper actions of discrete groups of affine transformations
Authors:
Jeffrey Danciger,
Todd A. Drumm,
William M. Goldman,
Ilia Smilga
Abstract:
In the early 1980's Margulis startled the world by showing the existence of proper affine actions of free groups on 3-space, answering a provocative and suggestive question Milnor posed in 1977. In this paper we discuss the historical background motivating this question, recent progress on this subject, and future directions inspired by this discovery.
In the early 1980's Margulis startled the world by showing the existence of proper affine actions of free groups on 3-space, answering a provocative and suggestive question Milnor posed in 1977. In this paper we discuss the historical background motivating this question, recent progress on this subject, and future directions inspired by this discovery.
△ Less
Submitted 21 February, 2020;
originally announced February 2020.
-
Online information of vaccines: information quality is an ethical responsibility of search engines
Authors:
Pietro Ghezzi,
Peter G Bannister,
Gonzalo Casino,
Alessia Catalani,
Michel Goldman,
Jessica Morley,
Marie Neunez,
Andreu Prados,
Mariarosaria Taddeo,
Tania Vanzolini,
Luciano Floridi
Abstract:
The fact that internet companies may record our personal data and track our online behavior for commercial or political purpose has emphasized aspects related to online privacy. This has also led to the development of search engines that promise no tracking and privacy. Search engines also have a major role in spreading low-quality health information such as that of anti-vaccine websites. This stu…
▽ More
The fact that internet companies may record our personal data and track our online behavior for commercial or political purpose has emphasized aspects related to online privacy. This has also led to the development of search engines that promise no tracking and privacy. Search engines also have a major role in spreading low-quality health information such as that of anti-vaccine websites. This study investigates the relationship between search engines' approach to privacy and the scientific quality of the information they return. We analyzed the first 30 webpages returned searching 'vaccines autism' in English, Spanish, Italian and French. The results show that alternative search engines (Duckduckgo, Ecosia, Qwant, Swisscows and Mojeek) may return more anti-vaccine pages (10 to 53 percent) than Google.com (zero). Some localized versions of Google, however, returned more anti-vaccine webpages (up to 10 percent) than Google.com. Our study suggests that designing a search engine that is privacy savvy and avoids issues with filter bubbles that can result from user tracking is necessary but insufficient; instead, mechanisms should be developed to test search engines from the perspective of information quality (particularly for health-related webpages), before they can be deemed trustworthy providers of public health information.
△ Less
Submitted 2 December, 2019;
originally announced December 2019.
-
Local regimes of turbulence in 3D magnetic reconnection
Authors:
G. Lapenta,
F. Pucci,
M. V. Goldman,
D. L. Newman
Abstract:
The process of magnetic reconnection when studied in Nature or when modeled in 3D simulations differs in one key way from the standard 2D paradigmatic cartoon: it is accompanied by much fluctuations in the electromagnetic fields and plasma properties. We developed a diagnostics to study the spectrum of fluctuations in the various regions around a reconnection site. We define the regions in terms o…
▽ More
The process of magnetic reconnection when studied in Nature or when modeled in 3D simulations differs in one key way from the standard 2D paradigmatic cartoon: it is accompanied by much fluctuations in the electromagnetic fields and plasma properties. We developed a diagnostics to study the spectrum of fluctuations in the various regions around a reconnection site. We define the regions in terms of the local value of the flux function that determines the distance form the reconnection site, with positive values in the outflow and negative values in the inflow. We find that fluctuations belong to two very different regimes depending on the local plasma beta (defined as the ratio of plasma and magnetic pressure). The first regime develops in the reconnection outflows where beta is high and is characterized by a strong link between plasma and electromagnetic fluctuations leading to momentum and energy exchanges via anomalous viscosity and resistivity. But there is a second, low beta regime: it develops in the inflow and in the region around the separatrix surfaces, including the reconnection electron diffusion region itself. It is remarkable that this low beta plasma, where the magnetic pressure dominates, remain laminar even though the electromagnetic fields are turbulent.
△ Less
Submitted 26 October, 2019;
originally announced October 2019.
-
Characterizing magnetic reconnection regions using Gaussian mixture models on particle velocity distributions
Authors:
Romain Dupuis,
Martin V. Goldman,
David L. Newman,
Jorge Amaya,
Giovanni Lapenta
Abstract:
We present a method based on unsupervised machine learning to identify regions of interest using particle velocity distributions as a signature pattern. An automatic density estimation technique is applied to particle distributions provided by PIC simulations to study magnetic reconnection. The key components of the method involve: i) a Gaussian mixture model determining the presence of a given nu…
▽ More
We present a method based on unsupervised machine learning to identify regions of interest using particle velocity distributions as a signature pattern. An automatic density estimation technique is applied to particle distributions provided by PIC simulations to study magnetic reconnection. The key components of the method involve: i) a Gaussian mixture model determining the presence of a given number of subpopulations within an overall population, and ii) a model selection technique with Bayesian Information Criterion to estimate the appropriate number of subpopulations. Thus, this method identifies automatically the presence of complex distributions, such as beams or other non-Maxwellian features, and can be used as a detection algorithm able to identify reconnection regions. The approach is demonstrated for specific double Harris sheet simulations but it can in principle be applied to any other type of simulation and observational data on the particle distribution function.
△ Less
Submitted 22 October, 2019;
originally announced October 2019.
-
The mapping class group action on SU(3)-character varieties
Authors:
William M. Goldman,
Sean Lawton,
Eugene Z. Xia
Abstract:
Let $Σ$ be a compact orientable surface of genus $g=1$ with $n=1$ boundary component. The mapping class group $Γ$ of $Σ$ acts on the SU(3)-character variety of $Σ$. We show that the action is ergodic with respect to the natural symplectic measure on the character variety.
Let $Σ$ be a compact orientable surface of genus $g=1$ with $n=1$ boundary component. The mapping class group $Γ$ of $Σ$ acts on the SU(3)-character variety of $Σ$. We show that the action is ergodic with respect to the natural symplectic measure on the character variety.
△ Less
Submitted 2 May, 2020; v1 submitted 24 September, 2019;
originally announced September 2019.
-
A variational approach to regularity theory in optimal transportation
Authors:
Michael Goldman
Abstract:
This paper describes recent results obtained in collaboration with M. Huesmann and F. Otto on the regularity of optimal transport maps. The main result is a quantitative version of the well-known fact that the linearization of the Monge-Amp{è}re equation around the identity is the Poisson equation. We present two applications of this result. The first one is a variational proof of the partial regu…
▽ More
This paper describes recent results obtained in collaboration with M. Huesmann and F. Otto on the regularity of optimal transport maps. The main result is a quantitative version of the well-known fact that the linearization of the Monge-Amp{è}re equation around the identity is the Poisson equation. We present two applications of this result. The first one is a variational proof of the partial regularity theorem of Figalli and Kim and the second is the rigorous validation of some predictions made by Carraciolo and al. on the structure of the optimal transport maps in matching problems.
△ Less
Submitted 12 July, 2019;
originally announced July 2019.
-
Matching on What Matters: A Pseudo-Metric Learning Approach to Matching Estimation in High Dimensions
Authors:
Gentry Johnson,
Brian Quistorff,
Matt Goldman
Abstract:
When pre-processing observational data via matching, we seek to approximate each unit with maximally similar peers that had an alternative treatment status--essentially replicating a randomized block design. However, as one considers a growing number of continuous features, a curse of dimensionality applies making asymptotically valid inference impossible (Abadie and Imbens, 2006). The alternative…
▽ More
When pre-processing observational data via matching, we seek to approximate each unit with maximally similar peers that had an alternative treatment status--essentially replicating a randomized block design. However, as one considers a growing number of continuous features, a curse of dimensionality applies making asymptotically valid inference impossible (Abadie and Imbens, 2006). The alternative of ignoring plausibly relevant features is certainly no better, and the resulting trade-off substantially limits the application of matching methods to "wide" datasets. Instead, Li and Fu (2017) recasts the problem of matching in a metric learning framework that maps features to a low-dimensional space that facilitates "closer matches" while still capturing important aspects of unit-level heterogeneity. However, that method lacks key theoretical guarantees and can produce inconsistent estimates in cases of heterogeneous treatment effects. Motivated by straightforward extension of existing results in the matching literature, we present alternative techniques that learn latent matching features through either MLPs or through siamese neural networks trained on a carefully selected loss function. We benchmark the resulting alternative methods in simulations as well as against two experimental data sets--including the canonical NSW worker training program data set--and find superior performance of the neural-net-based methods.
△ Less
Submitted 28 May, 2019;
originally announced May 2019.
-
Quantitative linearization results for the Monge-Ampère equation
Authors:
Michael Goldman,
Martin Huesmann,
Felix Otto
Abstract:
This paper is about quantitative linearization results for the Monge-Ampère equation with rough data. We develop a large-scale regularity theory and prove that if a measure $μ$ is close to the Lebesgue measure in Wasserstein distance at all scales, then the displacement of the macroscopic optimal coupling is quantitatively close at all scales to the gradient of the solution of the corresponding Po…
▽ More
This paper is about quantitative linearization results for the Monge-Ampère equation with rough data. We develop a large-scale regularity theory and prove that if a measure $μ$ is close to the Lebesgue measure in Wasserstein distance at all scales, then the displacement of the macroscopic optimal coupling is quantitatively close at all scales to the gradient of the solution of the corresponding Poisson equation. The main ingredient we use is a harmonic approximation result for the optimal transport plan between arbitrary measures. This is used in a Campanato iteration which transfers the information through the scales.
△ Less
Submitted 30 April, 2021; v1 submitted 23 May, 2019;
originally announced May 2019.
-
Non-convex functionals penalizing simultaneous oscillations along independent directions: rigidity estimates
Authors:
Michael Goldman,
B. Merlet
Abstract:
We study a family of non-convex functionals $\{\mathcal{E}\}$ on the space of measurable functions$u: Ω_1\timesΩ_2 \subset \mathbb{R}^{n_1}\times\mathbb{R}^{n_2} \to \mathbb{R}$. These functionals vanish on the non-convex subset $S(Ω_1\timesΩ_2)$ formed by functions of the form $u(x_1,x_2)=u_1(x_1)$ or $u(x_1,x_2)=u_2(x_2)$. We investigate under which conditions the converse implication "…
▽ More
We study a family of non-convex functionals $\{\mathcal{E}\}$ on the space of measurable functions$u: Ω_1\timesΩ_2 \subset \mathbb{R}^{n_1}\times\mathbb{R}^{n_2} \to \mathbb{R}$. These functionals vanish on the non-convex subset $S(Ω_1\timesΩ_2)$ formed by functions of the form $u(x_1,x_2)=u_1(x_1)$ or $u(x_1,x_2)=u_2(x_2)$. We investigate under which conditions the converse implication "$\mathcal{E}(u) = 0 \Rightarrow u \in S(Ω_1\timesΩ_2)$" holds. In particular, we show that the answer depends strongly on the smoothness of u. We also obtain quantitative versions of this implication by proving that (at least for some parameters) $\mathcal{E}(u)$ controls in a strong sense the distance of $u$ to $S(Ω_1\timesΩ_2)$.
△ Less
Submitted 20 May, 2019;
originally announced May 2019.
-
Electrolytic Conversion of Bicarbonate into CO in a Flow Cell
Authors:
Tengfei Li,
Eric W. Lees,
Maxwell Goldman,
Danielle A. Salvatore,
David M. Weekes,
Curtis P. Berlinguette
Abstract:
Electrochemical CO2 reduction offers a method to use renewable electricity to convert CO2 into CO and other carbon-based chemical building blocks. While nearly all studies rely on a CO2 feed, we show herein that aqueous bicarbonate solutions can also be electrochemically converted into CO gas at meaningful rates in a flow cell. We achieved this result in a flow cell containing a bipolar membrane (…
▽ More
Electrochemical CO2 reduction offers a method to use renewable electricity to convert CO2 into CO and other carbon-based chemical building blocks. While nearly all studies rely on a CO2 feed, we show herein that aqueous bicarbonate solutions can also be electrochemically converted into CO gas at meaningful rates in a flow cell. We achieved this result in a flow cell containing a bipolar membrane (BPM) and a silver nanoparticle catalyst on a porous carbon support. Electrolysis upon a N2-saturated 3.0-M potassium bicarbonate electrolyte solution yields CO with a faradaic efficiency (F.E.CO) of 81% at 25 mA cm-2 and 37% at 100 mA cm-2. This output is comparable to the analogous experiment where the electrolyte is saturated with gaseous CO2 (faradaic efficiency for CO is 78% at 25 mA cm-2 and 35% at 100 mA cm-2). The H+ flux from the BPM is critical to this chemistry in that it reacts with the bicarbonate feed to generate CO2, which is then reduced to CO at the gas diffusion electrode. These results are important in that they show that the addition of gaseous CO2 to bicarbonate electrolytes is not necessary in order to obtain reduced carbon products with a flow cell architecture. This process offers a means of using electrolysis to bypass the thermally-intensive step of extracting CO2 from bicarbonate solutions generated in carbon capture schemes.
△ Less
Submitted 11 May, 2019;
originally announced May 2019.
-
Learn Stereo, Infer Mono: Siamese Networks for Self-Supervised, Monocular, Depth Estimation
Authors:
Matan Goldman,
Tal Hassner,
Shai Avidan
Abstract:
The field of self-supervised monocular depth estimation has seen huge advancements in recent years. Most methods assume stereo data is available during training but usually under-utilize it and only treat it as a reference signal. We propose a novel self-supervised approach which uses both left and right images equally during training, but can still be used with a single input image at test time,…
▽ More
The field of self-supervised monocular depth estimation has seen huge advancements in recent years. Most methods assume stereo data is available during training but usually under-utilize it and only treat it as a reference signal. We propose a novel self-supervised approach which uses both left and right images equally during training, but can still be used with a single input image at test time, for monocular depth estimation. Our Siamese network architecture consists of two, twin networks, each learns to predict a disparity map from a single image. At test time, however, only one of these networks is used in order to infer depth. We show state-of-the-art results on the standard KITTI Eigen split benchmark as well as being the highest scoring self-supervised method on the new KITTI single view benchmark. To demonstrate the ability of our method to generalize to new data sets, we further provide results on the Make3D benchmark, which was not used during training.
△ Less
Submitted 1 May, 2019;
originally announced May 2019.
-
A violin sonata for reconnection
Authors:
G. Lapenta,
F. Pucci,
M. V. Goldman,
D. L. Newman
Abstract:
The process of magnetic reconnection when studied in Nature or when modeled in 3D simulations differs in one key way from the standard 2D paradigmatic cartoon: it is accompanied by much fluctuations in the electromagnetic fields and plasma properties. We developed a new diagnostics, the topographical fluctuations analysis (TFA) to study the spectrum of fluctuations in the various regions around a…
▽ More
The process of magnetic reconnection when studied in Nature or when modeled in 3D simulations differs in one key way from the standard 2D paradigmatic cartoon: it is accompanied by much fluctuations in the electromagnetic fields and plasma properties. We developed a new diagnostics, the topographical fluctuations analysis (TFA) to study the spectrum of fluctuations in the various regions around a reconnection site. We find that fluctuations belong to two very different regimes. The first regime is better known, it develops in the reconnection outflows and is characterized by a strong link between plasma and electromagnetic fluctuations leading to momentum and energy exchanges via anomalous viscosity and resistivity. But there is a second, new, regime: it develops in the inflow and in the region around the separatrix surfaces, including the reconnection diffusion region itself. In this new regime the plasma remains laminar but the electromagnetic fields fluctuates strongly. We present an analogy with the smooth continuous motion of the bow of a violin producing the vibrations of the strings to emit music.
△ Less
Submitted 3 April, 2019;
originally announced April 2019.
-
Properties of turbulence in the reconnection exhaust: numerical simulations compared with observations
Authors:
F. Pucci,
S. Servidio,
L. Sorriso-Valvo,
V. Olshevsky,
W. H. Matthaeus,
F. Malara,
M. V. Goldman,
D. L. Newman,
G. Lapenta
Abstract:
The properties of the turbulence which develops in the outflows of magnetic reconnection have been investigated using self-consistent plasma simulations, in three dimensions. As commonly observed in space plasmas, magnetic reconnection is characterized by the presence of turbulence. Here we provide a direct comparison of our simulations with reported observations of reconnection events in the magn…
▽ More
The properties of the turbulence which develops in the outflows of magnetic reconnection have been investigated using self-consistent plasma simulations, in three dimensions. As commonly observed in space plasmas, magnetic reconnection is characterized by the presence of turbulence. Here we provide a direct comparison of our simulations with reported observations of reconnection events in the magnetotail investigating the properties of the electromagnetic field and the energy conversion mechanisms. In particular, simulations show the development of a turbulent cascade consistent with spacecraft observations, statistics of the the dissipation mechanisms in the turbulent outflows similar to the one observed in reconnection jets in the magnetotail, and that the properties of turbulence vary as a function of the distance from the reconnecting X-line.
△ Less
Submitted 31 October, 2018;
originally announced November 2018.
-
Generation of turbulence in colliding reconnection jets
Authors:
F. Pucci,
W. H. Matthaeus,
A. Chasapis,
S. Servidio,
L. Sorriso-Valvo,
V. Olshevsky,
D. L. Newman,
M. V. Goldman,
G. Lapenta
Abstract:
The collision of magnetic reconnection jets is studied by means of a three dimensional numerical simulation at kinetic scale, in the presence of a strong guide field. We show that turbulence develops due to the jets collision producing several current sheets in reconnection outflows, aligned with the guide field direction. The turbulence is mainly two-dimensional, with stronger gradients in the pl…
▽ More
The collision of magnetic reconnection jets is studied by means of a three dimensional numerical simulation at kinetic scale, in the presence of a strong guide field. We show that turbulence develops due to the jets collision producing several current sheets in reconnection outflows, aligned with the guide field direction. The turbulence is mainly two-dimensional, with stronger gradients in the plane perpendicular to the guide field and a low wave-like activity in the parallel direction. First, we provide a numerical method to isolate the central turbulent region. Second, we analyze spatial second-order structure function and prove that turbulence is confined in this region. Finally, we compute local magnetic and electric frequency spectra, finding a trend in the sub-ion range that differs from typical cases for which the Taylor hypothesis is valid, as well as wave activity in the range between ion and electron cyclotron frequencies. Our results are relevant to understand observations of reconnection jets collisions in space plasmas.
△ Less
Submitted 31 October, 2018;
originally announced October 2018.
-
Infinite-randomness fixed point of the quantum superconductor-metal transitions in amorphous thin films
Authors:
Nicholas A. Lewellyn,
Ilana M. Percher,
JJ Nelson,
Javier Garcia-Barriocanal,
Irina Volotsenko,
Aviad Frydman,
Thomas Vojta,
Allen M. Goldman
Abstract:
The magnetic-field-tuned quantum superconductor-insulator transitions of disordered amorphous indium oxide films are a paradigm in the study of quantum phase transitions, and exhibit power-law scaling behavior. For superconducting indium oxide films with low disorder, such as the ones reported on here, the high-field state appears to be a quantum-corrected metal. Resistance data across the superco…
▽ More
The magnetic-field-tuned quantum superconductor-insulator transitions of disordered amorphous indium oxide films are a paradigm in the study of quantum phase transitions, and exhibit power-law scaling behavior. For superconducting indium oxide films with low disorder, such as the ones reported on here, the high-field state appears to be a quantum-corrected metal. Resistance data across the superconductor-metal transition in these films are shown here to obey an activated scaling form appropriate to a quantum phase transition controlled by an infinite randomness fixed point in the universality class of the random transverse-field Ising model. Collapse of the field-dependent resistance vs. temperature data is obtained using an activated scaling form appropriate to this universality class, using values determined through a modified form of power-law scaling analysis. This exotic behavior of films exhibiting a superconductor-metal transition is caused by the dissipative dynamics of superconducting rare regions immersed in a metallic matrix, as predicted by a recent renormalization group theory. The smeared crossing points of isotherms observed are due to corrections to scaling which are expected near an infinite randomness critical point, where the inverse disorder strength acts as an irrelevant scaling variable.
△ Less
Submitted 8 January, 2019; v1 submitted 7 September, 2018;
originally announced September 2018.
-
A two-point function approach to connectedness of drops in convex potentials
Authors:
Michael Goldman,
Guido De Philippis
Abstract:
We establish connectedness of volume constrained minimisers of energies involving surface tensions and convex potentials. By a previous result of McCann, this implies that minimisers are convex in dimension two. This positively answers an old question of Almgren. We also prove convexity of minimisers when the volume constraint is dropped. Our proof is based on the introduction of a new "two-point…
▽ More
We establish connectedness of volume constrained minimisers of energies involving surface tensions and convex potentials. By a previous result of McCann, this implies that minimisers are convex in dimension two. This positively answers an old question of Almgren. We also prove convexity of minimisers when the volume constraint is dropped. Our proof is based on the introduction of a new "two-point function" which measures the lack of convexity and which gives rise to a negative second variation of the energy.
△ Less
Submitted 5 September, 2018;
originally announced September 2018.
-
A large-scale regularity theory for the Monge-Ampere equation with rough data and application to the optimal matching problem
Authors:
Michael Goldman,
Martin Huesmann,
Felix Otto
Abstract:
The aim of this paper is to obtain quantitative bounds for solutions to the optimal matching problem in dimension two. These bounds show that up to a logarithmically divergent shift, the optimal transport maps are close to be the identity at every scale. These bounds allow us to pass to the limit as the system size goes to infinity and construct a locally optimal coupling between the Lebesgue meas…
▽ More
The aim of this paper is to obtain quantitative bounds for solutions to the optimal matching problem in dimension two. These bounds show that up to a logarithmically divergent shift, the optimal transport maps are close to be the identity at every scale. These bounds allow us to pass to the limit as the system size goes to infinity and construct a locally optimal coupling between the Lebesgue measure and the Poisson point process which retains the stationarity properties of the Poisson point process only at the level of second-order differences. Our quantitative bounds are obtained through a Campanato iteration scheme based on a deterministic and a stochastic ingredient. The deterministic part, which can be seen as our main contribution, is a regularity result for Monge-Ampère equations with rough right-hand side. Since we believe that it could be useful in other contexts, we prove it for general space dimensions. The stochastic part is a concentration result for the optimal matching problem which builds on previous work by Ambrosio, Stra and Trevisan.
△ Less
Submitted 28 August, 2018;
originally announced August 2018.
-
Non-linear Waves and Instabilities Leading to Secondary Reconnection in Reconnection Outflows
Authors:
Giovanni Lapenta,
Francesco Pucci,
Vyacheslav Olshevsky,
Sergio Servidio,
Luca Sorriso-Valvo,
David L. Newman,
Martin Goldman
Abstract:
Reconnection outflows are regions of intense recent scrutiny, from in situ observations and from simulations. These regions are host to a variety of instabilities and intense energy exchanges, often even superior to the main reconnection site. We report here a number of results drawn from investigation of simulations. First, the outflows are observed to become unstable to drift instabilities. Seco…
▽ More
Reconnection outflows are regions of intense recent scrutiny, from in situ observations and from simulations. These regions are host to a variety of instabilities and intense energy exchanges, often even superior to the main reconnection site. We report here a number of results drawn from investigation of simulations. First, the outflows are observed to become unstable to drift instabilities. Second, these instabilities lead to the formation of secondary reconnection sites. Third, the secondary processes are responsible for large energy exchanges and particle energization. Finally, the particle distribution function are modified to become non-Maxwellian and include multiple interpenetrating populations.
△ Less
Submitted 26 August, 2018;
originally announced August 2018.
-
Optical Control of a Single Nuclear Spin in the Solid State
Authors:
M. L. Goldman,
T. L. Patti,
D. Levonian,
S. F. Yelin,
M. D. Lukin
Abstract:
We demonstrate a novel method for coherent optical manipulation of individual nuclear spins in the solid state, mediated by the electronic states of a proximal quantum emitter. Specifically, using the nitrogen-vacancy (NV) color center in diamond, we demonstrate control of a proximal $^{14}$N nuclear spin via an all-optical Raman technique. We evaluate the extent to which the intrinsic physical pr…
▽ More
We demonstrate a novel method for coherent optical manipulation of individual nuclear spins in the solid state, mediated by the electronic states of a proximal quantum emitter. Specifically, using the nitrogen-vacancy (NV) color center in diamond, we demonstrate control of a proximal $^{14}$N nuclear spin via an all-optical Raman technique. We evaluate the extent to which the intrinsic physical properties of the NV center limit the performance of coherent control, and find that it is ultimately constrained by the relative rates of transverse hyperfine coupling and radiative decay in the NV center's excited state. Possible extensions and applications to other color centers are discussed.
△ Less
Submitted 13 August, 2018;
originally announced August 2018.
