Showing 1–50 of 112 results for author: Campos, C

Measurement of $C\!P$ violation observables in $D^+\rightarrow K^-K^+π^+$ decays

Authors: LHCb collaboration, R. Aaij, A. S. W. Abdelmotteleb, C. Abellan Beteta, F. Abudinén, T. Ackernley, A. A. Adefisoye, B. Adeva, M. Adinolfi, P. Adlarson, C. Agapopoulou, C. A. Aidala, Z. Ajaltouni, S. Akar, K. Akiba, P. Albicocco, J. Albrecht, F. Alessio, M. Alexander, Z. Aliouche, P. Alvarez Cartelle, R. Amalric, S. Amato, J. L. Amey, Y. Amhis , et al. (1109 additional authors not shown)

Abstract: A search for violation of the charge-parity $C\!P$ symmetry in the $D^+\rightarrow K^-K^+π^+$ decay is presented, with proton-proton collision data corresponding to an integrated luminosity of 5.4 fb$^{-1}$, collected at a center-of-mass energy of $13$ TeV with the LHCb detector. A novel model-independent technique is used to compare the $D^+$ and $D^-$ phase-space distributions, with instrumental… ▽ More A search for violation of the charge-parity $C\!P$ symmetry in the $D^+\rightarrow K^-K^+π^+$ decay is presented, with proton-proton collision data corresponding to an integrated luminosity of 5.4 fb$^{-1}$, collected at a center-of-mass energy of $13$ TeV with the LHCb detector. A novel model-independent technique is used to compare the $D^+$ and $D^-$ phase-space distributions, with instrumental asymmetries subtracted using the $D^+_{s}\rightarrow K^-K^+π^+$ decay as a control channel. The $p$-value for the hypothesis of $C\!P$ conservation is $8.1\%$. The $C\!P$ asymmetry observables $A_{C\!P|S}^{φπ^+} = (0.95 \pm 0.43_{stat} \pm 0.26_{syst})\times 10^{-3}$ and $A_{C\!P|S}^{\overline{K}^{*0}K^+} = (-0.26 \pm 0.56_{ stat} \pm 0.18_{syst})\times 10^{-3}$ are also measured. These results show no evidence of $C\!P$ violation and represent the most sensitive search performed through the phase space of a multibody decay. △ Less

Submitted 2 September, 2024; originally announced September 2024.

Comments: All figures and tables, along with machine-readable versions and any supplementary material and additional information, are available at https://lbfence.cern.ch/alcm/public/analysis/full-details/1616 (LHCb public pages)

Report number: CERN-EP-2024-204, LHCb-PAPER-2024-019

A New Control Law for TS Fuzzy Models: Less Conservative LMI Conditions by Using Membership Functions Derivative

Authors: Leonardo Amaral Mozelli, Victor Costa da Silva Campos

Abstract: This note proposes a new type of Parallel Distributed Controller (PDC) for Takagi-Sugeno (TS) fuzzy models. Our idea consists of using two control terms based on state feedback, one composed of a convex combination of linear gains weighted by the normalized membership grade, as in traditional PDC, and the other composed of linear gains weighted by the time-derivatives of the membership functions.… ▽ More This note proposes a new type of Parallel Distributed Controller (PDC) for Takagi-Sugeno (TS) fuzzy models. Our idea consists of using two control terms based on state feedback, one composed of a convex combination of linear gains weighted by the normalized membership grade, as in traditional PDC, and the other composed of linear gains weighted by the time-derivatives of the membership functions. We present the design conditions as Linear Matrix Inequalities, solvable through numerical optimization tools. Numerical examples are given to illustrate the advantages of the proposed approach, which contains the the traditional PDC as a special case. △ Less

Submitted 15 August, 2024; originally announced August 2024.

Comments: 20 pages, 4 figures

MSC Class: 93C42; 93D15; 37B25

Nerva: a Truly Sparse Implementation of Neural Networks

Authors: Wieger Wesselink, Bram Grooten, Qiao Xiao, Cassio de Campos, Mykola Pechenizkiy

Abstract: We introduce Nerva, a fast neural network library under development in C++. It supports sparsity by using the sparse matrix operations of Intel's Math Kernel Library (MKL), which eliminates the need for binary masks. We show that Nerva significantly decreases training time and memory usage while reaching equivalent accuracy to PyTorch. We run static sparse experiments with an MLP on CIFAR-10. On h… ▽ More We introduce Nerva, a fast neural network library under development in C++. It supports sparsity by using the sparse matrix operations of Intel's Math Kernel Library (MKL), which eliminates the need for binary masks. We show that Nerva significantly decreases training time and memory usage while reaching equivalent accuracy to PyTorch. We run static sparse experiments with an MLP on CIFAR-10. On high sparsity levels like $99\%$, the runtime is reduced by a factor of $4\times$ compared to a PyTorch model using masks. Similar to other popular frameworks such as PyTorch and Keras, Nerva offers a Python interface for users to work with. △ Less

Submitted 24 July, 2024; originally announced July 2024.

Comments: The Nerva library is available at https://github.com/wiegerw/nerva

Scaling Continuous Latent Variable Models as Probabilistic Integral Circuits

Authors: Gennaro Gala, Cassio de Campos, Antonio Vergari, Erik Quaeghebeur

Abstract: Probabilistic integral circuits (PICs) have been recently introduced as probabilistic models enjoying the key ingredient behind expressive generative models: continuous latent variables (LVs). PICs are symbolic computational graphs defining continuous LV models as hierarchies of functions that are summed and multiplied together, or integrated over some LVs. They are tractable if LVs can be analyti… ▽ More Probabilistic integral circuits (PICs) have been recently introduced as probabilistic models enjoying the key ingredient behind expressive generative models: continuous latent variables (LVs). PICs are symbolic computational graphs defining continuous LV models as hierarchies of functions that are summed and multiplied together, or integrated over some LVs. They are tractable if LVs can be analytically integrated out, otherwise they can be approximated by tractable probabilistic circuits (PC) encoding a hierarchical numerical quadrature process, called QPCs. So far, only tree-shaped PICs have been explored, and training them via numerical quadrature requires memory-intensive processing at scale. In this paper, we address these issues, and present: (i) a pipeline for building DAG-shaped PICs out of arbitrary variable decompositions, (ii) a procedure for training PICs using tensorized circuit architectures, and (iii) neural functional sharing techniques to allow scalable training. In extensive experiments, we showcase the effectiveness of functional sharing and the superiority of QPCs over traditional PCs. △ Less

Submitted 10 June, 2024; originally announced June 2024.

Momentum-based gradient descent methods for Lie groups

Authors: Cédric M. Campos, David Martín de Diego, José Torrente

Abstract: Polyak's Heavy Ball (PHB; Polyak, 1964), a.k.a. Classical Momentum, and Nesterov's Accelerated Gradient (NAG; Nesterov, 1983) are well know examples of momentum-descent methods for optimization. While the latter outperforms the former, solely generalizations of PHB-like methods to nonlinear spaces have been described in the literature. We propose here a generalization of NAG-like methods for Lie g… ▽ More Polyak's Heavy Ball (PHB; Polyak, 1964), a.k.a. Classical Momentum, and Nesterov's Accelerated Gradient (NAG; Nesterov, 1983) are well know examples of momentum-descent methods for optimization. While the latter outperforms the former, solely generalizations of PHB-like methods to nonlinear spaces have been described in the literature. We propose here a generalization of NAG-like methods for Lie group optimization based on the variational one-to-one correspondence between classical and accelerated momentum methods (Campos et al., 2023). Numerical experiments are shown. △ Less

Submitted 14 April, 2024; originally announced April 2024.

Comments: 24 pages, 2 algorithms, 5 figures

MSC Class: 65K10 (Primary) 70G45; 22E99 (Secondary)

A Semi-Lagrangian Approach for Time and Energy Path Planning Optimization in Static Flow Fields

Authors: Víctor C. da S. Campos, Armando A. Neto, Douglas G. Macharet

Abstract: Efficient path planning for autonomous mobile robots is a critical problem across numerous domains, where optimizing both time and energy consumption is paramount. This paper introduces a novel methodology that considers the dynamic influence of an environmental flow field and considers geometric constraints, including obstacles and forbidden zones, enriching the complexity of the planning problem… ▽ More Efficient path planning for autonomous mobile robots is a critical problem across numerous domains, where optimizing both time and energy consumption is paramount. This paper introduces a novel methodology that considers the dynamic influence of an environmental flow field and considers geometric constraints, including obstacles and forbidden zones, enriching the complexity of the planning problem. We formulate it as a multi-objective optimal control problem, propose a novel transformation called Harmonic Transformation, and apply a semi-Lagrangian scheme to solve it. The set of Pareto efficient solutions is obtained considering two distinct approaches: a deterministic method and an evolutionary-based one, both of which are designed to make use of the proposed Harmonic Transformation. Through an extensive analysis of these approaches, we demonstrate their efficacy in finding optimized paths. △ Less

Submitted 14 June, 2024; v1 submitted 25 March, 2024; originally announced March 2024.

Comments: 12 pages, initial paper submission; Preprint submitted to Journal of the Franklin Institute

Soft Learning Probabilistic Circuits

Authors: Soroush Ghandi, Benjamin Quost, Cassio de Campos

Abstract: Probabilistic Circuits (PCs) are prominent tractable probabilistic models, allowing for a range of exact inferences. This paper focuses on the main algorithm for training PCs, LearnSPN, a gold standard due to its efficiency, performance, and ease of use, in particular for tabular data. We show that LearnSPN is a greedy likelihood maximizer under mild assumptions. While inferences in PCs may use th… ▽ More Probabilistic Circuits (PCs) are prominent tractable probabilistic models, allowing for a range of exact inferences. This paper focuses on the main algorithm for training PCs, LearnSPN, a gold standard due to its efficiency, performance, and ease of use, in particular for tabular data. We show that LearnSPN is a greedy likelihood maximizer under mild assumptions. While inferences in PCs may use the entire circuit structure for processing queries, LearnSPN applies a hard method for learning them, propagating at each sum node a data point through one and only one of the children/edges as in a hard clustering process. We propose a new learning procedure named SoftLearn, that induces a PC using a soft clustering process. We investigate the effect of this learning-inference compatibility in PCs. Our experiments show that SoftLearn outperforms LearnSPN in many situations, yielding better likelihoods and arguably better samples. We also analyze comparable tractable models to highlight the differences between soft/hard learning and model querying. △ Less

Submitted 21 March, 2024; originally announced March 2024.

Probabilistic Circuits with Constraints via Convex Optimization

Authors: Soroush Ghandi, Benjamin Quost, Cassio de Campos

Abstract: This work addresses integrating probabilistic propositional logic constraints into the distribution encoded by a probabilistic circuit (PC). PCs are a class of tractable models that allow efficient computations (such as conditional and marginal probabilities) while achieving state-of-the-art performance in some domains. The proposed approach takes both a PC and constraints as inputs, and outputs a… ▽ More This work addresses integrating probabilistic propositional logic constraints into the distribution encoded by a probabilistic circuit (PC). PCs are a class of tractable models that allow efficient computations (such as conditional and marginal probabilities) while achieving state-of-the-art performance in some domains. The proposed approach takes both a PC and constraints as inputs, and outputs a new PC that satisfies the constraints. This is done efficiently via convex optimization without the need to retrain the entire model. Empirical evaluations indicate that the combination of constraints and PCs can have multiple use cases, including the improvement of model performance under scarce or incomplete data, as well as the enforcement of machine learning fairness measures into the model without compromising model fitness. We believe that these ideas will open possibilities for multiple other applications involving the combination of logics and deep probabilistic models. △ Less

Submitted 19 March, 2024; originally announced March 2024.

RR Lyrae Stars Belonging to the Candidate Globular Cluster Patchick 99

Authors: Evan Butler, Andrea Kunder, Zdenek Prudil, Kevin R. Covey, Macy Ball, Carlos Campos, Kaylen Gollnick, Julio Olivares Carvajal, Joanne Hughes, Kathryn Devine, Christian I. Johnson, A. Katherina Vivas, Michael R. Rich, Meridith Joyce, Iulia T. Simon, Tommaso Marchetti, Andreas J. Koch-Hansen, William I. Clarkson, Rebekah Kuss

Abstract: Patchick 99 is a candidate globular cluster located in the direction of the Galactic bulge, with a proper motion almost identical to the field and extreme field star contamination. A recent analysis suggests it is a low-luminosity globular cluster with a population of RR Lyrae stars. We present new spectra of stars in and around Patchick 99, targeting specifically the 3 RR Lyrae stars associated w… ▽ More Patchick 99 is a candidate globular cluster located in the direction of the Galactic bulge, with a proper motion almost identical to the field and extreme field star contamination. A recent analysis suggests it is a low-luminosity globular cluster with a population of RR Lyrae stars. We present new spectra of stars in and around Patchick 99, targeting specifically the 3 RR Lyrae stars associated with the cluster as well as the other RR Lyrae stars in the field. A sample of 53 giant stars selected from proper motions and a position on CMD are also observed. The three RR Lyrae stars associated with the cluster have similar radial velocities and distances, and two of the targeted giants also have radial velocities in this velocity regime and [Fe/H] metallicities that are slightly more metal-poor than the field. Therefore, if Patchick 99 is a bonafide globular cluster, it would have a radial velocity of -92+/-10 km s-1, a distance of 6.7+/-0.4 kpc (as determined from the RR Lyrae stars), and an orbit that confines it to the inner bulge. △ Less

Submitted 25 January, 2024; v1 submitted 24 January, 2024; originally announced January 2024.

Comments: Accepted to The Astrophysical Journal Letters. Replaced due to a typo in the title

Gemini: A Family of Highly Capable Multimodal Models

Authors: Gemini Team, Rohan Anil, Sebastian Borgeaud, Jean-Baptiste Alayrac, Jiahui Yu, Radu Soricut, Johan Schalkwyk, Andrew M. Dai, Anja Hauth, Katie Millican, David Silver, Melvin Johnson, Ioannis Antonoglou, Julian Schrittwieser, Amelia Glaese, Jilin Chen, Emily Pitler, Timothy Lillicrap, Angeliki Lazaridou, Orhan Firat, James Molloy, Michael Isard, Paul R. Barham, Tom Hennigan, Benjamin Lee , et al. (1325 additional authors not shown)

Abstract: This report introduces a new family of multimodal models, Gemini, that exhibit remarkable capabilities across image, audio, video, and text understanding. The Gemini family consists of Ultra, Pro, and Nano sizes, suitable for applications ranging from complex reasoning tasks to on-device memory-constrained use-cases. Evaluation on a broad range of benchmarks shows that our most-capable Gemini Ultr… ▽ More This report introduces a new family of multimodal models, Gemini, that exhibit remarkable capabilities across image, audio, video, and text understanding. The Gemini family consists of Ultra, Pro, and Nano sizes, suitable for applications ranging from complex reasoning tasks to on-device memory-constrained use-cases. Evaluation on a broad range of benchmarks shows that our most-capable Gemini Ultra model advances the state of the art in 30 of 32 of these benchmarks - notably being the first model to achieve human-expert performance on the well-studied exam benchmark MMLU, and improving the state of the art in every one of the 20 multimodal benchmarks we examined. We believe that the new capabilities of the Gemini family in cross-modal reasoning and language understanding will enable a wide variety of use cases. We discuss our approach toward post-training and deploying Gemini models responsibly to users through services including Gemini, Gemini Advanced, Google AI Studio, and Cloud Vertex AI. △ Less

Submitted 17 June, 2024; v1 submitted 18 December, 2023; originally announced December 2023.

The Milky Way Bulge extra-tidal star survey: BH 261 (AL 3)

Authors: Andrea Kunder, Zdenek Prudil, Kevin Covey, Joanne Hughes, Meridith Joyce, Iulia T. Simion, Rebekah Kuss, Carlos Campos, Christian I. Johnson, Catherine A. Pilachowski, Kristen A. Larson, Andreas J. Koch-Hansen, Tommaso Marchetti, Michael R. Rich, Evan Butler, William I. Clarkson, Michael J. Rivet, Kathryn Devine, A. Katherina Vivas, Gabriel I. Perren, Mario Soto, Erika Silva

Abstract: The Milky Way Bulge extra-tidal star survey (MWBest) is a spectroscopic survey with the goal of identifying stripped globular cluster stars from inner Galaxy clusters. In this way, an indication of the fraction of metal-poor bulge stars that originated from globular clusters can be determined. We observed and analyzed stars in and around BH 261, an understudied globular cluster in the bulge. From… ▽ More The Milky Way Bulge extra-tidal star survey (MWBest) is a spectroscopic survey with the goal of identifying stripped globular cluster stars from inner Galaxy clusters. In this way, an indication of the fraction of metal-poor bulge stars that originated from globular clusters can be determined. We observed and analyzed stars in and around BH 261, an understudied globular cluster in the bulge. From seven giants within the tidal radius of the cluster, we measured an average heliocentric radial velocity of <RV> = -61 +- 2.6 km/s with a radial velocity dispersion of σ= 6.1 +- 1.9 km/s. The large velocity dispersion may have arisen from tidal heating in the cluster's orbit about the Galactic center, or because BH 261 has a high dynamical mass as well as a high mass-to-light ratio. From spectra of five giants, we measure an average metallicity of <[Fe/H]> = -1.1 +- 0.2 dex. We also spectroscopically confirm an RR Lyrae star in BH 261, which yields a distance to the cluster of 7.1 +- 0.4~kpc. Stars with 3D velocities and metallicities consistent with BH 261 reaching to ~0.5 degrees from the cluster are identified. A handful of these stars are also consistent with the spatial distribution of that potential debris from models focussing on the most recent disruption of the cluster. △ Less

Submitted 27 October, 2023; originally announced October 2023.

Comments: accepted for publication in The Astronomical Journal

Probabilistic Integral Circuits

Authors: Gennaro Gala, Cassio de Campos, Robert Peharz, Antonio Vergari, Erik Quaeghebeur

Abstract: Continuous latent variables (LVs) are a key ingredient of many generative models, as they allow modelling expressive mixtures with an uncountable number of components. In contrast, probabilistic circuits (PCs) are hierarchical discrete mixtures represented as computational graphs composed of input, sum and product units. Unlike continuous LV models, PCs provide tractable inference but are limited… ▽ More Continuous latent variables (LVs) are a key ingredient of many generative models, as they allow modelling expressive mixtures with an uncountable number of components. In contrast, probabilistic circuits (PCs) are hierarchical discrete mixtures represented as computational graphs composed of input, sum and product units. Unlike continuous LV models, PCs provide tractable inference but are limited to discrete LVs with categorical (i.e. unordered) states. We bridge these model classes by introducing probabilistic integral circuits (PICs), a new language of computational graphs that extends PCs with integral units representing continuous LVs. In the first place, PICs are symbolic computational graphs and are fully tractable in simple cases where analytical integration is possible. In practice, we parameterise PICs with light-weight neural nets delivering an intractable hierarchical continuous mixture that can be approximated arbitrarily well with large PCs using numerical quadrature. On several distribution estimation benchmarks, we show that such PIC-approximating PCs systematically outperform PCs commonly learned via expectation-maximization or SGD. △ Less

Submitted 25 October, 2023; originally announced October 2023.

Modelling pyro-convection phenomenon during a mega-fire event in Portugal

Authors: Cátia Campos, Flavio Tiago Couto, Jean-Baptiste Filippi, Roberta Baggio, Rui Salgado

Abstract: The present study contributes to an increased understanding of pyro-convection phenomena by using a fire-atmosphere coupled simulation, and investigates in detail the large-scale meteorological conditions affecting Portugal during the occurrence of multiple mega-fires events on 15 October 2017. Two numerical simulations were performed using the MesoNH atmospheric model. The first simulation, was r… ▽ More The present study contributes to an increased understanding of pyro-convection phenomena by using a fire-atmosphere coupled simulation, and investigates in detail the large-scale meteorological conditions affecting Portugal during the occurrence of multiple mega-fires events on 15 October 2017. Two numerical simulations were performed using the MesoNH atmospheric model. The first simulation, was run for a large single domain (300 x 250 grid points) with a 15 km resolution. In the second one, the MesoNH was coupled to a fire propagation model (ForeFire) to study in detail the Quiaios's fire. To optimize both high resolution in the proximity of the fire region and computational efficiency, the simulation is set up using 3 nested domains (300 x 300 grid points) with horizontal resolution of 2000 m, 400 m, and 80 m respectively. The emission into the atmosphere of the heat and the water vapour fluxes caused by the evolving fire is managed by the ForeFire code. The fire spatio-temporal evolution is based on an assigned map, which follows what reported by public authorities. At the large scale, the simulation shows the evolution of the hurricane Ophelia, pointing out the influence of south/southwest winds on the rapid spread of active fires, as well as the subtropical moisture transport toward mainland Portugal in the early evening, when violent pyro-convective activity was observed in Central Portugal. The coupled simulation allowed to reproduce the formation of a PyroCu cloud inside the smoke plume. The convective updraughts caused by the fire led to the vertical transport of water vapour to higher levels and enhanced the development of a high-based cloud over a dry atmospheric layer within the smoke plume. △ Less

Submitted 25 October, 2023; originally announced October 2023.

Journal ref: Atmospheric Research, 2023, 290, pp.106776

Federated Self-Supervised Learning of Monocular Depth Estimators for Autonomous Vehicles

Authors: Elton F. de S. Soares, Carlos Alberto V. Campos

Abstract: Image-based depth estimation has gained significant attention in recent research on computer vision for autonomous vehicles in intelligent transportation systems. This focus stems from its cost-effectiveness and wide range of potential applications. Unlike binocular depth estimation methods that require two fixed cameras, monocular depth estimation methods only rely on a single camera, making them… ▽ More Image-based depth estimation has gained significant attention in recent research on computer vision for autonomous vehicles in intelligent transportation systems. This focus stems from its cost-effectiveness and wide range of potential applications. Unlike binocular depth estimation methods that require two fixed cameras, monocular depth estimation methods only rely on a single camera, making them highly versatile. While state-of-the-art approaches for this task leverage self-supervised learning of deep neural networks in conjunction with tasks like pose estimation and semantic segmentation, none of them have explored the combination of federated learning and self-supervision to train models using unlabeled and private data captured by autonomous vehicles. The utilization of federated learning offers notable benefits, including enhanced privacy protection, reduced network consumption, and improved resilience to connectivity issues. To address this gap, we propose FedSCDepth, a novel method that combines federated learning and deep self-supervision to enable the learning of monocular depth estimators with comparable effectiveness and superior efficiency compared to the current state-of-the-art methods. Our evaluation experiments conducted on Eigen's Split of the KITTI dataset demonstrate that our proposed method achieves near state-of-the-art performance, with a test loss below 0.13 and requiring, on average, only 1.5k training steps and up to 0.415 GB of weight data transfer per autonomous vehicle on each round. △ Less

Submitted 7 October, 2023; originally announced October 2023.

Comments: 16 pages, 8 figures, journal preprint

Probabilistic Multi-Dimensional Classification

Authors: Vu-Linh Nguyen, Yang Yang, Cassio de Campos

Abstract: Multi-dimensional classification (MDC) can be employed in a range of applications where one needs to predict multiple class variables for each given instance. Many existing MDC methods suffer from at least one of inaccuracy, scalability, limited use to certain types of data, hardness of interpretation or lack of probabilistic (uncertainty) estimations. This paper is an attempt to address all these… ▽ More Multi-dimensional classification (MDC) can be employed in a range of applications where one needs to predict multiple class variables for each given instance. Many existing MDC methods suffer from at least one of inaccuracy, scalability, limited use to certain types of data, hardness of interpretation or lack of probabilistic (uncertainty) estimations. This paper is an attempt to address all these disadvantages simultaneously. We propose a formal framework for probabilistic MDC in which learning an optimal multi-dimensional classifier can be decomposed, without loss of generality, into learning a set of (smaller) single-variable multi-class probabilistic classifiers and a directed acyclic graph. Current and future developments of both probabilistic classification and graphical model learning can directly enhance our framework, which is flexible and provably optimal. A collection of experiments is conducted to highlight the usefulness of this MDC framework. △ Less

Submitted 25 November, 2023; v1 submitted 10 June, 2023; originally announced June 2023.

Comments: Accepted for the 39th Conference on Uncertainty in Artificial Intelligence (UAI 2023)

Post-2000 Nonlinear Optical Materials and Measurements: Data Tables and Best Practices

Authors: Nathalie Vermeulen, Daniel Espinosa, Adam Ball, John Ballato, Philippe Boucaud, Georges Boudebs, Cecilia L. A. V. Campos, Peter Dragic, Anderson S. L. Gomes, Mikko J. Huttunen, Nathaniel Kinsey, Rich Mildren, Dragomir Neshev, Lazaro Padilha, Minhao Pu, Ray Secondo, Eiji Tokunaga, Dmitry Turchinovich, Jingshi Yan, Kresten Yvind, Ksenia Dolgaleva, Eric W. Van Stryland

Abstract: In its 60 years of existence, the field of nonlinear optics has gained momentum especially over the past two decades thanks to major breakthroughs in material science and technology. In this article, we present a new set of data tables listing nonlinear-optical properties for different material categories as reported in the literature since 2000. The papers included in the data tables are represen… ▽ More In its 60 years of existence, the field of nonlinear optics has gained momentum especially over the past two decades thanks to major breakthroughs in material science and technology. In this article, we present a new set of data tables listing nonlinear-optical properties for different material categories as reported in the literature since 2000. The papers included in the data tables are representative experimental works on bulk materials, solvents, 0D-1D-2D materials, metamaterials, fiber waveguiding materials, on-chip waveguiding materials, hybrid waveguiding systems, and materials suitable for nonlinear optics at THz frequencies. In addition to the data tables, we also provide best practices for performing and reporting nonlinear-optical experiments. These best practices underpin the selection process that was used for including papers in the tables. While the tables indeed show strong advancements in the field over the past two decades, we encourage the nonlinear-optics community to implement the identified best practices in future works. This will allow a more adequate comparison, interpretation and use of the published parameters, and as such further stimulate the overall progress in nonlinear-optical science and applications. △ Less

Submitted 21 May, 2023; v1 submitted 15 January, 2023; originally announced January 2023.

Comments: 209 pages, 8 figures, 14 tables. Journal of Physics:Photonics (2023)

Journal ref: J. Phys. Photonics, 5, 035001 (2023)

Continuous Mixtures of Tractable Probabilistic Models

Authors: Alvaro H. C. Correia, Gennaro Gala, Erik Quaeghebeur, Cassio de Campos, Robert Peharz

Abstract: Probabilistic models based on continuous latent spaces, such as variational autoencoders, can be understood as uncountable mixture models where components depend continuously on the latent code. They have proven to be expressive tools for generative and probabilistic modelling, but are at odds with tractable probabilistic inference, that is, computing marginals and conditionals of the represented… ▽ More Probabilistic models based on continuous latent spaces, such as variational autoencoders, can be understood as uncountable mixture models where components depend continuously on the latent code. They have proven to be expressive tools for generative and probabilistic modelling, but are at odds with tractable probabilistic inference, that is, computing marginals and conditionals of the represented probability distribution. Meanwhile, tractable probabilistic models such as probabilistic circuits (PCs) can be understood as hierarchical discrete mixture models, and thus are capable of performing exact inference efficiently but often show subpar performance in comparison to continuous latent-space models. In this paper, we investigate a hybrid approach, namely continuous mixtures of tractable models with a small latent dimension. While these models are analytically intractable, they are well amenable to numerical integration schemes based on a finite set of integration points. With a large enough number of integration points the approximation becomes de-facto exact. Moreover, for a finite set of integration points, the integration method effectively compiles the continuous mixture into a standard PC. In experiments, we show that this simple scheme proves remarkably effective, as PCs learnt this way set new state of the art for tractable models on many standard density estimation benchmarks. △ Less

Submitted 24 March, 2023; v1 submitted 21 September, 2022; originally announced September 2022.

Analyzing and Calibrating Risk Assessment by Software Developers

Authors: Yukasa Murakami, Masateru Tsunoda, Eduardo C. Campos

Abstract: In software project management, risk management is a critical factor. Project managers use existing lists of risk or perform brainstorming to identify the risks. However, it is not easy to perceive all the risks objectively. As a result, some risks are perceived based on subjective impression, which leads to risk biases. So, our goals are (i) We clarify the risk perception of developers to enhance… ▽ More In software project management, risk management is a critical factor. Project managers use existing lists of risk or perform brainstorming to identify the risks. However, it is not easy to perceive all the risks objectively. As a result, some risks are perceived based on subjective impression, which leads to risk biases. So, our goals are (i) We clarify the risk perception of developers to enhance the reliability of the brainstorming, and (ii) we calibrate the risk assessment based on a mathematical model to make more accurate risk list. In the analysis, we collected data concerning the risk perception of 69 professional software developers via a questionnaire. The average number of years of experience among these professionals was 18.3. Using the dataset, we applied factor analysis to clarify the factors that affect the evaluation of risk impact. The questionnaire was based on the risk perception theory established by Slovic, in which "dread" and "unknown" are the major factor of risk perception. The analysis result shows that (i) risk experience (i.e., whether a developer actually faced the risk or not) sometimes affects risk assessment (evaluation of risk impact), (ii) risk perception is considered to be based on unknown and dread factors, and (iii) risk assessment can be calibrated by a mathematical model (the average absolute error was 0.20). △ Less

Submitted 7 September, 2022; originally announced September 2022.

Comments: Japanese letter version is available at: https://www.jstage.jst.go.jp/article/jssst/35/4/35_37/_article/-char/en

Hypothesis tests for multiple responses regression: effect of probiotics on addiction and binge eating disorder

Authors: Lineu Alberto Cavazani de Freitas, Ligia de Oliveira Carlos, Antônio Carlos Ligocki Campos, Wagner Hugo Bonat

Abstract: Clinical trials are common in medical research where multiple non-Gaussian responses and time-dependent observations are frequent. The analysis of data from these studies requires statistical modeling techniques that take these characteristics into account. We propose a general strategy based on the Wald statistics to perform hypothesis tests like ANOVAs, MANOVAs and multiple comparison tests on r… ▽ More Clinical trials are common in medical research where multiple non-Gaussian responses and time-dependent observations are frequent. The analysis of data from these studies requires statistical modeling techniques that take these characteristics into account. We propose a general strategy based on the Wald statistics to perform hypothesis tests like ANOVAs, MANOVAs and multiple comparison tests on regression and dispersion parameters of multivariate covariance generalized linear models (McGLMs). McGLMs provide a general statistical modeling framework for normal and non-normal multivariate data analysis along with a wide range of correlation structures. We design different simulation scenarios to verify the properties of the proposed tests. The results are promising showing that the proposed tests present the levels of confidence close to the specified one for all simulation study scenarios. Complementary to the proposal, we developed implementations in the R language to carry out the tests presented, the codes are available in the supplementary material. The proposal is motivated by the analysis of a clinical trial that aims to evaluate the effect of the use of probiotics in the control of addiction and binge eating disorder in patients undergoing bariatric surgery. The subjects were separated into two groups (placebo and treatment) and evaluated at three different times. The results indicate that addiction and binge eating disorder reduce over time, but there is no difference between groups at each time point. △ Less

Submitted 29 July, 2022; originally announced August 2022.

Observation of well-defined Kohn-anomaly in high-quality graphene devices at room temperature

Authors: Andreij C. Gadelha, Rafael Nadas, Tiago C. Barbosa, Kenji Watanabe, Takashi Taniguchi, Leonardo C. Campos, Markus B. Raschke, Ado Jorio

Abstract: Due to its ultra-thin nature, the study of graphene quantum optoelectronics, like gate-dependent graphene Raman properties, is obscured by interactions with substrates and surroundings. For instance, the use of doped silicon with a capping thermal oxide layer limited the observation to low temperatures of a well-defined Kohn-anomaly behavior, related to the breakdown of the adiabatic Born-Oppenhei… ▽ More Due to its ultra-thin nature, the study of graphene quantum optoelectronics, like gate-dependent graphene Raman properties, is obscured by interactions with substrates and surroundings. For instance, the use of doped silicon with a capping thermal oxide layer limited the observation to low temperatures of a well-defined Kohn-anomaly behavior, related to the breakdown of the adiabatic Born-Oppenheimer approximation. Here, we design an optoelectronic device consisting of single-layer graphene electrically contacted with thin graphite leads, seated on an atomically flat hexagonal boron nitride (hBN) substrate and gated with an ultra-thin gold (Au) layer. We show that this device is optically transparent, has no background optical peaks and photoluminescence from the device components, and no generation of laser-induced electrostatic doping (photodoping). This allows for room-temperature gate-dependent Raman spectroscopy effects that have only been observed at cryogenic temperatures so far, above all the Kohn-anomaly phonon energy normalization. The new device architecture by decoupling graphene optoelectronic properties from the substrate effects, allows for the observation of quantum phenomena at room temperature. △ Less

Submitted 19 July, 2022; originally announced July 2022.

Thick-restarted joint Lanczos bidiagonalization for the GSVD

Authors: Fernando Alvarruiz, Carmen Campos, Jose E. Roman

Abstract: The computation of the partial generalized singular value decomposition (GSVD) of large-scale matrix pairs can be approached by means of iterative methods based on expanding subspaces, particularly Krylov subspaces. We consider the joint Lanczos bidiagonalization method, and analyze the feasibility of adapting the thick restart technique that is being used successfully in the context of other line… ▽ More The computation of the partial generalized singular value decomposition (GSVD) of large-scale matrix pairs can be approached by means of iterative methods based on expanding subspaces, particularly Krylov subspaces. We consider the joint Lanczos bidiagonalization method, and analyze the feasibility of adapting the thick restart technique that is being used successfully in the context of other linear algebra problems. Numerical experiments illustrate the effectiveness of the proposed method. We also compare the new method with an alternative solution via equivalent eigenvalue problems, considering accuracy as well as computational performance. The analysis is done using a parallel implementation in the SLEPc library. △ Less

Submitted 12 May, 2023; v1 submitted 8 June, 2022; originally announced June 2022.

PoVRPoint: Authoring Presentations in Mobile Virtual Reality

Authors: Verena Biener, Travis Gesslein, Daniel Schneider, Felix Kawala, Alexander Otte, Per Ola Kristensson, Michel Pahud, Eyal Ofek, Cuauhtli Campos, Matjaž Kljun, Klen Čopič Pucihar, Jens Grubert

Abstract: Virtual Reality (VR) has the potential to support mobile knowledge workers by complementing traditional input devices with a large three-dimensional output space and spatial input. Previous research on supporting VR knowledge work explored domains such as text entry using physical keyboards and spreadsheet interaction using combined pen and touch input. Inspired by such work, this paper probes the… ▽ More Virtual Reality (VR) has the potential to support mobile knowledge workers by complementing traditional input devices with a large three-dimensional output space and spatial input. Previous research on supporting VR knowledge work explored domains such as text entry using physical keyboards and spreadsheet interaction using combined pen and touch input. Inspired by such work, this paper probes the VR design space for authoring presentations in mobile settings. We propose PoVRPoint -- a set of tools coupling pen- and touch-based editing of presentations on mobile devices, such as tablets, with the interaction capabilities afforded by VR. We study the utility of extended display space to, for example, assist users in identifying target slides, supporting spatial manipulation of objects on a slide, creating animations, and facilitating arrangements of multiple, possibly occluded, shapes. Among other things, our results indicate that 1) the wide field of view afforded by VR results in significantly faster target slide identification times compared to a tablet-only interface for visually salient targets; and 2) the three-dimensional view in VR enables significantly faster object reordering in the presence of occlusion compared to two baseline interfaces. A user study further confirmed that the interaction techniques were found to be usable and enjoyable. △ Less

Submitted 17 January, 2022; originally announced January 2022.

Comments: IEEE VR 2022; to appear in IEEE transactions on visualization and computer graphics, 2022

ACM Class: I.3.7

Journal ref: In IEEE transactions on visualization and computer graphics, 2022

Raman spectra of twisted bilayer graphene close to the magic angle

Authors: Tiago C. Barbosa, Andreij C. Gadelha, Douglas A. A. Ohlberg, Kenji Watanabe, Takashi Taniguchi, Gilberto Medeiros-Ribeiro, Ado Jorio, Leonardo C. Campos

Abstract: In this work, we study the Raman spectra of twisted bilayer graphene samples as a function of their twist-angles ($θ$), ranging from 0.03$^\circ$ to 3.40$^\circ$, where local $θ$ are determined by analysis of their associated moire superlattices, as imaged by scanning microwave impedance microscopy. Three standard excitation laser lines are used (457, 532, and 633 nm wavelengths), and the main Ram… ▽ More In this work, we study the Raman spectra of twisted bilayer graphene samples as a function of their twist-angles ($θ$), ranging from 0.03$^\circ$ to 3.40$^\circ$, where local $θ$ are determined by analysis of their associated moire superlattices, as imaged by scanning microwave impedance microscopy. Three standard excitation laser lines are used (457, 532, and 633 nm wavelengths), and the main Raman active graphene bands (G and 2D) are considered. Our results reveal that electron-phonon interaction influences the G band's linewidth close to the magic angle regardless of laser excitation wavelength. Also, the 2D band lineshape in the $θ$ < 1$^\circ$ regime is dictated by crystal lattice and depends on both the Bernal (AB and BA) stacking bilayer graphene and strain soliton regions (SP). We propose a geometrical model to explain the 2D lineshape variations, and from it, we estimate the SP width when moving towards the magic angle. △ Less

Submitted 5 January, 2022; originally announced January 2022.

Picosecond Timing-Planes for Future Collider Detectors

Authors: P. W. Gorham, C. Campos, C. Hast, K. Jobe, C. Miki, L. Nguyen, M. Olmedo, R. Prechelt, C. Sonoda, G. S. Varner

Abstract: We report experimental test-beam results on dielectric-loaded waveguide detectors that utilize microwave Cherenkov signals to time and characterize high energy particle showers. These results are used to validate models and produce high-fidelity simulations of timing plane systems which yield picosecond time tags and millimeter spatial coordinates for the shower centroid. These timing planes, base… ▽ More We report experimental test-beam results on dielectric-loaded waveguide detectors that utilize microwave Cherenkov signals to time and characterize high energy particle showers. These results are used to validate models and produce high-fidelity simulations of timing plane systems which yield picosecond time tags and millimeter spatial coordinates for the shower centroid. These timing planes, based on the Askaryan effect in solid dielectrics, are most effective at the high center-of-momentum energies planned for the Future Circular Collider (FCC-hh), and are of particular interest in the forward region due to their high radiation immunity. We use our beam test results and GEANT4 simulations to validate a hybrid microwave detector model, which is used to simulate a reference timing plane design for the FCC forward calorimeters. Our results indicate that 0.5-3 ps particle timing is possible for a wide range of collision products in the reference FCC hadron collider detector, even with current technology. △ Less

Submitted 10 February, 2022; v1 submitted 1 December, 2021; originally announced December 2021.

Comments: 34 pages, 35 figures

Internal nonlinear transmission in an Yb mode-locked fiber laser through bifurcations

Authors: Cecília L. A. V. Campos, Lucas B. A. Mélo, Hugo L. D. de S. Cavalcante, Lúcio H. Acioli, Marcio H. G. de Miranda

Abstract: Mode-locked fiber lasers are rich dynamical systems that may present several different types of pulsed operating modes, depending on external control parameters such as pump power. A systematic experimental characterization of such regimes is a challenging problem. Here single pulse regimes of ytterbium mode-locked fiber laser are explored and related to the nonlinear amplitude modulation affectin… ▽ More Mode-locked fiber lasers are rich dynamical systems that may present several different types of pulsed operating modes, depending on external control parameters such as pump power. A systematic experimental characterization of such regimes is a challenging problem. Here single pulse regimes of ytterbium mode-locked fiber laser are explored and related to the nonlinear amplitude modulation affecting the pulses inside the laser cavity. To understand the full dynamics of this system, average power measurements (slow time scale) have been performed and compared to simultaneous time-resolved measurements (fast time scale). The average nonlinear transmission resulting from the amplitude modulation is measured, allowing to relate the transitions between dynamical regimes to the system's nonlinearity. △ Less

Submitted 1 December, 2021; originally announced December 2021.

Comments: 15 pages, 5 figures

Journal ref: Optics Communications, vol. 461 p. 125154 (2020)

Bayesian Modelling of Multivalued Power Curves from an Operational Wind Farm

Authors: L. A. Bull, P. A. Gardner, T. J. Rogers, N. Dervilis, E. J. Cross, E. Papatheou, A. E. Maguire, C. Campos, K. Worden

Abstract: Power curves capture the relationship between wind speed and output power for a specific wind turbine. Accurate regression models of this function prove useful in monitoring, maintenance, design, and planning. In practice, however, the measurements do not always correspond to the ideal curve: power curtailments will appear as (additional) functional components. Such multivalued relationships canno… ▽ More Power curves capture the relationship between wind speed and output power for a specific wind turbine. Accurate regression models of this function prove useful in monitoring, maintenance, design, and planning. In practice, however, the measurements do not always correspond to the ideal curve: power curtailments will appear as (additional) functional components. Such multivalued relationships cannot be modelled by conventional regression, and the associated data are usually removed during pre-processing. The current work suggests an alternative method to infer multivalued relationships in curtailed power data. Using a population-based approach, an overlapping mixture of probabilistic regression models is applied to signals recorded from turbines within an operational wind farm. The model is shown to provide an accurate representation of practical power data across the population. △ Less

Submitted 30 November, 2021; originally announced November 2021.

Journal ref: Mechanical Systems and Signal Processing (2021): 108530

Accuracy Evaluation of Touch Tasks in Commodity Virtual and Augmented Reality Head-Mounted Displays

Authors: Daniel Schneider, Verena Biener, Alexander Otte, Travis Gesslein, Philipp Gagel, Cuauhtli Campos, Klen Čopič Pucihar, Matjaž Kljun, Eyal Ofek, Michel Pahud, Per Ola Kristensson, Jens Grubert

Abstract: An increasing number of consumer-oriented head-mounted displays (HMD) for augmented and virtual reality (AR/VR) are capable of finger and hand tracking. We report on the accuracy of off-the-shelf VR and AR HMDs when used for touch-based tasks such as pointing or drawing. Specifically, we report on the finger tracking accuracy of the VR head-mounted displays Oculus Quest, Vive Pro and the Leap Moti… ▽ More An increasing number of consumer-oriented head-mounted displays (HMD) for augmented and virtual reality (AR/VR) are capable of finger and hand tracking. We report on the accuracy of off-the-shelf VR and AR HMDs when used for touch-based tasks such as pointing or drawing. Specifically, we report on the finger tracking accuracy of the VR head-mounted displays Oculus Quest, Vive Pro and the Leap Motion controller, when attached to a VR HMD, as well as the finger tracking accuracy of the AR head-mounted displays Microsoft HoloLens 2 and Magic Leap. We present the results of two experiments in which we compare the accuracy for absolute and relative pointing tasks using both human participants and a robot. The results suggest that HTC Vive has a lower spatial accuracy than the Oculus Quest and Leap Motion and that the Microsoft HoloLens 2 provides higher spatial accuracy than Magic Leap One. These findings can serve as decision support for researchers and practitioners in choosing which systems to use in the future. △ Less

Submitted 22 September, 2021; originally announced September 2021.

Comments: To appear in SUI 2021, November 09-10, Virtual Conference

ACM Class: I.3.7

Scale-aware direct monocular odometry

Authors: Carlos Campos, Juan D. Tardós

Abstract: We present a generic framework for scale-aware direct monocular odometry based on depth prediction from a deep neural network. In contrast with previous methods where depth information is only partially exploited, we formulate a novel depth prediction residual which allows us to incorporate multi-view depth information. In addition, we propose to use a truncated robust cost function which prevents… ▽ More We present a generic framework for scale-aware direct monocular odometry based on depth prediction from a deep neural network. In contrast with previous methods where depth information is only partially exploited, we formulate a novel depth prediction residual which allows us to incorporate multi-view depth information. In addition, we propose to use a truncated robust cost function which prevents considering inconsistent depth estimations. The photometric and depth-prediction measurements are integrated into a tightly-coupled optimization leading to a scale-aware monocular system which does not accumulate scale drift. Our proposal does not particularize for a concrete neural network, being able to work along with the vast majority of the existing depth prediction solutions. We demonstrate the validity and generality of our proposal evaluating it on the KITTI odometry dataset, using two publicly available neural networks and comparing it with similar approaches and the state-of-the-art for monocular and stereo SLAM. Experiments show that our proposal largely outperforms classic monocular SLAM, being 5 to 9 times more precise, beating similar approaches and having an accuracy which is closer to that of stereo systems. △ Less

Submitted 22 July, 2022; v1 submitted 21 September, 2021; originally announced September 2021.

Comments: This paper has been accepted for publication in the IROS2022 conference

A Semi-Lagrangian Approach for the Minimal Exposure Path Problem in Wireless Sensor Networks

Authors: Armando Alves Neto, Víctor C. da Silva Campos, Douglas G. Macharet

Abstract: A critical metric of the coverage quality in Wireless Sensor Networks (WSNs) is the Minimal Exposure Path (MEP), a path through the environment that least exposes an intruder to the sensor detecting nodes. Many approaches have been proposed in the last decades to solve this optimization problem, ranging from classic (grid-based and Voronoi-based) planners to genetic meta-heuristics. However, most… ▽ More A critical metric of the coverage quality in Wireless Sensor Networks (WSNs) is the Minimal Exposure Path (MEP), a path through the environment that least exposes an intruder to the sensor detecting nodes. Many approaches have been proposed in the last decades to solve this optimization problem, ranging from classic (grid-based and Voronoi-based) planners to genetic meta-heuristics. However, most of them are limited to specific sensing models and obstacle-free spaces. Still, none of them guarantee an optimal solution, and the state-of-the-art is expensive in terms of run-time. Therefore, in this paper, we propose a novel method that models the MEP as an Optimal Control problem and solves it by using a Semi-Lagrangian approach. This framework is shown to converge to the optimal MEP while also incorporates different homogeneous and heterogeneous sensor models and geometric constraints (obstacles). Experiments show that our method dominates the state-of-the-art, improving the results by approximately 10% with a relatively lower execution time. △ Less

Submitted 12 August, 2021; originally announced August 2021.

A Discrete Variational Derivation of Accelerated Methods in Optimization

Authors: Cédric M. Campos, Alejandro Mahillo, David Martín de Diego

Abstract: Many of the new developments in machine learning are connected with gradient-based optimization methods. Recently, these methods have been studied using a variational perspective. This has opened up the possibility of introducing variational and symplectic methods using geometric integration. In particular, in this paper, we introduce variational integrators which allow us to derive different meth… ▽ More Many of the new developments in machine learning are connected with gradient-based optimization methods. Recently, these methods have been studied using a variational perspective. This has opened up the possibility of introducing variational and symplectic methods using geometric integration. In particular, in this paper, we introduce variational integrators which allow us to derive different methods for optimization. Using both, Hamilton's and Lagrange-d'Alembert's principle, we derive two families of respective optimization methods in one-to-one correspondence that generalize Polyak's heavy ball and the well known Nesterov accelerated gradient method, the second of which mimics the behavior of the first reducing the oscillations of classical momentum methods. However, since the systems considered are explicitly time-dependent, the preservation of symplecticity of autonomous systems occurs here solely on the fibers. Several experiments exemplify the result. △ Less

Submitted 20 November, 2022; v1 submitted 4 June, 2021; originally announced June 2021.

Comments: 33 pages, 11 figures

Journal ref: Journal of Machine Learning Research 24 (2023) 25:1-33

Bayesian Kernelised Test of (In)dependence with Mixed-type Variables

Authors: Alessio Benavoli, Cassio de Campos

Abstract: A fundamental task in AI is to assess (in)dependence between mixed-type variables (text, image, sound). We propose a Bayesian kernelised correlation test of (in)dependence using a Dirichlet process model. The new measure of (in)dependence allows us to answer some fundamental questions: Based on data, are (mixed-type) variables independent? How likely is dependence/independence to hold? How high is… ▽ More A fundamental task in AI is to assess (in)dependence between mixed-type variables (text, image, sound). We propose a Bayesian kernelised correlation test of (in)dependence using a Dirichlet process model. The new measure of (in)dependence allows us to answer some fundamental questions: Based on data, are (mixed-type) variables independent? How likely is dependence/independence to hold? How high is the probability that two mixed-type variables are more than just weakly dependent? We theoretically show the properties of the approach, as well as algorithms for fast computation with it. We empirically demonstrate the effectiveness of the proposed method by analysing its performance and by comparing it with other frequentist and Bayesian approaches on a range of datasets and tasks with mixed-type variables. △ Less

Submitted 9 May, 2021; originally announced May 2021.

Blockchain Privacy Through Merge Avoidance and Mixing Services: a Hardness and an Impossibility Result

Authors: Jefferson E. Simoes, Eduardo Ferreira, Daniel S. Menasche, Carlos A. V. Campos

Abstract: Cryptocurrencies typically aim at preserving the privacy of their users. Different cryptocurrencies preserve privacy at various levels, some of them requiring users to rely on strategies to raise the privacy level to their needs. Among those strategies, we focus on two of them: merge avoidance and mixing services. Such strategies may be adopted on top of virtually any blockchain-based cryptocurren… ▽ More Cryptocurrencies typically aim at preserving the privacy of their users. Different cryptocurrencies preserve privacy at various levels, some of them requiring users to rely on strategies to raise the privacy level to their needs. Among those strategies, we focus on two of them: merge avoidance and mixing services. Such strategies may be adopted on top of virtually any blockchain-based cryptocurrency. In this paper, we show that whereas optimal merge avoidance leads to an NP-hard optimization problem, incentive-compatible mixing services are subject to a certain class of impossibility results. Together, our results contribute to the body of work on fundamental limits of privacy mechanisms in blockchain-based cryptocurrencies. △ Less

Submitted 4 January, 2021; originally announced January 2021.

Journal ref: Symposium on Cryptocurrency Analysis (SOCCA) 2020

DynaSLAM II: Tightly-Coupled Multi-Object Tracking and SLAM

Authors: Berta Bescos, Carlos Campos, Juan D. Tardós, José Neira

Abstract: The assumption of scene rigidity is common in visual SLAM algorithms. However, it limits their applicability in populated real-world environments. Furthermore, most scenarios including autonomous driving, multi-robot collaboration and augmented/virtual reality, require explicit motion information of the surroundings to help with decision making and scene understanding. We present in this paper Dyn… ▽ More The assumption of scene rigidity is common in visual SLAM algorithms. However, it limits their applicability in populated real-world environments. Furthermore, most scenarios including autonomous driving, multi-robot collaboration and augmented/virtual reality, require explicit motion information of the surroundings to help with decision making and scene understanding. We present in this paper DynaSLAM II, a visual SLAM system for stereo and RGB-D configurations that tightly integrates the multi-object tracking capability. DynaSLAM II makes use of instance semantic segmentation and of ORB features to track dynamic objects. The structure of the static scene and of the dynamic objects is optimized jointly with the trajectories of both the camera and the moving agents within a novel bundle adjustment proposal. The 3D bounding boxes of the objects are also estimated and loosely optimized within a fixed temporal window. We demonstrate that tracking dynamic objects does not only provide rich clues for scene understanding but is also beneficial for camera tracking. The project code will be released upon acceptance. △ Less

Submitted 15 October, 2020; originally announced October 2020.

ORB-SLAM3: An Accurate Open-Source Library for Visual, Visual-Inertial and Multi-Map SLAM

Authors: Carlos Campos, Richard Elvira, Juan J. Gómez Rodríguez, José M. M. Montiel, Juan D. Tardós

Abstract: This paper presents ORB-SLAM3, the first system able to perform visual, visual-inertial and multi-map SLAM with monocular, stereo and RGB-D cameras, using pin-hole and fisheye lens models. The first main novelty is a feature-based tightly-integrated visual-inertial SLAM system that fully relies on Maximum-a-Posteriori (MAP) estimation, even during the IMU initialization phase. The result is a syst… ▽ More This paper presents ORB-SLAM3, the first system able to perform visual, visual-inertial and multi-map SLAM with monocular, stereo and RGB-D cameras, using pin-hole and fisheye lens models. The first main novelty is a feature-based tightly-integrated visual-inertial SLAM system that fully relies on Maximum-a-Posteriori (MAP) estimation, even during the IMU initialization phase. The result is a system that operates robustly in real-time, in small and large, indoor and outdoor environments, and is 2 to 5 times more accurate than previous approaches. The second main novelty is a multiple map system that relies on a new place recognition method with improved recall. Thanks to it, ORB-SLAM3 is able to survive to long periods of poor visual information: when it gets lost, it starts a new map that will be seamlessly merged with previous maps when revisiting mapped areas. Compared with visual odometry systems that only use information from the last few seconds, ORB-SLAM3 is the first system able to reuse in all the algorithm stages all previous information. This allows to include in bundle adjustment co-visible keyframes, that provide high parallax observations boosting accuracy, even if they are widely separated in time or if they come from a previous mapping session. Our experiments show that, in all sensor configurations, ORB-SLAM3 is as robust as the best systems available in the literature, and significantly more accurate. Notably, our stereo-inertial SLAM achieves an average accuracy of 3.6 cm on the EuRoC drone and 9 mm under quick hand-held motions in the room of TUM-VI dataset, a setting representative of AR/VR scenarios. For the benefit of the community we make public the source code. △ Less

Submitted 23 April, 2021; v1 submitted 23 July, 2020; originally announced July 2020.

Towards Robust Classification with Deep Generative Forests

Authors: Alvaro H. C. Correia, Robert Peharz, Cassio de Campos

Abstract: Decision Trees and Random Forests are among the most widely used machine learning models, and often achieve state-of-the-art performance in tabular, domain-agnostic datasets. Nonetheless, being primarily discriminative models they lack principled methods to manipulate the uncertainty of predictions. In this paper, we exploit Generative Forests (GeFs), a recent class of deep probabilistic models th… ▽ More Decision Trees and Random Forests are among the most widely used machine learning models, and often achieve state-of-the-art performance in tabular, domain-agnostic datasets. Nonetheless, being primarily discriminative models they lack principled methods to manipulate the uncertainty of predictions. In this paper, we exploit Generative Forests (GeFs), a recent class of deep probabilistic models that addresses these issues by extending Random Forests to generative models representing the full joint distribution over the feature space. We demonstrate that GeFs are uncertainty-aware classifiers, capable of measuring the robustness of each prediction as well as detecting out-of-distribution samples. △ Less

Submitted 11 July, 2020; originally announced July 2020.

Comments: Presented at the ICML 2020 Workshop on Uncertainty and Robustness in Deep Learning

The limits of Near Field Immersion Microwave Microscopy evaluated by imaging bilayer graphene Moiré patterns

Authors: Douglas A. A. Ohlberg, Diego Tami, Andreij C. Gadelha, Eliel G. S. Neto, Fabiano C. Santana, Daniel Miranda, Wellington Avelino, Kenji Watanabe, Takashi Taniguchi, Leonardo C. Campos, Jhonattan C. Ramirez, Cássio Gonçalves do Rego, Ado Jorio, Gilberto Medeiros-Ribeiro

Abstract: Molecular and atomic imaging required the development of electron and scanning probe microscopies to surpass the physical limits dictated by diffraction. Nano-infrared experiments and pico-cavity tip-enhanced Raman spectroscopy imaging later demonstrated that radiation in the visible range can surpass this limit by using scanning probe tips to access the near-field regime. Here we show that ultima… ▽ More Molecular and atomic imaging required the development of electron and scanning probe microscopies to surpass the physical limits dictated by diffraction. Nano-infrared experiments and pico-cavity tip-enhanced Raman spectroscopy imaging later demonstrated that radiation in the visible range can surpass this limit by using scanning probe tips to access the near-field regime. Here we show that ultimate resolution can be obtained by using scanning microwave imaging microscopy to reveal structures with feature sizes down to 1~nm using a radiation of 0.1~m in wavelength. As a test material we use twisted bilayer graphene, which is not only a very important recent topic due to the discovery of correlated electron effects such as superconductivity, but also because it provides a sample where we can systematically tune a superstructure Moiré patterns modulation from below one up to tens of nanometers. By analyzing the tip-sample distance dynamics, we demonstrate that this ultimate 10$^8$ probe-to-pattern resolution can be achieved by using liquid immersion microscopy concepts and exquisite force control exerted on nanoscale water menisci. △ Less

Submitted 7 July, 2020; originally announced July 2020.

Comments: suppl. mat included, movies not included and available upon request by email to gilberto@dcc.ufmg.br

Joints in Random Forests

Authors: Alvaro H. C. Correia, Robert Peharz, Cassio de Campos

Abstract: Decision Trees (DTs) and Random Forests (RFs) are powerful discriminative learners and tools of central importance to the everyday machine learning practitioner and data scientist. Due to their discriminative nature, however, they lack principled methods to process inputs with missing features or to detect outliers, which requires pairing them with imputation techniques or a separate generative mo… ▽ More Decision Trees (DTs) and Random Forests (RFs) are powerful discriminative learners and tools of central importance to the everyday machine learning practitioner and data scientist. Due to their discriminative nature, however, they lack principled methods to process inputs with missing features or to detect outliers, which requires pairing them with imputation techniques or a separate generative model. In this paper, we demonstrate that DTs and RFs can naturally be interpreted as generative models, by drawing a connection to Probabilistic Circuits, a prominent class of tractable probabilistic models. This reinterpretation equips them with a full joint distribution over the feature space and leads to Generative Decision Trees (GeDTs) and Generative Forests (GeFs), a family of novel hybrid generative-discriminative models. This family of models retains the overall characteristics of DTs and RFs while additionally being able to handle missing features by means of marginalisation. Under certain assumptions, frequently made for Bayes consistency results, we show that consistency in GeDTs and GeFs extend to any pattern of missing input features, if missing at random. Empirically, we show that our models often outperform common routines to treat missing data, such as K-nearest neighbour imputation, and moreover, that our models can naturally detect outliers by monitoring the marginal probability of input features. △ Less

Submitted 19 November, 2020; v1 submitted 25 June, 2020; originally announced June 2020.

Journal ref: Advances in Neural Information Processing Systems 33 (2020)

Gate-tunable non-volatile photomemory effect in MoS$_2$ transistors

Authors: Andreij C. Gadelha, Alisson R. Cadore, Kenji Watanabe, Takashi Tanigushi, Ana M. de Paula, Leandro M. Malard, Rodrigo G. Lacerda, Leonardo C. Campos

Abstract: Non-volatile memory devices have been limited to flash architectures that are complex devices. Here, we present a unique photomemory effect in MoS$_2$ transistors. The photomemory is based on a photodoping effect - a controlled way of manipulating the density of free charges in monolayer MoS$_2$ using a combination of laser exposure and gate voltage application. The photodoping promotes changes on… ▽ More Non-volatile memory devices have been limited to flash architectures that are complex devices. Here, we present a unique photomemory effect in MoS$_2$ transistors. The photomemory is based on a photodoping effect - a controlled way of manipulating the density of free charges in monolayer MoS$_2$ using a combination of laser exposure and gate voltage application. The photodoping promotes changes on the conductance of MoS$_2$ leading to photomemory states with high memory on/off ratio. Such memory states are non-volatile with an expectation of retaining up to 50 % of the information for tens of years. Furthermore, we show that the photodoping is gate-tunable, enabling control of the recorded memory states. Finally, we propose a model to explain the photodoping, and we provide experimental evidence supporting such a phenomenon. In summary, our work includes the MoS$_2$ phototransistors in the non-volatile memory devices and expands the possibilities of memory application beyond conventional memory architectures. △ Less

Submitted 17 June, 2020; originally announced June 2020.

Journal ref: 2D Mater. 6 025036 (2019)

Local photodoping in monolayer MoS2

Authors: Andreij C. Gadelha, Alisson R. Cadore, Lucas Lafeta, Ana M. de Paula, Leandro M. Malard, Rodrigo G. Lacerda, Leonardo C. Campos

Abstract: Inducing electrostatic doping in 2D materials by laser exposure (photodoping effect) is an exciting route to tune optoelectronic phenomena. However, there is a lack of investigation concerning in what respect the action of photodoping in optoelectronic devices is local. Here, we employ scanning photocurrent microscopy (SPCM) techniques to investigate how a permanent photodoping modulates the photo… ▽ More Inducing electrostatic doping in 2D materials by laser exposure (photodoping effect) is an exciting route to tune optoelectronic phenomena. However, there is a lack of investigation concerning in what respect the action of photodoping in optoelectronic devices is local. Here, we employ scanning photocurrent microscopy (SPCM) techniques to investigate how a permanent photodoping modulates the photocurrent generation in MoS2 transistors locally. We claim that the photodoping fills the electronic states in MoS2 conduction band, preventing the photon-absorption and the photocurrent generation by the MoS2 sheet. Moreover, by comparing the persistent photocurrent (PPC) generation of MoS2 on top of different substrates, we elucidate that the interface between the material used for the gate and the insulator (gate-insulator interface) is essential for the photodoping generation. Our work gives a step forward to the understanding of the photodoping effect in MoS2 transistors and the implementation of such an effect in integrated devices. △ Less

Submitted 16 June, 2020; originally announced June 2020.

Journal ref: Nanotechnology 31 255701 (2020)

Lattice dynamics localization in low-angle twisted bilayer graphene

Authors: Andreij C. Gadelha, Douglas A. A. Ohlberg, Cassiano Rabelo, Eliel G. S. Neto, Thiago L. Vasconcelos, João L. Campos, Jessica S. Lemos, Vinícius Ornelas, Daniel Miranda, Rafael Nadas, Fabiano C. Santana, Kenji Watanabe, Takashi Taniguchi, Benoit van Troeye, Michael Lamparski, Vincent Meunier, Viet-Hung Nguyen, Dawid Paszko, Jean-Christophe Charlier, Leonardo C. Campos, Luiz G. Cançado, Gilberto Medeiros-Ribeiro, Ado Jorio

Abstract: A low twist angle between the two stacked crystal networks in bilayer graphene enables self-organized lattice reconstruction with the formation of a periodic domain. This superlattice modulates the vibrational and electronic structures, imposing new rules for electron-phonon coupling and the eventual observation of strong correlation and superconductivity. Direct optical images of the crystal supe… ▽ More A low twist angle between the two stacked crystal networks in bilayer graphene enables self-organized lattice reconstruction with the formation of a periodic domain. This superlattice modulates the vibrational and electronic structures, imposing new rules for electron-phonon coupling and the eventual observation of strong correlation and superconductivity. Direct optical images of the crystal superlattice in reconstructed twisted bilayer graphene are reported here, generated by the inelastic scattering of light in a nano-Raman spectroscope. The observation of the crystallographic structure with visible light is made possible due to lattice dynamics localization, the images resembling spectral variations caused by the presence of strain solitons and topological points. The results are rationalized by a nearly-free-phonon model and electronic calculations that highlight the relevance of solitons and topological points, particularly pronounced for structures with small twist angles. We anticipate our discovery to play a role in understanding Jahn-Teller effects and electronic Cooper pairing, among many other important phonon-related effects, and it may be useful for characterizing devices in the most prominent platform for the field of twistronics. △ Less

Submitted 16 June, 2020; originally announced June 2020.

Comments: 9 pages, 8 figures

MSC Class: 82D03 ACM Class: J.2

On hyperspherical associated Legendre functions: the extension of spherical harmonics to $N$ dimensions

Authors: L. M. B. C. Campos, M. J. S. Silva

Abstract: The solution in hyperspherical coordinates for $N$ dimensions is given for a general class of partial differential equations of mathematical physics including the Laplace, wave, heat and Helmholtz, Schrödinger, Klein-Gordon and telegraph equations and their combinations. The starting point is the Laplacian operator specified by the scale factors of hyperspherical coordinates. The general equation… ▽ More The solution in hyperspherical coordinates for $N$ dimensions is given for a general class of partial differential equations of mathematical physics including the Laplace, wave, heat and Helmholtz, Schrödinger, Klein-Gordon and telegraph equations and their combinations. The starting point is the Laplacian operator specified by the scale factors of hyperspherical coordinates. The general equation of mathematical physics is solved by separation of variables leading to the dependencies: (i) on time by the usual exponential function; (ii) on longitude by the usual sinusoidal function; (iii) on radius by Bessel functions of order generally distinct from cylindrical or spherical Bessel functions; (iv) on one latitude by associated Legendre functions; (v) on the remaining latitudes by an extension, namely the hyperspherical associated Legendre functions. The original associated Legendre functions are a particular case of the Gaussian hypergeometric functions, and the hyperspherical associated Legendre functions are also a more general particular case of the Gaussian hypergeometric functions so that it is not necessary to consider extended Gaussian hypergeometric functions. △ Less

Submitted 19 May, 2020; originally announced May 2020.

Comments: 52 pages, 2 figures, 1 table submitted to "Advances in Theoretical and Mathematical Physics"

The role of thickness inhomogeneities in hierarchical cortical folding

Authors: Lucas da Costa Campos, Raphael Hornung, Gerhard Gompper, Jens Elgeti, Svenja Caspers

Abstract: For long it has been known that specific patterns of folding are necessary for an optimally functioning brain. For instance, lissencephaly and polymicrogyria can lead to severe mental retardation, short life expectancy, epileptic seizures, and tetraplegia. The construction of a quantitative model on how and why these folds appear is the first step in understanding the cause of these conditions. In… ▽ More For long it has been known that specific patterns of folding are necessary for an optimally functioning brain. For instance, lissencephaly and polymicrogyria can lead to severe mental retardation, short life expectancy, epileptic seizures, and tetraplegia. The construction of a quantitative model on how and why these folds appear is the first step in understanding the cause of these conditions. In recent years, there have been various attempts to understand and model the mechanisms of brain folding. Previous works have shown that mechanical instabilities play a crucial role in the formation of brain folds, and that the geometry of the fetal brain is one of the main factors in dictating the folding characteristics. However, modeling higher-order folding, one of the main characteristics of the human brain, has not been fully tackled. The effects of thickness inhomogeneity in the gyrogenesis of the mammalian brain are studied through finite-element simulations of rectangular slabs. The slabs are divided into two distinct regions, where the outer layer mimics the gray matter, and the inner layer the underlying white matter. Differential growth is introduced by only growing the top layer. The brain tissue is modeled as a neo-Hookean hyperelastic material. Simulations are performed with both, homogeneous and inhomogeneous cortical thickness. The homogeneous cortex is shown to fold into a single wavelength, as common for bilayered materials, while the inhomogeneous cortex folds into more complex conformations: In their early stages of development, structures reminiscent of the deep sulci in the brain are obtained. As the cortex continues to develop, secondary undulations, which are shallower and more variable than the structures obtained in earlier gyrification stage emerge, reproducing well-known characteristics of higher-order folding in the mammalian, and particularly the human, brain. △ Less

Submitted 2 April, 2020; originally announced April 2020.

Comments: 13 pages,10 figures

Inertial-Only Optimization for Visual-Inertial Initialization

Authors: Carlos Campos, José M. M. Montiel, Juan D. Tardós

Abstract: We formulate for the first time visual-inertial initialization as an optimal estimation problem, in the sense of maximum-a-posteriori (MAP) estimation. This allows us to properly take into account IMU measurement uncertainty, which was neglected in previous methods that either solved sets of algebraic equations, or minimized ad-hoc cost functions using least squares. Our exhaustive initialization… ▽ More We formulate for the first time visual-inertial initialization as an optimal estimation problem, in the sense of maximum-a-posteriori (MAP) estimation. This allows us to properly take into account IMU measurement uncertainty, which was neglected in previous methods that either solved sets of algebraic equations, or minimized ad-hoc cost functions using least squares. Our exhaustive initialization tests on EuRoC dataset show that our proposal largely outperforms the best methods in the literature, being able to initialize in less than 4 seconds in almost any point of the trajectory, with a scale error of 5.3% on average. This initialization has been integrated into ORB-SLAM Visual-Inertial boosting its robustness and efficiency while maintaining its excellent accuracy. △ Less

Submitted 12 March, 2020; originally announced March 2020.

Comments: 2020 International Conference on Robotics and Automation

Probing the Electronic Properties of Monolayer MoS$_2$ via Interaction with Molecular Hydrogen

Authors: Natália P. Rezende, Alisson R. Cadore, Andreij C. Gadelha, Cíntia L. Pereira, Vinicius Ornelas, Kenji Watanabe, Takashi Taniguchi, André S. Ferlauto, Ângelo Malachias, Leonardo C. Campos, Rodrigo G. Lacerda

Abstract: This work presents a detailed experimental investigation of the interaction between molecular hydrogen (H$_2$) and monolayer MoS$_2$ field effect transistors (MoS$_2$ FET), aiming for sensing application. The MoS$_2$ FET exhibits a response to H$_2$ that covers a broad range of concentration (0.1 - 90%) at a relatively low operating temperature range (300-473 K). Most important, H$_2$ sensors base… ▽ More This work presents a detailed experimental investigation of the interaction between molecular hydrogen (H$_2$) and monolayer MoS$_2$ field effect transistors (MoS$_2$ FET), aiming for sensing application. The MoS$_2$ FET exhibits a response to H$_2$ that covers a broad range of concentration (0.1 - 90%) at a relatively low operating temperature range (300-473 K). Most important, H$_2$ sensors based on MoS$_2$ FETs show desirable properties such as full reversibility and absence of catalytic metal dopants (Pt or Pd). The experimental results indicate that the conductivity of MoS$_2$ monotonically increases as a function of the H$_2$ concentration due to a reversible charge transferring process. It is proposed that such process involves dissociative H$_2$ adsorption driven by interaction with sulfur vacancies in the MoS$_2$ surface (VS). This description is in agreement with related density functional theory studies about H$_2$ adsorption on MoS$_2$. Finally, measurements on partially defect-passivated MoS$_2$ FETs using atomic layer deposited aluminum oxide consist of an experimental indication that the VS plays an important role in the H$_2$ interaction with the MoS$_2$. These findings provide insights for futures applications in catalytic process between monolayer MoS$_2$ and H$_2$ and also introduce MoS$_2$ FETs as promising H$_2$ sensors. △ Less

Submitted 4 March, 2020; originally announced March 2020.

Journal ref: Advanced Electronic Materials, 2019

Reversible doping of graphene field effect transistors by molecular hydrogen: the role of the metal/graphene interface

Authors: C. L. Pereira, A. R. Cadore, N. P. Rezende, A. Gadelha, E. A. Soares, H. Chacham, L. C. Campos, R. G. Lacerda

Abstract: In this work, we present an investigation regarding how and why molecular hydrogen changes the electronic properties of graphene field effect transistors. We demonstrate that interaction with H2 leads to local doping of graphene near of the graphene-contact heterojunction. We also show that such interaction is strongly dependent on the characteristics of the metal-graphene interface. By changing t… ▽ More In this work, we present an investigation regarding how and why molecular hydrogen changes the electronic properties of graphene field effect transistors. We demonstrate that interaction with H2 leads to local doping of graphene near of the graphene-contact heterojunction. We also show that such interaction is strongly dependent on the characteristics of the metal-graphene interface. By changing the type metal in the contact, we observe that Ohmic contacts can be strongly or weakly electrostatically coupled with graphene. For strongly coupled contacts, the signature of the charge transfer effect promoted by the contacts results on an asymmetric ambipolar conduction, and such asymmetry can be tunable under interaction with H2. On the other hand, for contacts weakly coupled with graphene, the hydrogen interaction has a more profound effect. In such situation, the devices show a second charge neutrality point in graphene transistor transfer curves (a double-peak response) upon H2 exposure. We propose that this double-peak phenomenon arises from the decoupling of the work function of graphene and that of the metallic electrodes induced by the H2 molecules. We also show that the gas-induced modifications at the metal-graphene interface can be exploited to create a controlled graphene p-n junction, with considerable electron transfer to graphene layer and significant variation in the graphene resistance. These effects can pave the way for a suitable metallic contact engineering providing great potential for the application of such devices as gas sensors. △ Less

Submitted 4 March, 2020; originally announced March 2020.

Journal ref: 2D Materials, 2019

ELSI -- An Open Infrastructure for Electronic Structure Solvers

Authors: Victor Wen-zhe Yu, Carmen Campos, William Dawson, Alberto García, Ville Havu, Ben Hourahine, William P Huhn, Mathias Jacquelin, Weile Jia, Murat Keçeli, Raul Laasner, Yingzhou Li, Lin Lin, Jianfeng Lu, Jonathan Moussa, Jose E Roman, Álvaro Vázquez-Mayagoitia, Chao Yang, Volker Blum

Abstract: Routine applications of electronic structure theory to molecules and periodic systems need to compute the electron density from given Hamiltonian and, in case of non-orthogonal basis sets, overlap matrices. System sizes can range from few to thousands or, in some examples, millions of atoms. Different discretization schemes (basis sets) and different system geometries (finite non-periodic vs. infi… ▽ More Routine applications of electronic structure theory to molecules and periodic systems need to compute the electron density from given Hamiltonian and, in case of non-orthogonal basis sets, overlap matrices. System sizes can range from few to thousands or, in some examples, millions of atoms. Different discretization schemes (basis sets) and different system geometries (finite non-periodic vs. infinite periodic boundary conditions) yield matrices with different structures. The ELectronic Structure Infrastructure (ELSI) project provides an open-source software interface to facilitate the implementation and optimal use of high-performance solver libraries covering cubic scaling eigensolvers, linear scaling density-matrix-based algorithms, and other reduced scaling methods in between. In this paper, we present recent improvements and developments inside ELSI, mainly covering (1) new solvers connected to the interface, (2) matrix layout and communication adapted for parallel calculations of periodic and/or spin-polarized systems, (3) routines for density matrix extrapolation in geometry optimization and molecular dynamics calculations, and (4) general utilities such as parallel matrix I/O and JSON output. The ELSI interface has been integrated into four electronic structure code projects (DFTB+, DGDFT, FHI-aims, SIESTA), allowing us to rigorously benchmark the performance of the solvers on an equal footing. Based on results of a systematic set of large-scale benchmarks performed with Kohn-Sham density-functional theory and density-functional tight-binding theory, we identify factors that strongly affect the efficiency of the solvers, and propose a decision layer that assists with the solver selection process. Finally, we describe a reverse communication interface encoding matrix-free iterative solver strategies that are amenable, e.g., for use with planewave basis sets. △ Less

Submitted 4 July, 2020; v1 submitted 31 December, 2019; originally announced December 2019.

Journal ref: Computer Physics Communications 256 (2020) 107459

Numerical Optimal Control of HIV Transmission in Octave/MATLAB

Authors: Carlos Campos, Cristiana J. Silva, Delfim F. M. Torres

Abstract: We provide easy and readable GNU Octave/MATLAB code for the simulation of mathematical models described by ordinary differential equations and for the solution of optimal control problems through Pontryagin's maximum principle. For that, we consider a normalized HIV/AIDS transmission dynamics model based on the one proposed in our recent contribution (Silva, C.J.; Torres, D.F.M. A SICA compartment… ▽ More We provide easy and readable GNU Octave/MATLAB code for the simulation of mathematical models described by ordinary differential equations and for the solution of optimal control problems through Pontryagin's maximum principle. For that, we consider a normalized HIV/AIDS transmission dynamics model based on the one proposed in our recent contribution (Silva, C.J.; Torres, D.F.M. A SICA compartmental model in epidemiology with application to HIV/AIDS in Cape Verde. Ecol. Complex. 2017, 30, 70--75), given by a system of four ordinary differential equations. An HIV initial value problem is solved numerically using the ode45 GNU Octave function and three standard methods implemented by us in Octave/MATLAB: Euler method and second-order and fourth-order Runge-Kutta methods. Afterwards, a control function is introduced into the normalized HIV model and an optimal control problem is formulated, where the goal is to find the optimal HIV prevention strategy that maximizes the fraction of uninfected HIV individuals with the least HIV new infections and cost associated with the control measures. The optimal control problem is characterized analytically using the Pontryagin Maximum Principle, and the extremals are computed numerically by implementing a forward-backward fourth-order Runge-Kutta method. Complete algorithms, for both uncontrolled initial value and optimal control problems, developed under the free GNU Octave software and compatible with MATLAB are provided along the article. △ Less

Submitted 19 December, 2019; originally announced December 2019.

Comments: Submitted 1-Oct-2019; Revised 7-Dec-2019; Accepted 17-Dec-2019; Published 19-Dec-2019 in Math. Comput. Appl., EISSN 2297-8747 [https://doi.org/10.3390/mca25010001]

MSC Class: 34K28; 49N90; 92D30

Journal ref: Math. Comput. Appl. 25 (2020), no. 1, 20 pp

Collective motion of run-and-tumble particles drives aggregation in one-dimensional systems

Authors: C. Miguel Barriuso Gutierrez, Christian Vanhille Campos, Francisco Alarcon Oseguera, Ignacio Pagonabarraga, Ricardo Brito, Chantal Valeriani

Abstract: Active matter deals with systems whose particles consume energy at the individual level in order to move. To unravel features such as the emergence of collective structures several models have been suggested, such as the on-lattice model of run-and-tumble particles implemented via the Persistent Exclusion Process (PEP). In our work, we study a one dimensional system of run-and-tumble repulsive or… ▽ More Active matter deals with systems whose particles consume energy at the individual level in order to move. To unravel features such as the emergence of collective structures several models have been suggested, such as the on-lattice model of run-and-tumble particles implemented via the Persistent Exclusion Process (PEP). In our work, we study a one dimensional system of run-and-tumble repulsive or attractive particles, both on and off lattice. Additionally, we implement a cluster motility dynamics in the on-lattice case (since in the off-lattice case cluster motility arises from the individual particle dynamics). While we observe important differences between discrete and continuous dynamics, few common features are of particular importance. Increasing particle density drives aggregation across all different systems explored. For non-attractive particles, the effects of particle activity on aggregation are largely independent of the details of the dynamics. On the contrary, once attractive interactions are introduced, the steady state, which is completely determined by the interplay between these and the particles' activity, becomes highly dependent on the details of the dynamics. △ Less

Submitted 14 October, 2021; v1 submitted 3 December, 2019; originally announced December 2019.

NEP: a module for the parallel solution of nonlinear eigenvalue problems in SLEPc

Authors: Carmen Campos, Jose E. Roman

Abstract: SLEPc is a parallel library for the solution of various types of large-scale eigenvalue problems. In the last years we have been developing a module within SLEPc, called NEP, that is intended for solving nonlinear eigenvalue problems. These problems can be defined by means of a matrix-valued function that depends nonlinearly on a single scalar parameter. We do not consider the particular case of p… ▽ More SLEPc is a parallel library for the solution of various types of large-scale eigenvalue problems. In the last years we have been developing a module within SLEPc, called NEP, that is intended for solving nonlinear eigenvalue problems. These problems can be defined by means of a matrix-valued function that depends nonlinearly on a single scalar parameter. We do not consider the particular case of polynomial eigenvalue problems (which are implemented in a different module in SLEPc) and focus here on rational eigenvalue problems and other general nonlinear eigenproblems involving square roots or any other nonlinear function. The paper discusses how the NEP module has been designed to fit the needs of applications and provides a description of the available solvers, including some implementation details such as parallelization. Several test problems coming from real applications are used to evaluate the performance and reliability of the solvers. △ Less

Submitted 15 January, 2021; v1 submitted 25 October, 2019; originally announced October 2019.

Journal ref: ACM Trans. Math. Software, 47(3), Article 23, 2021

Fast and Robust Initialization for Visual-Inertial SLAM

Authors: Carlos Campos, J. M. M. Montiel, Juan D. Tardós

Abstract: Visual-inertial SLAM (VI-SLAM) requires a good initial estimation of the initial velocity, orientation with respect to gravity and gyroscope and accelerometer biases. In this paper we build on the initialization method proposed by Martinelli and extended by Kaiser et al. , modifying it to be more general and efficient. We improve accuracy with several rounds of visual-inertial bundle adjustment, a… ▽ More Visual-inertial SLAM (VI-SLAM) requires a good initial estimation of the initial velocity, orientation with respect to gravity and gyroscope and accelerometer biases. In this paper we build on the initialization method proposed by Martinelli and extended by Kaiser et al. , modifying it to be more general and efficient. We improve accuracy with several rounds of visual-inertial bundle adjustment, and robustify the method with novel observability and consensus tests, that discard erroneous solutions. Our results on the EuRoC dataset show that, while the original method produces scale errors up to 156%, our method is able to consistently initialize in less than two seconds with scale errors around 5%, which can be further reduced to less than 1% performing visual-inertial bundle adjustment after ten seconds. △ Less

Submitted 28 August, 2019; originally announced August 2019.

Comments: 2019 International Conference on Robotics and Automation

Journal ref: C. Campos, M. José M.M. and J. D. Tardós, "Fast and Robust Initialization for Visual-Inertial SLAM," 2019 International Conference on Robotics and Automation (ICRA), Montreal, QC, Canada, 2019, pp. 1288-1294