title: Modeling Continuous Stochastic Processes with Dynamic Normalizing Flows

abstract: Normalizing flows transform a simple base distribution into a complex target distribution and have proved to be powerful models for data generation and density estimation. In this work, we propose a novel type of normalizing flow driven by a differential deformation of the Wiener process. As a result, we obtain a rich time series model whose observable process inherits many of the appealing properties of its base process, such as efficient computation of likelihoods and marginals. Furthermore, our continuous treatment provides a natural framework for irregular time series with an independent arrival process, including straightforward interpolation. We illustrate the desirable properties of the proposed model on popular stochastic processes and demonstrate its superior flexibility to variational RNN and latent ODE baselines in a series of experiments on synthetic and real-world data.


title: A Taxonomy of Twists of Supersymmetric Yang--Mills Theory

abstract: We give a complete classification of twists of supersymmetric Yang--Mills theories in dimensions $2\leq n \leq 10$. We formulate supersymmetric Yang--Mills theory classically using the BV formalism, and then we construct an action of the supersymmetry algebra using the language of $L_\infty$ algebras. For each orbit in the space of square-zero supercharges in the supersymmetry algebra, under the action of the spin group and the group of R-symmetries, we give a description of the corresponding twisted theory. These twists can be described in terms of mixed holomorphic-topological versions of Chern--Simons and BF theory.


title: Nondecoupling in Multi-Higgs doublet models

abstract: We consider models with any number of Higgs doublets and study the conditions for decoupling. We show that, under very general circumstances, all the quadratic coefficients of the scalar potential must be present, except in special cases, which include terms related to directions of vanishing vacuum expectation values. We give a few examples. Moreover, we show that the decoupling of all charged scalars implies the decoupling of all extra neutral scalars and vanishing $\mathcal{CP}$ violation in scalar-pseudoscalar mixing.


title: Causal bounds for outcome-dependent sampling in observational studies

abstract: Outcome-dependent sampling designs are common in many different scientific fields including epidemiology, ecology, and economics. As with all observational studies, such designs often suffer from unmeasured confounding, which generally precludes the nonparametric identification of causal effects. Nonparametric bounds can provide a way to narrow the range of possible values for a nonidentifiable causal effect without making additional untestable assumptions. The nonparametric bounds literature has almost exclusively focused on settings with random sampling, and the bounds have often been derived with a particular linear programming method. We derive novel bounds for the causal risk difference, often referred to as the average treatment effect, in six settings with outcome-dependent sampling and unmeasured confounding for a binary outcome and exposure. Our derivations of the bounds illustrate two approaches that may be applicable in other settings where the bounding problem cannot be directly stated as a system of linear constraints. We illustrate our derived bounds in a real data example involving the effect of vitamin D concentration on mortality.


title: Thermodynamics and Cardy-like formula for nonminimally dressed, charged Lifshitz black holes in New Massive Gravity

abstract: In three dimensions, we consider the Einstein-Maxwell Lagrangian dressed by a nonminimally coupled scalar field in New Massive Gravity. For this theory, we provide two families of electrically charged Lifshitz black holes where their metric functions depend only on an integration constant. We calculate their masses using the quasilocal approach, as well as their entropy and electric charge. These charged configurations are interpreted as extremal in the sense that the mass vanishes identically while the entropy and electric charge are non zero thermodynamic quantities. Using these examples, we corroborate that the semiclassical entropy can be recovered through a charged Cardy-like formula, involving the corresponding magnetically charged solitons obtained by a double Wick rotation. Finally, the first law of thermodynamics, as well as the Smarr formula are also verified


title: Physics Constrained Learning for Data-driven Inverse Modeling from Sparse Observations

abstract: Deep neural networks (DNN) have been used to model nonlinear relations between physical quantities. Those DNNs are embedded in physical systems described by partial differential equations (PDE) and trained by minimizing a loss function that measures the discrepancy between predictions and observations in some chosen norm. This loss function often includes the PDE constraints as a penalty term when only sparse observations are available. As a result, the PDE is only satisfied approximately by the solution. However, the penalty term typically slows down the convergence of the optimizer for stiff problems. We present a new approach that trains the embedded DNNs while numerically satisfying the PDE constraints. We develop an algorithm that enables differentiating both explicit and implicit numerical solvers in reverse-mode automatic differentiation. This allows the gradients of the DNNs and the PDE solvers to be computed in a unified framework. We demonstrate that our approach enjoys faster convergence and better stability in relatively stiff problems compared to the penalty method. Our approach allows for the potential to solve and accelerate a wide range of data-driven inverse modeling, where the physical constraints are described by PDEs and need to be satisfied accurately.


title: MIDMod-OSN: A Microscopic-level Information Diffusion Model for Online Social Networks

abstract: As online social networks continue to be commonly used for the dissemination of information to the public, understanding the phenomena that govern information diffusion is crucial for many security and safety-related applications, such as maximizing information spread and misinformation containment during crises and natural disasters. In this study, we hypothesize that the features that contribute to information diffusion in online social networks are significantly influenced by the type of event being studied. We classify Twitter events as either informative or trending and then explore the node-to-node influence dynamics associated with information spread. We build a model based on Bayesian Logistic Regression for learning and prediction and Random Forests for feature selection. Experimental results from real-world data sets show that the proposed model outperforms state-of-the-art diffusion prediction models, achieving 93% accuracy in informative events and 86% in trending events. We observed that the models for informative and trending events differ significantly, both in the diffusion process and in the user features that govern the diffusion. Our findings show that followers play an important role in the diffusion process and it is possible to use the diffusion and OSN behavior of users for predicting the trending character of a message without having to count the number of reactions.


title: Co-VeGAN: Complex-Valued Generative Adversarial Network for Compressive Sensing MR Image Reconstruction

abstract: Compressive sensing (CS) is widely used to reduce the acquisition time of magnetic resonance imaging (MRI). Although state-of-the-art deep learning based methods have been able to obtain fast, high-quality reconstruction of CS-MR images, their main drawback is that they treat complex-valued MRI data as real-valued entities. Most methods either extract the magnitude from the complex-valued entities or concatenate them as two real-valued channels. In both the cases, the phase content, which links the real and imaginary parts of the complex-valued entities, is discarded. In order to address the fundamental problem of real-valued deep networks, i.e. their inability to process complex-valued data, we propose a novel framework based on a complex-valued generative adversarial network (Co-VeGAN). Our model can process complex-valued input, which enables it to perform high-quality reconstruction of the CS-MR images. Further, considering that phase is a crucial component of complex-valued entities, we propose a novel complex-valued activation function, which is sensitive to the phase of the input. Extensive evaluation of the proposed approach on different datasets using various sampling masks demonstrates that the proposed model significantly outperforms the existing CS-MRI reconstruction techniques in terms of peak signal-to-noise ratio as well as structural similarity index. Further, it uses significantly fewer trainable parameters to do so, as compared to the real-valued deep learning based methods.


title: Efficient exploration of zero-sum stochastic games

abstract: We investigate the increasingly important and common game-solving setting where we do not have an explicit description of the game but only oracle access to it through gameplay, such as in financial or military simulations and computer games. During a limited-duration learning phase, the algorithm can control the actions of both players in order to try to learn the game and how to play it well. After that, the algorithm has to produce a strategy that has low exploitability. Our motivation is to quickly learn strategies that have low exploitability in situations where evaluating the payoffs of a queried strategy profile is costly. For the stochastic game setting, we propose using the distribution of state-action value functions induced by a belief distribution over possible environments. We compare the performance of various exploration strategies for this task, including generalizations of Thompson sampling and Bayes-UCB to this new setting. These two consistently outperform other strategies.


title: Scalable Multi-Agent Inverse Reinforcement Learning via Actor-Attention-Critic

abstract: Multi-agent adversarial inverse reinforcement learning (MA-AIRL) is a recent approach that applies single-agent AIRL to multi-agent problems where we seek to recover both policies for our agents and reward functions that promote expert-like behavior. While MA-AIRL has promising results on cooperative and competitive tasks, it is sample-inefficient and has only been validated empirically for small numbers of agents -- its ability to scale to many agents remains an open question. We propose a multi-agent inverse RL algorithm that is more sample-efficient and scalable than previous works. Specifically, we employ multi-agent actor-attention-critic (MAAC) -- an off-policy multi-agent RL (MARL) method -- for the RL inner loop of the inverse RL procedure. In doing so, we are able to increase sample efficiency compared to state-of-the-art baselines, across both small- and large-scale tasks. Moreover, the RL agents trained on the rewards recovered by our method better match the experts than those trained on the rewards derived from the baselines. Finally, our method requires far fewer agent-environment interactions, particularly as the number of agents increases.


title: Asymptotic Analysis of Sampling Estimators for Randomized Numerical Linear Algebra Algorithms

abstract: The statistical analysis of Randomized Numerical Linear Algebra (RandNLA) algorithms within the past few years has mostly focused on their performance as point estimators. However, this is insufficient for conducting statistical inference, e.g., constructing confidence intervals and hypothesis testing, since the distribution of the estimator is lacking. In this article, we develop an asymptotic analysis to derive the distribution of RandNLA sampling estimators for the least-squares problem. In particular, we derive the asymptotic distribution of a general sampling estimator with arbitrary sampling probabilities. The analysis is conducted in two complementary settings, i.e., when the objective of interest is to approximate the full sample estimator or is to infer the underlying ground truth model parameters. For each setting, we show that the sampling estimator is asymptotically normally distributed under mild regularity conditions. Moreover, the sampling estimator is asymptotically unbiased in both settings. Based on our asymptotic analysis, we use two criteria, the Asymptotic Mean Squared Error (AMSE) and the Expected Asymptotic Mean Squared Error (EAMSE), to identify optimal sampling probabilities. Several of these optimal sampling probability distributions are new to the literature, e.g., the root leverage sampling estimator and the predictor length sampling estimator. Our theoretical results clarify the role of leverage in the sampling process, and our empirical results demonstrate improvements over existing methods.


title: Regularization by $\varepsilon$-metric. II. Limit $\varepsilon = + 0$

abstract: In a wide class of propagators regularized by the $\varepsilon$-metric [1], the $R$-operation is formulated. It is proved that the limit of renormalized Feynman integrals exists and is covariant. Possible applications in gravity are discussed. (The paper is an English translation of the second of two articles in Russian published by the author in 1987-88: V.D. Ivashchuk, Regularization by $\varepsilon$-metric. II. Limit $\varepsilon = +0$, Izvestiya Akademii Nauk Moldavskoy SSR, Ser. fiziko-tekhnicheskih i matematicheskih nauk, No. 1, p. 10-20 (1988) [in Russian] .)


title: Commutator subgroups of Sylow 2-subgroups of alternating group and Miller-Moreno groups as bases of new Key Exchange Protocol

abstract: The goal of this investigation is effective method of key exchange which based on non-commutative group $G$. The results of Ko et al. \cite{kolee} is improved and generalized. The size of a minimal generating set for the commutator subgroup of Sylow 2-subgroups of alternating group is found. The structure of the commutator subgroup of Sylow 2-subgroups of the alternating group ${A_{{2^{k}}}}$ is investigated and used in key exchange protocol which based on non-commutative group. We consider non-commutative generalization of CDH problem \cite{gu2013new, bohli2006towards} on base of metacyclic group of Miller-Moreno type (minimal non-abelian group). We show that conjugacy problem in this group is intractable. Effectivity of computation is provided due to using groups of residues by modulo $n$. The algorithm of generating (designing) common key in non-commutative group with 2 mutually commuting subgroups is constructed by us.


title: Scattered light noise characterisation at the Virgo interferometer with tvf-EMD adaptive algorithm

abstract: A methodology of adaptive time series analysis, based on Empirical Mode Decomposition (EMD), and on its time varying version tvf-EMD has been applied to strain data from the gravitational wave interferometer (IFO) Virgo in order to characterise scattered light noise affecting the sensitivity of the IFO in the detection frequency band. Data taken both during hardware injections, when a part of the IFO is put in oscillation for detector characterisation purposes, and during periods of science mode, when the IFO is fully locked and data are used for the detection of gravitational waves, were analysed. The adaptive nature of the EMD and tvf-EMD algorithms allows them to deal with nonlinear non-stationary data and hence they are particularly suited to characterise scattered light noise which is an intrinsically nonlinear and non-stationary noise. Obtained results show that tvf-EMD algorithm allows to obtain more precise results compared to the EMD algorithm, yielding higher cross-correlation values with the auxiliary channels that are the culprits of scattered light noise.


title: Cry Wolf: Toward an Experimentation Platform and Dataset for Human Factors in Cyber Security Analysis

abstract: Computer network defense is a partnership between automated systems and human cyber security analysts. The system behaviors, for example raising a high proportion of false alarms, likely impact cyber analyst performance. Experimentation in the analyst-system domain is challenging due to lack of access to security experts, the usability of attack datasets, and the training required to use security analysis tools. This paper describes Cry Wolf, an open source web application for user studies of cyber security analysis tasks. This paper also provides an open-access dataset of 73 true and false Intrusion Detection System (IDS) alarms derived from real-world examples of "impossible travel" scenarios. Cry Wolf and the impossible travel dataset were used in an experiment on the impact of IDS false alarm rate on analysts' abilities to correctly classify IDS alerts as true or false alarms. Results from that experiment are used to evaluate the quality of the dataset using difficulty and discrimination index measures drawn from classical test theory. Many alerts in the dataset provide good discrimination for participants' overall task performance.


title: Weak-strong uniqueness in weighted $L^2$ spaces and weak suitable solutions in local Morrey spaces for the MHD equations

abstract: We consider here the magneto-hydrodynamics (MHD) equations on the whole space. For the 3D case, in the setting of the weighted $L^2$ spaces we obtain a weak-strong uniqueness criterion provided that the velocity field and the magnetic field belong to a fairly general multipliers space. On the other hand, we study the local and global existence of weak suitable solutions for intermittent initial data, which is characterized through a local Morrey space. This large initial data space was also exhibit in a contemporary work [4] in the context of 3D Navier-Stokes equations. Finally, we make a discussion on the local and global existence problem in the 2D case.


title: Additive Manufacturing Approaches for Hydroxyapatite-Reinforced Composites

abstract: Additive manufacturing (AM) techniques have gained interest in the tissue engineering field thanks to their versatility and unique possibilities of producing constructs with complex macroscopic geometries and defined patterns. Recently, composite materials - namely heterogeneous biomaterials identified as continuous phase (matrix) and reinforcement (filler) - have been proposed as inks that can be processed by AM to obtain scaffolds with improved biomimetic and bioactive properties. Significant efforts have been dedicated to hydroxyapatite (HA)-reinforced composites, especially targeting bone tissue engineering, thanks to the chemical similarities of HA with respect to mineral components of native mineralized tissues. Here we review applications of AM techniques to process HA-reinforced composites and biocomposites for the production of scaffolds with biological matrices, including cellular tissues. The primary outcomes of recent investigations in terms of morphological, structural, and in vitro and in vivo biological properties of the materials are discussed. We classify the approaches based on the nature of the matrices employed to embed the HA reinforcements and produce the tissue substitutes and report a critical discussion on the presented state of the art as well as the future perspectives, to offer a comprehensive picture of the strategies investigated as well as challenges in this emerging field.


title: De novo topology optimization of Total Ossicular Replacement Prostheses

abstract: Conductive hearing loss, due to middle ear pathologies or traumas, affects more than 5% of the population worldwide. Passive prostheses to replace the ossicular chain mainly rely on piston-like titanium and/or hydroxyapatite devices, which in the long term suffer from extrusion. Although the basic shape of such devices always consists of a base for contact with the eardrum and a stem to have mechanical connection with the residual bony structures, a plethora of topologies have been proposed, mainly to help surgical positioning. In this work, we optimize the topology of a total ossicular replacement prosthesis, by maximizing the global stiffness and under the smallest possible volume constraint that ensures material continuity. This investigation optimizes the prosthesis topology in response to static displacement loads with amplitudes that normally occur during sound stimulation in a frequency range between 100 Hz and 10 kHz. Following earlier studies, we discuss how the presence and arrangement of holes on the surface of the prosthesis plate in contact with the umbo affect the overall geometry. Finally, we validate the designs through a finite-element model, in which we assess the prosthesis performance upon dynamic sound pressure loads by considering four different constitutive materials: titanium, cortical bone, silk, and collagen/hydroxyapatite. The results show that the selected prostheses present, almost independently of their constitutive material, a vibroacustic behavior close to that of the native ossicular chain, with a slight almost constant positive shift that reaches a maximum of 5 dB close to 1 kHz. This work represents a reference for the development of a new generation of middle ear prostheses with non-conventional topologies for fabrication via additive manufacturing technologies or ultraprecision machining in order to create patient-specific devices.


title: Evaluating Registration Without Ground Truth

abstract: We present a generic method for assessing the quality of non-rigid registration (NRR) algorithms, that does not depend on the existence of any ground truth, but depends solely on the data itself. The data is a set of images. The output of any NRR of such a set of images is a dense correspondence across the whole set. Given such a dense correspondence, it is possible to build various generative statistical models of appearance variation across the set. We show that evaluating the quality of the registration can be mapped to the problem of evaluating the quality of the resultant statistical model. The quality of the model entails a comparison between the model and the image data that was used to construct it. It should be noted that this approach does not depend on the specifics of the registration algorithm used (i.e., whether a groupwise or pairwise algorithm was used to register the set of images), or on the specifics of the modelling approach used. We derive an index of image model specificity that can be used to assess image model quality, and hence the quality of registration. This approach is validated by comparing our assessment of registration quality with that derived from ground truth anatomical labeling. We demonstrate that our approach is capable of assessing NRR reliably without ground truth. Finally, to demonstrate the practicality of our method, different NRR algorithms -- both pairwise and groupwise -- are compared in terms of their performance on 3D MR brain data.


title: Particle approximation of the BGK equation

abstract: In this paper we prove the convergence of a suitable particle system towards the BGK model. More precisely, we consider an interacting stochastic particle system in which each particle can instantaneously thermalize locally. We show that, under a suitable scaling limit, propagation of chaos does hold and that the one-particle distribution function converges to the solution of the BGK equation.


title: Enriched pre-Lie operads and freeness theorems

abstract: In this paper, we study the C-enriched pre-Lie operad defined by Calaque and Willwacher for any Hopf cooperad C to produce conceptual constructions of the operads acting on various deformation complexes. Maps between Hopf cooperads lead to maps between the corresponding enriched pre-Lie operads; we prove criteria for the module action of the domain on the codomain to be free, on the left and on the right. In particular, this implies a new functorial Poincar\'e--Birkhoff--Witt type theorem for universal enveloping brace algebras of pre-Lie algebras.


title: Video Monitoring Queries

abstract: Recent advances in video processing utilizing deep learning primitives achieved breakthroughs in fundamental problems in video analysis such as frame classification and object detection enabling an array of new applications. In this paper we study the problem of interactive declarative query processing on video streams. In particular we introduce a set of approximate filters to speed up queries that involve objects of specific type (e.g., cars, trucks, etc.) on video frames with associated spatial relationships among them (e.g., car left of truck). The resulting filters are able to assess quickly if the query predicates are true to proceed with further analysis of the frame or otherwise not consider the frame further avoiding costly object detection operations. We propose two classes of filters $IC$ and $OD$, that adapt principles from deep image classification and object detection. The filters utilize extensible deep neural architectures and are easy to deploy and utilize. In addition, we propose statistical query processing techniques to process aggregate queries involving objects with spatial constraints on video streams and demonstrate experimentally the resulting increased accuracy on the resulting aggregate estimation. Combined these techniques constitute a robust set of video monitoring query processing techniques. We demonstrate that the application of the techniques proposed in conjunction with declarative queries on video streams can dramatically increase the frame processing rate and speed up query processing by at least two orders of magnitude. We present the results of a thorough experimental study utilizing benchmark video data sets at scale demonstrating the performance benefits and the practical relevance of our proposals.


title: On an example of LaBuz

abstract: I acknowledge that the example of LaBuz does show that Proposition 9 in and therefore the proof of Proposition 10 is incorrect. I show what the correct "universal basis" is.


title: Efficient Rollout Strategies for Bayesian Optimization

abstract: Bayesian optimization (BO) is a class of sample-efficient global optimization methods, where a probabilistic model conditioned on previous observations is used to determine future evaluations via the optimization of an acquisition function. Most acquisition functions are myopic, meaning that they only consider the impact of the next function evaluation. Non-myopic acquisition functions consider the impact of the next $h$ function evaluations and are typically computed through rollout, in which $h$ steps of BO are simulated. These rollout acquisition functions are defined as $h$-dimensional integrals, and are expensive to compute and optimize. We show that a combination of quasi-Monte Carlo, common random numbers, and control variates significantly reduce the computational burden of rollout. We then formulate a policy-search based approach that removes the need to optimize the rollout acquisition function. Finally, we discuss the qualitative behavior of rollout policies in the setting of multi-modal objectives and model error.


title: Revisiting Saliency Metrics: Farthest-Neighbor Area Under Curve

abstract: Saliency detection has been widely studied because it plays an important role in various vision applications, but it is difficult to evaluate saliency systems because each measure has its own bias. In this paper, we first revisit the problem of applying the widely used saliency metrics on modern Convolutional Neural Networks(CNNs). Our investigation shows the saliency datasets have been built based on different choices of parameters and CNNs are designed to fit a dataset-specific distribution. Secondly, we show that the Shuffled Area Under Curve(S-AUC) metric still suffers from spatial biases. We propose a new saliency metric based on the AUC property, which aims at sampling a more directional negative set for evaluation, denoted as Farthest-Neighbor AUC(FN-AUC). We also propose a strategy to measure the quality of the sampled negative set. Our experiment shows FN-AUC can measure spatial biases, central and peripheral, more effectively than S-AUC without penalizing the fixation locations. Thirdly, we propose a global smoothing function to overcome the problem of few value degrees (output quantization) in computing AUC metrics. Comparing with random noise, our smooth function can create unique values without losing the relative saliency relationship.


title: PoisSolver: a Tool for Modelling Silicon Dangling Bond Clocking Networks

abstract: Advancements in the fabrication of silicon dangling bonds (SiDBs) reveal a potential platform for clocked field coupled nanocomputing structures. This work introduces PoisSolver, a finite element simulator for investigating clocked SiDB systems in the SiQAD design tool. Three clocking schemes borrowed from prior work on quantum-dot cellular automata are examined as potential building blocks for a general clocking framework for SiDB circuits. These clocking schemes are implemented in SiQAD, and power estimates are performed with geometrically agnostic methods to characterise each clocking scheme. Clocking schemes using a 14 nm technology node are found to dissipate 10-100 uW cm-2 at 1 GHz and 1-10 W cm-2 at 1 THz.


title: Stochastic Polyak Step-size for SGD: An Adaptive Learning Rate for Fast Convergence

abstract: We propose a stochastic variant of the classical Polyak step-size (Polyak, 1987) commonly used in the subgradient method. Although computing the Polyak step-size requires knowledge of the optimal function values, this information is readily available for typical modern machine learning applications. Consequently, the proposed stochastic Polyak step-size (SPS) is an attractive choice for setting the learning rate for stochastic gradient descent (SGD). We provide theoretical convergence guarantees for SGD equipped with SPS in different settings, including strongly convex, convex and non-convex functions. Furthermore, our analysis results in novel convergence guarantees for SGD with a constant step-size. We show that SPS is particularly effective when training over-parameterized models capable of interpolating the training data. In this setting, we prove that SPS enables SGD to converge to the true solution at a fast rate without requiring the knowledge of any problem-dependent constants or additional computational overhead. We experimentally validate our theoretical results via extensive experiments on synthetic and real datasets. We demonstrate the strong performance of SGD with SPS compared to state-of-the-art optimization methods when training over-parameterized models.


title: Variational Wasserstein Barycenters for Geometric Clustering

abstract: We propose to compute Wasserstein barycenters (WBs) by solving for Monge maps with variational principle. We discuss the metric properties of WBs and explore their connections, especially the connections of Monge WBs, to K-means clustering and co-clustering. We also discuss the feasibility of Monge WBs on unbalanced measures and spherical domains. We propose two new problems -- regularized K-means and Wasserstein barycenter compression. We demonstrate the use of VWBs in solving these clustering-related problems.


title: Provable Representation Learning for Imitation Learning via Bi-level Optimization

abstract: A common strategy in modern learning systems is to learn a representation that is useful for many tasks, a.k.a. representation learning. We study this strategy in the imitation learning setting for Markov decision processes (MDPs) where multiple experts' trajectories are available. We formulate representation learning as a bi-level optimization problem where the "outer" optimization tries to learn the joint representation and the "inner" optimization encodes the imitation learning setup and tries to learn task-specific parameters. We instantiate this framework for the imitation learning settings of behavior cloning and observation-alone. Theoretically, we show using our framework that representation learning can provide sample complexity benefits for imitation learning in both settings. We also provide proof-of-concept experiments to verify our theory.


title: Statistical inference for Axiom A attractors

abstract: From the climate system to the effect of the internet on society, chaotic systems appear to have a significant role in our future. Here a method of statistical learning for a class of chaotic systems is described along with underlying theory that can be used not only for predicting those systems a short time ahead, but also as a basis for statistical inference about their dynamics. The method is applied to prediction of 3 different systems. The statistical inference aspect can be applied to explore and enhance computer models of such systems which in turn can provide feedback for even better prediction and more precise inference.


title: Parsing Early Modern English for Linguistic Search

abstract: We investigate the question of whether advances in NLP over the last few years make it possible to vastly increase the size of data usable for research in historical syntax. This brings together many of the usual tools in NLP - word embeddings, tagging, and parsing - in the service of linguistic queries over automatically annotated corpora. We train a part-of-speech (POS) tagger and parser on a corpus of historical English, using ELMo embeddings trained over a billion words of similar text. The evaluation is based on the standard metrics, as well as on the accuracy of the query searches using the parsed data.


title: The Hills Mechanism and the Galactic Center S-stars

abstract: Our Galactic center contains young stars, including the few million year old clockwise disk between 0.05 and 0.5 pc from the Galactic center, and the S-star cluster of B-type stars at a galactocentric distance of ~0.01 pc. Recent observations suggest the S-stars are remnants of tidally disrupted binaries from the clockwise disk. In particular, Koposov et al. 2020 discovered a hypervelocity star that was ejected from the Galactic center 5 Myr ago, with a velocity vector consistent with the disk. We perform a detailed study of this binary disruption scenario. First, we quantify the plausible range of binary semimajor axes in the disk. Dynamical evaporation of such binaries is dominated by other disk stars rather than by the isotropic stellar population. For the expected range of semimajor axes in the disk, binary tidal disruptions can reproduce the observed S-star semimajor axis distribution. Reproducing the observed thermal eccentricity distribution of the S-stars requires an additional relaxation process. The flight time of the Koposov star suggests that this process must be effective within 10 Myr. We consider three possibilities: (i) scalar resonant relaxation from the observed isotropic star cluster, (ii) torques from the clockwise disk, and (iii) an intermediate-mass black hole. We conclude that the first and third mechanisms are fast enough to reproduce the observed S-star eccentricity distribution. Finally, we show that the primary star from an unequal-mass binary would be deposited at larger semimajor axes than the secondary, possibly explaining the dearth of O stars among the S-stars.


title: Adapting Planck's route to investigate the thermodynamics of the spin-half pyrochlore Heisenberg antiferromagnet

abstract: The spin-half pyrochlore Heisenberg antiferromagnet (PHAF) is one of the most challenging problems in the field of highly frustrated quantum magnetism. Stimulated by the seminal paper of M.~Planck [M.~Planck, Verhandl. Dtsch. phys. Ges. {\bf 2}, 202-204 (1900)] we calculate thermodynamic properties of this model by interpolating between the low- and high-temperature behavior. For that we follow ideas developed in detail by B.~Bernu and G.~Misguich and use for the interpolation the entropy exploiting sum rules [the ``entropy method'' (EM)]. We complement the EM results for the specific heat, the entropy, and the susceptibility by corresponding results obtained by the finite-temperature Lanczos method (FTLM) for a finite lattice of $N=32$ sites as well as by the high-temperature expansion (HTE) data. We find that due to pronounced finite-size effects the FTLM data for $N=32$ are not representative for the infinite system below $T \approx 0.7$. A similar restriction to $T \gtrsim 0.7$ holds for the HTE designed for the infinite PHAF. By contrast, the EM provides reliable data for the whole temperature region for the infinite PHAF. We find evidence for a gapless spectrum leading to a power-law behavior of the specific heat at low $T$ and for a single maximum in $c(T)$ at $T\approx 0.25$. For the susceptibility $\chi(T)$ we find indications of a monotonous increase of $\chi$ upon decreasing of $T$ reaching $\chi_0 \approx 0.1$ at $T=0$. Moreover, the EM allows to estimate the ground-state energy to $e_0\approx -0.52$.


title: Progressive Learning and Disentanglement of Hierarchical Representations

abstract: Learning rich representation from data is an important task for deep generative models such as variational auto-encoder (VAE). However, by extracting high-level abstractions in the bottom-up inference process, the goal of preserving all factors of variations for top-down generation is compromised. Motivated by the concept of "starting small", we present a strategy to progressively learn independent hierarchical representations from high- to low-levels of abstractions. The model starts with learning the most abstract representation, and then progressively grow the network architecture to introduce new representations at different levels of abstraction. We quantitatively demonstrate the ability of the presented model to improve disentanglement in comparison to existing works on two benchmark data sets using three disentanglement metrics, including a new metric we proposed to complement the previously-presented metric of mutual information gap. We further present both qualitative and quantitative evidence on how the progression of learning improves disentangling of hierarchical representations. By drawing on the respective advantage of hierarchical representation learning and progressive learning, this is to our knowledge the first attempt to improve disentanglement by progressively growing the capacity of VAE to learn hierarchical representations.


title: A codimension 2 component of the Gieseker-Petri locus

abstract: We show that the Brill-Noether locus M^3_{18,16} is an irreducible component of the Gieseker-Petri locus in genus 18 having codimension 2 in the moduli space of curves. This result disproves a conjecture predicting that the Gieseker-Petri locus is always divisorial.


title: Inversion of operator pencils on Banach space using Jordan chains when the generalized resolvent has an isolated essential singularity

abstract: We assume that the generalized resolvent for a bounded linear operator pencil mapping one Banach space onto another has an isolated essential singularity at the origin and is analytic on some annular region of the complex plane centred at the origin. In such cases the resolvent operator can be represented on the annulus by a convergent Laurent series and the spectral set has two components---a bounded component inside the inner boundary of the annulus and an unbounded component outside the outer boundary. In this paper we prove that the complementary spectral separation projections on the domain space are uniquely determined by the respective generating subspaces for the associated infinite-length generalized Jordan chains of vectors and that the domain space is the direct sum of these two subspaces. We show that the images of the generating subspaces under the mapping defined by the pencil provide a corresponding direct sum decomposition for the range space and that this is simply the decomposition defined by the complementary spectral separation projections on the range space. If the domain space has a Schauder basis we show that the separated systems of fundamental equations are reduced to two semi-infinite systems of matrix equations which can be solved recursively to obtain a basic solution and thereby determine the Laurent series coefficients for the resolvent operator on the given annular region.


title: Optimisation of Body-ground Contact for Augmenting Whole-Body Loco-manipulation of Quadruped Robots

abstract: Legged robots have great potential to perform loco-manipulation tasks, yet it is challenging to keep the robot balanced while it interacts with the environment. In this paper we study the use of additional contact points for maximising the robustness of loco-manipulation motions. Specifically, body-ground contact is studied for enhancing robustness and manipulation capabilities of quadrupedal robots. We propose to equip the robot with prongs: small legs rigidly attached to the body which ensure body-ground contact occurs in controllable point-contacts. The effect of these prongs on robustness is quantified by computing the Smallest Unrejectable Force (SUF), a measure of robustness related to Feasible Wrench Polytopes. We apply the SUF to assess the robustness of the system, and propose an effective approximation of the SUF that can be computed at near-real-time speed. We design a hierarchical quadratic programming based whole-body controller that controls stable interaction when the prongs are in contact with the ground. This novel concept of using prongs and the resulting control framework are all implemented on hardware to validate the effectiveness of the increased robustness and newly enabled loco-manipulation tasks, such as obstacle clearance and manipulation of a large object.


title: Neural Networks are Convex Regularizers: Exact Polynomial-time Convex Optimization Formulations for Two-layer Networks

abstract: We develop exact representations of training two-layer neural networks with rectified linear units (ReLUs) in terms of a single convex program with number of variables polynomial in the number of training samples and the number of hidden neurons. Our theory utilizes semi-infinite duality and minimum norm regularization. We show that ReLU networks trained with standard weight decay are equivalent to block $\ell_1$ penalized convex models. Moreover, we show that certain standard convolutional linear networks are equivalent semi-definite programs which can be simplified to $\ell_1$ regularized linear models in a polynomial sized discrete Fourier feature space.


title: Irreducibility of random polynomials of bounded degree

abstract: It is known that random monic integral polynomials of bounded degree $d$ and integral coefficients distributed uniformly and independently in $[-H,H]$ are irreducible over $\mathbb{Z}$ with probability tending to $1$ as $H\to \infty$. In this paper, we give a general criterion for guaranteeing the same conclusion under much more general coefficient distributions, allowing them to be nonuniformly and dependently distributed over arbitrary sets.


title: Coupling lattice instabilities across the interface in ultrathin oxide heterostructures

abstract: Oxide heterointerfaces constitute a rich platform for realizing novel functionalities in condensed matter. A key aspect is the strong link between structural and electronic properties, which can be modified by interfacing materials with distinct lattice symmetries. Here we determine the effect of the cubic-tetragonal distortion of $\text{SrTiO}_3$ on the electronic properties of thin films of $\text{SrIrO}_3$, a topological crystalline metal hosting a delicate interplay between spin-orbit coupling and electronic correlations. We demonstrate that below the transition temperature at 105 K, $\text{SrIrO}_3$ orthorhombic domains couple directly to tetragonal domains in $\text{SrTiO}_3$. This forces the in-phase rotational axis to lie in-plane and creates a binary domain structure in the $\text{SrIrO}_3$ film. The close proximity to the metal-insulator transition in ultrathin $\text{SrIrO}_3$ causes the individual domains to have strongly anisotropic transport properties, driven by a reduction of bandwidth along the in-phase axis. The strong structure-property relationships in perovskites make these compounds particularly suitable for static and dynamic coupling at interfaces, providing a promising route towards realizing novel functionalities in oxide heterostructures.


title: The lithium-rotation connection in the newly discovered young stellar stream Psc-Eri (Meingast 1)

abstract: Context. As a fragile element, lithium is a sensitive probe of physical processes occurring in stellar interiors. Aims. We aim at investigating the relationship between lithium abundance and rotation rate in low-mass members of the newly discovered 125~Myr-old Psc-Eri stellar stream. Methods. We obtained high resolution optical spectra and measure the equivalent width of the 607.8 nm LiI line for 40 members of the Psc-Eri stream, whose rotational periods have been derived by arXiv:1905.10588. Results. We show that a tight correlation exists between lithium content and rotation rate among the late-G to early K-type stars of the Psc-Eri stream. Fast rotators are systematically Li-rich, while slow rotators are Li-depleted. This trend mimics the one previously reported for the similar age Pleiades cluster. Conclusions. The lithium-rotation connection thus seems to be universal over a restricted effective temperature range for low-mass stars at or close to the zero-age main sequence, and does not depend on environmental conditions.


title: On the basic reproduction number in continuously structured populations

abstract: In the framework of population dynamics, the basic reproduction number R_0 is, by definition, the expected number of offspring that an individual has during its lifetime. In constant and time periodic environments it is calculated as the spectral radius of the so-called next-generation operator. In continuously structured populations defined in a Banach lattice X with concentrated states at birth one cannot define the next-generation operator in X. In the present paper we present an approach to compute the basic reproduction number of such models as the limit of the basic reproduction number of a sequence of models for which R_0 can be computed as the spectral radius of the next-generation operator. We apply these results to some examples: the (classical) size-dependent model, a size structured cell population model, a size structured model with diffusion in structure space (under some particular assumptions) and a (physiological) age-structured model with diffusion in structure space.


title: Physics-informed deep learning for incompressible laminar flows

abstract: Physics-informed deep learning has drawn tremendous interest in recent years to solve computational physics problems, whose basic concept is to embed physical laws to constrain/inform neural networks, with the need of less data for training a reliable model. This can be achieved by incorporating the residual of physics equations into the loss function. Through minimizing the loss function, the network could approximate the solution. In this paper, we propose a mixed-variable scheme of physics-informed neural network (PINN) for fluid dynamics and apply it to simulate steady and transient laminar flows at low Reynolds numbers. A parametric study indicates that the mixed-variable scheme can improve the PINN trainability and the solution accuracy. The predicted velocity and pressure fields by the proposed PINN approach are also compared with the reference numerical solutions. Simulation results demonstrate great potential of the proposed PINN for fluid flow simulation with a high accuracy.


title: Do nuclear rings in barred galaxies form at the shear minimum of the rotation curve?

abstract: It has been recently suggested that (i) nuclear rings in barred galaxies (including our own Milky Way) form at the radius where the shear parameter of the rotation curve reaches a minimum; (ii) the acoustic instability of Montenegro et al. is responsible for driving the turbulence and angular momentum transport in the central regions of barred galaxies. Here we test these suggestions by running simple hydrodynamical simulations in a logarithmic barred potential. Since the rotation curve of this potential is scale-free, the shear minimum theory predicts that no ring should form. We find that in contrast to this prediction, a ring does form in the simulation, with morphology consistent with that of nuclear rings in real barred galaxies. This proves that the presence of a shear-minimum is not a necessary condition for the formation of a ring. We also find that perturbations that are predicted to be acoustically unstable wind up and eventually propagate off to infinity, so that the system is actually stable. We conclude that (i) the shear-minimum theory is an unlikely mechanism for the formation of nuclear rings in barred galaxies; (ii) the acoustic instability is a spurious result and may not be able to drive turbulence in the interstellar medium, at least for the case without self-gravity. The question of the role of turbulent viscosity remains open.


title: Triplet Online Instance Matching Loss for Person Re-identification

abstract: Mining the shared features of same identity in different scene, and the unique features of different identity in same scene, are most significant challenges in the field of person re-identification (ReID). Online Instance Matching (OIM) loss function and Triplet loss function are main methods for person ReID. Unfortunately, both of them have drawbacks. OIM loss treats all samples equally and puts no emphasis on hard samples. Triplet loss processes batch construction in a complicated and fussy way and converges slowly. For these problems, we propose a Triplet Online Instance Matching (TOIM) loss function, which lays emphasis on the hard samples and improves the accuracy of person ReID effectively. It combines the advantages of OIM loss and Triplet loss and simplifies the process of batch construction, which leads to a more rapid convergence. It can be trained on-line when handle the joint detection and identification task. To validate our loss function, we collect and annotate a large-scale benchmark dataset (UESTC-PR) based on images taken from surveillance cameras, which contains 499 identities and 60,437 images. We evaluated our proposed loss function on Duke, Marker-1501 and UESTC-PR using ResNet-50, and the result shows that our proposed loss function outperforms the baseline methods by a maximum of 21.7%, including Softmax loss, OIM loss and Triplet loss.


title: Learning the mapping $\mathbf{x}\mapsto \sum_{i=1}^d x_i^2$: the cost of finding the needle in a haystack

abstract: The task of using machine learning to approximate the mapping $\mathbf{x}\mapsto\sum_{i=1}^d x_i^2$ with $x_i\in[-1,1]$ seems to be a trivial one. Given the knowledge of the separable structure of the function, one can design a sparse network to represent the function very accurately, or even exactly. When such structural information is not available, and we may only use a dense neural network, the optimization procedure to find the sparse network embedded in the dense network is similar to finding the needle in a haystack, using a given number of samples of the function. We demonstrate that the cost (measured by sample complexity) of finding the needle is directly related to the Barron norm of the function. While only a small number of samples is needed to train a sparse network, the dense network trained with the same number of samples exhibits large test loss and a large generalization gap. In order to control the size of the generalization gap, we find that the use of explicit regularization becomes increasingly more important as $d$ increases. The numerically observed sample complexity with explicit regularization scales as $\mathcal{O}(d^{2.5})$, which is in fact better than the theoretically predicted sample complexity that scales as $\mathcal{O}(d^{4})$. Without explicit regularization (also called implicit regularization), the numerically observed sample complexity is significantly higher and is close to $\mathcal{O}(d^{4.5})$.


title: The GAPS Programme at TNG XXI -- A GIARPS case-study of known young planetary candidates: confirmation of HD 285507 b and refutation of AD Leo b

abstract: The existence of hot Jupiters is still not well understood. Two main channels are thought to be responsible for their current location: a smooth planet migration through the proto-planetary disk or the circularization of an initial high eccentric orbit by tidal dissipation leading to a strong decrease of the semimajor axis. Different formation scenarios result in different observable effects, such as orbital parameters (obliquity/eccentricity), or frequency of planets at different stellar ages. In the context of the GAPS Young-Objects project, we are carrying out a radial velocity survey with the aim to search and characterize young hot-Jupiter planets. Our purpose is to put constraints on evolutionary models and establish statistical properties, such as the frequency of these planets from a homogeneous sample. Since young stars are in general magnetically very active, we performed multi-band (visible and near-infrared) spectroscopy with simultaneous GIANO-B + HARPS-N (GIARPS) observing mode at TNG. This helps to deal with stellar activity and distinguish the nature of radial velocity variations: stellar activity will introduce a wavelength-dependent radial velocity amplitude, whereas a Keplerian signal is achromatic. As a pilot study, we present here the cases of two already claimed hot Jupiters orbiting young stars: HD285507 b and AD Leo b. Our analysis of simultaneous high-precision GIARPS spectroscopic data confirms the Keplerian nature of HD285507's radial velocities variation and refines the orbital parameters of the hot Jupiter, obtaining an eccentricity consistent with a circular orbit. On the other hand, our analysis does not confirm the signal previously attributed to a planet orbiting AD Leo. This demonstrates the power of the multi-band spectroscopic technique when observing active stars.


title: A Double Q-Learning Approach for Navigation of Aerial Vehicles with Connectivity Constraint

abstract: This paper studies the trajectory optimization problem for an aerial vehicle with the mission of flying between a pair of given initial and final locations. The objective is to minimize the travel time of the aerial vehicle ensuring that the communication connectivity constraint required for the safe operation of the aerial vehicle is satisfied. We consider two different criteria for the connectivity constraint of the aerial vehicle which leads to two different scenarios. In the first scenario, we assume that the maximum continuous time duration that the aerial vehicle is out of the coverage of the ground base stations (GBSs) is limited to a given threshold. In the second scenario, however, we assume that the total time periods that the aerial vehicle is not covered by the GBSs is restricted. Based on these two constraints, we formulate two trajectory optimization problems. To solve these non-convex problems, we use an approach based on the double Q-learning method which is a model-free reinforcement learning technique and unlike the existing algorithms does not need perfect knowledge of the environment. Moreover, in contrast to the well-known Q-learning technique, our double Q-learning algorithm does not suffer from the over-estimation issue. Simulation results show that although our algorithm does not require prior information of the environment, it works well and shows near optimal performance.


title: OTSLM Toolbox for Structured Light Methods

abstract: We present a new Matlab toolbox for generating phase and amplitude patterns for digital micro-mirror device (DMD) and liquid crystal (LC) based spatial light modulators (SLMs). This toolbox consists of a collection of algorithms commonly used for generating patterns for these devices with a focus on optical tweezers beam shaping applications. In addition to the algorithms provided, we have put together a range of user interfaces for simplifying the use of these patterns. The toolbox currently has functionality to generate patterns which can be saved as a image or displayed on a device/screen using the supplied interface. We have only implemented interfaces for the devices our group currently uses but we believe that extending the code we provide to other devices should be fairly straightforward. The range of algorithms included in the toolbox is not exhaustive. However, by making the toolbox open sources and available on GitHub we hope that other researchers working with these devices will contribute their patterns/algorithms to the toolbox.


title: Long distance measurement using single soliton microcomb

abstract: Dispersive interferometry (DPI) takes a major interest in optical frequency comb (OFC) based long distance laser-based light detection and ranging (LIDAR) for the merits of strong anti-interference ability and long coherent length. However, the mismatch between the repetition rate of OFC and the resolution of optical spectrum acquisition system induces a large dead-zone which is a major obstacle for practical applications. Here, a new DPI LIDAR on the strength of high-repetition-rate soliton microcomb is demonstrated, which reaches a minimum Allan deviation of 27 nm for an outdoor 1179 m ranging experiment. The proposed scheme approaches a compact, high-accuracy, and none-dead-zone long distance ranging system, opening up new opportunities for emerging applications of frontier scientific researches and advanced manufacturing.


title: Forecasting the Intra-Day Spread Densities of Electricity Prices

abstract: Intra-day price spreads are of interest to electricity traders, storage and electric vehicle operators. This paper formulates dynamic density functions, based upon skewed-t and similar representations, to model and forecast the German electricity price spreads between different hours of the day, as revealed in the day-ahead auctions. The four specifications of the density functions are dynamic and conditional upon exogenous drivers, thereby permitting the location, scale and shape parameters of the densities to respond hourly to such factors as weather and demand forecasts. The best fitting and forecasting specifications for each spread are selected based on the Pinball Loss function, following the closed-form analytical solutions of the cumulative distribution functions.


title: Universal energy-accuracy tradeoffs in nonequilibrium cellular sensing

abstract: We combine stochastic thermodynamics, large deviation theory, and information theory to derive fundamental limits on the accuracy with which single cell receptors can estimate external concentrations. As expected, if estimation is performed by an ideal observer of the entire trajectory of receptor states, then no energy consuming non-equilibrium receptor that can be divided into bound and unbound states can outperform an equilibrium two-state receptor. However, when estimation is performed by a simple observer that measures the fraction of time the receptor is bound, we derive a fundamental limit on the accuracy of general nonequilibrium receptors as a function of energy consumption. We further derive and exploit explicit formulas to numerically estimate a Pareto-optimal tradeoff between accuracy and energy. We find this tradeoff can be achieved by nonuniform ring receptors with a number of states that necessarily increases with energy. Our results yield a novel thermodynamic uncertainty relation for the time a physical system spends in a pool of states, and generalize the classic 1977 Berg-Purcell limit on cellular sensing along multiple dimensions.

categories: physics.bio-ph cond-mat.stat-mech q-bio.CB q-bio.MN q-bio.SC


title: Unique Scales Preserve Self-Similar Integrate-and-Fire Functionality of Neuronal Clusters

abstract: Identifying the brain's neuronal cluster size to be presented as nodes in a network computation is critical to both neuroscience and artificial intelligence, as these define the cognitive blocks required for building intelligent computation. Experiments support many forms and sizes of neural clustering, while neural mass models (NMM) assume scale-invariant functionality. Here, we use computational simulations with brain-derived fMRI network to show that not only brain network stays structurally self-similar continuously across scales, but also neuron-like signal integration functionality is preserved at particular scales. As such, we propose a coarse-graining of network of neurons to ensemble-nodes, with multiple spikes making up its ensemble-spike, and time re-scaling factor defining its ensemble-time step. The fractal-like spatiotemporal structure and function that emerge permit strategic choice in bridging across experimental scales for computational modeling, while also suggesting regulatory constraints on developmental and/or evolutionary "growth spurts" in brain size, as per punctuated equilibrium theories in evolutionary biology.


title: Adaptive Multi-Receiver Coded Slotted ALOHA for Indoor Optical Wireless Communications

abstract: In this paper, we design a novel high-throughput random access scheme for an indoor optical wireless communication (OWC) massive Internet of Things (IoT) scenario. Due to the large density of both IoT devices and OWC access points (APs), we base the proposed scheme on multi-receiver coded slotted ALOHA. In this scenario, collisions at APs are resolved by a centralized interference cancellation decoder that may exploit both spatial and temporal diversity. By applying adaptive control of each OWC AP field of view (FOV), the proposed system is able to dynamically adapt to different IoT device activation rates, in order to maintain a high total throughput. Using illustrative simulation results, we demonstrate the design methodology and performance possibilities of the proposed method.


title: Real-time Fusion Network for RGB-D Semantic Segmentation Incorporating Unexpected Obstacle Detection for Road-driving Images

abstract: Semantic segmentation has made striking progress due to the success of deep convolutional neural networks. Considering the demands of autonomous driving, real-time semantic segmentation has become a research hotspot these years. However, few real-time RGB-D fusion semantic segmentation studies are carried out despite readily accessible depth information nowadays. In this paper, we propose a real-time fusion semantic segmentation network termed RFNet that effectively exploits complementary cross-modal information. Building on an efficient network architecture, RFNet is capable of running swiftly, which satisfies autonomous vehicles applications. Multi-dataset training is leveraged to incorporate unexpected small obstacle detection, enriching the recognizable classes required to face unforeseen hazards in the real world. A comprehensive set of experiments demonstrates the effectiveness of our framework. On Cityscapes, Our method outperforms previous state-of-the-art semantic segmenters, with excellent accuracy and 22Hz inference speed at the full 2048x1024 resolution, outperforming most existing RGB-D networks.


title: On Picent for blocks with normal defect group

abstract: We prove that if $b$ is a block of a finite group with normal abelian defect group and inertial quotient a direct product of elementary abelian groups, then $\operatorname{Picent}(b)$ is trivial. We also provide examples of blocks $b$ of finite groups with non-trivial $\operatorname{Picent}(b)$. We even have examples with normal abelian defect group and abelian inertial quotient.


title: Millimeter Wave Communications with an Intelligent Reflector: Performance Optimization and Distributional Reinforcement Learning

abstract: In this paper, a novel framework is proposed to optimize the downlink multi-user communication of a millimeter wave base station, which is assisted by a reconfigurable intelligent reflector (IR). In particular, a channel estimation approach is developed to measure the channel state information (CSI) in real-time. First, for a perfect CSI scenario, the precoding transmission of the BS and the reflection coefficient of the IR are jointly optimized, via an iterative approach, so as to maximize the sum of downlink rates towards multiple users. Next, in the imperfect CSI scenario, a distributional reinforcement learning (DRL) approach is proposed to learn the optimal IR reflection and maximize the expectation of downlink capacity. In order to model the transmission rate's probability distribution, a learning algorithm, based on quantile regression (QR), is developed, and the proposed QR-DRL method is proved to converge to a stable distribution of downlink transmission rate. Simulation results show that, in the error-free CSI scenario, the proposed approach yields over 30% and 2-fold increase in the downlink sum-rate, compared with a fixed IR reflection scheme and direct transmission scheme, respectively. Simulation results also show that by deploying more IR elements, the downlink sum-rate can be significantly improved. However, as the number of IR components increases, more time is required for channel estimation, and the slope of increase in the IR-aided transmission rate will become smaller. Furthermore, under limited knowledge of CSI, simulation results show that the proposed QR-DRL method, which learns a full distribution of the downlink rate, yields a better prediction accuracy and improves the downlink rate by 10% for online deployments, compared with a Q-learning baseline.


title: Multi Linear Regression applied to Communications systems Analysis

abstract: This paper develops a propagation model of electromagnetic signals emitted at frequencies of 20 and 40 MHz for the Ecuadorian jungle. It is expected that the results obtained at the end of this research will be applied to produce a complete coverage map for wireless communications technologies, which will optimize the radio spectrum in operations carried out by the Armed Forces in the Ecuadorian border jungle. The final expression found is an adjustment function that relates the Receiving Power (PRX) to factors that determine the geometry of the Fresnell Zone (Connectivity). The resulting model of the research improves the discrepancy between the simulated power (PRX) in commercial software and a sample of measured wireless transmissions in situ. The analysis was based on the results and methodology presented by Longley-Rice. It was determined the non-normality of the discrepancy between the losses (LLR) calculated by Longley-Rice Model (LMR) and the data obtained in the field, It was added correction coefficients on the expression of LMR. Subsequently, the mathematical expression was linearized to implement multivariate linear adjustment techniques. Alternative formulations to the Linear Regression model were sought and their goodness of fit was compared; all these techniques are introduced theoretically. To conclude, an analysis of the error of the found model is made. Mathematical modeling software such as MATLAB and SPSS were used for the formulation and numerical analysis, whose algorithms are also introduced. Finally, we propose future lines of research that allow a global understanding of the behavior of telecommunications technologies under hostile environments.


title: Generalized scale behavior and renormalization group for data analysis

abstract: Some recent results showed that renormalization group can be considered as a promising framework to address open issues in data analysis. In this work, we focus on one of these aspects, closely related to principal component analysis for the case of large dimensional data sets with covariance having a nearly continuous spectrum. In this case, the distinction between "noise-like" and "non-noise" modes becomes arbitrary and an open challenge for standard methods. Observing that both renormalization group and principal component analysis search for simplification for systems involving many degrees of freedom, we aim to use the renormalization group argument to clarify the turning point between noise and information modes. The analogy between coarse-graining renormalization and principal component analysis has been investigated in [Journal of Statistical Physics,167, Issue 3-4, pp 462-475, (2017)], from a perturbative framework, and the implementation with real sets of data by the same authors showed that the procedure may reflect more than a simple formal analogy. In particular, the separation of sampling noise modes may be controlled by a non-Gaussian fixed point, reminiscent of the behaviour of critical systems. In our analysis, we go beyond the perturbative framework using nonperturbative techniques to investigate non-Gaussian fixed points and propose a deeper formalism allowing going beyond power-law assumptions for explicit computations.


title: Oscillation of high-energy neutrinos from choked jets in stellar and merger ejecta

abstract: We present a comprehensive study on oscillation of high-energy neutrinos from two different environments: blue supergiant progenitors that may harbor low-power gamma-ray burst (GRB) jets and neutron star merger ejecta that would be associated with short gamma-ray bursts. We incorporate the radiation constraint that gives a necessary condition for nonthermal neutrino production, and account for the time evolution of the jet, which allows us to treat neutrino oscillation in matter more accurately. For massive star progenitors, neutrino injection inside the star can lead to nonadiabatic oscillation patterns in the 1 TeV - 10 TeV and is also visible in the flavor ratio. For neutron star merger ejecta, we find a similar behavior in the 100 GeV - 10 TeV region and the oscillation may result in a $\nu_e$ excess around 1 TeV. These features, which enable us to probe the progenitors of long and short GRBs, could be seen by future neutrino detectors with precise flavor ratio measurements. We also discuss potential contributions to the diffuse neutrino flux measured by IceCube, and find parameter sets allowing choked low-power GRB jets to account for the neutrino flux in the 10 TeV - 100 TeV range without violating the existing constraints.


title: The Rise of Active Galactic Nuclei in the GAlaxy Evolution and Assembly semi-analytic model

abstract: We present a new implementation of the GAlaxy Evolution and Assembly (GAEA) semi-analytic model, that features an improved modelling of the process of cold gas accretion onto supermassive black hole (SMBHs), derived from both analytic arguments and high-resolution simulations. We consider different scenarios for the loss of angular momentum required for the available cold gas to be accreted onto the central SMBHs, and we compare different combinations of triggering mechanisms, including galaxy mergers and disc instabilities in star forming discs. We compare our predictions with the luminosity function (LF) observed for Active Galactic Nuclei (AGN) and we confirm that a non-instantaneous accretion timescale (either in the form of a low-angular momentum reservoir or as an assumed light curve evolution) is needed in order to reproduce the measured evolution of the AGN-LF and the so-called AGN-downsizing trend. Moreover, we also study the impact of AGN feedback, in the form of AGN-driven outflows, on the SF properties of model galaxies, using prescriptions derived both from empirical studies or from numerical experiments. We show that AGN-driven outflows are effective in suppressing the residual star formation rate in massive galaxies ($> 10^{11} M_\odot$) without changing their overall assembly history. These winds also affect the SFR of lower mass galaxies, resulting in a too large fraction of passive galaxies at $< 10^{10} M_\odot$. Finally, we study the Eddington ratio distribution as a function of SMBH mass, showing that only objects more massive than $10^8 M_\odot$ are already in a self-regulated state as inferred from observations.


title: Dynamic Power Allocation and Virtual Cell Formation for Throughput-Optimal Vehicular Edge Networks in Highway Transportation

abstract: This paper investigates highly mobile vehicular networks from users' perspectives in highway transportation. Particularly, a centralized software-defined architecture is introduced in which centralized resources can be assigned, programmed, and controlled using the anchor nodes (ANs) of the edge servers. Unlike the legacy networks, where a typical user is served from only one access point (AP), in the proposed system model, a vehicle user is served from multiple APs simultaneously. While this increases the reliability and the spectral efficiency of the assisted users, it also necessitates an accurate power allocation in all transmission time slots. As such, a joint user association and power allocation problem is formulated to achieve enhanced reliability and weighted user sum rate. However, the formulated problem is a complex combinatorial problem, remarkably hard to solve. Therefore, fine-grained machine learning algorithms are used to efficiently optimize joint user associations and power allocations of the APs in a highly mobile vehicular network. Furthermore, a distributed single-agent reinforcement learning algorithm, namely SARL-MARL, is proposed which obtains nearly identical genie-aided optimal solutions within a nominal number of training episodes than the baseline solution. Simulation results validate that our solution outperforms existing schemes and can attain genie-aided optimal performances.


title: Observation of the semi-muonic decay $D^+\to \omega\mu^+\nu_\mu$

abstract: We report the first observation of the semi-muonic decay $D^+\to \omega \mu^+\nu_\mu$ using an $e^+e^-$ collision data sample corresponding to an integrated luminosity of $2.93~\mathrm{fb}^{-1}$ collected with the BESIII detector at a center-of-mass energy of 3.773 GeV. The absolute branching fraction of the $D^+\to \omega \mu^+\nu_\mu$ decay is measured to be ${\mathcal B}_{D^+\to \omega\mu^+\nu_\mu}=(17.7\pm1.8_{\rm stat}\pm1.1_{\rm syst})\times 10^{-4}$. Its ratio with the world average value of the branching fraction of the $D^+\to \omega e^+\nu_e$ decay probes lepton flavor universality and it is determined to be ${\mathcal B}_{D^+\to \omega \mu^+\nu_\mu}/{\mathcal B}^{\rm PDG}_{D^+\to \omega e^+\nu_e}=1.05\pm0.14$, in agreement with the standard model expectation within one standard deviation.


title: Exceptional coupling in extreme skin-depth waveguides for extremely low waveguide crosstalk

abstract: Photonic chips can miniaturize complicate optical systems very tiny and portable, providing versatile functionalities for many optical applications. Increasing the photonic chip integration density is highly desired as it provides more functionalities, low cost, and lower power consumption. However, photonic chip integration density is limited by the waveguide crosstalk, which is caused by the evanescent waves in the cladding. Here we show that the waveguide crosstalk can be suppressed completely with the exceptional coupling in extreme skin-depth (eskid) waveguides. The anisotropic dielectric perturbations in the coupled eskid waveguides cause such an exceptional coupling, resulting in infinitely long coupling length. We demonstrate the extreme suppression of waveguide crosstalk via exceptional coupling on a silicon-on-insulator (SOI) platform, which is compatible with a complementary metal-oxide-semiconductor (CMOS) process. The idea of exceptional coupling in eskid waveguides can be applied to many other photonic devices as well, significantly reducing entire chip footprints.


title: Nematic and time-reversal breaking superconductivities coexisting with quadrupole order in a $\Gamma_3$ system

abstract: We discuss superconductivity in a model on a cubic lattice for a $\Gamma_3$ non-Kramers system. In previous studies, it is revealed that $d$-wave superconductivity with $E_g$ symmetry occurs in a wide parameter range in a $\Gamma_3$ system. Such anisotropic superconductivity can break the cubic symmetry of the lattice. In a $\Gamma_3$ system, the quadrupole degrees of freedom are active and the effect of the cubic symmetry breaking should be important. Here, we investigate the coexisting states of the $d$-wave superconductivity and quadrupole order by a mean-field theory. In particular, we discuss possible competition and cooperation between the superconductivity and quadrupole order depending on types of them. We find nematic superconductivity breaking the cubic symmetry and coexisting with quadrupole order. In the present model, we also find $d+id$ superconductivity, which breaks time-reversal symmetry but retains the cubic symmetry. We also discuss the effects of uniaxial stress on these superconducting states.


title: Spin-orbit magnetic state readout in scaled ferromagnetic/heavy metal nanostructures

abstract: Efficient detection of the magnetic state at nanoscale dimensions is an important step to utilize spin logic devices for computing. Magnetoresistance effects have been hitherto used in magnetic state detection, but they suffer from energetically unfavorable scaling and do not generate an electromotive force that can be used to drive a circuit element for logic device applications. Here, we experimentally show that a favorable miniaturization law is possible via the use of spin-Hall detection of the in-plane magnetic state of a magnet. This scaling law allows us to obtain a giant signal by spin Hall effect in CoFe/Pt nanostructures and quantify an effective spin-to-charge conversion rate for the CoFe/Pt system. The spin-to-charge conversion can be described as a current source with an internal resistance, i.e., it generates an electromotive force that can be used to drive computing circuits. We predict that the spin-orbit detection of magnetic states can reach high efficiency at reduced dimensions, paving the way for scalable spin-orbit logic devices and memories.


title: Automating Discovery of Dominance in Synchronous Computer-Mediated Communication

abstract: With the advent of electronic interaction, dominance (or the assertion of control over others) has acquired new dimensions. This study investigates the dynamics and characteristics of dominance in virtual interaction by analyzing electronic chat transcripts of groups solving a hidden profile task. We investigate computer-mediated communication behavior patterns that demonstrate dominance and identify a number of relevant variables. These indicators are calculated with automatic and manual coding of text transcripts. A comparison of both sets of variables indicates that automatic text analysis methods yield similar conclusions than manual coding. These findings are encouraging to advance research in text analysis methods in general, and in the study of virtual team dominance in particular.


title: Scheduled Restart Momentum for Accelerated Stochastic Gradient Descent

abstract: Stochastic gradient descent (SGD) with constant momentum and its variants such as Adam are the optimization algorithms of choice for training deep neural networks (DNNs). Since DNN training is incredibly computationally expensive, there is great interest in speeding up the convergence. Nesterov accelerated gradient (NAG) improves the convergence rate of gradient descent (GD) for convex optimization using a specially designed momentum; however, it accumulates error when an inexact gradient is used (such as in SGD), slowing convergence at best and diverging at worst. In this paper, we propose Scheduled Restart SGD (SRSGD), a new NAG-style scheme for training DNNs. SRSGD replaces the constant momentum in SGD by the increasing momentum in NAG but stabilizes the iterations by resetting the momentum to zero according to a schedule. Using a variety of models and benchmarks for image classification, we demonstrate that, in training DNNs, SRSGD significantly improves convergence and generalization; for instance in training ResNet200 for ImageNet classification, SRSGD achieves an error rate of 20.93% vs. the benchmark of 22.13%. These improvements become more significant as the network grows deeper. Furthermore, on both CIFAR and ImageNet, SRSGD reaches similar or even better error rates with significantly fewer training epochs compared to the SGD baseline.


title: Cepheid distances based on Gaia and VMC@VISTA observations

abstract: We present new accurate Period-Luminosity (PL) and Period-Wesenheit (PW) relations in the V,J,Ks bands based on a sample of more than 4500 Cepheids in the Large Magellanic Cloud (LMC) whose photometry was obtained in the context of the VISTA Magellanic Clouds (VMC) Survey. The excellent precision of these data allows us to study the geometry of the LMC and to establish a solid baseline for extra-galactic distance scale studies. To calibrate the zero points of these PL/PW relations, we adopted Gaia Data Release 2 parallaxes for more than 2000 Milky Way Cepheids. The implications for the measurement of $H_0$ are briefly discussed.


title: Backpropamine: training self-modifying neural networks with differentiable neuromodulated plasticity

abstract: The impressive lifelong learning in animal brains is primarily enabled by plastic changes in synaptic connectivity. Importantly, these changes are not passive, but are actively controlled by neuromodulation, which is itself under the control of the brain. The resulting self-modifying abilities of the brain play an important role in learning and adaptation, and are a major basis for biological reinforcement learning. Here we show for the first time that artificial neural networks with such neuromodulated plasticity can be trained with gradient descent. Extending previous work on differentiable Hebbian plasticity, we propose a differentiable formulation for the neuromodulation of plasticity. We show that neuromodulated plasticity improves the performance of neural networks on both reinforcement learning and supervised learning tasks. In one task, neuromodulated plastic LSTMs with millions of parameters outperform standard LSTMs on a benchmark language modeling task (controlling for the number of parameters). We conclude that differentiable neuromodulation of plasticity offers a powerful new framework for training neural networks.


title: Is The Leader Robot an Adequate Sensor for Posture Estimation and Ergonomic Assessment of A Human Teleoperator?

abstract: Ergonomic assessment of human posture plays a vital role in understanding work-related safety and health. Current posture estimation approaches face occlusion challenges in teleoperation and physical human-robot interaction. We investigate if the leader robot is an adequate sensor for posture estimation in teleoperation and we introduce a new probabilistic approach that relies solely on the trajectory of the leader robot for generating observations. We model the human using a redundant, partially-observable dynamical system and we infer the posture using a standard particle filter. We compare our approach with postures from a commercial motion capture system and also two least-squares optimization approaches for human inverse kinematics. The results reveal that the proposed approach successfully estimates human postures and ergonomic risk scores comparable to those estimates from gold-standard motion capture.


title: Asymmetrical nebula of the M33 variable GR290 (WR/LBV)

abstract: Context: GR290 (M33 V0532=Romano's star) is a luminous M33 object undergoing photometric variability typical for luminous blue variable (LBV) stars. It lies inside Wolf-Rayet region in the Hertzsprung-Russell diagram and possesses a WN8 type spectrum at the light minima. Analysis of Gran Telescopio Canarias (GTC) spectra obtained in 2016 led to the conclusion that it is surrounded by an unresolved HII region formed mostly of ejected material from the central star, and disclosed the presence of a second, more extended asymmetrical emission region. Aims: The aim of this paper is to further explore the structure of the nearby environment of GR290. Methods: Long-slit spectra of GR290 were obtained with three slit orientations in the visual and red spectral regions. The emission-line distribution for each slit was analyzed. Results: We confirm the presence of an asymmetric HII region that extends ~50 pc to the south; ~30pc to the north and southeast; ~20 pc to the east and northwest and ~10pc to the west. We also present the first spectrum to be acquired of a star belonging to the neighboring OB88 association, J013501.87+304157.3, which we classify as a B-type supergiant with a possible binary companion.


title: TESS Photometry of the Eclipsing $\delta$ Scuti Star AI Hydrae

abstract: AI Hya has been known as an eclipsing binary with a monoperiodic $\delta$ Sct pulsator. We present the results from its {\it TESS} photometry observed during Sector 7. Including our five minimum epochs, the eclipse timing diagram displays the apsidal motion with a rate of $\dot{\omega}$ = 0.075$\pm$0.031 deg year$^{-1}$, which corresponds to an apsidal period of U = 4800$\pm$2000 years. The binary star model represents that the smaller, less massive primary component is 427 K hotter than the pulsating secondary, and our distance of 612$\pm$36 pc is in good agreement with the $Gaia$ distance of 644$\pm$26 pc. We subtracted the binary effects from the observed {\it TESS} data and applied a multifrequency analysis to these residuals. The result reveals that AI Hya is multiperiodic in its pulsation. Of 14 signals detected, four ($f_1$, $f_2$, $f_3$, $f_6$) may be considered independent pulsation frequencies. The period ratios of $P_{\rm pul}/P_{\rm orb}$ = 0.012$-$0.021 and the pulsation constants of $Q$ = 0.30$-$0.52 days correspond to $\delta$ Sct pulsations in binaries. We found that the secondary component of AI Hya pulsates in both radial fundamental $F$ modes ($f_2$ and $f_3$) and non-radial $g_1$ modes with a low degree of $\ell$ = 2 ($f_1$ and $f_6$).


title: Trivial 2-cocycles for invariants of mod p homology spheres and Perron's conjecture

abstract: The main target of this thesis is to solve the Perron's conjecture. This conjecture affirms that some function on the mod p Torelli group, with values in Z/p, is an invariant of mod p homology 3-spheres. In order to solve this conjecture, in this thesis we first study the mod p homology 3-spheres, the rational homology 3-spheres and those that can be realized as a Heegaard splitting with gluing map an element of the mod p Torelli group. In particular we give a criterion to determine whenever a rational homology 3-sphere has a Heegaard splitting with gluing map an element of the Torelli group mod p, and using this criterion we prove that not all mod p homology 3-spheres can be realized in such way. Next, we extend the results of the article ''Trivial cocycles and invariants of homology 3-spheres'' obtaining a construction of invariants with values to an abelian group without restrictions, from a suitable family of 2-cocycles on the Torelli group. In particular, we explain the influence of the invariant of Rohlin in the lost of uniqueness in such construction. Later, using the same tools, we obtain a construction of invariants of rational homology spheres that have a Heegaard splitting with gluing map an element of the mod p Torelli group, from a suitable family of 2-cocycles on the mod p Torelli group where appears an invariant of mod p homology spheres which does not appear in the literature, who plays the same role that Rohlin invariant in the lost of uniqueness of our construction. Finally, we prove that Perrron's conjecture is false providing a cohomological obstruction that is given by the fact that the first characteristic class of surface bundles reduced modulo p does not vanish.


title: Revised Estimates of Ocean Surface Drag in Strong Winds

abstract: Air-sea drag governs the momentum transfer between the atmosphere and the ocean, and remains largely unknown in hurricane winds. We revisit the momentum budget and eddy-covariance methods to estimate the surface drag coefficient in the laboratory. Our drag estimates agree with field measurements in low-to-moderate winds, and previous laboratory measurements in hurricane-force winds. The drag coefficient saturates at $2.6 \times 10^{-3}$ and $U_{10} \approx 25\ m\ s^{-1}$, in agreement with previous laboratory results by Takagaki et al. (2012). During our analysis, we discovered an error in the original source code used by Donelan et al. (2004). We present the corrected data and describe the correction procedure. Although the correction to the data does not change the key finding of drag saturation in strong winds, its magnitude and wind speed threshold are significantly changed. Our findings emphasize the need for an updated and unified drag parameterization based on field and laboratory data.


title: Deep learning predicts total knee replacement from magnetic resonance images

abstract: Knee Osteoarthritis (OA) is a common musculoskeletal disorder in the United States. When diagnosed at early stages, lifestyle interventions such as exercise and weight loss can slow OA progression, but at later stages, only an invasive option is available: total knee replacement (TKR). Though a generally successful procedure, only 2/3 of patients who undergo the procedure report their knees feeling ''normal'' post-operation, and complications can arise that require revision. This necessitates a model to identify a population at higher risk of TKR, particularly at less advanced stages of OA, such that appropriate treatments can be implemented that slow OA progression and delay TKR. Here, we present a deep learning pipeline that leverages MRI images and clinical and demographic information to predict TKR with AUC $0.834 \pm 0.036$ (p < 0.05). Most notably, the pipeline predicts TKR with AUC $0.943 \pm 0.057$ (p < 0.05) for patients without OA. Furthermore, we develop occlusion maps for case-control pairs in test data and compare regions used by the model in both, thereby identifying TKR imaging biomarkers. As such, this work takes strides towards a pipeline with clinical utility, and the biomarkers identified further our understanding of OA progression and eventual TKR onset.


title: Efficient Quantum Circuit Decompositions via Intermediate Qudits

abstract: Many quantum algorithms make use of ancilla, additional qubits used to store temporary information during computation, to reduce the total execution time. Quantum computers will be resource-constrained for years to come so reducing ancilla requirements is crucial. In this work, we give a method to generate ancilla out of idle qubits by placing some in higher-value states, called qudits. We show how to take a circuit with many $O(n)$ ancilla and design an ancilla-free circuit with the same asymptotic depth. Using this, we give a circuit construction for an in-place adder and a constant adder both with $O(\log n)$ depth using temporary qudits and no ancilla.


title: The Effects Of Technology Driven Information Categories On Performance In Electronic Trading Markets

abstract: Electronic trading markets have evolved rapidly with continued adoption of new technologies and growing in-formation acquisition and processing capabilities. Traditional perspectives on trading performance adopted a mono-lithic view of information. Past research and practitioner heuristics posit that adopting new technologies and incorpo-rating more information should increase price efficiency and trading performance uniformity. However, along with technological change, information dynamics have evolved significantly resulting in immense growth in data volumes, and increased complexity of information categories. The present research explores behavioral trading performance under varying information category conditions and argues that unfettered technological developments and information consumption will not necessarily lead to consistent improvement in uniformity of trading performance. In this study, we employ an artificial stock market based economic experiment to examine the role of technol-ogy driven information categories in influencing trading decisions in electronic markets. Financial electronic markets are used as an information-rich mature markets representation to analyze information category driven trading perfor-mance. The results show that a variation of information categories can influence trading performance. The findings provide a basis to better understand behavioral phenomena in electronic markets and can be used to explain anomalies as well as to manage trading performance in electronic markets.


title: On-Orbit Operations Simulator for Workload Measurement during Telerobotic Training

abstract: Training for telerobotic systems often makes heavy use of simulated platforms, which ensure safe operation during the learning process. Outer space is one domain in which such a simulated training platform would be useful, as On-Orbit Operations (O3) can be costly, inefficient, or even dangerous if not performed properly. In this paper, we present a new telerobotic training simulator for the Canadarm2 on the International Space Station (ISS), which is able to modulate workload through the addition of confounding factors such as latency, obstacles, and time pressure. In addition, multimodal physiological data is collected from subjects as they perform a task from the simulator under these different conditions. As most current workload measures are subjective, we analyse objective measures from the simulator and EEG data that can provide a reliable measure. ANOVA of task data revealed which simulator-based performance measures could predict the presence of latency and time pressure. Furthermore, EEG classification using a Riemannian classifier and Leave-One-Subject-Out cross-validation showed promising classification performance and allowed for comparison of different channel configurations and preprocessing methods. Additionally, Riemannian distance and beta power of EEG data were investigated as potential cross-trial and continuous workload measures.


title: Is there NaI in the atmosphere of HD 209458b? Effect of the centre-to-limb variation and Rossiter-McLaughlin effect in transmission spectroscopy studies

abstract: HD 209458b was the first transiting planet discovered, and the first for which an atmosphere, in particular NaI, was detected. With time, it has become one of the most frequently studied planets, with a large diversity of atmospheric studies using low- and high-resolution spectroscopy. Here, we present transit spectroscopy observations of HD 209458b using the HARPS-N and CARMENES spectrographs. We fit the Rossiter-McLaughlin effect by combining radial velocity data from both instruments (nine transits in total), measuring a projected spin-orbit angle of $-1.6\pm0.3~{\rm deg}$. We also present the analysis of high-resolution transmission spectroscopy around the NaI region at $590~{\rm nm}$, using a total of five transit observations. In contrast to previous studies where atmospheric NaI absorption is detected, we find that for all of the nights, whether individually or combined, the transmission spectra can be explained by the combination of the centre-to-limb variation and the Rossiter-McLaughlin effect. This is also observed in the time-evolution maps and transmission light curves, but at lower signal-to-noise ratio. Other strong lines such as H$\alpha$, CaII IRT, the MgI triplet region, and KI D1 are analysed, and are also consistent with the modelled effects, without considering any contribution from the exoplanet atmosphere. Thus, the transmission spectrum reveals no detectableNaI absorption in HD 209458b. We discuss how previous pioneering studies of this benchmark object may have overlooked these effects. While for some star-planet systems these effects are small, for other planetary atmospheres the results reported in the literature may require revision.


title: Dynamical decoupling of a geometric qubit

abstract: Quantum bits or qubits naturally decohere by becoming entangled with uncontrollable environments. Dynamical decoupling is thereby required to disentangle qubits from an environment by periodically reversing the qubit bases, but this causes rotation error to accumulate. Whereas a conventional qubit is rotated within the SU(2) two-level system, a geometric qubit defined in the degenerate subspace of a V-shaped SU(3) three-level system is geometrically rotated via the third ancillary level to acquire a geometric phase. We here demonstrate that, simply by introducing detuning, the dynamical decoupling of the geometric qubit on a spin triplet electron in a nitrogen-vacancy center in diamond can be made to spontaneously suppress error accumulation. The geometric dynamical decoupling extends the coherence time of the geometric qubit up to 1.9 ms, limited by the relaxation time, with 128 decoupling gates at room temperature. Our technique opens a route to holonomic quantum memory for use in various quantum applications requiring sequential operations


title: Statistical Adaptive Stochastic Gradient Methods

abstract: We propose a statistical adaptive procedure called SALSA for automatically scheduling the learning rate (step size) in stochastic gradient methods. SALSA first uses a smoothed stochastic line-search procedure to gradually increase the learning rate, then automatically switches to a statistical method to decrease the learning rate. The line search procedure ``warms up'' the optimization process, reducing the need for expensive trial and error in setting an initial learning rate. The method for decreasing the learning rate is based on a new statistical test for detecting stationarity when using a constant step size. Unlike in prior work, our test applies to a broad class of stochastic gradient algorithms without modification. The combined method is highly robust and autonomous, and it matches the performance of the best hand-tuned learning rate schedules in our experiments on several deep learning tasks.


title: Formulas in connection with parameters related to convexity of paths on three vertices: caterpillars and unit interval graphs

abstract: We present formulas to compute the P3-geodetic number, the P3-hull number and the percolation time for a caterpillar, in terms of certain sequences associated with it. In addition, we find a connection between the percolation time of a unit interval graph and a parameter involving the diameter of a unit interval graph related to it. Finally, we present a hereditary graph class, defined by forbidden induced subgraphs, such that its percolation time is equal to one.