title: Compressive Embedding and Visualization using Graphs

abstract: Visualizing high-dimensional data has been a focus in data analysis communities for decades, which has led to the design of many algorithms, some of which are now considered references (such as t-SNE for example). In our era of overwhelming data volumes, the scalability of such methods have become more and more important. In this work, we present a method which allows to apply any visualization or embedding algorithm on very large datasets by considering only a fraction of the data as input and then extending the information to all data points using a graph encoding its global similarity. We show that in most cases, using only $\mathcal{O}(\log(N))$ samples is sufficient to diffuse the information to all $N$ data points. In addition, we propose quantitative methods to measure the quality of embeddings and demonstrate the validity of our technique on both synthetic and real-world datasets.


title: The isoperimetric problem for Lens spaces

abstract: We solve the isoperimetric problem in the Lens spaces with large fundamental group. Namely, we prove that the isoperimetric surfaces are geodesic spheres or tori of revolution about geodesics. We also show that the isoperimetric problem in L(3,1) and L(3,2) follows from the proof of the Willmore conjecture by Marques and Neves.


title: Formation of Laves Phases in Repulsive and Attractive Hard Sphere Suspensions

abstract: Colloidal Laves phases (LPs) of $MgCu_2$ type are promising precursors for diamond structure photonic materials. They have been predicted for hard sphere binary mixtures, but not yet observed. We here report a time resolved static light scattering study on their formation in a binary mixture of buoyant experimental hard sphere approximants (size ratio $\Gamma=0.77$, molar fraction of small spheres $x_S = 0.76$) for volume fractions between melting and the glass transition. In line with theoretical expectation, all samples form LPs of $MgZn_2$ structure on the time scale of weeks to months. $MgNi_2$ structures are absent, $MgCu_2$ structures and randomly stacked LPs prevail at elevated volume fraction. The addition of small amounts of non-adsorbing polymer switches the interaction to depletion attractive and results in significantly accelerated crystallization kinetics and improved crystal quality.


title: Two definitions of the Gell-Mann channels - a comparative analysis

abstract: Our goal is to compare two unequivalent definitions of the Gell-Mann channels. It turns out that both definitions coincide for qubits and qutrits. In higher dimensions, there exist some constraints under which the channels describe the same dynamics. Finally, we find the GKSL time-local generators for a class of the Gell-Mann channels.


title: The 2-adic valuation of generalized Fibonacci sequences with an application to certain Diophantine equations

abstract: In this paper we focus on finding all the factorials expressible as a product of a fixed number of $2k$-nacci numbers with $k \geq 2$. We derive the 2-adic valuation of the $2k$-nacci sequence and use it to establish bounds on the solutions of the initial equation. In addition, we specify a more general family of sequences, for which we can perform a similar procedure. We also investigate a possible connection of these results with $p$-regular sequences.


title: Spin relaxation in corrugated graphene

abstract: In graphene, out-of-plane (flexural) vibrations and static ripples imposed by the substrate relax the electron spin, intrinsically protected by mirror symmetry. We calculate the relaxation times in different scenarios, accounting for all the possible spin-phonon couplings allowed by the hexagonal symmetry of the lattice. Scattering by flexural phonons imposes the ultimate bound to the spin lifetimes, in the ballpark of hundreds of nano-seconds at room temperature. This estimate and the behavior as a function of the carrier concentration are substantially altered by the presence of tensions or the pinning with the substrate. Static ripples also influence the spin transport in the diffusive regime, dominated by motional narrowing. We find that the D'yakonov-Perel' mechanism saturates when the mean free path is comparable to the correlation length of the heights profile. In this regime, the spin-relaxation times are exclusively determined by the geometry of the corrugations. Simple models for typical corrugations lead to lifetimes of the order of tens of micro-seconds.


title: Post-edit Analysis of Collective Biography Generation

abstract: Text generation is increasingly common but often requires manual post-editing where high precision is critical to end users. However, manual editing is expensive so we want to ensure this effort is focused on high-value tasks. And we want to maintain stylistic consistency, a particular challenge in crowd settings. We present a case study, analysing human post-editing in the context of a template-based biography generation system. An edit flow visualisation combined with manual characterisation of edits helps identify and prioritise work for improving end-to-end efficiency and accuracy.


title: Environmentally mediated coherent control of a spin qubit in diamond

abstract: The coherent control of spin qubits forms the basis of many applications in quantum information processing and nanoscale sensing, imaging and spectroscopy. Such control is conventionally achieved by direct driving of the qubit transition with a resonant global field, typically at microwave frequencies. Here we introduce an approach that relies on the resonant driving of nearby environment spins, whose localised magnetic field in turn drives the qubit when the environmental spin Rabi frequency matches the qubit resonance. This concept of environmentally mediated resonance (EMR) is explored experimentally using a qubit based on a single nitrogen-vacancy (NV) centre in diamond, with nearby electronic spins serving as the environmental mediators. We demonstrate EMR driven coherent control of the NV spin-state, including the observation of Rabi oscillations, free induction decay, and spin-echo. This technique also provides a way to probe the nanoscale environment of spin qubits, which we illustrate by acquisition of electron spin resonance spectra of single NV centres in various settings.


title: Improved lower bound for the number of unimodular zeros of self-reciprocal polynomials with coefficients in a finite set

abstract: Let $n_1 < n_2 < \cdots < n_N$ be non-negative integers. In a private communication Brian Conrey asked how fast the number of real zeros of the trigonometric polynomials $T_N(\theta) = \sum_{j=1}^N {\cos (n_j\theta)}$ tends to $\infty$ as a function of $N$. Conrey's question in general does not appear to be easy.Let ${\mathcal P}_n(S)$ be the set of all algebraic polynomials of degree at most $n$ with each of their coefficients in $S$. For a finite set $S \subset {\mathbb C}$ let $M = M(S) := \max\{|z|: z \in S\}$. It has been shown recently that if $S \subset {\mathbb R}$ is a finite set and $(P_n)$ is a sequence of self-reciprocal polynomials $P_n \in {\mathcal P}_n(S)$ with $|P_n(1)|$ tending to $\infty$, then the number of zeros of $P_n$ on the unit circle also tends to $\infty$. In this paper we show that if $S \subset {\mathbb Z}$ is a finite set, then every self-reciprocal polynomial $P \in {\mathcal P}_n(S)$ has at least $$c(\log\log\log|P(1)|)^{1-\varepsilon}-1$$ zeros on the unit circle of ${\mathbb C}$ with a constant $c > 0$ depending only on $\varepsilon > 0$ and$M = M(S)$. Our new result improves the exponent $1/2 - \varepsilon$ in a recent result by Julian Sahasrabudhe to $1 - \varepsilon$. Sahasrabudhe's new idea [66] is combined with the approach used in [34] offering an essencially simplified way to achieve our improvement. We note that in both Sahasrabudhe's paper and our paper the assumption that the finite set $S$ contains only integers is deeply exploited.


title: Large deviations of a tracer in the symmetric exclusion process

abstract: The one-dimensional symmetric exclusion process, the simplest interacting particle process, is a lattice-gas made of particles that hop symmetrically on a discrete line respecting hard-core exclusion. The system is prepared on the infinite lattice with a step initial profile with average densities $\rho_{+}$ and $\rho_{-}$ on the right and on the left of the origin. When $\rho_{+} = \rho_{-}$, the gas is at equilibrium and undergoes stationary fluctuations. When these densities are unequal, the gas is out of equilibrium and will remain so forever. A tracer, or a tagged particle, is initially located at the boundary between the two domains; its position $X_t$ is a random observable in time, that carries information on the non-equilibrium dynamics of the whole system. We derive an exact formula for the cumulant generating function and the large deviation function of $X_t$, in the long time limit, and deduce the full statistical properties of the tracer's position. The equilibrium fluctuations of the tracer's position, when the density is uniform, are obtained as an important special case.


title: Sustainable Fair Division

abstract: In this paper, I summarize our work on online fair division. In particular, I present two models for online fair division: (1) one existing model for fair division in food banks and (2) one new model for fair division of deceased organs to patients. I further discuss simple mechanisms for these models that allocate the resources as they arrive to agents. In practice, agents are often risk-averse having imperfect information. Within this assumption, I report several interesting axiomatic and complexity results for these mechanisms and conclude with future work.


title: Network Systems and String Stability

abstract: Network systems and their control are highly important and appear in a variety of applications, including vehicle platooning and formation con- trol. Especially vehicle platoons are highly investigated and an interesting problem that arises in this area is string stability, which broadly spoken means that a input signal amplifies unbounded as it travels through the vehicle string. However, various definitions are commonly used. In this paper, we aim to formalise the notion of string stability and illustrate the importance of those distinctions on simulation examples. A second goal is to generalise the found definitions for general network systems.


title: Improved lower bounds for the Mahler measure of the Fekete polynomials

abstract: We show that there is an absolute constant $c > 1/2$ such that the Mahler measure of the Fekete polynomials $f_p$ of the form $$f_p(z) := \sum_{k=1}^{p-1}{\left( \frac kp \right)z^k}\,,$$ (where the coefficients are the usual Legendre symbols) is at least $c\sqrt{p}$ for all sufficiently large primes $p$. This improves the lower bound $\left(\frac 12 - \varepsilon\right)\sqrt{p}$ known before for the Mahler measure of the Fekete polynomials $f_p$ for all sufficiently large primes $p \geq c_{\varepsilon}$. Our approach is based on the study of the zeros of the Fekete polynomials on the unit circle.


title: Failure and Scaling of Graphene Nanocomposites

abstract: This work proposes an investigation on the scaling of the structural strength of polymer/graphene nanocomposites. To this end, fracture tests on geometrically scaled Single Edge Notch Bending (SENB) specimens with varying contents of graphene were conducted to study the effects of nanomodification on the scaling. It is shown that, while the strength of the pristine polymer scales according to Linear Elastic Fracture Mechanics (LEFM), this is not the case for nanocomposites, even for very low graphene contents. In fact, small specimens exhibited a more pronounced ductility with limited scaling and a significant deviation from LEFM whereas larger specimens behaved in a more brittle way, with scaling of nominal strength closer to the one predicted by LEFM. This behavior, due to the significant size of the Fracture Process Zone (FPZ) compared to the specimen size, needs to be taken into serious consideration. In facts, it is shown that, for the specimen sizes investigated in this work, neglecting the non-linear effects of the FPZ can lead to an underestimation of the fracture energy as high as 113%, this error decreasing for increasing specimen sizes.


title: Copula-based piecewise regression

abstract: Most common parametric families of copulas are totally ordered, and in many cases they are also positively or negatively regression dependent and therefore they lead to monotone regression functions, which makes them not suitable for dependence relationships that imply or suggest a non-monotone regression function. A gluing copula approach is proposed to decompose the underlying copula into totally ordered copulas that combined may lead to a non-monotone regression function.


title: A one-dimensional mathematical model of collecting lymphatics coupled with an electro-fluid-mechanical contraction model and valve dynamics

abstract: We propose a one-dimensional model for collecting lymphatics coupled with a novel Electro-Fluid-Mechanical Contraction (EFMC) model for dynamical contractions, based on a modified FitzHugh-Nagumo model for action potentials. The one-dimensional model for a compliant lymphatic vessel is a set of hyperbolic Partial Differential Equations (PDEs). The EFMC model combines the electrical activity of lymphangions (action potentials) with fluid-mechanical feedback (stretch of the lymphatic wall and wall shear stress) and the mechanical variation of the lymphatic wall properties (contractions). The EFMC model is governed by four Ordinary Differential Equations (ODEs) and phenomenologically relies on: (1) environmental calcium influx, (2) stretch-activated calcium influx, and (3) contraction inhibitions induced by wall shear stresses. We carried out a complete mathematical analysis of the stability of the stationary state of the EFMC model. Overall, the EFMC model allows imitating the influence of pressure and wall shear stress on the frequency of contractions observed experimentally. Lymphatic valves are modelled using a well-established lumped-parameter model which allows simulating stenotic and regurgitant valves. We analysed several lymphodynamical indexes of a single lymphangion for a wide range of upstream and downstream pressure combinations. Stenotic and regurgitant valves were modelled, and their effects are here quantified. Results for stenotic valves showed in the downstream lymphangion that for low frequencies of contractions the Calculated Pump Flow (CPF) index remained almost unaltered, while for high frequencies the CPF dramatically decreased depending on the severity of the stenosis (up to 93% for a severe stenosis). Results for incompetent valves showed that the net flow during a lymphatic cycle tends to zero as the degree of incompetence increases.


title: Radio Emission from Pulsar Wind Nebulae without Surrounding Supernova Ejecta: Application to FRB 121102

abstract: In this paper, we propose a new scenario in which a rapidly-rotating strongly-magnetized pulsar without any surrounding supernova ejecta produces fast radio bursts (FRBs) repeatedly via some mechanisms, and meanwhile, an ultra-relativistic electron/positron pair wind from the pulsar sweeps up its ambient dense interstellar medium, giving rise to a non-relativistic pulsar wind nebula (PWN). We show that the synchrotron radio emission from such a PWN is bright enough to account for the recently-discovered persistent radio source associated with the repeating FRB 121102 in reasonable ranges of the model parameters. In addition, our PWN scenario is consistent with the non-evolution of the dispersion measure inferred from all the repeating bursts observed in four years.


title: Robust Hierarchical Bayes Small Area Estimation for Nested Error Regression Model

abstract: National statistical institutes in many countries are now mandated to produce reliable statistics for important variables such as population, income, unemployment, health outcomes, etc. for small areas, defined by geography and/or demography. Due to small samples from these areas, direct sample-based estimates are often unreliable. Model-based small area estimation is now extensively used to generate reliable statistics by "borrowing strength" from other areas and related variables through suitable models. Outliers adversely influence standard model-based small area estimates. To deal with outliers, Sinha and Rao (2009) proposed a robust frequentist approach. In this article, we present a robust Bayesian alternative to the nested error regression model for unit-level data to mitigate outliers. We consider a two-component scale mixture of normal distributions for the unit-level error to model outliers and present a computational approach to produce Bayesian predictors of small area means under a noninformative prior for model parameters. A real example and extensive simulations convincingly show robustness of our Bayesian predictors to outliers. Simulations comparison of these two procedures with Bayesian predictors by Datta and Ghosh (1991) and M-quantile estimators by Chambers et al. (2014) shows that our proposed procedure is better than the others in terms of bias, variability, and coverage probability of prediction intervals, when there are outliers. The superior frequentist performance of our procedure shows its dual (Bayes and frequentist) dominance, and makes it attractive to all practitioners, both Bayesian and frequentist, of small area estimation.


title: IEEE 802.11ad-based Radar: An Approach to Joint Vehicular Communication-Radar System

abstract: Millimeter-wave (mmWave) radar is widely used in vehicles for applications such as adaptive cruise control and collision avoidance. In this paper, we propose an IEEE 802.11ad-based radar for long-range radar (LRR) applications at the 60 GHz unlicensed band. We exploit the preamble of a single-carrier (SC) physical layer (PHY) frame, which consists of Golay complementary sequences with good correlation properties, as a radar waveform. This system enables a joint waveform for automotive radar and a potential mmWave vehicular communication system based on IEEE 802.11ad, allowing hardware reuse. To formulate an integrated framework of vehicle-to-vehicle (V2V) communication and LRR based on a mmWave consumer wireless local area network (WLAN) standard, we make typical assumptions for LRR applications and incorporate the full duplex radar assumption due to the possibility of sufficient isolation and self-interference cancellation. We develop single- and multi-frame radar receiver algorithms for target detection as well as range and velocity estimation within a coherent processing interval. Our proposed radar processing algorithms leverage channel estimation and time-frequency synchronization techniques used in a conventional IEEE 802.11ad receiver with minimal modifications. Analysis and simulations show that in a single target scenario, a Gbps data rate is achieved simultaneously with cm-level range accuracy and cm/s-level velocity accuracy. The target vehicle is detected with a high probability of detection ($>$99.9$\%$) at a low false alarm of 10$^{-6}$ for an equivalent isotropically radiated power (EIRP) of 43 dBm up to a vehicle separation distance of 200 m.


title: Symmetric structure for the endomorphism algebra of projective-injective module in parabolic category

abstract: We show that for any singular dominant integral weight $\lambda$ of a complex semisimple Lie algebra $\mathfrak{g}$, the endomorphism algebra $B$ of any projective-injective module of the parabolic BGG category $\mathcal{O}_\lambda^{\mathfrak{p}}$ is a symmetric algebra (as conjectured by Khovanov) extending the results of Mazorchuk and Stroppel for the regular dominant integral weight. Moreover, the endomorphism algebra $B$ is equipped with a homogeneous (non-degenerate) symmetrizing form. In the appendix, there is a short proof due to K. Coulembier and V. Mazorchuk showing that the endomorphism algebra $B_\lambda^{\mathfrak{p}}$ of the basic projective-injective module of $\mathcal{O}_\lambda^{\mathfrak{p}}$ is a symmetric algebra.


title: Spatial-Temporal Imaging of Anisotropic Photocarrier Dynamics in Black Phosphorus

abstract: As an emerging single elemental layered material with a low symmetry in-plane crystal lattice, black phosphorus (BP) has attracted significant research interest owing to its unique electronic and optoelectronic properties, including its widely tunable bandgap, polarization dependent photoresponse and highly anisotropic in-plane charge transport. Despite extensive study of the steady-state charge transport in BP, there has not been direct characterization and visualization of the hot carriers dynamics in BP immediately after photoexcitation, which is crucial to understanding the performance of BP-based optoelectronic devices. Here we use the newly developed scanning ultrafast electron microscopy (SUEM) to directly visualize the motion of photo-excited hot carriers on the surface of BP in both space and time. We observe highly anisotropic in-plane diffusion of hot holes, with a 15-times higher diffusivity along the armchair (x-) direction than that along the zigzag (y-) direction. Our results provide direct evidence of anisotropic hot carrier transport in BP and demonstrate the capability of SUEM to resolve ultrafast hot carrier dynamics in layered two-dimensional materials.


title: The socle filtrations of principal series representations of $SL(3,\mathbb{R})$ and $Sp(2,\mathbb{R})$

abstract: We study the structure of the $(\mathfrak{g},K)$-modules of the principal series representations of $SL(3,\mathbb{R})$ and $Sp(2,\mathbb{R})$ induced from minimal parabolic subgroups, in the case when the infinitesimal character is nonsingular. The composition factors of these modules are known by Kazhdan-Lusztig-Vogan conjecture. In this paper, we give complete descriptions of the socle filtrations of these modules.


title: Search for the Radiative Leptonic Decay $D^{+}\to \gamma e^{+} {\nu}_{e}$

abstract: Using an electron-positron collision data sample of 2.93 fb$^{-1}$ collected at a center-of-mass energy of $\sqrt{s}=3.773$ GeV with the BESIII detector, we present the first search for the radiative leptonic decay $D^{+} \rightarrow \gamma e^{+}{\nu}_{e}$. The analysis is performed with a double tag method. We do not observe a significant $D^{+} \rightarrow \gamma e^{+}{\nu}_{e}$ signal, and obtain an upper limit on the branching fraction of $D^{+} \rightarrow \gamma e^{+}{\nu}_{e}$ decay with the energy of radiative photon larger than 10 MeV of $3.0\times10^{-5}$ at the 90% confidence level.


title: Temporal Observables and Entangled Histories

abstract: We demonstrate that temporal observables, which are sensitive to a system's history (as opposed to its state), implicate entangled histories. We exemplify protocols for measuring such observables, and algorithms for predicting the (stochastic) outcomes of such measurements. Temporal observables allow us to define, and potentially to measure, precise mathematical consequences of intrinsically disjoint, yet mutually accessible, branches within the evolution of a pure quantum state.


title: Progressively Diffused Networks for Semantic Image Segmentation

abstract: This paper introduces Progressively Diffused Networks (PDNs) for unifying multi-scale context modeling with deep feature learning, by taking semantic image segmentation as an exemplar application. Prior neural networks, such as ResNet, tend to enhance representational power by increasing the depth of architectures and driving the training objective across layers. However, we argue that spatial dependencies in different layers, which generally represent the rich contexts among data elements, are also critical to building deep and discriminative representations. To this end, our PDNs enables to progressively broadcast information over the learned feature maps by inserting a stack of information diffusion layers, each of which exploits multi-dimensional convolutional LSTMs (Long-Short-Term Memory Structures). In each LSTM unit, a special type of atrous filters are designed to capture the short range and long range dependencies from various neighbors to a certain site of the feature map and pass the accumulated information to the next layer. From the extensive experiments on semantic image segmentation benchmarks (e.g., ImageNet Parsing, PASCAL VOC2012 and PASCAL-Part), our framework demonstrates the effectiveness to substantially improve the performances over the popular existing neural network models, and achieves state-of-the-art on ImageNet Parsing for large scale semantic segmentation.


title: Understanding the Cancelation of Double Poles in the Pfaffian of CHY-formulism

abstract: For a physical field theory, the tree-level amplitudes should possess only single poles. However, when computing amplitudes with Cachazo-He-Yuan (CHY) formulation, individual terms in the intermediate steps will contribute higher-order poles. In this paper, we investigate the cancelation of higher-order poles in CHY formula with Pfaffian as the building block. We develop a diagrammatic rule for expanding the reduced Pfaffian. Then by organizing diagrams in appropriate groups and applying the cross-ratio identities, we show that all potential contributions to higher-order poles in the reduced Pfaffian are canceled out, i.e., only single poles survive in Yang-Mills theory and gravity. Furthermore, we show the cancelations of higher-order poles in other field theories by introducing appropriate truncations, based on the single pole structure of Pfaffian.


title: Continuation Methods for Computing Z-/H-eigenpairs of Nonnegative Tensors

abstract: In this paper, a homotopy continuation method for the computation of nonnegative Z-/H-eigenpairs of a nonnegative tensor is presented. We show that the homotopy continuation method is guaranteed to compute a nonnegative eigenpair. Additionally, using degree analysis, we show that the number of positive Z-eigenpairs of an irreducible nonnegative tensor is odd. A novel homotopy continuation method is proposed to compute an odd number of positive Z-eigenpairs, and some numerical results are presented.


title: Scale-invariance in soft gamma repeaters

abstract: The statistical properties of the soft gamma repeater SGR J1550--5418 are investigated carefully. We find that the cumulative distributions of fluence, peak flux and duration can be well fitted by a bent power law, while the cumulative distribution of waiting time follows a simple power law. In particular, the probability density functions of fluctuations of fluence, peak flux, and duration have a sharp peak and fat tails, which can be well fitted by a $q$-Gaussian function. The $q$ values keep approximately steady for different scale intervals, indicating a scale-invariant structure of soft gamma repeaters. Those results support that the origin of soft gamma repeaters is crustquakes of neutron stars with extremely strong magnetic fields.


title: Chaos Engineering

abstract: Modern software-based services are implemented as distributed systems with complex behavior and failure modes. Many large tech organizations are using experimentation to verify the reliability of such systems. We use the term "Chaos Engineering" to refer to this approach, and discuss the underlying principles and how to use it to run experiments.


title: Hyperfine Wave Functions and Force Densities for the Hydrogen Atom

abstract: This study addresses the effect of the magnetic hyperfine interaction on the relativistic H1s wave functions. These are used to calculate the electric, magnetic, and confinement force densities acting on the 1s electron. The magnetic field couples Dirac equations for different angular momenta. These are solved numerically for the hyperfine singlet and triplet, as well as for a classical magnetic dipole. In the singlet ground state the hyperfine interaction shifts the electron density toward the proton. A similar shift is found for the classical dipole, and an opposite shift for the triplet. The cross-over between charge accumulation and depletion occurs at 1.325 times the Bohr radius. The behavior of the wave functions is investigated down to distances smaller than the proton radius, including the incorporation of virtual positrons. The force densities are determined and balanced against each other.


title: Intertwining operators among twisted modules associated to not-necessarily-commuting automorphisms

abstract: We introduce intertwining operators among twisted modules or twisted intertwining operators associated to not-necessarily-commuting automorphisms of a vertex operator algebra. Let $V$ be a vertex operator algebra and let $g_{1}$, $g_{2}$ and $g_{3}$ be automorphisms of $V$. We prove that for $g_{1}$-, $g_{2}$- and $g_{3}$-twisted $V$-modules $W_{1}$, $W_{2}$ and $W_{3}$, respectively, such that the vertex operator map for $W_{3}$ is injective, if there exists a twisted intertwining operator of type ${W_{3}\choose W_{1}W_{2}}$ such that the images of its component operators span $W_{3}$, then $g_{3}=g_{1}g_{2}$. We also construct what we call the skew-symmetry and contragredient isomorphisms between spaces of twisted intertwining operators among twisted modules of suitable types. The proofs of these results involve careful analysis of the analytic extensions corresponding to the actions of the not-necessarily-commuting automorphisms of the vertex operator algebra.


title: Novel Techniques to Derive Capacity Results for Multi-User Interference Channels

abstract: The Interference Channels (ICs) represent fundamental building blocks of wireless communication networks. Despite considerable progress in network information theory, available capacity results for ICs, specifically those with more than two users, are still very limited. One of the main difficulties in the analysis of these networks is how to establish useful capacity outer bounds for them. In this paper, novel techniques requiring subtle sequential application of Csiszar-Korner identity are developed to establish efficient single-letter outer bounds on the sum-rate capacity of interference networks. By using the outer bounds, a full characterization of the sum-rate capacity is then derived for various multi-user ICs under specific conditions. Our capacity results hold for both discrete and Gaussian networks.


title: Constant-roll inflation: confrontation with recent observational data

abstract: The previously proposed class of phenomenological inflationary models in which the assumption of inflaton slow-roll is replaced by the more general, constant-roll condition is compared with the most recent cosmological observational data, mainly the Planck ones. Models in this two-parametric class which remain viable appear to be close to the slow-roll ones, and their inflaton potentials are close to (but still different from) that of the natural inflation model. Permitted regions for the two model parameters are presented.


title: On the Weight Hierarchy of Locally Repairable Codes

abstract: An $(n,k,r)$ \emph{locally repairable code} (LRC) is an $[n,k,d]$ linear code where every code symbol can be repaired from at most $r$ other code symbols. An LRC is said to be optimal if the minimum distance attains the Singleton-like bound $d \le n-k-\lceil k/r \rceil +2$. The \emph{generalized Hamming weights} (GHWs) of linear codes are fundamental parameters which have many useful applications. Generally it is difficult to determine the GHWs of linear codes. In this paper, we study the GHWs of LRCs. Firstly, we obtain a generalized Singleton-like bound on the $i$-th $(1\le i \le k)$ GHWs of general $(n,k,r)$ LRCs. Then, it is shown that for an optimal $(n,k,r)$ LRC with $r \mid k$, its weight hierarchy can be completely determined, and the $i$-th GHW of an optimal $(n,k,r)$ LRC with $r \mid k$ attains the proposed generalized Singleton-like bound for all $1\le i \le k$. For an optimal $(n,k,r)$ LRC with $r \nmid k$, we give lower bounds on the GHWs of the LRC and its dual code. Finally, two general bounds on linear codes in terms of the GHWs are presented. Moreover, it is also shown that some previous results on the bounds of minimum distances of linear codes can also be explained or refined in terms of the GHWs.


title: A Platform for Automating Chaos Experiments

abstract: The Netflix video streaming system is composed of many interacting services. In such a large system, failures in individual services are not uncommon. This paper describes the Chaos Automation Platform, a system for running failure injection experiments on the production system to verify that failures in non-critical services do not result in system outages.


title: Semicontinuous Banach spaces for Schr\"odinger's Eq. with Dirac-$\delta'$ potential

abstract: Schr\'{o}dinger's equation with distributional $\delta$, or $\delta'$ potentials has been well studied in the past. There are challenges in simultaneously addressing some of the inherent issues of the system: The functional operator cannot exist entirely within the standard $L^2$ Hilbert spaces. On differentiable manifolds, the domain of the free kinetic energy operator is in the space of harmonic forms. Locally, by the Hodge decomposition theorem and the standard distributional calculus, the space of functionals of a $\delta$ or $\delta'$ potential must be orthogonal to the free kinetic energy operator. Restricting to semicontinuous topologies presents opportunities to address these, and other issues. We develop, in great detail, a formalism of Banach spaces with semicontinuous topologies, and their properties are extensively defined and studied. For $C(\overline{\mathbb{R}})$ functions, the spaces are indistinguishable. The semicontinuous analogs of the $L^p$ spaces, are nontrivial and result in a dense topologically continuous embedding of the semicontinuous $L^p$ spaces into the semicontinuous $C(\overline{\mathbb{R}})$ spaces. Here, certain classes of distributions may be inverted in terms of their primitive functions. Also many operators are inherently self adjoint. We define equivalence relations between the cohomology classes of distributions and derivatives of their associated primitives on local sections of $\overline{\mathbb{R}}$. Here Hamilton's equations are canonical, and define a connection on the fibers of the base space. Semicontinuity provides a resolution to the above domain and interaction problems, and easily integrable Feynman functional. We arrive at a compatible domain which is Krein ($\mathfrak{H}$) over disjoint components of $\overline{\mathbb{R}}$. The subspaces of $\mathfrak{H}$ are isomorphic to the semicontinuous Hilbert spaces of the Hamiltonian.


title: Parallel Optimization of Polynomials for Large-scale Problems in Stability and Control

abstract: In this thesis, we focus on some of the NP-hard problems in control theory. Thanks to the converse Lyapunov theory, these problems can often be modeled as optimization over polynomials. To avoid the problem of intractability, we establish a trade off between accuracy and complexity. We develop a sequence of tractable optimization problems - in the form of LPs and SDPs - whose solutions converge to the exact solution of the NP-hard problem. However, the computational and memory complexity of these LPs and SDPs grow exponentially with the progress of the sequence - meaning that improving the accuracy of the solutions requires solving SDPs with tens of thousands of decision variables and constraints. Setting up and solving such problems is a significant challenge. The existing optimization algorithms and software are only designed to use desktop computers or small cluster computers - machines which do not have sufficient memory for solving such large SDPs. This in fact is the reason we seek parallel algorithms for setting-up and solving large SDPs on supercomputers. We propose parallel algorithms for stability analysis of two classes of systems: 1) Linear systems with a large number of uncertain parameters; 2) Nonlinear systems defined by polynomial vector fields. First, we develop a distributed parallel algorithm which applies Polya's and Handelman's theorems to some variants of parameter-dependent Lyapunov inequalities with parameters defined over the simplex. The result is a sequence of SDPs which possess a block-diagonal structure. We then develop a parallel SDP solver which exploits this structure to map the computation, memory and communication to a distributed parallel environment. Numerical tests on a supercomputer demonstrate the ability of the algorithm to efficiently utilize hundreds and potentially thousands of processors and analyze systems with 100+ dimensional state-space.


title: Some asymptotic results for nonlinear Hawkes processes

abstract: Hawkes process is a class of simple point processes with self-exciting and clustering properties. Hawkes process has been widely applied in finance, neuroscience, social networks, criminology, seismology, and many other fields. In this paper, we study fluctuations, large deviations and moderate deviations nonlinear Hawkes processes in a new asymptotic regime, the large intensity function and the small exciting function regime. It corresponds to the large baseline intensity asymptotics for the linear case, and can also be interpreted as the asymptotics for the mean process of Hawkes processes on a large network.


title: Relay-Aided MIMO Cellular Networks Using Opposite Directional Interference Alignment

abstract: In this paper, we first propose interference alignment (IA) scheme for uplink transmission of multiple-input-mulitple-output (MIMO) cellular network with a help of relay which operates in halfduplex mode. The proposed scheme only requires global channel state information (CSI) knowledge at relay and no transmitter beamforming and time extension is required at user equipment (UE), which differs from the conventional IA schemes for cellular network. We derive the feasibility condition of the proposed scheme for the general network configuration and analyze the degrees-of-freedom (DoF) performance of the proposed IA scheme. Extension of proposed scheme for downlink and full-duplex network are further described in this paper. It is also shown that the same advantage as the uplink case can be obtained for downlink case through relay induced interfering multiple access channel (IMAC) and interfering broadcast channel (IBC) duality. Furthermore, full-duplex network is shown to have same advantages with half-duplex cases.


title: Blocking Self-avoiding Walks Stops Cyber-epidemics: A Scalable GPU-based Approach

abstract: Cyber-epidemics, the widespread of fake news or propaganda through social media, can cause devastating economic and political consequences. A common countermeasure against cyber-epidemics is to disable a small subset of suspected social connections or accounts to effectively contain the epidemics. An example is the recent shutdown of 125,000 ISIS-related Twitter accounts. Despite many proposed methods to identify such subset, none are scalable enough to provide high-quality solutions in nowadays billion-size networks. To this end, we investigate the Spread Interdiction problems that seek most effective links (or nodes) for removal under the well-known Linear Threshold model. We propose novel CPU-GPU methods that scale to networks with billions of edges, yet, possess rigorous theoretical guarantee on the solution quality. At the core of our methods is an $O(1)$-space out-of-core algorithm to generate a new type of random walks, called Hitting Self-avoiding Walks (HSAWs). Such a low memory requirement enables handling of big networks and, more importantly, hiding latency via scheduling of millions of threads on GPUs. Comprehensive experiments on real-world networks show that our algorithms provides much higher quality solutions and are several order of magnitude faster than the state-of-the art. Comparing to the (single-core) CPU counterpart, our GPU implementations achieve significant speedup factors up to 177x on a single GPU and 338x on a GPU pair.


title: On generalization of Bailey's identity involving product of generalized hypergeometric series

abstract: The aim of this research paper is to obtain explicit expressions of (i) $ {}_1F_1 \left[\begin{array}{c} \alpha \\ 2\alpha + i \end{array} ; x \right]. {}_1F_1\left[ \begin{array}{c} \beta \\ 2\beta + j \end{array} ; x \right]$ (ii) ${}_1F_1 \left[ \begin{array}{c} \alpha \\ 2\alpha - i \end{array} ; x \right] . {}_1F_1 \left[ \begin{array}{c} \beta \\ 2\beta - j \end{array} ; x \right]$ (iii) ${}_1F_1 \left[ \begin{array}{c} \alpha \\ 2\alpha + i \end{array} ; x \right] . {}_1F_1 \left[\begin{array}{c} \beta \\ 2\beta - j \end{array} ; x \right]$ in the most general form for any $i,j=0,1,2,\ldots$ For $i=j=0$, we recover well known and useful identity due to Bailey. The results are derived with the help of a well known Bailey's formula involving products of generalized hypergeometric series and generalization of Kummer's second transformation formulas available in the literature. A few interesting new as well as known special cases have also been given.


title: Cold collisions of heavy $^2\Sigma$ molecules with alkali-metal atoms in a magnetic field: Ab initio analysis and prospects for sympathetic cooling of SrOH$(^2\Sigma)$ by Li($^2$S)

abstract: We use accurate ab initio and quantum scattering calculations to explore the prospects for sympathetic cooling of the heavy molecular radical SrOH($^2\Sigma$) by ultracold Li atoms in a magnetic trap. A two-dimensional potential energy surface (PES) for the triplet electronic state of Li-SrOH is calculated ab initio using the partially spin-restricted coupled cluster method with single, double and perturbative triple excitations and a large correlation-consistent basis set. The highly anisotropic PES has a deep global minimum in the skewed Li-HOSr geometry with $D_e=4932$ cm$^{-1}$ and saddle points in collinear configurations. Our quantum scattering calculations predict low spin relaxation rates in fully spin-polarized Li+SrOH collisions with the ratios of elastic to inelastic collision rates well in excess of 100 over a wide range of magnetic fields (1-1000 G) and collision energies (10$^{-5}-0.1$~K) suggesting favorable prospects for sympathetic cooling of SrOH molecules with spin-polarized Li atoms in a magnetic trap. We find that spin relaxation in Li+SrOH collisions occurs via a direct mechanism mediated by the magnetic dipole-dipole interaction between the electron spins of Li and SrOH, and that the indirect (spin-rotation) mechanism is strongly suppressed. The upper limit to the Li+SrOH reaction rate coefficient calculated for the singlet PES using adiabatic capture theory is found to decrease from $4\times 10^{-10}$~cm$^3$/s to a limiting value of $3.5\times 10^{-10}$ cm$^3$/s with decreasing temperature from 0.1 K to 1 $\mu$K.


title: Approaching the Intrinsic Photoluminescence Linewidth in Transition Metal Dichalcogenide Monolayers

abstract: Excitonic states in monolayer transition metal dichalcogenides (TMDCs) have been the subject of extensive recent interest. Their intrinsic properties can, however, be obscured due to the influence of inhomogeneity in the external environment. Here we report methods for fabricating high quality TMDC monolayers with narrow photoluminescence (PL) linewidth approaching the intrinsic limit. We find that encapsulation in hexagonal boron nitride (h-BN) sharply reduces the PL linewidth, and that passivation of the oxide substrate by an alkyl monolayer further decreases the linewidth and also minimizes the charged exciton (trion) peak. The combination of these sample preparation methods results in much reduced spatial variation in the PL emission, with a full-width-at-half-maximum as low as 1.7 meV. Analysis of the PL line shape yields a homogeneous width of 1.43$\pm$0.08 meV and inhomogeneous broadening of 1.1$\pm$0.3 meV.


title: Multi-Sensor Control for Multi-Object Bayes Filters

abstract: Sensor management in multi-object stochastic systems is a theoretically and computationally challenging problem. This paper presents a novel approach to the multi-target multi-sensor control problem within the partially observed Markov decision process (POMDP) framework. We model the multi-object state as a labeled multi-Bernoulli random finite set (RFS), and use the labeled multi-Bernoulli filter in conjunction with minimizing a task-driven control objective function: posterior expected error of cardinality and state (PEECS). A major contribution is a guided search for multi-dimensional optimization in the multi-sensor control command space, using coordinate descent method. In conjunction with the Generalized Covariance Intersection method for multi-sensor fusion, a fast multi-sensor algorithm is achieved. Numerical studies are presented in several scenarios where numerous controllable (mobile) sensors track multiple moving targets with different levels of observability. The results show that our method works significantly faster than the approach taken by a state of art method, with similar tracking errors.


title: Data-driven polynomial ridge approximation using variable projection

abstract: Inexpensive surrogates are useful for reducing the cost of science and engineering studies involving large-scale, complex computational models with many input parameters. A ridge approximation is one class of surrogate that models a quantity of interest as a nonlinear function of a few linear combinations of the input parameters. When used in parameter studies (e.g., optimization or uncertainty quantification), ridge approximations allow the low dimensional structure to be exploited, reducing the effective dimension. We introduce a new, fast algorithm for constructing a ridge approximation where the nonlinear function is a polynomial. This polynomial ridge approximation is chosen to minimize least squared mismatch between the surrogate and the quantity of interest on a given set of inputs. Naively, this would require optimizing both the polynomial coefficients and the linear combination of weights; the latter of which define a low-dimensional subspace of the input space. However, given a fixed subspace the optimal polynomial can be found by solving a linear least-squares problem, and hence by using variable projection the polynomial can be implicitly found leaving an optimization problem over the subspace alone. We provide an algorithm that finds this polynomial ridge approximation by minimizing over the Grassmann manifold of low-dimensional subspaces using a Gauss-Newton method. We provide details of this optimization algorithm and demonstrate its performance on several numerical examples. Our Gauss-Newton method has superior theoretical guarantees and faster convergence than the alternating approach for polynomial ridge approximation earlier proposed by Constantine, Eftekhari, Hokanson, and Ward [https://doi.org/10.1016/j.cma.2017.07.038] that alternates between (i) optimizing the polynomial coefficients given the subspace and (ii) optimizing the subspace given the coefficients.


title: Robust Sparse Estimation Tasks in High Dimensions

abstract: In this paper we initiate the study of whether or not sparse estimation tasks can be performed efficiently in high dimensions, in the robust setting where an $\eps$-fraction of samples are corrupted adversarially. We study the natural robust version of two classical sparse estimation problems, namely, sparse mean estimation and sparse PCA in the spiked covariance model. For both of these problems, we provide the first efficient algorithms that provide non-trivial error guarantees in the presence of noise, using only a number of samples which is similar to the number required for these problems without noise. In particular, our sample complexities are sublinear in the ambient dimension $d$. Our work also suggests evidence for new computational-vs-statistical gaps for these problems (similar to those for sparse PCA without noise) which only arise in the presence of noise.


title: The Business of Height Pairings

abstract: In algebraic geometry there is the notion of a height pairing of algebraic cycles, which lies at the confluence of arithmetic, Hodge theory and topology. After explaining a motivating example situation, we introduce new directions in this subject.


title: Event-by-Event Simulations of Jet Modification Using the MATTER Event Generator

abstract: The modification of hard jets in the Quark Gluon Plasma (QGP) is studied using the MATTER event generator. Based on the higher twist formalism of energy loss, the MATTER event generator simulates the evolution of highly virtual partons through a medium. These partons sampled from an underlying PYTHIA kernel undergo splitting through a combination of vacuum and medium induced emission. The momentum exchange with the medium is simulated via the jet transport coefficient $\hat{q}$, which is assumed to scale with the entropy density at a given location in the medium. The entropy density is obtained from a relativistic viscous fluid dynamics simulation (VISH2+1D) in 2+1 space time dimensions. Results for jet and hadron observables are presented using an independent fragmentation model. These proceedings will focus on the physics input and simulation details of the MATTER event generator as compared to a variety of test observables.


title: Semantic Compression for Edge-Assisted Systems

abstract: A novel semantic approach to data selection and compression is presented for the dynamic adaptation of IoT data processing and transmission within "wireless islands", where a set of sensing devices (sensors) are interconnected through one-hop wireless links to a computational resource via a local access point. The core of the proposed technique is a cooperative framework where local classifiers at the mobile nodes are dynamically crafted and updated based on the current state of the observed system, the global processing objective and the characteristics of the sensors and data streams. The edge processor plays a key role by establishing a link between content and operations within the distributed system. The local classifiers are designed to filter the data streams and provide only the needed information to the global classifier at the edge processor, thus minimizing bandwidth usage. However, the better the accuracy of these local classifiers, the larger the energy necessary to run them at the individual sensors. A formulation of the optimization problem for the dynamic construction of the classifiers under bandwidth and energy constraints is proposed and demonstrated on a synthetic example.


title: An elliptic theory of indicial weights and applications to non-linear geometry problems

abstract: Given an elliptic operator $P$ on a non-compact manifold (with proper asymptotic conditions), there is a discrete set of numbers called indicial roots. It's known that $P$ is Fredholm between weighted Sobolev spaces if and only if the weight is not indicial. We show that an elliptic theory exists even when the weight is indicial. We also discuss some simple applications to Yang-Mills theory and minimal surfaces.


title: Hemingway: Modeling Distributed Optimization Algorithms

abstract: Distributed optimization algorithms are widely used in many industrial machine learning applications. However choosing the appropriate algorithm and cluster size is often difficult for users as the performance and convergence rate of optimization algorithms vary with the size of the cluster. In this paper we make the case for an ML-optimizer that can select the appropriate algorithm and cluster size to use for a given problem. To do this we propose building two models: one that captures the system level characteristics of how computation, communication change as we increase cluster sizes and another that captures how convergence rates change with cluster sizes. We present preliminary results from our prototype implementation called Hemingway and discuss some of the challenges involved in developing such a system.


title: Pulsar Wind Nebulae inside Supernova Remnants as Cosmic-Ray PeVatrons

abstract: We propose that cosmic-ray PeVatrons are pulsar wind nebulae (PWNe) inside supernova remnants (SNRs). The PWN initially expands into the freely expanding stellar ejecta. Then, the PWN catches up with the shocked region of the SNR, where particles can be slightly accelerated by the back and forth motion between the PWN and the SNR, and some particles diffuse into the PWN. Afterwards the PWN is compressed by the SNR, where the particles in the PWN are accelerated by the adiabatic compression. Using a Monte Carlo simulation, we show that particles accelerated by the SNR to 0.1 PeV can be reaccelerated to 1 PeV until the end of the PWN compression.


title: Multiplicities of Character Values of Binary Sidel'nikov-Lempel-Cohn-Eastman Sequences

abstract: Binary Sidel'nikov-Lempel-Cohn-Eastman sequences (or SLCE sequences) over F 2 have even period and almost perfect autocorrelation. However, the evaluation of the linear complexity of these sequences is really difficult. In this paper, we continue the study of [1]. We first express the multiple roots of character polynomials of SLCE sequences into certain kinds of Jacobi sums. Then by making use of Gauss sums and Jacobi sums in the "semiprimitive" case, we derive new divisibility results for SLCE sequences.


title: A variant of Gromov's problem on H\"older equivalence of Carnot groups

abstract: It is unknown if there exists a locally $\alpha$-H\"older homeomorphism $f:\mathbb{R}^3\to \mathbb{H}^1$ for any $\frac{1}{2}< \alpha\le \frac{2}{3}$, although the identity map $\mathbb{R}^3\to \mathbb{H}^1$ is locally $\frac{1}{2}$-H\"older. More generally, Gromov asked: Given $k$ and a Carnot group $G$, for which $\alpha$ does there exist a locally $\alpha$-H\"older homeomorphism $f:\mathbb{R}^k\to G$? Here, we equip a Carnot group $G$ with the Carnot-Carath\'eodory metric. In 2014, Balogh, Hajlasz, and Wildrick considered a variant of this problem. These authors proved that if $k>n$, there does not exist an injective, $(\frac{1}{2}+)$-H\"older mapping $f:\mathbb{R}^k\to \mathbb{H}^n$ that is also locally Lipschitz as a mapping into $\mathbb{R}^{2n+1}$. For their proof, they use the fact that $\mathbb{H}^n$ is purely $k$-unrectifiable for $k>n$. In this paper, we will extend their result from the Heisenberg group to model filiform groups and Carnot groups of step at most three. We will now require that the Carnot group is purely $k$-unrectifiable. The main key to our proof will be showing that $(\frac{1}{2}+)$-H\"older maps $f:\mathbb{R}^k\to G$ that are locally Lipschitz into Euclidean space, are weakly contact. Proving weak contactness in these two settings requires understanding the relationship between the algebraic and metric structures of the Carnot group. We will use coordinates of the first and second kind for Carnot groups.


title: Comments on Entanglement Propagation

abstract: We extend our work on entanglement propagation following a local quench in 2+1 dimensional holographic conformal field theories. We find that entanglement propagates along an emergent lightcone, whose speed of propagation $v_E$ seems distinct from other measures of quantum information spreading. We compare the relations we find to information and hydrodynamic velocities in strongly coupled 2+1 dimensional theories. While early-time entanglement velocities corresponding to small entangling regions are numerically close to the butterfly velocity, late-time entanglement velocities for large regions show less regularity. We also generalize and extend our previous results regarding the late-time decay of the entanglement entropy back to its equilibrium value.


title: Cosine Normalization: Using Cosine Similarity Instead of Dot Product in Neural Networks

abstract: Traditionally, multi-layer neural networks use dot product between the output vector of previous layer and the incoming weight vector as the input to activation function. The result of dot product is unbounded, thus increases the risk of large variance. Large variance of neuron makes the model sensitive to the change of input distribution, thus results in poor generalization, and aggravates the internal covariate shift which slows down the training. To bound dot product and decrease the variance, we propose to use cosine similarity or centered cosine similarity (Pearson Correlation Coefficient) instead of dot product in neural networks, which we call cosine normalization. We compare cosine normalization with batch, weight and layer normalization in fully-connected neural networks as well as convolutional networks on the data sets of MNIST, 20NEWS GROUP, CIFAR-10/100 and SVHN. Experiments show that cosine normalization achieves better performance than other normalization techniques.


title: Equation of state for hybrid stars with strangeness

abstract: Considering the mass constraint from the resent pulsar observations, we study the properties of neutron stars including hyperons and quarks explicitly. Using the chiral quark-meson coupling model with relativistic Hartree-Fock approximation, the equation of state (EoS) for hadronic matter is calculated by taking into account the strange ($\sigma^{\ast}$ and $\phi$) mesons as well as the light non-strange ($\sigma$, $\omega$, $\vec{\rho}$, and $\vec{\pi}$) mesons in SU(3) flavor symmetry. On the other hand, the EoS for quark matter is constructed with the simple MIT bag or the flavor-SU(3) Nambu-Jona-Lasinio model, and we investigate the effect of the hadron-quark coexistence on the neutron-star properties, imposing smooth crossover or Gibbs criterion for chemical equilibrium. The mass-radius relation of a neutron star, as well as physical quantities such as EoSs, particle fractions, and the speed of sound in matter are presented. We find that, in order to prevent the quark appearance at very low densities, the stiff hadronic EoS should be required under both of the hadron-quark crossover and the first-order phase transition.


title: Acoustic metasurfaces for scattering-free anomalous reflection and refraction

abstract: Manipulation of acoustic wavefronts by thin and planar devices, known as metasurfaces, has been extensively studied, in view of many important applications. Reflective and refractive metasurfaces are designed using the generalized reflection and Snell's laws, which tell that local phase shifts at the metasurface supply extra momentum to the wave, presumably allowing arbitrary control of reflected or transmitted waves. However, as it has been recently shown for the electromagnetic counterpart, conventional metasurfaces based on the generalized laws of reflection and refraction have important drawbacks in terms of power efficiency. This work presents a new synthesis method of acoustic metasurfaces for anomalous reflection and transmission that overcomes the fundamental limitations of conventional designs, allowing full control of acoustic energy flow. The results show that different mechanisms are necessary in the reflection and transmission scenarios for ensuring perfect performance. Metasurfaces for anomalous reflection require non-local response, which allows energy channeling along the metasurface. On other hand, for perfect manipulation of anomalously transmitted waves, local and non-symmetric response is required. These conclusions are interpreted through appropriate surface impedance models which are used to find possible physical implementations of perfect metasurfaces in each scenario. We hope that this advance in the design of acoustic metasurfaces opens new avenues not only for perfect anomalous reflection and transmission but also for realizing more complex functionalities, such as focusing, self-bending or vortex generation.


title: Characterization of 1-Tough Graphs using Factors

abstract: For a graph $G$, let $odd(G)$ and $\omega(G)$ denote the number of odd components and the number of components of $G$, respectively. Then it is well-known that $G$ has a 1-factor if and only if $odd(G-S)\le |S|$ for all $S\subset V(G)$. Also it is clear that $odd(G-S) \le \omega(G-S)$. In this paper we characterize a 1-tough graph $G$, which satisfies $\omega(G-S) \le |S|$ for all $\emptyset \ne S \subset V(G)$, using an $H$-factor of a set-valued function $H:V(G) \to \{ \{1\}, \{0,2\} \}$. Moreover, we generalize this characterization to a graph that satisfies $\omega(G-S) \le f(S)$ for all $\emptyset \ne S \subset V(G)$, where $f:V(G) \to \{1,3,5, \ldots\}$.


title: The Complexity of All $(g,f)$-Factor Problem

abstract: Let $G$ be a graph with vertex set $V$ and let $g, f : V\rightarrow \mathbb{Z}^+$ be two functions such that $g\le f$. We say that $G$ has all $(g, f )$-factors if $G$ has an $h$-factor for every $h: V\rightarrow \mathbb{Z}^+$ such that $g(v)\le h(v)\le f (v)$ for every $v\in V$ and $\sum_{v\in{V}}h(v)\equiv 0\pmod 2$. Two decades ago, Niessen derived from Tutte's $f$-factor theorem a similar characterization for the property of graphs having all $(g, f )$-factors and asked whether there is a polynomial time algorithm for testing whether a graph $G$ has all $(g, f )$-factors (A characterization of graphs having all $(g, f )$-Factors, \emph{J. Combin. Theory, Ser. B}, \textbf{72} (1998), 152--156). In this paper, we show that it is NP-hard to determine whether a graph $G$ has all $(g,f)$-factors, which gives a negative answer to the question of Niessen.


title: Spin-accumulation capacitance and its application to magnetoimpedance

abstract: It has been known that spin-dependent capacitances usually coexist with geometric capacitances in a magnetic multilayer. However, the charge and energy storage of the capacitance due to spin accumulation has not been fully understood. Here, we resolve this problem starting from the charge storage in the spin degree of freedom: spin accumulation manifests itself as an excess of electrons in one spin channel and an equal deficiency in the other under the quasi-neutrality condition. This enables us to model the two spin channels as the two plates of a capacitor. Taking a ferromagnet/nonmagnet junction as an example and using a method similar to that for treating quantum capacitance, we find that a spin-accumulation (SA) capacitance can be introduced for each layer to measure its ability to store spins. A spatial charge storage is not essential for the SA capacitor and the energy stored in it is the splitting energy of the spin-dependent chemical potentials instead of the electrostatic energy. The SA capacitance is essentially a quantum capacitance due to spin accumulation on the scale of the spin-diffusion length. The SA capacitances can be used to reinterpret the imaginary part of the low-frequency magnetoimpedance.


title: Gate Tunable Magneto-resistance of Ultra-Thin WTe2 Devices

abstract: In this work, the magneto-resistance (MR) of ultra-thin WTe2/BN heterostructures far away from electron-hole equilibrium is measured. The change of MR of such devices is found to be determined largely by a single tunable parameter, i.e. the amount of imbalance between electrons and holes. We also found that the magnetoresistive behavior of ultra-thin WTe2 devices is well-captured by a two-fluid model. According to the model, the change of MR could be as large as 400,000%, the largest potential change of MR among all materials known, if the ultra-thin samples are tuned to neutrality when preserving the mobility of 167,000 cm2V-1s-1 observed in bulk samples. Our findings show the prospects of ultra-thin WTe2 as a variable magnetoresistance material in future applications such as magnetic field sensors, information storage and extraction devices, and galvanic isolators. The results also provide important insight into the electronic structure and the origin of the large MR in ultra-thin WTe2 samples.


title: Production of 82Se enriched Zinc Selenide (ZnSe) crystals for the study of neutrinoless double beta decay

abstract: High purity Zinc Selenide (ZnSe) crystals are produced starting from elemental Zn and Se to be used for the search of the neutrinoless double beta decay (0{\nu}DBD) of 82Se. In order to increase the number of emitting nuclides, enriched 82Se is used. Dedicated production lines for the synthesis and conditioning of the Zn82Se powder in order to make it suitable for crystal growth were assembled compliant with radio-purity constraints specific to rare event physics experiments. Besides routine check of impurities concentration, high sensitivity measurements are made for radio-isotope concentrations in raw materials, reactants, consumables, ancillaries and intermediary products used for ZnSe crystals production. Indications are given on the crystals perfection and how it is achieved. Since very expensive isotopically enriched material (82Se) is used, a special attention is given for acquiring the maximum yield in the mass balance of all production stages. Production and certification protocols are presented and resulting ready-to-use Zn82Se crystals are described.


title: From Photo Streams to Evolving Situations

abstract: Photos are becoming spontaneous, objective, and universal sources of information. This paper develops evolving situation recognition using photo streams coming from disparate sources combined with the advances of deep learning. Using visual concepts in photos together with space and time information, we formulate the situation detection into a semi-supervised learning framework and propose new graph-based models to solve the problem. To extend the method for unknown situations, we introduce a soft label method which enables the traditional semi-supervised learning framework to accurately predict predefined labels as well as effectively form new clusters. To overcome the noisy data which degrades graph quality, leading to poor recognition results, we take advantage of two kinds of noise-robust norms which can eliminate the adverse effects of outliers in visual concepts and improve the accuracy of situation recognition. Finally, we demonstrate the idea and the effectiveness of the proposed model on Yahoo Flickr Creative Commons 100 Million.


title: Complexity of Possibly-gapped Histogram and Analysis of Histogram (ANOHT)

abstract: Without unrealistic continuity and smoothness assumptions on a distributional density of one dimensional dataset, constructing an authentic possibly-gapped histogram becomes rather complex. The candidate ensemble is described via a two-layer Ising model, and its size is shown to grow exponentially. This exponential complexity makes any exhaustive search in-feasible and all boundary parameters local. For data compression via Uniformity, the decoding error criterion is nearly independent of sample size. These characteristics nullify statistical model selection techniques, such as Minimum Description Length (MDL). Nonetheless practical and nearly optimal solutions are algorithmically computable. A data-driven algorithm is devised to construct such histograms along the branching hierarchy of a Hierarchical Clustering tree. Such resultant histograms naturally manifest data's physical information contents: deterministic structures of bin-boundaries coupled with stochastic structures of Uniformity within each bin. Without enforcing unrealistic Normality and constant variance assumptions, an application of possibly-gapped histogram is devised, called analysis of Histogram (ANOHT), to replace Analysis of Variance (ANOVA). Its potential applications are foreseen in digital re-normalization schemes and associative pattern extraction among features of heterogeneous data types. Thus constructing possibly-gapped histograms becomes a prerequisite for knowledge discovery, via exploratory data analysis and unsupervised Machine Learning.


title: Mobility Increases the Data Offloading Ratio in D2D Caching Networks

abstract: Caching at mobile devices, accompanied by device-to-device (D2D) communications, is one promising technique to accommodate the exponentially increasing mobile data traffic. While most previous works ignored user mobility, there are some recent works taking it into account. However, the duration of user contact times has been ignored, making it difficult to explicitly characterize the effect of mobility. In this paper, we adopt the alternating renewal process to model the duration of both the contact and inter-contact times, and investigate how the caching performance is affected by mobility. The data offloading ratio, i.e., the proportion of requested data that can be delivered via D2D links, is taken as the performance metric. We first approximate the distribution of the communication time for a given user by beta distribution through moment matching. With this approximation, an accurate expression of the data offloading ratio is derived. For the homogeneous case where the average contact and inter-contact times of different user pairs are identical, we prove that the data offloading ratio increases with the user moving speed, assuming that the transmission rate remains the same. Simulation results are provided to show the accuracy of the approximate result, and also validate the effect of user mobility.


title: The system of ionized gas dynamics

abstract: The aim of this paper is to study a system of three equations for ionized gas dynamics at high temperature, in one spatial dimension. In addition to the mass density, pressure and particle velocity, a further quantity is needed, namely, the degree of ionization. The system is supplemented by the first and second law of thermodynamics and by an equation of state; all of them involve the degree of ionization. At last, under the assumption of thermal equilibrium, the system is closed by requiring Saha's ionization equation. The geometric properties of the system are rather complicated: in particular, we prove the loss of convexity (genuine nonlinearity) for both forward and backward characteristic fields, and hence the loss of concavity of the physical entropy. This takes place in a small bounded region, which we are able to characterize by numerical estimates on the state functions. The structure of shock waves is also studied by a detailed analysis of the Hugoniot locus, which will be used in a forthcoming paper to study the shock tube problem.


title: Phase Diagram of Restricted Boltzmann Machines and Generalised Hopfield Networks with Arbitrary Priors

abstract: Restricted Boltzmann Machines are described by the Gibbs measure of a bipartite spin glass, which in turn corresponds to the one of a generalised Hopfield network. This equivalence allows us to characterise the state of these systems in terms of retrieval capabilities, both at low and high load. We study the paramagnetic-spin glass and the spin glass-retrieval phase transitions, as the pattern (i.e. weight) distribution and spin (i.e. unit) priors vary smoothly from Gaussian real variables to Boolean discrete variables. Our analysis shows that the presence of a retrieval phase is robust and not peculiar to the standard Hopfield model with Boolean patterns. The retrieval region is larger when the pattern entries and retrieval units get more peaked and, conversely, when the hidden units acquire a broader prior and therefore have a stronger response to high fields. Moreover, at low load retrieval always exists below some critical temperature, for every pattern distribution ranging from the Boolean to the Gaussian case.

categories: cond-mat.dis-nn cs.LG physics.data-an stat.ML


title: Resonant Scattering Characteristics of Homogeneous Dielectric Sphere

abstract: In the present article the classical problem of electromagnetic scattering by a single homogeneous sphere is revisited. Main focus is the study of the scattering behavior as a function of the material contrast and the size parameters for all electric and magnetic resonances of a dielectric sphere. Specifically, the Pad\'e approximants are introduced and utilized as an alternative system expansion of the Mie coefficients. Low order Pad\'e approximants can give compact and physically insightful expressions for the scattering system and the enabled dynamic mechanisms. Higher order approximants are used for predicting accurately the resonant pole spectrum. These results are summarized into general pole formulae, covering up to fifth order magnetic and forth order electric resonances of a small dielectric sphere. Additionally, the connection between the radiative damping process and the resonant linewidth is investigated. The results obtained reveal the fundamental connection of the radiative damping mechanism with the maximum width occurring for each resonance. Finally, the suggested system ansatz is used for studying the resonant absorption maximum through a circuit-inspired perspective.


title: Studies on spaces of initial conditions for nonautonomous mappings of the plane

abstract: We study nonautonomous mappings of the plane by means of spaces of initial conditions. First we introduce the notion of a space of initial conditions for nonautonomous systems and we study the basic properties of general equations that have spaces of initial conditions. Then, we consider the minimization of spaces of initial conditions for nonautonomous systems and we show that if a nonautonomous mapping of the plane with a space of initial conditions, and unbounded degree growth, has zero algebraic entropy, then it must be one of the discrete Painlev\'{e} equations in the Sakai classification.


title: Signatures of monic polynomials

abstract: To a univariate monic polynomial is attached a special planar forest that is called the picture of the polynomial. Isotopy classes of pictures are called signatures. All combinatorially possible signatures are realized and spaces of polynomials realizing a given signature are contractible. A finite cell complex for the cohomology of the braid groups is obtained.


title: $L^p$-Analysis of the Hodge--Dirac operator associated with Witten Laplacians on complete Riemannian manifolds

abstract: We prove $R$-bisectoriality and boundedness of the $H^\infty$-functional calculus in $L^p$ for all $1<p<\infty$ for the Hodge-Dirac operator associated with Witten Laplacians on complete Riemannian manifolds with non-negative Bakry-Emery Ricci curvature on $k$-forms.


title: X-ray study of extended emission around M86 observed with Suzaku

abstract: We analyzed the Suzaku data of M86 and its adjacent regions to study the extended emission around it. The M86 core, the plume, and the tail extending toward the northwest were clearly detected, as well as the extended halo around them. From the position angle $\sim45\deg$ to $\sim275\deg$, the surface brightness distribution of the core and the extended halo was represented relatively well with a single $\beta$-model of $\beta\sim0.5$ up to 15'-20'. The X-ray spectra of the core was represented with a two-temperature model of $kT\sim 0.9$ keV and $\sim0.6$ keV. The temperatures of the core and the halo have a positive gradient in the center, and reach the maximum of $kT\sim1.0$ keV at $r\sim7'$, indicating that the halo gas is located in a larger scale potential structure than that of the galaxy. The temperatures of the plume and the tail were $0.86\pm0.01$ keV and $1.00\pm0.01$ keV. We succeeded in determining the abundances of $\alpha$-element separately, for the core, the plume, the tail, and the halo for the first time. Abundance ratios with respect to Fe were consistent with the solar ratios everywhere, except for Ne. The abundance of Fe was $\sim0.7$ in the core and in the plume, while that in the tail was $\sim1.0$, but the difference was not significant considering the uncertainties of the ICM. The abundance of the halo was almost the same up to $r\sim10'$, and then it becomes significantly smaller (0.2-0.3) at $r\gtrsim10'$, indicating the gas with low metal abundance still remains in the outer halo. From the surface brightness distribution, we estimated the gas mass ($\sim3\times10^{10}M_\odot$) and the dynamical mass ($\sim3\times10^{12}M_\odot$) in $r<100$~kpc. The gas mass to the dynamical mass ratio was $10^{-3}$-$10^{-2}$, suggesting a significant fraction of the halo gas has been stripped.


title: Memory Efficient Max Flow for Multi-label Submodular MRFs

abstract: Multi-label submodular Markov Random Fields (MRFs) have been shown to be solvable using max-flow based on an encoding of the labels proposed by Ishikawa, in which each variable $X_i$ is represented by $\ell$ nodes (where $\ell$ is the number of labels) arranged in a column. However, this method in general requires $2\,\ell^2$ edges for each pair of neighbouring variables. This makes it inapplicable to realistic problems with many variables and labels, due to excessive memory requirement. In this paper, we introduce a variant of the max-flow algorithm that requires much less storage. Consequently, our algorithm makes it possible to optimally solve multi-label submodular problems involving large numbers of variables and labels on a standard computer.


title: Existence and stability of circular orbits in general static and spherically symmetric spacetimes

abstract: The existence and stability of circular orbits (CO) in static and spherically symmetric (SSS) spacetime are important because of their practical and potential usefulness. In this paper, using the fixed point method, we first prove a necessary and sufficient condition on the metric function for the existence of timelike COs in SSS spacetimes. After analyzing the asymptotic behavior of the metric, we then show that asymptotic flat SSS spacetime that corresponds to a negative Newtonian potential at large $r$ will always allow the existence of CO. The stability of the CO in a general SSS spacetime is then studied using the Lyapunov exponent method. Two sufficient conditions on the (in)stability of the COs are obtained. For null geodesics, a sufficient condition on the metric function for the (in)stability of null CO is also obtained. We then illustrate one powerful application of these results by showing that an SU(2) Yang-Mills-Einstein SSS spacetime whose metric function is not known, will allow the existence of timelike COs. We also used our results to assert the existence and (in)stabilities of a number of known SSS metrics.


title: The HTL resumed propagators in the light cone gauge

abstract: The expression of the HTL resumed gluon propagator in the light cone gauge is derived. In the real time mechanism, using the Mandelstam Leibbrant prescription of $(n\cdot K)^{-1}$, we calculate the transverse and longitudinal parts of the gluon HTL self-energy and prove the transverse and longitudinal parts do not have divergence. We also calculate the quark self energy in the HTL approximation, and find it gauge independent. We analytically calculate the damping rates of the hard quark and gluon with this HTL resumed gluon propagator.


title: Learning Spatial Regularization with Image-level Supervisions for Multi-label Image Classification

abstract: Multi-label image classification is a fundamental but challenging task in computer vision. Great progress has been achieved by exploiting semantic relations between labels in recent years. However, conventional approaches are unable to model the underlying spatial relations between labels in multi-label images, because spatial annotations of the labels are generally not provided. In this paper, we propose a unified deep neural network that exploits both semantic and spatial relations between labels with only image-level supervisions. Given a multi-label image, our proposed Spatial Regularization Network (SRN) generates attention maps for all labels and captures the underlying relations between them via learnable convolutions. By aggregating the regularized classification results with original results by a ResNet-101 network, the classification performance can be consistently improved. The whole deep neural network is trained end-to-end with only image-level annotations, thus requires no additional efforts on image annotations. Extensive evaluations on 3 public datasets with different types of labels show that our approach significantly outperforms state-of-the-arts and has strong generalization capability. Analysis of the learned SRN model demonstrates that it can effectively capture both semantic and spatial relations of labels for improving classification performance.


title: On the rate of convergence in the central limit theorem for hierarchical Laplacian

abstract: Let $(X,d)$ be a proper ultrametric space. Given a measure $m$ on $X$ and a function $C(B)$ defined on the set of all non-singleton balls $B$ we consider the hierarchical Laplacian $L=L_{C}$. Choosing a sequence $\{\varepsilon (B)\}$ of i.i.d. random variables we define the perturbed function $C(B,\omega )$ and the perturbed hierarchical Laplacian $L^{\omega }=L_{C(\omega )}.$ We study the arithmetic means $\overline{\lambda }(\omega )$ of the $L^{\omega }$-eigenvalues. Under some mild assumptions the normalized arithmetic means $\big( \overline{\lambda }-\mathbb{E}\overline{\lambda }\big) /\sigma \big( \overline{\lambda }\big) $ converge in law to the standard normal distribution. In this note we study convergence in the total variation distance and estimate the rate of convergence.


title: Non-emptiness of Brill-Noether Loci over very general quintic hypersurface

abstract: In this article we study Brill-Noether loci of moduli space of stable bundles over smooth surfaces. We define Petri map as an analogy with the case of curves. We show the non-emptiness of certain Brill-Noether loci over very general quintic hypersurface in $\mathbb{P}^3$, and use the Petri map to produce components of expected dimension.


title: Schoenberg Representations and Gramian Matrices of Mat\'ern Functions

abstract: We represent Mat\'ern functions in terms of Schoenberg's integrals which ensure the positive definiteness and prove the systems of translates of Mat\'ern functions form Riesz sequences in $L^2(\R^n)$ or Sobolev spaces. Our approach is based on a new class of integral transforms that generalize Fourier transforms for radial functions. We also consider inverse multi-quadrics and obtain similar results.


title: Reply to the "Comment on Fermion production in a magnetic field in a de Sitter Universe"

abstract: In this paper we study the problem of gauge invariance of the first order transition amplitudes in de Sitter QED in Coulomb gauge. We consider the gauge transformations which preserve the Coulomb gauge, that contain the gradient of the gauge function. The final results prove that the first order transition amplitudes do not change at a gradient transformation of the vector potential because the only allowed transformation is $\Lambda=0$. Our results suggest that the remarks made in the comment by Nicolaevici and Farkas {\em Phys. Rev. D} \textbf{95}, are not directly applicable to the results in our paper since their proposed gauge transformations do not preserve the Coulomb gauge.


title: Compactly supported reproducing kernels for $L^2$-based Sobolev spaces and Hankel-Schoenberg transforms

abstract: We exhibit three classes of compactly supported functions which provide reproducing kernels for the Sobolev spaces $H^\delta(\R^d)$ of arbitrary order $\,\delta>d/2.\,$ Our method of construction is based on a new class of oscillatory integral transforms that incorporate radial Fourier transforms and Hankel transforms.


title: A dynamical stability limit for the charge density wave in K0.3MoO3

abstract: We study the response of the one-dimensional charge density wave in K0.3MoO3 to different types of excitation with femtosecond optical pulses. We compare the response to direct excitation of the lattice at mid-infrared frequencies with that to the injection of quasi-particles across the low-energy charge density wave gap and to charge transfer excitations in the near infrared. For all three cases, we observe a fluence threshold above which the amplitude-mode oscillation frequency is softened and the mode becomes increasingly damped. We show that all the data can be collapsed onto a universal curve in which the melting of the charge density wave occurs abruptly at a critical lattice excursion. These data highlight the existence of a universal stability limit for a charge density wave, reminiscent of the empirical Lindemann criterion for the stability of a crystal lattice.


title: Flows and Shocks: Some Recent Developments in Symbiotic Star and Nova Research

abstract: There have been several surprising developments in our understanding of symbiotic binary stars and nova eruptions over the last decade or so based on multiwavelength data. For example, symbiotic stars without shell burning on their white dwarfs have been revealed through their X-ray emission, UV excess, and UV variability. These purely accretion powered symbiotic stars have much weaker optical emission lines than those with shell burning, and are therefore harder to discover. Yet they may be as numerous as symbiotic stars with shell burning. Interestingly, both types of symbiotic stars are capable of driving strong outflows, leading to colliding wind X-ray emission and spatially resolved jets. For nova eruptions, the most surprising discovery has been that they are capable of particle acceleration as evidenced by Fermi detection of novae as transient GeV gamma-ray sources. For nova eruptions in cataclysmic variables, this implicates internal shocks, between a slow, dense outflow and a fast outflow or wind. Other signatures of shocks include thermal X-rays and non-thermal radio emissions, and a substantial fraction of optical emission may be shock-powered in the early phase of novae. Radio (V959 Mon) and HST (V959 Mon and T Pyx) images of nova shells within a few years of their respective eruptions suggest that the ejecta from a nova may commonly consist of an equatorial ring and a bipolar outflow.


title: An Examination of the Benefits of Scalable TTI for Heterogeneous Traffic Management in 5G Networks

abstract: The rapid growth in the number and variety of connected devices requires 5G wireless systems to cope with a very heterogeneous traffic mix. As a consequence, the use of a fixed TTI during transmission is not necessarily the most efficacious method when heterogeneous traffic types need to be simultaneously serviced.This work analyzes the benefits of scheduling based on exploiting scalable TTI, where the channel assignment and the TTI duration are adapted to the deadlines and requirements of different services. We formulate an optimization problem by taking individual service requirements into consideration. We then prove that the optimization problem is NP-hard and provide a heuristic algorithm, which provides an effective solution to the problem. Numerical results show that our proposed algorithm is capable of finding near-optimal solutions to meet the latency requirements of mission critical communication services, while providing a good throughput performance for mobile broadband services.


title: $\delta$-Greedy $t$-spanner

abstract: We introduce a new geometric spanner, $\delta$-Greedy, whose construction is based on a generalization of the known Path-Greedy and Gap-Greedy spanners. The $\delta$-Greedy spanner combines the most desirable properties of geometric spanners both in theory and in practice. More specifically, it has the same theoretical and practical properties as the Path-Greedy spanner: a natural definition, small degree, linear number of edges, low weight, and strong $(1+\varepsilon)$-spanner for every $\varepsilon>0$. The $\delta$-Greedy algorithm is an improvement over the Path-Greedy algorithm with respect to the number of shortest path queries and hence with respect to its construction time. We show how to construct such a spanner for a set of $n$ points in the plane in $O(n^2 \log n)$ time. The $\delta$-Greedy spanner has an additional parameter, $\delta$, which indicates how close it is to the Path-Greedy spanner on the account of the number of shortest path queries. For $\delta = t$ the output spanner is identical to the Path-Greedy spanner, while the number of shortest path queries is, in practice, linear. Finally, we show that for a set of $n$ points placed independently at random in a unit square the expected construction time of the $\delta$-Greedy algorithm is $O(n \log n)$. Our analysis indicates that the $\delta$-Greedy spanner gives the best results among the known spanners of expected $O(n \log n)$ time for random point sets. Moreover, the analysis implies that by setting $\delta = t$, the $\delta$-Greedy algorithm provides a spanner identical to the Path-Greedy spanner in expected $O(n \log n)$ time.


title: Reducing the Computational Complexity of Multicasting in Large-Scale Antenna Systems

abstract: In this paper, we study the physical layer multicasting to multiple co-channel groups in large-scale antenna systems. The users within each group are interested in a common message and different groups have distinct messages. In particular, we aim at designing the precoding vectors solving the so-called quality of service (QoS) and weighted max-min fairness (MMF) problems, assuming that the channel state information is available at the base station (BS). To solve both problems, the baseline approach exploits the semidefinite relaxation (SDR) technique. Considering a BS with $N$ antennas, the SDR complexity is more than $\mathcal{O}(N^{6})$, which prevents its application in large-scale antenna systems. To overcome this issue, we present two new classes of algorithms that, not only have significantly lower computational complexity than existing solutions, but also largely outperform the SDR based methods. Moreover, we present a novel duality between transformed versions of the QoS and the weighted MMF problems. The duality explicitly determines the solution to the weighted MMF problem given the solution to the QoS problem, and vice versa. Numerical results are used to validate the effectiveness of the proposed solutions and to make comparisons with existing alternatives under different operating conditions.


title: Excellent extensions and homological conjectures

abstract: In this paper, we introduce the notion of excellent extension of rings. Let $\Gamma$ be an excellent extension of an artin algebra $\Lambda$, we prove that $\Lambda$ satisfies the Gorenstein symmetry conjecture (resp. finitistic dimension conjecture, Auslander-Gorenstein conjecture, Nakayama conjecture) if and only if so does $\Gamma$. As a special case of excellent extensions, if $G$ is a finite group whose order is invertible in $\Lambda$ acting on $\Lambda$ and $\Lambda$ is $G$-stable, we prove that if the skew group algebras $\Lambda G$ satisfies strong Nakayama conjecture (resp. generalized Nakayama conjecture), then so does $\Lambda$.


title: Quantum criticality of spinons

abstract: The free fermion nature of interacting spins in one dimensional (1D) spin chains still lacks a rigorous study. In this letter we show that the length-$1$ spin strings significantly dominate critical properties of spinons, magnons and free fermions in the 1D antiferromagnetic spin-1/2 chain. Using the Bethe ansatz solution we analytically calculate exact scaling functions of thermal and magnetic properties of the model, providing a rigorous understanding of the quantum criticality of spinons. It turns out that the double peaks in specific heat elegantly mark two crossover temperatures fanning out from the critical point, indicating three quantum phases: the Tomonaga-Luttinger liquid (TLL), quantum critical and fully polarized ferromagnetic phases. For the TLL phase, the Wilson ratio $R_W=4K_s$ remains almost temperature-independent, here $K_s$ is the Luttinger parameter. Furthermore, applying our results we precisely determine the quantum scalings and critical exponents of all magnetic properties in the ideal 1D spin-1/2 antiferromagnet Cu(C${}_4$H${}_4$N${}_2$)(NO${}_3$)${}_2$ recently studied in Phys. Rev. Lett. {\bf 114}, 037202 (2015)]. We further find that the magnetization peak used in experiments is not a good quantity to map out the finite temperature TLL phase boundary.


title: Classical molecular dynamics simulations of fusion and fragmentation in fullerene-fullerene collisions

abstract: We present the results of classical molecular dynamics simulations of collision-induced fusion and fragmentation of C$_{60}$ fullerenes, performed by means of the MBN Explorer software package. The simulations provide information on structural differences of the fused compound depending on kinematics of the collision process. The analysis of fragmentation dynamics at different initial conditions shows that the size distributions of produced molecular fragments are peaked for dimers, which is in agreement with a well-established mechanism of C$_{60}$ fragmentation via preferential C$_2$ emission. Atomic trajectories of the colliding particles are analyzed and different fragmentation patterns are observed and discussed. On the basis of the performed simulations, characteristic time of C$_2$ emission is estimated as a function of collision energy. The results are compared with experimental time-of-flight distributions of molecular fragments and with earlier theoretical studies. Considering the widely explored case study of C$_{60}$--C$_{60}$ collisions, we demonstrate broad capabilities of the MBN Explorer software, which can be utilized for studying collisions of a broad variety of nanoscale and biomolecular systems by means of classical molecular dynamics.


title: An Early Diagnostics of the Geoeffectiveness of Solar Eruptions from Photospheric Magnetic Flux Observations: The Transition from SOHO to SDO

abstract: In our previous articles (Chertok et al.: 2013, Solar Phys. 282, 175, and 2015, Solar Phys. 290, 627), we presented a preliminary tool for the early diagnostics of the geoeffectiveness of solar eruptions based on the estimate of the total unsigned line-of-sight photospheric magnetic flux in accompanying extreme-ultraviolet arcades and dimmings. This tool was based on the analysis of eruptions observed in 1996-2005 with SOHO/EIT and MDI. Empirical relationships were obtained to estimate the probable importance of upcoming space weather disturbances caused by an eruption, which just occurred, without data on the associated coronal mass ejections. It was possible to estimate the intensity of a non-recurrent geomagnetic storm (GMS) and Forbush decrease (FD), as well as their onset and peak times. After 2010-2011, data on solar eruptions are obtained with SDO/AIA and HMI. We use relatively short intervals of overlapping EIT-AIA and MDI-HMI detailed observations and a number of large eruptions over the next five years with the 12-hour cadence EIT images to adapt the SOHO diagnostic tool to SDO data. The adopted brightness thresholds select from the EIT 195 \AA\ and AIA 193 \AA\ image practically the same areas of arcades and dimmings with a cross-calibration factor of 3.6-5.8 (5.0-8.2) for the AIA exposure time of 2.0 s (2.9 s). For the same photospheric areas, the MDI magnetic flux systematically exceeds the HMI flux by a factor of 1.4. Based on these results, the empirical diagnostic relationships obtained from SOHO data are adjusted to SDO instruments. Examples of a post-diagnostics based on SDO data are presented. As before, the tool is applicable to non-recurrent GMSs and FDs caused by nearly central eruptions from active regions, provided that the southern component of the interplanetary magnetic field near the Earth is predominantly negative, which is not predicted by this tool.


title: Virtual photon polarization in ultrarelativistic heavy-ion collisions

abstract: The polarization of direct photons produced in an ultrarelativistic heavy-ion collision reflects the momentum anisotropy of the quark-gluon plasma created in the collision. This paper presents a general framework, based on the photon spectral functions in the plasma, for analyzing the angular distribution and thus the polarization of dileptons in terms of the plasma momentum anisotropies. The rates of dilepton production depend, in general, on four independent spectral functions, corresponding to two transverse polarizations, one longitudinal polarization, and -- in plasmas in which the momentum anisotropy is not invariant under parity in the local rest frame of the matter -- a new spectral function, $\rho_n$, related to the anisotropy direction in the collision. The momentum anisotropy appears in the difference of the two transverse spectral functions, as well as in $\rho_n$. As an illustration, we delineate the spectral functions for dilepton pairs produced in the lowest order Drell-Yan process of quark-antiquark annihilation to a virtual photon.


title: Growth of two-dimensional Au patches in graphene pores: a density-functional study

abstract: Inspired by recent studies of various two-dimensional (2D) metals such as Au, Fe and Ag, we study the growth of two-dimensional gold patches in graphene pores by density-functional theory. We find that at room temperature gold atoms diffuse readily on top of both graphene and two-dimensional gold with energy barriers less than $0.5$ eV. Furthermore, gold atoms move without barriers from the top of graphene to its edge and from the top of 2D gold to its edge. The energy barriers are absent even at the interface of 2D gold and graphene, so that the gold atoms move effortlessly across the interface. We hope our demonstration for the propensity of diffusing gold atoms to grow 2D gold patches in graphene pores will inspire the fabrication of these patches experimentally.


title: An approximate global solution of Einstein's equations for differentially rotating compact body

abstract: We obtain an approximate global stationary and axisymmetric solution of Einstein's equations which can be thought as a simple star model: a self-gravitating perfect fluid ball with a differential rotation motion pattern. Using the post--Minkowskian formalism (weak field approximation) and considering rotation as a perturbation (slow rotation approximation), we find approximate interior and exterior (asymptotically flat) solutions to this problem in harmonic coordinates. Interior and exterior solutions are matched, in the sense of Lichnerowicz, on the surface of zero pressure to get a global solution. The resulting metric depends on four arbitrary constants: mass density, rotational velocity at $r=0$, a parameter which takes into account of the change in the rotational velocity through the star and the star radius in the non-rotation limit. The mass, angular momentum, quadrupole moment and other constants of the exterior metric are determined in terms of these four parameters.