-
Stability Properties of the Impulsive Goodwin's Oscillator in 1-cycle
Authors:
Anton V. Proskurnikov,
Alexander Medvedev
Abstract:
The Impulsive Goodwin's Oscillator (IGO) is a mathematical model of a hybrid closed-loop system. It arises by closing a special kind of continuous linear positive time-invariant system with impulsive feedback, which employs both amplitude and frequency pulse modulation. The structure of IGO precludes the existence of equilibria, and all its solutions are oscillatory. With its origin in mathematica…
▽ More
The Impulsive Goodwin's Oscillator (IGO) is a mathematical model of a hybrid closed-loop system. It arises by closing a special kind of continuous linear positive time-invariant system with impulsive feedback, which employs both amplitude and frequency pulse modulation. The structure of IGO precludes the existence of equilibria, and all its solutions are oscillatory. With its origin in mathematical biology, the IGO also presents a control paradigm useful in a wide range of applications, in particular dosing of chemicals and medicines. Since the pulse modulation feedback mechanism introduces significant nonlinearity and non-smoothness in the closedloop dynamics, conventional controller design methods fail to apply. However, the hybrid dynamics of IGO reduce to a nonlinear, time-invariant discrete-time system, exhibiting a one-to-one correspondence between periodic solutions of the original IGO and those of the discrete-time system. The paper proposes a design approach that leverages the linearization of the equivalent discrete-time dynamics in the vicinity of a fixed point. A simple and efficient local stability condition of the 1-cycle in terms of the characteristics of the amplitude and frequency modulation functions is obtained.
△ Less
Submitted 27 March, 2024;
originally announced March 2024.
-
Impulsive feedback control for dosing applications
Authors:
Alexander Medvedev,
Anton Proskurnikov,
Zhanybai T. Zhusubaliyev
Abstract:
This paper addresses a design procedure of pulse-modulated feedback control solving a dosing problem defined for implementation in a manual mode. Discrete dosing, as a control strategy, is characterized by exerting control action on the plant in impulsive manner at certain time instants. Dosing applications appear primarily in chemical industry and medicine where the control signal constitutes a s…
▽ More
This paper addresses a design procedure of pulse-modulated feedback control solving a dosing problem defined for implementation in a manual mode. Discrete dosing, as a control strategy, is characterized by exerting control action on the plant in impulsive manner at certain time instants. Dosing applications appear primarily in chemical industry and medicine where the control signal constitutes a sequence of (chemically or pharmacologically) active substance quantities (doses) administered to achieve a desired result. When the doses and the instants of their administration are adjusted as functions of some measured variable, a feedback control loop exhibiting nonlinear dynamics arises. The impulsive character of the interaction between the controller and the plant makes the resulting closed-loop system non-smooth. Limitations of the control law with respect to control goals are discussed. An application of the approach at hand to neuromuscular blockade in closed-loop anesthesia is considered in a numerical example.
△ Less
Submitted 10 November, 2023;
originally announced November 2023.
-
Output corridor control via design of impulsive Goodwin's oscillator
Authors:
Alexander Medvedev,
Anton V. Proskurnikov,
Zhanybai T. Zhusubaliyev
Abstract:
In the Impulsive Goodwin's oscillator (IGO), a continuous positive linear time-invariant (LTI) plant is controlled by an amplitude- and frequency-modulated feedback into an oscillating solution. Self-sustained oscillations in the IGO model have been extensively used to portray periodic rhythms in endocrine systems, whereas the potential of the concept as a controller design approach still remains…
▽ More
In the Impulsive Goodwin's oscillator (IGO), a continuous positive linear time-invariant (LTI) plant is controlled by an amplitude- and frequency-modulated feedback into an oscillating solution. Self-sustained oscillations in the IGO model have been extensively used to portray periodic rhythms in endocrine systems, whereas the potential of the concept as a controller design approach still remains mainly unexplored. This paper proposes an algorithm to design the feedback of the IGO so that the output of the continuous plant is kept (at stationary conditions) within a pre-defined corridor, i.e. within a bounded interval of values. The presented framework covers single-input single-output LTI plants as well as positive Wiener and Hammerstein models that often appear in process and biomedical control. A potential application of the developed impulsive control approach to a minimal Wiener model of pharmacokinetics and pharmacodynamics of a muscle relaxant used in general anesthesia is discussed.
△ Less
Submitted 18 October, 2023;
originally announced October 2023.
-
Neural Fiber Activation in Unipolar vs Bipolar Deep Brain Stimulation
Authors:
Anna Franziska Frigge,
Alexander Medvedev,
Elena Jiltsova,
Dag Nyholm
Abstract:
Deep Brain Stimulation (DBS) is an established and powerful treatment method in various neurological disorders. It involves chronically delivering electrical pulses to a certain stimulation target in the brain in order to alleviate the symptoms of a disease. Traditionally, the effect of DBS on neural tissue has been modeled based on the geometrical intersection of the static Volume of Tissue Activ…
▽ More
Deep Brain Stimulation (DBS) is an established and powerful treatment method in various neurological disorders. It involves chronically delivering electrical pulses to a certain stimulation target in the brain in order to alleviate the symptoms of a disease. Traditionally, the effect of DBS on neural tissue has been modeled based on the geometrical intersection of the static Volume of Tissue Activated (VTA) and the stimulation target. Recent studies suggest that the Dentato-Rubro-Thalamic Tract (DRTT) may serve as a potential common underlying stimulation target for tremor control in Essential Tremor (ET). However, clinical observations highlight that the therapeutic effect of DBS, especially in ET, is strongly influenced by the dynamic DBS parameters such as pulse width and frequency, as well as stimulation polarity. This study introduces a computational model to elucidate the effect of the stimulation signal shape on the DRTT under neural input. The simulation results suggest that achieving a specific pulse amplitude threshold is necessary before eliciting the therapeutic effect through adjustments in pulse widths and frequencies becomes feasible. Longer pulse widths proved more likely to induce firing, thus requiring a lower stimulation amplitude. Additionally, the modulation effect of bipolar configurations on neural traffic was found to vary significantly depending on the chosen stimulation polarity and the direction of neural traffic. Further, bipolar configurations demonstrated the ability to selectively influence firing patterns in different fiber tracts.
△ Less
Submitted 2 October, 2023;
originally announced October 2023.
-
Security Allocation in Networked Control Systems under Stealthy Attacks
Authors:
Anh Tung Nguyen,
André M. H. Teixeira,
Alexander Medvedev
Abstract:
This paper considers the problem of security allocation in a networked control system under stealthy attacks. The system is comprised of interconnected subsystems represented by vertices. A malicious adversary selects a single vertex on which to conduct a stealthy data injection attack with the purpose of maximally disrupting a distant target vertex while remaining undetected. Defense resources ag…
▽ More
This paper considers the problem of security allocation in a networked control system under stealthy attacks. The system is comprised of interconnected subsystems represented by vertices. A malicious adversary selects a single vertex on which to conduct a stealthy data injection attack with the purpose of maximally disrupting a distant target vertex while remaining undetected. Defense resources against the adversary are allocated by a defender on several selected vertices. First, the objectives of the adversary and the defender with uncertain targets are formulated in a probabilistic manner, resulting in an expected worst-case impact of stealthy attacks. Next, we provide a graph-theoretic necessary and sufficient condition under which the cost for the defender and the expected worst-case impact of stealthy attacks are bounded. This condition enables the defender to restrict the admissible actions to dominating sets of the graph representing the network. Then, the security allocation problem is solved through a Stackelberg game-theoretic framework. Finally, the obtained results are validated through a numerical example of a 50-vertex networked control system.
△ Less
Submitted 2 April, 2024; v1 submitted 31 August, 2023;
originally announced August 2023.
-
Closed-form expressions for the pure time delay in terms of the input and output Laguerre spectra
Authors:
Alexander Medvedev
Abstract:
The pure time delay operator is considered in continuous and discrete time under the assumption of the input signal being integrable (summable) with square. By making use of a discrete convolution operator with polynomial Markov parameters, a common framework for handling the continuous and discrete case is set. Closed-form expressions for the delay value are derived in terms of the Laguerre spect…
▽ More
The pure time delay operator is considered in continuous and discrete time under the assumption of the input signal being integrable (summable) with square. By making use of a discrete convolution operator with polynomial Markov parameters, a common framework for handling the continuous and discrete case is set. Closed-form expressions for the delay value are derived in terms of the Laguerre spectra of the output and input signals. The expressions hold for any feasible value of the Laguerre parameter and can be utilized for e.g. building time-delay estimators that allow for non-persistent input. A simulation example is provided to illustrate the principle of Laguerre-domain time delay modeling and analysis.
△ Less
Submitted 20 June, 2023;
originally announced June 2023.
-
Design of the Impulsive Goodwin's Oscillator: A Case Study
Authors:
Alexander Medvedev,
Anton V. Proskurnikov,
Zhanybai T. Zhusubaliyev
Abstract:
The impulsive Goodwin's oscillator (IGO) is a hybrid model composed of a third-order continuous linear part and a pulse-modulated feedback. This paper introduces a design problem of the IGO to admit a desired periodic solution. The dynamics of the continuous states represent the plant to be controlled, whereas the parameters of the impulsive feedback constitute design degrees of freedom. The desig…
▽ More
The impulsive Goodwin's oscillator (IGO) is a hybrid model composed of a third-order continuous linear part and a pulse-modulated feedback. This paper introduces a design problem of the IGO to admit a desired periodic solution. The dynamics of the continuous states represent the plant to be controlled, whereas the parameters of the impulsive feedback constitute design degrees of freedom. The design objective is to select the free parameters so that the IGO exhibits a stable 1-cycle with desired characteristics. The impulse-to-impulse map of the oscillator is demonstrated to always possess a positive fixed point that corresponds to the desired periodic solution; the closed-form expressions to evaluate this fixed point are provided. Necessary and sufficient conditions for orbital stability of the 1-cycle are presented in terms of the oscillator parameters and exhibit similarity to the problem of static output control. An IGO design procedure is proposed and validated by simulation. The nonlinear dynamics of the designed IGO are reviewed by means of bifurcation analysis. Applications of the design procedure to dosing problems in chemical industry and biomedicine are envisioned.
△ Less
Submitted 18 May, 2023;
originally announced May 2023.
-
On the continuity assumption of "Finite Adaptability in Multistage Linear Optimization'' by Bertsimas and Caramanis
Authors:
Safia Kedad-Sidhoum,
Anton Medvedev,
Frédéric Meunier
Abstract:
Two-stage robust optimization is a fundamental paradigm for modeling and solving optimization problems with uncertain parameters. A now classical method within this paradigm is finite-adaptability, introduced by Bertsimas and Caramanis (IEEE Transactions on Automatic Control, 2010). In this note, we point out that the continuity assumption they stated to ensure the convergence of the method is not…
▽ More
Two-stage robust optimization is a fundamental paradigm for modeling and solving optimization problems with uncertain parameters. A now classical method within this paradigm is finite-adaptability, introduced by Bertsimas and Caramanis (IEEE Transactions on Automatic Control, 2010). In this note, we point out that the continuity assumption they stated to ensure the convergence of the method is not correct, and we propose an alternative assumption for which we prove the desired convergence.
△ Less
Submitted 17 May, 2024; v1 submitted 9 May, 2023;
originally announced May 2023.
-
Robust One-Step Estimation of Impulsive Time Series
Authors:
Håkan Runvik,
Alexander Medvedev
Abstract:
The paper deals with the estimation of a signal model in the form of the output of a continuous linear time-invariant system driven by a sequence of instantaneous impulses, i.e. an impulsive time series. This modeling concept arises in, e.g., endocrinology when episodic hormone secretion events and elimination rates are simultaneously estimated from sampled hormone concentration measurements. The…
▽ More
The paper deals with the estimation of a signal model in the form of the output of a continuous linear time-invariant system driven by a sequence of instantaneous impulses, i.e. an impulsive time series. This modeling concept arises in, e.g., endocrinology when episodic hormone secretion events and elimination rates are simultaneously estimated from sampled hormone concentration measurements. The pulsatile secretion is modeled with a train of Dirac impulses constituting the input to a linear plant, which represents stimulated hormone secretion and elimination. A previously developed one-step estimation algorithm effectively resolves the trade-off between data fit and impulsive input sparsity. The present work improves the algorithm so that it requires less manual tuning and produces more accurate results through the use of an information criterion. It is also extended to handle outliers and unknown basal levels that are commonly recognized issues in biomedical data. The algorithm performance is evaluated both theoretically and experimentally on synthetic and clinical data.
△ Less
Submitted 26 April, 2023;
originally announced April 2023.
-
Jupiter Science Enabled by ESA's Jupiter Icy Moons Explorer
Authors:
Leigh N. Fletcher,
Thibault Cavalié,
Davide Grassi,
Ricardo Hueso,
Luisa M. Lara,
Yohai Kaspi,
Eli Galanti,
Thomas K. Greathouse,
Philippa M. Molyneux,
Marina Galand,
Claire Vallat,
Olivier Witasse,
Rosario Lorente,
Paul Hartogh,
François Poulet,
Yves Langevin,
Pasquale Palumbo,
G. Randall Gladstone,
Kurt D. Retherford,
Michele K. Dougherty,
Jan-Erik Wahlund,
Stas Barabash,
Luciano Iess,
Lorenzo Bruzzone,
Hauke Hussmann
, et al. (25 additional authors not shown)
Abstract:
ESA's Jupiter Icy Moons Explorer (JUICE) will provide a detailed investigation of the Jovian system in the 2030s, combining a suite of state-of-the-art instruments with an orbital tour tailored to maximise observing opportunities. We review the Jupiter science enabled by the JUICE mission, building on the legacy of discoveries from the Galileo, Cassini, and Juno missions, alongside ground- and spa…
▽ More
ESA's Jupiter Icy Moons Explorer (JUICE) will provide a detailed investigation of the Jovian system in the 2030s, combining a suite of state-of-the-art instruments with an orbital tour tailored to maximise observing opportunities. We review the Jupiter science enabled by the JUICE mission, building on the legacy of discoveries from the Galileo, Cassini, and Juno missions, alongside ground- and space-based observatories. We focus on remote sensing of the climate, meteorology, and chemistry of the atmosphere and auroras from the cloud-forming weather layer, through the upper troposphere, into the stratosphere and ionosphere. The Jupiter orbital tour provides a wealth of opportunities for atmospheric and auroral science: global perspectives with its near-equatorial and inclined phases, sampling all phase angles from dayside to nightside, and investigating phenomena evolving on timescales from minutes to months. The remote sensing payload spans far-UV spectroscopy (50-210 nm), visible imaging (340-1080 nm), visible/near-infrared spectroscopy (0.49-5.56 $μ$m), and sub-millimetre sounding (near 530-625\,GHz and 1067-1275\,GHz). This is coupled to radio, stellar, and solar occultation opportunities to explore the atmosphere at high vertical resolution; and radio and plasma wave measurements of electric discharges in the Jovian atmosphere and auroras. Cross-disciplinary scientific investigations enable JUICE to explore coupling processes in giant planet atmospheres, to show how the atmosphere is connected to (i) the deep circulation and composition of the hydrogen-dominated interior; and (ii) to the currents and charged particle environments of the external magnetosphere. JUICE will provide a comprehensive characterisation of the atmosphere and auroras of this archetypal giant planet.
△ Less
Submitted 26 October, 2023; v1 submitted 20 April, 2023;
originally announced April 2023.
-
Neural Network Algorithm for Intercepting Targets Moving Along Known Trajectories by a Dubins' Car
Authors:
Ivan Nasonov,
Andrey Galyaev,
Andrey Medvedev
Abstract:
The task of intercepting a target moving along a rectilinear or circular trajectory by a Dubins' car is formulated as a time-optimal control problem with an arbitrary direction of the car's velocity at the interception moment. To solve this problem and to synthesize interception trajectories, neural network methods of unsupervised learning based on the Deep Deterministic Policy Gradient algorithm…
▽ More
The task of intercepting a target moving along a rectilinear or circular trajectory by a Dubins' car is formulated as a time-optimal control problem with an arbitrary direction of the car's velocity at the interception moment. To solve this problem and to synthesize interception trajectories, neural network methods of unsupervised learning based on the Deep Deterministic Policy Gradient algorithm are used. The analysis of the obtained control laws and interception trajectories in comparison with the analytical solutions of the interception problem is performed. The mathematical modeling for the parameters of the target movement that the neural network had not seen before during training is carried out. Model experiments are conducted to test the stability of the neural solution. The effectiveness of using neural network methods for the synthesis of interception trajectories for given classes of target movements is shown.
△ Less
Submitted 12 April, 2023;
originally announced April 2023.
-
Optimal Detector Placement in Networked Control Systems under Cyber-attacks with Applications to Power Networks
Authors:
Anh Tung Nguyen,
Sribalaji C. Anand,
André M. H. Teixeira,
Alexander Medvedev
Abstract:
This paper proposes a game-theoretic method to address the problem of optimal detector placement in a networked control system under cyber-attacks. The networked control system is composed of interconnected agents where each agent is regulated by its local controller over unprotected communication, which leaves the system vulnerable to malicious cyber-attacks. To guarantee a given local performanc…
▽ More
This paper proposes a game-theoretic method to address the problem of optimal detector placement in a networked control system under cyber-attacks. The networked control system is composed of interconnected agents where each agent is regulated by its local controller over unprotected communication, which leaves the system vulnerable to malicious cyber-attacks. To guarantee a given local performance, the defender optimally selects a single agent on which to place a detector at its local controller with the purpose of detecting cyber-attacks. On the other hand, an adversary optimally chooses a single agent on which to conduct a cyber-attack on its input with the aim of maximally worsening the local performance while remaining stealthy to the defender. First, we present a necessary and sufficient condition to ensure that the maximal attack impact on the local performance is bounded, which restricts the possible actions of the defender to a subset of available agents. Then, by considering the maximal attack impact on the local performance as a game payoff, we cast the problem of finding optimal actions of the defender and the adversary as a zero-sum game. Finally, with the possible action sets of the defender and the adversary, an algorithm is devoted to determining the Nash equilibria of the zero-sum game that yield the optimal detector placement. The proposed method is illustrated on an IEEE benchmark for power systems.
△ Less
Submitted 12 April, 2023;
originally announced April 2023.
-
Context Query Simulation for Smart Carparking Scenarios in the Melbourne CDB
Authors:
Shakthi Weerasinghe,
Arkaday Zaslavsky,
Alireza Hassani,
Seng W. Loke,
Alexey Medvedev,
Amin Abken
Abstract:
The rapid growth in Internet of Things (IoT) has ushered in the way for better context-awareness enabling more smarter applications. Although for the growth in the number of IoT devices, Context Management Platforms (CMPs) that integrate different domains of IoT to produce context information lacks scalability to cater to a high volume of context queries. Research in scalability and adaptation in…
▽ More
The rapid growth in Internet of Things (IoT) has ushered in the way for better context-awareness enabling more smarter applications. Although for the growth in the number of IoT devices, Context Management Platforms (CMPs) that integrate different domains of IoT to produce context information lacks scalability to cater to a high volume of context queries. Research in scalability and adaptation in CMPs are of significant importance due to this reason. However, there is limited methods to benchmarks and validate research in this area due to the lack of sizable sets of context queries that could simulate real-world situations, scenarios, and scenes. Commercially collected context query logs are not publicly accessible and deploying IoT devices, and context consumers in the real-world at scale is expensive and consumes a significant effort and time. Therefore, there is a need to develop a method to reliably generate and simulate context query loads that resembles real-world scenarios to test CMPs for scale. In this paper, we propose a context query simulator for the context-aware smart car parking scenario in Melbourne Central Business District in Australia. We present the process of generating context queries using multiple real-world datasets and publicly accessible reports, followed by the context query execution process. The context query generator matches the popularity of places with the different profiles of commuters, preferences, and traffic variations to produce a dataset of context query templates containing 898,050 records. The simulator is executable over a seven-day profile which far exceeds the simulation time of any IoT system simulator. The context query generation process is also generic and context query language independent.
△ Less
Submitted 13 February, 2023;
originally announced February 2023.
-
Cycles in Impulsive Goodwin's Oscillators of Arbitrary Order
Authors:
Anton V. Proskurnikov,
Håkan Runvik,
Alexander Medvedev
Abstract:
Existence of periodical solutions, i.e. cycles, in the Impulsive Goodwin's Oscillator (IGO) with the continuous part of an arbitrary order m is considered. The original IGO with a third-order continuous part is a hybrid model that portrays a chemical or biochemical system composed of three substances represented by their concentrations and arranged in a cascade. The first substance in the chain is…
▽ More
Existence of periodical solutions, i.e. cycles, in the Impulsive Goodwin's Oscillator (IGO) with the continuous part of an arbitrary order m is considered. The original IGO with a third-order continuous part is a hybrid model that portrays a chemical or biochemical system composed of three substances represented by their concentrations and arranged in a cascade. The first substance in the chain is introduced via an impulsive feedback where both the impulse frequency and weights are modulated by the measured output of the continuous part. It is shown that, under the standard assumptions on the IGO, a positive periodic solution with one firing of the pulse-modulated feedback in the least period also exists in models with any m >= 1. Furthermore, the uniqueness of this 1-cycle is proved for the IGO with m <= 10 whereas, for m > 10, the uniqueness can still be guaranteed under mild assumptions on the frequency modulation function.
△ Less
Submitted 2 February, 2023;
originally announced February 2023.
-
From Traditional Adaptive Data Caching to Adaptive Context Caching: A Survey
Authors:
Shakthi Weerasinghe,
Arkady Zaslavsky,
Seng W. Loke,
Alireza Hassani,
Amin Abken,
Alexey Medvedev
Abstract:
Context information is in demand more than ever with the rapid increase in the number of context-aware Internet of Things applications developed worldwide. Research in context and context-awareness is being conducted to broaden its applicability in light of many practical and technical challenges. One of the challenges is improving performance when responding to a large number of context queries.…
▽ More
Context information is in demand more than ever with the rapid increase in the number of context-aware Internet of Things applications developed worldwide. Research in context and context-awareness is being conducted to broaden its applicability in light of many practical and technical challenges. One of the challenges is improving performance when responding to a large number of context queries. Context Management Platforms that infer and deliver context to applications measure this problem using Quality of Service (QoS) parameters. Although caching is a proven way to improve QoS, transiency of context and features such as variability and heterogeneity of context queries pose an additional real-time cost management problem. This paper presents a critical survey of the state-of-the-art in adaptive data caching with the objective of developing a body of knowledge in cost- and performance-efficient adaptive caching strategies. We comprehensively survey a large number of research publications and evaluate, compare, and contrast different techniques, policies, approaches, and schemes in adaptive caching. Our critical analysis is motivated by the focus on adaptively caching context as a core research problem. A formal definition for adaptive context caching is then proposed, followed by identified features and requirements of a well-designed, objective optimal adaptive context caching strategy.
△ Less
Submitted 9 February, 2023; v1 submitted 21 November, 2022;
originally announced November 2022.
-
Continuous Time-Delay Estimation From Sampled Measurements
Authors:
Mohamed Abdalmoaty,
Alexander Medvedev
Abstract:
An algorithm for continuous time-delay estimation from sampled output data and known input of finite energy is presented. The continuous time-delay modeling allows for the estimation of subsample delays. The proposed estimation algorithm consists of two steps. First, the continuous Laguerre spectrum of the output signal is estimated from discrete-time (sampled) noisy measurements. Second, an estim…
▽ More
An algorithm for continuous time-delay estimation from sampled output data and known input of finite energy is presented. The continuous time-delay modeling allows for the estimation of subsample delays. The proposed estimation algorithm consists of two steps. First, the continuous Laguerre spectrum of the output signal is estimated from discrete-time (sampled) noisy measurements. Second, an estimate of the delay value is obtained in Laguerre domain given a continuous-time description of the input. The second step of the algorithm is shown to be intrinsically biased, the bias sources are established, and the bias itself is modeled. The proposed delay estimation approach is compared in a Monte-Carlo simulation with state-of-the-art methods implemented in time, frequency, and Laguerre domain demonstrating comparable or higher accuracy for the considered case.
△ Less
Submitted 20 November, 2022;
originally announced November 2022.
-
Characterization of the Thermospheric Mean Winds and Circulation during Solstice using ICON/MIGHTI Observations
Authors:
Erdal Yiğit,
Manbharat Dhadly,
Alexander S. Medvedev,
Brian J. Harding,
Christoph R. Englert,
Qian Wu,
Thomas J. Immel
Abstract:
Using the horizontal neutral wind observations from the MIGHTI instrument onboard NASA's ICON (Ionospheric Connection Explorer) spacecraft with continuous coverage, we determine the climatology of the mean zonal and meridional winds and the associated mean circulation at low- to middle latitudes ($10^\circ$S-40$^{\circ}$N) for Northern Hemisphere {summer} solstice conditions between 90 km and 200…
▽ More
Using the horizontal neutral wind observations from the MIGHTI instrument onboard NASA's ICON (Ionospheric Connection Explorer) spacecraft with continuous coverage, we determine the climatology of the mean zonal and meridional winds and the associated mean circulation at low- to middle latitudes ($10^\circ$S-40$^{\circ}$N) for Northern Hemisphere {summer} solstice conditions between 90 km and 200 km altitudes, specifically on 20 June 2020 solstice as well as for a one-month period from 8 June-7 July 2020 {and for Northern winter season from 16 December 2019-31 January 2020, which spans a 47-day period, providing full local time coverage}. The data are averaged within appropriate altitude, longitude, latitude, solar zenith angle, and local time bins to produce mean wind distributions. The geographical distributions and local time variations of the mean horizontal circulation are evaluated. The instantaneous horizontal winds exhibit a significant degree of spatiotemporal variability often exceeding $\pm 150 $ m s$^{-1}$. The daily averaged zonal mean winds demonstrate day-to-day variability. Eastward zonal winds and northward (winter-to-summer) meridional winds are prevalent in the lower thermosphere, which provides indirect observational evidence of the eastward momentum deposition by small-scale gravity waves. The mean neutral winds and circulation exhibit smaller scale structures in the lower thermosphere (90-120 km), while they are more homogeneous in the upper thermosphere, indicating the increasingly dissipative nature of the thermosphere. The mean wind and circulation patterns inferred from ICON/MIGHTI measurements can be used to constrain and validate general circulation models, as well as input for numerical wave models.
△ Less
Submitted 17 October, 2022;
originally announced October 2022.
-
Noise reduction in Laguerre-domain discrete delay estimation
Authors:
Mohamed Abdalmoaty,
Alexander Medvedev
Abstract:
This paper introduces a stochastic framework for a recently proposed discrete-time delay estimation method in Laguerre-domain, i.e. with the delay block input and output signals being represented by the corresponding Laguerre series. A novel Laguerre domain disturbance model is devised, which allows the involved signals to be square-summable sequences and is suitable in a number of important appli…
▽ More
This paper introduces a stochastic framework for a recently proposed discrete-time delay estimation method in Laguerre-domain, i.e. with the delay block input and output signals being represented by the corresponding Laguerre series. A novel Laguerre domain disturbance model is devised, which allows the involved signals to be square-summable sequences and is suitable in a number of important applications. The relation to two commonly used time-domain disturbance models is clarified. Furthermore, by forming the input signal in a certain way, the signal shape of an additive output disturbance can be estimated and utilized for noise reduction. It is demonstrated that a significant improvement in the delay estimation error is achieved when the noise sequence is correlated. The noise reduction approach is applicable to other Laguerre-domain problems than pure delay estimation.
△ Less
Submitted 26 July, 2022;
originally announced July 2022.
-
Input Sequence and Parameter Estimation in Impulsive Biomedical Models
Authors:
Håkan Runvik,
Alexander Medvedev
Abstract:
A hybrid model for biomedical time series comprising a continuous second-order linear time-invariant system driven by an input sequence of positively weighted Dirac delta-functions is considered. The problem of the joint estimation of the input sequence and the continuous system parameters from output measurements is investigated. A solution that builds upon and refines a previously published leas…
▽ More
A hybrid model for biomedical time series comprising a continuous second-order linear time-invariant system driven by an input sequence of positively weighted Dirac delta-functions is considered. The problem of the joint estimation of the input sequence and the continuous system parameters from output measurements is investigated. A solution that builds upon and refines a previously published least-squares formulation is proposed. Based on a thorough analysis of the properties of the least-squares solution, improvements in terms of accuracy and ease of use are achieved on synthetic data, compared to the original algorithm.
△ Less
Submitted 4 July, 2022;
originally announced July 2022.
-
A Single-Adversary-Single-Detector Zero-Sum Game in Networked Control Systems
Authors:
Anh Tung Nguyen,
André M. H. Teixeira,
Alexander Medvedev
Abstract:
This paper proposes a game-theoretic approach to address the problem of optimal sensor placement for detecting cyber-attacks in networked control systems. The problem is formulated as a zero-sum game with two players, namely a malicious adversary and a detector. Given a protected target vertex, the detector places a sensor at a single vertex to monitor the system and detect the presence of the adv…
▽ More
This paper proposes a game-theoretic approach to address the problem of optimal sensor placement for detecting cyber-attacks in networked control systems. The problem is formulated as a zero-sum game with two players, namely a malicious adversary and a detector. Given a protected target vertex, the detector places a sensor at a single vertex to monitor the system and detect the presence of the adversary. On the other hand, the adversary selects a single vertex through which to conduct a cyber-attack that maximally disrupts the target vertex while remaining undetected by the detector. As our first contribution, for a given pair of attack and monitor vertices and a known target vertex, the game payoff function is defined as the output-to-output gain of the respective system. Then, the paper characterizes the set of feasible actions by the detector that ensures bounded values of the game payoff. Finally, an algebraic sufficient condition is proposed to examine whether a given vertex belongs to the set of feasible monitor vertices. The optimal sensor placement is then determined by computing the mixed-strategy Nash equilibrium of the zero-sum game through linear programming. The approach is illustrated via a numerical example of a 10-vertex networked control system with a given target vertex.
△ Less
Submitted 27 May, 2022;
originally announced May 2022.
-
Bayesian Monitoring of COVID-19 in Sweden
Authors:
Robin Marin,
Håkan Runvik,
Alexander Medvedev,
Stefan Engblom
Abstract:
In an effort to provide regional decision support for the public healthcare, we design a data-driven compartment-based model of COVID-19 in Sweden. From national hospital statistics we derive parameter priors, and we develop linear filtering techniques to drive the simulations given data in the form of daily healthcare demands. We additionally propose a posterior marginal estimator which provides…
▽ More
In an effort to provide regional decision support for the public healthcare, we design a data-driven compartment-based model of COVID-19 in Sweden. From national hospital statistics we derive parameter priors, and we develop linear filtering techniques to drive the simulations given data in the form of daily healthcare demands. We additionally propose a posterior marginal estimator which provides for an improved temporal resolution of the reproduction number estimate as well as supports robustness checks via a parametric bootstrap procedure.
From our computational approach we obtain a Bayesian model of predictive value which provides important insight into the progression of the disease, including estimates of the effective reproduction number, the infection fatality rate, and the regional-level immunity. We successfully validate our posterior model against several different sources, including outputs from extensive screening programs. Since our required data in comparison is easy and non-sensitive to collect, we argue that our approach is particularly promising as a tool to support monitoring and decisions within public health.
△ Less
Submitted 28 July, 2023; v1 submitted 2 May, 2022;
originally announced May 2022.
-
Skew-invariant curves and the algebraic independence of Mahler functions
Authors:
Alice Medvedev,
Khoa Dang Nguyen,
Thomas Scanlon
Abstract:
For $p \in \mathbb{Q}_+ \smallsetminus \{ 1 \}$ a positive rational number different from one, we say that the Puisseux series $f \in \mathbb{C}((t))^\text{alg}$ is $p$-Mahler of non-exceptional polynomial type if there is a polynomial $P \in \mathbb{C}(t)^\text{alg}[X]$ of degree at least two which is not conjugate to either a monomial or to plus or minus a Chebyshev polynomial for which the equa…
▽ More
For $p \in \mathbb{Q}_+ \smallsetminus \{ 1 \}$ a positive rational number different from one, we say that the Puisseux series $f \in \mathbb{C}((t))^\text{alg}$ is $p$-Mahler of non-exceptional polynomial type if there is a polynomial $P \in \mathbb{C}(t)^\text{alg}[X]$ of degree at least two which is not conjugate to either a monomial or to plus or minus a Chebyshev polynomial for which the equation $f(t^p) = P(f(t))$ holds. We show that if $p$ and $q$ are multiplicatively independent and $f$ and $g$ are $p$-Mahler and $q$-Mahler, respectively, of non-exceptional polynomial type, then $f$ and $g$ are algebraically independent over $\mathbb{C}(t)$. This theorem is proven as a consequence of a more general theorem that if $f$ is $p$-Mahler of non-exceptional polynomial type, and $g_1, \ldots, g_n$ each satisfy some difference equation with respect to the substitution $t \mapsto t^q$, then $f$ is algebraically independent from $g_1, \ldots, g_n$. These theorems are themselves consequences of a refined classification of skew-invariant curves for split polynomial dynamical systems on $\mathbb{A}^2$.
△ Less
Submitted 9 March, 2022;
originally announced March 2022.
-
Modeling the hydrological cycle in the atmosphere of Mars: Influence of a bimodal size distribution of aerosol nucleation particles
Authors:
Dmitry S. Shaposhnikov,
Alexander V. Rodin,
Alexander S. Medvedev,
Anna A. Fedorova,
Takeshi Kurod,
Paul Hartogh
Abstract:
We present a new implementation of the hydrological cycle scheme into a general circulation model of the Martian atmosphere. The model includes a semi-Lagrangian transport scheme for water vapor and ice, and accounts for microphysics of phase transitions between them. The hydrological scheme includes processes of saturation, nucleation, particle growth, sublimation and sedimentation under the assu…
▽ More
We present a new implementation of the hydrological cycle scheme into a general circulation model of the Martian atmosphere. The model includes a semi-Lagrangian transport scheme for water vapor and ice, and accounts for microphysics of phase transitions between them. The hydrological scheme includes processes of saturation, nucleation, particle growth, sublimation and sedimentation under the assumption of a variable size distribution. The scheme has been implemented into the Max Planck Institute Martian general circulation model (MPI--MGCM) and tested assuming mono- and bimodal log-normal distributions of ice condensation nuclei. We present a comparison of the simulated annual variations, horizontal and vertical distributions of water vapor and ice clouds with the available observations from instruments onboard Mars orbiters. The accounting for bi-modality of aerosol particle distribution improves the simulations of the annual hydrological cycle, including predicted ice clouds mass, opacity, number density, particle radii. The increased number density and lower nucleation rates brings the simulated cloud opacities closer to observations. Simulations show a weak effect of the excess of small aerosol particles on the simulated water vapor distributions.
△ Less
Submitted 13 January, 2022;
originally announced January 2022.
-
Martian Dust Storms and Gravity Waves: Disentangling Water Transport to the Upper Atmosphere
Authors:
Dmitry S. Shaposhnikov,
Alexander S. Medvedev,
Alexander V. Rodin,
Erdal Yiğit,
Paul Hartogh
Abstract:
Simulations with the Max Planck Institute Martian general circulation model for Martian years 28 and 34 reveal details of the water "pump" mechanism and the role of gravity wave (GW) forcing. Water is advected to the upper atmosphere mainly by upward branches of the meridional circulation: in low latitudes during equinoxes and over the south pole during solstices. Molecular diffusion plays little…
▽ More
Simulations with the Max Planck Institute Martian general circulation model for Martian years 28 and 34 reveal details of the water "pump" mechanism and the role of gravity wave (GW) forcing. Water is advected to the upper atmosphere mainly by upward branches of the meridional circulation: in low latitudes during equinoxes and over the south pole during solstices. Molecular diffusion plays little role in water transport in the middle atmosphere and across the mesopause. GWs modulate the circulation and temperature during global dust storms, thus changing the timing and intensity of the transport. At equinoxes, they facilitate water accumulation in the polar warming regions in the middle atmosphere followed by stronger upwelling over the equator. As equinoctial storms decay, GWs tend to accelerate the reduction of water in the thermosphere. GWs delay the onset of the transport during solstitial storms and change the globally averaged amount of water in the upper atmosphere by 10-25%.
△ Less
Submitted 13 January, 2022;
originally announced January 2022.
-
Exotic Magnetic and Electronic Properties of Layered CrI3 Single Crystals Under High Pressure
Authors:
Anirudha Ghosh,
D. Singh,
Qingge Mu,
Y. Kvashnin,
G. Haider,
M. Jonak,
D. Chareev,
T. Aramaki,
S. A. Medvedev,
R. Klingeler,
M. Mito,
E. H. Abdul-Hafidh,
J. Vejpravova,
M. Kalbac,
R. Ahuja,
Olle Eriksson,
Mahmoud Abdel-Hafiez
Abstract:
Through advanced experimental techniques on CrI$_{3}$ single crystals, we derive a previously not discussed pressure-temperature phase diagram. We find that $T_{c}$ increases to $\sim$ 66\,K with pressure up to $\sim$ 3\,GPa followed by a decrease to $\sim$ 10\,K at 21.2\,GPa. The experimental results are reproduced by theoretical calculations based on density functional theory where electron-elec…
▽ More
Through advanced experimental techniques on CrI$_{3}$ single crystals, we derive a previously not discussed pressure-temperature phase diagram. We find that $T_{c}$ increases to $\sim$ 66\,K with pressure up to $\sim$ 3\,GPa followed by a decrease to $\sim$ 10\,K at 21.2\,GPa. The experimental results are reproduced by theoretical calculations based on density functional theory where electron-electron interactions are treated by a static on-site Hubbard U on Cr 3$d$ orbitals. The origin of the pressure induced reduction of the ordering temperature is associated with a decrease of the calculated bond angle, from 95$^{\circ}$ at ambient pressure to $\sim$ 85$^{\circ}$ at 25\,GPa. Above 22\,GPa, the magnetically ordered state is essentially quenched, possibly driving the system to a Kitaev spin-liquid state at low temperature, thereby opening up the possibility of further exploration of long-range quantum entanglement between spins. The pressure-induced semiconductor-to-metal phase transition was revealed by high-pressure resistivity that is accompanied by a transition from a robust ferromagnetic state to gradually more dominating anti-ferromagnetic interactions and was consistent with theoretical modeling.
△ Less
Submitted 31 July, 2021;
originally announced August 2021.
-
Variations of the Martian Thermospheric Gravity Wave Activity during the Recent Solar Minimum as Observed by MAVEN
Authors:
Erdal Yiğit,
Alexander S. Medvedev,
Paul Hartogh
Abstract:
Atmospheric gravity (buoyancy) waves (GWs) are of great importance for the energy and momentum budget of all planetary atmospheres. Propagating upward waves carry energy and momentum from the lower atmosphere to thermospheric altitudes and re-distribute them there. On Mars, GWs dominate the variability of the thermosphere and ionosphere. We provide a comprehensive climatology of Martian thermosphe…
▽ More
Atmospheric gravity (buoyancy) waves (GWs) are of great importance for the energy and momentum budget of all planetary atmospheres. Propagating upward waves carry energy and momentum from the lower atmosphere to thermospheric altitudes and re-distribute them there. On Mars, GWs dominate the variability of the thermosphere and ionosphere. We provide a comprehensive climatology of Martian thermospheric GW activity at solar minimum (end of Solar Cycle 24) inferred from measurements by Neutral Gas and Ions Mass Spectrometer on board Mars Atmosphere and Volatile EvolutioN (NGIMS/MAVEN). The results are compared and interpreted using a one-dimensional spectral nonlinear GW model. Monthly mean GW activity varies strongly as a function of altitude (150-230 km) between 6-25%, reaching a maximum at $\sim$170 km. GW activity systematically exhibits a local time variability with nighttime values exceeding those during daytime, in accordance with previous studies. The analysis suggests that the day-night difference is primarily caused by a competition between dissipation due to molecular diffusion and wave growth due to decreasing background density. Thus, convective instability mechanism is likely to play a less important role in limiting GW amplitudes in the upper thermosphere, which explains their local time behavior.
△ Less
Submitted 22 July, 2021; v1 submitted 21 July, 2021;
originally announced July 2021.
-
Robust Weyl Semimetal and Pressure-induced Superconductivity in Quasi-One-Dimensional (NbSe4)2I
Authors:
Cuiying Pei,
Wujun Shi,
Yi Zhao,
Lingling Gao,
Jiacheng Gao,
Yiwei Li,
Haiyin Zhu,
Qing Zhang,
Na Yu,
Changhua Li,
Weizheng Cao,
Sergey A. Medvedev,
Claudia Felser,
Binghai Yan,
Zhongkai Liu,
Yulin Chen,
Zhijun Wang,
Yanpeng Qi
Abstract:
A search for the single material system that simultaneously exhibits topological phase and intrinsic superconductivity has been largely limited, although such a system is far more favorable especially for the quantum device applications. Except artificially engineered topological superconductivity in heterostructure systems, another alternative is to have superconductivity arising from the topolog…
▽ More
A search for the single material system that simultaneously exhibits topological phase and intrinsic superconductivity has been largely limited, although such a system is far more favorable especially for the quantum device applications. Except artificially engineered topological superconductivity in heterostructure systems, another alternative is to have superconductivity arising from the topological materials by pressure or other clean technology. Here, based on first-principles calculations, we first show that quasi-one-dimensional compound (NbSe4)2I represents a rare example of a chiral Weyl semimetal in which the set of symmetry-related Weyl points (WPs) exhibit the same chiral charge at a certain energy. The net chiral charge (NCC) of the below Fermi level EF (or a certain energy) can be tuned by pressure. In addition, a partial disorder induced by pressure accompanied with superconductivity emerges. Although amorphization of the iodine sub-lattice under high pressure, the one-dimensional NbSe4 chains in (NbSe4)2I remain intact and provide a superconducting channel in one dimension. Our combined theoretical and experimental research provide critical insight into a new phase of the one-dimensional system, in which distinctive phase transitions and correlated topological states emerge upon compression.
△ Less
Submitted 27 May, 2021;
originally announced May 2021.
-
Modelling railway delay propagation as diffusion-like spreading
Authors:
Mark M. Dekker,
Alexey N. Medvedev,
Jan Rombouts,
Grzegorz Siudem,
Liubov Tupikina
Abstract:
Railway systems form an important means of transport across the world. However, congestions or disruptions may significantly decrease these systems' efficiencies, making predicting and understanding the resulting train delays a priority for railway organisations. Delays are studied in a wide variety of models, which usually simulate trains as discrete agents carrying delays. In contrast, in this p…
▽ More
Railway systems form an important means of transport across the world. However, congestions or disruptions may significantly decrease these systems' efficiencies, making predicting and understanding the resulting train delays a priority for railway organisations. Delays are studied in a wide variety of models, which usually simulate trains as discrete agents carrying delays. In contrast, in this paper, we define a novel model for studying delays, where they spread across the railway network via a diffusion-like process. This type of modelling has various advantages such as quick computation and ease of applying various statistical tools like spectral methods, but it also comes with limitations related to the directional and discrete nature of delays and the trains carrying them. We apply the model to the Belgian railways and study its performance in simulating the delay propagation in severely disrupted railway situations. In particular, we discuss the role of spatial aggregation by proposing to cluster the Belgian railway system into sets of stations and adapt the model accordingly. We find that such aggregation significantly increases the model's performance. For some particular situations, a non-trivial optimal level of spatial resolution is found on which the model performs best. Our results show the potential of this type of delay modelling to understand large-scale properties of railway systems.
△ Less
Submitted 13 May, 2021;
originally announced May 2021.
-
XAMG: A library for solving linear systems with multiple right-hand side vectors
Authors:
Boris Krasnopolsky,
Alexey Medvedev
Abstract:
This paper presents the XAMG library for solving large sparse systems of linear algebraic equations with multiple right-hand side vectors. The library specializes but is not limited to the solution of linear systems obtained from the discretization of elliptic differential equations. A corresponding set of numerical methods includes Krylov subspace, algebraic multigrid, Jacobi, Gauss-Seidel, and C…
▽ More
This paper presents the XAMG library for solving large sparse systems of linear algebraic equations with multiple right-hand side vectors. The library specializes but is not limited to the solution of linear systems obtained from the discretization of elliptic differential equations. A corresponding set of numerical methods includes Krylov subspace, algebraic multigrid, Jacobi, Gauss-Seidel, and Chebyshev iterative methods. The parallelization is implemented with MPI+POSIX shared memory hybrid programming model, which introduces a three-level hierarchical decomposition using the corresponding per-level synchronization and communication primitives. The code contains a number of optimizations, including the multilevel data segmentation, compression of indices, mixed-precision floating-point calculations, vector status flags, and others. The XAMG library uses the program code of the well-known hypre library to construct the multigrid matrix hierarchy. The XAMG's own implementation for the solve phase of the iterative methods provides up to a twofold speedup compared to hypre for the tests performed. Additionally, XAMG provides extended functionality to solve systems with multiple right-hand side vectors.
△ Less
Submitted 12 March, 2021;
originally announced March 2021.
-
Dust storm-enhanced gravity wave activity in the Martian thermosphere observed by MAVEN and implication for atmospheric escape
Authors:
Erdal Yiğit,
Alexander S. Medvedev,
Mehdi Benna,
Bruce Jakosky
Abstract:
Lower atmospheric global dust storms affect the small- and large-scale weather and variability of the whole Martian atmosphere. Analysis of the CO$_2$ density data from the Neutral Gas and Ion Mass Spectrometer instrument (NGIMS) on board NASA's Mars Atmosphere Volatile EvolutioN (MAVEN) spacecraft show a remarkable increase of GW-induced density fluctuations in the thermosphere during the 2018 ma…
▽ More
Lower atmospheric global dust storms affect the small- and large-scale weather and variability of the whole Martian atmosphere. Analysis of the CO$_2$ density data from the Neutral Gas and Ion Mass Spectrometer instrument (NGIMS) on board NASA's Mars Atmosphere Volatile EvolutioN (MAVEN) spacecraft show a remarkable increase of GW-induced density fluctuations in the thermosphere during the 2018 major dust storm with distinct latitude and local time variability. The mean thermospheric GW activity increases by a factor of two during the storm event. The magnitude of relative density perturbations is around 20% on average and 40% locally. One and a half months later, the GW activity gradually decreases. Enhanced temperature disturbances in the Martian thermosphere can facilitate atmospheric escape. For the first time, we estimate that, for a 20% and 40% GW-induced disturbances, the net increase of Jeans escape flux of hydrogen is a factor of 1.3 and 2, respectively.
△ Less
Submitted 19 January, 2021;
originally announced January 2021.
-
Effects of latitude-dependent gravity wave source variations on the middle and upper atmosphere
Authors:
Erdal Yiğit,
Alexander S. Medvedev,
Manfred Ern
Abstract:
Atmospheric gravity waves (GWs) are generated in the lower atmosphere by various weather phenomena. They propagate upward, carry energy and momentum to higher altitudes, and appreciably influence the general circulation upon depositing them in the middle and upper atmosphere. We use a three-dimensional first-principle general circulation model (GCM) with an implemented nonlinear whole atmosphere G…
▽ More
Atmospheric gravity waves (GWs) are generated in the lower atmosphere by various weather phenomena. They propagate upward, carry energy and momentum to higher altitudes, and appreciably influence the general circulation upon depositing them in the middle and upper atmosphere. We use a three-dimensional first-principle general circulation model (GCM) with an implemented nonlinear whole atmosphere GW parameterization to study the global climatology of wave activity and produced effects at altitudes up to the upper thermosphere. The numerical experiments were guided by the GW momentum fluxes and temperature variances as measured in 2010 by the SABER (Sounding of the Atmosphere using Broadband Emission Radiometry) instrument onboard NASA's TIMED (Thermosphere Ionosphere Mesosphere Energetics Dynamics) satellite. This includes the latitudinal dependence and magnitude of GW activity in the lower stratosphere for the boreal summer season. The modeling results were compared to the SABER temperature and total absolute momentum flux, and Upper Atmosphere Research Satellite (UARS) data in the mesosphere and lower thermosphere. Simulations suggest that, in order to reproduce the observed circulation and wave activity in the middle atmosphere, smaller than the measured GW fluxes have to be used at the source level in the lower atmosphere. This is because observations contain a broader spectrum of GWs, while parameterizations capture only a portion relevant to the middle and upper atmosphere dynamics. Accounting for the latitudinal variations of the source appreciably improves simulations.
△ Less
Submitted 23 December, 2020;
originally announced December 2020.
-
Laser scattering by submicron droplets originated during the electrical explosion of thin metal wires
Authors:
V. M. Romanova,
G. V. Ivanenkov,
E. V. Parkevich,
I. N. Tilikin,
M. A. Medvedev,
T. A. Shelkovenko,
S. A. Pikuz,
A. S. Selyukov
Abstract:
This paper presents the results of studying of dispersed media formation during the electrical explosion of thin metal wires in vacuum by using low-current generators ($\sim 1$-$10$~kA). Particular attention is paid to the analysis of the composition and structure of the corresponding explosion products as well as to the problem of their visualization using simultaneous laser interferometry and sh…
▽ More
This paper presents the results of studying of dispersed media formation during the electrical explosion of thin metal wires in vacuum by using low-current generators ($\sim 1$-$10$~kA). Particular attention is paid to the analysis of the composition and structure of the corresponding explosion products as well as to the problem of their visualization using simultaneous laser interferometry and shadow imaging at two wavelengths (1.064~$μ$m and 0.532~$μ$m). Our findings point to the fact that the important role in the visualization of the explosion products belongs to multiple scattering by submicron droplets of dense condensed matter, which are mixed with metal vapor. The hypothesis on the existence of submicron droplets in the products of exploding metal wires correlates with the results obtained by soft x-ray radiography combined with a laser probing technique. Taking into account the multiple scattering by submicron droplets, it is possible to significantly clarify the parameters of the explosion products visualized via laser probing techniques as well as to gain a deeper insight into the physics behind the electrical wire explosion.
△ Less
Submitted 3 December, 2020; v1 submitted 27 October, 2020;
originally announced November 2020.
-
NGC 5474 X-1: a neutron star ULX in an old stellar cluster?
Authors:
Kirill Atapin,
Alexander Vinokurov,
Arkadiy Sarkisyan,
Yulia Solovyeva,
Alexander Kostenkov,
Aleksei Medvedev,
Sergei Fabrika
Abstract:
We present the optical and X-ray study of the ultraluminous X-ray source (ULX) NGC 5474 X-1. The X-ray spectrum taken during the bright state of the source ($L_X \sim 2\times10^{40}$ egs/s) shows signatures of a broad absorption line at $\simeq 8$ keV which may be a cyclotron resonant scattering feature. This implies that this system may host a neutron star with a magnetic field $\sim 10^{12}$ G.…
▽ More
We present the optical and X-ray study of the ultraluminous X-ray source (ULX) NGC 5474 X-1. The X-ray spectrum taken during the bright state of the source ($L_X \sim 2\times10^{40}$ egs/s) shows signatures of a broad absorption line at $\simeq 8$ keV which may be a cyclotron resonant scattering feature. This implies that this system may host a neutron star with a magnetic field $\sim 10^{12}$ G. The first observation of this area with the Hubble Space Telescope (HST) carried out 14 months later revealed that the source was bright in the optical range as well. The subsequent observations have shown that the source faded in both ranges (more than $2.8^m$ in the U band and by a factor of 50-100 in X-rays) and has never become bright again. Deeper HST observations made it possible to impose constraints on the donor star spectral class and mass ($<\,7$ M$_\odot$), as well as to identify a stellar cluster of about 1 Gyr, the centre of which is located at a projected distance of $\simeq 2$ pc from NGC 5474 X-1. The ULX can be a member of this old cluster, however, the presence of stars with ages of $\sim10$ Myr within 300 pc around the ULX does not completely exclude the possibility that it is just an accidental projection.
△ Less
Submitted 29 November, 2023; v1 submitted 2 November, 2020;
originally announced November 2020.
-
Coexistence of Superconductivity and Charge Density Wave in Tantalum Disulfide: Experiment and Theory
Authors:
Y. Kvashnin,
D. VanGennep,
M. Mito,
S. A. Medvedev,
R. Thiyagarajan,
O. Karis,
A. N. Vasiliev,
O. Eriksson,
M. Abdel-Hafiez
Abstract:
The coexistence of charge density wave (CDW) and superconductivity in tantalum disulfide (2H-TaS$_2$) at ambient pressure, is boosted by applying hydrostatic pressures up to 30GPa, thereby inducing a typical dome-shaped superconducting phase. The ambient pressure CDW ground state which begins at TCDW = 76 K, with critically small Fermi surfaces, was found to be fully suppressed at Pc = 8.7GPa. Aro…
▽ More
The coexistence of charge density wave (CDW) and superconductivity in tantalum disulfide (2H-TaS$_2$) at ambient pressure, is boosted by applying hydrostatic pressures up to 30GPa, thereby inducing a typical dome-shaped superconducting phase. The ambient pressure CDW ground state which begins at TCDW = 76 K, with critically small Fermi surfaces, was found to be fully suppressed at Pc = 8.7GPa. Around Pc, we observe a superconducting dome with a maximum superconducting transition temperature Tc = 9.1 K. First-principles calculations of the electronic structure predict that, under ambient conditions, the undistorted structure is characterized by a phonon instability at finite momentum close to the experimental CDW wave vector. Upon compression, this instability is found to disappear, indicating the suppression of CDW order. The calculations reveal an electronic topological transition (ETT), which occurs before the suppression of the phonon instability, suggesting that the ETT alone is not directly causing the structural change in the system. The temperature dependence of the first vortex penetration field has been experimentally obtained by two independent methods and the corresponding lower critical field H$_{c1}$ was deduced. While a d wave and single-gap BCS prediction cannot describe our H$_{c1}$ experiments, the temperature dependence of the H$_{c1}$ can be well described by a single-gap anisotropic s-wave order parameter.
△ Less
Submitted 29 July, 2020;
originally announced July 2020.
-
Initialization of a Disease Transmission Model
Authors:
Håkan Runvik,
Alexander Medvedev,
Robin Eriksson,
Stefan Engblom
Abstract:
Approaches to the calculation of the full state vector of a larger epidemiological model for the spread of COVID-19 in Sweden at the initial time instant from available data and with a simplified dynamical model are proposed and evaluated. The larger epidemiological model is based on a continuous Markov chain and captures the demographic composition of and the transport flows between the counties…
▽ More
Approaches to the calculation of the full state vector of a larger epidemiological model for the spread of COVID-19 in Sweden at the initial time instant from available data and with a simplified dynamical model are proposed and evaluated. The larger epidemiological model is based on a continuous Markov chain and captures the demographic composition of and the transport flows between the counties of Sweden. Its intended use is to predict the outbreak development in temporal and spatial coordinates as well as across the demographic groups. It can also support evaluating and comparing of prospective intervention strategies in terms of e.g. lockdown in certain areas or isolation of specific age groups. The simplified model is a discrete time-invariant linear system that has cumulative infectious incidence, infected population, asymptomatic population, exposed population, and infectious pressure as the state variables. Since the system matrix of the model depends on a number transition rates, structural properties of the model are investigated for suitable parameter ranges. It is concluded that the model becomes unobservable for some parameter values. Two contrasting approaches to the initial state estimation are considered. One is a version of Rauch-Tung-Striebel smoother and another is based on solving a batch nonlinear optimization problem. The benefits and shortcomings of the considered estimation techniques are analyzed and compared on synthetic data for several Swedish counties.
△ Less
Submitted 17 July, 2020;
originally announced July 2020.
-
Cartan connections for stochastic developments on sub-Riemannian manifolds
Authors:
Ivan Beschastnyi,
Karen Habermann,
Alexandr Medvedev
Abstract:
Analogous to the characterisation of Brownian motion on a Riemannian manifold as the development of Brownian motion on a Euclidean space, we construct sub-Riemannian diffusions on equinilpotentisable sub-Riemannian manifolds by developing a canonical stochastic process arising as the lift of Brownian motion to an associated model space. The notion of stochastic development we introduce for equinil…
▽ More
Analogous to the characterisation of Brownian motion on a Riemannian manifold as the development of Brownian motion on a Euclidean space, we construct sub-Riemannian diffusions on equinilpotentisable sub-Riemannian manifolds by developing a canonical stochastic process arising as the lift of Brownian motion to an associated model space. The notion of stochastic development we introduce for equinilpotentisable sub-Riemannian manifolds uses Cartan connections, which take the place of the Levi-Civita connection in Riemannian geometry. We first derive a general expression for the generator of the stochastic process which is the stochastic development with respect to a Cartan connection of the lift of Brownian motion to the model space. We further provide a necessary and sufficient condition for the existence of a Cartan connection which develops the canonical stochastic process to the sub-Riemannian diffusion associated with the sub-Laplacian defined with respect to the Popp volume. We illustrate the construction of a suitable Cartan connection for free sub-Riemannian structures with two generators and we discuss an example where the condition is not satisfied.
△ Less
Submitted 7 December, 2020; v1 submitted 29 June, 2020;
originally announced June 2020.
-
Simplex2Vec embeddings for community detection in simplicial complexes
Authors:
Jacob Charles Wright Billings,
Mirko Hu,
Giulia Lerda,
Alexey N. Medvedev,
Francesco Mottes,
Adrian Onicas,
Andrea Santoro,
Giovanni Petri
Abstract:
Topological representations are rapidly becoming a popular way to capture and encode higher-order interactions in complex systems. They have found applications in disciplines as different as cancer genomics, brain function, and computational social science, in representing both descriptive features of data and inference models. While intense research has focused on the connectivity and homological…
▽ More
Topological representations are rapidly becoming a popular way to capture and encode higher-order interactions in complex systems. They have found applications in disciplines as different as cancer genomics, brain function, and computational social science, in representing both descriptive features of data and inference models. While intense research has focused on the connectivity and homological features of topological representations, surprisingly scarce attention has been given to the investigation of the community structures of simplicial complexes. To this end, we adopt recent advances in symbolic embeddings to compute and visualize the community structures of simplicial complexes. We first investigate the stability properties of embedding obtained for synthetic simplicial complexes to the presence of higher order interactions. We then focus on complexes arising from social and brain functional data and show how higher order interactions can be leveraged to improve clustering detection and assess the effect of higher order interaction on individual nodes. We conclude delineating limitations and directions for extension of this work.
△ Less
Submitted 21 June, 2019;
originally announced June 2019.
-
Signal detection algorithms for interferometric sensors with harmonic phase modulation: distortion analysis and suppression
Authors:
Leonid Liokumovich,
Andrei Medvedev,
Konstantin Muravyov,
Philipp Skliarov,
Nikolai Ushakov
Abstract:
In current paper, distortions in digital demodulation schemes with harmonic phase modulation for interferometric optical sensors are considered. In particular, the influence of target signal variations on phase demodulation errors is theoretically evaluated. An analytical expression, describing the phase error magnitude dependence on first derivative and mean value of measured signal and amplitude…
▽ More
In current paper, distortions in digital demodulation schemes with harmonic phase modulation for interferometric optical sensors are considered. In particular, the influence of target signal variations on phase demodulation errors is theoretically evaluated. An analytical expression, describing the phase error magnitude dependence on first derivative and mean value of measured signal and amplitude of the phase modulation in case of simple 4-point demodulation algorithm is derived. After that, an approach for synthesizing algorithms with suppressed sensitivity to target signal variations is developed. Based on this approach, a novel 4+1 demodulation algorithm is proposed. It is shown analytically that the demodulation error of new 4+1 algorithm is proportional to the second derivative of target signal, and therefore, is typically several orders of magnitude smaller than in case of 4-point algorithm. The correspondence between analytical expressions and real phase errors, induced by target signal variations is verified by means of numeric simulation.
△ Less
Submitted 30 April, 2019;
originally announced April 2019.
-
Seasonal Water "Pump" in the Atmosphere of Mars: Vertical Transport to the Thermosphere
Authors:
Dmitry S. Shaposhnikov,
Alexander S. Medvedev,
Alexander V. Rodin,
Paul Hartogh
Abstract:
We present results of simulations with the Max Planck Institute general circulation model (MPI-MGCM) implementing a hydrological cycle scheme. The simulations reveal a seasonal water "pump" mechanism responsible for the upward transport of water vapor. This mechanism occurs in high latitudes above 60$^\circ$ of the southern hemisphere at perihelion, when the upward branch of the meridional circula…
▽ More
We present results of simulations with the Max Planck Institute general circulation model (MPI-MGCM) implementing a hydrological cycle scheme. The simulations reveal a seasonal water "pump" mechanism responsible for the upward transport of water vapor. This mechanism occurs in high latitudes above 60$^\circ$ of the southern hemisphere at perihelion, when the upward branch of the meridional circulation is particularly strong. A combination of the mean vertical flux with variations induced by solar tides facilitates penetration of water across the "bottleneck" at approximately 60 km. The meridional circulation then transports water across the globe to the northern hemisphere. Since the intensity of the meridional cell is tightly controlled by airborne dust, the water abundance in the thermosphere strongly increases during dust storms.
△ Less
Submitted 12 April, 2019;
originally announced April 2019.
-
Small cycles, generalized prisms and Hamiltonian cycles in the Bubble-sort graph
Authors:
Elena V. Konstantinova,
Alexey N. Medvedev
Abstract:
The Bubble-sort graph $BS_n,\,n\geqslant 2$, is a Cayley graph over the symmetric group $Sym_n$ generated by transpositions from the set $\{(1 2), (2 3),\ldots, (n-1 n)\}$. It is a bipartite graph containing all even cycles of length $\ell$, where $4\leqslant \ell\leqslant n!$. We give an explicit combinatorial characterization of all its $4$- and $6$-cycles. Based on this characterization, we def…
▽ More
The Bubble-sort graph $BS_n,\,n\geqslant 2$, is a Cayley graph over the symmetric group $Sym_n$ generated by transpositions from the set $\{(1 2), (2 3),\ldots, (n-1 n)\}$. It is a bipartite graph containing all even cycles of length $\ell$, where $4\leqslant \ell\leqslant n!$. We give an explicit combinatorial characterization of all its $4$- and $6$-cycles. Based on this characterization, we define generalized prisms in $BS_n,\,n\geqslant 5$, and present a new approach to construct a Hamiltonian cycle based on these generalized prisms.
△ Less
Submitted 4 April, 2021; v1 submitted 12 January, 2019;
originally announced January 2019.
-
Variations on the Feferman-Vaught Theorem, with applications to $\prod_p \mathbb{F}_p$
Authors:
Alice Medvedev,
Alexander Van Abel
Abstract:
Using the Feferman-Vaught Theorem, we prove that a definable subset of a product structure must be a Boolean combination of open sets, in the product topology induced by giving each factor structure the discrete topology. We prove a converse of the Feferman-Vaught theorem for families of structures with certain properties, including families of integral domains. We use these results to obtain char…
▽ More
Using the Feferman-Vaught Theorem, we prove that a definable subset of a product structure must be a Boolean combination of open sets, in the product topology induced by giving each factor structure the discrete topology. We prove a converse of the Feferman-Vaught theorem for families of structures with certain properties, including families of integral domains. We use these results to obtain characterizations of the definable subsets of $\prod_p \mathbb{F}_p$ -- in particular, every formula is equivalent to a Boolean combination of $\exists \forall \exists$ formulae.
△ Less
Submitted 6 December, 2018;
originally announced December 2018.
-
Influence of gravity waves on the climatology of high-altitude Martian carbon dioxide ice clouds
Authors:
Erdal Yiğit,
Alexander S. Medvedev,
Paul Hartogh
Abstract:
Carbon dioxide (CO$_2$) ice clouds have been routinely observed in the middle atmosphere of Mars. However, there are still uncertainties concerning physical mechanisms that control their altitude, geographical, and seasonal distributions. Using the Max Planck Institute Martian General Circulation Model (MPI-MGCM), incorporating a state-of-the-art whole atmosphere subgrid-scale gravity wave paramet…
▽ More
Carbon dioxide (CO$_2$) ice clouds have been routinely observed in the middle atmosphere of Mars. However, there are still uncertainties concerning physical mechanisms that control their altitude, geographical, and seasonal distributions. Using the Max Planck Institute Martian General Circulation Model (MPI-MGCM), incorporating a state-of-the-art whole atmosphere subgrid-scale gravity wave parameterization [Yiğit et al., 2008], we demonstrate that internal gravity waves generated by lower atmospheric weather processes have wide reaching impact on the Martian climate. Globally, GWs cool the upper atmosphere of Mars by $\sim$10 % and facilitate high-altitude CO$_2$ ice cloud formation. CO$_2$ ice cloud seasonal variations in the mesosphere and the mesopause region appreciably coincide with the spatio-temporal variations of GW effects, providing insight into the observed distribution of clouds. Our results suggest that GW propagation and dissipation constitute a necessary physical mechanism for CO$_2$ ice cloud formation in the Martian upper atmosphere during all seasons.
△ Less
Submitted 2 December, 2018;
originally announced December 2018.
-
The anatomy of Reddit: An overview of academic research
Authors:
Alexey N. Medvedev,
Renaud Lambiotte,
Jean-Charles Delvenne
Abstract:
Online forums provide rich environments where users may post questions and comments about different topics. Understanding how people behave in online forums may shed light on the fundamental mechanisms by which collective thinking emerges in a group of individuals, but it has also important practical applications, for instance to improve user experience, increase engagement or automatically identi…
▽ More
Online forums provide rich environments where users may post questions and comments about different topics. Understanding how people behave in online forums may shed light on the fundamental mechanisms by which collective thinking emerges in a group of individuals, but it has also important practical applications, for instance to improve user experience, increase engagement or automatically identify bullying. Importantly, the datasets generated by the activity of the users are often openly available for researchers, in contrast to other sources of data in computational social science. In this survey, we map the main research directions that arose in recent years and focus primarily on the most popular platform, Reddit. We distinguish and categorise research depending on their focus on the posts or on the users, and point to different types of methodologies to extract information from the structure and dynamics of the system. We emphasize the diversity and richness of the research in terms of questions and methods, and suggest future avenues of research.
△ Less
Submitted 4 June, 2020; v1 submitted 25 October, 2018;
originally announced October 2018.
-
Identifying exogenous and endogenous activity in social media
Authors:
Kazuki Fujita,
Alexey Medvedev,
Shinsuke Koyama,
Renaud Lambiotte,
Shigeru Shinomoto
Abstract:
The occurrence of new events in a system is typically driven by external causes and by previous events taking place inside the system. This is a general statement, applying to a range of situations including, more recently, to the activity of users in Online social networks (OSNs). Here we develop a method for extracting from a series of posting times the relative contributions of exogenous, e.g.…
▽ More
The occurrence of new events in a system is typically driven by external causes and by previous events taking place inside the system. This is a general statement, applying to a range of situations including, more recently, to the activity of users in Online social networks (OSNs). Here we develop a method for extracting from a series of posting times the relative contributions of exogenous, e.g. news media, and endogenous, e.g. information cascade. The method is based on the fitting of a generalized linear model (GLM) equipped with a self-excitation mechanism. We test the method with synthetic data generated by a nonlinear Hawkes process, and apply it to a real time series of tweets with a given hashtag. In the empirical dataset, the estimated contributions of exogenous and endogenous volumes are close to the amounts of original tweets and retweets respectively. We conclude by discussing the possible applications of the method, for instance in online marketing.
△ Less
Submitted 2 August, 2018;
originally announced August 2018.
-
Classification of homogeneous affine tube domains with large automorphism groups in arbitrary dimensions
Authors:
Vladimir Ezhov,
Alexandr Medvedev,
Gerd Schmalz
Abstract:
We classify tube domains in $\mathbb C^{n+1}$ ($n\ge 1$) with affinely homogeneous base of their boundary and a.) with positive definite Levi form and b.) with Lorentzian type Levi form and affine isotropy of dimension at least $\frac{(n-2)(n-3)}2$.
We classify tube domains in $\mathbb C^{n+1}$ ($n\ge 1$) with affinely homogeneous base of their boundary and a.) with positive definite Levi form and b.) with Lorentzian type Levi form and affine isotropy of dimension at least $\frac{(n-2)(n-3)}2$.
△ Less
Submitted 6 April, 2018;
originally announced April 2018.
-
The peculiarities of near-cathode processes in air discharge at atmospheric pressure
Authors:
E V Parkevich,
M A Medvedev,
A V Agafonov,
S I Tkachenko,
A V Oginov,
A I Khirianova,
A R Mingaleev,
T A Shelkovenko,
S A Pikuz
Abstract:
Formation of near-cathode plasma at the instant of breakdown of the air gap was studied by the methods of picosecond laser probing. It was demonstrated that 1-2 ns after a sharp rise of the current through the discharge gap dense plasma clots with Ne~10^20 cm^-3 and dNe/dx ~ 10^24 cm^-4 are formed on the cathode surface. It was shown that these highly ionized regions lead to initiation and develop…
▽ More
Formation of near-cathode plasma at the instant of breakdown of the air gap was studied by the methods of picosecond laser probing. It was demonstrated that 1-2 ns after a sharp rise of the current through the discharge gap dense plasma clots with Ne~10^20 cm^-3 and dNe/dx ~ 10^24 cm^-4 are formed on the cathode surface. It was shown that these highly ionized regions lead to initiation and development of the spark channel originating from the cathode. We propose that the observed formations correspond to erosive plasma formed from the cathode material.
△ Less
Submitted 4 April, 2018;
originally announced April 2018.
-
Modelling structure and predicting dynamics of discussion threads in online boards
Authors:
Alexey N. Medvedev,
Jean-Charles Delvenne,
Renaud Lambiotte
Abstract:
Internet boards are platforms for online discussions about a variety of topics. On these boards, individuals may start a new thread on a specific matter, or leave comments in an existing discussion. The resulting collective process leads to the formation of `discussion trees', where nodes represent a post and comments, and an edge represents a `reply-to' relation. The structure of discussion trees…
▽ More
Internet boards are platforms for online discussions about a variety of topics. On these boards, individuals may start a new thread on a specific matter, or leave comments in an existing discussion. The resulting collective process leads to the formation of `discussion trees', where nodes represent a post and comments, and an edge represents a `reply-to' relation. The structure of discussion trees has been analysed in previous works, but only from a static perspective. In this paper, we focus on their structural and dynamical properties by modelling their formation as a self-exciting Hawkes process. We first study a Reddit dataset to show that the structure of the trees resemble those produced by a Galton-Watson process with a special root offspring distribution. The dynamical aspect of the model is then used to predict future commenting activity and the final size of a discussion tree. We compare the efficiency of our approach with previous works and show its superiority for the prediction of the dynamics of discussions.
△ Less
Submitted 4 June, 2020; v1 submitted 30 January, 2018;
originally announced January 2018.
-
Constant Curvature Models in Sub-Riemannian Geometry
Authors:
D. Alekseevsky,
A. Medvedev,
J. Slovak
Abstract:
Each sub-Riemannian geometry with bracket generating distribution enjoys a background structure determined by the distribution itself. At the same time, those geometries with constant sub-Riemannian symbols determine a unique Cartan connection leading to their principal invariants. We provide cohomological description of the structure of these curvature invariants in the cases where the background…
▽ More
Each sub-Riemannian geometry with bracket generating distribution enjoys a background structure determined by the distribution itself. At the same time, those geometries with constant sub-Riemannian symbols determine a unique Cartan connection leading to their principal invariants. We provide cohomological description of the structure of these curvature invariants in the cases where the background structure is one of the parabolic geometries. As an illustration, constant curvature models are discussed for certain sub-Riemannian geometries.
△ Less
Submitted 4 October, 2018; v1 submitted 29 December, 2017;
originally announced December 2017.
-
Homogeneous Levi non-degenerate hypersurfaces in $\mathbb{C}^3$
Authors:
Boris Doubrov,
Alexandr Medvedev,
Dennis The
Abstract:
We classify all (locally) homogeneous Levi non-degenerate real hypersurfaces in $\mathbb{C}^3$ with symmetry algebra of dimension $\geq 6$.
We classify all (locally) homogeneous Levi non-degenerate real hypersurfaces in $\mathbb{C}^3$ with symmetry algebra of dimension $\geq 6$.
△ Less
Submitted 7 November, 2017;
originally announced November 2017.
-
Pressure-induced Lifshitz transition in NbP: Raman, x-ray diffraction, electrical transport and density functional theory
Authors:
Satyendra Nath Gupta,
Anjali Singh,
Koushik Pal,
D. V. S. Muthu,
C. Shekhar,
Yanpeng Qi,
Pavel G. Naumov,
Sergey A. Medvedev,
C. Felser,
U. V. Waghmare,
A. K. Sood
Abstract:
We report high pressure Raman, synchrotron x-ray diffraction and electrical transport studies on Weyl semimetals NbP and TaP along with first-principles density functional theoretical (DFT) analysis. The frequencies of first-order Raman modes of NbP harden with increasing pressure and exhibit a slope change at P$_c$ $\sim$ 9 GPa, and its resistivity exhibits a minimum at P$_c$. The pressure-depend…
▽ More
We report high pressure Raman, synchrotron x-ray diffraction and electrical transport studies on Weyl semimetals NbP and TaP along with first-principles density functional theoretical (DFT) analysis. The frequencies of first-order Raman modes of NbP harden with increasing pressure and exhibit a slope change at P$_c$ $\sim$ 9 GPa, and its resistivity exhibits a minimum at P$_c$. The pressure-dependent volume of NbP exhibits a change in its bulk modulus from 207 GPa to 243 GPa at P$_c$. Using DFT calculations, we show that these anomalies are associated with pressure induced Lifshitz transition which involves appearance of electron and hole pockets in its electronic structure. In contrast, results of Raman and synchrotron x-ray diffraction experiments on TaP and DFT calculations show that TaP is quite robust under pressure and does not undergo any phase transition.
△ Less
Submitted 27 September, 2017;
originally announced September 2017.