-
Communication-efficient Vertical Federated Learning via Compressed Error Feedback
Authors:
Pedro Valdeira,
João Xavier,
Cláudia Soares,
Yuejie Chi
Abstract:
Communication overhead is a known bottleneck in federated learning (FL). To address this, lossy compression is commonly used on the information communicated between the server and clients during training. In horizontal FL, where each client holds a subset of the samples, such communication-compressed training methods have recently seen significant progress. However, in their vertical FL counterpar…
▽ More
Communication overhead is a known bottleneck in federated learning (FL). To address this, lossy compression is commonly used on the information communicated between the server and clients during training. In horizontal FL, where each client holds a subset of the samples, such communication-compressed training methods have recently seen significant progress. However, in their vertical FL counterparts, where each client holds a subset of the features, our understanding remains limited. To address this, we propose an error feedback compressed vertical federated learning (EFVFL) method to train split neural networks. In contrast with previous communication-compressed methods for vertical FL, EFVFL does not require a vanishing compression error for the gradient norm to converge to zero for smooth nonconvex problems. By leveraging error feedback, our method can achieve a $\mathcal{O}(1/T)$ convergence rate in the full-batch case, improving over the state-of-the-art $\mathcal{O}(1/\sqrt{T})$ rate under $\mathcal{O}(1/\sqrt{T})$ compression error, and matching the rate of uncompressed methods. Further, when the objective function satisfies the Polyak-Łojasiewicz inequality, our method converges linearly. In addition to improving convergence rates, our method also supports the use of private labels. Numerical experiments show that EFVFL significantly improves over the prior art, confirming our theoretical results.
△ Less
Submitted 20 June, 2024;
originally announced June 2024.
-
Outlier-resilient model fitting via percentile losses: Methods for general and convex residuals
Authors:
João Domingos,
João Xavier
Abstract:
We consider the problem of robustly fitting a model to data that includes outliers by formulating a percentile optimization problem. This problem is non-smooth and non-convex, hence hard to solve. We derive properties that the minimizers of such problems must satisfy. These properties lead to methods that solve the percentile formulation both for general residuals and for convex residuals. The met…
▽ More
We consider the problem of robustly fitting a model to data that includes outliers by formulating a percentile optimization problem. This problem is non-smooth and non-convex, hence hard to solve. We derive properties that the minimizers of such problems must satisfy. These properties lead to methods that solve the percentile formulation both for general residuals and for convex residuals. The methods fit the model to subsets of the data, and then extract the solution of the percentile formulation from these partial fits. As illustrative simulations show, such methods endure higher outlier percentages, when compared with standard robust estimates. Additionally, the derived properties provide a broader and alternative theoretical validation for existing robust methods, whose validity was previously limited to specific forms of the residuals.
△ Less
Submitted 15 May, 2024;
originally announced May 2024.
-
Enhanced Whispering Gallery Mode Phase Shift using Indistinguishable Photon Pairs
Authors:
Callum Jones,
Antonio Vidiella-Barranco,
Jolly Xavier,
Frank Vollmer
Abstract:
We present a theoretical investigation of a whispering gallery mode (WGM) resonator coupled to a Mach-Zehnder interferometer (MZI) and show a bimodal coincidence transmission spectrum when the input state is an indistinguishable photon pair. This is due to the doubled WGM phase shift experienced by the path-entangled state in the interferometer. Further, we model the noise in a WGM resonance shift…
▽ More
We present a theoretical investigation of a whispering gallery mode (WGM) resonator coupled to a Mach-Zehnder interferometer (MZI) and show a bimodal coincidence transmission spectrum when the input state is an indistinguishable photon pair. This is due to the doubled WGM phase shift experienced by the path-entangled state in the interferometer. Further, we model the noise in a WGM resonance shift measurement comparing photon pairs with a coherent state. At least a four-fold improvement in the signal-to-noise ratio (SNR) is possible, with clear implications for quantum-enhanced WGM sensing.
△ Less
Submitted 26 April, 2024;
originally announced April 2024.
-
Assessing Engraftment Following Fecal Microbiota Transplant
Authors:
Chloe Herman,
Bridget M. Barker,
Thais F. Bartelli,
Vidhi Chandra,
Rosa Krajmalnik-Brown,
Mary Jewell,
Le Li,
Chen Liao,
Florencia McAllister,
Khemlal Nirmalkar,
Joao B. Xavier,
J. Gregory Caporaso
Abstract:
Fecal Microbiota Transplant (FMT) is an FDA approved treatment for recurrent Clostridium difficile infections, and is being explored for other clinical applications, from alleviating digestive and neurological disorders, to priming the microbiome for cancer treatment, and restoring microbiomes impacted by cancer treatment.
Quantifying the extent of engraftment following an FMT is important in de…
▽ More
Fecal Microbiota Transplant (FMT) is an FDA approved treatment for recurrent Clostridium difficile infections, and is being explored for other clinical applications, from alleviating digestive and neurological disorders, to priming the microbiome for cancer treatment, and restoring microbiomes impacted by cancer treatment.
Quantifying the extent of engraftment following an FMT is important in determining if a recipient didn't respond because the engrafted microbiome didn't produce the desired outcomes (a successful FMT, but negative treatment outcome), or the microbiome didn't engraft (an unsuccessful FMT and negative treatment outcome). The lack of a consistent methodology for quantifying FMT engraftment extent hinders the assessment of FMT success and its relation to clinical outcomes, and presents challenges for comparing FMT results and protocols across studies.
Here we review 46 studies of FMT in humans and model organisms and group their approaches for assessing the extent to which an FMT engrafts into three criteria: 1) Chimeric Asymmetric Community Coalescence investigates microbiome shifts following FMT engraftment. 2) Donated Microbiome Indicator Features tracks donated microbiome features as a signal of engraftment with methods such as differential abundance testing based on the current sample collection, or tracking changes in feature abundances that have been previously identified. 3) Temporal Stability examines how resistant post-FMT recipient's microbiomes are to reverting back to their baseline microbiome. Investigated together, these criteria provide a clear assessment of microbiome engraftment.
We discuss the pros and cons of each of these criteria, providing illustrative examples of their application. We also introduce key terminology and recommendations on how FMT studies can be analyzed for rigorous engraftment extent assessment.
△ Less
Submitted 10 April, 2024;
originally announced April 2024.
-
Secure quantum imaging with decoy state heralded single photons
Authors:
Siddhant Vernekar,
Jolly Xavier
Abstract:
Weak coherent source (WCS) and spontaneous parametric down converted heralded single photon pairs have found applications in quantum key distribution (QKD) and quantum imaging (QI) experiments. Decoy state methods have also been used to enhance the security for QKD and QI. We study quantum secured imaging with the decoy state heralded single photon source (HSPS). The HSPSs superior performance in…
▽ More
Weak coherent source (WCS) and spontaneous parametric down converted heralded single photon pairs have found applications in quantum key distribution (QKD) and quantum imaging (QI) experiments. Decoy state methods have also been used to enhance the security for QKD and QI. We study quantum secured imaging with the decoy state heralded single photon source (HSPS). The HSPSs superior performance in low photon number regimes makes it an ideal candidate for integrating quantum key distribution protocols to reduce measurement uncertainty and ensure secure QI. Furthermore, our results also infer the influence of the decoy state WCS, due to its higher operating speed than decoy state HSPS, would be effective in conditions that allow higher mean photon numbers for quantum secured imaging.
△ Less
Submitted 18 February, 2024;
originally announced February 2024.
-
Real-time imaging of standing-wave patterns in microresonators
Authors:
Haochen Yan,
Alekhya Ghosh,
Arghadeep Pal,
Hao Zhang,
Toby Bi,
George Ghalanos,
Shuangyou Zhang,
Lewis Hill,
Yaojing Zhang,
Yongyong Zhuang,
Jolly Xavier,
Pascal DelHaye
Abstract:
Real-time characterization of microresonator dynamics is important for many applications. In particular it is critical for near-field sensing and understanding light-matter interactions. Here, we report camera-facilitated imaging and analysis of standing wave patterns in optical ring resonators. The standing wave pattern is generated through bi-directional pumping of a microresonator and the scatt…
▽ More
Real-time characterization of microresonator dynamics is important for many applications. In particular it is critical for near-field sensing and understanding light-matter interactions. Here, we report camera-facilitated imaging and analysis of standing wave patterns in optical ring resonators. The standing wave pattern is generated through bi-directional pumping of a microresonator and the scattered light from the microresonator is collected by a short-wave infrared (SWIR) camera. The recorded scattering patterns are wavelength dependent, and the scattered intensity exhibits a linear relation with the circulating power within the microresonator. By modulating the relative phase between the two pump waves, we can control the generated standing waves movements and characterize the resonator with the SWIR camera. The visualized standing wave enables subwavelength distance measurements of scattering targets with nanometer-level accuracy. This work opens new avenues for applications in on-chip near-field (bio-)sensing, real time characterization of photonic integrated circuits and backscattering control in telecom systems.
△ Less
Submitted 15 January, 2024;
originally announced January 2024.
-
Hot Electron Generation and Manipulation in Nano-spiked Plasmonic Cavity Arrays
Authors:
Siddhartha Banerjee,
Jolly Xavier
Abstract:
The generation of the localized surface plasmon resonance (LSPR) on the surface of plasmonic structures in the nanoscale has paved the way for advanced biosensing, surpassing the conventional detection limits. The electric field enhancement (electromagnetic hot spots) between two plasmonic nano structures at close quarters produces hot electrons with a change in the electron density of the materia…
▽ More
The generation of the localized surface plasmon resonance (LSPR) on the surface of plasmonic structures in the nanoscale has paved the way for advanced biosensing, surpassing the conventional detection limits. The electric field enhancement (electromagnetic hot spots) between two plasmonic nano structures at close quarters produces hot electrons with a change in the electron density of the material. Techniques such as photoinjection are suitable to inject hot electrons from the metal crossing the Schottky barrier to the semiconductor leading to the development of photocurrent. sea urchin looking spiked nanoparticles serve a great interest in research of plasmonic materials as they can generate hot electrons at higher rate compared to other structures. By means of finite element method (FEM) analysis, we investigate the electric field enhancement generated between the cavity of nano star spike pairs and extend the study to an array of them. We further do a parametric study and identify the importance of using such an array to generate hot electrons whose applications are envisaged for light harvesting, enhanced photodetection, photo catalysis etc.
△ Less
Submitted 10 December, 2023;
originally announced December 2023.
-
A Multi-Token Coordinate Descent Method for Semi-Decentralized Vertical Federated Learning
Authors:
Pedro Valdeira,
Yuejie Chi,
Cláudia Soares,
João Xavier
Abstract:
Communication efficiency is a major challenge in federated learning (FL). In client-server schemes, the server constitutes a bottleneck, and while decentralized setups spread communications, they do not necessarily reduce them due to slower convergence. We propose Multi-Token Coordinate Descent (MTCD), a communication-efficient algorithm for semi-decentralized vertical federated learning, exploiti…
▽ More
Communication efficiency is a major challenge in federated learning (FL). In client-server schemes, the server constitutes a bottleneck, and while decentralized setups spread communications, they do not necessarily reduce them due to slower convergence. We propose Multi-Token Coordinate Descent (MTCD), a communication-efficient algorithm for semi-decentralized vertical federated learning, exploiting both client-server and client-client communications when each client holds a small subset of features. Our multi-token method can be seen as a parallel Markov chain (block) coordinate descent algorithm and it subsumes the client-server and decentralized setups as special cases. We obtain a convergence rate of $\mathcal{O}(1/T)$ for nonconvex objectives when tokens roam over disjoint subsets of clients and for convex objectives when they roam over possibly overlapping subsets. Numerical results show that MTCD improves the state-of-the-art communication efficiency and allows for a tunable amount of parallel communications.
△ Less
Submitted 18 September, 2023;
originally announced September 2023.
-
What it takes to solve the Origin(s) of Life: An integrated review of techniques
Authors:
OoLEN,
Silke Asche,
Carla Bautista,
David Boulesteix,
Alexandre Champagne-Ruel,
Cole Mathis,
Omer Markovitch,
Zhen Peng,
Alyssa Adams,
Avinash Vicholous Dass,
Arnaud Buch,
Eloi Camprubi,
Enrico Sandro Colizzi,
Stephanie Colón-Santos,
Hannah Dromiack,
Valentina Erastova,
Amanda Garcia,
Ghjuvan Grimaud,
Aaron Halpern,
Stuart A Harrison,
Seán F. Jordan,
Tony Z Jia,
Amit Kahana,
Artemy Kolchinsky,
Odin Moron-Garcia
, et al. (13 additional authors not shown)
Abstract:
Understanding the origin(s) of life (OoL) is a fundamental challenge for science in the 21st century. Research on OoL spans many disciplines, including chemistry, physics, biology, planetary sciences, computer science, mathematics and philosophy. The sheer number of different scientific perspectives relevant to the problem has resulted in the coexistence of diverse tools, techniques, data, and sof…
▽ More
Understanding the origin(s) of life (OoL) is a fundamental challenge for science in the 21st century. Research on OoL spans many disciplines, including chemistry, physics, biology, planetary sciences, computer science, mathematics and philosophy. The sheer number of different scientific perspectives relevant to the problem has resulted in the coexistence of diverse tools, techniques, data, and software in OoL studies. This has made communication between the disciplines relevant to the OoL extremely difficult because the interpretation of data, analyses, or standards of evidence can vary dramatically. Here, we hope to bridge this wide field of study by providing common ground via the consolidation of tools and techniques rather than positing a unifying view on how life emerges. We review the common tools and techniques that have been used significantly in OoL studies in recent years. In particular, we aim to identify which information is most relevant for comparing and integrating the results of experimental analyses into mathematical and computational models. This review aims to provide a baseline expectation and understanding of technical aspects of origins research, rather than being a primer on any particular topic. As such, it spans broadly -- from analytical chemistry to mathematical models -- and highlights areas of future work that will benefit from a multidisciplinary approach to tackling the mystery of life's origin. Ultimately, we hope to empower a new generation of OoL scientists by reviewing how they can investigate life's origin, rather than dictating how to think about the problem.
△ Less
Submitted 24 August, 2023; v1 submitted 22 August, 2023;
originally announced August 2023.
-
Effect of long-range hopping on dynamic quantum phase transitions of an exactly solvable free-fermion model: Nonanalyticities at almost all times
Authors:
J. C. Xavier,
José A. Hoyos
Abstract:
In this work, we investigate quenches in a free-fermion chain with long-range hopping which decay with the distance with an exponent $ν$ and has range $D$. By exploring the exact solution of the model, we found that the dynamic free energy is non-analytical, in the thermodynamic limit, whenever the sudden quench crosses the equilibrium quantum critical point. We were able to determine the non-anal…
▽ More
In this work, we investigate quenches in a free-fermion chain with long-range hopping which decay with the distance with an exponent $ν$ and has range $D$. By exploring the exact solution of the model, we found that the dynamic free energy is non-analytical, in the thermodynamic limit, whenever the sudden quench crosses the equilibrium quantum critical point. We were able to determine the non-analyticities of dynamic free energy $f(t)$ at some critical times $t^{c}$ by solving nonlinear equations. We also show that the Yang-Lee-Fisher (YLF) zeros cross the real-time axis at those critical times. We found that the number of nontrivial critical times, $N_{s},$ depends on $ν$ and $D$. In particular, we show that for small $ν$ and large $D$ the dynamic free energy presents non-analyticities in any time interval $Δt\sim1/D\ll1$, i.e., there are \emph{non-analyticities at almost all times}. For the spacial case $ν=0$, we obtain the critical times in terms of a simple expression of the model parameters and also show that $f(t)$ is non-analytical even for finite system under anti-periodic boundary condition, when we consider some special values of quench parameters. We also show that, generically, the first derivative of the dynamic free energy is discontinuous at the critical time instant when the YLF zeros are non-degenerate. On the other hand, when they become degenerate, all derivatives of $f(t)$ exist at the associated critical instant.
△ Less
Submitted 5 December, 2023; v1 submitted 9 August, 2023;
originally announced August 2023.
-
Mapping and maneuvering long-term natural orbits around Titania, a satellite of Uranus
Authors:
Silvia Giuliatti Winter,
Jadilene Xavier,
Antônio Bertachini Prado,
Andre Amarante
Abstract:
In this work, we present the results of a set of numerical simulations carried out to obtain long-duration orbits for a probe around Titania, Uranus' largest satellite. We also propose orbital maneuvers to extend the lifetime of some orbits. Titania's $J_2$ and $C_{22}$ gravitational coefficients and Uranus' gravitational perturbation are considered. The analysis of lifetime sensitivity due to pos…
▽ More
In this work, we present the results of a set of numerical simulations carried out to obtain long-duration orbits for a probe around Titania, Uranus' largest satellite. We also propose orbital maneuvers to extend the lifetime of some orbits. Titania's $J_2$ and $C_{22}$ gravitational coefficients and Uranus' gravitational perturbation are considered. The analysis of lifetime sensitivity due to possible errors in $J_2$ and $C_{22}$ values is investigated using multiple regression models. Simulations were performed for eccentricity equal 10-4, and lifetime maps were constructed. The results show that low-altitude orbits have longer lifetimes due to the balance between the disturbance of Uranus and the gravitational coefficients of Titania. The results also show that non-zero values of periapsis longitude ($ω$) and ascending node longitude ($Ω$) are essential to increase lifespan. Furthermore, the results indicate that the most economical maneuver occurs for final orbits of radius equal to 1050 km, this is observed for all inclination values.
△ Less
Submitted 13 July, 2023;
originally announced July 2023.
-
Robust Target Localization in 2D: A Value-at-Risk Approach
Authors:
João Domingos,
João Xavier
Abstract:
This paper consider considers the problem of locating a two dimensional target from range-measurements containing outliers. Assuming that the number of outlier is known, we formulate the problem of minimizing inlier losses while ignoring outliers. This leads to a combinatorial, non-convex, non-smooth problem involving the percentile function. Using the framework of risk analysis from Rockafellar e…
▽ More
This paper consider considers the problem of locating a two dimensional target from range-measurements containing outliers. Assuming that the number of outlier is known, we formulate the problem of minimizing inlier losses while ignoring outliers. This leads to a combinatorial, non-convex, non-smooth problem involving the percentile function. Using the framework of risk analysis from Rockafellar et al., we start by interpreting this formulation as a Value-at-risk (VaR) problem from portfolio optimization. To the best of our knowledge, this is the first time that a localization problem was formulated using risk analysis theory. To study the VaR formulation, we start by designing a majorizer set that contains any solution of a general percentile problem. This set is useful because, when applied to a localization scenario in 2D, it allows to majorize the solution set in terms of singletons, circumferences, ellipses and hyperbolas. Using know parametrization of these curves, we propose a grid method for the original non-convex problem. So we reduce the task of optimizing the VaR objective to that of efficiently sampling the proposed majorizer set. We compare our algorithm with four benchmarks in target localization. Numerical simulations show that our method is fast while, on average, improving the accuracy of the best benchmarks by at least 100m in a 1 Km$^2$ area.
△ Less
Submitted 5 July, 2023; v1 submitted 2 July, 2023;
originally announced July 2023.
-
Quotients of torus endomorphisms have parabolic orbifolds
Authors:
Sofía Llavayol,
Juliana Xavier
Abstract:
In this work we show that every quotient of a torus endomorphism has a parabolic orbifold. This answers a question of Mario Bonk and Daniel Meyer posed in their book "Expanding Thurston maps".
In this work we show that every quotient of a torus endomorphism has a parabolic orbifold. This answers a question of Mario Bonk and Daniel Meyer posed in their book "Expanding Thurston maps".
△ Less
Submitted 22 May, 2023;
originally announced May 2023.
-
Distributed detection of ARMA signals
Authors:
João Domingos,
João Xavier
Abstract:
This paper considers a distributed detection setup where agents in a network want to detect a time-varying signal embedded in temporally correlated noise. The signal of interest is the impulse response of an ARMA (auto-regressive moving average) filter, and the noise is the output of yet another ARMA filter which is fed white Gaussian noise. For this extended problem setup, which can prompt novel…
▽ More
This paper considers a distributed detection setup where agents in a network want to detect a time-varying signal embedded in temporally correlated noise. The signal of interest is the impulse response of an ARMA (auto-regressive moving average) filter, and the noise is the output of yet another ARMA filter which is fed white Gaussian noise. For this extended problem setup, which can prompt novel behaviour, we propose a comprehensive solution. First, we extend the well-known running consensus detector (RCD) to this correlated setup; then, we design an efficient implementation of the RCD by exploiting the underlying ARMA structures; and, finally, we derive the theoretical asymptotic performance of the RCD in this ARMA setup. It turns out that the error probability at each agent exhibits one of two regimes: either (a) the error probability decays exponentially fast to zero or (b) it converges to a strictly positive error floor. While regime (a) spans staple results in large deviation theory, regime (b) is new in distributed detection and is elicited by the ARMA setup. We fully characterize these two scenarios: we give necessary and sufficient conditions, phrased in terms of the zero and poles of the underlying ARMA models, for the emergence of each regime, and provide closed-form expressions for both the decay rates of regime (a) and the positive error floors of regime (b). Our analysis also shows that the ARMA setup leads to two novel features: (1) the threshold level used in RCD can influence the asymptotics of the error probabilities and (2) some agents might be weakly informative, in the sense that their observations do not improve the asymptotic performance of RCD and, as such, can be safely muted to save sensing resources. Numerical simulations illustrate and confirm the theoretical findings.
△ Less
Submitted 14 April, 2023;
originally announced April 2023.
-
Time-dependent Mandel Q parameter analysis for a hexagonal boron nitride single photon source
Authors:
Callum Jones,
Jolly Xavier,
Samir Vartabi Kashanian,
Minh Nguyen,
Igor Aharonovich,
Frank Vollmer
Abstract:
The time-dependent Mandel Q parameter, Q(T), provides a measure of photon number variance for a light source as a function of integration time. Here, we use Q(T) to characterise single photon emission from a quantum emitter in hexagonal boron nitride (hBN). Under pulsed excitation a negative Q parameter was measured, indicating photon antibunching at an integration time of 100 ns. For larger integ…
▽ More
The time-dependent Mandel Q parameter, Q(T), provides a measure of photon number variance for a light source as a function of integration time. Here, we use Q(T) to characterise single photon emission from a quantum emitter in hexagonal boron nitride (hBN). Under pulsed excitation a negative Q parameter was measured, indicating photon antibunching at an integration time of 100 ns. For larger integration times Q is positive and the photon statistics become super-Poissonian, and we show by comparison with a Monte Carlo simulation for a three-level emitter that this is consistent with the effect of a metastable shelving state. Looking towards technological applications for hBN single photon sources, we propose that Q(T) provides valuable information on the intensity stability of single photon emission. This is useful in addition to the commonly used $g^{(2)}(τ)$ function for the complete characterisation of a hBN emitter.
△ Less
Submitted 24 January, 2023;
originally announced January 2023.
-
The growth rate inequality for Thurston maps with non hyperbolic orbifolds
Authors:
J. Iglesias,
A. Portela,
A. Rovella,
J. Xavier
Abstract:
Let $f: S^2 \to S^2$ be a continuous map of degree $d$, $|d|>1$, and let $N_nf$ denote the number of fixed points of $f^n$. We show that if $f$ is a Thurston map with non hyperbolic orbifold, then either the growth rate inequality $\limsup \frac{1}{n} \log N_nf\geq \log |d|$ holds for $f$ or $f$ has exactly two critical points which are fixed and totally invariant.
Let $f: S^2 \to S^2$ be a continuous map of degree $d$, $|d|>1$, and let $N_nf$ denote the number of fixed points of $f^n$. We show that if $f$ is a Thurston map with non hyperbolic orbifold, then either the growth rate inequality $\limsup \frac{1}{n} \log N_nf\geq \log |d|$ holds for $f$ or $f$ has exactly two critical points which are fixed and totally invariant.
△ Less
Submitted 7 November, 2022;
originally announced November 2022.
-
Mapping Long-Term Natural Orbits about Titania, a Satellite of Uranus
Authors:
J. Xavier,
A. B. Prado,
S. M. Giuliatti Winter,
A. Amarante
Abstract:
Close polar and circular orbits are of great interest for the exploration of natural satellites. There are still no studies in the literature investigating orbits around Titania, the largest satellite of Uranus. In this work, we present results of a set of numerical simulations carried out to obtain long-duration orbits for a probe around Titania. Through an expansion of the gravitational potentia…
▽ More
Close polar and circular orbits are of great interest for the exploration of natural satellites. There are still no studies in the literature investigating orbits around Titania, the largest satellite of Uranus. In this work, we present results of a set of numerical simulations carried out to obtain long-duration orbits for a probe around Titania. Through an expansion of the gravitational potential up to second order, the asymmetry of the gravitational field due to Titania's coefficient $C_{22}$, the zonal coefficient $J_2$, and the gravitational perturbation of Uranus is considered. The analysis of lifetime sensitivity due to possible errors in the values of $J_2$ and $C_{22}$ is investigated using multiple regression models. Simulations were performed for different eccentricity values, and lifetime maps were constructed. The results show that low-altitude and near-circular orbits have longer lifetimes due to the balance between the disturbance of Uranus and the gravitational coefficients of Titania. The results also show that non-zero values of the longitude of periapsis ($ω$) and longitude of the ascending node ($Ω$) are essential to increase the lifetime up to eight times compared to cases where $ω= Ω=0^\circ$. We also show that an orbit with eccentricity $10^{-3}$ is the most affected by errors in the values of $J_2$ and $C_{22}$.
△ Less
Submitted 27 March, 2022;
originally announced March 2022.
-
Decentralized EM to Learn Gaussian Mixtures from Datasets Distributed by Features
Authors:
Pedro Valdeira,
Cláudia Soares,
João Xavier
Abstract:
Expectation Maximization (EM) is the standard method to learn Gaussian mixtures. Yet its classic, centralized form is often infeasible, due to privacy concerns and computational and communication bottlenecks. Prior work dealt with data distributed by examples, horizontal partitioning, but we lack a counterpart for data scattered by features, an increasingly common scheme (e.g. user profiling with…
▽ More
Expectation Maximization (EM) is the standard method to learn Gaussian mixtures. Yet its classic, centralized form is often infeasible, due to privacy concerns and computational and communication bottlenecks. Prior work dealt with data distributed by examples, horizontal partitioning, but we lack a counterpart for data scattered by features, an increasingly common scheme (e.g. user profiling with data from multiple entities). To fill this gap, we provide an EM-based algorithm to fit Gaussian mixtures to Vertically Partitioned data (VP-EM). In federated learning setups, our algorithm matches the centralized EM fitting of Gaussian mixtures constrained to a subspace. In arbitrary communication graphs, consensus averaging allows VP-EM to run on large peer-to-peer networks as an EM approximation. This mismatch comes from consensus error only, which vanishes exponentially fast with the number of consensus rounds. We demonstrate VP-EM on various topologies for both synthetic and real data, evaluating its approximation of centralized EM and seeing that it outperforms the available benchmark.
△ Less
Submitted 24 January, 2022;
originally announced January 2022.
-
Disorder-induced dynamical Griffiths singularities after certain quantum quenches
Authors:
José A. Hoyos,
R. F. P. Costa,
J. C. Xavier
Abstract:
We demonstrate that in a class of disordered quantum systems the dynamical partition function is not an analytical function in a time window after certain quantum quenches. We related this behavior to rare and large regions with atypical inhomogeneity configurations. We also quantify the strength of the associated singularities and their signatures in experiments and numerical studies.
We demonstrate that in a class of disordered quantum systems the dynamical partition function is not an analytical function in a time window after certain quantum quenches. We related this behavior to rare and large regions with atypical inhomogeneity configurations. We also quantify the strength of the associated singularities and their signatures in experiments and numerical studies.
△ Less
Submitted 10 October, 2022; v1 submitted 24 January, 2022;
originally announced January 2022.
-
Coexistence of spontaneous dimerization and magnetic order in a transverse-field Ising ladder with four-spin interactions
Authors:
J. C. Xavier,
R. G. Pereira,
M. E. S. Nunes,
J. A. Plascak
Abstract:
The spin-1/2 transverse field two-leg Ising ladder with nearest-neighbor exchange and plaquette four-spin interaction $J_{4}$ is studied analytically and numerically with the density matrix renormalization group approach. The quantum phase diagram in the transverse field $B$ versus $J_{4}$ plane has been obtained. There are three different phases: a paramagnetic (PM) phase for high values of the t…
▽ More
The spin-1/2 transverse field two-leg Ising ladder with nearest-neighbor exchange and plaquette four-spin interaction $J_{4}$ is studied analytically and numerically with the density matrix renormalization group approach. The quantum phase diagram in the transverse field $B$ versus $J_{4}$ plane has been obtained. There are three different phases: a paramagnetic (PM) phase for high values of the transverse field and, for low values of $B$, a ferromagnetic (FM) ordered phase for small $J_{4}$ and a dimerized-rung (DR) phase for large negative values of $J_{4}$. All phases are separated by quantum phase transition lines meeting at a multicritical point. The critical lines have been obtained by exploring the entanglement entropy. The results show that along the critical lines the central charge is $c=1/2$, while at the multicritical point one has $c=1$. The scaling dimension of the energy operator is $X_ε=1$, in agreement with the universality class of the critical behavior of the quantum Ising chain. An effective field theory for the multicritical point is also discussed. The FM and the DR order parameters have also been computed and we found a region where the FM and the DR phases coexist.
△ Less
Submitted 22 November, 2021;
originally announced November 2021.
-
Unambiguous signature of exchange interactions between nanoparticles in a collective system
Authors:
C. A. Iglesias,
J. C. R. de Araújo,
J. Xavier,
R. B. da Silva,
J. M. Soares,
L. M. Silva,
J. H. de Araújo,
S. N. de Medeiros,
P. B. Souza,
C. C. Plá Cid,
M. Gamino,
M. A. Correa,
F. Bohn
Abstract:
We report a simple, efficient manner to assess magnetic interactions, more specifically the exchange interaction, in systems as blocked magnetic nanoparticles. Specifically, we investigate what is the theoretical limit for the external magnetic susceptibility in a system described by the plain old Stoner-Wohlfarth model. We go beyond and introduce a general mean field theory for interacting system…
▽ More
We report a simple, efficient manner to assess magnetic interactions, more specifically the exchange interaction, in systems as blocked magnetic nanoparticles. Specifically, we investigate what is the theoretical limit for the external magnetic susceptibility in a system described by the plain old Stoner-Wohlfarth model. We go beyond and introduce a general mean field theory for interacting systems, thus estimating how the magnetic susceptibility is affected due to the dipolar and exchange interactions inside the system. We disclose a fundamental inequality for the magnetic susceptibility and show its violation is an unambiguous signature of the existence of exchange interactions between nanoparticles in a collective system. To test the robustness of our theoretical achievements, we examine magnetization measurements and Henkel plots for blocked magnetic nanopowders. The agreement between experiment and theory provides evidence to confirm the validity of our findings.
△ Less
Submitted 8 October, 2021;
originally announced October 2021.
-
Robust Localization with Bounded Noise: Creating a Superset of the Possible Target Positions via Linear-Fractional Representations
Authors:
João Domingos,
Cláudia Soares,
João Xavier
Abstract:
Locating a target is key in many applications, namely in high-stakes real-world scenarios, like detecting humans or obstacles in vehicular networks. In scenarios where precise statistics of the measurement noise are unavailable, applications require localization methods that assume minimal knowledge on the noise distribution. We present a scalable algorithm delimiting a tight superset of all possi…
▽ More
Locating a target is key in many applications, namely in high-stakes real-world scenarios, like detecting humans or obstacles in vehicular networks. In scenarios where precise statistics of the measurement noise are unavailable, applications require localization methods that assume minimal knowledge on the noise distribution. We present a scalable algorithm delimiting a tight superset of all possible target locations, assuming range measurements to known landmarks, contaminated with bounded noise and unknown distributions. This superset is of primary interest in robust statistics since it is a tight majorizer of the set of Maximum-Likelihood (ML) estimates parametrized by noise densities respecting two main assumptions: (1) the noise distribution is supported on a ellipsoidal uncertainty region and (2) the measurements are non-negative with probability one. We create the superset through convex relaxations that use Linear Fractional Representations (LFRs), a well-known technique in robust control. For low noise regimes the supersets created by our method double the accuracy of a standard semidefinite relaxation. For moderate to high noise regimes our method still improves the benchmark but the benefit tends to be less significant, as both supersets tend to have the same size (area).
△ Less
Submitted 12 June, 2022; v1 submitted 6 October, 2021;
originally announced October 2021.
-
COCO Denoiser: Using Co-Coercivity for Variance Reduction in Stochastic Convex Optimization
Authors:
Manuel Madeira,
Renato Negrinho,
João Xavier,
Pedro M. Q. Aguiar
Abstract:
First-order methods for stochastic optimization have undeniable relevance, in part due to their pivotal role in machine learning. Variance reduction for these algorithms has become an important research topic. In contrast to common approaches, which rarely leverage global models of the objective function, we exploit convexity and L-smoothness to improve the noisy estimates outputted by the stochas…
▽ More
First-order methods for stochastic optimization have undeniable relevance, in part due to their pivotal role in machine learning. Variance reduction for these algorithms has become an important research topic. In contrast to common approaches, which rarely leverage global models of the objective function, we exploit convexity and L-smoothness to improve the noisy estimates outputted by the stochastic gradient oracle. Our method, named COCO denoiser, is the joint maximum likelihood estimator of multiple function gradients from their noisy observations, subject to co-coercivity constraints between them. The resulting estimate is the solution of a convex Quadratically Constrained Quadratic Problem. Although this problem is expensive to solve by interior point methods, we exploit its structure to apply an accelerated first-order algorithm, the Fast Dual Proximal Gradient method. Besides analytically characterizing the proposed estimator, we show empirically that increasing the number and proximity of the queried points leads to better gradient estimates. We also apply COCO in stochastic settings by plugging it in existing algorithms, such as SGD, Adam or STRSAGA, outperforming their vanilla versions, even in scenarios where our modelling assumptions are mismatched.
△ Less
Submitted 7 September, 2021;
originally announced September 2021.
-
Distributed Banach-Picard Iteration: Application to Distributed EM and Distributed PCA
Authors:
Francisco L. Andrade,
Mário A. T. Figueiredo,
João Xavier
Abstract:
In recent work, we proposed a distributed Banach-Picard iteration (DBPI) that allows a set of agents, linked by a communication network, to find a fixed point of a locally contractive (LC) map that is the average of individual maps held by said agents. In this work, we build upon the DBPI and its local linear convergence (LLC) guarantees to make several contributions. We show that Sanger's algorit…
▽ More
In recent work, we proposed a distributed Banach-Picard iteration (DBPI) that allows a set of agents, linked by a communication network, to find a fixed point of a locally contractive (LC) map that is the average of individual maps held by said agents. In this work, we build upon the DBPI and its local linear convergence (LLC) guarantees to make several contributions. We show that Sanger's algorithm for principal component analysis (PCA) corresponds to the iteration of an LC map that can be written as the average of local maps, each map known to each agent holding a subset of the data. Similarly, we show that a variant of the expectation-maximization (EM) algorithm for parameter estimation from noisy and faulty measurements in a sensor network can be written as the iteration of an LC map that is the average of local maps, each available at just one node. Consequently, via the DBPI, we derive two distributed algorithms - distributed EM and distributed PCA - whose LLC guarantees follow from those that we proved for the DBPI. The verification of the LC condition for EM is challenging, as the underlying operator depends on random samples, thus the LC condition is of probabilistic nature.
△ Less
Submitted 26 January, 2022; v1 submitted 20 June, 2021;
originally announced June 2021.
-
Distributed Banach-Picard Iteration for Locally Contractive Maps
Authors:
Francisco L. Andrade,
Mário A. T. Figueiredo,
João Xavier
Abstract:
The Banach-Picard iteration is widely used to find fixed points of locally contractive (LC) maps. This paper extends the Banach-Picard iteration to distributed settings; specifically, we assume the map of which the fixed point is sought to be the average of individual (not necessarily LC) maps held by a set of agents linked by a communication network. An additional difficulty is that the LC map is…
▽ More
The Banach-Picard iteration is widely used to find fixed points of locally contractive (LC) maps. This paper extends the Banach-Picard iteration to distributed settings; specifically, we assume the map of which the fixed point is sought to be the average of individual (not necessarily LC) maps held by a set of agents linked by a communication network. An additional difficulty is that the LC map is not assumed to come from an underlying optimization problem, which prevents exploiting strong global properties such as convexity or Lipschitzianity. Yet, we propose a distributed algorithm and prove its convergence, in fact showing that it maintains the linear rate of the standard Banach-Picard iteration for the average LC map. As another contribution, our proof imports tools from perturbation theory of linear operators, which, to the best of our knowledge, had not been used before in the theory of distributed computation.
△ Less
Submitted 28 December, 2021; v1 submitted 31 March, 2021;
originally announced April 2021.
-
Specific loss power of magnetic nanoparticles (fluid) hyperthermia in non-adiabatic conditions
Authors:
C. A. M. Iglesias,
J. C. R. de Araújo,
J. Xavier,
R. B. da Silva,
J. M. Soares,
E. L. Brito,
L. Streck,
J. L. C. Fonseca,
C. C. Plá Cid,
M. Gamino,
E. F. Silva,
C. Chesman,
M. A. Correa,
S. N. de Medeiros,
F. Bohn
Abstract:
We investigate the magnetic nanoparticles (fluid) hyperthermia in non-adiabatic conditions through the calorimetric method. Specifically, we propose a theoretical approach to magnetic hyperthermia from a thermodynamic point of view. To test the robustness of the approach, we perform hyperthermia experiments and analyze the thermal behavior of magnetite and magnesium ferrite magnetic nanoparticles…
▽ More
We investigate the magnetic nanoparticles (fluid) hyperthermia in non-adiabatic conditions through the calorimetric method. Specifically, we propose a theoretical approach to magnetic hyperthermia from a thermodynamic point of view. To test the robustness of the approach, we perform hyperthermia experiments and analyze the thermal behavior of magnetite and magnesium ferrite magnetic nanoparticles dispersed in water submitted to an alternating magnetic field. From our findings, besides estimating the specific loss power value from a non-adiabatic process, thus enhancing the accuracy in the determination of this quantity, we provide physical meaning to parameters found in literature that still remained not fully understood, and bring to light how they can be obtained experimentally.
△ Less
Submitted 9 March, 2021;
originally announced March 2021.
-
Branched coverings of the sphere having a completely invariant continuum with infinitely many Wada Lakes
Authors:
J. Iglesias,
A. Portela,
A. Rovella,
J. Xavier
Abstract:
We construct a family of smooth branched coverings of degree $2$ of the sphere $S^2$ having a completely invariant indecomposable continuum $K$ and infinitely many Wada Lakes.
We construct a family of smooth branched coverings of degree $2$ of the sphere $S^2$ having a completely invariant indecomposable continuum $K$ and infinitely many Wada Lakes.
△ Less
Submitted 6 December, 2022; v1 submitted 18 June, 2020;
originally announced June 2020.
-
Confinement and bound states of bound states in a transverse-field two-leg Ising ladder
Authors:
Flavia B. Ramos,
Mate Lencses,
J. C. Xavier,
Rodrigo G. Pereira
Abstract:
Weakly coupled Ising chains provide a condensed-matter realization of confinement. In these systems, kinks and antikinks bind into mesons due to an attractive interaction potential that increases linearly with the distance between the particles. While single mesons have been directly observed in experiments, the role of the multiparticle continuum and bound states of mesons in the excitation spect…
▽ More
Weakly coupled Ising chains provide a condensed-matter realization of confinement. In these systems, kinks and antikinks bind into mesons due to an attractive interaction potential that increases linearly with the distance between the particles. While single mesons have been directly observed in experiments, the role of the multiparticle continuum and bound states of mesons in the excitation spectrum is far less clear. Using time-dependent density matrix renormalization group methods, we study the dynamical structure factors of one- and two-spin operators in a transverse-field two-leg Ising ladder in the ferromagnetic phase. The propagation of time-dependent correlations and the two-spin excitation spectrum reveal the existence of interchain bound states, which are absent in the one-spin dynamical structure factor. We also identify two-meson bound states that appear at higher energies, above the thresholds of several two-meson continua.
△ Less
Submitted 17 July, 2020; v1 submitted 6 May, 2020;
originally announced May 2020.
-
StochaLM: a Stochastic alternate Linearization Method for distributed optimization
Authors:
Inês Almeida,
João Xavier
Abstract:
We present the Stochastic alternate Linearization Method (StochaLM), a token-based method for distributed optimization. This algorithm finds the solution of a consensus optimization problem by solving a sequence of subproblems where some components of the cost are linearized around specific anchor points. StochaLM can be interpreted as a dual block coordinate ascent method whose block components a…
▽ More
We present the Stochastic alternate Linearization Method (StochaLM), a token-based method for distributed optimization. This algorithm finds the solution of a consensus optimization problem by solving a sequence of subproblems where some components of the cost are linearized around specific anchor points. StochaLM can be interpreted as a dual block coordinate ascent method whose block components are selected using the state of an ergodic Markov chain. This sampling process is neither essentially cyclic nor independent over time, preventing us from using proofs of convergence of dual block coordinate ascent methods done in previous works. The proof of convergence of our method is, therefore, also novel. We show that, if the cost is strongly convex and the network is fully connected, then, with probability one, the primal sequence generated by StochaLM converges to the solution. Our method is application-friendly, as it has no global hyperparameters to tune; any hyperparameters can be tuned locally by the agents, using information regarding their private cost and neighbourhood. Our method is, therefore, decentralized in the truest sense. Numerical experiments evidence that our method converges to the solution faster than other token-based methods, even when these methods' hyperparameters are tuned for optimal performance.
△ Less
Submitted 27 December, 2021; v1 submitted 23 April, 2020;
originally announced April 2020.
-
Entanglement and boundary entropy in quantum spin chains with arbitrary direction of the boundary magnetic fields
Authors:
J. C. Xavier,
M. A. Rajabpour
Abstract:
We calculate the entanglement and the universal boundary entropy (BE) in the critical quantum spin chains, such as the transverse field Ising chain and the XXZ chain, with arbitrary direction of the boundary magnetic field (ADBMF). We determine the boundary universality class that an ADBMF induces. In particular, we show that the induced boundary conformal field theory (BCFT) depends on the point…
▽ More
We calculate the entanglement and the universal boundary entropy (BE) in the critical quantum spin chains, such as the transverse field Ising chain and the XXZ chain, with arbitrary direction of the boundary magnetic field (ADBMF). We determine the boundary universality class that an ADBMF induces. In particular, we show that the induced boundary conformal field theory (BCFT) depends on the point on the Bloch sphere where the boundary magnetic field directs. We show that the classification of the directions boils down to the simple fact that the boundary field breaks the bulk symmetry or does not. We present a procedure to estimate the universal BE, based on the finite-size corrections of the entanglement entropy, that apply to the ADBMF. To calculate the universal BE in the XXZ chain, we use the density matrix renormalization group (DMRG). The transverse field XY chain with ADBMF after Jordan-Wigner (JW) transformation is not a quadratic free fermion Hamiltonian. We map this model to a quadratic free fermion chain by introducing two extra ancillary spins coupled to the main chain at the boundaries, which makes the problem {\it{integrable}}. The eigenstates of the transverse field XY chain can be obtained by proper projection in the enlarged chain. Using this mapping, we are able to calculate the entanglement entropy of the transverse field XY chain using the usual correlation matrix technique up to relatively large sizes.
△ Less
Submitted 15 June, 2020; v1 submitted 28 February, 2020;
originally announced March 2020.
-
Modeling microbial cross-feeding at intermediate scale portrays community dynamics and species coexistence
Authors:
Chen Liao,
Tong Wang,
Sergei Maslov,
Joao B. Xavier
Abstract:
Social interaction between microbes can be described at many levels of details, ranging from the biochemistry of cell-cell interactions to the ecological dynamics of populations. Choosing the best level to model microbial communities without losing generality remains a challenge. Here we propose to model cross-feeding interactions at an intermediate level between genome-scale metabolic models of i…
▽ More
Social interaction between microbes can be described at many levels of details, ranging from the biochemistry of cell-cell interactions to the ecological dynamics of populations. Choosing the best level to model microbial communities without losing generality remains a challenge. Here we propose to model cross-feeding interactions at an intermediate level between genome-scale metabolic models of individual species and consumer-resource models of ecosystems, which is suitable to empirical data. We applied our method to three published examples of multi-strain Escherichia coli communities with increasing complexity consisting of uni-, bi-, and multi-directional cross-feeding of either substitutable metabolic byproducts or essential nutrients. The intermediate-scale model accurately described empirical data and could quantify exchange rates elusive by other means, such as the byproduct secretions, even for a complex community of 14 amino acid auxotrophs. We used the three models to study each community's limits of robustness to perturbations such as variations in resource supply, antibiotic treatments and invasion by other "cheaters" species. Our analysis provides a foundation to quantify cross-feeding interactions from experimental data, and highlights the importance of metabolic exchanges in the dynamics and stability of microbial communities.
△ Less
Submitted 19 February, 2020;
originally announced February 2020.
-
Primal-dual methods for large-scale and distributed convex optimization and data analytics
Authors:
Dusan Jakovetic,
Dragana Bajovic,
Joao Xavier,
Jose M. F. Moura
Abstract:
The augmented Lagrangian method (ALM) is a classical optimization tool that solves a given "difficult" (constrained) problem via finding solutions of a sequence of "easier"(often unconstrained) sub-problems with respect to the original (primal) variable, wherein constraints satisfaction is controlled via the so-called dual variables. ALM is highly flexible with respect to how primal sub-problems c…
▽ More
The augmented Lagrangian method (ALM) is a classical optimization tool that solves a given "difficult" (constrained) problem via finding solutions of a sequence of "easier"(often unconstrained) sub-problems with respect to the original (primal) variable, wherein constraints satisfaction is controlled via the so-called dual variables. ALM is highly flexible with respect to how primal sub-problems can be solved, giving rise to a plethora of different primal-dual methods. The powerful ALM mechanism has recently proved to be very successful in various large scale and distributed applications. In addition, several significant advances have appeared, primarily on precise complexity results with respect to computational and communication costs in the presence of inexact updates and design and analysis of novel optimal methods for distributed consensus optimization. We provide a tutorial-style introduction to ALM and its variants for solving convex optimization problems in large scale and distributed settings. We describe control-theoretic tools for the algorithms' analysis and design, survey recent results, and provide novel insights in the context of two emerging applications: federated learning and distributed energy trading.
△ Less
Submitted 14 April, 2020; v1 submitted 18 December, 2019;
originally announced December 2019.
-
Linearization of topologically Anosov homeomorphisms of non compact surfaces
Authors:
Gonzalo Cousillas,
Jorge Groisman,
Juliana Xavier
Abstract:
We study the dynamics of Topologically Anosov homeomorphisms of non compact surfaces. In the case of surfaces of genus zero and finite type, we classify them. We prove that if $f:S \to S$, is a Topologically Anosov homeomorphism where $S$ is a non-compact surface of genus zero and finite type, then $S= \R ^ 2$ and $f$ is conjugate to a homothety or reverse homothety (depending on wether $f$ preser…
▽ More
We study the dynamics of Topologically Anosov homeomorphisms of non compact surfaces. In the case of surfaces of genus zero and finite type, we classify them. We prove that if $f:S \to S$, is a Topologically Anosov homeomorphism where $S$ is a non-compact surface of genus zero and finite type, then $S= \R ^ 2$ and $f$ is conjugate to a homothety or reverse homothety (depending on wether $f$ preserves or reverses orientation). A weaker version of this result was conjectured in a previous work.
△ Less
Submitted 24 September, 2019;
originally announced September 2019.
-
Distortion analysis of crystalline and locally quasicrystalline 2D photonic structures with GISAXS
Authors:
Mika Pflüger,
Victor Soltwisch,
Jolly Xavier,
Jürgen Probst,
Frank Scholze,
Christiane Becker,
Michael Krumrey
Abstract:
In this study, grazing incidence small-angle X-ray scattering (GISAXS) is used to collect statistical information on dimensional parameters in an area of 20 mm x 15 mm on photonic structures produced by nanoimprint lithography. The photonic structures are composed of crystalline and locally quasicrystalline two-dimensional patterns with structure sizes between about 100 nm and 10 $μ$m to enable br…
▽ More
In this study, grazing incidence small-angle X-ray scattering (GISAXS) is used to collect statistical information on dimensional parameters in an area of 20 mm x 15 mm on photonic structures produced by nanoimprint lithography. The photonic structures are composed of crystalline and locally quasicrystalline two-dimensional patterns with structure sizes between about 100 nm and 10 $μ$m to enable broadband visible light absorption for use in solar energy harvesting. These first GISAXS measurements on locally quasicrystalline samples demonstrate that GISAXS is capable of showing the locally quasicrystalline nature of the samples while at the same time revealing the long-range periodicity introduced due to the lattice design. We describe the scattering qualitatively in the framework of the distorted wave Born approximation using a hierarchical model mirroring the sample design, which consists of a rectangular and locally quasicrystalline supercell which is repeated periodically to fill the whole surface.
The nanoimprinted samples are compared to a sample manufactured using electron beam lithography and the distortions of the periodic and locally quasiperiodic samples are quantified statistically. Due to the high sensitivity of GISAXS to deviations from the perfect lattice, the misalignment of the crystallographic axes was measured with a resolution of 0.015°, showing distortions up to +/- 0.15° in the investigated samples.
△ Less
Submitted 19 March, 2019;
originally announced March 2019.
-
Adaptive Density Matrix Renormalization Group for Disordered Systems
Authors:
J. C. Xavier,
J. A. Hoyos,
E. Miranda
Abstract:
We propose a simple modification of the density matrix renormalization group (DMRG) method in order to tackle strongly disordered quantum spin chains. Our proposal, akin to the idea of the adaptive time-dependent DMRG, enables us to reach larger system sizes in the strong disorder limit by avoiding most of the metastable configurations which hinder the performance of the standard DMRG method. We b…
▽ More
We propose a simple modification of the density matrix renormalization group (DMRG) method in order to tackle strongly disordered quantum spin chains. Our proposal, akin to the idea of the adaptive time-dependent DMRG, enables us to reach larger system sizes in the strong disorder limit by avoiding most of the metastable configurations which hinder the performance of the standard DMRG method. We benchmark our adaptive method by revisiting the random antiferromagnetic XXZ spin-1/2 chain for which we compute the random-singlet ground-state average spin-spin correlation functions and von Neumann entanglement entropy. We then apply our method to the bilinear-biquadratic random antiferromagnetic spin-1 chain tuned to the antiferromagnet and gapless highly symmetric SU(3) point. We find the new result that the mean correlation function decays algebraically with the same universal exponent $φ=2$ as the spin-1/2 chain. We then perform numerical and analytical strong-disorder renormalization-group calculations which confirm this finding and generalize it for any highly symmetric SU($N$) random-singlet state.
△ Less
Submitted 12 November, 2018; v1 submitted 11 September, 2018;
originally announced September 2018.
-
Latent Space Temporal Model of Microbial Abundance to Predict Domination and Bacteremia
Authors:
Ruiqi Zhong,
Tyler Joseph,
Joao B Xavier,
Itsik Pe'er
Abstract:
Gut microbial composition has been linked to multiple health outcomes. Yet, temporal analysis of this composition had been limited to deterministic models. In this paper, we introduce a probabilistic model for the dynamics of intestinal microbiomes that takes into account interaction among bacteria as well as external effects such as antibiotics. The model successfully deals with pragmatic issues…
▽ More
Gut microbial composition has been linked to multiple health outcomes. Yet, temporal analysis of this composition had been limited to deterministic models. In this paper, we introduce a probabilistic model for the dynamics of intestinal microbiomes that takes into account interaction among bacteria as well as external effects such as antibiotics. The model successfully deals with pragmatic issues such as random measurement error and varying time intervals between measurements through latent space modeling. We demonstrate utility of the model by using latent state features to predict the clinical events of intestinal domination and bacteremia, improving accuracy over existing methods. We further leverage this framework to validate known links between antibiotics and clinical outcomes, while discovering new ones.
△ Less
Submitted 31 August, 2018;
originally announced August 2018.
-
FADE: Fast and Asymptotically efficient Distributed Estimator for dynamic networks
Authors:
António Simões,
João Xavier
Abstract:
Consider a set of agents that wish to estimate a vector of parameters of their mutual interest. For this estimation goal, agents can sense and communicate. When sensing, an agent measures (in additive gaussian noise) linear combinations of the unknown vector of parameters. When communicating, an agent can broadcast information to a few other agents, by using the channels that happen to be randomly…
▽ More
Consider a set of agents that wish to estimate a vector of parameters of their mutual interest. For this estimation goal, agents can sense and communicate. When sensing, an agent measures (in additive gaussian noise) linear combinations of the unknown vector of parameters. When communicating, an agent can broadcast information to a few other agents, by using the channels that happen to be randomly at its disposal at the time.
To coordinate the agents towards their estimation goal, we propose a novel algorithm called FADE (Fast and Asymptotically efficient Distributed Estimator), in which agents collaborate at discrete time-steps; at each time-step, agents sense and communicate just once, while also updating their own estimate of the unknown vector of parameters.
FADE enjoys five attractive features: first, it is an intuitive estimator, simple to derive; second, it withstands dynamic networks, that is, networks whose communication channels change randomly over time; third, it is strongly consistent in that, as time-steps play out, each agent's local estimate converges (almost surely) to the true vector of parameters; fourth, it is both asymptotically unbiased and efficient, which means that, across time, each agent's estimate becomes unbiased and the mean-square error (MSE) of each agent's estimate vanishes to zero at the same rate of the MSE of the optimal estimator at an almighty central node; fifth, and most importantly, when compared with a state-of-art consensus+innovation (CI) algorithm, it yields estimates with outstandingly lower mean-square errors, for the same number of communications -- for example, in a sparsely connected network model with 50 agents, we find through numerical simulations that the reduction can be dramatic, reaching several orders of magnitude.
△ Less
Submitted 31 July, 2018;
originally announced July 2018.
-
Distributed Estimation Via a Roaming Token
Authors:
Lucas Balthazar,
João Xavier,
Bruno Sinopoli
Abstract:
We present an algorithm for the problem of linear distributed estimation of a parameter in a network where a set of agents are successively taking measurements. The approach considers a roaming token in a network that carries the estimate, and jumps from one agent to another in its vicinity according to the probabilities of a Markov chain. When the token is at an agent it records the agent's local…
▽ More
We present an algorithm for the problem of linear distributed estimation of a parameter in a network where a set of agents are successively taking measurements. The approach considers a roaming token in a network that carries the estimate, and jumps from one agent to another in its vicinity according to the probabilities of a Markov chain. When the token is at an agent it records the agent's local information. We analyze the proposed algorithm and show that it is consistent and asymptotically optimal, in the sense that its mean-square-error (MSE) rate of decay approaches the centralized one as the number of iterations increases. We show these results for a scenario where the network changes over time, and we consider two different set of assumptions on the network instantiations: they are i.i.d. and connected on the average, or they are deterministic and strongly connected for every finite time window of a fixed size. Simulations show our algorithm is competitive with consensus+innovations type algorithms, achieving a smaller MSE at each iteration in all considered scenarios.
△ Less
Submitted 4 July, 2018;
originally announced July 2018.
-
Topologically Anosov plane homeomorphisms
Authors:
Gonzalo Cousillas,
Jorge Groisman,
Juliana Xavier
Abstract:
This paper deals with classifying the dynamics of {\it Topologically Anosov} plane homeomorphisms. We prove that a Topologically Anosov homeomorphism $f:\mathbb{R}^2 \to \mathbb{R}^2$ is conjugate to a homothety if it is the time one map of a flow. We also obtain results for the cases when the nonwandering set of $f$ reduces to a fixed point, or if there exists an open, connected, simply connected…
▽ More
This paper deals with classifying the dynamics of {\it Topologically Anosov} plane homeomorphisms. We prove that a Topologically Anosov homeomorphism $f:\mathbb{R}^2 \to \mathbb{R}^2$ is conjugate to a homothety if it is the time one map of a flow. We also obtain results for the cases when the nonwandering set of $f$ reduces to a fixed point, or if there exists an open, connected, simply connected proper subset $U$ such that $U \subset \mathrm{Int}(\overline {f(U)})$, and such that $ \cup_{n\geq 0} f^n (U)= \mathbb{R}^2$. In the general case, we prove a structure theorem for the $α$-limits of orbits with empty $ω$-limit (or the $ω$-limits of orbits with empty $α$-limit), and we show that any basin of attraction (or repulsion) must be unbounded.
△ Less
Submitted 7 May, 2018;
originally announced May 2018.
-
Equipartition of the Entanglement Entropy
Authors:
J. C. Xavier,
F. C. Alcaraz,
G. Sierra
Abstract:
The entanglement in a quantum system that possess an internal symmetry, characterized by the Sz-magnetization or U(1)-charge, is distributed among different sectors. The aim of this letter is to gain a deeper understanding of the contribution to the entanglement entropy in each of those sectors for the ground state of conformal invariant critical one dimensional systems. We find surprisingly that…
▽ More
The entanglement in a quantum system that possess an internal symmetry, characterized by the Sz-magnetization or U(1)-charge, is distributed among different sectors. The aim of this letter is to gain a deeper understanding of the contribution to the entanglement entropy in each of those sectors for the ground state of conformal invariant critical one dimensional systems. We find surprisingly that the entanglement entropy is equally distributed among the different magnetization sectors. Its value is given by the standard area law violating logarithmic term, that depends on the central charge c, minus a double logarithmic correction related to the zero temperature susceptibility. This result provides a new method to estimate simultaneously the central charge c and the critical exponents of U(1)-symmetric quantum chains. The method is numerically simple and gives precise results for the spin-1/2 quantum XXZ chain. We also compute the probability distribution of the magnetization in contiguous sublattices.
△ Less
Submitted 17 July, 2018; v1 submitted 17 April, 2018;
originally announced April 2018.
-
DJAM: distributed Jacobi asynchronous method for learning personal models
Authors:
Inês Almeida,
João Xavier
Abstract:
Processing data collected by a network of agents often boils down to solving an optimization problem. The distributed nature of these problems calls for methods that are, themselves, distributed. While most collaborative learning problems require agents to reach a common (or consensus) model, there are situations in which the consensus solution may not be optimal. For instance, agents may want to…
▽ More
Processing data collected by a network of agents often boils down to solving an optimization problem. The distributed nature of these problems calls for methods that are, themselves, distributed. While most collaborative learning problems require agents to reach a common (or consensus) model, there are situations in which the consensus solution may not be optimal. For instance, agents may want to reach a compromise between agreeing with their neighbors and minimizing a personal loss function. We present DJAM, a Jacobi-like distributed algorithm for learning personalized models. This method is implementation-friendly: it has no hyperparameters that need tuning, it is asynchronous, and its updates only require single-neighbor interactions. We prove that DJAM converges with probability one to the solution, provided that the personal loss functions are strongly convex and have Lipschitz gradient. We then give evidence that DJAM is on par with state-of-the-art methods: our method reaches a solution with error similar to the error of a carefully tuned ADMM in about the same number of single-neighbor interactions.
△ Less
Submitted 20 July, 2018; v1 submitted 26 March, 2018;
originally announced March 2018.
-
Evolution at the edge of expanding populations
Authors:
Maxime Deforet,
Carlos Carmona-Fontaine,
Kirill S. Korolev,
Joao B. Xavier
Abstract:
Predicting evolution of expanding populations is critical to control biological threats such as invasive species and cancer metastasis. Expansion is primarily driven by reproduction and dispersal, but nature abounds with examples of evolution where organisms pay a reproductive cost to disperse faster. When does selection favor this 'survival of the fastest?' We searched for a simple rule, motivate…
▽ More
Predicting evolution of expanding populations is critical to control biological threats such as invasive species and cancer metastasis. Expansion is primarily driven by reproduction and dispersal, but nature abounds with examples of evolution where organisms pay a reproductive cost to disperse faster. When does selection favor this 'survival of the fastest?' We searched for a simple rule, motivated by evolution experiments where swarming bacteria evolved into an hyperswarmer mutant which disperses $ \sim 100\%$ faster but pays a growth cost of $\sim 10 \%$ to make many copies of its flagellum. We analyzed a two-species model based on the Fisher equation to explain this observation: the population expansion rate ($v$) results from an interplay of growth ($r$) and dispersal ($D$) and is independent of the carrying capacity: $v=2\sqrt{rD}$. A mutant can take over the edge only if its expansion rate ($v_2$) exceeds the expansion rate of the established species ($v_1$); this simple condition ($v_2 > v_1$) determines the maximum cost in slower growth that a faster mutant can pay and still be able to take over. Numerical simulations and time-course experiments where we tracked evolution by imaging bacteria suggest that our findings are general: less favorable conditions delay but do not entirely prevent the success of the fastest. Thus, the expansion rate defines a traveling wave fitness, which could be combined with trade-offs to predict evolution of expanding populations.
△ Less
Submitted 25 April, 2019; v1 submitted 21 November, 2017;
originally announced November 2017.
-
On the growth rate inequality for periodic points in the two sphere
Authors:
G. Honorato,
J. Iglesias,
A. Portela,
A. Rovella,
F. Valenzuela,
J. Xavier
Abstract:
Let $f:S^2\to S^2$ be a continuous map such that $deg f = d, |d|>1$. Suppose $f$ has two attracting fixed points denoted $N$ and $S$ and let $A=S^2\setminus \{N,S\}$. Assume that if a loop $γ\subset f^{-1}(A)$ is homotopically trivial in $A$, then $f(γ)$ is also homotopically trivial in $A$. Then, for all $n$, $f$ has at least $|d^n -1|$ fixed points.
Let $f:S^2\to S^2$ be a continuous map such that $deg f = d, |d|>1$. Suppose $f$ has two attracting fixed points denoted $N$ and $S$ and let $A=S^2\setminus \{N,S\}$. Assume that if a loop $γ\subset f^{-1}(A)$ is homotopically trivial in $A$, then $f(γ)$ is also homotopically trivial in $A$. Then, for all $n$, $f$ has at least $|d^n -1|$ fixed points.
△ Less
Submitted 19 July, 2017;
originally announced July 2017.
-
Planar Baumslag-Solitar actions
Authors:
Juan Alonso,
Cristobal Rivas,
Juliana Xavier
Abstract:
Following previous work, where representations of $BS(1,n)$ by planar orientation preserving homeomorphisms with linear diagonalizable conjugating element was studied, we consider the elliptic and parabolic cases. As an application, we prove that there are no faithful representations of $BS(1,n)$ by toral homeomorphisms with conjugating element the Dehn twist map.
Following previous work, where representations of $BS(1,n)$ by planar orientation preserving homeomorphisms with linear diagonalizable conjugating element was studied, we consider the elliptic and parabolic cases. As an application, we prove that there are no faithful representations of $BS(1,n)$ by toral homeomorphisms with conjugating element the Dehn twist map.
△ Less
Submitted 17 March, 2022; v1 submitted 18 July, 2017;
originally announced July 2017.
-
One some planar Baumslag-Solitar actions
Authors:
Juan Alonso,
Nancy Guelman,
Cristobal Rivas,
Juliana Xavier
Abstract:
Let $BS(1,n)= \langle a,b : a b a ^{-1} = b ^n\rangle$ be the solvable Baumslag-Solitar group for $n \geq 2$. We study representations of $BS(1, n)$ on the plane by orientation preserving homeomorphisms, assuming that $a$ acts as a linear map and $b$ as a map with bounded displacement. We find that the possibilities for a faithful action depend greatly on the Jordan canonical form of the map $h$ d…
▽ More
Let $BS(1,n)= \langle a,b : a b a ^{-1} = b ^n\rangle$ be the solvable Baumslag-Solitar group for $n \geq 2$. We study representations of $BS(1, n)$ on the plane by orientation preserving homeomorphisms, assuming that $a$ acts as a linear map and $b$ as a map with bounded displacement. We find that the possibilities for a faithful action depend greatly on the Jordan canonical form of the map $h$ defined by the action of $a$. In case $h$ is diagonalizable over $\mathbb R$, we shall give examples or prove rigidity theorems depending on the eigenvalues. We also show some rigidity in the cases where $h$ is elliptic or parabolic. Then we give applications to the actions of $BS(1, n)$ by homeomorphisms of the torus.
△ Less
Submitted 17 March, 2022; v1 submitted 27 March, 2017;
originally announced March 2017.
-
Preorder Construct on Simple Undirected Graphs
Authors:
Augusto Almeida Santos,
José M. F. Moura,
João Xavier
Abstract:
We construct a novel preorder on the set of nodes of a simple undirected graph. We prove that the preorder (induced by the topology of the graph) is preserved, e.g., by the logistic dynamical system (both in discrete and continuous time). Moreover, the underlying equivalence relation of the preorder corresponds to the coarsest equitable partition (CEP). This will further imply that the logistic dy…
▽ More
We construct a novel preorder on the set of nodes of a simple undirected graph. We prove that the preorder (induced by the topology of the graph) is preserved, e.g., by the logistic dynamical system (both in discrete and continuous time). Moreover, the underlying equivalence relation of the preorder corresponds to the coarsest equitable partition (CEP). This will further imply that the logistic dynamical system on a graph preserves its coarsest equitable partition. The results provide a nontrivial invariant set for the logistic and the like dynamical systems, as we show. We note that our construct provides a functional characterization for the CEP as an alternative to the pure set theoretical iterated degree sequences characterization. The construct and results presented might have independent interest for analysis on graphs or qualitative analysis of dynamical systems over networks.
△ Less
Submitted 9 March, 2017;
originally announced March 2017.
-
Thermodynamic Limit of Interacting Particle Systems over Time-varying Sparse Random Networks
Authors:
Augusto Almeida Santos,
Soummya Kar,
José M. F. Moura,
João Xavier
Abstract:
We establish a functional weak law of large numbers for observable macroscopic state variables of interacting particle systems (e.g., voter and contact processes) over fast time-varying sparse random networks of interactions. We show that, as the number of agents $N$ grows large, the proportion of agents $\left(\overline{Y}_{k}^{N}(t)\right)$ at a certain state $k$ converges in distribution -- or,…
▽ More
We establish a functional weak law of large numbers for observable macroscopic state variables of interacting particle systems (e.g., voter and contact processes) over fast time-varying sparse random networks of interactions. We show that, as the number of agents $N$ grows large, the proportion of agents $\left(\overline{Y}_{k}^{N}(t)\right)$ at a certain state $k$ converges in distribution -- or, more precisely, weakly with respect to the uniform topology on the space of \emph{càdlàg} sample paths -- to the solution of an ordinary differential equation over any compact interval $\left[0,T\right]$. Although the limiting process is Markov, the prelimit processes, i.e., the normalized macrostate vector processes $\left(\mathbf{\overline{Y}}^{N}(t)\right)=\left(\overline{Y}_{1}^{N}(t),\ldots,\overline{Y}_{K}^{N}(t)\right)$, are non-Markov as they are tied to the \emph{high-dimensional} microscopic state of the system, which precludes the direct application of standard arguments for establishing weak convergence. The techniques developed in the paper for establishing weak convergence might be of independent interest.
△ Less
Submitted 26 February, 2017;
originally announced February 2017.
-
Sphere branched coverings and the growth rate inequality
Authors:
J. Iglesias,
A. Portela,
A. Rovella,
J. Xavier
Abstract:
We show that the growth inequality rate $$\limsup \frac{1}{n} \log (\# Fix (f^n))\geq \log d$$ holds for branched coverings of degree $d$ of the sphere $S^2$ having a completely invariant simply connected region $R$ with locally connected boundary, except in some degenerate cases with known couterexamples.
We show that the growth inequality rate $$\limsup \frac{1}{n} \log (\# Fix (f^n))\geq \log d$$ holds for branched coverings of degree $d$ of the sphere $S^2$ having a completely invariant simply connected region $R$ with locally connected boundary, except in some degenerate cases with known couterexamples.
△ Less
Submitted 7 December, 2016;
originally announced December 2016.
-
Robustness Properties in Fictitious-Play-Type Algorithms
Authors:
Brian Swenson,
Soummya Kar,
João Xavier,
David S. Leslie
Abstract:
Fictitious play (FP) is a canonical game-theoretic learning algorithm which has been deployed extensively in decentralized control scenarios. However standard treatments of FP, and of many other game-theoretic models, assume rather idealistic conditions which rarely hold in realistic control scenarios. This paper considers a broad class of best response learning algorithms, that we refer to as FP-…
▽ More
Fictitious play (FP) is a canonical game-theoretic learning algorithm which has been deployed extensively in decentralized control scenarios. However standard treatments of FP, and of many other game-theoretic models, assume rather idealistic conditions which rarely hold in realistic control scenarios. This paper considers a broad class of best response learning algorithms, that we refer to as FP-type algorithms. In such an algorithm, given some (possibly limited) information about the history of actions, each individual forecasts the future play and chooses a (myopic) best action given their forecast. We provide a unifed analysis of the behavior of FP-type algorithms under an important class of perturbations, thus demonstrating robustness to deviations from the idealistic operating conditions that have been previously assumed. This robustness result is then used to derive convergence results for two control-relevant relaxations of standard game-theoretic applications: distributed (network-based) implementation without full observability and asynchronous deployment (including in continuous time). In each case the results follow as a direct consequence of the main robustness result.
△ Less
Submitted 28 September, 2016;
originally announced September 2016.
-
Deterministic aperiodic composite lattice-structured silicon thin films for photon management
Authors:
Jolly Xavier,
Juergen Probst,
Christiane Becker
Abstract:
Exotic manipulation of the flow of photons in nanoengineered semiconductor materials with an aperiodic distribution of nanostructures plays a key role in efficiency-enhanced and industrially viable broadband photonic technologies. Through a generic deterministic nanotechnological route, in addition to periodic, transversely quasicrystallographic or disordered random photonic lattices, here we show…
▽ More
Exotic manipulation of the flow of photons in nanoengineered semiconductor materials with an aperiodic distribution of nanostructures plays a key role in efficiency-enhanced and industrially viable broadband photonic technologies. Through a generic deterministic nanotechnological route, in addition to periodic, transversely quasicrystallographic or disordered random photonic lattices, here we show scalable nanostructured semiconductor thin films on large area nanoimprinted substrates up to 4cm^2 with advanced functional features of aperiodic composite nanophotonic lattices having tailorable supercell tiles. The richer Fourier spectra of the presented artificially nanostructured materials with well-defined lattice point morphologies are designed functionally akin to two-dimensional incommensurate intergrowth aperiodic lattices-comprising periodic photonic crystals and in-plane quasicrystals as subgroups. The composite photonic lattice-structured crystalline silicon thin films with tapered nanoholes or nanocone-nanoholes are presented showing up to +30 % achievable short circuit current density in comparison to a periodic counterpart where as it is up to +190 % in comparison to a reference unstructured silicon thin film of 300 nanometers of thickness. In view of scalable bottom-up integrated device fabrication processes, the structural analysis is further extended to liquid phase crystallized double-side-textured deterministic aperiodic lattice-structured 10 micrometers thick large area crystalline silicon film.
△ Less
Submitted 19 May, 2016;
originally announced May 2016.