title: Separable nonlinear least-squares parameter estimation for complex dynamic systems

abstract: Nonlinear dynamic models are widely used for characterizing functional forms of processes that govern complex biological pathway systems. Over the past decade, validation and further development of these models became possible due to data collected via high-throughput experiments using methods from molecular biology. While these data are very beneficial, they are typically incomplete and noisy, so that inferring parameter values for complex dynamic models is associated with serious computational challenges. Fortunately, many biological systems have embedded linear mathematical features, which may be exploited, thereby improving fits and leading to better convergence of optimization algorithms. In this paper, we explore options of inference for dynamic models using a novel method of {\it separable nonlinear least-squares optimization}, and compare its performance to the traditional nonlinear least-squares method. The numerical results from extensive simulations suggest that the proposed approach is at least as accurate as the traditional nonlinear least-squares, but usually superior, while also enjoying a substantial reduction in computational time.


title: Efficient Three-party Computation: An Information-theoretic Approach from Cut-and-Choose

abstract: As far as we know, the literature on secure computation from cut-and-choose has focused on achieving computational security against malicious adversaries. It is unclear whether the idea of cut-and-choose can be adapted to secure computation with information-theoretic security. In this work we explore the possibility of using cut-and-choose in information theoretic setting for secure three-party computation (3PC). Previous work on 3PC has mainly focus on the semi-honest case, and is motivated by the observation that real-word deployments of multi-party computation (MPC) seem to involve few parties. We propose a new protocol for information-theoretically secure 3PC tolerating one malicious party with cheating probability $2^{-s}$ using $s$ runs of circuit computation in the cut-and-choose paradigm. The computational cost of our protocol is essentially only a small constant worse than that of state-of-the-art 3PC protocols against a semi-honest corruption, while its communication round is greatly reduced compared to other maliciously secure 3PC protocols in information-theoretic setting.


title: Introduction to the 35th International Conference on Logic Programming Special Issue

abstract: We are proud to introduce this special issue of Theory and Practice of Logic Programming (TPLP), dedicated to the regular papers accepted for the 35th International Conference on Logic Programming (ICLP). The ICLP meetings started in Marseille in 1982 and since then constitute the main venue for presenting and discussing work in the area of logic programming. Under consideration for acceptance in TPLP.


title: The Spectral Energy Distributions of Active Galactic Nuclei

abstract: We present spectral energy distributions (SEDs) of 41 active galactic nuclei, derived from multiwavelength photometry and archival spectroscopy. All of the SEDs span at least 0.09 to 30 micron, but in some instances wavelength coverage extends into the X-ray, far-infrared and radio. For some AGNs we have fitted the measured far-infrared photometry with greybody models, while radio flux density measurements have been approximated by power-laws or polynomials. We have been able to fill some of the gaps in the spectral coverage using interpolation or extrapolation of simple models. In addition to the 41 individual AGN SEDs, we have produced 72 Seyfert SEDs by mixing SEDs of the central regions of Seyferts with galaxy SEDs. Relative to the literature, our templates have broader wavelength coverage and/or higher spectral resolution. We have tested the utility of our SEDs by using them to generate photometric redshifts for 0 < z < 6.12 AGNs in the Bootes field (selected with X-ray, IR and optical criteria) and, relative to SEDs from the literature, they produce comparable or better photometric redshifts with reduced flux density residuals.


title: Magnetic helicity in multiply connected domains

abstract: Magnetic helicity is a fundamental quantity of magnetohydrodynamics that carries topological information about the magnetic field. By `topological information', we usually refer to the linkage of magnetic field lines. For domains that are not simply connected, however, helicity also depends on the topology of the domain. In this paper, we expand the standard definition of magnetic helicity in simply connected domains to multiply connected domains in $\mathbb{R}^3$ of arbitrary topology. We also discuss how using the classic Biot-Savart operator simplifies the expression for helicity and how domain topology affects the physical interpretation of helicity.


title: Analysis of Electric Dipole Moment of $^{225}$Ra Atom using the Relativistic Normal Coupled-cluster Theory

abstract: In view of the large differences in the previous calculations of enhancement factors to the parity and time-reversal violating (P,T-odd) electric dipole moment (EDM) of $^{225}$Ra due to nuclear Schiff moment (NSM) and tensor-pseudotensor (T-PT) electron-nucleus (e-N) interactions between the relativistic coupled-cluster (RCC) theory and other many-body methods, we employ the relativistic normal coupled-cluster (RNCC) theory to explain the discrepancies. The normalization of the wave function in the RNCC theory becomes unity by construction. This feature removes the ambiguity associated with the uncertainties in calculations that could arise due to mismatch in cancellation of the normalization factor of the wave function in a truncated RCC method. Moreover, all the terms in the expression for EDM using the RNCC method naturally terminate, in contrast to the RCC approach. By taking an average of the results from two variants each of both the RCC and RNCC methods, we recommend enhancement factors to the EDM of 225Ra due to NSM as $-$6.29(1) $\times 10^{-17} |e| $cm $( |e| fm^3)$ and due to T-PT e-N coupling constant as $-$12.66(14) $\times {10^{-20} \langle \sigma_N \rangle | e | }$cm, for the nuclear Pauli spinor, $\sigma_N$. This is corroborated by analyzing the dipole polarizability ($\alpha_d$) value of $^{225}$Ra, which is obtained as 244(13) $ea_0^3$. We also compare our results for all three properties with previous calculations that employ different many-body methods. Our $\alpha_d$ value agrees very well with the results that are obtained by carrying out rigorous analyses using other variants of RCC methods.


title: Radiative parton energy loss and baryon stopping in $AA$ collisions

abstract: We study the radiative energy loss contribution to proton stopping in $AA$ collisions. The analyses is performed within the light-cone path integral approach to induced gluon emission. We have found that the radiative correction can fill in partly the midrapidity dip in the net proton rapidity distribution in $AA$ collisions at $\sqrt{s}\sim 10$ GeV. We argue that at $\sqrt{s}\sim 10$ GeV the net proton fluctuations at midrapidity may be dominated by the initial fluctuations of the proton flow, which, to a good accuracy, should be binomial.


title: Slide Reduction, Revisited---Filling the Gaps in SVP Approximation

abstract: We show how to generalize Gama and Nguyen's slide reduction algorithm [STOC '08] for solving the approximate Shortest Vector Problem over lattices (SVP). As a result, we show the fastest provably correct algorithm for $\delta$-approximate SVP for all approximation factors $n^{1/2+\varepsilon} \leq \delta \leq n^{O(1)}$. This is the range of approximation factors most relevant for cryptography.


title: Loops versus lines and the compression stiffening of cells

abstract: Both animal and plant tissue exhibit a nonlinear rheological phenomenon known as compression stiffening, or an increase in moduli with increasing uniaxial compressive strain. Does such a phenomenon exist in single cells, which are the building blocks of tissues? One expects an individual cell to compression soften since the semiflexible biopolymer-based cytoskeletal network maintains the mechanical integrity of the cell and in vitro semiflexible biopolymer networks typically compression soften. To the contrary, we find that mouse embryonic fibroblasts (mEFs) compression stiffen under uniaxial compression via atomic force microscopy (AFM) studies. To understand this finding, we uncover several potential mechanisms for compression stiffening. First, we study a single semiflexible polymer loop modeling the actomyosin cortex enclosing a viscous medium modeled as an incompressible fluid. Second, we study a two-dimensional semiflexible polymer/fiber network interspersed with area-conserving loops, which are a proxy for vesicles and fluid-based organelles. Third, we study two-dimensional fiber networks with angular-constraining crosslinks, i.e. semiflexible loops on the mesh scale. In the latter two cases, the loops act as geometric constraints on the fiber network to help stiffen it via increased angular interactions. We find that the single semiflexible polymer loop model agrees well with our AFM experiments until approximately 35% compressive strain. We also find for the fiber network with area-conserving loops model that the stress-strain curves are sensitive to the packing fraction and size distribution of the area-conserving loops, thereby creating a mechanical fingerprint across different cell types. Finally, we make comparisons between this model and experiments on fibrin networks interlaced with beads as well as discuss the tissue-scale implications of cellular compression stiffening.


title: Generalized random matrix model with additional interactions

abstract: We introduce a log-gas model that is a generalization of a random matrix ensemble with an additional interaction, whose strength depends on a parameter $\gamma$. The equilibrium density is computed by numerically solving the Riemann-Hilbert problem associated with the ensemble. The effect of the additional parameter $\gamma$ associated with the two-body interaction can be understood in terms of an effective $\gamma$-dependent single-particle confining potential.


title: Multi-phonon interactions between nitrogen-vacancy centers and nanomechanical resonators

abstract: We investigate the multi-phonon interactions in a hybrid system composed of a nitrogen-vacancy center and a mechanical resonator.We show that, through appropriate sideband engineering analogy to the Mollow or Lamb-Dicke dynamics, the enhanced nonlinear interactions can dominate the coupled system. As an example, we show the preparation of nonclassical states of the mechanical motion and explore the quantum correlation of $n$-phonon with the engineered dissipation,based on the large multi-phonon coupling strength attainable in this sideband structure. This work takes full advantage of the structure and coherence features of defect centers in diamond, and may be useful for quantum information processing.


title: A Nash-Type Fictitious Game Framework to Time-Inconsistent Stochastic Control Problems

abstract: In this paper, a Nash-type fictitious game framework is introduced to handle a time-inconsistent linear-quadratic optimal control. The Nash-type game in this framework is called fictitious as it is between the decision maker (called real player) and an auxiliary control variable (called fictitious player) with the real player and fictitious player looking for time-consistent policy and precomitted optimal policy, respectively. Namely, the fictitious-game framework is actually an auxiliary-variable-based mechanism where the fictitious player is our particular design. Noting that the real player's cost functional is revised in accordance with that of fictitious player, the equilibrium policy of real player is called an open-loop self-coordination control of original linear-quadratic problem. As a generalization, a time-inconsistent nonzero-sum stochastic linear-quadratic dynamic game is investigated, where one player is to look for precommitted optimal policy and the other player is to search time-consistent policy. Necessary and sufficient conditions are presented to ensure the existence of open-loop equilibrium of the nonzero-sum game, which resort to a set of Riccati-like equations and linear equations. By applying the developed theory of nonzero-sum game, open-loop self-coordination control of the linear-quadratic optimal control is fully characterized, and multi-period mean-variance portfolio selection is also investigated. Finally, numerical simulations are presented, which show the efficiency of the proposed fictitious-game framework.


title: Evolution equations on time-dependent intervals

abstract: We study initial boundary value problems for linear evolution partial differential equations (PDEs) posed on a time-dependent interval $l_1(t)<x<l_2(t)$, $0<t<T$, where $l_1(t)$ and $l_2(t)$ are given, real, differentiable functions, and $T$ is an arbitrary constant. For such problems, we show how to characterise the unknown boundary values in terms of the given initial and boundary conditions. As illustrative examples we consider the heat equation and the linear Schr\"{o}dinger equation. In the first case, the unknown Neumann boundary values are expressed in terms of the Dirichlet boundary values and of the initial value through the unique solution of a system of two linear integral equations with explicit kernels. In the second case, a similar result can be proved but only for a more restrictive class of boundary curves.}


title: On the integrability of the Abel and of the extended Li\'{e}nard equations

abstract: We present some exact integrability cases of the extended Li\'{e}nard equation $y^{\prime \prime }+f\left( y\right) \left(y^{\prime }\right)^{n}+k\left( y\right) \left(y^{\prime }\right)^{m}+g\left(y\right) y^{\prime }+h\left( y\right) =0$, with $n>0$ and $m>0$ arbitrary constants, while $f(y)$, $k(y)$, $g(y)$, and $h(y)$ are arbitrary functions. The solutions are obtained by transforming the equation Li\'{e}nard equation to an equivalent first kind first order Abel type equation given by $\frac{dv}{dy} =f\left( y\right) v^{3-n}+k\left( y\right) v^{3-m}+g\left( y\right) v^{2}+h\left( y\right) v^{3}$, with $v=1/y^{\prime }$. As a first step in our study we obtain three integrability cases of the extended quadratic-cubic Li\'{e}nard equation, corresponding to $n=2$ and $m=3$, by assuming that particular solutions of the associated Abel equation are known. Under this assumption the general solutions of the Abel and Li\'{e}nard equations with coefficients satisfying some differential conditions can be obtained in an exact closed form. With the use of the Chiellini integrability condition, we show that if a particular solution of the Abel equation is known, the general solution of the extended quadratic cubic Li\'{e}nard equation can be obtained by quadratures. The Chiellini integrability condition is extended to generalized Abel equations with $g(y)\equiv 0$ and $h(y)\equiv 0$, and arbitrary $n$ and $m$, thus allowing to obtain the general solution of the corresponding Li\'{e}nard equation. The application of the generalized Chiellini condition to the case of the reduced Riccati equation is also considered.


title: Learning to Explore in Motion and Interaction Tasks

abstract: Model free reinforcement learning suffers from the high sampling complexity inherent to robotic manipulation or locomotion tasks. Most successful approaches typically use random sampling strategies which leads to slow policy convergence. In this paper we present a novel approach for efficient exploration that leverages previously learned tasks. We exploit the fact that the same system is used across many tasks and build a generative model for exploration based on data from previously solved tasks to improve learning new tasks. The approach also enables continuous learning of improved exploration strategies as novel tasks are learned. Extensive simulations on a robot manipulator performing a variety of motion and contact interaction tasks demonstrate the capabilities of the approach. In particular, our experiments suggest that the exploration strategy can more than double learning speed, especially when rewards are sparse. Moreover, the algorithm is robust to task variations and parameter tuning, making it beneficial for complex robotic problems.


title: Improving cosmological parameter estimation with the future 21 cm observation from SKA

abstract: Future observations of 21 cm emission from neutral hydrogen survey will become a promising approach to probe the large scale structure of Universe. In this paper, we investigate the impacts of Square Kilometer Array (SKA) 21 cm observation on the estimation of cosmological parameters. We use the simulated data of the baryonic acoustic oscillation (BAO) measurements based on the future SKA experiment with the intensity mapping (IM) technique to do the analysis. For the current observations, we use the latest cosmic microwave background (CMB) observation from {\it Planck} 2018, the optical BAO measurements, and the Type Ia supernovae (SN) observation (Pantheon compilation). We find that the SKA mock data could break the degeneracy between the matter density and the Hubble constant, further improving the cosmological constraints to a great extent. We also find that the constraint on the equation of state parameters of dark energy could be significantly improved by including the SKA mock data into the cosmological global fit.


title: Subdivisions of digraphs in tournaments

abstract: We show that for every positive integer $k$, any tournament with minimum out-degree at least $(2+o(1))k^2$ contains a subdivision of the complete directed graph on $k$ vertices, which is best possible up to a factor of $8$. This may be viewed as a directed analogue of a theorem proved by Bollob\'as and Thomason, and independently by Koml\'os and Szemer\'edi, concerning subdivisions of cliques in graphs with sufficiently high average degree. We also consider the following problem: given $k$, what is the smallest positive integer $f(k)$ such that any $f(k)$-vertex tournament contains a $1$-subdivision of the transitive tournament on $k$ vertices? We show that $f(k)= O\left (k^2\log^3 k\right)$ which is best possible up to the logarithmic factors.


title: Unsupervised Stemming based Language Model for Telugu Broadcast News Transcription

abstract: In Indian Languages , native speakers are able to understand new words formed by either combining or modifying root words with tense and / or gender. Due to data insufficiency, Automatic Speech Recognition system (ASR) may not accommodate all the words in the language model irrespective of the size of the text corpus. It also becomes computationally challenging if the volume of the data increases exponentially due to morphological changes to the root word. In this paper a new unsupervised method is proposed for a Indian language: Telugu, based on the unsupervised method for Hindi, to generate the Out of Vocabulary (OOV) words in the language model. By using techniques like smoothing and interpolation of pre-processed data with supervised and unsupervised stemming, different issues in language model for Indian language: Telugu has been addressed. We observe that the smoothing techniques Witten-Bell and Kneser-Ney perform well when compared to other techniques on pre-processed data from supervised learning. The ASRs accuracy is improved by 0.76% and 0.94% with supervised and unsupervised stemming respectively.


title: Automatic acute ischemic stroke lesion segmentation using semi-supervised learning

abstract: Ischemic stroke is a common disease in the elderly population, which can cause long-term disability and even death. However, the time window for treatment of ischemic stroke in its acute stage is very short. To fast localize and quantitively evaluate the acute ischemic stroke (AIS) lesions, many deep-learning-based lesion segmentation methods have been proposed in the literature, where a deep convolutional neural network (CNN) was trained on hundreds of fully labeled subjects with accurate annotations of AIS lesions. Despite that high segmentation accuracy can be achieved, the accurate labels should be annotated by experienced clinicians, and it is therefore very time-consuming to obtain a large number of fully labeled subjects. In this paper, we propose a semi-supervised method to automatically segment AIS lesions in diffusion weighted images and apparent diffusion coefficient maps. By using a large number of weakly labeled subjects and a small number of fully labeled subjects, our proposed method is able to accurately detect and segment the AIS lesions. In particular, our proposed method consists of three parts: 1) a double-path classification net (DPC-Net) trained in a weakly-supervised way is used to detect the suspicious regions of AIS lesions; 2) a pixel-level K-Means clustering algorithm is used to identify the hyperintensive regions on the DWIs; and 3) a region-growing algorithm combines the outputs of the DPC-Net and the K-Means to obtain the final precise lesion segmentation. In our experiment, we use 460 weakly labeled subjects and 15 fully labeled subjects to train and fine-tune the proposed method. By evaluating on a clinical dataset with 150 fully labeled subjects, our proposed method achieves a mean dice coefficient of 0.642, and a lesion-wise F1 score of 0.822.


title: Fallback Strategies in Operation Control of Microgrids with Communication Failures

abstract: This paper proposes a model predictive control (MPC)-based energy management system to address communication failures in an islanded microgrid (MG). The energy management system relies on a communication network to monitor and control power generation in the units. A communication failure to a unit inhibits the ability of the MPC to change the power of the unit. However, this unit is electrically connected to the network and ignoring this aspect could have adverse effect in the operation of the microgrid. Therefore, this paper considers an MPC design that includes the electrical connectivity of the unit during communication failures. This paper also highlights the benefit of including the lower layer control behaviour of the MG to withstand communication failures. The proposed approaches are validated with a case study.


title: Deep Triplet Neural Networks with Cluster-CCA for Audio-Visual Cross-modal Retrieval

abstract: Cross-modal retrieval aims to retrieve data in one modality by a query in another modality, which has been a very interesting research issue in the field of multimedia, information retrieval, and computer vision, and database. Most existing works focus on cross-modal retrieval between text-image, text-video, and lyrics-audio.Little research addresses cross-modal retrieval between audio and video due to limited audio-video paired datasets and semantic information. The main challenge of audio-visual cross-modal retrieval task focuses on learning joint embeddings from a shared subspace for computing the similarity across different modalities, where generating new representations is to maximize the correlation between audio and visual modalities space. In this work, we propose a novel deep triplet neural network with cluster canonical correlation analysis(TNN-C-CCA), which is an end-to-end supervised learning architecture with audio branch and video branch.We not only consider the matching pairs in the common space but also compute the mismatching pairs when maximizing the correlation. In particular, two significant contributions are made: i) a better representation by constructing deep triplet neural network with triplet loss for optimal projections can be generated to maximize correlation in the shared subspace. ii) positive examples and negative examples are used in the learning stage to improve the capability of embedding learning between audio and video. Our experiment is run over 5-fold cross-validation, where average performance is applied to demonstrate the performance of audio-video cross-modal retrieval. The experimental results achieved on two different audio-visual datasets show the proposed learning architecture with two branches outperforms existing six CCA-based methods and four state-of-the-art based cross-modal retrieval methods.


title: Personalized Music Recommendation with Triplet Network

abstract: Since many online music services emerged in recent years so that effective music recommendation systems are desirable. Some common problems in recommendation system like feature representations, distance measure and cold start problems are also challenges for music recommendation. In this paper, I proposed a triplet neural network, exploiting both positive and negative samples to learn the representation and distance measure between users and items, to solve the recommendation task.


title: Permutation Matrices, Their Discrete Derivatives and Extremal Properties

abstract: For a permutation $\pi$, and the corresponding permutation matrix, we introduce the notion of {\em discrete derivative}, obtained by taking differences of successive entries in $\pi$. We characterize the possible derivatives of permutations, and consider questions for permutations with certain properties satisfied by the derivative. For instance, we consider permutations with distinct derivatives, and the relationship to so-called Costas arrays.


title: Permutation Matrix Representation Quantum Monte Carlo

abstract: We present a quantum Monte Carlo algorithm for the simulation of general quantum and classical many-body models within a single unifying framework. The algorithm builds on a power series expansion of the quantum partition function in its off-diagonal terms and is both parameter-free and Trotter error-free. In our approach, the quantum dimension consists of products of elements of a permutation group. As such, it allows for the study of a very wide variety of models on an equal footing. To demonstrate the utility of our technique, we use it to clarify the emergence of the sign problem in the simulations of non-stoquastic physical models. We showcase the flexibility of our algorithm and the advantages it offers over existing state-of-the-art by simulating transverse-field Ising model Hamiltonians and comparing the performance of our technique against that of the stochastic series expansion algorithm. We also study a transverse-field Ising model augmented with randomly chosen two-body transverse-field interactions.


title: Congruences in character tables of symmetric groups

abstract: If $\lambda$ and $\mu$ are two non-empty Young diagrams with the same number of squares, and $\boldsymbol\lambda$ and $\boldsymbol\mu$ are obtained by dividing each square into $d^2$ congruent squares, then the corresponding character value $\chi_{\boldsymbol\lambda}(\boldsymbol\mu)$ is divisible by $d!$.


title: Baryogenesis in the paradigm of quintessential inflation

abstract: We explore the possibility of baryogenesis in the framework of quintessential inflation. We focus on the model independent features of the underlying paradigm and demonstrate that the required baryon asymmetry can successfully be generated in this scenario. To this effect, we use the effective field theory framework with desired terms in the Lagrangian necessary to mimic baryon number violation \textit{\`{a} la} spontaneous baryogenesis which can successfully evade Sakharov's requirement allowing us to generate the observed baryon asymmetry in the equilibrium process. Our estimates are independent of the underlying physical process responsible for baryon number violation. The underlying framework of quintessential inflation essentially includes the presence of kinetic regime after inflation which gives rise to blue spectrum of gravitational wave background at high frequencies. In addition to baryogenesis, we discuss the prospects of detection of relic gravitational wave background, in the future proposed missions, sticking to model independent treatment.


title: Non-particulate Klein-Gordon modes formed by inflation

abstract: In a full solution for a scalar quantum field coupled to an accelerating isotropic universe, all constituent non-autonomous modes of elementary excitation cease to oscillate and become unstable at a discrete sequence of times. After canonical quantization the time-frozen Hamiltonian has eigenstates that can be viewed as neither particles nor oscillatory radiation. Under standard canonical quantization, the Hamiltonian has a natural time dependent partitioning into a light component and a dark component. Under equipartition of energy, the energy density of the dark component remains constant. The dark component consists of a finite number of low wave -number repulsive units with time varying force constant. Although these unstable modes have no ground state of minimum energy there exist only finitely many of them and so their quantum Hilbert space representation, by the Stone-von Neumann theorem, is unique up to unitary equivalence. The remaining infinite number of stable modes still have a unique Fock-Cook representation and so overall there is still a preferred physical representation.


title: Audio-Visual Embedding for Cross-Modal MusicVideo Retrieval through Supervised Deep CCA

abstract: Deep learning has successfully shown excellent performance in learning joint representations between different data modalities. Unfortunately, little research focuses on cross-modal correlation learning where temporal structures of different data modalities, such as audio and video, should be taken into account. Music video retrieval by given musical audio is a natural way to search and interact with music contents. In this work, we study cross-modal music video retrieval in terms of emotion similarity. Particularly, audio of an arbitrary length is used to retrieve a longer or full-length music video. To this end, we propose a novel audio-visual embedding algorithm by Supervised Deep CanonicalCorrelation Analysis (S-DCCA) that projects audio and video into a shared space to bridge the semantic gap between audio and video. This also preserves the similarity between audio and visual contents from different videos with the same class label and the temporal structure. The contribution of our approach is mainly manifested in the two aspects: i) We propose to select top k audio chunks by attention-based Long Short-Term Memory (LSTM)model, which can represent good audio summarization with local properties. ii) We propose an end-to-end deep model for cross-modal audio-visual learning where S-DCCA is trained to learn the semantic correlation between audio and visual modalities. Due to the lack of music video dataset, we construct 10K music video dataset from YouTube 8M dataset. Some promising results such as MAP and precision-recall show that our proposed model can be applied to music video retrieval.


title: Noble gases in carbonate melts: constraints on the solubility and the surface tension by molecular dynamics simulation

abstract: Although they are rare elements in the Earth's mantle, noble gases (NG) owe to their strongly varying masses contrasting physical behaviors making them important geochemical tracers. When partial melting occurs at depth, the partitioning of NGs between phases is controlled by a distribution coefficient that can be determined from the solubility of the NGs in each phase. Here we report quantitative calculations of the solubility of He, Ne, Ar and Xe in carbonate melts based on molecular dynamics simulations. The NG solubilities are first calculated in K2CO3-CaCO3 mixtures at 1 bar and favorably compared to the only experimental data available to date. Then we investigate the effect of pressure (up to 6 GPa), focusing on two melt compositions: a dolomitic one and a natrocarbonatitic one (modeling the lava of Ol Doinyo Lengai). The solubility decreases with the amount of alkaline-earth cation in the melt and with the size of the noble gas. In the natrocarbonatitic melt, Henry's law is fulfilled at low pressures (up to ~ 0.1 GPa). At higher pressures the solubility levels off or even starts to diminish smoothly. In contrast, in molten dolomite the effect of pressure is negligible on the studied P range (3-6 GPa). At the pressures of the Earth's mantle, the solubilities of noble gases in carbonate melts are still of the same order of magnitude as the ones in molten silicates (~ 1-10 mol%). This suggests that carbonatitic melts at depth are not preferential carriers of noble gases, even if the dependence with the melt composition is not negligible and has to be evaluated on a case-by-case basis. Finally we evaluate the surface tension at the interface between carbonate melts and noble gases and its evolution with pressure. Whatever the composition of the melt and of the NG phase, the surface tension increases when P increases from 0 to 6 GPa.


title: Intrinsic area near the origin for self-similar growth-fragmentations and related random surfaces

abstract: We study the behaviour of a natural measure defined on the leaves of the genealogical tree of some branching processes, namely self-similar growth-fragmentation processes. Each particle, or cell, is attributed a positive mass that evolves in continuous time according to a positive self-similar Markov process and gives birth to children at negative jumps events. We are interested in the asymptotics of the mass of the ball centered at the root, as its radius decreases to $0$. We obtain the almost sure behaviour of this mass when the Eve cell starts with a strictly positive size. This differs from the situation where the Eve cell grows indefinitely from size 0. In this case, we show that, when properly rescaled, the mass of the ball converges in distribution towards a non-degenerate random variable. We then derive bounds describing the almost sure behaviour of the rescaled mass. Those results are applied to certain random surfaces, exploiting the connection between growth-fragmentations and random planar maps obtained in Bertoin et al. [6]. This allows us to extend a result of Le Gall [24] on the volume of a free Brownian disk close to its boundary, to a larger family of stable disks. The upper bound of the mass of a typical ball in the Brownian map is refined, and we obtain a lower bound as well.


title: Bi-cross validation for estimating spectral clustering hyper parameters

abstract: One challenge impeding the analysis of terabyte scale x-ray scattering data from the Linac Coherent Light Source LCLS, is determining the number of clusters required for the execution of traditional clustering algorithms. Here we demonstrate that previous work using bi-cross validation (BCV) to determine the number of singular vectors directly maps to the spectral clustering problem of estimating both the number of clusters and hyper parameter values. These results indicate that the process of estimating the number of clusters should not be divorced from the process of estimating other hyper parameters. Applying this method to LCLS x-ray scattering data enables the identification of dropped shots without manually setting boundaries on detector fluence and provides a path towards identifying rare and anomalous events.


title: Show Me Your Account: Detecting MMORPG Game Bot Leveraging Financial Analysis with LSTM

abstract: With the rapid growth of MMORPG market, game bot detection has become an essential task for maintaining stable in-game ecosystem. To classify bots from normal users, detection methods are proposed in both game client and server-side. Among various classification methods, data mining method in server-side captured unique characteristics of bots efficiently. For features used in data mining, behavioral and social actions of character are analyzed with numerous algorithms. However, bot developers can evade the previous detection methods by changing bot's activities continuously. Eventually, overall maintenance cost increases because the selected features need to be updated along with the change of bot's behavior. To overcome this limitation, we propose improved bot detection method with financial analysis. As bot's activity absolutely necessitates the change of financial status, analyzing financial fluctuation effectively captures bots as a key feature. We trained and tested model with actual data of Aion, a leading MMORPG in Asia. Leveraging that LSTM efficiently recognizes time-series movement of data, we achieved meaningful detection performance. Further on this model, we expect sustainable bot detection system in the near future.


title: Groups with Spanier-Whitehead duality

abstract: Building on work by Kasparov, we study the notion of Spanier-Whitehead K-duality for a discrete group. It is defined as duality in the KK-category between two C*-algebras which are naturally attached to the group, namely the reduced group C*-algebra and the crossed product for the group action on the universal example for proper actions. We compare this notion to the Baum-Connes conjecture by constructing duality classes based on two methods: the standard "gamma element" technique, and the more recent approach via cycles with property gamma. As a result of our analysis, we prove Spanier-Whitehead duality for a large class of groups, including Bieberbach's space groups, groups acting on trees, and lattices in Lorentz groups.


title: Weyl fermions in a non-abelian gauge background and trace anomalies

abstract: We study the trace and chiral anomalies of Weyl fermions in a non-abelian gauge background in four dimensions. Using a Pauli-Villars regularization we identify the trace anomaly, proving that it can be cast in a gauge invariant form, even in the presence of the non-abelian chiral anomaly, that we rederive to check the consistency of our methods. In particular, we find that the trace anomaly does not contain any parity-odd topological contribution, whose presence has been debated in the recent literature.


title: Polynomial analogues of restricted multicolor b-ary partition functions

abstract: Given an integer base $b\geq 2$, a number $\rho\geq 1$ of colors, and a finite sequence $\Lambda=(\lambda_1,\ldots,\lambda_\rho)$ of positive integers, we introduce the concept of a $\Lambda$-restricted $\rho$-colored $b$-ary partition of an integer $n\geq 1$. We also define a sequence of polynomials in $\lambda_1+\cdots+\lambda_\rho$ variables, and prove that the $n$th polynomial characterizes all $\Lambda$-restricted $\rho$-colored $b$-ary partitions of $n$. In the process we define a recurrence relation for the polynomials in question, obtain explicit formulas and identify a factorization theorem.


title: Factoring Catalan numbers

abstract: The paper describes a prime factorization of the Catalan numbers. Odd prime factors are distributed in layers in accordance with Legendre's formula. The content of each layer is a network of the intervals, Chebyshev's Segments. The primes of Segment are not calculated and are selected on the basis of its bounds. Layers contain non-repeated primes. Repeated factors are formed when primes are duplicated among different layers. The paper slightly modifies Kummer's theorem for the selection of individual prime factors, also starting from the boundaries of Segments. In conclusion, the reader is offered a software service for factorization of the Catalan number with index up to $10^8$


title: Optimization-based Settingless Algorithm Combining Protection and Fault Identification

abstract: The demand for faster protection algorithms is growing due to the increasingly faster dynamics in the system. The majority of existing algorithms require empirically selected set-points, which may reduce sensitivity to internal faults and cause security problems. This paper addresses these challenges by proposing a settingless time-domain unit protection algorithm for medium-voltage lines. The main idea of the algorithm is to identify which model of a protected line, i.e. healthy or with an internal fault, is more consistent with the input measurements. This is done by solving a number of small-scale convex optimization problems, which at the same time determine the characteristics of an internal fault that best fit the measurements. Thus, the proposed algorithm merges protection, fault location and fault type identification functionalities. The algorithm's performance is extensively tested on a grid model in MATLAB Simulink for different types of generation and grid operating conditions. The results demonstrate that the algorithm can operate quickly and reliably, and accurately estimate fault characteristics even in the presence of noisy measurements and uncertain line parameters.


title: Strong-coupling limit of the driven dissipative light-matter interaction

abstract: We approach the strong-coupling thermodynamic limit in the response of the open driven Jaynes-Cummings (JC) oscillator. We do so by highlighting the role of quantum fluctuations against the semiclassical response in three distinct regimes of operation. We begin by demonstrating the persistence of photon blockade, predicted in [H. J. Carmichael, Phys. Rev. X 5, 031028 (2015)], as a manifestation of the inherently-quantum and nonlinear JC spectrum revealed for vanishing dissipation. We then proceed to discuss the importance of bistability, which, despite being present in photon blockade, is able to provide an alignment between the semiclassical nonlinearity and quantum dynamics only for a driving amplitude having the same order of magnitude as the light-matter coupling strength. This resolution brings us to the critical point of the well-known quantum phase transition of second order on resonance, where the quantum and semiclassical pictures are once more contrasted for a varying participation of the two coherent interactions when going through the collapse of the quasi-energy spectrum.


title: A $\mu$-$\tau$-philic Higgs doublet confronted with the muon g-2, $\tau$ decays and LHC data

abstract: We propose a two-Higgs-doublet model in which one Higgs doublet has the same interactions with fermions as the SM, and another Higgs doublet only has the $\mu$-$\tau$ LFV interactions. Assuming that the Yukawa matrices are real and symmetric, we impose various relevant theoretical and experimental constraints, and find that the excesses of muon $g-2$ and lepton flavour universality in the $\tau$ decays can be simultaneously explained in the region of small mass splittings between the heavy CP-even Higgs and the CP-odd Higgs ($m_A > m_H$). The multi-lepton event searches at the LHC can sizably reduce the mass ranges of extra Higgses, and $m_H$ is required to be larger than 560 GeV.


title: Studies of Star-forming Complexes in the Galaxies NGC 628, NGC 2976, and NGC 3351

abstract: We analyze parameters of the interstellar matter emission in star-forming complexes in the high metallicity galaxies NGC~628, NGC~2976, and NGC~3351, which have different morphological types. The relation between H$\alpha$ emission and emission in CO and HI lines is considered along with the relation between H$\alpha$ emission and dust emission in the infrared range (IR). The fluxes and surface brightnesses in the UV and IR correlate well with H$\alpha$ emission. The HI emission also correlates well with H$\alpha$, while the correlation between the CO and H$\alpha$ emission is much less prominent. The ratio of the fluxes at 8 and 24 $\mu$m decreases with increasing H$\alpha$ flux. This may be due to changes in the properties of the dust ensemble (a decrease in the mass fraction of polycyclic aromatic hydrocarbons) or to changes in excitation conditions. Analysis of the kinematics of the CO lines shows that the CO flux grows with increasing velocity scatter $\Delta V$ when $\Delta V\lesssim70$~km/s. Preliminary evidence for the existence of star-forming complexes with higher values of $\Delta V$ is presented, and the increase in the velocity scatter is accompanied by a decrease in the CO line luminosity of the complex.


title: Generalised thresholding of hidden variable network models with scale-free property

abstract: The hidden variable formalism (based on the assumption of some intrinsic node parameters) turned out to be a remarkably efficient and powerful approach in describing and analyzing the topology of complex networks. Owing to one of its most advantageous property - namely proven to be able to reproduce a wide range of different degree distribution forms - it has become a standard tool for generating networks having the scale-free property. One of the most intensively studied version of this model is based on a thresholding mechanism of the exponentially distributed hidden variables associated to the nodes (intrinsic vertex weights), which give rise to the emergence of a scale-free network where the degree distribution $p(k)\sim k^{-\gamma}$ is decaying with an exponent of $\gamma =2$. Here we propose a generalization and modification of this model by extending the set of connection probabilities and hidden variable distributions that lead to the aforementioned degree distribution, and analyze the conditions leading to the above behavior analytically. In addition, we propose a relaxation of the hard threshold in the connection probabilities, which opens up the possibility for obtaining sparse scale free networks with arbitrary scaling exponent.


abstract: New families of non-parity-time-symmetric complex potentials with all-real spectra are derived by the supersymmetry method and the pseudo-Hermiticity method. With the supersymmetry method, we find families of non-parity-time-symmetric complex partner potentials which share the same spectrum as base potentials with known real spectra, such as the (complex) Wadati potentials. Different from previous supersymmetry derivations of potentials with real spectra, our derivation does not utilize discrete eigenmodes of base potentials. As a result, our partner potentials feature explicit analytical expressions which contain free functions. With the pseudo-Hermiticity method, we derive a new class of non-parity-time-symmetric complex potentials with free functions and constants, whose eigenvalues appear as conjugate pairs. This eigenvalue symmetry forces the spectrum to be all-real for a wide range of choices of these functions and constants in the potential. Tuning these free functions and constants, phase transition can also be induced, where conjugate pairs of complex eigenvalues emerge in the spectrum.


title: Hardware-efficient quantum algorithm for the simulation of open-system dynamics and thermalisation

abstract: The quantum open-system simulation is an important category of quantum simulation. By simulating the thermalisation process at the zero temperature, we can solve the ground-state problem of quantum systems. To realise the open-system evolution on the quantum computer, we need to encode the environment using qubits. However, usually the environment is much larger than the system, i.e. numerous qubits are required if the environment is directly encoded. In this paper, we propose a way to simulate open-system dynamics by reproducing reservoir correlation functions using a minimised Hilbert space. In this way, we only need a small number of qubits to represent the environment. To simulate the $n$-th-order expansion of the time-convolutionless master equation by reproducing up to $n$-time correlation functions, the number of qubits representing the environment is $\sim \lfloor \frac{n}{2} \rfloor \log_2(N_\omega N_\beta)$. Here, $N_\omega$ is the number of frequencies in the discretised environment spectrum, and $N_\beta$ is the number of terms in the system-environment interaction. By reproducing two-time correlation functions, i.e. taking $n = 2$, we can simulate the Markovian quantum master equation. In our algorithm, the environment on the quantum computer could be even smaller than the system.


title: A note on the topological slice genus of satellite knots

abstract: This paper presents evidence supporting the surprising conjecture that in the topological category the slice genus of a satellite knot $P(K)$ is bounded above by the sum of the slice genera of $K$ and $P(U)$. Our main result establishes this conjecture for a variant of the topological slice genus, the $\mathbb{Z}$-slice genus. As an application, we show that the $(n,1)$-cable of any 3-genus 1 knot (e.g. the figure 8 or trefoil knot) has topological slice genus at most 1. Further, we show that the lower bounds on the slice genus coming from the Tristram-Levine and Casson-Gordon signatures cannot be used to disprove the conjecture. Notably, the conjectured upper bound does not involve the algebraic winding number of the pattern $P$. This stands in stark contrast with the smooth category, where for example there are many genus 1 knots whose $(n,1)$-cables have arbitrarily large smooth 4-genera.


title: Large-Scale Traffic Signal Control Using a Novel Multi-Agent Reinforcement Learning

abstract: Finding the optimal signal timing strategy is a difficult task for the problem of large-scale traffic signal control (TSC). Multi-Agent Reinforcement Learning (MARL) is a promising method to solve this problem. However, there is still room for improvement in extending to large-scale problems and modeling the behaviors of other agents for each individual agent. In this paper, a new MARL, called Cooperative double Q-learning (Co-DQL), is proposed, which has several prominent features. It uses a highly scalable independent double Q-learning method based on double estimators and the UCB policy, which can eliminate the over-estimation problem existing in traditional independent Q-learning while ensuring exploration. It uses mean field approximation to model the interaction among agents, thereby making agents learn a better cooperative strategy. In order to improve the stability and robustness of the learning process, we introduce a new reward allocation mechanism and a local state sharing method. In addition, we analyze the convergence properties of the proposed algorithm. Co-DQL is applied on TSC and tested on a multi-traffic signal simulator. According to the results obtained on several traffic scenarios, Co- DQL outperforms several state-of-the-art decentralized MARL algorithms. It can effectively shorten the average waiting time of the vehicles in the whole road system.


title: Moderate deviations of density-dependent Markov chains

abstract: The density-dependent Markov chain (DDMC) introduced in \cite{Kurtz1978} is a continuous time Markov process applied in fields such as epidemics, chemical reactions and so on. In this paper, we give moderate deviation principles of paths of DDMC under some generally satisfied assumptions. The proofs for the lower and upper bounds of our main result utilize an exponential martingale and a generalized version of Girsanov's theorem. The exponential martingale is defined according to the generator of DDMC.


title: Using a 3D printer for 2D beam profile measurements in proton radiotherapy

abstract: We have used an inexpensive 3D printer mounting an inexpensive pre-irradiated silicon diode to measure the 2D dose distribution of a small (approximately 10mm FWHM) proton beam. z was measured with high resolution whereas x was changed on a 2mm grid for a total of 12 scans and 2 repeats in 20 minutes. The peak dose to the diode was approximately 530 cGy/s, and no degradation in diode sensitivity was observed. We present beam profiles and a 2D beam contour derived from these data. Most of the scan time was spent changing x by hand, and a more appropriate 3D printer along with an improved setup and printer control software would reduce that to about 1 minute for an end-to-end absolute 2D dose measurement. We propose extending this to 3D (2D transverse plus depth) for appropriate small beams by adding a small water tank.


title: Comparing Convective Self-Aggregation in Idealized Models to Observed Moist Static Energy Variability near the Equator

abstract: Idealized convection-permitting simulations of radiative-convective equilibrium (RCE) have become a popular tool for understanding the physical processes leading to horizontal variability of tropical water vapor and rainfall. However, the applicability of idealized simulations to nature is still unclear given that important processes are typically neglected, such as lateral vapor advection by extratropical intrusions, or interactive ocean coupling. Here, we exploit spectral analysis to compactly summarize the multi-scale processes supporting convective aggregation. By applying this framework to high-resolution reanalysis data and satellite observations in addition to idealized simulations, we compare convective-aggregation processes across horizontal scales and data sets. The results affirm the validity of the RCE simulations as an analogy to the real world. Column moist static energy tendencies share similar signs and scale-selectivity in convection-permitting models and observations: Radiation increases variance at wavelengths above 1,000km, while advection damps variance across wavelengths, and surface fluxes mostly reduce variance between 1,000km and 10,000km.


title: On the monomial reduction number of a monomial ideal in $K[x,y]$

abstract: The reduction number of monomial ideals in the polynomial $K[x,y]$ is studied. We focus on ideals $I$ for which $J=(x^a,y^b)$ is a reduction ideal. The computation of the reduction number amounts to solve linear inequalities. In some special cases the reduction number can be explicitly computed.


title: Wavelet transforms associated with the Index Whittaker transform

abstract: In this paper, we exploit the theory of convolution of index Whittaker transform for study of continuous and discrete Index Whittaker wavelet transform and discuss some of its basic properties. Certain boundedness, Plancherel as well as reconstruction formula for the continuous Index Whittaker Wavelet Transform (CIWWT) are obtain and Finally we discuss the discrete version of index Whittaker wavelet transform.


title: Frequency selection in a gravitationally stretched capillary jet in the jetting regime

abstract: A capillary jet falling under the effect of gravity continuously stretches while thinning downstream. We report here the effect of external periodic forcing on such a spatially varying jet in the jetting regime. Surprisingly, the optimal forcing frequency producing the most unstable jet is found to be highly dependent on the forcing amplitude. Taking benefit of the one-dimensional Eggers & Dupont (J. Fluid Mech., vol. 262, 1994, 205-221) equations, we investigate the case through nonlinear simulations and linear stability analysis. In the local framework the WKBJ formalism, established for weakly non-parallel flows, fails to capture the nonlinear simulation results quantitatively. However in the global framework, the resolvent analysis supplemented by a simple approximation of the required response norm inducing breakup, is shown to correctly predict the optimal forcing frequency at a given forcing amplitude and the resulting jet breakup length. The results of the resolvent analysis are found to be in good agreement with those of the nonlinear simulations.


title: Implementation of Morse-Witten theory for a polydisperse wet 2D foam simulation

abstract: The Morse-Witten theory (D. Morse and T. Witten, EPL 22 (1993) 549-555) provides a formulation for the inter-bubble forces and corresponding deformations in a liquid foam, accurate in the limit of high liquid fraction. Here we show how the theory may be applied in practice, including allowing for polydispersity in the bubble sizes. The resulting equilibrated 2D structures are consistent with direct calculations, within the limitations of the theory. The path to developing a 3D model is outlined for future work.


title: The edge ideal of a graph and its splitting graphs

abstract: We introduce and study the concept which we call the splitting of a graph and compare algebraic properties of the edge ideals of graphs and those of their splitting graphs.


title: Modeling Engagement Dynamics of Online Discussions using Relativistic Gravitational Theory

abstract: Online discussions are valuable resources to study user behaviour on a diverse set of topics. Unlike previous studies which model a discussion in a static manner, in the present study, we model it as a time-varying process and solve two inter-related problems -- predict which user groups will get engaged with an ongoing discussion, and forecast the growth rate of a discussion in terms of the number of comments. We propose RGNet (Relativistic Gravitational Nerwork), a novel algorithm that uses Einstein Field Equations of gravity to model online discussions as `cloud of dust' hovering over a user spacetime manifold, attracting users of different groups at different rates over time. We also propose GUVec, a global user embedding method for an online discussion, which is used by RGNet to predict temporal user engagement. RGNet leverages different textual and network-based features to learn the dust distribution for discussions. We employ four baselines -- first two using LSTM architecture, third one using Newtonian model of gravity, and fourth one using a logistic regression adopted from a previous work on engagement prediction. Experiments on Reddit dataset show that RGNet achieves 0.72 Micro F1 score and 6.01% average error for temporal engagement prediction of user groups and growth rate forecasting, respectively, outperforming all the baselines significantly. We further employ RGNet to predict non-temporal engagement -- whether users will comment to a given post or not. RGNet achieves 0.62 AUC for this task, outperforming existing baseline by 8.77% AUC.


title: Predicting Critical Transitions in Multiscale Dynamical Systems Using Reservoir Computing

abstract: We study the problem of predicting rare critical transition events for a class of slow-fast nonlinear dynamical systems. The state of the system of interest is described by a slow process, whereas a faster process drives its evolution and induces critical transitions. By taking advantage of recent advances in reservoir computing, we present a data-driven method to predict the future evolution of the state. We show that our method is capable of predicting a critical transition event at least several numerical time steps in advance. We demonstrate the success as well as the limitations of our method using numerical experiments on three examples of systems, ranging from low dimensional to high dimensional. We discuss the mathematical and broader implications of our results.


title: Multistate ring polymer instantons and nonadiabatic reaction rates

abstract: We present two multistate ring polymer instanton (RPI) formulations, both obtained from an exact path integral representation of the quantum canonical partition function for multistate systems. The two RPIs differ in their treatment of the electronic degrees of freedom; whereas the Mean-Field (MF)-RPI averages over the electronic state contributions, the Mapping Variable (MV)-RPI employs explicit continuous Cartesian variables to represent the electronic states. We compute both RPIs for a series of model two-state systems coupled to a single nuclear mode with electronic coupling values chosen to describe dynamics in both adiabatic and nonadiabatic regimes. We show that the MF-RPI for symmetric systems are in good agreement with previous literature, and we show that our numerical techniques are robust for systems with non-zero driving force. The nuclear MF-RPI and the nuclear MV-RPI are similar, but the MV-RPI uniquely reports on the changes in the electronic state populations along the instanton path. In both cases, we analytically demonstrate the existence of a zero-mode and we numerically and that these solutions are true instantons with a single unstable mode as expected for a first order saddle point. Finally, we use the MF-RPI to accurately calculate rate constants for adiabatic and nonadiabatic model systems with the coupling strength varying over three orders of magnitude.


title: Approximation of the Lagrange and Markov spectra

abstract: The (classical) Lagrange spectrum is a closed subset of the positive real numbers defined in terms of diophantine approximation. Its structure is quite involved. This article describes a polynomial time algorithm to approximate it in Hausdorff distance. It also extends to approximate the Markov spectrum related to infimum of binary quadratic forms.


title: Uniform logical proofs for Riesz representation theorem, Daniell-Stone theorem and Stone's representation theorem for probability algebras

abstract: Riesz representation theorem, Daniell-Stone theorem for Daniell integrals and Stone's representation theorem for probability and measure algebras are three important classical results in analysis concerning existence of measures with certain properties. Many proofs of these theorems can be found in the literature of analysis, from elementary ones which use ordinary techniques from measure theory, to more sophisticated ones, such as those employing techniques from nonstandard analysis, in particular for Riesz representation theorem. In this paper, as the first goal, we give new proofs for all these three theorems. Our proofs have a mild logical flavor and are uniform in the sense that they are all based on the same general idea and rely on the application of the same technical tool from logic to measure theory, namely logical compactness theorem. In fact, as the second goal of the paper, we try to reveal more the power of logical methods in analysis in particular measure theory, and make stronger connections between analysis and logic. We use the setting of "integration logic" which is a logical framework (and one of the forms of probability logics) for studying measure and probability structures by logical means. Indeed, we elaborate this setting and use its expressive power and a version of compactness theorem holding in it to show its application in measure theory by giving new proofs for the above-mentioned measure existence theorems. As mentioned, an advantage of these proofs is that they are all given in a uniform way since they are all based on the logical compactness theorem. The paper is mostly written for general mathematicians, in particular the people active in analysis or logic as the main audience. So it is self-contained and the reader does not need to have any advanced prerequisite knowledge from logic or measure theory.


title: Lightweight and Scalable Particle Tracking and Motion Clustering of 3D Cell Trajectories

abstract: Tracking cell particles in 3D microscopy videos is a challenging task but is of great significance for modeling the motion of cells. Proper characterization of the cell's shape, evolution, and their movement over time is crucial to understanding and modeling the mechanobiology of cell migration in many diseases. One in particular, toxoplasmosis is the disease caused by the parasite Toxoplasma gondii. Roughly, one-third of the world's population tests positive for T. gondii. Its virulence is linked to its lytic cycle, predicated on its motility and ability to enter and exit nucleated cells; therefore, studies elucidating its motility patterns are critical to the eventual development of therapeutic strategies. Here, we present a computational framework for fast and scalable detection, tracking, and identification of T. gondii motion phenotypes in 3D videos, in a completely unsupervised fashion. Our pipeline consists of several different modules including preprocessing, sparsification, cell detection, cell tracking, trajectories extraction, parametrization of the trajectories; and finally, a clustering step. Additionally, we identified the computational bottlenecks, and developed a lightweight and highly scalable pipeline through a combination of task distribution and parallelism. Our results prove both the accuracy and performance of our method.


title: Lifting methods for manifold-valued variational problems

abstract: Lifting methods allow to transform hard variational problems such as segmentation and optical flow estimation into convex problems in a suitable higher-dimensional space. The lifted models can then be efficiently solved to a global optimum, which allows to find approximate global minimizers of the original problem. Recently, these techniques have also been applied to problems with values in a manifold. We provide a review of such methods in a refined framework based on a finite element discretization of the range, which extends the concept of sublabel-accurate lifting to manifolds. We also generalize existing methods for total variation regularization to support general convex regularization.


title: On the quenched functional CLT in 2d random sceneries, examples

abstract: We prove a quenched functional central limit theorem (quenched FCLT) for the sums of a random field (r.f.) along a 2d-random walk in different situations: when the r.f. is iid with a second order moment (random sceneries), or when it is generated by the action of commuting automorphisms of a torus. We consider also a quenched version of the FCLT when the random walk is replaced by a Lorentz process in the random scenery.


title: Spin-layer locking of interlayer excitons trapped in moir\'e potentials

abstract: Van der Waals heterostructures offer attractive opportunities to design quantum materials. For instance, transition metal dichalcogenides (TMDs) possess three quantum degrees of freedom: spin, valley index, and layer index. Further, twisted TMD heterobilayers can form moir\'e patterns that modulate the electronic band structure according to atomic registry, leading to spatial confinement of interlayer exciton (IXs). Here we report the observation of spin-layer locking of IXs trapped in moir\'e potentials formed in a heterostructure of bilayer 2H-MoSe$_2$ and monolayer WSe$_2$. The phenomenon of locked electron spin and layer index leads to two quantum-confined IX species with distinct spin-layer-valley configurations. Furthermore, we observe that the atomic registries of the moir\'e trapping sites in the three layers are intrinsically locked together due to the 2H-type stacking characteristic of bilayer TMDs. These results identify the layer index as a useful degree of freedom to engineer tunable few-level quantum systems in two-dimensional heterostructures.


title: The Polarized Three-Loop Anomalous Dimensions from On-Shell Massive Operator Matrix Elements

abstract: We calculate all contributions $\propto T_F$ to the polarized three-loop anomalous dimensions in the M-scheme using massive operator matrix elements and compare to results in the literature. This includes the complete anomalous dimensions $\gamma_{qq}^{(2),\rm PS}$ and $\gamma_{qg}^{(2)}$. We also obtain the complete two-loop polarized anomalous dimensions in an independent calculation. While for most of the anomalous dimensions the usual direct computation methods in Mellin $N$-space can be applied since all recurrences factorize at first order, this is not the case for $\gamma_{qg}^{(2)}$. Due to the necessity of deeper expansions of the master integrals in the dimensional parameter $\varepsilon = D-4$, we had to use the method of arbitrary high moments to eliminate elliptic contributions in intermediate steps. 4000 moments were generated to determine this anomalous dimension and 2640 moments turned out to be sufficient. As an aside, we also recalculate the contributions $\propto T_F$ to the three-loop QCD $\beta$-function.


title: Non-Hermitian coupled cluster method for non-stationary systems and its interaction-picture reinterpretation

abstract: The interaction picture in a non-Hermitian realization is discussed in detail and considered for its practical use in many-body quantum physics. The resulting non-Hermitian interaction-picture (NHIP) description of dynamics, in which both the wave functions and operators belonging to physical observables cease to remain constant in time, is a non-Hermitian generalization of the traditional Dirac picture of standard quantum mechanics, which itself is widely used in quantum field theory calculations. Particular attention is paid here to the variational (or, better, bivariational) and dynamical (i.e., non-stationary) aspects that are characteristic of the coupled cluster method (CCM) techniques that nowadays form one of the most versatile and most accurate of all available formulations of quantum many-body theory. In so doing we expose and exploit multiple parallels between the NHIP and the CCM in its time-dependent versions.


title: A theory of incremental compression

abstract: The ability to find short representations, i.e. to compress data, is crucial for many intelligent systems. We present a theory of incremental compression showing that arbitrary data strings, that can be described by a set of features, can be compressed by searching for those features incrementally, which results in a partition of the information content of the string into a complete set of pairwise independent pieces. The description length of this partition turns out to be close to optimal in terms of the Kolmogorov complexity of the string. Exploiting this decomposition, we introduce ALICE - a computable ALgorithm for Incremental ComprEssion - and derive an expression for its time complexity. Finally, we show that our concept of a feature is closely related to Martin-L\"of randomness tests, thereby formalizing the meaning of "property" for computable objects.


title: A Critical Note on the Evaluation of Clustering Algorithms

abstract: Experimental evaluation is a major research methodology for investigating clustering algorithms and many other machine learning algorithms. For this purpose, a number of benchmark datasets have been widely used in the literature and their quality plays a key role on the value of the research work. However, in most of the existing studies, little attention has been paid to the properties of the datasets and they are often regarded as black-box problems. For example, it is common to use datasets intended for classification in clustering research and assume class la-bels as the ground truth for judging the quality of cluster-ing. In our work, with the help of advanced visualization and dimension reduction techniques, we show that this practice may seriously compromise the research quality and produce misleading results. We suggest that the applicability of existing benchmark datasets should be carefully revisited and significant efforts need to be devoted to improving the current practice of experimental evaluation of clustering algorithms to ensure an essential match between algorithms and problems.


title: A note on degenerate Euler and Bernoulli polynomials of complex variable

abstract: In this paper, we study the degenerate version of the new type Euler polynomials, namely degenerate cosine-Euler polynomials and sime-Euler polynomials and also corresponding ones for Bernoulli polynomials, namely degenerate cosine Bernoulli polynomials and degenerate sine-Bernoulli polynomials by considering the degenerate Euler polynomials of complex variable and the degenerate Bernoulli polynomials of complex variable.


title: High-resistance YBa2Cu3O7-x grain-boundary Josephson junctions fabricated by electromigration

abstract: [100]-tilt grain-boundary YBa2Cu3O7-x (YBCO) junctions are promising for investigation of macroscopic quantum phenomena in high-Tc Josephson junctions. However, fabrication of the [100]-tilt grain-boundary YBCO junctions with a high resistance, which are required to study quantum effects, is difficult because of a high transparency of a tunnel barrier in this type of junctions. Here, we demonstrate a modification of grain-boundary barrier properties with a new approach to an oxygen electromigration in the YBCO grain-boundary junctions when the oxygen diffuses under an applied electric field from the grain-boundary to a BaTbO3 layer deposited atop of an YBCO film. Using this approach, we changed the normal-state resistance of the junctions from tens to several hundred Ohms without a degradation of their characteristic voltage IcRn and determined a barrier height and thickness by measuring the quasiparticle tunnelling current.


title: Pressure-Thresholded Response in Cylindrically Shocked Cyclotrimethylene Trinitramine (RDX)

abstract: We demonstrate a strongly thresholded response in cyclotrimethylene trinitramine (RDX) when it is cylindrically shocked using a novel waveguide geometry. Using ultrafast single-shot multi-frame imaging, we demonstrate that <100-{\mu}m diameter single crystals of RDX embedded in a polymer host deform along preferential planes for >100 ns after the shock first arrives in the crystal. We use in-situ imaging and time-resolved photoemission to demonstrate that short-lived chemistry occurs with complex deformation pathways. Using scanning electron microscopy and ultrasmall-angle X-ray scattering, we demonstrate that the shock-induced dynamics leave behind porous crystals, with pore shapes and sizes that change significantly with shock energy. A threshold pressure of ~ 12 GPa at the center of convergence separated the single-mode planar crystal deformations from the chemistry-coupled multi-plane dynamics at higher pressures. Our observations indicate preferential directions for deformation for our cylindrically shocked system, despite the applied stress along many different crystallographic planes.


title: An approximate factorization method for inverse acoustic scattering with phaseless total-field data

abstract: This paper is concerned with the inverse acoustic scattering problem with phaseless total-field data at a fixed frequency. An approximate factorization method is developed to numerically reconstruct both the location and shape of the unknown scatterer from the phaseless total-field data generated by incident plane waves at a fixed frequency and measured on the circle $\partial B_R$ with a sufficiently large radius $R$. The theoretical analysis of our method is based on the asymptotic property in the operator norm from $H^{1/2}(\mathbb{S}^1)$ to $H^{-1/2}(\mathbb{S}^1)$ of the phaseless total-field operator defined in terms of the phaseless total-field data measured on $\partial B_R$ with large enough $R$, where $H^s(\mathbb{S}^1)$ is a Sobolev space on the unit circle $\mathbb{S}^1$ for real number $s$, together with the factorization of a modified far-field operator. The asymptotic property of the phaseless total-field operator is also established in this paper with the theory of oscillatory integrals. The unknown scatterer can be either an impenetrable obstacle of sound-soft, sound-hard or impedance type or an inhomogeneous medium with a compact support, and the proposed inversion algorithm does not need to know the boundary condition of the unknown obstacle in advance. Numerical examples are also carried out to demonstrate the effectiveness of our inversion method. To the best of our knowledge, it is the first attempt to develop a factorization type method for inverse scattering problems with phaseless data.


title: Steady waves in flows over periodic bottoms

abstract: We study the formation of steady waves in two-dimensional fluids under a current with mean velocity $c$ flowing over a periodic bottom. Using a formulation based on the Dirichlet-Neumann operator, we establish the unique continuation of a steady solution from the trivial solution when a flat bottom is perturbed, except for a sequence of velocities $c_{k}$. The main contribution is the proof that at least two steady solutions exist close to a non-degenerate $S^{1}$-orbit of non-constant steady waves when a flat bottom is perturbed. Consequently, we obtain persistence of at least two steady waves close to a non-degenerate $S^{1}$-orbit of Stokes waves bifurcating from the velocities $c_{k}$.


title: Avoidable paths in graphs

abstract: We prove a recent conjecture of Beisegel et al. that for every positive integer k, every graph containing an induced P_k also contains an avoidable P_k. Avoidability generalises the notion of simpliciality best known in the context of chordal graphs. The conjecture was only established for k in {1,2} (Ohtsuki et al. 1976, and Beisegel et al. 2019, respectively). Our result also implies a result of Chv\'atal et al. 2002, which assumed cycle restrictions. We provide a constructive and elementary proof, relying on a single trick regarding the induction hypothesis. In the line of previous works, we discuss conditions for multiple avoidable paths to exist.


title: Coexistence of hidden attractors and self-excited attractors through breaking heteroclinic-like orbits of switched systems

abstract: In this paper an approach to generate hidden attractors based on piecewise linear (PWL) systems is studied. The approach consists of the coexistence of self-excited attrators and hidden attractors, i.e., the equilibria of the system are immersed in the basin of attraction of the seft-excited attractors, so hidden attractors appear around these self-excited attractors. The approach starts by generating a double-scroll attractor based on two equilibria presented heteroclinic orbits. Then, two equilibria are added to the system, so biestability is generated by displaying two self-excited attractor. In this paper we show that hidden attractors arise as a consequence of the rupture of trajectories that resemble heteroclinic orbits at a larger scale. Therefore, multistability appears naturally as an interesting phenomenon present in this class of dynamical systems via hidden attractors and self-excited attractors. The study suggests the feasibility of the geometric design of new classes of multistable systems with coexistence of the two classes of attractors or even multistable systems without equilibria at all.


title: Towards Online Observability-Aware Trajectory Optimization for Landmark-based Estimators

abstract: As autonomous systems increasingly rely on onboard sensing for localization and perception, the parallel tasks of motion planning and state estimation become more strongly coupled. This coupling is well-captured by augmenting the planning objective with a posterior-covariance penalty -- however, prediction of the estimator covariance is challenging when the observation model depends on unknown landmarks, as is the case in Simultaneous Localization and Mapping (SLAM). This paper addresses these challenges in the case of landmark- and SLAM-based estimators, enabling efficient prediction (and ultimately minimization) of this performance metric. First, we provide an interval-based filtering approximation of the SLAM inference process which allows for recursive propagation of the ego-covariance while avoiding the quadratic complexity of explicitly tracking landmark uncertainty. Secondly, we introduce a Lie-derivative measurement bundling scheme that simplifies the recursive "bundled" update, representing significant computational savings for high-rate sensors such as cameras. Finally, we identify a large class of measurement models (which includes orthographic camera projection) for which the contributions from each landmark can be directly combined, making evaluation of the information gained at each timestep (nearly) independent of the number of landmarks. This also enables the generalization from finite sets of landmarks $\{\ell^{(n)} \}$ to distributions, foregoing the need for fully-specified linearization points at planning time and allowing for new landmarks to be anticipated. Taken together, these contributions allow SLAM performance to be accurately and efficiently predicted, paving the way for online, observability-aware trajectory optimization in unknown space.


title: Low-Frequency Electronic Noise in Superlattice and Random-Packed Thin Films of Colloidal Quantum Dots

abstract: We report measurements of low-frequency electronic noise in ordered superlattice, weakly-ordered and random-packed thin films of 6.5 nm PbSe quantum dots prepared using several different ligand chemistries. For all samples, the normalized noise spectral density of the dark current revealed a Lorentzian component, reminiscent of the generation-recombination noise, superimposed on the 1/f background (f is the frequency). An activation energy of 0.3 eV was extracted from the temperature dependence of the noise spectra. The noise level in the ordered films was lower than that in the weakly-ordered and random-packed films. A large variation in the magnitude of the noise spectral density was also observed in samples with different ligand treatments. The obtained results are important for application of colloidal quantum dot films in photodetectors.


title: Object-Aware Instance Labeling for Weakly Supervised Object Detection

abstract: Weakly supervised object detection (WSOD), where a detector is trained with only image-level annotations, is attracting more and more attention. As a method to obtain a well-performing detector, the detector and the instance labels are updated iteratively. In this study, for more efficient iterative updating, we focus on the instance labeling problem, a problem of which label should be annotated to each region based on the last localization result. Instead of simply labeling the top-scoring region and its highly overlapping regions as positive and others as negative, we propose more effective instance labeling methods as follows. First, to solve the problem that regions covering only some parts of the object tend to be labeled as positive, we find regions covering the whole object focusing on the context classification loss. Second, considering the situation where the other objects contained in the image can be labeled as negative, we impose a spatial restriction on regions labeled as negative. Using these instance labeling methods, we train the detector on the PASCAL VOC 2007 and 2012 and obtain significantly improved results compared with other state-of-the-art approaches.


title: The rise and rise of interdisciplinary research: Understanding the interaction dynamics of three major fields -- Physics, Mathematics & Computer Science

abstract: The distinction between sciences is becoming increasingly more artificial -- an approach from one area can be easily applied to the other. More exciting research nowadays is happening perhaps at the interfaces of disciplines like Physics, Mathematics and Computer Science. How do these interfaces emerge and interact? For instance, is there a specific pattern in which these fields cite each other? In this article, we investigate a collection of more than 1.2 million papers from three different scientific disciplines -- Physics, Mathematics, and Computer Science. We show how over a timescale the citation patterns from the core science fields (Physics, Mathematics) to the applied and fast-growing field of Computer Science have drastically increased. Further, we observe how certain subfields in these disciplines are shrinking while others are becoming tremendously popular. For instance, an intriguing observation is that citations from Mathematics to the subfield of machine learning in Computer Science in recent times are exponentially increasing.


title: Quantitative Assessment of the Toner and Tu Theory of Polar Flocks

abstract: We present a quantitative assessment of the Toner and Tu theory describing the universal scaling of fluctuations in polar phases of dry active matter. Using large scale simulations of the Vicsek model in two and three dimensions, we find the overall phenomenology and generic algebraic scaling predicted by Toner and Tu, but our data on density correlations reveal some qualitative discrepancies. The values of the associated scaling exponents we estimate differ significantly from those conjectured in 1995. In particular, we identify a large crossover scale beyond which flocks are only weakly anisotropic. We discuss the meaning and consequences of these results.


title: Eigenvectors from eigenvalues: A survey of a basic identity in linear algebra

abstract: If $A$ is an $n \times n$ Hermitian matrix with eigenvalues $\lambda_1(A),\dots,\lambda_n(A)$ and $i,j = 1,\dots,n$, then the $j^{\mathrm{th}}$ component $v_{i,j}$ of a unit eigenvector $v_i$ associated to the eigenvalue $\lambda_i(A)$ is related to the eigenvalues $\lambda_1(M_j),\dots,\lambda_{n-1}(M_j)$ of the minor $M_j$ of $A$ formed by removing the $j^{\mathrm{th}}$ row and column by the formula $$ |v_{i,j}|^2\prod_{k=1;k\neq i}^{n}\left(\lambda_i(A)-\lambda_k(A)\right)=\prod_{k=1}^{n-1}\left(\lambda_i(A)-\lambda_k(M_j)\right)\,.$$ We refer to this identity as the \emph{eigenvector-eigenvalue identity} and show how this identity can also be used to extract the relative phases between the components of any given eigenvector. Despite the simple nature of this identity and the extremely mature state of development of linear algebra, this identity was not widely known until very recently. In this survey we describe the many times that this identity, or variants thereof, have been discovered and rediscovered in the literature (with the earliest precursor we know of appearing in 1834). We also provide a number of proofs and generalizations of the identity.


title: Single-Exciton Gain and Stimulated Emission Across the Infrared Optical Telecom Band from Robust Heavily-doped PbS Colloidal Quantum Dots

abstract: Materials with optical gain in the infrared are of paramount importance for optical communications, medical diagnostics1 and silicon photonics2,3 . The current technology is based either on costly III-V semiconductors that are not monolithic to silicon CMOS technology or Er-doped fiber technology that does not make use of the full fiber transparency window. Colloidal quantum dots (CQD) offer a unique opportunity as an optical gain medium4 in view of their tunable bandgap, solution processability and CMOS compatibility. Their potential for narrower linewidths5 and the lower-than-bulk degeneracy6 has led to dramatic progress towards successful demonstration of optical gain4, stimulated emission7 and lasing8,9,10 in the visible part of spectrum utilizing CdSe-based CQDs. Infrared Pb-chalcogenide colloidal quantum dots however exhibit higher state degeneracy and as a result the demonstration of optical gain has imposed very high thresholds.11,12 Here we demonstrate room-temperature, infrared stimulated emission, tunable across the optical communication band, based on robust electronically doped PbS CQDs, that reach gain threshold at the single exciton regime, representing a four-fold reduction from the theoretical limit of an eight-fold degenerate system and two orders of magnitude lower than prior reports.


title: A New Definition of the Dimension of Graphs

abstract: Enlighted by the recent work of Hao Huang on sensitivity conjecture [arXiv:1907.00847], we propose a new definition of the dimension of graphs and establish a relationship between the chromatic number and the dimension.


title: Newtonian orbits of nanoparticles interacting with structured light beams

abstract: We perform numerical analysis to study the orbits described by subwavelength size particles interacting with structured light beams. Our solution to the particle dynamics considers: (i) the gradient force, (ii) the radiation pressure, and (iii) the force from the curl of the spin angular momentum. The last two terms, (ii) and (iii), constitute the scattering forces. The optical structures of interest are vortices, vector, and Full-Poincar\'{e} beams. From our numerical results, we show that the particle is expelled from the beam, independent of the gradient force intensity if some of the scattering forces have cylindrical symmetry. Furthermore, we found spiral orbits for some particular conditions and types of Full-Poincare beams.


title: Radial Anharmonic Oscillator: Perturbation Theory, New Semiclassical Expansion, Approximating Eigenfunctions. I. Generalities, Cubic Anharmonicity Case

abstract: For the general $D$-dimensional radial anharmonic oscillator with potential $V(r)= \frac{1}{g^2}\,\hat{V}(gr)$ the Perturbation Theory (PT) in powers of coupling constant $g$ (weak coupling regime) and in inverse, fractional powers of $g$ (strong coupling regime) is developed constructively in $r$-space and in $(gr)$ space, respectively. The Riccati-Bloch (RB) equation and Generalized Bloch (GB) equation are introduced as ones which govern dynamics in coordinate $r$-space and in $(gr)$-space, respectively, exploring the logarithmic derivative of wavefunction $y$. It is shown that PT in powers of $g$ developed in RB equation leads to Taylor expansion of $y$ at small $r$ while being developed in GB equation leads to a new form of semiclassical expansion at large $(g r)$: it coincides with loop expansion in path integral formalism. In complementary way PT for large $g$ developed in RB equation leads to an expansion of $y$ at large $r$ and developed in GB equation leads to an expansion at small $(g r)$. Interpolating all four expansions for $y$ leads to a compact function (called the {\it Approximant}), which should uniformly approximate the exact eigenfunction at $r \in [0, \infty)$ for any coupling constant $g \geq 0$ and dimension $D > 0$. Free parameters of the Approximant are fixed by taking it as a trial function in variational calculus. As a concrete application the low-lying states of the cubic anharmonic oscillator $V=r^2+gr^3$ are considered. It is shown that the relative deviation of the Approximant from the exact ground state eigenfunction is $\lesssim 10^{-4}$ for $r \in [0, \infty)$ for coupling constant $g \geq 0$ and dimension $D=1,2,\ldots$. In turn, the variational energies of the low-lying states are obtained with unprecedented accuracy 7-8 s.d. for $g \geq 0$ and $D=1,2,\ldots$.


title: Large Scale Geometries of Infinite Strings

abstract: We introduce geometric consideration into the theory of formal languages. We aim to shed light on our understanding of global patterns that occur on infinite strings. We utilise methods of geometric group theory. Our emphasis is on large scale geometries. Two infinite strings have the same large scale geometry if there are colour preserving bi-Lipschitz maps with distortions between the strings. Call these maps quasi-isometries. Introduction of large scale geometries poses several questions. The first question asks to study the partial order induced by quasi-isometries. This partial order compares large scale geometries; as such it presents an algebraic tool for classification of global patterns. We prove there is a greatest large scale geometry and infinitely many minimal large scale geometries. The second question is related to understanding the quasi-isometric maps on various classes of strings. The third question investigates the sets of large scale geometries of strings accepted by computational models, e.g. B\"uchi automata. We provide an algorithm that describes large scale geometries of strings accepted by B\"uchi automata. This links large scale geometries with automata theory. The fourth question studies the complexity of the quasi-isometry problem. We show the problem is $\Sigma_3^0$-complete thus providing a bridge with computability theory. Finally, the fifth question asks to build algebraic structures that are invariants of large scale geometries. We invoke asymptotic cones, a key concept in geometric group theory, defined via model-theoretic notion of ultra-product. Partly, we study asymptotic cones of algorithmically random strings thus connecting the topic with algorithmic randomness.


title: Some Orbits of Free Words that are Determined by Measures on Finite Groups

abstract: Every word in a free group $F$ induces a probability measure on every finite group in a natural manner. It is an open problem whether two words that induce the same measure on every finite group, necessarily belong to the same orbit of $\mathrm{Aut}F$. A special case of this problem, when one of the words is the primitive word $x$, was settled positively by the third author and Parzanchevski [arXiv:1202.3269]. Here we extend this result to the case where one of the words is $x^d$ or $\left[x,y\right]^{d}$ for an arbitrary $d\in\mathbb{Z}$.


title: Almost-rigidity of frameworks

abstract: We extend the mathematical theory of rigidity of frameworks (graphs embedded in $d$-dimensional space) to consider nonlocal rigidity and flexibility properties. We provide conditions on a framework under which (I) as the framework flexes continuously it must remain inside a small ball, a property we call "almost-rigidity"; (II) any other framework with the same edge lengths must lie outside a much larger ball; (III) if the framework deforms by some given amount, its edge lengths change by a minimum amount; (IV) there is a nearby framework that is prestress stable, and thus rigid. The conditions can be tested efficiently using semidefinite programming. The test is a slight extension of the test for prestress stability of a framework, and gives analytic expressions for the radii of the balls and the edge length changes. Examples illustrate how the theory may be applied in practice, and we provide an algorithm to test for rigidity or almost-rigidity. We briefly discuss how the theory may be applied to tensegrities.


title: Cloud-based Management of Energy-Efficient Dense IEEE 802.11ax Networks

abstract: During the last decade, the number of devices connected to the Internet by Wi-Fi has grown significantly. A high density of both the client devices and the hot spots posed new challenges related to providing the desired quality of service in the current and emerging scenarios. To cope with the negative effects caused by network densification, modern Wi-Fi is becoming more and more centralized. To improve network efficiency, today many new Wi-Fi deployments are under control of management systems that optimize network parameters in a centralized manner. In the paper, for such a cloud management system, we develop an algorithm which aims at maximizing energy efficiency and also keeps fairness among clients. For that, we design an objective function and solve an optimization problem using the branch and bound approach. To evaluate the efficiency of the developed solution, we implement it in the NS-3 simulator and compare with existing solutions and legacy behavior.


title: New Constructions of Subspace Codes Using Subsets of MRD codes in Several Blocks

abstract: A basic problem for the constant dimension subspace coding is to determine the maximal possible size A_q (n, d, k) of a set of k-dimensional subspaces in Fnq such that the subspace distance satisfies d(U, V )> or =d for any two different subspaces U andV in this set. We present two new constructions of constant dimension subspace codes using subsets of maximal rank-distance (MRD) codes in several blocks. This method is firstly applied to the linkage construction and secondly to arbitrary number of blocks of lifting MRD codes. In these two constructions, subsets of MRD codes with bounded ranks play an essential role. The Delsarte theorem of the rank distribution of MRD codes is an important ingredient to count codewords in our constructed constant dimension subspace codes. We give many new lower bounds for A_q (n, d, k). More than 110 new constant dimension subspace codes better than previously best known codes are constructed.


title: Anderson localization for multi-frequency quasi-periodic operators on $\mathbb{Z}^d$

abstract: We establish Anderson localization for general analytic $k$-frequency quasi-periodic operators on $\mathbb{Z}^d$ for \textit{arbitrary} $k,d$.


title: Theory of Optical Nonlocality in Polar Dielectrics

abstract: Sub-wavelength confinement of mid-infrared light can be achieved exploiting the metal-like optical response of polar dielectric crystals in their Reststrahlen spectral region, where they support evanescent modes termed surface phonon polaritons. In the past few years the investigation of phonon polaritons localised in nanoresonators and layered heterostructures has enjoyed remarkable success, highlighting them as a promising platform for mid-infrared nanophotonic applications. Here we prove that the standard local dielectric description of phonon polaritons in nanometric objects fails due to the nonlocal nature of the phonon response and we develop the corresponding nonlocal theory. Application of our general theory to both dielectric nanospheres and thin films demonstrates that polar dielectrics exhibit a rich nonlocal phenomenology, qualitatively different from the one of plasmonic systems, due to the negative dispersion of phononic optical modes.


title: Above-barrier heavy-ion fusion cross-sections using the relativistic mean-field approach: case of spherical colliding nuclei

abstract: The double folding (DF) approach is one of the widely used methods for finding nucleus-nucleus interaction potential. In the present work, the influence of the nuclear matter density on the DF potential and on the Coulomb barrier parameters is studied systematically for collisions of spherical nuclei. The value of the parameter $B_Z=Z_P Z_T/(A_P^{1/3}+A_T^{1/3})$ (estimating the Coulomb barrier height) varies in these calculations from 10 MeV up to 150 MeV. The novel feature of this study is that the nuclear densities came from the Relativistic Mean Field approach (RMF). For the nucleon-nucleon effective interaction, the M3Y forces with the finite range exchange term and density dependence are employed. The above barrier fusion cross sections are calculated within the framework of the trajectory model with surface friction. Results are compared with the previous study in which the nuclear density came from the Hartree-Fock (HF) calculations and with the high precision experimental data. This comparison demonstrates that i) agreement between the theoretical and experimental cross sections obtained with RMF and HF densities is of the same quality and ii) the values of the only adjustable parameter (friction strength) obtained with RMF and HF densities strongly correlate.


title: Noncompact complete Riemannian manifolds with singular continuous spectrum embedded into the essential spectrum of the Laplacian, I. The hyperbolic case

abstract: We construct Riemannian manifolds with singular continuous spectrum embedded in the absolutely continuous spectrum of the Laplacian. Our manifolds are asymptotically hyperbolic with sharp curvature bounds.