Obfuscation resilient Android Malware Analysis Based on Contrastive Learning. Due to its open source nature, Android operating system has been the main target of attackers to exploit. Malware creators always perform different code obfuscations on their apps to hide malicious activities. Features extracted from these obfuscated samples through program analysis contain many useless and disguised features, which leads to many false negatives. To address the issue, in this paper, we demonstrate that obfuscation resilient malware analysis can be achieved through contrastive learning. We take the Android malware classification as an example to demonstrate our analysis. The key insight behind our analysis is that contrastive learning can be used to reduce the difference introduced by obfuscation while amplifying the difference between malware and benign apps (or other types of malware). Based on the proposed analysis, we design a system that can achieve robust and interpretable classification of Android malware. To achieve robust classification, we perform contrastive learning on malware samples to learn an encoder that can automatically extract robust features from malware samples. To achieve interpretable classification, we transform the function call graph of a sample into an image by centrality analysis. Then the corresponding heatmaps are obtained by visualization techniques. These heatmaps can help users understand why the malware is classified as this family. We implement IFDroid and perform extensive evaluations on two widely used datasets. Experimental results show that IFDroid is superior to state of the art Android malware familial classification systems. Moreover, IFDroid is capable of maintaining 98.2 true positive rate on classifying 8,112 obfuscated malware samples.


On the Use and Misuse of Absorbing States in Multi agent Reinforcement Learning. The creation and destruction of agents in cooperative multi agent reinforcement learning (MARL) is a critically under explored area of research. Current MARL algorithms often assume that the number of agents within a group remains fixed throughout an experiment. However, in many practical problems, an agent may terminate before their teammates. This early termination issue presents a challenge the terminated agent must learn from the group's success or failure which occurs beyond its own existence. We refer to propagating value from rewards earned by remaining teammates to terminated agents as the Posthumous Credit Assignment problem. Current MARL methods handle this problem by placing these agents in an absorbing state until the entire group of agents reaches a termination condition. Although absorbing states enable existing algorithms and APIs to handle terminated agents without modification, practical training efficiency and resource use problems exist. In this work, we first demonstrate that sample complexity increases with the quantity of absorbing states in a toy supervised learning task for a fully connected network, while attention is more robust to variable size input. Then, we present a novel architecture for an existing state of the art MARL algorithm which uses attention instead of a fully connected layer with absorbing states. Finally, we demonstrate that this novel architecture significantly outperforms the standard architecture on tasks in which agents are created or destroyed within episodes as well as standard multi agent coordination tasks.


Modelling the Dissipation Range of von K 'arm 'an Turbulence Spectrum. Noise pollution caused by inflow turbulence is a major problem in many applications, such as propellers and fans. Leading edge noise models, e.g., Amiet's model, are widely applied to predict the noise produced by these applications. This noise prediction model relies on the accuracy of the turbulence spectrum, which is usually assumed to be the von K 'arm 'an energy spectrum for isotropic turbulence. However, the von K 'arm 'an spectrum does not model accurately the dissipation range of the turbulent energy, resulting in incorrect far field noise predictions for the high frequency range. An exponential correction can be applied to the spectrum to model the dissipation range. This correction depends on the dissipation frequency, which is the frequency where the energy spectrum changes from a dependence of k 5 3 to an exponential decay. This study experimentally investigates nearly isotropic inflow turbulence and determines the flow field characteristics that affect the dissipation frequency in order to model this frequency and the dissipation range. Experiments have been conducted with two passive grids and hot wire anemometry in the Aeroacoustic Wind Tunnel of the University of Twente. The turbulence uniformity in the test section and the turbulence development in the streamwise direction were analyzed, showing that the grid generated turbulence was mostly uniform and nearly isotropic. The dissipation frequency was observed to depend on the turbulence intensity, the free stream velocity, and the turbulence length scale. An expression to compute this frequency is proposed, as well as a formula to predict the dissipation range. The predicted leading edge noise is affected by the dissipation range modelling in the high frequency range, presenting a decrease in level up to 17 dB for the highest frequencies.


A note on the sum of finite multiple harmonic q series on r text (r 1) indices. We study the sum of the finite multiple harmonic q series on r text (r 1) indices at roots of unity with r 1,2,3 . And we give the equivalent conditions of two conjectures regarding cyclic sums of finite multiple harmonic q series on 1 text 2 text 3 indices at roots of unity, posed recently by Kh. Pilehrood, T. Pilehrood and R. Tauraso.


A Meta embedding based Ensemble Approach for ICD Coding Prediction. International Classification of Diseases (ICD) are the de facto codes used globally for clinical coding. These codes enable healthcare providers to claim reimbursement and facilitate efficient storage and retrieval of diagnostic information. The problem of automatically assigning ICD codes has been approached in literature as a multilabel classification, using neural models on unstructured data. Our proposed approach enhances the performance of neural models by effectively training word vectors using routine medical data as well as external knowledge from scientific articles. Furthermore, we exploit the geometric properties of the two sets of word vectors and combine them into a common dimensional space, using meta embedding techniques. We demonstrate the efficacy of this approach for a multimodal setting, using unstructured and structured information. We empirically show that our approach improves the current state of the art deep learning architectures and benefits ensemble models.


Convergence rate for a class of supercritical superprocesses. Suppose X X t, t ge 0 is a supercritical superprocess. Let phi be the non negative eigenfunction of the mean semigroup of X corresponding to the principal eigenvalue lambda 0 . Then M t( phi) e lambda t langle phi, X t rangle, t geq 0, is a non negative martingale with almost sure limit M infty( phi) . In this paper we study the rate at which M t( phi) M infty( phi) converges to 0 as t to infty when the process may not have finite variance. Under some conditions on the mean semigroup, we provide sufficient and necessary conditions for the rate in the almost sure sense. Some results on the convergence rate in L p with p in(1, 2) are also obtained.


Orbital varieties in types B and C . We correct the proof of the main result of an earlier paper, parametrizing orbital varieties in a complex simple Lie algebra of type B or C in terms of domino tableaux and showing how to compute the orbital variety attached to an element of the Weyl group in either of these types.


Recovering the Potential and Order in One Dimensional Time Fractional Diffusion with Unknown Initial Condition and Source. This paper is concerned with an inverse problem of recovering a potential term and fractional order in a one dimensional subdiffusion problem, which involves a Djrbashian Caputo fractional derivative of order alpha in(0,1) in time, from the lateral Cauchy data. In the model, we do not assume a full knowledge of the initial data and the source term, since they might be unavailable in some practical applications. We prove the unique recovery of the spatially dependent potential coefficient and the order alpha of the derivation simultaneously from the measured trace data at one end point, when the model is equipped with a boundary excitation with a compact support away from t 0 . One of the initial data and the source can also be uniquely determined, provided that the other is known. The analysis employs a representation of the solution and the time analyticity of the associated function. Further, we discuss a two stage procedure, directly inspired by the analysis, for the numerical identification of the order and potential coefficient, and illustrate the feasibility of the recovery with several numerical experiments.


Temporal Nuances of Coordination Networks. Current network based methods for detecting coordinated inauthentic behaviour on social media focus primarily on inferring links between accounts based on "behavioural traces" 1 , such as retweeting the same tweet or posting the same URL. Assuming the goal of coordination is amplification, i.e., boosting a message within a constrained period, most approaches use a temporal window to ensure the co activity occurs within a specific timeframe 1 4 . These methods could all exploit temporal elements of coordinated activity. We describe preliminary research regarding coordination network semantics, coordination network construction, relevant observations in three political Twitter datasets and the role of 'cheerleaders' to find social bots. See https github.com weberdc find hccs for code and data.


Distributed Deep Learning Using Volunteer Computing Like Paradigm. Use of Deep Learning (DL) in commercial applications such as image classification, sentiment analysis and speech recognition is increasing. When training DL models with large number of parameters and or large datasets, cost and speed of training can become prohibitive. Distributed DL training solutions that split a training job into subtasks and execute them over multiple nodes can decrease training time. However, the cost of current solutions, built predominantly for cluster computing systems, can still be an issue. In contrast to cluster computing systems, Volunteer Computing (VC) systems can lower the cost of computing, but applications running on VC systems have to handle fault tolerance, variable network latency and heterogeneity of compute nodes, and the current solutions are not designed to do so. We design a distributed solution that can run DL training on a VC system by using a data parallel approach. We implement a novel asynchronous SGD scheme called VC ASGD suited for VC systems. In contrast to traditional VC systems that lower cost by using untrustworthy volunteer devices, we lower cost by leveraging preemptible computing instances on commercial cloud platforms. By using preemptible instances that require applications to be fault tolerant, we lower cost by 70 90 and improve data security.


Malicious Selling Strategies in Livestream Shopping A Case Study of Alibaba's Taobao and ByteDance's Douyin. Livestream shopping is getting more and more popular as a new shopping form. Also, due to the COVID 19 pandemic, people have shifted to online shopping platforms. However, the broader user adoption comes with a cost many streamers' malicious selling incidents have recently been reported. In this study, we aim to explore streamers' malicious selling strategies and how viewers may perceive these strategies. First, we collected 40 livestream shopping sessions from two popular livestream platforms in China Taobao and Douyin (TikTok of Chinese version). We identified three categories of malicious selling strategies (i.e., Compulsive, Restrictive, and Designing) and found platform designs mostly enhanced these malicious selling strategies. Second, through an interview study with 13 end users, we provide a rich description of users' awareness of malicious selling strategies and challenges to counter malicious selling. We conclude the paper by discussing the policy and design implications to counter malicious selling.


Why Can You Lay Off Heads? Investigating How BERT Heads Transfer. The huge size of the widely used BERT family models has led to recent efforts about model distillation. The main goal of distillation is to create a task agnostic pre trained model that can be fine tuned on downstream tasks without fine tuning its full sized version. Despite the progress of distillation, to what degree and for what reason a task agnostic model can be created from distillation has not been well studied. Also, the mechanisms behind transfer learning of those BERT models are not well investigated either. Therefore, this work focuses on analyzing the acceptable deduction when distillation for guiding the future distillation procedure. Specifically, we first inspect the prunability of the Transformer heads in RoBERTa and ALBERT using their head importance estimation proposed by Michel et al. (2019), and then check the coherence of the important heads between the pre trained task and downstream tasks. Hence, the acceptable deduction of performance on the pre trained task when distilling a model can be derived from the results, and we further compare the behavior of the pruned model before and after fine tuning. Our studies provide guidance for future directions about BERT family model distillation.


140 GHz Urban Microcell Propagation Measurements for Spatial Consistency Modeling. Sub Terahertz frequencies (frequencies above 100 GHz) have the potential to satisfy the unprecedented demand on data rate on the order of hundreds of Gbps for sixth generation (6G) wireless communications and beyond. Accurate beam tracking and rapid beam selection are increasingly important since antenna arrays with more elements generate narrower beams to compensate for additional path loss within the first meter of propagation distance at sub THz frequencies. Realistic channel models for above 100 GHz are needed, and should include spatial consistency to model the spatial and temporal channel evolution along the user trajectory. This paper introduces recent outdoor urban microcell (UMi) propagation measurements at 142 GHz along a 39 m times 12 m rectangular route (102 m long), where each consecutive and adjacent receiver location is 3 m apart from each other. The measured power delay profiles and angular power spectrum at each receiver location are used to study spatial autocorrelation properties of various channel parameters such as shadow fading, delay spread, and angular spread along the track. Compared to the correlation distances reported in the 3GPP TR 38.901 for frequencies below 100 GHz, the measured correlation distance of shadow fading at 142 GHz (3.8 m) is much shorter than the 10 13 m as specified in 3GPP the measured correlation distances of delay spread and angular spread at 142 GHz (both 12 m) are comparable to the 7 10 m as specified in 3GPP. This result may guide the development of a statistical spatially consistent channel model for frequencies above 100 GHz in the UMi street canyon environment.


Retrieving linear and nonlinear optical dispersions of matter combined experiment numerical ellipsometry in Silicon, Gold and Indium Tin Oxide. The predominant methods currently used to determine nonlinear optical constants like the nonlinear refractive index n2 or chi3 rely mostly on experimental, open and closed z scan techniques and beam deflection methods. While these methods work well when the linear absorption is relatively small or negligible, the retrieval process is more complicated for a strongly scattering, dispersive or absorbing medium. The study of optics at the nanoscale in the ps or fs regimes demands the development of new theoretical tools experimental approaches, to extract and verify both linear and nonlinear optical dispersions exhibited by matter, especially when material constituents are fashioned into nanostructures of arbitrary shape. We present a practical, combined experimental and theoretical approach based on a hydrodynamic model that uses experimental results of harmonic generation conversion efficiencies to retrieve complex, nonlinear dispersion curves, not necessarily only for third order processes. We provide examples for materials that are of special interest to nanophotonics, silicon, gold, and indium tin oxide, which displays nonlocal effects and a zero crossing of the real part of the dielectric constant. The results for silicon and gold compare well with analytical predictions based on the nonlinear oscillator model. Based on our assessment of THG conversion efficiencies in silicon, we predict chi3(w) and chi3(3w) are of order 10 ( 17)(m V) 2, in the visible and IR ranges, with respective peaks of 10 ( 14) and 10 ( 16)(m V) 2 in the UV range. Similarly, gold's chi3(w) and chi3(3w) are of order 10 ( 17) and 10 ( 16)(m V) 2, and predict chi3(w) 10 ( 17)(m V) 2 and chi3(3w) 10 ( 18)(m V) 2 for ITO. These results suggest that judicious exploitation of the nonlinear dispersion of ordinary semiconductors can transform device physics in spectral regions that extend well into the UV range.


On Capital Allocation for a Risk Measure Derived from Ruin Theory. This paper addresses allocation methodologies for a risk measure inherited from ruin theory. Specifically, we consider a dynamic value at risk (VaR) measure defined as the smallest initial capital needed to ensure that the ultimate ruin probability is less than a given threshold. We introduce an intuitively appealing, novel allocation method, with a focus on its application to capital reserves which are determined through the dynamic value at risk (VaR) measure. Various desirable properties of the presented approach are derived including a limit result when considering a large time horizon and the comparison with the frequently used gradient allocation method. In passing, we introduce a second allocation method and discuss its relation to the other allocation approaches. A number of examples illustrate the applicability and performance of the allocation approaches.


DIVA Dataset Derivative of a Learning Task. We present a method to compute the derivative of a learning task with respect to a dataset. A learning task is a function from a training set to the validation error, which can be represented by a trained deep neural network (DNN). The "dataset derivative" is a linear operator, computed around the trained model, that informs how perturbations of the weight of each training sample affect the validation error, usually computed on a separate validation dataset. Our method, DIVA (Differentiable Validation) hinges on a closed form differentiable expression of the leave one out cross validation error around a pre trained DNN. Such expression constitutes the dataset derivative. DIVA could be used for dataset auto curation, for example removing samples with faulty annotations, augmenting a dataset with additional relevant samples, or rebalancing. More generally, DIVA can be used to optimize the dataset, along with the parameters of the model, as part of the training process without the need for a separate validation dataset, unlike bi level optimization methods customary in AutoML. To illustrate the flexibility of DIVA, we report experiments on sample auto curation tasks such as outlier rejection, dataset extension, and automatic aggregation of multi modal data.


Fast and Secure Key Generation with Channel Obfuscation in Slowly Varying Environments. The physical layer secret key generation has emerged as a promising solution for establishing cryptographic keys by leveraging reciprocal and time varying wireless channels. However, existing approaches suffer from low key generation rates and vulnerabilities under various attacks in slowly varying environments. We propose a new physical layer secret key generation approach with channel obfuscation, which improves the dynamic property of channel parameters based on random filtering and random antenna scheduling. Our approach makes one party obfuscate the channel to allow the legitimate party to obtain similar dynamic channel parameters yet prevents a third party from inferring the obfuscation information. Our approach allows more random bits to be extracted from the obfuscated channel parameters by a joint design of the K L transform and adaptive quantization. A testbed implementation shows that our approach, compared to the existing ones that we evaluate, performs the best in generating high entropy bits at a fast rate and a high security level in slowly varying environments. Specifically, our approach can achieve a significantly faster secret bit generation rate at about 67 bit pkt, and the key sequences can pass the randomness tests of the NIST test suite.


Data Driven VRP A Neural Network Model to Learn Hidden Preferences for VRP. The traditional Capacitated Vehicle Routing Problem (CVRP) minimizes the total distance of the routes under the capacity constraints of the vehicles. But more often, the objective involves multiple criteria including not only the total distance of the tour but also other factors such as travel costs, travel time, and fuel consumption.Moreover, in reality, there are numerous implicit preferences ingrained in the minds of the route planners and the drivers. Drivers, for instance, have familiarity with certain neighborhoods and knowledge of the state of roads, and often consider the best places for rest and lunch breaks. This knowledge is difficult to formulate and balance when operational routing decisions have to be made. This motivates us to learn the implicit preferences from past solutions and to incorporate these learned preferences in the optimization process. These preferences are in the form of arc probabilities, i.e., the more preferred a route is, the higher is the joint probability. The novelty of this work is the use of a neural network model to estimate the arc probabilities, which allows for additional features and automatic parameter estimation. This first requires identifying suitable features, neural architectures and loss functions, taking into account that there is typically few data available. We investigate the difference with a prior weighted Markov counting approach, and study the applicability of neural networks in this setting.


Spatial Temporal Transformer for Dynamic Scene Graph Generation. Dynamic scene graph generation aims at generating a scene graph of the given video. Compared to the task of scene graph generation from images, it is more challenging because of the dynamic relationships between objects and the temporal dependencies between frames allowing for a richer semantic interpretation. In this paper, we propose Spatial temporal Transformer (STTran), a neural network that consists of two core modules (1) a spatial encoder that takes an input frame to extract spatial context and reason about the visual relationships within a frame, and (2) a temporal decoder which takes the output of the spatial encoder as input in order to capture the temporal dependencies between frames and infer the dynamic relationships. Furthermore, STTran is flexible to take varying lengths of videos as input without clipping, which is especially important for long videos. Our method is validated on the benchmark dataset Action Genome (AG). The experimental results demonstrate the superior performance of our method in terms of dynamic scene graphs. Moreover, a set of ablative studies is conducted and the effect of each proposed module is justified. Code available at https github.com yrcong STTran.


An a posteriori strategy for adaptive schemes in time and space. A nonlinear adaptive procedure for optimising both the schemes in time and space is proposed in view of increasing the numerical efficiency and reducing the computational time. The method is based on a four parameter family of schemes we shall tune in function of the physical data (velocity, diffusion), the characteristic size in time and space, and the local regularity of the function leading to a nonlinear procedure. The textit a posteriori strategy we adopt consists in, given the solution at time t n , computing a candidate solution with the highest accurate schemes in time and space for all the nodes. Then, for the nodes that present some instabilities, both the schemes in time and space are modified and adapted in order to preserve the stability with a large time step. The updated solution is computed with node dependent schemes both in time and space. For the sake of simplicity, only convection diffusion problems are addressed as a prototype with a two parameters five points finite difference method for the spatial discretisation together with an explicit time two parameters four stages Runge Kutta method. We prove that we manage to obtain an optimal time step algorithm that produces accurate numerical approximations exempt of non physical oscillations.


On the Computation of the Algebraic Closure of Finitely Generated Groups of Matrices. We investigate the complexity of computing the Zariski closure of a finitely generated group of matrices. The Zariski closure was previously shown to be computable by Derksen, Jeandel, and Koiran, but the termination argument for their algorithm appears not to yield any complexity bound. In this paper we follow a different approach and obtain a bound on the degree of the polynomials that define the closure. Our bound shows that the closure can be computed in elementary time. We also obtain upper bounds on the length of chains of linear algebraic groups, where all the groups are generated over a fixed number field.


Variations of Ramanujan's Euler Products. We study the meromorphy of various Euler products of degree two attached to cusp forms including Ramanujan's Delta function.


A Note on Projection Based Recovery of Clusters in Markov Chains. Let T 0 be the transition matrix of a purely clustered Markov chain, i.e. a direct sum of k geq 2 irreducible stochastic matrices. Given a perturbation T(x) T 0 xE of T 0 such that T(x) is also stochastic, how small must x be in order for us to recover the indices of the direct summands of T 0 ? We give a simple algorithm based on the orthogonal projection matrix onto the left or right singular subspace corresponding to the k smallest singular values of I T(x) which allows for exact recovery all clusters when x O left( frac sigma n k E 2 sqrt n 1 right) and approximate recovery of a single cluster when x O left( frac sigma n k E 2 right) , where n 1 is the size of the largest cluster and sigma n k the (k 1) st smallest singular value of T 0 .


Consistent SPH simulations of the anisotropic dispersion of a contaminant plume. Solute transport through heterogeneous porous media is governed by fuid advection, molecular diffusion and anisotropic dispersion. The dispersion is assumed to obey Fick's law and the dispersion coefficient is defined as a second rank tensor. However, this problem has revealed to be a very difficult one because independently of the numerical methods employed, the solutions are seen to exhibit artificial oscillations and negative concentrations when the dispersivity becomes anisotropic. Here we report consistent SPH simulations of the anisotropic dispersion of a Gaussian contaminant plume in porous media using the open source code DualSPHysics. Consistency of the SPH method is restored by increasing the spatial resolution along with the number of neighbours within the compact support of the interpolating kernel. The solution shows that as the number of neighbours is increased with resolution, full convergence of the numerical solutions is guaranteed regardless of the dispersivity. However, despite the restored consistency, negative concentrations, albeit at a lower level, are still present. This suggests that a compromise between the number of neighbours and the size of the smoothing length must be guaranteed such that sufficient implicit numerical diffusion remains to damp out unphysical oscillations.


Beta families arising from a v 2 9 self map on S (3,v 1 8) . We show that v 2 9 is a permanent cycle in the 3 primary Adams Novikov spectral sequence computing pi (S (3,v 1 8)) , and use this to conclude that the families beta 9t 3 i for i 1,2 , beta 9t 6 i for i 1,2,3 , beta 9t 9 i for i 1, dots,8 , alpha 1 beta 9t 3 3 , and alpha 1 beta 9t 7 are permanent cycles in the 3 primary Adams Novikov spectral sequence for the sphere for all t geq 0 . We use a computer program by Wang to determine the additive and partial multiplicative structure of the Adams Novikov E 2 page for the sphere in relevant degrees. The i 1 cases recover previously known results of Behrens Pemmaraju and the second author. The results about beta 9t 3 3 , beta 9t 6 3 and beta 9t 8 9 were previously claimed by the second author the computer calculations allow us to give a more direct proof. As an application, we determine the image of the Hurewicz map pi S to pi tmf at p 3 .


Optimal Sequential Stochastic Deployment of Multiple Passenger Robots. We present a new algorithm for deploying passenger robots in marsupial robot systems. A marsupial robot system consists of a carrier robot (e.g., a ground vehicle), which is highly capable and has a long mission duration, and at least one passenger robot (e.g., a short duration aerial vehicle) transported by the carrier. We optimize the performance of passenger robot deployment by proposing an algorithm that reasons over uncertainty by exploiting information about the prior probability distribution of features of interest in the environment. Our algorithm is formulated as a solution to a sequential stochastic assignment problem (SSAP). The key feature of the algorithm is a recurrence relationship that defines a set of observation thresholds that are used to decide when to deploy passenger robots. Our algorithm computes the optimal policy in O(NR) time, where N is the number of deployment decision points and R is the number of passenger robots to be deployed. We conducted drone deployment exploration experiments on real world data from the DARPA Subterranean challenge to test the SSAP algorithm. Our results show that our deployment algorithm outperforms other competing algorithms, such as the classic secretary approach and baseline partitioning methods, and is comparable to an offline oracle algorithm.


Byzantine robust Federated Learning through Spatial temporal Analysis of Local Model Updates. Federated Learning (FL) enables multiple distributed clients (e.g., mobile devices) to collaboratively train a centralized model while keeping the training data locally on the client. Compared to traditional centralized machine learning, FL offers many favorable features such as offloading operations which would usually be performed by a central server and reducing risks of serious privacy leakage. However, Byzantine clients that send incorrect or disruptive updates due to system failures or adversarial attacks may disturb the joint learning process, consequently degrading the performance of the resulting model. In this paper, we propose to mitigate these failures and attacks from a spatial temporal perspective. Specifically, we use a clustering based method to detect and exclude incorrect updates by leveraging their geometric properties in the parameter space. Moreover, to further handle malicious clients with time varying behaviors, we propose to adaptively adjust the learning rate according to momentum based update speculation. Extensive experiments on 4 public datasets demonstrate that our algorithm achieves enhanced robustness comparing to existing methods under both cross silo and cross device FL settings with faulty malicious clients.


Analysis of a semi augmented mixed finite element method for double diffusive natural convection in porous media. In this paper we study a stationary double diffusive natural convection problem in porous media given by a Navier Stokes Darcy type system, for describing the velocity and the pressure, coupled to a vector advection diffusion equation describing the heat and substance concentration, of a viscous fluid in a porous media with physical boundary conditions. The model problem is rewritten in terms of a first order system, without the pressure, based on the introduction of the strain tensor and a nonlinear pseudo stress tensor in the fluid equations. After a variational approach, the resulting weak model is then augmented using appropriate redundant penalization terms for the fluid equations along with a standard primal formulation for the heat and substance concentration. Then, it is rewritten as an equivalent fixed point problem. Well posedness and uniqueness results for both the continuous and the discrete schemes are stated, as well as the respective convergence result under certain regularity assumptions combined with the Lax Milgram theorem, and the Banach and Brouwer fixed point theorems. In particular, Raviart Thomas elements of order k are used for approximating the pseudo stress tensor, piecewise polynomials of degree leq k and leq k 1 are utilized for approximating the strain tensor and the velocity, respectively, and the heat and substance concentration are approximated by means of Lagrange finite elements of order leq k 1 . Optimal a priori error estimates are derived and confirmed through some numerical examples that illustrate the performance of the proposed semi augmented mixed primal scheme.


Diffusion free scaling in rotating spherical Rayleigh B 'enard convection. Direct numerical simulations are employed to reveal three distinctly different flow regions in rotating spherical Rayleigh B 'enard convection. In the low latitude region mathrm I vertical (parallel to the axis of rotation) convective columns are generated between the hot inner and the cold outer sphere. The mid latitude region mathrm II is dominated by vertically aligned convective columns formed between the Northern and Southern hemispheres of the outer sphere. The diffusion free scaling, which indicates bulk dominated convection, originates from this mid latitude region. In the equator region mathrm III the vortices are affected by the outer spherical boundary and are much shorter than in region mathrm II . Thermally driven turbulence with background rotation in spherical Rayleigh B 'enard convection is found to be characterized by three distinctly different flow regions. The diffusion free scaling, which indicates the heat transfer is bulk dominated, originates from the mid latitude region in which vertically aligned vortices are stretched between the Northern and Southern hemispheres of the outer sphere. These results show that the flow physics in rotating convection are qualitatively different in planar and spherical geometries. This finding underlines that it is crucial to study convection in spherical geometries to better understand geophysical and astrophysical flow phenomena.


Integrability properties of quasi regular representations of NA groups. Let G N rtimes A , where N is a graded Lie group and A mathbb R acts on N via homogeneous dilations. The quasi regular representation pi mathrm ind A G (1) of G can be realised to act on L 2 (N) . It is shown that for a class of analysing vectors the associated wavelet transform defines an isometry from L 2 (N) into L 2 (G) and that the integral kernel of the corresponding orthogonal projector has polynomial off diagonal decay. The obtained reproducing formula is instrumental for proving decomposition theorems for function spaces on nilpotent Lie groups.


Global Index on Financial Losses due to Crime in the United States. Crime can have a volatile impact on investments. Despite the potential importance of crime rates in investments, there are no indices dedicated to evaluating the financial impact of crime in the United States. As such, this paper presents an index based insurance portfolio for crime in the United States by utilizing the financial losses reported by the Federal Bureau of Investigation for property crimes and cybercrimes. Our research intends to help investors envision risk exposure in our portfolio, gauge investment risk based on their desired risk level, and hedge strategies for potential losses due to economic crashes. Underlying the index, we hedge the investments by issuing marketable European call and put options and providing risk budgets (diversifying risk to each type of crime). We find that real estate, ransomware, and government impersonation are the main risk contributors. We then evaluate the performance of our index to determine its resilience to economic crisis. The unemployment rate potentially demonstrates a high systemic risk on the portfolio compared to the economic factors used in this study. In conclusion, we provide a basis for the securitization of insurance risk from certain crimes that could forewarn investors to transfer their risk to capital market investors.


Predicting Driver Self Reported Stress by Analyzing the Road Scene. Several studies have shown the relevance of biosignals in driver stress recognition. In this work, we examine something important that has been less frequently explored We develop methods to test if the visual driving scene can be used to estimate a drivers' subjective stress levels. For this purpose, we use the AffectiveROAD video recordings and their corresponding stress labels, a continuous human driver provided stress metric. We use the common class discretization for stress, dividing its continuous values into three classes low, medium, and high. We design and evaluate three computer vision modeling approaches to classify the driver's stress levels (1) object presence features, where features are computed using automatic scene segmentation (2) end to end image classification and (3) end to end video classification. All three approaches show promising results, suggesting that it is possible to approximate the drivers' subjective stress from the information found in the visual scene. We observe that the video classification, which processes the temporal information integrated with the visual information, obtains the highest accuracy of 0.72 , compared to a random baseline accuracy of 0.33 when tested on a set of nine drivers.


Enumerating chambers of hyperplane arrangements with symmetry. We introduce a new algorithm for enumerating chambers of hyperplane arrangements which exploits their underlying symmetry groups. Our algorithm counts the chambers of an arrangement as a byproduct of computing its characteristic polynomial. We showcase our julia implementation, based on OSCAR, on examples coming from hyperplane arrangements with applications to physics and computer science.


Modelling Art Interpretation and Meaning. A Data Model for Describing Iconology and Iconography. Iconology is a branch of art history that investigates the meaning of artworks in relation to their social and cultural background. Nowadays, several interdisciplinary research fields leverage theoretical frameworks close to iconology to pursue quantitative Art History with data science methods and Semantic Web technologies. However, while Iconographic studies have been recently addressed in ontologies, a complete description of aspects relevant to iconological studies is still missing. In this article, we present a preliminary study on eleven case studies selected from the literature and we envision new terms for extending existing ontologies. We validate new terms according to a common evaluation method and we discuss our results in the light of the opportunities that such an extended ontology would arise in the community of Digital Art History.


On an extremal problem for locally sparse multigraphs. A multigraph G is an (s,q) graph if every s set of vertices in G supports at most q edges of G , counting multiplicities. Mubayi and Terry posed the problem of determining the maximum of the product of the edge multiplicities in an (s,q) graph on n vertices. We give an asymptotic solution to this problem for the family (s,q) (2r, a binom 2r 2 mathrm ex (2r, K r 1 ) 1 ) with r, a in mathbb Z geq 2 . This greatly generalises previous results on the problem due to Mubayi and Terry and to Day, Treglown and the author, who between them had resolved the special case r 2 , and asymptotically confirms an infinite family of cases in a conjecture of Day, Treglown and the author.


Displacement pseudostress formulation for the linear elasticity spectral problem. In this paper we analyze a mixed displacement pseudostress formulation for the elasticity eigenvalue problem. We propose a finite element method to approximate the pseudostress tensor with Raviart Thomas elements and the displacement with piecewise polynomials. With the aid of the classic theory for compact operators, we prove that our method is convergent and does not introduce spurious modes. Also, we obtain error estimates for the proposed method. Finally, we report some numerical tests supporting the theoretical results.


Yetter Drinfeld modules for group cograded Hopf quasigroups. Let H be a crossed group cograded Hopf quasigroup. We first introduce the notion of p Yetter Drinfeld quasimodule over H . If the antipode of H is bijective, we show that the category mathscr Y mathscr D mathscr Q(H) of Yetter Drinfeld quasimodules over H is a crossed category, and the subcategory mathscr Y mathscr D(H) of Yetter Drinfeld modules is a braided crossed category.


Relative velocity based reward functions for crowd navigation of robots. The four wheeled Mecanum robot is widely used in various industries due to its maneuverability and strong load capacity, which is suitable for performing precise transportation tasks in a narrow environment, but while the Mecanum wheel robot has mobility, it also consumes more energy than ordinary robots. The power consumed by the Mecanum wheel mobile robot varies enormously depending on their operating regimes and environments. Therefore, only knowing the working environment of the robot and the accurate power consumption model can we accurately predict the power consumption of the robot. In order to increase the appli cable scenarios of energy consumption modeling for Mecanum wheel robots and improve the accuracy of energy consumption modeling, this paper focuses on various factors that affect the energy consumption of the Mecanum wheel robot, such as motor temperature, terrain, the center of gravity position, etc. The model is derived from the kinematic and kinetic model combined with electrical engineering and energy flow principles. The model has been simulated in MATLAB and experimentally validated with the four wheeled Mecanum robot platform in our lab. Experimental results show that the model is 90 accurate. The results of energy consumption modeling can help robots to save energy by helping them to perform rational path planning and task planning.


Potent preservers of incidence algebras. Let X be a finite connected poset, F a field and I(X,F) the incidence algebra of X over F . We describe the bijective linear idempotent preservers varphi I(X,F) to I(X,F) . Namely, we prove that, whenever mathrm char (F) ne 2 , varphi is either an automorphism or an anti automorphism of I(X,F) . If mathrm char (F) 2 and F 2 , then varphi is a (in general, non proper) Lie automorphism of I(X,F) . Finally, if F mathbb Z 2 , then varphi is the composition of a bijective shift map and a Lie automorphism of I(X,F) . Under certain restrictions on the characteristic of F we also obtain descriptions of the bijective linear maps which preserve tripotents and, more generally, k potents of I(X,F) for k ge 3 .


Stop Building Castles on a Swamp! The Crisis of Reproducing Automatic Search in Evidence based Software Engineering. The evidence based approach has increasingly been employed to synthesize empirical findings from the primary research in software engineering. Nevertheless, the reproducibility of evidence based software engineering (EBSE) studies seems to be underemphasized. In our investigation into the automatic search of 311 sample studies, more than 50 of the search strings are not reusable about 87.5 of the search activities (e.g., search field settings) are unrepeatable and more than 95 of the whole automatic search implementations are unreproducible. Considering that searching is a cornerstone of an EBSE study, we are afraid that the reproducibility of the current secondary research could be worse than we can imagine. By analyzing and reporting the root causes of the aforementioned observations, we urge collaboration and cooperation among all the stakeholders in our community to improve the research reproducibility in EBSE.


Randomized block Krylov methods for approximating extreme eigenvalues. Randomized block Krylov subspace methods form a powerful class of algorithms for computing the extreme eigenvalues of a symmetric matrix or the extreme singular values of a general matrix. The purpose of this paper is to develop new theoretical bounds on the performance of randomized block Krylov subspace methods for these problems. For matrices with polynomial spectral decay, the randomized block Krylov method can obtain an accurate spectral norm estimate using only a constant number of steps (that depends on the decay rate and the accuracy). Furthermore, the analysis reveals that the behavior of the algorithm depends in a delicate way on the block size. Numerical evidence confirms these predictions.


A homotopic approach to policy gradients for linear quadratic regulators with nonlinear controls. We study the convergence of deterministic policy gradient algorithms in continuous state and action space for the prototypical Linear Quadratic Regulator (LQR) problem when the search space is not limited to the family of linear policies. We first provide a counterexample showing that extending the policy class to piecewise linear functions results in local minima of the policy gradient algorithm. To solve this problem, we develop a new approach that involves sequentially increasing a discount factor between iterations of the original policy gradient algorithm. We finally prove that this homotopic variant of policy gradient methods converges to the global optimum of the undiscounted Linear Quadratic Regulator problem for a large class of Lipschitz, non linear policies.


Okutsu frames of irreducible polynomials over henselian fields. For a henselian valued field (K,v) we establish a complete parallelism between the arithmetic properties of irreducible polynomials F in K x , encoded by their Okutsu frames, and the valuation theoretic properties of their induced valuations v F on K x , encoded by their MacLane Vaqui 'e chains. This parallelism was only known for defectless irreducible polynomials.


Pythagorean Triples, Complex Numbers, Abelian Groups and Prime Numbers. It is well known that pythagorean triples can be represented by points of the unit circle with rational coordinates. These points form an abelian group, and we describe its structure. This structural description yields, almost immediately, an enumeration of the normalized pythagorean triples with a given hypotenuse, and also to an effective method for producing all such triples. This effective method seems to be new. This paper is intended for the general mathematical audience, including undergraduate mathematics students, and therefore it contains plenty of background material, some history and several examples and exercises.


Querying Linked Data how to ensure user's quality requirements. In the distributed and dynamic framework of the Web, data quality is a big challenge. The Linked Open Data (LOD) provides an enormous amount of data, the quality of which is difficult to control. Quality is intrinsically a matter of usage, so consumers need ways to specify quality rules that make sense for their use, in order to get only data conforming to these rules. We propose a user side query framework equipped with a checker of constraints and confidence levels on data resulting from LOD providers ' query evaluations. We detail its theoretical foundations and we provide experimental results showing that the check additional cost is reasonable and that integrating the constraints in the queries further improves it significantly.


Scientific Machine Learning Benchmarks. The breakthrough in Deep Learning neural networks has transformed the use of AI and machine learning technologies for the analysis of very large experimental datasets. These datasets are typically generated by large scale experimental facilities at national laboratories. In the context of science, scientific machine learning focuses on training machines to identify patterns, trends, and anomalies to extract meaningful scientific insights from such datasets. With a new generation of experimental facilities, the rate of data generation and the scale of data volumes will increasingly require the use of more automated data analysis. At present, identifying the most appropriate machine learning algorithm for the analysis of any given scientific dataset is still a challenge for scientists. This is due to many different machine learning frameworks, computer architectures, and machine learning models. Historically, for modelling and simulation on HPC systems such problems have been addressed through benchmarking computer applications, algorithms, and architectures. Extending such a benchmarking approach and identifying metrics for the application of machine learning methods to scientific datasets is a new challenge for both scientists and computer scientists. In this paper, we describe our approach to the development of scientific machine learning benchmarks and review other approaches to benchmarking scientific machine learning.


Equidistribution of torsion points in abelian varieties. We prove that the torsion points of an abelian variety are equidistributed over the corresponding berkovich space with respect to the canonical measure.


Unsupervised Representation Learning for 3D Point Cloud Data. Though a number of point cloud learning methods have been proposed to handle unordered points, most of them are supervised and require labels for training. By contrast, unsupervised learning of point cloud data has received much less attention to date. In this paper, we propose a simple yet effective approach for unsupervised point cloud learning. In particular, we identify a very useful transformation which generates a good contrastive version of an original point cloud. They make up a pair. After going through a shared encoder and a shared head network, the consistency between the output representations are maximized with introducing two variants of contrastive losses to respectively facilitate downstream classification and segmentation. To demonstrate the efficacy of our method, we conduct experiments on three downstream tasks which are 3D object classification (on ModelNet40 and ModelNet10), shape part segmentation (on ShapeNet Part dataset) as well as scene segmentation (on S3DIS). Comprehensive results show that our unsupervised contrastive representation learning enables impressive outcomes in object classification and semantic segmentation. It generally outperforms current unsupervised methods, and even achieves comparable performance to supervised methods. Our source codes will be made publicly available.


Analytical Process Scheduling Optimization for Heterogeneous Multi core Systems. In this paper, we propose the first optimum process scheduling algorithm for an increasingly prevalent type of heterogeneous multicore (HEMC) system that combines high performance big cores and energy efficient small cores with the same instruction set architecture (ISA). Existing algorithms are all heuristics based, and the well known IPC driven approach essentially tries to schedule high scaling factor processes on big cores. Our analysis shows that, for optimum solutions, it is also critical to consider placing long running processes on big cores. Tests of SPEC 2006 cases on various big small core combinations show that our proposed optimum approach is up to 34 faster than the IPC driven heuristic approach in terms of total workload completion time. The complexity of our algorithm is O(NlogN) where N is the number of processes. Therefore, the proposed optimum algorithm is practical for use.


Uncertainty aware Safe Exploratory Planning using Gaussian Process and Neural Control Contraction Metric. In this paper, we consider the problem of using a robot to explore an environment with an unknown, state dependent disturbance function while avoiding some forbidden areas. The goal of the robot is to safely collect observations of the disturbance and construct an accurate estimate of the underlying disturbance function. We use Gaussian Process (GP) to get an estimate of the disturbance from data with a high confidence bound on the regression error. Furthermore, we use neural Contraction Metrics to derive a tracking controller and the corresponding high confidence uncertainty tube around the nominal trajectory planned for the robot, based on the estimate of the disturbance. From the robustness of the Contraction Metric, error bound can be pre computed and used by the motion planner such that the actual trajectory is guaranteed to be safe. As the robot collects more and more observations along its trajectory, the estimate of the disturbance becomes more and more accurate, which in turn improves the performance of the tracking controller and enlarges the free space that the robot can safely explore. We evaluate the proposed method using a carefully designed environment with a ground vehicle. Results show that with the proposed method the robot can thoroughly explore the environment safely and quickly.


Offensive Language Identification in Low resourced Code mixed Dravidian languages using Pseudo labeling. Social media has effectively become the prime hub of communication and digital marketing. As these platforms enable the free manifestation of thoughts and facts in text, images and video, there is an extensive need to screen them to protect individuals and groups from offensive content targeted at them. Our work intends to classify codemixed social media comments posts in the Dravidian languages of Tamil, Kannada, and Malayalam. We intend to improve offensive language identification by generating pseudo labels on the dataset. A custom dataset is constructed by transliterating all the code mixed texts into the respective Dravidian language, either Kannada, Malayalam, or Tamil and then generating pseudo labels for the transliterated dataset. The two datasets are combined using the generated pseudo labels to create a custom dataset called CMTRA. As Dravidian languages are under resourced, our approach increases the amount of training data for the language models. We fine tune several recent pretrained language models on the newly constructed dataset. We extract the pretrained language embeddings and pass them onto recurrent neural networks. We observe that fine tuning ULMFiT on the custom dataset yields the best results on the code mixed test sets of all three languages. Our approach yields the best results among the benchmarked models on Tamil English, achieving a weighted F1 Score of 0.7934 while scoring competitive weighted F1 Scores of 0.9624 and 0.7306 on the code mixed test sets of Malayalam English and Kannada English, respectively.


Choose a Transformer Fourier or Galerkin. In this paper, we apply the self attention from the state of the art Transformer in Attention Is All You Need for the first time to a data driven operator learning problem related to partial differential equations. An effort is put together to explain the heuristics of, and to improve the efficacy of the attention mechanism. By employing the operator approximation theory in Hilbert spaces, it is demonstrated for the first time that the softmax normalization in the scaled dot product attention is sufficient but not necessary. Without softmax, the approximation capacity of a linearized Transformer variant can be proved to be comparable to a Petrov Galerkin projection layer wise, and the estimate is independent with respect to the sequence length. A new layer normalization scheme mimicking the Petrov Galerkin projection is proposed to allow a scaling to propagate through attention layers, which helps the model achieve remarkable accuracy in operator learning tasks with unnormalized data. Finally, we present three operator learning experiments, including the viscid Burgers' equation, an interface Darcy flow, and an inverse interface coefficient identification problem. The newly proposed simple attention based operator learner, Galerkin Transformer, shows significant improvements in both training cost and evaluation accuracy over its softmax normalized counterparts.


Further study of Squeezing function Corresponding to Polydisk. In the present article, we further investigate the properties of squeezing function corresponding to polydisk. We work out some explicit expressions of squeezing function corresponding to polydisk.


Combining the hybrid mimetic mixed method with the Scharfetter Gummel scheme for magnetised transport in plasmas. In this paper, we propose a numerical scheme for fluid models of magnetised plasmas. One important feature of the numerical scheme is that it should be able to handle the anisotropy induced by the magnetic field. In order to do so, we propose the use of the hybrid mimetic mixed (HMM) scheme for diffusion. This is combined with a hybridised variant of the Scharfetter Gummel (SG) scheme for advection. The proposed hybrid scheme can be implemented very efficiently via static condensation. Numerical tests are then performed to show the applicability of the combined HMM SG scheme, even for highly anisotropic magnetic fields.


Counting rational points close to p adic integers and applications in Diophantine approximation. We find upper and lower bounds on the number of rational points that are psi approximations of some n dimensional p adic integer. Lattice point counting techniques are used to find the upper bound result, and a Pigeon hole principle style argument is used to find the lower bound result. We use these results to find the Hausdorff dimension for the set of p adic weighted simultaneously approximable points intersected with p adic coordinate hyperplanes. For the lower bound result we show that the set of rational points that tau approximate a p adic integer form a set of resonant points that can be used to construct a local ubiquitous system of rectangles.


Lotka Volterra competition diffusion system the critical competition case. We consider the reaction diffusion competition system in the so called it critical competition case . The associated ODE system then admits infinitely many equilibria, which makes the analysis intricate. We first prove the non existence of it ultimately monotone traveling waves by applying the phase plane analysis. Next, we study the large time behavior of the solution of the Cauchy problem with a compactly supported initial datum. We not only reveal that the "faster" species excludes the "slower" one (with a known it spreading speed ), but also provide a sharp description of the profile of the solution, thus shedding light on a new it bump phenomenon .


VizLinter A Linter and Fixer Framework for Data Visualization. Despite the rising popularity of automated visualization tools, existing systems tend to provide direct results which do not always fit the input data or meet visualization requirements. Therefore, additional specification adjustments are still required in real world use cases. However, manual adjustments are difficult since most users do not necessarily possess adequate skills or visualization knowledge. Even experienced users might create imperfect visualizations that involve chart construction errors. We present a framework, VizLinter, to help users detect flaws and rectify already built but defective visualizations. The framework consists of two components, (1) a visualization linter, which applies well recognized principles to inspect the legitimacy of rendered visualizations, and (2) a visualization fixer, which automatically corrects the detected violations according to the linter. We implement the framework into an online editor prototype based on Vega Lite specifications. To further evaluate the system, we conduct an in lab user study. The results prove its effectiveness and efficiency in identifying and fixing errors for data visualizations.


The full renormalization horseshoe for multimodal maps. In this paper, we consider the renormalization operator mathcal R for multimodal maps. We prove the renormalization operator mathcal R is a self homeomorphism on any totally mathcal R invariant set. As a corollary, we prove the existence of the full renormalization horseshoe for multimodal maps.


On regular systems of finite classical polar spaces. Let cal P be a finite classical polar space of rank d . An m regular system with respect to (k 1) dimensional projective spaces of cal P , 1 le k le d 1 , is a set cal R of generators of cal P with the property that every (k 1) dimensional projective space of cal P lies on exactly m generators of cal R . Regular systems of polar spaces are investigated. Some non existence results about certain 1 regular systems of polar spaces with low rank are proved and a procedure to obtain m' regular systems from a given m regular system is described. Finally, three different construction methods of regular systems w.r.t. points of various polar spaces are discussed.


Learning disentangled representations via product manifold projection. We propose a novel approach to disentangle the generative factors of variation underlying a given set of observations. Our method builds upon the idea that the (unknown) low dimensional manifold underlying the data space can be explicitly modeled as a product of submanifolds. This definition of disentanglement gives rise to a novel weakly supervised algorithm for recovering the unknown explanatory factors behind the data. At training time, our algorithm only requires pairs of non i.i.d. data samples whose elements share at least one, possibly multidimensional, generative factor of variation. We require no knowledge on the nature of these transformations, and do not make any limiting assumption on the properties of each subspace. Our approach is easy to implement, and can be successfully applied to different kinds of data (from images to 3D surfaces) undergoing arbitrary transformations. In addition to standard synthetic benchmarks, we showcase our method in challenging real world applications, where we compare favorably with the state of the art.


SOON Scenario Oriented Object Navigation with Graph based Exploration. The ability to navigate like a human towards a language guided target from anywhere in a 3D embodied environment is one of the 'holy grail' goals of intelligent robots. Most visual navigation benchmarks, however, focus on navigating toward a target from a fixed starting point, guided by an elaborate set of instructions that depicts step by step. This approach deviates from real world problems in which human only describes what the object and its surrounding look like and asks the robot to start navigation from anywhere. Accordingly, in this paper, we introduce a Scenario Oriented Object Navigation (SOON) task. In this task, an agent is required to navigate from an arbitrary position in a 3D embodied environment to localize a target following a scene description. To give a promising direction to solve this task, we propose a novel graph based exploration (GBE) method, which models the navigation state as a graph and introduces a novel graph based exploration approach to learn knowledge from the graph and stabilize training by learning sub optimal trajectories. We also propose a new large scale benchmark named From Anywhere to Object (FAO) dataset. To avoid target ambiguity, the descriptions in FAO provide rich semantic scene information includes object attribute, object relationship, region description, and nearby region description. Our experiments reveal that the proposed GBE outperforms various state of the arts on both FAO and R2R datasets. And the ablation studies on FAO validates the quality of the dataset.


Distributed stochastic inertial accelerated methods with delayed derivatives for nonconvex problems. Stochastic gradient methods (SGMs) are predominant approaches for solving stochastic optimization. On smooth nonconvex problems, a few acceleration techniques have been applied to improve the convergence rate of SGMs. However, little exploration has been made on applying a certain acceleration technique to a stochastic subgradient method (SsGM) for nonsmooth nonconvex problems. In addition, few efforts have been made to analyze an (accelerated) SsGM with delayed derivatives. The information delay naturally happens in a distributed system, where computing workers do not coordinate with each other. In this paper, we propose an inertial proximal SsGM for solving nonsmooth nonconvex stochastic optimization problems. The proposed method can have guaranteed convergence even with delayed derivative information in a distributed environment. Convergence rate results are established to three classes of nonconvex problems weakly convex nonsmooth problems with a convex regularizer, composite nonconvex problems with a nonsmooth convex regularizer, and smooth nonconvex problems. For each problem class, the convergence rate is O(1 K frac 1 2 ) in the expected value of the gradient norm square, for K iterations. In a distributed environment, the convergence rate of the proposed method will be slowed down by the information delay. Nevertheless, the slow down effect will decay with the number of iterations for the latter two problem classes. We test the proposed method on three applications. The numerical results clearly demonstrate the advantages of using the inertial based acceleration. Furthermore, we observe higher parallelization speed up in asynchronous updates over the synchronous counterpart, though the former uses delayed derivatives. Our source code is released at https github.com RPI OPT Inertial SsGM


Triharmonic Curves in f Kenmotsu Manifolds. The aim of this paper is to study triharmonic curves in three dimensional f Kenmotsu manifolds. We investigate necessary and sufficient conditions for Frenet curves, and specifically for slant and Legendre curves to be triharmonic. Then we prove that triharmonic Frenet curves with constant curvature are Frenet helices in three dimensional f Kenmotsu manifolds. Next, we give a nonexistence theorem that there is no triharmonic Legendre curve in three dimensional f Kenmotsu manifolds.


On the equivalence between weak BMO and the space of derivatives of the Zygmund class. In this paper, we will discuss the space of functions of weak bounded mean oscillation. In particular, we will show that this space is the dual space of the special atom space, whose dual space was already known to be the space of derivative of functions (in the sense of distribution) belonging to the Zygmund class of functions. We show in particular that this proves that the Hardy space H 1 strictly contains the space special atom space.


Atomistic insights into hydrogen bonds of water around detonation nanodiamonds effects of surface chemistry and dissolved ions. Prior studies on the effects of surface chemistry of detonation nanodiamonds (DNDs) on Hydrogen Bonds (HBs) of water have not unambiguously distinguished among three broad categories of HBs water site HBs (HB (s) ), water water HBs in hydration layers of DNDs (HB (h) ), and water water HBs in the bulk region (HB (b) ). Furthermore, the effects of dissolved ions in the solution on the aforementioned HBs have not been studied yet. In this study, we addressed these issues via Molecular Dynamics simulations. Each system under study contained a single DND that was functionalized either with purely H or with a mixture of H and one of NH 2 , COOH, or OH groups. Then, the product was dissolved in 0.1M aqueous solution of any of KCl, NaCl, CaCl 2 , or MgCl 2 salts. All DNDs except for DND OH, which was charge neutral, existed as either a charged or a neutral particle in the solution. Whereas DND H with 56 net charges induced the strongest HB (h) , the addition of 28 charges to DND H led to a drop in HB (h) strength and became identical to what we observed for HB (h) strength around DND OH. Considering only neutral DNDs, HB (h) strength increased in the sequence as DND H approx DND COOH DND OH DND NH 2 . For the case of the negatively charged DND COOH, while HB (h) in either NaCl or KCl solution was slightly weaker than HB (b) , HB (h) strength corresponding to CaCl 2 or MgCl 2 solutions were enhanced compared to HB (b) . Interestingly, Mg 2 and to a lesser extent Ca 2 around the charged DND COOH were associated with the strongest HB (s) in this study.


Relative generalized Hamming weights of evaluation codes. The aim of this work is to algebraically describe the relative generalized Hamming weights of evaluation codes. We give a lower bound for these weights in terms of a footprint bound. We prove that this bound can be sharp. We compute the next to minimal weight of toric codes over hypersimplices of degree 1.


LogDP Combining Dependency and Proximity for Log based Anomaly Detection. Log analysis is an important technique that engineers use for troubleshooting faults of large scale service oriented systems. In this study, we propose a novel semi supervised log based anomaly detection approach, LogDP, which utilizes the dependency relationships among log events and proximity among log sequences to detect the anomalies in massive unlabeled log data. LogDP divides log events into dependent and independent events, then learns normal patterns of dependent events using dependency and independent events using proximity. Events violating any normal pattern are identified as anomalies. By combining dependency and proximity, LogDP is able to achieve high detection accuracy. Extensive experiments have been conducted on real world datasets, and the results show that LogDP outperforms six state of the art methods.


Constraints on the Dimensionality of Space. Complex structures can only form in a universe that allows for bound states. While this is clearly observed in three dimensions, added degrees of freedom in a higher dimensional space preclude the immediate assumption that binding potentials can in fact exist. In this paper, we derive a constraint on the dimensionality of a universe in the presence of an arbitrary set of forces. We then apply this constraint to systems with several example potentials. In doing so, we find that bound states in higher than 3 dimensions are in fact possible under specific circumstances which we characterize.


Translates of S arithmetic orbits and applications. We prove that certain sequences of periodic orbits of the diagonal group in the space of lattices equidistribute. As an application we obtain new information regarding the sequence of best approximations to certain vectors with algebraic coordinates. In order to prove these results we generalize the seminal work of Eskin Mozes and Shah about the equidistribution of translates of periodic measures from the real case to the S arithmetic case.


Data driven turbulence modeling in separated flows considering physical mechanism analysis. Accurate simulation of turbulent flow with separation is an important but challenging problem. In this paper, a data driven Reynolds averaged turbulence modeling approach, field inversion and machine learning is implemented to modify the Spalart Allmaras model separately on three cases, namely, the S809 airfoil, a periodic hill and the GLC305 airfoil with ice shape 944. Field inversion based on a discrete adjoint method is used to quantify the model form uncertainty with limited experimental data. An artificial neural network is trained to predict the model corrections with local flow features to extract generalized modeling knowledge. Physical knowledge of the nonequilibrium turbulence in the separating shear layer is considered when setting the prior model uncertainty. The results show that the model corrections from the field inversion demonstrate strong consistency with the underlying physical mechanism of nonequilibrium turbulence. The quantity of interest from the observation data can be reproduced with relatively high accuracy by the augmented model. In addition, the validation in similar flow conditions shows a certain extent of generalization ability.


Resilient Active Information Acquisition with Teams of Robots. Emerging applications of collaborative autonomy, such as Multi Target Tracking, Unknown Map Exploration, and Persistent Surveillance, require robots plan paths to navigate an environment while maximizing the information collected via on board sensors. In this paper, we consider such information acquisition tasks but in adversarial environments, where attacks may temporarily disable the robots' sensors. We propose the first receding horizon algorithm, aiming for robust and adaptive multi robot planning against any number of attacks, which we call Resilient Active Information acquisitioN (RAIN). RAIN calls, in an online fashion, a Robust Trajectory Planning (RTP) subroutine which plans attack robust control inputs over a look ahead planning horizon. We quantify RTP's performance by bounding its suboptimality. We base our theoretical analysis on notions of curvature introduced in combinatorial optimization. We evaluate RAIN in three information acquisition scenarios Multi Target Tracking, Occupancy Grid Mapping, and Persistent Surveillance. The scenarios are simulated in C and a Unity based simulator. In all simulations, RAIN runs in real time, and exhibits superior performance against a state of the art baseline information acquisition algorithm, even in the presence of a high number of attacks. We also demonstrate RAIN's robustness and effectiveness against varying models of attacks (worst case and random), as well as, varying replanning rates.


Higman Thompson groups from self similar groupoid actions. Given a self similar groupoid action (G,E) on a finite directed graph, we prove some properties of the corresponding ample groupoid of germs mathcal G(G,E) . We study the analogue of the Higman Thompson group associated to (G,E) using G tables and relate it to the topological full group of mathcal G(G,E) , which is isomorphic to a subgroup of unitaries in the algebra C (G,E) . After recalling some concepts in groupoid homology, we discuss the Matui's AH conjecture for mathcal G(G,E) in some particular cases.


How Metaphors Impact Political Discourse A Large Scale Topic Agnostic Study Using Neural Metaphor Detection. Metaphors are widely used in political rhetoric as an effective framing device. While the efficacy of specific metaphors such as the war metaphor in political discourse has been documented before, those studies often rely on small number of hand coded instances of metaphor use. Larger scale topic agnostic studies are required to establish the general persuasiveness of metaphors as a device, and to shed light on the broader patterns that guide their persuasiveness. In this paper, we present a large scale data driven study of metaphors used in political discourse. We conduct this study on a publicly available dataset of over 85K posts made by 412 US politicians in their Facebook public pages, up until Feb 2017. Our contributions are threefold we show evidence that metaphor use correlates with ideological leanings in complex ways that depend on concurrent political events such as winning or losing elections we show that posts with metaphors elicit more engagement from their audience overall even after controlling for various socio political factors such as gender and political party affiliation and finally, we demonstrate that metaphoricity is indeed the reason for increased engagement of posts, through a fine grained linguistic analysis of metaphorical vs. literal usages of 513 words across 70K posts.


PROPm Allocations of Indivisible Goods to Multiple Agents. We study the classic problem of fairly allocating a set of indivisible goods among a group of agents, and focus on the notion of approximate proportionality known as PROPm. Prior work showed that there exists an allocation that satisfies this notion of fairness for instances involving up to five agents, but fell short of proving that this is true in general. We extend this result to show that a PROPm allocation is guaranteed to exist for all instances, independent of the number of agents or goods. Our proof is constructive, providing an algorithm that computes such an allocation and, unlike prior work, the running time of this algorithm is polynomial in both the number of agents and the number of goods.


Modeling the interplay between seasonal flu outcomes and individual vaccination decisions. Seasonal influenza presents an ongoing challenge to public health. The rapid evolution of the flu virus necessitates annual vaccination campaigns, but the decision to get vaccinated or not in a given year is largely voluntary, at least in the United States, and many people decide against it. In early attempts to model these yearly flu vaccine decisions, it was often assumed that individuals behave rationally, and do so with perfect information assumptions that allowed the techniques of classical economics and game theory to be applied. However, the usual assumptions are contradicted by the emerging empirical evidence about human decision making behavior in this context. We develop a simple model of coupled disease spread and vaccination dynamics that instead incorporates experimental observations from social psychology to model annual vaccine decision making more realistically. We investigate population level effects of these new decision making assumptions, with the goal of understanding whether the population can self organize into a state of herd immunity, and if so, under what conditions. Our model agrees with established results while also revealing more subtle population level behavior, including biennial oscillations about the herd immunity threshold.


Perspective taking and Pragmatics for Generating Empathetic Responses Focused on Emotion Causes. Empathy is a complex cognitive ability based on the reasoning of others' affective states. In order to better understand others and express stronger empathy in dialogues, we argue that two issues must be tackled at the same time (i) identifying which word is the cause for the other's emotion from his or her utterance and (ii) reflecting those specific words in the response generation. However, previous approaches for recognizing emotion cause words in text require sub utterance level annotations, which can be demanding. Taking inspiration from social cognition, we leverage a generative estimator to infer emotion cause words from utterances with no word level label. Also, we introduce a novel method based on pragmatics to make dialogue models focus on targeted words in the input during generation. Our method is applicable to any dialogue models with no additional training on the fly. We show our approach improves multiple best performing dialogue agents on generating more focused empathetic responses in terms of both automatic and human evaluation.


Latexify Math Mathematical Formula Markup Revision to Assist Collaborative Editing in Math Q A Sites. Collaborative editing questions and answers plays an important role in quality control of Mathematics Stack Exchange which is a math Q A Site. Our study of post edits in Mathematics Stack Exchange shows that there is a large number of math related edits about latexifying formulas, revising LaTeX and converting the blurred math formula screenshots to LaTeX sequence. Despite its importance, manually editing one math related post especially those with complex mathematical formulas is time consuming and error prone even for experienced users. To assist post owners and editors to do this editing, we have developed an edit assistance tool, MathLatexEdit for formula latexification, LaTeX revision and screenshot transcription. We formulate this formula editing task as a translation problem, in which an original post is translated to a revised post. MathLatexEdit implements a deep learning based approach including two encoder decoder models for textual and visual LaTeX edit recommendation with math specific inference. The two models are trained on large scale historical original edited post pairs and synthesized screenshot formula pairs. Our evaluation of MathLatexEdit not only demonstrates the accuracy of our model, but also the usefulness of MathLatexEdit in editing real world posts which are accepted in Mathematics Stack Exchange.


A self consistent hybrid model of kinetic striations in low current Argon discharges. A self consistent hybrid model of standing and moving striations was developed for low current DC discharges in noble gases. We introduced the concept of surface diffusion in phase space described by a tensor diffusion in the nonlocal Fokker Planck kinetic equation for electrons in the collisional plasma. Electrons diffuse along surfaces of constant total energy between energy jumps in inelastic collisions with atoms. Numerical solutions of the 1d1u kinetic equation for electrons were obtained by two methods and coupled to ion transport and Poisson solver. We studied the dynamics of striation formation in Townsend and glow discharges in Argon gas at low discharge currents using a two level excitation ionization model and a full chemistry model, which includes stepwise and Penning ionization. Standing striations appeared in Townsend and glow discharges at low currents. Moving striations were obtained for glow discharges at currents exceeding a critical value. They originate at the anode and propagate towards the cathode. We have obtained two types of moving striations with the two level and full chemistry models, which resemble the s and p striations previously observed in the experiments. Simulations indicate that processes in the anode region could control moving striations in the positive column plasma. The developed model helps clarify the nature of standing and moving striations in DC discharges of noble gases at low discharge currents and low gas pressures.


Femtosecond Pulse Generation via an Integrated Electro Optic Time Lens. Integrated femtosecond pulse and frequency comb sources are critical components for a wide range of applications. The leading approaches for on chip pulse generation rely on mode locking inside microresonator with either third order nonlinearity or with semiconductor gain. These approaches, however, are limited in noise performance, wavelength tunability and repetition rates. Alternatively, sub picosecond pulses can be synthesized without mode locking, by modulating a continuous wave (CW) single frequency laser using a cascade of electro optic (EO) modulators. This method is particularly attractive due to its simplicity, robustness, and frequency agility but has been realized only on a tabletop using multiple discrete EO modulators and requiring optical amplifiers (to overcome large insertion losses), microwave amplifiers, and phase shifters. Here we demonstrate a chip scale femtosecond pulse source implemented on an integrated lithium niobate (LN) photonic platform18, using cascaded low loss electro optic amplitude and phase modulators and chirped Bragg grating, forming a time lens system. The device is driven by a CW distributed feedback (DFB) chip laser and controlled by a single CW microwave source without the need for any stabilization or locking. We measure femtosecond pulse trains (520 fs duration) with a 30 GHz repetition rate, flat top optical spectra with a 10 dB optical bandwidth of 12.6 nm, individual comb line powers above 0.1 milliwatt, and pulse energies of 0.54 picojoule. Our results represent a tunable, robust and low cost integrated pulsed light source with CW to pulse conversion efficiencies an order of magnitude higher than achieved with previous integrated sources. Our pulse generator can find applications from ultrafast optical measurement to networks of distributed quantum computers.


(Im)balance in the Representation of News? An Extensive Study on a Decade Long Dataset from India. (Im)balance in the representation of news has always been a topic of debate in political circles. The concept of balance has often been discussed and studied in the context of the social responsibility theory and the prestige press in the USA. While various qualitative, as well as quantitative measures of balance, have been suggested in the literature, a comprehensive analysis of all these measures across a large dataset of the post truth era comprising different popular news media houses and over a sufficiently long temporal scale in a non US democratic setting is lacking. We use this concept of balance to measure and understand the evolution of imbalance in Indian media on various journalistic metrics on a month by month basis. For this study, we amass a huge dataset of over four million political articles from India for 9 years and analyze the extent and quality of coverage given to issues and political parties in the context of contemporary influential events for three leading newspapers. We use several state of the art NLP tools to effectively understand political polarization (if any) manifesting in these articles over time. We find that two out of the three news outlets are more strongly clustered in their imbalance metrics. We also observe that only a few locations are extensively covered across all the news outlets and the situation is only slightly getting better for one of the three news outlets. Cloze tests show that the changing landscape of events get reflected in all the news outlets with border and terrorism issues dominating in around 2010 while economic aspects like unemployment, GST, demonetization, etc. became more dominant in the period 2014 2018. Further, cloze tests clearly portray the changing popularity profile of the political parties over time.


Developing and experimenting with LEO satellite constellations in OMNeT . In this paper, we present our work in designing and implementing a LEO satellite constellation simulation model1 within OMNeT and INET, which is validated by comparing the results with existing work. Our model builds upon the fundamentals of the Open Source Satellite Simulator (OS 3 ), which was ported to INET 4.3. We describe how the model was integrated on top of the INET and ported OS 3 Framework. We then experiment with the simulation model to demonstrate its viability in simulating LEO satellite constellations. This involved simulating both outdated and more recent satellite constellations, using FCC filing information, to validate latency results.


Distributed Deep Reinforcement Learning for Intelligent Traffic Monitoring with a Team of Aerial Robots. This paper studies the traffic monitoring problem in a road network using a team of aerial robots. The problem is challenging due to two main reasons. First, the traffic events are stochastic, both temporally and spatially. Second, the problem has a non homogeneous structure as the traffic events arrive at different locations of the road network at different rates. Accordingly, some locations require more visits by the robots compared to other locations. To address these issues, we define an uncertainty metric for each location of the road network and formulate a path planning problem for the aerial robots to minimize the network's average uncertainty. We express this problem as a partially observable Markov decision process (POMDP) and propose a distributed and scalable algorithm based on deep reinforcement learning to solve it. We consider two different scenarios depending on the communication mode between the agents (aerial robots) and the traffic management center (TMC). The first scenario assumes that the agents continuously communicate with the TMC to send receive real time information about the traffic events. Hence, the agents have global and real time knowledge of the environment. However, in the second scenario, we consider a challenging setting where the observation of the aerial robots is partial and limited to their sensing ranges. Moreover, in contrast to the first scenario, the information exchange between the aerial robots and the TMC is restricted to specific time instances. We evaluate the performance of our proposed algorithm in both scenarios for a real road network topology and demonstrate its functionality in a traffic monitoring system.


Characterization and Modeling of Native MOSFETs Down to 4.2 K. The extremely low threshold voltage (Vth) of native MOSFETs (Vth 0V 300K) is conducive to the design of cryogenic circuits. Previous research on cryogenic MOSFETs mainly focused on the standard threshold voltage (SVT) and low threshold voltage (LVT) MOSFETs. In this paper, we characterize native MOSFETs within the temperature range from 300K to 4.2K. The cryogenic Vth increases up to 0.25V (W L 10um 10um) and the improved subthreshold swing (SS) 14.30mV dec 4.2K. The off state current (Ioff) and the gate induced drain leakage (GIDL) effect are ameliorated greatly. The step up effect caused by the substrate charge and the transconductance peak effect caused by the energy quantization in different sub bands are also discussed. Based on the EKV model, we modified the mobility calculation equations and proposed a compact model of large size native MOSFETs suitable for the range of 300K to 4.2K. The mobility related parameters are extracted via a machine learning approach and the temperature dependences of the scattering mechanisms are analyzed. This work is beneficial to both the research on cryogenic MOSFETs modeling and the design of cryogenic CMOS circuits for quantum chips.


MCDAL Maximum Classifier Discrepancy for Active Learning. Recent state of the art active learning methods have mostly leveraged Generative Adversarial Networks (GAN) for sample acquisition however, GAN is usually known to suffer from instability and sensitivity to hyper parameters. In contrast to these methods, we propose in this paper a novel active learning framework that we call Maximum Classifier Discrepancy for Active Learning (MCDAL) which takes the prediction discrepancies between multiple classifiers. In particular, we utilize two auxiliary classification layers that learn tighter decision boundaries by maximizing the discrepancies among them. Intuitively, the discrepancies in the auxiliary classification layers' predictions indicate the uncertainty in the prediction. In this regard, we propose a novel method to leverage the classifier discrepancies for the acquisition function for active learning. We also provide an interpretation of our idea in relation to existing GAN based active learning methods and domain adaptation frameworks. Moreover, we empirically demonstrate the utility of our approach where the performance of our approach exceeds the state of the art methods on several image classification and semantic segmentation datasets in active learning setups.


Self regulating soliton domain walls in microresonators. Dissipative soliton Kerr frequency combs in microresonators have recently been demonstrated with the self injection locking process. They have the advantage of turnkey deterministic comb generation and simplifying dark soliton generation in the normal dispersion regime. Here, the formation process of dark pulses triggered by self injection locking is studied by regarding them as a pair of domain walls that connect domains having different intracavity powers. The self injection locking mechanism allows the domain walls to self regulate their position so that a wide range of dark comb states can be accessed, and the duty cycle is controlled by the feedback phase. Direct imaging of the dark pulse shape using the electro optic sampling technique is used to verify the theory. The results provide new physical insights as well as a new operational modality for this important class of nonlinear waves.


Learning multiple gaits of quadruped robot using hierarchical reinforcement learning. There is a growing interest in learning a velocity command tracking controller of quadruped robot using reinforcement learning due to its robustness and scalability. However, a single policy, trained end to end, usually shows a single gait regardless of the command velocity. This could be a suboptimal solution considering the existence of optimal gait according to the velocity for quadruped animals. In this work, we propose a hierarchical controller for quadruped robot that could generate multiple gaits (i.e. pace, trot, bound) while tracking velocity command. Our controller is composed of two policies, each working as a central pattern generator and local feedback controller, and trained with hierarchical reinforcement learning. Experiment results show 1) the existence of optimal gait for specific velocity range 2) the efficiency of our hierarchical controller compared to a controller composed of a single policy, which usually shows a single gait. Codes are publicly available.


Memory Augmented Reinforcement Learning for Image Goal Navigation. In this work, we present a memory augmented approach for image goal navigation. Earlier attempts, including RL based and SLAM based approaches have either shown poor generalization performance, or are heavily reliant on pose depth sensors. Our method is based on an attention based end to end model that leverages an episodic memory to learn to navigate. First, we train a state embedding network in a self supervised fashion, and then use it to embed previously visited states into the agent's memory. Our navigation policy takes advantage of this information through an attention mechanism. We validate our approach with extensive evaluations, and show that our model establishes a new state of the art on the challenging Gibson dataset. Furthermore, we achieve this impressive performance from RGB input alone, without access to additional information such as position or depth, in stark contrast to related work.


Accelerating Serverless Computing by Harvesting Idle Resources. Serverless computing automates fine grained resource scaling and simplifies the development and deployment of online services with stateless functions. However, it is still non trivial for users to allocate appropriate resources due to various function types, dependencies, and input sizes. Misconfiguration of resource allocations leaves functions either under provisioned or over provisioned and leads to continuous low resource utilization. This paper presents Freyr, a new resource manager (RM) for serverless platforms that maximizes resource efficiency by dynamically harvesting idle resources from over provisioned functions to under provisioned functions. Freyr monitors each function's resource utilization in real time, detects over provisioning and under provisioning, and learns to harvest idle resources safely and accelerates functions efficiently by applying deep reinforcement learning algorithms along with a safeguard mechanism. We have implemented and deployed a Freyr prototype in a 13 node Apache OpenWhisk cluster. Experimental results show that 38.8 of function invocations have idle resources harvested by Freyr, and 39.2 of invocations are accelerated by the harvested resources. Freyr reduces the 99th percentile function response latency by 32.1 compared to the baseline RMs.


GAN Inversion A Survey. GAN inversion aims to invert a given image back into the latent space of a pretrained GAN model, for the image to be faithfully reconstructed from the inverted code by the generator. As an emerging technique to bridge the real and fake image domains, GAN inversion plays an essential role in enabling the pretrained GAN models such as StyleGAN and BigGAN to be used for real image editing applications. Meanwhile, GAN inversion also provides insights on the interpretation of GAN's latent space and how the realistic images can be generated. In this paper, we provide an overview of GAN inversion with a focus on its recent algorithms and applications. We cover important techniques of GAN inversion and their applications to image restoration and image manipulation. We further elaborate on some trends and challenges for future directions.


Guaranteed epsilon optimal solutions with the linear optimizer ART3 O. The linear optimization algorithm ART3 O introduced by Chen et al. (2010) can efficiently solve large scale inverse planning problems encountered in radiation therapy by iterative projection. Its major weakness is that it cannot guarantee epsilon optimality of the final solution due to an arbitrary stopping criterion. We propose an improvement to ART3 O where the stopping criterion is based on Farkas' lemma. The same theory can be used to detect inconsistency in other projection methods as well. The proposed algorithm guarantees to find an epsilon optimal solution in finite time. The algorithm is demonstrated on numerical examples in radiation therapy.


Wind Farm Icing Loss Forecast Pertinent to Winter Extremes. The 2021 Texas power crisis has highlighted the vulnerability of the power system under wind extremes, particularly with the increasing penetration of energy resources that depend on weather conditions (e.g., wind energy). The current wind power forecast models do not effectively consider the impact of such extreme weather events. In the present study, we provide a fast and robust statistical model developed using ten years of utility scale turbine data at the Eolos Wind Energy Research Station to forecast the icing losses under such weather conditions. This model covers different cold climate impacts, including precipitation icing, frost contamination, and low temperature effect. This model has been assessed using three large scale (larger than 100 MW) wind farm data involving turbines with different capacities and from different manufacturers across multiple geographic regions. Notably, the model has been used to predict the wind power losses in entire Texas (larger than 91 of total wind installation in Texas) during 2021 Texas power crisis. The proposed model can be easily integrated into the existing wind farm and power grid operations, allowing the power system operators can develop more appropriate and pinpointed plans to balance the severe and sudden energy deficits and increase the system integrity under winter extremes.