title: SNR-adaptive OCT angiography enabled by statistical characterization of intensity and decorrelation with multi-variate time series model

abstract: In OCT angiography (OCTA), decorrelation computation has been widely used as a local motion index to identify dynamic flow from static tissues, but its dependence on SNR severely degrades the vascular visibility, particularly in low- SNR regions. To mathematically characterize the decorrelation-SNR dependence of OCT signals, we developed a multi-variate time series (MVTS) model. Based on the model, we derived a universal asymptotic linear relation of decorrelation to inverse SNR (iSNR), with the variance in static and noise regions determined by the average kernel size. Accordingly, with the population distribution of static and noise voxels being explicitly calculated in the iSNR and decorrelation (ID) space, a linear classifier is developed by removing static and noise voxels at all SNR, to generate a SNR-adaptive OCTA, termed as ID-OCTA. Then, flow phantom and human skin experiments were performed to validate the proposed ID-OCTA. Both qualitative and quantitative assessments demonstrated that ID-OCTA offers a superior visibility of blood vessels, particularly in the deep layer. Finally, implications of this work on both system design and hemodynamic quantification are further discussed.


title: Another approach to Hom-Lie bialgebras via Manin triples

abstract: In this paper, we study Hom-Lie bialgebras by a new notion of the dual representation of a representation of a Hom-Lie algebra. Motivated by the essential connection between Lie bialgebras and Manin triples, we introduce the notion of a Hom-Lie bialgebra with emphasis on its compatibility with a Manin triple of Hom-Lie algebras associated to a nondegenerate symmetric bilinear form satisfying a new invariance condition. With this notion, coboundary Hom-Lie bialgebras can be studied without a skew-symmetric condition of $r\in\mathfrak{g}\otimes \mathfrak{g}$, naturally leading to the classical Hom-Yang-Baxter equation whose solutions are used to construct coboundary Hom-Lie bialgebras. In particular, they are used to obtain a canonical Hom-Lie bialgebra structure on the double space of a Hom-Lie bialgebra. We also derive solutions of the classical Hom-Yang-Baxter equation from $\mathcal{O}$-operators and Hom-left-symmetric algebras.


title: PLMP -- Point-Line Minimal Problems in Complete Multi-View Visibility

abstract: We present a complete classification of all minimal problems for generic arrangements of points and lines completely observed by calibrated perspective cameras. We show that there are only 30 minimal problems in total, no problems exist for more than 6 cameras, for more than 5 points, and for more than 6 lines. We present a sequence of tests for detecting minimality starting with counting degrees of freedom and ending with full symbolic and numeric verification of representative examples. For all minimal problems discovered, we present their algebraic degrees, i.e. the number of solutions, which measure their intrinsic difficulty. It shows how exactly the difficulty of problems grows with the number of views. Importantly, several new minimal problems have small degrees that might be practical in image matching and 3D reconstruction.


title: Uniform dimension results for the inverse images of symmetric L\'evy processes

abstract: We prove uniform Hausdorff and packing dimension results for the inverse images of a large class of real-valued symmetric L\'evy processes. Our main result for the Hausdorff dimension extends that of Kaufman (1985) for Brownian motion and that of Song, Xiao, and Yang (2018) for $\alpha$-stable L\'evy processes with $1<\alpha<2$. Along the way, we also prove an upper bound for the uniform modulus of continuity of the local times of these processes.


title: Bipartite Graphs as Polynomials, and Polynomials as Bipartite Graphs (with a view towards dividing in $\mathbb{N}[x],$ $\mathbb{N}[x,y]$)

abstract: The aim of this paper is to show that any finite undirected bipartite graph can be considered as a polynomial $p \in \mathbb{N}[x]$, and any directed finite bipartite graph can be considered as a polynomial $p\in\mathbb{N}[x,y]$, and vise verse. We also show that the multiplication in semirings $\mathbb{N}[x]$, $\mathbb{N}[x,y]$ correspondences to a operations of the corresponding graphs which looks like a ``perturbed'' products of graphs. As an application, we give a new point of view to dividing in semirings $\mathbb{N}[x]$, $\mathbb{N}[x,y]$. Finally, we endow the set of all bipartite graphs with the Zariski topology.


title: Optical spectroscopic observations of gamma-ray blazar candidates VIII: the 2016-2017 follow up campaign carried out at SPM, NOT, KPNO and SOAR telescopes

abstract: The third Fermi source catalog lists 3033 gamma-ray sources above $4\sigma$ significance. More than 30% are classified as either unidentified/unassociated Gamma-ray sources (UGSs), with about 20% classified as Blazar candidates of uncertain types (BCUs). To confirm the blazar-like nature of candidate counterparts of UGSs and BCUs, we started in 2012 an optical spectroscopic follow up campaign. We report here the spectra of 36 targets with observations from the Observatorio Astron\'omico Nacional San Pedro M\'artir, the Southern Astrophysical Research Observatory, the Kitt Peak National Observatory and the Northern Optical Telescope, between 2016 and 2017. We confirm the BL Lac nature of 23 sources, and the flat spectrum radio quasar nature of other 7 ones. We also provide redshift estimates for 19 out of these 30 confirmations, with only one being a lower limit due to spectral features ascribable to intervening systems along the line of sight. As in previous analyses, the largest fraction of now-classified BCUs belong to the class of BL Lac objects, that appear to be the most elusive class of active galactic nuclei. One of the BL Lacs identified in this work, associated with 3FGL J2213.6-4755, lies at a redshift of $z>$1.529, making it one of the few distant gamma-ray BL Lac objects.


title: A mixture of experts model for predicting persistent weather patterns

abstract: Weather and atmospheric patterns are often persistent. The simplest weather forecasting method is the so-called persistence model, which assumes that the future state of a system will be similar (or equal) to the present state. Machine learning (ML) models are widely used in different weather forecasting applications, but they need to be compared to the persistence model to analyse whether they provide a competitive solution to the problem at hand. In this paper, we devise a new model for predicting low-visibility in airports using the concepts of mixture of experts. Visibility level is coded as two different ordered categorical variables: cloud height and runway visual height. The underlying system in this application is stagnant approximately in 90% of the cases, and standard ML models fail to improve on the performance of the persistence model. Because of this, instead of trying to simply beat the persistence model using ML, we use this persistence as a baseline and learn an ordinal neural network model that refines its results by focusing on learning weather fluctuations. The results show that the proposal outperforms persistence and other ordinal autoregressive models, especially for longer time horizon predictions and for the runway visual height variable.


title: The Encoding and Decoding Complexities of Entanglement-Assisted Quantum Stabilizer Codes

abstract: Quantum error-correcting codes are used to protect quantum information from decoherence. A raw state is mapped, by an encoding circuit, to a codeword so that the most likely quantum errors from a noisy quantum channel can be removed after a decoding process. A good encoding circuit should have some desired features, such as low depth, few gates, and so on. In this paper, we show how to practically implement an encoding circuit of gate complexity $O(n(n-k+c)/\log n)$ for an $[[n,k;c]]$ quantum stabilizer code with the help of $c$ pairs of maximally-entangled states. For the special case of an $[[n,k]]$ stabilizer code with $c=0$, the encoding complexity is $O(n(n-k)/\log n)$, which is previously known to be $O(n^2/\log n)$. For $c>0,$ this suggests that the benefits from shared entanglement come at an additional cost of encoding complexity. Finally we discuss decoding of entanglement-assisted quantum stabilizer codes and extend previously known computational hardness results on decoding quantum stabilizer codes.


title: Optical Spectroscopic Observations of Gamma-Ray Blazar Candidates. VII. Follow-up Campaign in the Southern Hemisphere

abstract: Searching for low energy counterparts of gamma-ray sources is one of the major challenges in modern gamma-ray astronomy. In the third Fermi source catalog about 30 % of detected sources are unidentified/unassociated Gamma-ray Sources (UGSs). We recently started an optical spectroscopic follow up campaign to confirm the blazar-like nature of candidates counterparts of UGSs. Here we report the spectra of 61 targets collected with the Southern Astrophysical Research Telescope (SOAR) between 2014 and the 2017. Our sample includes 33 potential counterparts of UGSs, selected on the basis of WISE colors, and 27 blazar candidates of uncertain type associated with gamma-ray sources of the last release of the Fermi catalog. We confirm the BZB nature of 20 sources lying within the positional uncertainty region of the UGSs. All the observed BCUs show blazar-like spectra, classified as 2 BZQs and 25 BZBs, for which we obtained 6 redshift estimates. Within the BCUs observations we report the redshift estimate for the BZB associated with, 3FGL J1106.4-3643 that is the second most distant BL Lac known to date, at z>1.084.


title: The Sloan Digital Sky Survey Reverberation Mapping Project: Low-Ionization Broad-Line Widths and Implications for Virial Black Hole Mass Estimation

abstract: The width of the broad emission lines in quasars is commonly characterized either by the full-width-at-half-maximum (FWHM) or the square root of the second moment of the line profile ($\sigma_{\rm line}$), and used as an indicator of the virial velocity of the broad-line region (BLR) in the estimation of black hole (BH) mass. We measure FWHM and $\sigma_{\rm line}$ for H$\alpha$, H$\beta$ and Mg II broad lines in both the mean and root-mean-square (rms) spectra of a large sample of quasars from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project. We introduce a new quantitative recipe to measure $\sigma_{\rm line}$ that is reproducible, less susceptible to noise and blending in the wings, and scales with the intrinsic width of the line. We compare the four definitions of line width (FWHM and $\sigma_{\rm line}$ in mean and rms spectra, respectively) for each of the three broad lines and among different lines. There are strong correlations among different width definitions for each line, providing justification for using the line width measured in single-epoch spectroscopy as a virial velocity indicator. There are also strong correlations among different lines, suggesting alternative lines to H$\beta$ can be used to estimate virial BH masses. We further investigate the correlations between virial BH masses using different line width definitions and the stellar velocity dispersion of the host galaxies, and the dependence of line shape (characterized by the ratio FWHM/$\sigma_{\rm line}$) on physical properties of the quasar. Our results provide further evidence that FWHM is more sensitive to the orientation of a flattened BLR geometry than $\sigma_{\rm line}$, but the overall comparison between the virial BH mass and host stellar velocity dispersion does not provide conclusive evidence that one particular width definition is significantly better than the others.


title: Kinematics of disk galaxies in (proto-)clusters at z=1.5

abstract: We observed star-forming galaxies at z~1.5 selected from the HyperSuprimeCam Subaru Strategic Program. The galaxies are part of two significant overdensities of [OII] emitters identified via narrow-band imaging and photometric redshifts from grizy photometry. We used VLT/KMOS to carry out Halpha integral field spectroscopy of 46 galaxies in total. Ionized gas maps, star formation rates and velocity fields were derived from the Halpha emission line. We quantified morphological and kinematical asymmetries to test for potential gravitational (e.g. galaxy-galaxy) or hydrodynamical (e.g. ram-pressure) interactions. Halpha emission was detected in 36 targets. 34 of the galaxies are members of two (proto-)clusters at z=1.47, confirming our selection strategy to be highly efficient. By fitting model velocity fields to the observed ones, we determined the intrinsic maximum rotation velocity Vmax of 14 galaxies. Utilizing the luminosity-velocity (Tully-Fisher) relation, we find that these galaxies are more luminous than their local counterparts of similar mass by up to ~4 mag in the rest-frame B-band. In contrast to field galaxies at z<1, the offsets of the z~1.5 (proto-)cluster galaxies from the local Tully-Fisher relation are not correlated with their star formation rates but with the ratio between Vmax and gas velocity dispersion sigma_g. This probably reflects that, as is observed in the field at similar redshifts, fewer disks have settled to purely rotational kinematics and high Vmax/sigma_g ratios. Due to relatively low galaxy velocity dispersions (sigma_v < 400 km/s) of the (proto-)clusters, gravitational interactions likely are more efficient, resulting in higher kinematical asymmetries, than in present-day clusters. (abbr.)


title: Ground-Based Radial Velocity as Critical Support for Future NASA Earth-Finding Missions

abstract: Future space-based direct imaging missions are poised to search for biosignatures in the atmospheres of potentially habitable planets orbiting nearby AFGKM stars. Although these missions could conduct a survey of high-priority target stars to detect candidate Earth-like planets, conducting a precursor radial velocity (RV) survey will benefit future direct imaging missions in four ways. First, an RV survey capable of detecting signals as small as 8 cm/s over timescales of a few years could discover potentially habitable Earth-mass planets orbiting dozens of nearby GKM stars accessible to space-based direct imaging. Second, RVs will improve scheduling efficiency by reducing the required number of revisits for orbit determination, and revealing when a planet of interest is most observable. Third, RV observations will reveal the masses of gas and ice giants that could be mistaken for Earth-mass planets, thereby reducing the time spent identifying false positives. Fourth, mass measurements from RVs will provide the surface gravities necessary for interpreting atmospheric spectra and potential biosignatures.


title: The Meeting of Acquaintances: A Cost-efficient Authentication Scheme for Light-weight Objects with Transient Trust Level and Plurality Approach

abstract: Wireless sensor networks consist of a large number of distributed sensor nodes so that potential risks are becoming more and more unpredictable. The new entrants pose the potential risks when they move into the secure zone. To build a door wall that provides safe and secured for the system, many recent research works applied the initial authentication process. However, the majority of the previous articles only focused on the Central Authority (CA) since this leads to an increase in the computation cost and energy consumption for the specific cases on the Internet of Things (IoT). Hence, in this article, we will lessen the importance of these third parties through proposing an enhanced authentication mechanism that includes key management and evaluation based on the past interactions to assist the objects joining a secured area without any nearby CA. We refer to a mobility dataset from CRAWDAD collected at the University Politehnica of Bucharest and rebuild into a new random dataset larger than the old one. The new one is an input for a simulated authenticating algorithm to observe the communication cost and resource usage of devices. Our proposal helps the authenticating flexible, being strict with unknown devices into the secured zone. The threshold of maximum friends can modify based on the optimization of the symmetric-key algorithm to diminish communication costs (our experimental results compare to previous schemes less than 2000 bits) and raise flexibility in resource-constrained environments.


title: Dynamic Spatial Verification for Large-Scale Object-Level Image Retrieval

abstract: Images from social media can reflect diverse viewpoints, heated arguments, and expressions of creativity, adding new complexity to retrieval tasks. Researchers working onContent-Based Image Retrieval (CBIR) have traditionally tuned their algorithms to match filtered results with user search intent. However, we are now bombarded with composite images of unknown origin, authenticity, and even meaning. With such uncertainty, users may not have an initial idea of what the results of a search query should look like. For instance, hidden people, spliced objects, and subtly altered scenes can be difficult for a user to detect initially in a meme image, but may contribute significantly to its composition. We propose a new approach for spatial verification that aims at modeling object-level regions dynamically clustering keypoints in a 2D Hough space, which are then used to accurately weight small contributing objects within the results, without the need for costly object detection steps. We call this method Objects in Scene to Objects in Scene (OS2OS) score, and it is optimized for fast matrix operations on CPUs. OS2OS performs comparably to state-of-the-art methods in classic CBIR problems, on the Oxford5K, Paris 6K, and Google-Landmarks datasets, without the need for bounding boxes. It also succeeds in emerging retrieval tasks such as image composite matching in the NIST MFC2018 dataset and meme-style composite imagery fromReddit.


title: Self-similar spreading in a merging-splitting model of animal group size

abstract: In a recent study of certain merging-splitting models of animal-group size (Degond et al., J. Nonl. Sci. 27 (2017) 379), it was shown that an initial size distribution with infinite first moment leads to convergence to zero in weak sense, corresponding to unbounded growth of group size. In the present paper we show that for any such initial distribution with a power-law tail, the solution approaches a self-similar spreading form. A one-parameter family of such self-similar solutions exists, with densities that are completely monotone, having power-law behavior in both small and large size regimes, with different exponents.


title: Effect of mounting strut and cavitator shape on the ventilation demand for ventilated supercavitation

abstract: The present work reports behaviors regarding the formation and collapse of a ventilated supercavity while varying the cavitator shapes, including triangle, disk, and cone and varying mounting struts. Three cavitators with the same frontal area are fabricated with 3D printing and mounted on a forward facing model (FFM). The ventilation requirements to generate (C_Qf) and sustain (C_Qc) a supercavity are tested over a wide range of Froude number (Fr) for each cavitator and compared with backward facing model(BFM). Compared to the triangle and disk cavitators, the cone-shaped cavitator requires the least amount of air to generate a supercavity in nearly all of the tested flow regime except very high Fr. The C_Qc of disk FFM is lower than that of its BFM counterpart at small Fr and exceeds the BFM C_Qc with further increase of Fr. The cone cavitator has the smallest C_Qc among all the cavitators across the range of Fr in our experiments. Simultaneous internal pressure and cavity dimension measurements are conducted to elucidate the cavity sustenance behaviors. The cone-generated cavity yields a significantly smaller maximum diameter and a shorter half-length. Cavity geometric information and cavity pressure measurements with high-speed imaging of re-entrant jet are employed to estimate the re-entrant jet momentum under different Fr for disk and cone cavitators. The estimated re-entrant jet momentum shows reasonable match with the ventilation air momentum under C_Qc in lower Fr for both cavitator cases, with the disk cavitator case yielding significantly stronger re-entrant jet, providing support to the re-entrant jet mechanism governing on the cavity collapse. Our study sheds some light on the cavitator design and ventilation strategy for a supercavitating vehicle in practice.


title: Exploiting Synthetically Generated Data with Semi-Supervised Learning for Small and Imbalanced Datasets

abstract: Data augmentation is rapidly gaining attention in machine learning. Synthetic data can be generated by simple transformations or through the data distribution. In the latter case, the main challenge is to estimate the label associated to new synthetic patterns. This paper studies the effect of generating synthetic data by convex combination of patterns and the use of these as unsupervised information in a semi-supervised learning framework with support vector machines, avoiding thus the need to label synthetic examples. We perform experiments on a total of 53 binary classification datasets. Our results show that this type of data over-sampling supports the well-known cluster assumption in semi-supervised learning, showing outstanding results for small high-dimensional datasets and imbalanced learning problems.


title: Adaptive wavefront correction of dynamic multimode beam based on modal decomposition

abstract: We propose and demonstrate a method for the adaptive wavefront correction of dynamic multimode fiber beams for the first time. The wavefront of incident beam is reconstructed in real-time based on the complete modal information, which obtained by using the modal decomposition of correlation filter method. For the proof of principle, both of the modal decomposition and the wavefront correction are implemented using the same computer-generated hologram, which encoded into a phase-only spatial light modulator. We demonstrate the wavefront correction of dynamic multimode beam at a rate of 5Hz and achieve a 1.73-fold improvement on the average power-in-the-bucket. The experimental results indicate the feasibility of the real-time wavefront correction for the large mode area fiber laser by adaptive optics.


title: Nucleosynthesis in heavy-ion collisions at the LHC via the Saha equation

abstract: The production of light (anti-)(hyper-)nuclei in heavy-ion collisions at the LHC is considered in the framework of the Saha equation, making use of the analogy between the evolution of the early universe after the Big Bang and that of "Little Bangs" created in the lab. Assuming that disintegration and regeneration reactions involving light nuclei proceed in relative chemical equilibrium after the chemical freeze-out of hadrons, their abundances are determined through the famous cosmological Saha equation of primordial nucleosynthesis and show no exponential dependence on the temperature typical for the thermal model. A quantitative analysis, performed using the hadron resonance gas model in partial chemical equilibrium, shows agreement with experimental data of the ALICE collaboration on d, $^3$He, $^3_\Lambda$H, and $^4$He yields for a very broad range of temperatures at $T \lesssim 155$ MeV. The presented picture is supported by the observed suppression of resonance yields in central Pb-Pb collisions at the LHC.


title: Generalization of k-means Related Algorithms

abstract: This article briefly introduced Arthur and Vassilvitshii's work on \textbf{k-means++} algorithm and further generalized the center initialization process. It is found that choosing the most distant sample point from the nearest center as new center can mostly have the same effect as the center initialization process in the \textbf{k-means++} algorithm.


title: One Proximate Kitaev Spin Liquid in the $K$-$J$-$\Gamma$ Model on the Honeycomb Lattice

abstract: In addition to the Kitaev ($K$) interaction, candidate Kitaev materials also possess Heisenberg ($J$) and off-diagonal symmetric ($\Gamma$) couplings. We investigate the quantum ($S = 1/2$) $K$-$J$-$\Gamma$ model on the honeycomb lattice by a variational Monte Carlo (VMC) method. In addition to the "generic" Kitaev spin liquid (KSL), we find that there is just one proximate KSL (PKSL) phase, while the rest of the phase diagram contains different magnetically ordered states. The PKSL is a gapless Z$_2$ state with 14 Majorana cones, which in contrast to the KSL has a gapless spin response. In a magnetic field applied normal to the honeycomb plane, it realizes two of Kitaev's gapped chiral spin-liquid phases, of which one is non-Abelian with Chern number $\nu = 5$ and the other is Abelian with $\nu = 4$. These two phases could be distinguished by their thermal Hall conductance.


title: Core Access Hybridization in 5G

abstract: The fifth generation (5G) cellular communication system will have several advanced features, majority of which will come into existence for the first time. According to the ITU recommendations, 5G will provide minimum 3 GB/s end to end data rates under the static conditions and in the mobile condition the minimum data rate is 100 Mb/s. The latency will be brought down to 1 ms and the device densities will increase by many folds. These features cannot be provided in the pure wireless domain. Certainly, the support of the optical networks is essential to provide these advanced features in 5G communication. These features need special architectures. In this article, we provide advanced optical wireless hybrid architectures for 5G networks. We explain the importance of this architecture and also show the nodal configuration of these hybrid networks.


title: Navigated interventions in the head and neck area: standardized assessment of a new handy field generator

abstract: Electromagnetic (EM) tracking enables localization of surgical instruments within the magnetic field emitted by an EM field generator (FG). Usually, the larger a FG is, the larger its tracking volume is. However, the company NDI (Northern Digital Inc., Waterloo, ON, Canada) recently introduced the Planar 10-11 FG, which combines a compact construction (97mm x 112mm x 31mm) with a relatively large, cylindrical tracking volume (diameter: 340mm, height: 340mm). Using the standardized assessment protocol of Hummel et al., the FG was tested with regard to its tracking accuracy and to its robustness with respect to external sources of disturbance. The mean positional error (5cm distance metric according to Hummel protocol) was 0.59mm, with a mean jitter of 0.26mm in the standard setup. The mean orientational error was found to be 0.10{\deg}. The highest positional error (4.82mm) due to metallic sources of disturbance was caused by the steel SST 303. In contrast, steel SST 416 caused the lowest positional error (0.10mm). Overall, the Planar 10-11 FG tends to achieve better tracking accuracy results compared to other NDI FGs. Due to its compact construction and portability, the FG could contribute to increased clinical use of EM tracking systems.


title: Enhancement of superluminal weak values under Lorentz boost

abstract: We study the local group velocity defined as the weak value of the velocity operator in the (1+1) dimensional Klein-Gordon as well as Dirac theory. It was shown by Berry [ J. Phys. A 45, 185308 (2012)] that when the pre- and post-selected states for evaluating the weak value are chosen at random from an ensemble of available states, the local group velocity has a universal probability distribution which can have both subluminal and superluminal components. In this work, we numerically explore the role of Lorentz boost and its impact on the superluminal fraction of the total probability distribution. We show that the dependence (enhancement) of the superluminal fraction on Lorentz boost of the total probability distribution differs both qualitatively and quantitatively for the Klein-Gordon waves and Dirac waves. For the Klein-Gordan waves, the asymmetry in the distribution of group velocities around the zero velocity point in the laboratory frame is entirely responsible for the observation of relative enhancement in the boosted frame. On the other hand, for the Dirac waves, we observe an enhancement irrespective of whether the laboratory frame velocity distribution is symmetric or not.


title: On symmetries of generalized Calogero model and Polychronakos-Frahm chain

abstract: The symmetry of the generalized Polychronakos-Frahm chain is obtained from the Dunkl-operator deformation of the unitary algebra, which describes the symmetry of the generalized Calogero model.

categories: hep-th cond-mat.str-el math-ph math.MP nlin.SI


title: Panchromatic patterns by paths

abstract: Let $H=(V_H,A_H)$ be a digraph, possibly with loops, and let $D=(V_D, A_D)$ be a loopless multidigraph with a colouring of its arcs $c: A_D \rightarrow V_H$. An $H$-path of $D$ is a path $(v_0, \dots, v_n)$ of $D$ such that $(c(v_{i-1}, v_i), c(v_i,v_{i+1}))$ is an arc of $H$ for every $1 \le i \le n-1$. For $u, v \in V_D$, we say that $u$ reaches $v$ by $H$-paths if there exists an $H$-path from $u$ to $v$ in $D$. A subset $S \subseteq V_D$ is $H$-absorbent of $D$ if every vertex in $V_D-S$ reaches by $H$-paths some vertex in $S$, and it is $H$-independent if no vertex in $S$ can reach another (different) vertex in $S$ by $H$-pahts. An $H$-kernel is an independent by $H$-paths and absorbent by $H$-paths subset of $V_D$. We define $\tilde{\mathscr{B}}_1$ as the set of digraphs $H$ such that any $H$-arc-coloured tournament has an $H$-absorbent by paths vertex; the set $\tilde{\mathscr{B}}_2$ consists of the digraphs $H$ such that any $H$-arc-coloured digraph $D$ has an independent, $H$-absorbent by paths set; analogously, the set $\tilde{\mathscr{B}}_3$ is the set of digraphs $H$ such that every $H$-arc-coloured digraph $D$ contains an $H$-kernel by paths. In this work, we present a characterization of $\tilde{\mathscr{B}}_2$, and provide structural properties of the digraphs in $\tilde{\mathscr{B}}_3$ which settle up its characterization except for the analysis of a single digraph on three vertices.


title: A Fast Particle-Based Approach for Calibrating a 3-D Model of the Antarctic Ice Sheet

abstract: We consider the scientifically challenging and policy-relevant task of understanding the past and projecting the future dynamics of the Antarctic ice sheet. The Antarctic ice sheet has shown a highly nonlinear threshold response to past climate forcings. Triggering such a threshold response through anthropogenic greenhouse gas emissions would drive drastic and potentially fast sea level rise with important implications for coastal flood risks. Previous studies have combined information from ice sheet models and observations to calibrate model parameters. These studies have broken important new ground but have either adopted simple ice sheet models or have limited the number of parameters to allow for the use of more complex models. These limitations are largely due to the computational challenges posed by calibration as models become more computationally intensive or when the number of parameters increases. Here we propose a method to alleviate this problem: a fast sequential Monte Carlo method that takes advantage of the massive parallelization afforded by modern high performance computing systems. We use simulated examples to demonstrate how our sample-based approach provides accurate approximations to the posterior distributions of the calibrated parameters. The drastic reduction in computational times enables us to provide new insights into important scientific questions, for example, the impact of Pliocene era data and prior parameter information on sea level projections. These studies would be computationally prohibitive with other computational approaches for calibration such as Markov chain Monte Carlo or emulation-based methods. We also find considerable differences in the distributions of sea level projections when we account for a larger number of uncertain parameters.


title: A Formalization of Robustness for Deep Neural Networks

abstract: Deep neural networks have been shown to lack robustness to small input perturbations. The process of generating the perturbations that expose the lack of robustness of neural networks is known as adversarial input generation. This process depends on the goals and capabilities of the adversary, In this paper, we propose a unifying formalization of the adversarial input generation process from a formal methods perspective. We provide a definition of robustness that is general enough to capture different formulations. The expressiveness of our formalization is shown by modeling and comparing a variety of adversarial attack techniques.


title: What is the spectral category?

abstract: For a category $\mathcal{C}$ with finite limits and a class $\mathcal{S}$ of monomorphisms in $\mathcal{C}$ that is pullback stable, contains all isomorphisms, is closed under composition, and has the strong left cancellation property, we use pullback stable $\mathcal{S}$-essential monomorphisms in $\mathcal{C}$ to construct a spectral category $\mathrm{Spec}(\mathcal{C},\mathcal{S})$. We show that it has finite limits and that the canonical functor $\mathcal{C}\to \mathrm{Spec}(\mathcal{C},\mathcal{S})$ preserves finite limits. When $\mathcal{C}$ is a normal category, assuming for simplicity that $\mathcal{S}$ is the class of all monomorphisms in $\mathcal{C}$, we show that pullback stable $\mathcal{S}$-essential monomorphisms are the same as what we call subobject-essential monomorphisms.


title: Automated Classification of Histopathology Images Using Transfer Learning

abstract: There is a strong need for automated systems to improve diagnostic quality and reduce the analysis time in histopathology image processing. Automated detection and classification of pathological tissue characteristics with computer-aided diagnostic systems are a critical step in the early diagnosis and treatment of diseases. Once a pathology image is scanned by a microscope and loaded onto a computer, it can be used for automated detection and classification of diseases. In this study, the DenseNet-161 and ResNet-50 pre-trained CNN models have been used to classify digital histopathology patches into the corresponding whole slide images via transfer learning technique. The proposed pre-trained models were tested on grayscale and color histopathology images. The DenseNet-161 pre-trained model achieved a classification accuracy of 97.89% using grayscale images and the ResNet-50 model obtained the accuracy of 98.87% for color images. The proposed pre-trained models outperform state-of-the-art methods in all performance metrics to classify digital pathology patches into 24 categories.


title: Lock-Free Transactional Adjacency List

abstract: Adjacency lists are frequently used in graphing or map based applications. Although efficient concurrent linked-list algorithms are well known, it can be difficult to adapt these approaches to build a high-performance adjacency list. Furthermore, it can often be desirable to execute operations in these data structures transactionally, or perform a sequence of operations in one atomic step. In this paper, we present a lock-free transactional adjacency list based on a multi-dimensional list (MDList). We are able to combine known linked list strategies with the capability of the MDList in order to efficiently organize graph vertexes and their edges. We design our underlying data structure to be node-based and linearizable, then use the Lock-Free Transactional Transformation (LFTT) methodology to efficiently enable transactional execution. In our performance evaluation, our lock-free transactional adjacency list achieves an average of 50% speedup over a transactional boosting implementation.


title: Representations of the Vertex Reinforced Jump Process as a mixture of Markov processes on $\mathbb{Z}^d$ and infinite trees

abstract: This paper concerns the Vertex Reinforced Jump Process (VRJP) and its representations as a Markov process in random environment. We show that all possible representations of the VRJP as a mixture of Markov processes can be expressed in a similar form, using a random potential and harmonic functions for an associated operator. This allows to show that the VRJP on $\mathbb{Z}^d$ (with certain initial conditions) has a unique representation, by proving that an associated Martin boundary is trivial. Moreover, on infinite trees, we construct a family of representations, that are all different when the VRJP is transient and the tree is $d$-regular (with $d\geq 3$).


title: Consequences of Modified Cosmologies in DM abundance and PeV IceCube signals

abstract: To explain the high-energy astrophysical neutrino flux with energies $\sim PeV$ recently observed by IceCube collaboration and the relic abundance of DM in the Universe, we consider the minimal 4-dimensional operator $\sim y_{\alpha \chi}\overline{{L_{L_{\alpha}}}}\, H\, \chi$. The cosmological background is assumed to evolve according to $f({\cal T})$ cosmology, where ${\cal T}$ is the scalar torsion.


title: A Strongly Consistent Sparse $k$-means Clustering with Direct $l_1$ Penalization on Variable Weights

abstract: We propose the Lasso Weighted $k$-means ($LW$-$k$-means) algorithm as a simple yet efficient sparse clustering procedure for high-dimensional data where the number of features ($p$) can be much larger compared to the number of observations ($n$). In the $LW$-$k$-means algorithm, we introduce a lasso-based penalty term, directly on the feature weights to incorporate feature selection in the framework of sparse clustering. $LW$-$k$-means does not make any distributional assumption of the given dataset and thus, induces a non-parametric method for feature selection. We also analytically investigate the convergence of the underlying optimization procedure in $LW$-$k$-means and establish the strong consistency of our algorithm. $LW$-$k$-means is tested on several real-life and synthetic datasets and through detailed experimental analysis, we find that the performance of the method is highly competitive against some state-of-the-art procedures for clustering and feature selection, not only in terms of clustering accuracy but also with respect to computational time.


title: Safe Learning-Based Control of Stochastic Jump Linear Systems: a Distributionally Robust Approach

abstract: We consider the problem of designing control laws for stochastic jump linear systems where the disturbances are drawn randomly from a finite sample space according to an unknown distribution, which is estimated from a finite sample of i.i.d. observations. We adopt a distributionally robust approach to compute a mean-square stabilizing feedback gain with a given probability. The larger the sample size, the less conservative the controller, yet our methodology gives stability guarantees with high probability, for any number of samples. Using tools from statistical learning theory, we estimate confidence regions for the unknown probability distributions (ambiguity sets) which have the shape of total variation balls centered around the empirical distribution. We use these confidence regions in the design of appropriate distributionally robust controllers and show that the associated stability conditions can be cast as a tractable linear matrix inequality (LMI) by using conjugate duality. The resulting design procedure scales gracefully with the size of the probability space and the system dimensions. Through a numerical example, we illustrate the superior sample complexity of the proposed methodology over the stochastic approach.


title: Parallel ADMM for robust quadratic optimal resource allocation problems

abstract: An alternating direction method of multipliers (ADMM) solver is described for optimal resource allocation problems with separable convex quadratic costs and constraints and linear coupling constraints. We describe a parallel implementation of the solver on a graphics processing unit (GPU) using a bespoke quartic function minimizer. An application to robust optimal energy management in hybrid electric vehicles is described, and the results of numerical simulations comparing the computation times of the parallel GPU implementation with those of an equivalent serial implementation are presented.


title: BpForms and BcForms: Tools for concretely describing non-canonical polymers and complexes to facilitate comprehensive biochemical networks

abstract: Although non-canonical residues, caps, crosslinks, and nicks play an important role in the function of many DNA, RNA, proteins, and complexes, we do not fully understand how networks of non-canonical macromolecules generate behavior. One barrier is our limited formats, such as IUPAC, for abstractly describing macromolecules. To overcome this barrier, we developed BpForms and BcForms, a toolkit of ontologies, grammars, and software for abstracting the primary structure of polymers and complexes as combinations of residues, caps, crosslinks, and nicks. The toolkit can help quality control, exchange, and integrate information about the primary structure of macromolecules into fine-grained global networks of intracellular biochemistry.


title: Directed immersions for complex structures

abstract: We analyze the differential relation corresponding to integrability of almost complex structures, reformulated as a directed immersion relation by Demailly and Gaussier. Combining results of Clemente [3], we show that applying h-principle techniques yields the following statement: for an almost complex manifold with arbitrary metric $(X, J, g)$, and for $\epsilon > 0$, there exists a smooth function $f : X \rightarrow \mathbb{R}$ and almost complex structure $J'$ on $X$ such that $J$ and $J'$ are $C^0$-close on the graph of $f$ with respect to the extended metric on $X \times \mathbb{R}$, and such that the Nijenhuis tensor of $J'$ on the graph has pointwise sup norm less than $C\epsilon$, where $C$ is a constant depending only on $J$ and $g$. This is an updated version of a previous preprint titled "Almost complex manifolds are (almost) complex".


title: Spectral analysis of mixing in 2D high-Reynolds flows

abstract: We use spectral analysis of Eulerian and Lagrangian dynamics to study the advective mixing in an incompressible 2D bounded cavity flow. A significant property of such a rotational flow at high Reynolds numbers is that mixing in its core is slower than wall-adjacent areas and corner eddies. We explain this property by considering the resonance between frequencies of unsteady motion -- captured by the Koopman spectral analysis of the velocity field -- and the circulation frequency of Lagrangian tracers in the mean flow. In high-Reynolds rotational 2D flows, the vorticity in the rotational core is uniformly distributed, which leads to uniform distribution of circulation periods in the mean flow, i.e., the kinematics in the core of mean flow is like rigid-body rotation. When this ``rigid" core is exposed to velocity fluctuations arising from bifurcations at high Reynolds, it shows more resilience toward resonance in Lagrangian motion and hence mixes more slowly compared to other areas. We also show how our qualitative resonance argument extends to chaotic flows where the classical tools of dynamical systems are not applicable.


title: Quantum critical behaviour and Lifshitz transition in intermediate valence $\alpha$-YbAlB$_4$

abstract: Intermetallic compounds containing $f$-electron elements have been prototypical materials for investigating strong electron correlations and quantum criticality (QC). Their heavy fermion ground state evoked by the magnetic $f$-electrons is susceptible to the onset of quantum phases, such as magnetism or superconductivity, due to the enhanced effective mass ($m^{*}$) and a corresponding decrease of the Fermi temperature. However, the presence of $f$-electron valence fluctuations to a non-magnetic state is regarded an anathema to QC, as it usually generates a paramagnetic Fermi-liquid state with quasiparticles of moderate $m^{*}$. Such systems are typically isotropic, with a characteristic energy scale $T_0$ of the order of hundreds of kelvins that require large magnetic fields or pressures to promote a valence or magnetic instability. Here we show that the intermediate valence compound $\alpha$-YbAlB$_4$ surprisingly exhibits both quantum critical behaviour and a Lifshitz transition under low magnetic field, which is attributed to the anisotropy of the hybridization between the conduction and localized $f$-electrons. These findings suggest a new route to bypass the large valence energy scale in developing the QC.


title: Downlink Training in Cell-Free Massive MIMO: A Blessing in Disguise

abstract: Cell-free Massive MIMO (multiple-input multiple-output) refers to a distributed Massive MIMO system where all the access points (APs) cooperate to coherently serve all the user equipments (UEs), suppress inter-cell interference and mitigate the multiuser interference. Recent works demonstrated that, unlike co-located Massive MIMO, the \textit{channel hardening} is, in general, less pronounced in cell-free Massive MIMO, thus there is much to benefit from estimating the downlink channel. In this study, we investigate the gain introduced by the downlink beamforming training, extending the previously proposed analysis to non-orthogonal uplink and downlink pilots. Assuming single-antenna APs, conjugate beamforming and independent Rayleigh fading channel, we derive a closed-form expression for the per-user achievable downlink rate that addresses channel estimation errors and pilot contamination both at the AP and UE side. The performance evaluation includes max-min fairness power control, greedy pilot assignment methods, and a comparison between achievable rates obtained from different capacity-bounding techniques. Numerical results show that downlink beamforming training, although increases pilot overhead and introduces additional pilot contamination, improves significantly the achievable downlink rate. Even for large number of APs, it is not fully efficient for the UE relying on the statistical channel state information for data decoding.


title: Approximation and Non-parametric Estimation of ResNet-type Convolutional Neural Networks

abstract: Convolutional neural networks (CNNs) have been shown to achieve optimal approximation and estimation error rates (in minimax sense) in several function classes. However, previous analyzed optimal CNNs are unrealistically wide and difficult to obtain via optimization due to sparse constraints in important function classes, including the H\"older class. We show a ResNet-type CNN can attain the minimax optimal error rates in these classes in more plausible situations -- it can be dense, and its width, channel size, and filter size are constant with respect to sample size. The key idea is that we can replicate the learning ability of Fully-connected neural networks (FNNs) by tailored CNNs, as long as the FNNs have \textit{block-sparse} structures. Our theory is general in a sense that we can automatically translate any approximation rate achieved by block-sparse FNNs into that by CNNs. As an application, we derive approximation and estimation error rates of the aformentioned type of CNNs for the Barron and H\"older classes with the same strategy.


title: Approximating local properties by tensor network states with constant bond dimension

abstract: Suppose we would like to approximate all local properties of a quantum many-body state to accuracy $\delta$. In one dimension, we prove that an area law for the Renyi entanglement entropy $R_\alpha$ with index $\alpha<1$ implies a matrix product state representation with bond dimension $\mathrm{poly}(1/\delta)$. For (at most constant-fold degenerate) ground states of one-dimensional gapped Hamiltonians, it suffices that the bond dimension is almost linear in $1/\delta$. In two dimensions, an area law for $R_\alpha(\alpha<1)$ implies a projected entangled pair state representation with bond dimension $e^{O(1/\delta)}$. In the presence of logarithmic corrections to the area law, similar results are obtained in both one and two dimensions.

categories: quant-ph cond-mat.str-el math-ph math.MP


title: Diagonal reduction of matrices over Bezout rings of stable range 1 with the Kazimirsky condition

abstract: We constuct the theory of diagonalizability for matrices over Bezout rings of stable range 1 with the Kazimirsky condition. It is shown that a ring of stable range 1 with the right (left) Kazimirsky condition is an elementary divisor ring if and only if it is a duo ring.


title: QED effects are negligible for neutron-star spin-down

abstract: The energy loss of a rotationally powered pulsar is primarily carried away as electromagnetic radiation and a particle wind. Considering that the magnetic field strength of pulsars ranges from about $10^8$ to $10^{15}$ G, one could expect quantum electrodynamics (QED) to play a role in their spin-down, especially for strongly magnetized ones (magnetars). In fact several authors have argued that QED corrections will dominate the spin-down for slowly rotating stars. They called this effect quantum vacuum friction (QVF). However, QVF was originally derived using a problematic self-torque technique, which leads to a dramatic overestimation of this spin-down effect. Here, instead of using QVF, we explicitly calculate the energy loss from rotating neutron stars using the Poynting vector and a model for a particle wind, and we include the QED one-loop corrections. We express the excess emission as QED one-loop corrections to the radiative magnetic moment of a neutron star. We do find a small component of the spin-down luminosity that originates from the vacuum polarization. However, it never exceeds one percent of the classical magnetic dipole radiation in neutron stars for all physically interesting field strengths. Therefore, we find that the radiative corrections of QED are irrelevant in the energetics of neutron-star spin-down.


title: Dynamics of shadow system of a singular Gierer-Meinhardt system on an evolving domain

abstract: The main purpose of the current paper is to contribute towards the comprehension of the dynamics of the shadow system of a singular Gierer-Meinhardt model on an isotropically evolving domain. In the case where the inhibitor's response to the activator's growth is rather weak, then the shadow system of the Gierer-Meinhardt model is reduced to a single though non-local equation whose dynamics is thoroughly investigated throughout the manuscript. The main focus is on the derivation of blow-up results for this non-local equation, which can be interpreted as instability patterns of the shadow system. In particular, a {\it diffusion-driven instability (DDI)}, or {\it Turing instability}, in the neighbourhood of a constant stationary solution, which then is destabilised via diffusion-driven blow-up, is observed. The latter indicates the formation of some unstable patterns, whilst some stability results of global-in-time solutions towards non-constant steady states guarantee the occurrence of some stable patterns. Most of the derived results are confirmed numerically and also compared with the ones in the case of a stationary domain.


title: Measurement of the production cross section for a Higgs boson in association with a vector boson in the $H \rightarrow WW^{\ast} \rightarrow \ell\nu\ell\nu$ channel in $pp$ collisions at $\sqrt{s}$ = 13 TeV with the ATLAS detector

abstract: A measurement of the Higgs boson production cross sections via associated $WH$ and $ZH$ production using $H \rightarrow WW^{\ast} \rightarrow \ell\nu\ell\nu$ decays, where $\ell$ stands for either an electron or a muon, is presented. Results for combined $WH$ and $ZH$ production are also presented. The analysis uses events produced in proton-proton collisions collected with the ATLAS detector at the Large Hadron Collider in 2015 and 2016. The data correspond to an integrated luminosity of 36.1 fb$ ^{-1}$ recorded at a centre-of-mass energy of 13 TeV. The products of the $H \rightarrow WW^{\ast}$ branching fraction times the $WH$ and $ZH$ cross sections are measured to be $0.67^{+0.31}_{-0.27}$(stat.)$^{+0.18}_{-0.14}$(syst.) pb and $0.54^{+0.31}_{-0.24}$(stat.)$^{+0.15}_{-0.07}$(syst.) pb respectively, in agreement with the Standard Model predictions.


title: Measuring light with light dependent resistors: an easy approach for optics experiments

abstract: We entertain the use of light dependent resistors as a viable option as measuring sensors in optics laboratory experiments or classroom demonstrations. The main advantages of theses devices are essentially very low cost, easy handling and commercial availability which can make them interesting for instructors with limited resources. Simple calibration procedures were developed indicating a precision of $\sim 5\% $ for illuminance measurements. Optical experiments were carried out as proof of feasibility for measurements of reflected and transmitted light and its quality results are presented. In particular, the sensor measurements allowed to verify the angular distribution of a Lambertian reflective material, to observe transmitted and reflected specular light on a glass slab as function of the incoming angle of a light beam, and to estimate glass refractive index with values averaging $1.51\pm0.06$ in satisfactory agreement with the expected 1.52 value.


title: Nonlocal statistical field theory of dipolar particles in electrolyte solutions

abstract: We present a nonlocal statistical field theory of a dilute electrolyte solution with small additive of dipolar particles. We postulate that every dipolar particle is associated with an arbitrary probability distribution function (PDF) of distance between its charge centers. Using the standard Hubbard-Stratonovich transformation, we represent the configuration integral of the system in the functional integral form. We show that in the limit of a small permanent dipole moment, the functional in integrand exponent takes the well known form of the Poisson-Boltzmann-Langevin (PBL) functional. In the mean-field approximation we obtain a non-linear integro-differential equation with respect to the mean-field electrostatic potential, generalizing the PBL equation for the point-like dipoles obtained first by Abrashkin et al. We apply the obtained equation in its linearized form to derivation of the expressions for the mean-field electrostatic potential of the point-like test ion and its solvation free energy in salt-free solution, as well as in solution with salt ions. For the 'Yukawa'-type PDF we obtain analytic relations for both the electrostatic potential and the solvation free energy of the point-like test ion. We obtain a general expression for the bulk electrostatic free energy of the solution within the random phase approximation (RPA). For the salt-free solution of the dipolar particles for the Yukawa-type PDF we obtain an analytic relation for the electrostatic free energy.


title: Nonlocal statistical field theory of dipolar particles forming chain-like clusters

abstract: We present a nonlocal statistical field theory of a diluted solution of dipolar particles which are capable of forming chain-like clusters in accordance with the 'head-to-tail' mechanism. As in our previous study [Yu.A. Budkov 2018 J. Phys.: Condens. Matter 30 344001], we model dipolar particles as dimers comprised of oppositely charged point-like groups, separated by fluctuating distance. For the special case of the Yukawa-type distribution function of distance between the charged groups of dipolar particles we obtain an analytical expression for the electrostatic free energy of solution within the random phase approximation. We show that an increase in the association constant leads to a decrease in the absolute value of the electrostatic free energy of solution, preventing its phase separation which is in agreement with the former computer simulations and theoretical results. We obtain a non-linear integro-differential equation for the self-consistent field potential created by the fixed external charges in a solution medium, taking into account the association of dipolar particles. As a consequence of the derived self-consistent field equation, in regime of weak electrostatic interactions, we obtain an analytical expression for the electrostatic potential of the point-like test ion, surrounded by the chain-like clusters of the dipolar particles. We show that in the mean-field approximation the association does not change the bulk dielectric permittivity of the solution, but increases the solvation radius of the point-like charge, relative to the theory of non-associating dipolar particles.


title: Invariant connections, Lie algebra actions, and foundations of numerical integration on manifolds

abstract: Motivated by numerical integration on manifolds, we relate the algebraic properties of invariant connections to their geometric properties. Using this perspective, we generalize some classical results of Cartan and Nomizu to invariant connections on algebroids. This has fundamental consequences for the theory of numerical integrators, giving a characterization of the spaces on which Butcher and Lie-Butcher series methods, which generalize Runge-Kutta methods, may be applied.


title: Middleware Building Blocks for Workflow Systems

abstract: This paper describes a building blocks approach to the design of scientific workflow systems. We discuss RADICAL-Cybertools as one implementation of the building blocks concept, showing how they are designed and developed in accordance with this approach. This paper offers three main contributions: (i) showing the relevance of the design principles underlying the building blocks approach to support scientific workflows on high performance computing platforms; (ii) illustrating a set of building blocks that enable multiple points of integration, "unifying" conceptual reasoning across otherwise very different tools and systems; and (iii) case studies discussing how RADICAL-Cybertools are integrated with existing workflow, workload, and general purpose computing systems and used to develop domain-specific workflow systems.


title: Correct power for cluster-randomized difference-in-difference trials with loss to follow-up

abstract: Cluster randomized trials with measurements at baseline can improve power over post-test only designs by using difference in difference designs. However, subjects may be lost to follow-up between the baseline and follow-up periods. While equations for sample size and variance have been developed assuming no loss to follow-up ("cohort") and completely different subjects at baseline and follow-up ("cross-sectional") difference in difference designs, equations have yet to be developed when some subjects are observed in both periods ("mixture" designs). We present a general equation for calculating the variance in difference in difference designs and derive special cases assuming loss to follow-up with replacement of lost subjects and assuming loss to follow-up with no replacement but retaining the baseline measurements of all subjects. Relative efficiency plots, plots of variance against subject autocorrelation, and plots of variance by follow-up rate and subject autocorrelation are used to compare cohort, cross-sectional, and mixture approaches. Results indicate that when loss to follow-up to uncommon, mixture designs are almost as efficient as cohort designs with a given initial sample size. When loss to follow-up is common, mixture designs with full replacement maintain efficiency relative to cohort designs. Finally, our results provide guidance on whether to replace lost subjects during trial design and analysis.


title: Gravitoelectromagnetism in metric f(R) and Brans-Dicke theories with a potential

abstract: A Gravitoelectromagnetism formalism in the context of metric f(R) theory is presented and the analogue Lorentz force law is derived. Some interesting results such as the dependence of the deviation from General Relativity on the absolute value of the scalar potential are found, it is also found that the f(R) effects are only relevant at a shorter distance or when the distance is much less than the compton wavelength, and that the effects are attractive in nature. An investigation of gravitational time delay in the context of metric f(R) is also presented showing that the Ricci scalar alone is responsible for the time delay effect which seems to suggest that the extra scalar degree of freedom associated to f(R) does not provide any modification. Also, to generalise our results, the Lorentz force law and gravitational time delay in the case of Brans-Dicke theories with a potential are derived; it is shown that the results are consistent with those obtained in the case of metric f(R) and General Relativity in the appropriate limits.


title: The matrix Dyson equation and its applications for random matrices

abstract: These lecture notes are a concise introduction of recent techniques to prove local spectral universality for a large class of random matrices. The general strategy is presented following the recent book with H.T. Yau. We extend the scope of this book by focusing on new techniques developed to deal with generalizations of Wigner matrices that allow for non-identically distributed entries and even for correlated entries. This requires to analyze a system of nonlinear equations, or more generally a nonlinear matrix equation called the Matrix Dyson Equation (MDE). We demonstrate that stability properties of the MDE play a central role in random matrix theory. The analysis of MDE is based upon joint works with J. Alt, O. Ajanki, D. Schr\"oder and T. Kr\"uger that are supported by the ERC Advanced Grant, RANMAT 338804 of the European Research Council. The lecture notes were written for the 27th Annual PCMI Summer Session on Random Matrices held in 2017. The current edited version will appear in the IAS/Park City Mathematics Series, Vol. 26.


title: A simplified general relativistic model to analyze the structure of spiral galaxies

abstract: The dark matter hypothesis, which is not called into question here, explains why typical rotation curves of spiral galaxies do not follow a Keplerian profile. It is however not sufficient in itself to explain why the whole matter distribution in spiral galaxies is such that the rotation curve generally presents a flat profile in the disk region. To understand this property, a model considering general relativistic effects is developed. It is stressed that the aim is not to explain the flat rotation curve of spiral galaxies without dark matter. More specifically, the analytical stationary solution of an axisymmetric rotating pressureless fluid for the linearized equations of the theory of general relativity is determined. It is demonstrated that this solution leads to some constraints on the rotation curve, by looking to its limit behavior when neglecting general relativistic effects. In particular, the positiveness of the density imposes the rotation curve to be flat in the regions where the density and general relativistic effects are small, such as in the disk region. General relativistic effects hence remain negligible in the disk region, but their consideration proved to be necessary to establish the aforementioned constraint. Such constraint cannot be derived from a Newtonian approach. The model is finally applied on two specific cases to demonstrate its ability to predict the rotation curve from a typical density profile along the galactic plane. These examples suggest that for some galaxies, general relativistic effects can be significant close to the bulge region and should be taken into account to have a proper understanding of their rotation curve.


title: Blow-up profile of neutron stars in the Hartree-Fock-Bogoliubov theory

abstract: We consider the gravitational collapse for neutron stars in the Hartree-Fock-Bogoliubov theory. We prove that when the number particle becomes large and the gravitational constant is small such that the attractive interaction strength approaches the Chandrasekhar limit mass slowly, the minimizers develop a universal blow-up profile. It is given by the Lane-Emden solution.


title: Optimal Linear Discriminators For The Discrete Choice Model In Growing Dimensions

abstract: Manski's celebrated maximum score estimator for the discrete choice model, which is an optimal linear discriminator, has been the focus of much investigation in both the econometrics and statistics literatures, but its behavior under growing dimension scenarios largely remains unknown. This paper addresses that gap. Two different cases are considered: $p$ grows with $n$ but at a slow rate, i.e. $p/n \rightarrow 0$; and $p \gg n$ (fast growth). In the binary response model, we recast Manski's score estimation as empirical risk minimization for a classification problem, and derive the $\ell_2$ rate of convergence of the score estimator under a \emph{transition condition} in terms of our margin parameter that calibrates the level of difficulty of the estimation problem. We also establish upper and lower bounds for the minimax $\ell_2$ error in the binary choice model that differ by a logarithmic factor, and construct a minimax-optimal estimator in the slow growth regime. Some extensions to the general case -- the multinomial response model -- are also considered. Last but not least, we use a variety of learning algorithms to compute the maximum score estimator in growing dimensions.


title: Omnipotent Virtual Giant for Remote Human-Swarm Interaction

abstract: This paper proposes an intuitive human-swarm interaction framework inspired by our childhood memory in which we interacted with living ants by changing their positions and environments as if we were omnipotent relative to the ants. In virtual reality, analogously, we can be a super-powered virtual giant who can supervise a swarm of mobile robots in a vast and remote environment by flying over or resizing the world and coordinate them by picking and placing a robot or creating virtual walls. This work implements this idea by using Virtual Reality along with Leap Motion, which is then validated by proof-of-concept experiments using real and virtual mobile robots in mixed reality. We conduct a usability analysis to quantify the effectiveness of the overall system as well as the individual interfaces proposed in this work. The results revealed that the proposed method is intuitive and feasible for interaction with swarm robots, but may require appropriate training for the new end-user interface device.


title: Dynamic intertemporal utility optimization by means of Riccati transformation of Hamilton-Jacobi Bellman equation

abstract: In this paper we investigate a dynamic stochastic portfolio optimization problem involving both the expected terminal utility and intertemporal utility maximization. We solve the problem by means of a solution to a fully nonlinear evolutionary Hamilton-Jacobi-Bellman (HJB) equation. We propose the so-called Riccati method for transformation of the fully nonlinear HJB equation into a quasi-linear parabolic equation with non-local terms involving the intertemporal utility function. As a numerical method we propose a semi-implicit scheme in time based on a finite volume approximation in the spatial variable. By analyzing an explicit traveling wave solution we show that the numerical method is of the second experimental order of convergence. As a practical application we compute optimal strategies for a portfolio investment problem motivated by market financial data of German DAX 30 Index and show the effect of considering intertemporal utility on optimal portfolio selection.


title: Photonic topological phase transition on demand

abstract: On-demand, switchable phase transitions between topologically non-trivial and trivial photonic states are demonstrated. Specifically, it is shown that integration of a 2D array of coupled ring resonators within a thermal heater array enables unparalleled control over topological protection of photonic modes. Importantly, auxiliary control over spatial phase modulation opens up a way to guide topologically protected edge modes along generated virtual boundaries. The proposed approach can lead to practical realizations of topological phase transitions in many photonic applications, including topologically protected photonic memory/logic devices, robust optical modulators, and switches.


title: An Analytical Model for Performance and Lifetime Estimation of Hybrid DRAM-NVM Main Memories

abstract: NVMs have promising advantages (e.g., lower idle power, higher density) over the existing predominant main memory technology, DRAM. Yet, NVMs also have disadvantages (e.g., limited endurance). System architects are therefore examining hybrid DRAM-NVM main memories to enable the advantages of NVMs while avoiding the disadvantages as much as possible. Unfortunately, the hybrid memory design space is very large and complex due to the existence of very different types of NVMs and their rapidly-changing characteristics. Therefore, optimization of performance and lifetime of hybrid memory based computing platforms and their experimental evaluation using traditional simulation methods can be very time-consuming and sometimes even impractical. As such, it is necessary to develop a fast and flexible analytical model to estimate the performance and lifetime of hybrid memories on various workloads. This paper presents an analytical model for hybrid memories based on Markov decision processes. The proposed model estimates the hit ratio and lifetime for various configurations of DRAM-NVM hybrid main memories. Our model also provides accurate estimation of the effect of data migration policies on the hybrid memory hit ratio, one of the most important factors in hybrid memory performance and lifetime. Such an analytical model can aid designers to tune hybrid memory configurations to improve performance and/or lifetime. We present several optimizations that make our model more efficient while maintaining its accuracy. Our experimental evaluations show that the proposed model (a) accurately predicts the hybrid memory hit ratio with an average error of 4.61% on a commodity server, (b) accurately estimates the NVM lifetime with an average error of 2.93%, and (c) is on average 4x faster than conventional state-of-the-art simulation platforms for hybrid memories.


title: Diophantine Problem in Some Metabelian Groups

abstract: In this paper we show that Diophantine problem for quadratic equations in Baumslag-Solitar groups $BS(1,k)$ and in wreath products $A \wr \mathbb{Z}$, where $A$ is a finitely generated abelian group and $\mathbb{Z}$ is an infinite cyclic group, is decidable.


title: Equivariant Degenerations of Plane Curve Orbits

abstract: In a series of papers, Aluffi and Faber computed the degree of the $GL_3$ orbit closure of an arbitrary plane curve. We attempt to generalize this to the equivariant setting by studying how orbits degenerate under some natural specializations, yielding a fairly complete picture in the case of plane quartics.


title: Double Sums of Kloosterman Sums in Finite Fields

abstract: We bound double sums of Kloosterman sums over a finite field ${\mathbb F}_{q}$, with one or both parameters ranging over an affine space over its prime subfield ${\mathbb F}_p \subseteq {\mathbb F}_{q} $. These are finite fields analogues of a series of recent results by various authors in finite fields and residue rings. Our results are based on recent advances in additive combinatorics in arbitrary finite field.


title: On the Performance of Mobility-Aware D2D Caching Networks

abstract: The increase in demand for spectrum-based services forms a bottleneck in wireless networks. Device-to-Device (D2D) caching networks tackle this problem by exploiting user's behavior predictability and the possibility of sharing data between them to alleviate the network congestion. However, capturing mobility statistics allows Service Providers (SPs) to enhance their caching strategies. In this work, we introduce a mobility-aware D2D caching network where SP harnesses user demand and mobility statistics to minimize the incurred service cost through an optimal caching policy. We investigate two caching schemes: centralized and decentralized caching schemes. In the centralized caching scheme, SP makes the caching decision towards its cost minimization to increase its profit. However, the complexity of optimal caching policy grows exponentially with the number of users. Therefore, we discuss a greedy caching algorithm which has a polynomial order complexity. We also use this greedy algorithm to establish upper and lower bounds on the proactive service gain achieved by the optimal caching policy. In the decentralized caching scheme, users take over and make their caching decisions, in a distributed fashion affected by the SP pricing policy, towards their payment minimization. We formulated the tension between the SP and users as a Stackelberg game. The Best response analysis was used to identify a sub-game perfect Nash equilibrium (SPNE) between users. The optimal solution of proposed model was found to depend on the SP reward preference, which affects the assigned memory in users' devices. We found some regimes for the reward value where the SPNE was non-unique. A fair allocation caching policy was adopted to choose one of these SPNEs. To understand the impact of user behavior, we investigated some special cases to explore how user's mobility statistics affect their caching decision.


title: On the adjoint Markov policies in stochastic differential games

abstract: We consider time-homogeneous uniformly nondegenerate stochastic differential games in domains and propose constructing $\varepsilon$-optimal strategies and policies by using adjoint Markov strategies and adjoint Markov policies which are actually time-homogeneous Markov, however, relative not to the original process but to a couple of processes governed by a system consisting of the main original equation and of an adjoint stochastic equations of the same type as the main one. We show how to find $\varepsilon$-optimal strategies and policies in these classes by using the solvability in Sobolev spaces of not the original Isaacs equation but of its appropriate modification. We also give an example of a uniformly nondegenerate game where our assumptions are not satisfied and where we conjecture that there are no not only optimal Markov but even $\varepsilon$-optimal adjoint (time-homogeneous) Markov strategies for one of the players.


title: One Bit Spectrum Sensing in Cognitive Radio Sensor Networks

abstract: This paper proposes a spectrum sensing algorithm from one bit measurements in a cognitive radio sensor network. A likelihood ratio test (LRT) for the one bit spectrum sensing problem is derived. Different from the one bit spectrum sensing research work in the literature, the signal is assumed to be a discrete random correlated Gaussian process, where the correlation is only available within immediate successive samples of the received signal. The employed model facilitates the design of a powerful detection criteria with measurable analytical performance. One bit spectrum sensing criterion is derived for one sensor which is then generalized to multiple sensors. Performance of the detector is analyzed by obtaining closed-form formulas for the probability of false alarm and the probability of detection. Simulation results corroborate the theoretical findings and confirm the efficacy of the proposed detector in the context of highly correlated signals and large number of sensors.


title: Zero supermode-based multipartite entanglement in $\chi^{(2)}$ nonlinear waveguides arrays

abstract: We show that arrays of $\chi^{(2)}$ nonlinear waveguides in the second harmonic generation regime are a promising source of continuous-variable entanglement. We indeed demonstrate analytically that optical arrays with odd number of waveguides injected with the zero-eigenvalue fundamental supermode entangle this fundamental supermode with a collective harmonic field. Moreover the fundamental individual modes are multipartite entangled and their entanglement grows with propagation length. The device is scalable, robust to losses, does not rely on specific values of nonlinearity and coupling and is easily realized with current technology. It thus stands as an unprecedented candidate for generation of multipartite continuous-variable entanglement for optical quantum information processing.


title: Machine Learning Methods Economists Should Know About

abstract: We discuss the relevance of the recent Machine Learning (ML) literature for economics and econometrics. First we discuss the differences in goals, methods and settings between the ML literature and the traditional econometrics and statistics literatures. Then we discuss some specific methods from the machine learning literature that we view as important for empirical researchers in economics. These include supervised learning methods for regression and classification, unsupervised learning methods, as well as matrix completion methods. Finally, we highlight newly developed methods at the intersection of ML and econometrics, methods that typically perform better than either off-the-shelf ML or more traditional econometric methods when applied to particular classes of problems, problems that include causal inference for average treatment effects, optimal policy estimation, and estimation of the counterfactual effect of price changes in consumer choice models.


title: Parametric Fokker-Planck equation

abstract: We derive the Fokker-Planck equation on the parametric space. It is the Wasserstein gradient flow of relative entropy on the statistical manifold. We pull back the PDE to a finite dimensional ODE on parameter space. Some analytical example and numerical examples are presented.


title: Truly Batch Apprenticeship Learning with Deep Successor Features

abstract: We introduce a novel apprenticeship learning algorithm to learn an expert's underlying reward structure in off-policy model-free \emph{batch} settings. Unlike existing methods that require a dynamics model or additional data acquisition for on-policy evaluation, our algorithm requires only the batch data of observed expert behavior. Such settings are common in real-world tasks---health care, finance or industrial processes ---where accurate simulators do not exist or data acquisition is costly. To address challenges in batch settings, we introduce Deep Successor Feature Networks(DSFN) that estimate feature expectations in an off-policy setting and a transition-regularized imitation network that produces a near-expert initial policy and an efficient feature representation. Our algorithm achieves superior results in batch settings on both control benchmarks and a vital clinical task of sepsis management in the Intensive Care Unit.


title: Besov class via heat semigroup on Dirichlet spaces III: BV functions and sub-Gaussian heat kernel estimates

abstract: With a view toward fractal spaces, by using a Korevaar-Schoen space approach, we introduce the class of bounded variation (BV) functions in a general framework of strongly local Dirichlet spaces with a heat kernel satisfying sub-Gaussian estimates. Under a weak Bakry-\'Emery curvature type condition, which is new in this setting, this BV class is identified with a heat semigroup based Besov class. As a consequence of this identification, properties of BV functions and associated BV measures are studied in detail. In particular, we prove co-area formulas, global $L^1$ Sobolev embeddings and isoperimetric inequalities. It is shown that for nested fractals or their direct products the BV class we define is dense in $L^1$. The examples of the unbounded Vicsek set, unbounded Sierpinski gasket and unbounded Sierpinski carpet are discussed.

categories: math.MG math-ph math.AP math.FA math.MP math.PR


title: Ensemble Methods for Causal Effects in Panel Data Settings

abstract: This paper studies a panel data setting where the goal is to estimate causal effects of an intervention by predicting the counterfactual values of outcomes for treated units, had they not received the treatment. Several approaches have been proposed for this problem, including regression methods, synthetic control methods and matrix completion methods. This paper considers an ensemble approach, and shows that it performs better than any of the individual methods in several economic datasets. Matrix completion methods are often given the most weight by the ensemble, but this clearly depends on the setting. We argue that ensemble methods present a fruitful direction for further research in the causal panel data setting.


title: Review of human decision-making during computer security incident analysis

abstract: We review practical advice on decision-making during computer security incident response. Scope includes standards from the IETF, ISO, FIRST, and the US intelligence community. To focus on human decision-making, the scope is the evidence collection, analysis, and reporting phases of response. The results indicate both strengths and gaps. A strength is available advice on how to accomplish many specific tasks. However, there is little guidance on how to prioritize tasks in limited time or how to interpret, generalize, and convincingly report results. Future work should focus on these gaps in explication and specification of decision-making during incident analysis.


title: Determining satisfiability of 3-SAT in polynomial time

abstract: In this paper, we provide a deterministic polynomial time algorithm that determines satisfiability of 3-SAT. The complexity analysis for the algorithm takes into account no efficiency and yet provides a low enough bound, that efficient versions are practical with respect to today's hardware. We accompany this paper with a serial version of the algorithm without non-trivial efficiencies (link: polynomial3sat.org).


title: Residual Non-local Attention Networks for Image Restoration

abstract: In this paper, we propose a residual non-local attention network for high-quality image restoration. Without considering the uneven distribution of information in the corrupted images, previous methods are restricted by local convolutional operation and equal treatment of spatial- and channel-wise features. To address this issue, we design local and non-local attention blocks to extract features that capture the long-range dependencies between pixels and pay more attention to the challenging parts. Specifically, we design trunk branch and (non-)local mask branch in each (non-)local attention block. The trunk branch is used to extract hierarchical features. Local and non-local mask branches aim to adaptively rescale these hierarchical features with mixed attentions. The local mask branch concentrates on more local structures with convolutional operations, while non-local attention considers more about long-range dependencies in the whole feature map. Furthermore, we propose residual local and non-local attention learning to train the very deep network, which further enhance the representation ability of the network. Our proposed method can be generalized for various image restoration applications, such as image denoising, demosaicing, compression artifacts reduction, and super-resolution. Experiments demonstrate that our method obtains comparable or better results compared with recently leading methods quantitatively and visually.


title: A Higher-Order Kolmogorov-Smirnov Test

abstract: We present an extension of the Kolmogorov-Smirnov (KS) two-sample test, which can be more sensitive to differences in the tails. Our test statistic is an integral probability metric (IPM) defined over a higher-order total variation ball, recovering the original KS test as its simplest case. We give an exact representer result for our IPM, which generalizes the fact that the original KS test statistic can be expressed in equivalent variational and CDF forms. For small enough orders ($k \leq 5$), we develop a linear-time algorithm for computing our higher-order KS test statistic; for all others ($k \geq 6$), we give a nearly linear-time approximation. We derive the asymptotic null distribution for our test, and show that our nearly linear-time approximation shares the same asymptotic null. Lastly, we complement our theory with numerical studies.


title: Probing the linear polarization of photons in ultraperipheral heavy ion collisions

abstract: We propose to measure the linear polarization of the external electromagnetic fields of a relativistic heavy ion through azimuthal asymmetries in dilepton production in ultraperipheral collisions. The asymmetries estimated with the equivalent photon approximation are shown to be sizable.


title: On resonance contribution to balance functions

abstract: It is known that resonance decays influence the shape of the charge-balance functions measured in hadronic collisions. That is reflected in their rapidity and azimuthal widths and the integral, and therefore has consequences for different model interpretations. In this paper, we show that the contribution from neutral resonance decays can be removed from the balance function in an analytical way, and test the performance of the removal procedure with PYTHIA events. Prospects for applications of the procedure to real data analysis of balance functions are also discussed.


title: The $B$ Anomalies and New Physics in $b \to s e^+ e^-$

abstract: We investigate the implications of the latest LHCb measurement of $R_K$ for NP explanations of the $B$ anomalies. The previous data could be explained if the $b \to s \mu^+ \mu^-$ NP is in (I) $C_{9,{\rm NP}}^{\mu\mu}$ or (II) $C_{9,{\rm NP}}^{\mu\mu} = -C_{10,{\rm NP}}^{\mu\mu}$, with scenario (I) providing a better explanation than scenario (II). This continues to hold with the new measurement of $R_K$. However, for both scenarios, this measurement leads to a slight tension of $O(1\sigma)$ between separate fits to the $b \to s \mu^+ \mu^-$ and $R_{K^{(*)}}$ data. In this paper, we investigate whether this tension can be alleviated with the addition of NP in $b \to s e^+ e^-$. In particular, we examine the effect of adding such NP to scenarios (I) and (II). We find several scenarios in which this leads to improvements in the fits. $Z'$ and LQ models with contributions to both $b \to s \mu^+ \mu^-$ and $b \to s e^+ e^-$ can reproduce the data, but only within scenarios based on (II). If the tension persists in future measurements, it may be necessary to consider NP models with more than one particle contributing to $b \to s \ell^+ \ell^-$.


title: Optimizing the trade-off between number of cops and capture time in Cops and Robbers

abstract: The cop throttling number $th_c(G)$ of a graph $G$ for the game of Cops and Robbers is the minimum of $k + capt_k(G)$, where $k$ is the number of cops and $capt_k(G)$ is the minimum number of rounds needed for $k$ cops to capture the robber on $G$ over all possible games in which both players play optimally. In this paper, we construct a family of graphs having $th_c(G)= \Omega(n^{2/3})$, establish a sublinear upper bound on the cop throttling number, and show that the cop throttling number of chordal graphs is $O(\sqrt{n})$. We also introduce the product cop throttling number $th_c^{\times}(G)$ as a parameter that minimizes the person-hours used by the cops. This parameter extends the notion of speed-up that has been studied in the context of parallel processing and network decontamination. We establish bounds on the product cop throttling number in terms of the cop throttling number, characterize graphs with low product cop throttling number, and show that for a chordal graph $G$, $th_c^{\times}=1+rad(G)$.


title: Fast Statistical Iterative Reconstruction for MVCT in TomoTherapy

abstract: Statistical iterative reconstruction is expected to improve the image quality of megavoltage computed tomography (MVCT). However, one of the challenges of iterative reconstruction is its large computational cost. The purpose of this work is to develop a fast iterative reconstruction algorithm by combining several iterative techniques and by optimizing reconstruction parameters. Megavolt projection data was acquired from a TomoTherapy system and reconstructed using our statistical iterative reconstruction. Total variation was used as the regularization term and the weight of the regularization term was determined by evaluating signal-to-noise ratio (SNR), contrast-to-noise ratio (CNR), and visual assessment of spatial resolution using Gammex and Cheese phantoms. Gradient decent with an adaptive convergence parameter, ordered subset expectation maximization (OSEM), and CPU/GPU parallelization were applied in order to accelerate the present reconstruction algorithm. The SNR and CNR of the iterative reconstruction were several times better than that of filtered back projection (FBP). The GPU parallelization code combined with the OSEM algorithm reconstructed an image several hundred times faster than a CPU calculation. With 500 iterations, which provided good convergence, our method produced a 512$\times$512 pixel image within a few seconds. The image quality of the present algorithm was much better than that of FBP for patient data. An image from the iterative reconstruction in TomoTherapy can be obtained within few seconds by fine-tuning the parameters. The iterative reconstruction with GPU was fast enough for clinical use, and largely improve the MVCT images.


title: Feedback-Induced Quantum Phase Transitions Using Weak Measurements

abstract: We show that applying feedback and weak measurements to a quantum system induces phase transitions beyond the dissipative ones. Feedback enables controlling essentially quantum properties of the transition, i.e., its critical exponent, as it is driven by the fundamental quantum fluctuations due to measurement. Feedback provides the non-Markovianity and nonlinearity to the hybrid quantum-classical system, and enables simulating effects similar to spin-bath problems and Floquet time crystals with tunable long-range (long-memory) interactions.

categories: quant-ph cond-mat.quant-gas physics.atom-ph physics.optics


title: Travelling wave solutions in a negative nonlinear diffusion-reaction model

abstract: We use a geometric approach to prove the existence of smooth travelling wave solutions of a nonlinear diffusion-reaction equation with logistic kinetics and a convex nonlinear diffusivity function which changes sign twice in our domain of interest. We determine the minimum wave speed, c*, and investigate its relation to the spectral stability of the travelling wave solutions.


title: Neutral surface topology

abstract: Neutral surfaces, along which most of the mixing in the ocean occurs, are notoriously difficult objects: they do not exist as well-defined surfaces, and as such can only be approximated. In a hypothetical ocean where neutral surfaces are well-defined, the in-situ density on the surface is a multivalued function of the pressure on the surface, $\underset{\sim}{p}$. The surface is decomposed into geographic regions where there is one connected pressure contour per pressure value, making this function single-valued in each region. The regions are represented by arcs of the Reeb graph of $\underset{\sim}{p}$. The regions meet at saddles of $\underset{\sim}{p}$ which are represented by internal nodes of the Reeb graph. Leaf nodes represent extrema of $\underset{\sim}{p}$. Cycles in the Reeb graph are created by islands and other holes in the neutral surface. This topological theory of neutral surfaces is used to create a new class of approximately neutral surfaces in the real ocean, called topobaric surfaces, which are very close to neutral and fast to compute. Topobaric surfaces are the topologically correct extension of orthobaric density surfaces to be geographically dependent, which is fundamental to neutral surfaces. Also considered is the possibility that helical neutral trajectories might have a larger pitch around islands than in the open ocean.


title: An Exact No Free Lunch Theorem for Community Detection

abstract: A precondition for a No Free Lunch theorem is evaluation with a loss function which does not assume a priori superiority of some outputs over others. A previous result for community detection by Peel et al. (2017) relies on a mismatch between the loss function and the problem domain. The loss function computes an expectation over only a subset of the universe of possible outputs; thus, it is only asymptotically appropriate with respect to the problem size. By using the correct random model for the problem domain, we provide a stronger, exact No Free Lunch theorem for community detection. The claim generalizes to other set-partitioning tasks including core/periphery separation, $k$-clustering, and graph partitioning. Finally, we review the literature of proposed evaluation functions and identify functions which (perhaps with slight modifications) are compatible with an exact No Free Lunch theorem.


title: Laws of large numbers in the Raise and Peel model

abstract: We establish the exact laws of large numbers for two time additive quantities in the Raise and Peel model, the number of tiles removed by avalanches and the number of global avalanches happened by given time. The validity of conjectures for the related stationary state correlation functions then follow. The proof is based on the technique of Baxter's T-Q equation applied to the associated XXZ chain and on its solution at $\Delta=-1/2$ obtained by Fridkin, Stroganov and Zagier.

categories: math-ph cond-mat.stat-mech math.MP nlin.SI


title: Discrete para-product operators on variable Hardy spaces

abstract: Let $p(\cdot):\mathbb R^n\rightarrow(0,\infty)$ be a variable exponent function satisfying the globally log-H\"older continuous condition. In this paper, we obtain the boundedness of para-product operators $\pi_b$ on variable Hardy spaces $H^{p(\cdot)}(\mathbb R^n)$, where $b\in BMO(\mathbb R^n)$. As an application, we show that non-convolution type Calder\'on-Zygmund operators $T$ are bounded on $H^{p(\cdot)}(\mathbb R^n)$ if and only if $T^\ast1=0$, where $\frac{n}{n+\epsilon}<\mbox{essinf}_{x\in\mathbb R^n} p\le \mbox{esssup}_{x\in\mathbb R^n} p\le 1$, $\epsilon$ is the regular exponent of kernel of $T$. Our approach relies on the discrete version of Calder\'on's reproducing formula, discrete Littlewood-Paley-Stein theory and almost orthogonal estimates. These results still hold for variable Hardy space on spaces of homogeneous type by using our methods.


abstract: McDougall (1989) proved that neutral surfaces possess an exact geostrophic streamfunction, but its form has remained unknown. The exact geostrophic streamfunction for neutral surfaces is derived here. It involves a path integral of the specific volume along a neutral trajectory. On a neutral surface, the specific volume is a multivalued function of the pressure on the surface, $\underset{\sim}{p}$. By decomposing the neutral surface into regions where the specific volume is a single-valued function of $\underset{\sim}{p}$, the path integral is simply a sum of integrals of these single-valued functions. The regions are determined by the Reeb graph of $\underset{\sim}{p}$, and the neutral trajectory is encoded by a walk on this graph. Islands, and other holes in the neutral surface, can create cycles in the Reeb graph, causing the exact geostrophic streamfunction on a neutral surface to be multivalued. Moreover, neutral surfaces are ill-defined in the real ocean. Hence, the topobaric geostrophic streamfunction is presented: a single-valued approximation of the exact geostrophic streamfunction for neutral surfaces, for use on any well-defined, approximately neutral surface. Numerical tests on several approximately neutral surfaces reveal that the topobaric geostrophic streamfunction estimates the geostrophic velocity with an error that is about an order of magnitude smaller than that for any previously known geostrophic streamfunction. Also, the Montgomery potential is generalized, leading to an alternate form of the exact geostrophic streamfunction for neutral surfaces. This form is used to construct the orthobaric Montgomery potential, an easily computable geostrophic streamfunction that estimates the geostrophic velocity more accurately than any previously known geostrophic streamfunction, but less so than the topobaric geostrophic streamfunction.


title: Convergence of generalized Collatz problem in k-adic field

abstract: In this article, we define a new k-adic series transformation called Z-transformation and probe into its fixed point and periodicity. We extend the number field of the transform period problem to a wider k-adic field. Different constraints are imposed on k, then different periodic columns are formed after finite Z-transformations. We obtain that their periodic sequences are M1= {1,2} and M2={1,2}\cup {n_0}\cup {n'} respectively after derivation. As an application, it can provide a reference for C problems in more complex algebraic ystems.