Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
17works
0followers
20topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

17 published item(s)

preprint2024arXiv

Approximation Algorithms for the Weighted Nash Social Welfare via Convex and Non-Convex Programs

In an instance of the weighted Nash Social Welfare problem, we are given a set of $m$ indivisible items, $\mathscr{G}$, and $n$ agents, $\mathscr{A}$, where each agent $i \in \mathscr{A}$ has a valuation $v_{ij}\geq 0$ for each item $j\in \mathscr{G}$. In addition, every agent $i$ has a non-negative weight $w_i$ such that the weights collectively sum up to $1$. The goal is to find an assignment $σ:\mathscr{G}\rightarrow \mathscr{A}$ that maximizes $\prod_{i\in \mathscr{A}} \left(\sum_{j\in σ^{-1}(i)} v_{ij}\right)^{w_i}$, the product of the weighted valuations of the players. When all the weights equal $\frac1n$, the problem reduces to the classical Nash Social Welfare problem, which has recently received much attention. In this work, we present a $5\cdot\exp\left(2\cdot D_{\text{KL}}(\mathbf{w}\, ||\, \frac{\vec{\mathbf{1}}}{n})\right) = 5\cdot\exp\left(2\log{n} + 2\sum_{i=1}^n w_i \log{w_i}\right)$-approximation algorithm for the weighted Nash Social Welfare problem, where $D_{\text{KL}}(\mathbf{w}\, ||\, \frac{\vec{\mathbf{1}}}{n})$ denotes the KL-divergence between the distribution induced by $\mathbf{w}$ and the uniform distribution on $[n]$. We show a novel connection between the convex programming relaxations for the unweighted variant of Nash Social Welfare presented in \cite{cole2017convex, anari2017nash}, and generalize the programs to two different mathematical programs for the weighted case. The first program is convex and is necessary for computational efficiency, while the second program is a non-convex relaxation that can be rounded efficiently. The approximation factor derives from the difference in the objective values of the convex and non-convex relaxation.

preprint2022arXiv

$λ_\infty$ & Maximum Variance Embedding: Measuring and Optimizing Connectivity of A Graph Metric

Bobkov, Houdré, and the last author [2000] introduced a Poincaré-type functional parameter, $λ_\infty$, of a graph and related it to connectivity of the graph via Cheeger-type inequalities. A work by the second author, Raghavendra, and Vempala [2013] related the complexity of $λ_\infty$ to the so-called small-set expansion (SSE) problem and further set forth the desiderata for NP-hardness of this optimization problem. We confirm the conjecture that computing $λ_\infty$ is NP-hard for weighted trees. Beyond measuring connectivity in many applications we want to optimize it. This, via convex duality, leads to a problem in machine learning known as the Maximum Variance Embedding (MVE). The output is a function from vertices to a low dim Euclidean space, subject to bounds on Euclidean distances between neighbors. The objective is to maximize output variance. Special cases of MVE into $n$ and $1$ dims lead to absolute algebraic connectivity [1990] and spread constant [1998], that measure connectivity of the graph and its Cartesian $n$-power, respectively. MVE has other applications in measuring diffusion speed and robustness of networks, clustering, and dimension reduction. We show that computing MVE in tree-width dims is NP-hard, while only one additional dim beyond width of a given tree-decomposition makes the problem in P. We show that MVE of a tree in 2 dims defines a non-convex yet benign optimization landscape, i.e., local=global optima. We further develop a linear time combinatorial algorithm for this case. Finally, we denote approximate Maximum Variance Embedding is tractable in significantly lower dims. For trees and general graphs, for which Maximum Variance Embedding cannot be solved in less than $2$ and $Ω(n)$ dims, we provide $1+\varepsilon$ approximation algorithms for embedding into $1$ and $O(\log n /\varepsilon^2)$ dims, respectively.

preprint2022arXiv

Bistable colloidal orientation near a charged surface

Anisotropic particles oriented in a specific direction can act as artificial atoms and molecules, and their controlled assembly can result in a wide variety of ordered structures. Towards this, we demonstrate the orientation transitions of uncharged peanut-shaped polystyrene colloids, suspended in a non-ionic aprotic polar solvent, near a flat surface whose potential is static or time-varying. The charged surface is coated with an insulating dielectric layer to suppress electric currents. The transition between several orientation states such as random, normal or parallel orientation with respect to the surface, is examined for two different colloid sizes at low-frequency ($\sim 10-350$ kHz) or static fields, and at small electric potentials. In time-varying (AC) field, a detailed phase diagram in the potential-frequency plane indicating the transition between particles parallel or normal to the surface is reported. We next present the first study of orientation switching in static (DC) fields, where no electro-osmotic or other flow is present. A reversible change between the two colloidal states is explained by a theory showing that the sum of electrostatic and gravitational energies of the colloid is bistable. The number of colloids in each of the two states depends on the external potential, particle and solvent permittivities, particle aspect ratio, and distance from the electrode.

preprint2022arXiv

Constant-Factor Approximation Algorithms for Socially Fair $k$-Clustering

We study approximation algorithms for the socially fair $(\ell_p, k)$-clustering problem with $m$ groups, whose special cases include the socially fair $k$-median ($p=1$) and socially fair $k$-means ($p=2$) problems. We present (1) a polynomial-time $(5+2\sqrt{6})^p$-approximation with at most $k+m$ centers (2) a $(5+2\sqrt{6}+ε)^p$-approximation with $k$ centers in time $n^{2^{O(p)}\cdot m^2}$, and (3) a $(15+6\sqrt{6})^p$ approximation with $k$ centers in time $k^{m}\cdot\text{poly}(n)$. The first result is obtained via a refinement of the iterative rounding method using a sequence of linear programs. The latter two results are obtained by converting a solution with up to $k+m$ centers to one with $k$ centers using sparsification methods for (2) and via an exhaustive search for (3). We also compare the performance of our algorithms with existing bicriteria algorithms as well as exactly $k$ center approximation algorithms on benchmark datasets, and find that our algorithms also outperform existing methods in practice.

preprint2022arXiv

Determinant Maximization via Matroid Intersection Algorithms

Determinant maximization problem gives a general framework that models problems arising in as diverse fields as statistics \cite{pukelsheim2006optimal}, convex geometry \cite{Khachiyan1996}, fair allocations\linebreak \cite{anari2016nash}, combinatorics \cite{AnariGV18}, spectral graph theory \cite{nikolov2019proportional}, network design, and random processes \cite{kulesza2012determinantal}. In an instance of a determinant maximization problem, we are given a collection of vectors $U=\{v_1,\ldots, v_n\} \subset \RR^d$, and a goal is to pick a subset $S\subseteq U$ of given vectors to maximize the determinant of the matrix $\sum_{i\in S} v_i v_i^\top $. Often, the set $S$ of picked vectors must satisfy additional combinatorial constraints such as cardinality constraint $\left(|S|\leq k\right)$ or matroid constraint ($S$ is a basis of a matroid defined on the vectors). In this paper, we give a polynomial-time deterministic algorithm that returns a $r^{O(r)}$-approximation for any matroid of rank $r\leq d$. This improves previous results that give $e^{O(r^2)}$-approximation algorithms relying on $e^{O(r)}$-approximate \emph{estimation} algorithms \cite{NikolovS16,anari2017generalization,AnariGV18,madan2020maximizing} for any $r\leq d$. All previous results use convex relaxations and their relationship to stable polynomials and strongly log-concave polynomials. In contrast, our algorithm builds on combinatorial algorithms for matroid intersection, which iteratively improve any solution by finding an \emph{alternating negative cycle} in the \emph{exchange graph} defined by the matroids. While the $\det(.)$ function is not linear, we show that taking appropriate linear approximations at each iteration suffice to give the improved approximation algorithm.

preprint2022arXiv

Tropicalization of Graph Profiles

A graph profile records all possible densities of a fixed finite set of graphs. Profiles can be extremely complicated; for instance the full profile of any triple of connected graphs is not known, and little is known about hypergraph profiles. We introduce the tropicalization of graph and hypergraph profiles. Tropicalization is a well-studied operation in algebraic geometry, which replaces a variety (the set of real or complex solutions to a finite set of algebraic equations) with its "combinatorial shadow". We prove that the tropicalization of a graph profile is a closed convex cone, which still captures interesting combinatorial information. We explicitly compute these tropicalizations for arbitrary sets of complete and star hypergraphs. We show they are rational polyhedral cones even though the corresponding profiles are not even known to be semialgebraic in some of these cases. We then use tropicalization to prove strong restrictions on the power of the sums of squares method, equivalently Cauchy-Schwarz calculus, to test (which is weaker than certification) the validity of graph density inequalities. In particular, we show that sums of squares cannot test simple binomial graph density inequalities, or even their approximations. Small concrete examples of such inequalities are presented, and include the famous Blakley-Roy inequalities for paths of odd length. As a consequence, these simple inequalities cannot be written as a rational sum of squares of graph densities.

preprint2021arXiv

Adaptive Bin Packing with Overflow

Motivated by bursty bandwidth allocation and by the allocation of virtual machines to servers in the cloud, we consider the online problem of packing items with random sizes into unit-capacity bins. Items arrive sequentially, but upon arrival an item's actual size is unknown; only its probabilistic information is available to the decision maker. Without knowing this size, the decision maker must irrevocably pack the item into an available bin or place it in a new bin. Once packed in a bin, the decision maker observes the item's actual size, and overflowing the bin is a possibility. An overflow incurs a large penalty cost and the corresponding bin is unusable for the rest of the process. In practical terms, this overflow models delayed services, failure of servers, and/or loss of end-user goodwill. The objective is to minimize the total expected cost given by the sum of the number of opened bins and the overflow penalty cost. We present an online algorithm with expected cost at most a constant factor times the cost incurred by the optimal packing policy when item sizes are drawn from an i.i.d. sequence of unknown length. We give a similar result when item size distributions are exponential with arbitrary rates. We also study the offline model, where distributions are known in advance but must be packed sequentially. We construct a soft-capacity PTAS for this problem, and show that the complexity of computing the optimal offline cost is $\#\mathbf{P}$-hard. Finally, we provide an empirical study of our online algorithm's performance.

preprint2021arXiv

Dynamic Resource Allocation in the Cloud with Near-Optimal Efficiency

Cloud computing has motivated renewed interest in resource allocation problems with new consumption models. A common goal is to share a resource, such as CPU or I/O bandwidth, among distinct users with different demand patterns as well as different quality of service requirements. To ensure these service requirements, cloud offerings often come with a service level agreement (SLA) between the provider and the users. An SLA specifies the amount of a resource a user is entitled to utilize. In many cloud settings, providers would like to operate resources at high utilization while simultaneously respecting individual SLAs. There is typically a tradeoff between these two objectives; for example, utilization can be increased by shifting away resources from idle users to "scavenger" workload, but with the risk of the former then becoming active again. We study this fundamental tradeoff by formulating a resource allocation model that captures basic properties of cloud computing systems, including SLAs, highly limited feedback about the state of the system, and variable and unpredictable input sequences. Our main result is a simple and practical algorithm that achieves near-optimal performance on the above two objectives. First, we guarantee nearly optimal utilization of the resource even if compared to the omniscient offline dynamic optimum. Second, we simultaneously satisfy all individual SLAs up to a small error. The main algorithmic tool is a multiplicative weight update algorithm, and a primal-dual argument to obtain its guarantees. We also provide numerical validation on real data to demonstrate the performance of our algorithm in practical applications.

preprint2021arXiv

Effect of discrete breathers on the specific heat of a nonlinear chain

A nonlinear chain with six-order polynomial on-site potential is used to analyze the evolution of the total to kinetic energy ratio during development of modulational instability of extended nonlinear vibrational modes. For the on-site potential of hard-type (soft-type) anharmonicity, the instability of $q=π$ mode ($q=0$ mode) results in the appearance of long-living discrete breathers (DBs) that gradually radiate their energy and eventually the system approaches thermal equilibrium with spatially uniform and temporally constant temperature. In the hard-type (soft-type) anharmonicity case, the total to kinetic energy ratio is minimal (maximal) in the regime of maximal energy localization by DBs. It is concluded that DBs affect specific heat of the nonlinear chain and for the case of hard-type (soft-type) anharmonicity they reduce (increase) the specific heat.

preprint2020arXiv

Effect of discrete breathers on macroscopic properties of the Fermi-Pasta-Ulam chain

The effect of discrete breathers (DBs) on macroscopic properties of the Fermi-Pasta-Ulam chain with symmetric and asymmetric potentials is investigated. The total to kinetic energy ratio (related to specific heat), stress (related to thermal expansion), and Young's modulus are monitored during the development of modulational instability of the zone boundary mode. The instability results in the formation of chaotic DBs followed by the transition to thermal equilibrium when DBs disappear due to energy radiation in the form of small-amplitude phonons. It is found that DBs reduce the specific heat for all the considered chain parameters. They increase the thermal expansion when the potential is asymmetric and, as expected, thermal expansion is not observed in the case of symmetric potential. The Young's modulus in the presence of DBs is smaller than in thermal equilibrium for the symmetric potential and for the potential with a small asymmetry, but it is larger than in thermal equilibrium for the potential with greater asymmetry. Our results can be useful for setting experiments on the identification of DBs in crystals by measuring their macroscopic properties.

preprint2020arXiv

Experimental observation of end-pinching in the Rayleigh breakup of an electrodynamically levitated charged droplet

The experimental time-lapse images of the breakup phenomenon of a charged droplet (diameter ~100-300 micro-m) levitated in an electrodynamic (ED) balance is reported. During the breakup process, a levitated charged droplet undergoes evaporation leading to a reduction in droplet size and increase in the corresponding surface charge density. As the surface charge density reaches to a critical value, known as the Rayleigh limit, the droplet undergoes breakup by forming a jet which further ejects highly charged progeny droplets. All the successive events of the droplet breakup process such as drop deformation, breakup, and relaxation of the drop back to spherical shape after ejection of progeny droplets have been recorded using the high-speed camera at 1.3 hundred thousand frames per second. The droplet is observed to eject 3-5 progeny droplets from a jet indicating end pinching mode of breakup. The jet is then observed to relax back after ejecting 31% of the total charge and ~3\% of the total mass. A suitable theory is provided to supplement the experimental observations, and a reasonable agreement is observed. Additionally, the theory is extended for the prediction of total mass loss and the entire lifetime of a charged droplet.

preprint2020arXiv

Maximizing Determinants under Matroid Constraints

Given vectors $v_1,\dots,v_n\in\mathbb{R}^d$ and a matroid $M=([n],I)$, we study the problem of finding a basis $S$ of $M$ such that $\det(\sum_{i \in S}v_i v_i^\top)$ is maximized. This problem appears in a diverse set of areas such as experimental design, fair allocation of goods, network design, and machine learning. The current best results include an $e^{2k}$-estimation for any matroid of rank $k$ and a $(1+ε)^d$-approximation for a uniform matroid of rank $k\ge d+\frac dε$, where the rank $k\ge d$ denotes the desired size of the optimal set. Our main result is a new approximation algorithm with an approximation guarantee that depends only on the dimension $d$ of the vectors and not on the size $k$ of the output set. In particular, we show an $(O(d))^{d}$-estimation and an $(O(d))^{d^3}$-approximation for any matroid, giving a significant improvement over prior work when $k\gg d$. Our result relies on the existence of an optimal solution to a convex programming relaxation for the problem which has sparse support; in particular, no more than $O(d^2)$ variables of the solution have fractional values. The sparsity results rely on the interplay between the first-order optimality conditions for the convex program and matroid theory. We believe that the techniques introduced to show sparsity of optimal solutions to convex programs will be of independent interest. We also give a randomized algorithm that rounds a sparse fractional solution to a feasible integral solution to the original problem. To show the approximation guarantee, we utilize recent works on strongly log-concave polynomials and show new relationships between different convex programs studied for the problem. Finally, we use the estimation algorithm and sparsity results to give an efficient deterministic approximation algorithm with an approximation guarantee that depends solely on the dimension $d$.

preprint2020arXiv

Multi-Criteria Dimensionality Reduction with Applications to Fairness

Dimensionality reduction is a classical technique widely used for data analysis. One foundational instantiation is Principal Component Analysis (PCA), which minimizes the average reconstruction error. In this paper, we introduce the "multi-criteria dimensionality reduction" problem where we are given multiple objectives that need to be optimized simultaneously. As an application, our model captures several fairness criteria for dimensionality reduction such as our novel Fair-PCA problem and the Nash Social Welfare (NSW) problem. In Fair-PCA, the input data is divided into $k$ groups, and the goal is to find a single $d$-dimensional representation for all groups for which the minimum variance of any one group is maximized. In NSW, the goal is to maximize the product of the individual variances of the groups achieved by the common low-dimensional space. Our main result is an exact polynomial-time algorithm for the two-criterion dimensionality reduction problem when the two criteria are increasing concave functions. As an application of this result, we obtain a polynomial time algorithm for Fair-PCA for $k=2$ groups and a polynomial time algorithm for NSW objective for $k=2$ groups. We also give approximation algorithms for $k>2$. Our technical contribution in the above results is to prove new low-rank properties of extreme point solutions to semi-definite programs. We conclude with experiments indicating the effectiveness of algorithms based on extreme point solutions of semi-definite programs on several real-world data sets.

preprint2020arXiv

On the Unreasonable Effectiveness of the Greedy Algorithm: Greedy Adapts to Sharpness

Submodular maximization has been widely studied over the past decades, mostly because of its numerous applications in real-world problems. It is well known that the standard greedy algorithm guarantees a worst-case approximation factor of 1-1/e when maximizing a monotone submodular function under a cardinality constraint. However, empirical studies show that its performance is substantially better in practice. This raises a natural question of explaining this improved performance of the greedy algorithm. In this work, we define sharpness for submodular functions as a candidate explanation for this phenomenon. The sharpness criterion is inspired by the concept of strong convexity in convex optimization. We show that the greedy algorithm provably performs better as the sharpness of the submodular function increases. This improvement ties in closely with the faster convergence rates of first order methods for sharp functions in convex optimization. Finally, we perform a computational study to empirically support our theoretical results and show that sharpness explains the greedy performance better than other justifications in the literature.

preprint2020arXiv

Rank one tensor completion problem

In this paper, we consider the rank-one tensor completion problem. We address the question of existence and uniqueness of the rank-one solution. In particular we show that the global uniqueness over the field of real numbers can be verified in a polynomial time. We give examples showing that there is an essential difference between the question of global uniqueness over the fields of real and complex numbers. Finally we briefly discuss the rank-one approximation problem for noisy observations.

preprint2019arXiv

Effect of trap potential on the Rayleigh breakup of a levitated charged droplet

Rayleigh instability that results in the breakup of a charged droplet, levitated in a quadrupole trap, has been investigated in the literature, but only scarcely. We report here asymmetric breakup of a charged drop, levitated in a loose trap, wherein, the droplet is stabilized at an off-center location in the trap. This aspect of levitation leads to an asymmetric breakup of the charged drop, predominantly in a direction opposite to that of gravity. In a first of its kind of study, we capture the successive events of the droplet deformation, breakup and relaxation of the drop after jet ejection using high speed imaging at a couple of hundred thousand frames per second. A pertinent question of the effect of the electrodynamic trap parameters such as applied voltage as well as physical parameters such as the size of the drop, gravity and conductivity on the characteristics of droplet breakup is also explored. A clear effect of the trap strength on the deformation (both symmetric and asymmetric) is observed. Moreover, the cone angle at the pole undergoing asymmetric breakup is almost independent of the applied field investigated in the experiments. All the experimental observations are compared with numerical simulations carried out using the boundary element method (BEM) in the Stokes flow limit. The BEM simulations are also extended to other experimentally achievable parameters. It is observed that the breakup is mostly field influenced, and not field induced. A plausible theory for the observations is reported, and a sensitive role of the sign of the charge on the droplet and the sign of the end cap potential, as well as the off-center location of the droplet in the trap.

preprint2019arXiv

Stability of electrodynamically levitated one or many charged droplets in the presence of noise

The theory of the effect of external fluctuation force on the stability and spatial distribution of mutually interacting and slowly evaporating charged drops, levitated in an electrodynamic balance, is presented using classical pseudo-potential approach. The theory is supplemented with numerical simulations where the non-homogeneous modified Mathieu equation is solved for single droplet as well as many droplets. The transition from the well ordered Coulombic crystal to randomly distributed liquid like structure is observed above a threshold value of the order parameter. The theory and simulations are found to be in fair agreement with each other. The simulation is aimed at studying the stability of structures for capturing the pollutant particle form the air streams using contactless membrane.