-
Duality of Hoffman constants
Authors:
Javier F. Pena,
Juan C. Vera,
Luis F. Zuluaga
Abstract:
Suppose $A\in \mathbb{R}^{m\times n}$ and consider the following canonical systems of inequalities defined by $A$: $$ \begin{array}{l} Ax=b\\ x \ge 0 \end{array} \qquad \text{ and }\qquad A^T y - c \le 0. $$ We establish some novel duality relationships between the Hoffman constants for the above constraint systems of linear inequalities provided some suitable Slater condition holds. The crux of o…
▽ More
Suppose $A\in \mathbb{R}^{m\times n}$ and consider the following canonical systems of inequalities defined by $A$: $$ \begin{array}{l} Ax=b\\ x \ge 0 \end{array} \qquad \text{ and }\qquad A^T y - c \le 0. $$ We establish some novel duality relationships between the Hoffman constants for the above constraint systems of linear inequalities provided some suitable Slater condition holds. The crux of our approach is a Hoffman duality inequality for polyhedral systems of constraints. The latter in turn yields an interesting duality identity between the Hoffman constants of the following box-constrained systems of inequalities: $$ \begin{array}{l} Ax=b\\ \ell \le x \le u \end{array}\qquad \text{ and }\qquad \ell \le A^T y - c \le u $$ for $\ell, u\in \mathbb{R}^n$ with $\ell < u.$
△ Less
Submitted 15 December, 2023;
originally announced December 2023.
-
Optimization hierarchies for distance-avoiding sets in compact spaces
Authors:
Bram Bekker,
Olga Kuryatnikova,
Fernando Mário de Oliveira Filho,
Juan C. Vera
Abstract:
Witsenhausen's problem asks for the maximum fraction $α_n$ of the $n$-dimensional unit sphere that can be covered by a measurable set containing no pairs of orthogonal points. The best upper bounds for $α_n$ are given by extensions of the Lovász theta number. In this paper, optimization hierarchies based on the Lovász theta number, like the Lasserre hierarchy, are extended to Witsenhausen's proble…
▽ More
Witsenhausen's problem asks for the maximum fraction $α_n$ of the $n$-dimensional unit sphere that can be covered by a measurable set containing no pairs of orthogonal points. The best upper bounds for $α_n$ are given by extensions of the Lovász theta number. In this paper, optimization hierarchies based on the Lovász theta number, like the Lasserre hierarchy, are extended to Witsenhausen's problem and similar problems. These hierarchies are shown to converge and are used to compute the best upper bounds for $α_n$ in low dimensions.
△ Less
Submitted 11 April, 2023;
originally announced April 2023.
-
Linear convergence of the Douglas-Rachford algorithm via a generic error bound condition
Authors:
Javier Peña,
Juan C. Vera,
Luis F. Zuluaga
Abstract:
We provide new insight into the convergence properties of the Douglas-Rachford algorithm for the problem $\min_x \{f(x)+g(x)\}$, where $f$ and $g$ are convex functions. Our approach relies on and highlights the natural primal-dual symmetry between the above problem and its Fenchel dual $\min_{u} \{ f^*(u) + g_*(u)\}$ where $g_*(u):=g^*(-u)$. Our main development is to show the linear convergence o…
▽ More
We provide new insight into the convergence properties of the Douglas-Rachford algorithm for the problem $\min_x \{f(x)+g(x)\}$, where $f$ and $g$ are convex functions. Our approach relies on and highlights the natural primal-dual symmetry between the above problem and its Fenchel dual $\min_{u} \{ f^*(u) + g_*(u)\}$ where $g_*(u):=g^*(-u)$. Our main development is to show the linear convergence of the algorithm when a natural error bound condition on the Douglas-Rachford operator holds. We leverage our error bound condition approach to show and estimate the algorithm's linear rate of convergence for three special classes of problems. The first one is when $f$ or$g$ and $f^*$ or $g_*$ are strongly convex relative to the primal and dual optimal sets respectively. The second one is when~$f$ and~$g$ are piecewise linear-quadratic functions. The third one is when~$f$ and~$g$ are the indicator functions of closed convex cones. In all three cases the rate of convergence is determined by a suitable measure of well-posedness of the problem. In the conic case, if the two closed convex cones are a linear subspace $L$ and $\mathbb{R}^n_+$, we establish the following stronger {\em finite termination} result: the Douglas-Rachford algorithm identifies the {\em maximum support sets} for $L\cap \mathbb{R}^n_+$ and $L^{\perp}\cap\mathbb{R}^n_+$ in finitely many steps. Our developments have straightforward extensions to the more general linearly constrained problem $\min_{x,y} \{f(x) + g( y):Ax + By = b\}$ thereby highlighting a direct and straightforward relationship between the Douglas-Rachford algorithm and the alternating direction method of multipliers (ADMM).
△ Less
Submitted 11 November, 2021;
originally announced November 2021.
-
Sparse non-SOS Putinar-type Positivstellensätze
Authors:
Lorenz M. Roebers,
Juan C. Vera,
Luis F. Zuluaga
Abstract:
Recently, non-SOS Positivstellensätze for polynomials on compact semialgebraic sets, following the general form of Schmüdgen's Positivstellensatz, have been derived by appropriately replacing the SOS polynomials with other classes of polynomials. An open question in the literature is how to obtain similar results following the general form of Putinar's Positivstellensatz. Extrapolating the algebra…
▽ More
Recently, non-SOS Positivstellensätze for polynomials on compact semialgebraic sets, following the general form of Schmüdgen's Positivstellensatz, have been derived by appropriately replacing the SOS polynomials with other classes of polynomials. An open question in the literature is how to obtain similar results following the general form of Putinar's Positivstellensatz. Extrapolating the algebraic geometry tools used to obtain this type of result in the SOS case fails to answer this question, because algebraic geometry methods strongly use hallmark properties of the class of SOS polynomials, such as closure under multiplication and closure under composition with other polynomials. In this article, using a new approach, we show the existence of Putinar-type Positivstellensätze that are constructed using non-SOS classes of non-negative polynomials, such as SONC, SDSOS and DSOS polynomials. Even not necessarily non-negative classes of polynomials such as sums of arithmetic-mean/geometric-mean polynomials could be used. Furthermore, we show that these certificates can be written with inherent sparsity characteristics. Such characteristics can be further exploited when the sparsity structure of both the polynomial whose non-negativity is being certified and the polynomials defining the semialgebraic set of interest are known. In contrast with related literature focused on exploiting sparsity in SOS Positivstellensätze, these latter results show how to exploit sparsity in a more general setting in which non-SOS polynomials are used to construct the Positivstellensätze.
△ Less
Submitted 12 October, 2021;
originally announced October 2021.
-
Computing near-optimal Value-at-Risk portfolios using Integer Programming techniques
Authors:
Onur Babat,
Juan C. Vera,
Luis F. Zuluaga
Abstract:
Value-at-Risk (VaR) is one of the main regulatory tools used for risk management purposes. However, it is difficult to compute optimal VaR portfolios; that is, an optimal risk-reward portfolio allocation using VaR as the risk measure. This is due to VaR being non-convex and of combinatorial nature. In particular, it is well known that the VaR portfolio problem can be formulated as a mixed integer…
▽ More
Value-at-Risk (VaR) is one of the main regulatory tools used for risk management purposes. However, it is difficult to compute optimal VaR portfolios; that is, an optimal risk-reward portfolio allocation using VaR as the risk measure. This is due to VaR being non-convex and of combinatorial nature. In particular, it is well known that the VaR portfolio problem can be formulated as a mixed integer linear program (MILP) that is difficult to solve with current MILP solvers for medium to large-scale instances of the problem. Here, we present an algorithm to compute near-optimal VaR portfolios that takes advantage of this MILP formulation and provides a guarantee of the solution's near-optimality. As a byproduct, we obtain an algorithm to compute tight lower bounds on the VaR portfolio problem that outperform related algorithms proposed in the literature for this purpose. The near-optimality guarantee provided by the proposed algorithm is obtained thanks to the relation between minimum risk portfolios satisfying a reward benchmark and the corresponding maximum reward portfolios satisfying a risk benchmark. These alternate formulations of the portfolio allocation problem have been frequently studied in the case of convex risk measures and concave reward functions. Here, this relationship is considered for general risk measures and reward functions. To illustrate the efficiency of the presented algorithm, numerical results are presented using historical asset returns from the US financial market.
△ Less
Submitted 22 June, 2021;
originally announced July 2021.
-
Investigation of the correlation patterns and the Compton dominance variability of Mrk 421 in 2017
Authors:
MAGIC Collaboration,
V. A. Acciari,
S. Ansoldi,
L. A. Antonelli,
A. Arbet Engels,
M. Artero,
K. Asano,
A. Babić,
A. Baquero,
U. Barres de Almeida,
J. A. Barrio,
I. Batković,
J. Becerra González,
W. Bednarek,
L. Bellizzi,
E. Bernardini,
M. Bernardos,
A. Berti,
J. Besenrieder,
W. Bhattacharyya,
C. Bigongiari,
O. Blanch,
Ž. Bošnjak,
G. Busetto,
R. Carosi
, et al. (263 additional authors not shown)
Abstract:
We present a detailed characterisation and theoretical interpretation of the broadband emission of the paradigmatic TeV blazar Mrk 421, with special focus on the multi-band flux correlations. The dataset has been collected through an extensive multiwavelength campaign organised between 2016 December and 2017 June. The instruments involved are MAGIC, FACT, Fermi-LAT, Swift, GASP-WEBT, OVRO, Medicin…
▽ More
We present a detailed characterisation and theoretical interpretation of the broadband emission of the paradigmatic TeV blazar Mrk 421, with special focus on the multi-band flux correlations. The dataset has been collected through an extensive multiwavelength campaign organised between 2016 December and 2017 June. The instruments involved are MAGIC, FACT, Fermi-LAT, Swift, GASP-WEBT, OVRO, Medicina and Metsähovi. Additionally, four deep exposures (several hours long) with simultaneous MAGIC and NuSTAR observations allowed a precise measurement of the falling segments of the two spectral components. The very-high-energy (VHE; E > 100 GeV) gamma rays and X-rays are positively correlated at zero time lag, but the strength and characteristics of the correlation change substantially across the various energy bands probed. The VHE versus X-ray fluxes follow different patterns, partly due to substantial changes in the Compton dominance during a few days without a simultaneous increase in the X-ray flux (i.e. orphan gamma-ray activity). Studying the broadband spectral energy distribution (SED) during the days including NuSTAR observations, we show that these changes can be explained within a one-zone leptonic model with a blob that increases its size over time. Our multi-band correlation study also hints at an anti-correlation between UV/optical and X-ray at a significance higher than 3 sigmas. A VHE flare observed on 2017 February 4 shows gamma-ray variability on multi-hour timescales, with a factor 10 increase in the TeV flux but only a moderate increase in the keV flux. The related broadband SED is better described by a two-zone leptonic scenario rather than by a one-zone scenario. We find that the flare can be produced by the appearance of a compact second blob populated by high energetic electrons spanning a narrow range of Lorentz factors.
△ Less
Submitted 10 June, 2021;
originally announced June 2021.
-
VHE gamma-ray detection of FSRQ QSO B1420+326 and modeling of its enhanced broadband state in 2020
Authors:
V. A. Acciari,
S. Ansoldi,
L. A. Antonelli,
A. Arbet Engels,
M. Artero,
K. Asano,
D. Baack,
A. Babić,
A. Baquero,
U. Barres de Almeida,
J. A. Barrio,
J. Becerra González,
W. Bednarek,
L. Bellizzi,
E. Bernardini,
M. Bernardos,
A. Berti,
J. Besenrieder,
W. Bhattacharyya,
C. Bigongiari,
A. Biland,
O. Blanch,
G. Bonnoli,
Ž. Bošnjak,
G. Busetto
, et al. (209 additional authors not shown)
Abstract:
Context. QSO B1420+326 is a blazar classified as a Flat Spectrum Radio Quasar (FSRQ). In the beginning of 2020 it underwent an enhanced flux state. An extensive multiwavelength campaign allowed us to trace the evolution of the flare. Aims. We search for VHE gamma-ray emission from QSO B1420+326 during this flaring state. We aim to characterize and model the broadband emission of the source over di…
▽ More
Context. QSO B1420+326 is a blazar classified as a Flat Spectrum Radio Quasar (FSRQ). In the beginning of 2020 it underwent an enhanced flux state. An extensive multiwavelength campaign allowed us to trace the evolution of the flare. Aims. We search for VHE gamma-ray emission from QSO B1420+326 during this flaring state. We aim to characterize and model the broadband emission of the source over different phases of the flare. Methods. The source was observed with a number of instruments in radio, near infrared, optical (including polarimetry and spectroscopy), ultra-violet, X-ray and gamma-ray bands. We use dedicated optical spectroscopy results to estimate the accretion disk and the dust torus luminosity. We perform spectral energy distribution modeling in the framework of combined Synchrotron-Self-Compton and External Compton scenario in which the electron energy distribution is partially determined from acceleration and cooling processes. Results. During the enhanced state the flux of both SED components drastically increased and the peaks were shifted to higher energies. Follow up observations with the MAGIC telescopes led to the detection of very-high-energy gamma-ray emission from this source, making it one of only a handful of FSRQs known in this energy range. Modeling allows us to constrain the evolution of the magnetic field and electron energy distribution in the emission region. The gamma-ray flare was accompanied by a rotation of the optical polarization vector during a low polarization state. Also, a new, superluminal radio knot contemporaneously appeared in the radio image of the jet. The optical spectroscopy shows a prominent FeII bump with flux evolving together with the continuum emission and a MgII line with varying equivalent width.
△ Less
Submitted 21 December, 2020;
originally announced December 2020.
-
Multiwavelength variability and correlation studies of Mrk 421 during historically low X-ray and $γ$-ray activity in 2015$-$2016
Authors:
MAGIC Collaboration,
V. A. Acciari,
S. Ansoldi,
L. A. Antonelli,
K. Asano,
A. Babić,
B. Banerjee,
A. Baquero,
U. Barres de Almeida,
J. A. Barrio,
J. Becerra González,
W. Bednarek,
L. Bellizzi,
E. Bernardini,
M. Bernardos,
A. Berti,
J. Besenrieder,
W. Bhattacharyya,
C. Bigongiari,
O. Blanch,
G. Bonnoli,
Ž. Bošnjak,
G. Busetto,
R. Carosi,
G. Ceribella
, et al. (205 additional authors not shown)
Abstract:
We report a characterization of the multi-band flux variability and correlations of the nearby (z=0.031) blazar Markarian 421 (Mrk 421) using data from Metsähovi, Swift, Fermi-LAT, MAGIC, FACT and other collaborations and instruments from November 2014 till June 2016. Mrk 421 did not show any prominent flaring activity, but exhibited periods of historically low activity above 1 TeV (F…
▽ More
We report a characterization of the multi-band flux variability and correlations of the nearby (z=0.031) blazar Markarian 421 (Mrk 421) using data from Metsähovi, Swift, Fermi-LAT, MAGIC, FACT and other collaborations and instruments from November 2014 till June 2016. Mrk 421 did not show any prominent flaring activity, but exhibited periods of historically low activity above 1 TeV (F$_{>1\mathrm{TeV}}<$ 1.7$\times$10$^{-12}$ ph cm$^{-2}$ s$^{-1}$) and in the 2-10 keV (X-ray) band (F$_{2-10 \mathrm{keV}}<$3.6$\times$10$^{-11}$ erg cm$^{-2}$ s$^{-1}$), during which the Swift-BAT data suggests an additional spectral component beyond the regular synchrotron emission. The highest flux variability occurs in X-rays and very-high-energy (E$>$0.1 TeV) $γ$-rays, which, despite the low activity, show a significant positive correlation with no time lag. The HR$_\mathrm{keV}$ and HR$_\mathrm{TeV}$ show the harder-when-brighter trend observed in many blazars, but the trend flattens at the highest fluxes, which suggests a change in the processes dominating the blazar variability. Enlarging our data set with data from years 2007 to 2014, we measured a positive correlation between the optical and the GeV emission over a range of about 60 days centered at time lag zero, and a positive correlation between the optical/GeV and the radio emission over a range of about 60 days centered at a time lag of $43^{+9}_{-6}$ days.This observation is consistent with the radio-bright zone being located about 0.2 parsec downstream from the optical/GeV emission regions of the jet. The flux distributions are better described with a LogNormal function in most of the energy bands probed, indicating that the variability in Mrk 421 is likely produced by a multiplicative process.
△ Less
Submitted 2 December, 2020;
originally announced December 2020.
-
Static Hedging of Weather and Price Risks in Electricity Markets
Authors:
Javier Pantoja Robayo,
Juan C. Vera
Abstract:
We present the closed-form solution to the problem of hedging price and quantity risks for energy retailers (ER), using financial instruments based on electricity price and weather indexes. Our model considers an ER who is intermediary in a regulated electricity market. ERs buy a fixed quantity of electricity at a variable cost and must serve a variable demand at a fixed cost. Thus ERs are subject…
▽ More
We present the closed-form solution to the problem of hedging price and quantity risks for energy retailers (ER), using financial instruments based on electricity price and weather indexes. Our model considers an ER who is intermediary in a regulated electricity market. ERs buy a fixed quantity of electricity at a variable cost and must serve a variable demand at a fixed cost. Thus ERs are subject to both price and quantity risks. To hedge such risks, an ER could construct a portfolio of financial instruments based on price and weather indexes. We construct the closed form solution for the optimal portfolio for the mean-Var model in the discrete setting. Our model does not make any distributional assumption.
△ Less
Submitted 17 November, 2020;
originally announced November 2020.
-
Contingent Capital with Stock Price Triggers in Interbank Networks
Authors:
Anne G. Balter,
Nikolaus Schweizer,
Juan C. Vera
Abstract:
This paper studies existence and uniqueness of equilibrium prices in a model of the banking sector in which banks trade contingent convertible bonds with stock price triggers among each other. This type of financial product was proposed as an instrument for stabilizing the global banking system after the financial crisis. Yet it was recognized early on that these products may create circularity pr…
▽ More
This paper studies existence and uniqueness of equilibrium prices in a model of the banking sector in which banks trade contingent convertible bonds with stock price triggers among each other. This type of financial product was proposed as an instrument for stabilizing the global banking system after the financial crisis. Yet it was recognized early on that these products may create circularity problems in the definition of stock prices - even in the absence of trade. We find that if conversion thresholds are such that bond holders are indifferent about marginal conversions, there exists a unique equilibrium irrespective of the network structure. When thresholds are lower, existence of equilibrium breaks down while higher thresholds may lead to multiplicity of equilibria. Moreover, there are complex network effects. One bank's conversion may trigger further conversions - or prevent them, depending on the constellations of asset values and conversion triggers.
△ Less
Submitted 12 November, 2020;
originally announced November 2020.
-
Long term variability and correlation study of the blazar 3C 454.3 in radio, NIR and optical wavebands
Authors:
A. Sarkar,
V. R. Chitnis,
A. C. Gupta,
H. Gaur,
S. R. Patel,
P. J. Wiita,
A. E. Volvach,
M. Tornikoski,
W. Chamani,
S. Enestam,
A. Lähteenmäki,
J. Tammi,
R. J. C Vera,
L. N. Volvach
Abstract:
We performed a long-term optical (B, V, R bands), infra-red (J and K bands) and radio band (15, 22, 37 GHz band) study on the flat spectrum radio quasar, 3C 454.3, using the data collected over a period of more than 8 years (MJD 54500--57500). The temporal variability, spectral properties and inter-waveband correlations were studied by dividing the available data into smaller segments with more re…
▽ More
We performed a long-term optical (B, V, R bands), infra-red (J and K bands) and radio band (15, 22, 37 GHz band) study on the flat spectrum radio quasar, 3C 454.3, using the data collected over a period of more than 8 years (MJD 54500--57500). The temporal variability, spectral properties and inter-waveband correlations were studied by dividing the available data into smaller segments with more regular sampling. This helped us constrain the size and the relative locations of the emission regions for different wavebands. Spectral analysis of the source revealed the interplay between the accretion disk and jet emission. The source predominantly showed a redder-when-brighter trend, though we observed a bluer-when-brighter trend at high flux levels which could be signatures of particle acceleration and radiative cooling. Significant correlations with near-zero lag were seen between various optical/infra-red bands, indicating that these emission regions are co-spatial. Correlations with a time lag of about 10--100 days are seen between optical/infra-red and radio bands indicating these emissions arise from different regions. We also observe the DCF peak lag change from year to year. We try to explain these differences using a curved jet model where the different emission regions have different viewing angles resulting in a frequency dependent Doppler factor. This variable Doppler factor model explains the variability timescales and the variation in DCF peak lag between the radio and optical emissions in different segments. Lags of 6-180 days are seen between emissions in various radio bands, indicating a core-shift effect
△ Less
Submitted 31 October, 2019;
originally announced October 2019.
-
Reducing non-negativity over general semialgebraic sets to non-negativity over simple sets
Authors:
Olga Kuryatnikova,
Juan C. Vera,
Luis F. Zuluaga
Abstract:
A non-negativity certificate (NNC) is a way to write a polynomial so that its non-negativity on a semialgebraic set becomes evident. Positivstellensätze (Psätze) guarantee the existence of NNCs. Both, NNCs and Psätze underlie powerful algorithmic techniques for optimization. This paper proposes a universal approach to derive new Psätze for general semialgebraic sets from ones developed for simpler…
▽ More
A non-negativity certificate (NNC) is a way to write a polynomial so that its non-negativity on a semialgebraic set becomes evident. Positivstellensätze (Psätze) guarantee the existence of NNCs. Both, NNCs and Psätze underlie powerful algorithmic techniques for optimization. This paper proposes a universal approach to derive new Psätze for general semialgebraic sets from ones developed for simpler sets, such as a box, a simplex, or the non-negative orthant. We provide several results illustrating the approach. First, by considering Handelman's Positivstellensatz (Psatz) over a box, we construct non-SOS Schmüdgen-type Psätze over any compact semialgebraic set. That is, a family of Psätze that follow the structure of the fundamental Schmüdgen's Psatz, but where instead of SOS polynomials, any class of polynomials containing the non-negative constants can be used, such as SONC, DSOS/SDSOS, hyperbolic or sums of AM/GM polynomials. Secondly, by considering the simplex as the simple set, we derive a sparse Psatz over general compact sets, which does not require any structural assumptions of the set. Finally, by considering Pólya's Psatz over the non-negative orthant, we derive a new non-SOS Psatz over unbounded sets which satisfy some generic conditions. All these results contribute to the literature regarding the use of non-SOS polynomials and sparse NNCs to derive Psätze over compact and unbounded sets. Throughout the article, we illustrate our results with relevant examples and numerical experiments.
△ Less
Submitted 22 March, 2024; v1 submitted 14 September, 2019;
originally announced September 2019.
-
Equivalence and invariance of the chi and Hoffman constants of a matrix
Authors:
Javier F. Pena,
Juan C. Vera,
Luis F. Zuluaga
Abstract:
We show that the following two condition measures of a full column rank matrix $A \in \mathbb{R}^{m\times n}$ are identical: the chi constant and a signed Hoffman constant. This identity is naturally suggested by the evident invariance of the chi constant under sign changes of the rows of $A$. We also show that similar equivalence and invariance properties extend to variants of the chi and Hoffman…
▽ More
We show that the following two condition measures of a full column rank matrix $A \in \mathbb{R}^{m\times n}$ are identical: the chi constant and a signed Hoffman constant. This identity is naturally suggested by the evident invariance of the chi constant under sign changes of the rows of $A$. We also show that similar equivalence and invariance properties extend to variants of the chi and Hoffman constants that depend only on the linear subspace $A(\mathbb{R}^n):=\{Ax: x\in\mathbb{R}^n\} \subseteq \mathbb{R}^m$. Finally, we show similar identities between the chi constants and signed versions of Renegar's and Grassmannian condition measures.
△ Less
Submitted 18 May, 2020; v1 submitted 15 May, 2019;
originally announced May 2019.
-
Generalizations of Schoenberg's theorem on positive definite kernels
Authors:
Olga Kuryatnikova,
Juan C. Vera
Abstract:
The seminal theorem of I.J. Schoenberg characterizes positive definite (p.d.) kernels on the unit sphere $S^{n-1}$ invariant under the automorphisms of the sphere. We obtain two generalizations of this theorem for p.d. kernels on fiber bundles. Our first theorem characterizes invariant p.d. kernels on bundles whose fiber is a product of a compact set and the unit sphere. This result implies, in pa…
▽ More
The seminal theorem of I.J. Schoenberg characterizes positive definite (p.d.) kernels on the unit sphere $S^{n-1}$ invariant under the automorphisms of the sphere. We obtain two generalizations of this theorem for p.d. kernels on fiber bundles. Our first theorem characterizes invariant p.d. kernels on bundles whose fiber is a product of a compact set and the unit sphere. This result implies, in particular, a characterization of invariant under the automorphisms of the sphere, p.d. kernels on a product of $S^{n-1}$ and a compact set. Our second result characterizes invariant p.d. kernels on the bundle whose fiber is $S^{n-1}$, base space is $(S^{n-1})^{r}$ and map is the projection on the base space. This set of kernels is isomorphic to the set of invariant under the automorphisms of the sphere continuous functions $F$ on $(S^{n-1})^{r+2}$ such that $F(\cdot,\cdot ,Z)$ is positive definite for every $Z \in (S^{n-1})^{r}$. When $Z$ is fixed, this class reduces to the class of p.d. kernels invariant under the stabilizer of $Z$ in the automorphism group of the sphere. For $r=1$ these kernels have been used to obtain upper bounds for the spherical codes problem. Our extension for $r>1$ can be used to construct new upper bounds on the size of spherical codes.
△ Less
Submitted 2 April, 2019;
originally announced April 2019.
-
Positive semidefinite approximations to the cone of copositive kernels
Authors:
Olga Kuryatnikova,
Juan C. Vera
Abstract:
It has been shown that the maximum stable set problem in some infinite graphs, and the kissing number problem in particular, reduces to a minimization problem over the cone of copositive kernels. Optimizing over this infinite dimensional cone is not tractable, and approximations of this cone have been hardly considered in literature. We propose two convergent hierarchies of subsets of copositive k…
▽ More
It has been shown that the maximum stable set problem in some infinite graphs, and the kissing number problem in particular, reduces to a minimization problem over the cone of copositive kernels. Optimizing over this infinite dimensional cone is not tractable, and approximations of this cone have been hardly considered in literature. We propose two convergent hierarchies of subsets of copositive kernels, in terms of non-negative and positive definite kernels. We use these hierarchies and representation theorems for invariant positive definite kernels on the sphere to construct new SDP-based bounds on the kissing number. This results in fast-to-compute upper bounds on the kissing number that lie between the currently existing LP and SDP bounds.
△ Less
Submitted 1 December, 2018;
originally announced December 2018.
-
Using Column Generation to Solve Extensions to the Markowitz Model
Authors:
Lorenz M. Roebers,
Aras Selvi,
Juan C. Vera
Abstract:
We introduce a solution scheme for portfolio optimization problems with cardinality constraints. Typical portfolio optimization problems are extensions of the classical Markowitz mean-variance portfolio optimization model. We solve such type of problems using a method similar to column generation. In this scheme, the original problem is restricted to a subset of the assets resulting in a master co…
▽ More
We introduce a solution scheme for portfolio optimization problems with cardinality constraints. Typical portfolio optimization problems are extensions of the classical Markowitz mean-variance portfolio optimization model. We solve such type of problems using a method similar to column generation. In this scheme, the original problem is restricted to a subset of the assets resulting in a master convex quadratic problem. Then the dual information of the master problem is used in a sub-problem to propose more assets to consider. We also consider other extensions to the Markowitz model to diversify the portfolio selection within the given intervals for active weights.
△ Less
Submitted 21 June, 2019; v1 submitted 30 November, 2018;
originally announced December 2018.
-
X-ray, UV, and Radio Timing Observations of the Radio Galaxy 3C 120
Authors:
Alan P. Marscher,
Svetlana G. Jorstad,
Karen E. Williamson,
Anne Lähteenmäki,
Merja Tornikoski,
John M. Hunter,
Katya A. Leidig,
Muhammad Zain Mobeen,
Rafael J. C. Vera,
Wara Chamani
Abstract:
We report the results of monitoring of the radio galaxy 3C 120 with the Neil Gehrels Swift Observatory, Very Long Baseline Array, and Metsähovi Radio Observatory. The UV-optical continuum spectrum and R-band polarization can be explained by a superposition of an inverted-spectrum source with a synchrotron component containing a disordered magnetic field. The UV-optical and X-ray light curves inclu…
▽ More
We report the results of monitoring of the radio galaxy 3C 120 with the Neil Gehrels Swift Observatory, Very Long Baseline Array, and Metsähovi Radio Observatory. The UV-optical continuum spectrum and R-band polarization can be explained by a superposition of an inverted-spectrum source with a synchrotron component containing a disordered magnetic field. The UV-optical and X-ray light curves include dips and flares, while several superluminal knots appear in the parsec-scale jet. The recovery time of the second dip was longer at UV-optical wavelengths, in conflict with a model in which the inner accretion disk (AD) is disrupted during a dip and then refilled from outer to inner radii. We favor an alternative scenario in which occasional polar alignments of the magnetic field in the disk and corona cause the flux dips and formation of shocks in the jet. Similar to observations of Seyfert galaxies, intra-band time lags of flux variations are longer than predicted by the standard AD model. This suggests that scattering or some other reprocessing occurs. The 37 GHz light curve is well correlated with the optical-UV variations, with a ~20-day delay. A radio flare in the jet occurred in a superluminal knot 0.14 milliarcseconds downstream of the 43 GHz "core," which places the site of the preceding X-ray/UV/optical flare within the core 0.5-1.3 pc from the black hole. The inverted UV-optical flare spectrum can be explained by a nearly mono-energetic electron distribution with energy similar to the minimum energy inferred in the TeV gamma-ray emitting regions of some BL Lacertae objects.
△ Less
Submitted 25 September, 2018;
originally announced September 2018.
-
Deciding Robust Feasibility and Infeasibility Using a Set Containment Approach: An Application to Stationary Passive Gas Network Operations
Authors:
Denis Aßmann,
Frauke Liers,
Michael Stingl,
Juan C. Vera
Abstract:
In this paper we study feasibility and infeasibility of nonlinear two-stage fully adjustable robust feasibility problems with an empty first stage. This is equivalent to deciding whether the uncertainty set is contained within the projection of the feasible region onto the uncertainty-space. Moreover, the considered sets are assumed to be described by polynomials. For answering this question, two…
▽ More
In this paper we study feasibility and infeasibility of nonlinear two-stage fully adjustable robust feasibility problems with an empty first stage. This is equivalent to deciding whether the uncertainty set is contained within the projection of the feasible region onto the uncertainty-space. Moreover, the considered sets are assumed to be described by polynomials. For answering this question, two very general approaches using methods from polynomial optimization are presented - one for showing feasibility and one for showing infeasibility. The developed methods are approximated through sum of squares polynomials and solved using semidefinite programs. Deciding robust feasibility and infeasibility is important for gas network operations, which is a nonconvex feasibility problem where the feasible set is described by a composition of polynomials with the absolute value function. Concerning the gas network problem, different topologies are considered. It is shown that a tree structured network can be decided exactly using linear programming. Furthermore, a method is presented to reduce a tree network with one additional arc to a single cycle network. In this case, the problem can be decided by eliminating the absolute value functions and solving the resulting linearly many polynomial optimization problems. Lastly, the effectivity of the methods is tested on a variety of small cyclic networks. It turns out that for instances where robust feasibility or infeasibility can be decided successfully, level 2 or level 3 of the Lasserre relaxation hierarchy typically is sufficient.
△ Less
Submitted 30 August, 2018;
originally announced August 2018.
-
The broad-band properties of the intermediate synchrotron peaked BL Lac S2 0109+22 from radio to VHE gamma rays
Authors:
MAGIC Collaboration,
S. Ansoldi,
L. A. Antonelli,
C. Arcaro,
D. Baack,
A. Babić,
B. Banerjee,
P. Bangale,
U. Barres de Almeida,
J. A. Barrio,
J. Becerra González,
W. Bednarek,
E. Bernardini,
R. Ch. Berse,
A. Berti,
J. Besenrieder,
W. Bhattacharyya,
C. Bigongiari,
A. Biland,
O. Blanch,
G. Bonnoli,
R. Carosi,
G. Ceribella,
A. Chatterjee,
S. M. Colak
, et al. (141 additional authors not shown)
Abstract:
The MAGIC telescopes observed S2 0109+22 in 2015 July during its flaring activity in high energy gamma rays observed by Fermi-LAT. We analyse the MAGIC data to characterise the very high energy (VHE) gamma-ray emission of S2 0109+22, which belongs to the subclass of intermediate synchrotron peak (ISP) BL Lac objects. We study the multi-frequency emission in order to investigate the source classifi…
▽ More
The MAGIC telescopes observed S2 0109+22 in 2015 July during its flaring activity in high energy gamma rays observed by Fermi-LAT. We analyse the MAGIC data to characterise the very high energy (VHE) gamma-ray emission of S2 0109+22, which belongs to the subclass of intermediate synchrotron peak (ISP) BL Lac objects. We study the multi-frequency emission in order to investigate the source classification. Finally, we compare the source long-term behaviour to other VHE gamma-ray emitting (TeV) blazars. We performed a temporal and spectral analysis of the data centred around the MAGIC interval of observation (MJD 57225-57231). Long-term radio and optical data have also been investigated using the discrete correlation function. The redshift of the source is estimated through optical host-galaxy imaging and also using the amount of VHE gamma-ray absorption. The quasi-simultaneous multi-frequency spectral energy distribution (SED) is modelled with the conventional one-zone synchrotron self-Compton (SSC) model. MAGIC observations resulted in the detection of the source at a significance level of $5.3\,σ$. The VHE gamma-ray emission of S2 0109+22 is variable on a daily time scale. VHE gamma-ray luminosity of the source is lower than the average of TeV BL Lacs. The optical polarization, and long-term optical/radio behaviour of the source are different from the general population of TeV blazars. All these findings agree with the classification of the source as an ISP BL Lac object. We estimate the source redshift as $z = 0.36 \pm 0.07$. The SSC parameters describing the SED are rather typical for blazars.
△ Less
Submitted 5 July, 2018;
originally announced July 2018.
-
Radio jets and gamma-ray emission in radio-silent narrow-line Seyfert 1 galaxies
Authors:
A. Lähteenmäki,
E. Järvelä,
V. Ramakrishnan,
M. Tornikoski,
J. Tammi,
R. J. C. Vera,
W. Chamani
Abstract:
We have detected six narrow-line Seyfert 1 (NLS1) galaxies at 37 GHz that were previously classified as radio silent and two that were classified as radio quiet. These detections reveal the presumption that NLS1 galaxies labelled radio quiet or radio silent and hosted by spiral galaxies are unable to launch jets to be incorrect. The detections are a plausible indicator of the presence of a powerfu…
▽ More
We have detected six narrow-line Seyfert 1 (NLS1) galaxies at 37 GHz that were previously classified as radio silent and two that were classified as radio quiet. These detections reveal the presumption that NLS1 galaxies labelled radio quiet or radio silent and hosted by spiral galaxies are unable to launch jets to be incorrect. The detections are a plausible indicator of the presence of a powerful, most likely relativistic jet because this intensity of emission at 37 GHz cannot be explained by, for example, radiation from supernova remnants. Additionally, one of the detected NLS1 galaxies is a newly discovered source of gamma rays and three others are candidates for future detections.
△ Less
Submitted 6 June, 2018;
originally announced June 2018.
-
New bounds for truthful scheduling on two unrelated selfish machines
Authors:
Olga Kuryatnikova,
Juan C. Vera
Abstract:
We consider the minimum makespan problem for $n$ tasks and two unrelated parallel selfish machines. Let $R_n$ be the best approximation ratio of randomized monotone scale-free algorithms. This class contains the most efficient algorithms known for truthful scheduling on two machines. We propose a new $Min-Max$ formulation for $R_n$, as well as upper and lower bounds on $R_n$ based on this formulat…
▽ More
We consider the minimum makespan problem for $n$ tasks and two unrelated parallel selfish machines. Let $R_n$ be the best approximation ratio of randomized monotone scale-free algorithms. This class contains the most efficient algorithms known for truthful scheduling on two machines. We propose a new $Min-Max$ formulation for $R_n$, as well as upper and lower bounds on $R_n$ based on this formulation. For the lower bound, we exploit pointwise approximations of cumulative distribution functions (CDFs). For the upper bound, we construct randomized algorithms using distributions with piecewise rational CDFs. Our method improves upon the existing bounds on $R_n$ for small $n$. In particular, we obtain almost tight bounds for $n=2$ showing that $|R_2-1.505996|<10^{-6}$.
△ Less
Submitted 25 September, 2019; v1 submitted 13 November, 2017;
originally announced November 2017.
-
Positive polynomials on unbounded domains
Authors:
Javer Pena,
Juan C. Vera,
Luis F. Zuluaga
Abstract:
Certificates of non-negativity such as Putinar's Positivstellensatz have been used to obtain powerful numerical techniques to solve polynomial optimization (PO) problems. Putinar's certificate uses sum-of-squares (sos) polynomials to certify the non-negativity of a given polynomial over a domain defined by polynomial inequalities. This certificate assumes the Archimedean property of the associated…
▽ More
Certificates of non-negativity such as Putinar's Positivstellensatz have been used to obtain powerful numerical techniques to solve polynomial optimization (PO) problems. Putinar's certificate uses sum-of-squares (sos) polynomials to certify the non-negativity of a given polynomial over a domain defined by polynomial inequalities. This certificate assumes the Archimedean property of the associated quadratic module, which in particular implies compactness of the domain. In this paper we characterize the existence of a certificate of non-negativity for polynomials over a possibly unbounded domain, without the use of the associated quadratic module. Next, we show that the certificate can be used to convergent linear matrix inequality (LMI) hierarchies for PO problems with unbounded feasible sets. Furthermore, by using copositive polynomials to certify non-negativity, instead of sos polynomials, the certificate allows the use of a very rich class of convergent LMI hierarchies to approximate the solution of general PO problems. Throughout the article we illustrate our results with various examples certifying the non-negativity of polynomials over possibly unbounded sets defined by polynomial equalities or inequalities.
△ Less
Submitted 11 September, 2017;
originally announced September 2017.
-
On the convergence rate of grid search for polynomial optimization over the simplex
Authors:
Etienne de Klerk,
Monique Laurent,
Zhao Sun,
Juan C. Vera
Abstract:
We consider the approximate minimization of a given polynomial on the standard simplex, obtained by taking the minimum value over all rational grid points with given denominator ${r} \in \mathbb{N}$. It was shown in [De Klerk, E., Laurent, M., Sun, Z.: An error analysis for polynomial optimization over the simplex based on the multivariate hypergeometric distribution. {\em SIAM J. Optim.} 25(3) 14…
▽ More
We consider the approximate minimization of a given polynomial on the standard simplex, obtained by taking the minimum value over all rational grid points with given denominator ${r} \in \mathbb{N}$. It was shown in [De Klerk, E., Laurent, M., Sun, Z.: An error analysis for polynomial optimization over the simplex based on the multivariate hypergeometric distribution. {\em SIAM J. Optim.} 25(3) 1498--1514 (2015)] that the relative accuracy of this approximation depends on $r$ as $O(1/r^2)$ if there exists a rational global minimizer. In this note we show that the rational minimizer condition is not necessary to obtain the $O(1/r^2)$ bound.
△ Less
Submitted 9 March, 2016; v1 submitted 5 October, 2015;
originally announced October 2015.
-
Improved Bounds on the Phase Transition for the Hard-Core Model in 2-Dimensions
Authors:
Juan C. Vera,
Eric Vigoda,
Linji Yang
Abstract:
For the hard-core lattice gas model defined on independent sets weighted by an activity $λ$, we study the critical activity $λ_c(\mathbb{Z}^2)$ for the uniqueness/non-uniqueness threshold on the 2-dimensional integer lattice $\mathbb{Z}^2$. The conjectured value of the critical activity is approximately $3.796$. Until recently, the best lower bound followed from algorithmic results of Weitz (2006)…
▽ More
For the hard-core lattice gas model defined on independent sets weighted by an activity $λ$, we study the critical activity $λ_c(\mathbb{Z}^2)$ for the uniqueness/non-uniqueness threshold on the 2-dimensional integer lattice $\mathbb{Z}^2$. The conjectured value of the critical activity is approximately $3.796$. Until recently, the best lower bound followed from algorithmic results of Weitz (2006). Weitz presented an FPTAS for approximating the partition function for graphs of constant maximum degree $Δ$ when $λ<λ_c(\mathbb{T}_Δ)$ where $\mathbb{T}_Δ$ is the infinite, regular tree of degree $Δ$. His result established a certain decay of correlations property called strong spatial mixing (SSM) on $\mathbb{Z}^2$ by proving that SSM holds on its self-avoiding walk tree $T_{\mathrm{saw}}^σ(\mathbb{Z}^2)$ where $σ=(σ_v)_{v\in \mathbb{Z}^2}$ and $σ_v$ is an ordering on the neighbors of vertex $v$. As a consequence he obtained that $λ_c(\mathbb{Z}^2)\geqλ_c( \mathbb{T}_4) = 1.675$. Restrepo et al. (2011) improved Weitz's approach for the particular case of $\mathbb{Z}^2$ and obtained that $λ_c(\mathbb{Z}^2)>2.388$. In this paper, we establish an upper bound for this approach, by showing that, for all $σ$, SSM does not hold on $T_{\mathrm{saw}}^σ(\mathbb{Z}^2)$ when $λ>3.4$. We also present a refinement of the approach of Restrepo et al. which improves the lower bound to $λ_c(\mathbb{Z}^2)>2.48$.
△ Less
Submitted 9 July, 2014; v1 submitted 3 June, 2013;
originally announced June 2013.
-
Phase Transition for Glauber Dynamics for Independent Sets on Regular Trees
Authors:
Ricardo Restrepo,
Daniel Stefankovic,
Juan C. Vera,
Eric Vigoda,
Linji Yang
Abstract:
We study the effect of boundary conditions on the relaxation time of the Glauber dynamics for the hard-core model on the tree. The hard-core model is defined on the set of independent sets weighted by a parameter $λ$, called the activity. The Glauber dynamics is the Markov chain that updates a randomly chosen vertex in each step. On the infinite tree with branching factor $b$, the hard-core model…
▽ More
We study the effect of boundary conditions on the relaxation time of the Glauber dynamics for the hard-core model on the tree. The hard-core model is defined on the set of independent sets weighted by a parameter $λ$, called the activity. The Glauber dynamics is the Markov chain that updates a randomly chosen vertex in each step. On the infinite tree with branching factor $b$, the hard-core model can be equivalently defined as a broadcasting process with a parameter $ω$ which is the positive solution to $λ=ω(1+ω)^b$, and vertices are occupied with probability $ω/(1+ω)$ when their parent is unoccupied. This broadcasting process undergoes a phase transition between the so-called reconstruction and non-reconstruction regions at $ω_r\approx \ln{b}/b$. Reconstruction has been of considerable interest recently since it appears to be intimately connected to the efficiency of local algorithms on locally tree-like graphs, such as sparse random graphs. In this paper we show that the relaxation time of the Glauber dynamics on regular $b$-ary trees $T_h$ of height $h$ and $n$ vertices, undergoes a phase transition around the reconstruction threshold. In particular, we construct a boundary condition for which the relaxation time slows down at the reconstruction threshold. More precisely, for any $ω\le \ln{b}/b$, for $T_h$ with any boundary condition, the relaxation time is $Ω(n)$ and $O(n^{1+o_b(1)})$. In contrast, above the reconstruction threshold we show that for every $δ>0$, for $ω=(1+δ)\ln{b}/b$, the relaxation time on $T_h$ with any boundary condition is $O(n^{1+δ+ o_b(1)})$, and we construct a boundary condition where the relaxation time is $Ω(n^{1+δ/2 - o_b(1)})$.
△ Less
Submitted 14 July, 2010;
originally announced July 2010.
-
Phase transition for the mixing time of the Glauber dynamics for coloring regular trees
Authors:
Prasad Tetali,
Juan C. Vera,
Eric Vigoda,
Linji Yang
Abstract:
We prove that the mixing time of the Glauber dynamics for random k-colorings of the complete tree with branching factor b undergoes a phase transition at $k=b(1+o_b(1))/\ln{b}$. Our main result shows nearly sharp bounds on the mixing time of the dynamics on the complete tree with n vertices for $k=Cb/\ln{b}$ colors with constant C. For $C\geq1$ we prove the mixing time is $O(n^{1+o_b(1)}\ln{n})$.…
▽ More
We prove that the mixing time of the Glauber dynamics for random k-colorings of the complete tree with branching factor b undergoes a phase transition at $k=b(1+o_b(1))/\ln{b}$. Our main result shows nearly sharp bounds on the mixing time of the dynamics on the complete tree with n vertices for $k=Cb/\ln{b}$ colors with constant C. For $C\geq1$ we prove the mixing time is $O(n^{1+o_b(1)}\ln{n})$. On the other side, for $C<1$ the mixing time experiences a slowing down; in particular, we prove it is $O(n^{1/C+o_b(1)}\ln{n})$ and $Ω(n^{1/C-o_b(1)})$. The critical point C=1 is interesting since it coincides (at least up to first order) with the so-called reconstruction threshold which was recently established by Sly. The reconstruction threshold has been of considerable interest recently since it appears to have close connections to the efficiency of certain local algorithms, and this work was inspired by our attempt to understand these connections in this particular setting.
△ Less
Submitted 27 November, 2012; v1 submitted 18 August, 2009;
originally announced August 2009.
-
Randomly coloring planar graphs with fewer colors than the maximum degree
Authors:
Thomas P. Hayes,
Juan C. Vera,
Eric Vigoda
Abstract:
We study Markov chains for randomly sampling $k$-colorings of a graph with maximum degree $Δ$. Our main result is a polynomial upper bound on the mixing time of the single-site update chain known as the Glauber dynamics for planar graphs when $k=Ω(Δ/\logΔ)$. Our results can be partially extended to the more general case where the maximum eigenvalue of the adjacency matrix of the graph is at most…
▽ More
We study Markov chains for randomly sampling $k$-colorings of a graph with maximum degree $Δ$. Our main result is a polynomial upper bound on the mixing time of the single-site update chain known as the Glauber dynamics for planar graphs when $k=Ω(Δ/\logΔ)$. Our results can be partially extended to the more general case where the maximum eigenvalue of the adjacency matrix of the graph is at most $Δ^{1-\eps}$, for fixed $\eps > 0$.
The main challenge when $k \le Δ+ 1$ is the possibility of "frozen" vertices, that is, vertices for which only one color is possible, conditioned on the colors of its neighbors. Indeed, when $Δ= O(1)$, even a typical coloring can have a constant fraction of the vertices frozen. Our proofs rely on recent advances in techniques for bounding mixing time using "local uniformity" properties.
△ Less
Submitted 31 August, 2011; v1 submitted 11 June, 2007;
originally announced June 2007.
