title: Simulation Study of Hemispheric Phase-Asymmetry in the Solar Cycle

abstract: Observations of the sun suggest that solar activities systematically create north-south hemispheric asymmetries. For instance, the hemisphere in which the sunspot activity is more active tends to switch after the early half of each solar cycle. Svalgaard & Kamide (2013) recently pointed out that the time gaps of polar field reversal between the north and south hemispheres are simply consequences of the asymmetry of sunspot activity. However, the mechanism underlying the asymmetric feature in solar cycle activities is not yet well understood. In this paper, in order to explain the cause of the asymmetry from the theoretical point of view, we investigate the relationship between the dipole- and quadrupole-type components of the magnetic field in the solar cycle using the mean-field theory based on the flux transport dynamo model. As a result, we found that there are two different attractors of the solar cycle, in which either the north or the south polar field is first reversed, and that the flux transport dynamo model well explains the phase-asymmetry of sunspot activity and the polar field reversal without any ad hoc source of asymmetry.


title: Local Uniqueness and Refined Spike Profiles of Ground States for Two-Dimensional Attractive Bose-Einstein Condensates

abstract: We consider ground states of two-dimensional Bose-Einstein condensates in a trap with attractive interactions, which can be described equivalently by positive minimizers of the $L^2-$critical constraint Gross-Pitaevskii energy functional. It is known that ground states exist if and only if $a< a^*:= \|w\|_2^2$, where $a$ denotes the interaction strength and $w$ is the unique positive solution of $\Delta w-w+w^3=0$ in $R^2$. In this paper, we prove the local uniqueness and refined spike profiles of ground states as $a\nearrow a^*$, provided that the trapping potential $h(x)$ is homogeneous and $H(y)=\int_{R^2} h(x+y)w^2(x)dx$ admits a unique and non-degenerate critical point.


title: Long-range dispersion effects on the water/vapor interface simulated using the most common models

abstract: The long-range contribution to dispersion forces is known to have a major impact on the properties of inhomogeneous fluids, and its correct treatment is increasingly recognized as being a necessary requirement to avoid cutoff-related artefacts. Although analytical corrections for quantities like the surface tension are known, these can not take into account the structural changes induced by the long-range contributions. Here, we analyze the interfacial properties of seven popular water models, comparing the results with the cut-off version of the dispersion potential. The differences in surface tension estimates are in all cases found to be less than 2 mN/m.


title: On warped product manifolds satisfying some pseudosymmetric type conditions

abstract: The object of the present paper is to study the characterization of warped product manifolds satisfying some pseudosymmetric type conditions, especially, due to projective curvature tensor. For this purpose we consider a warped product manifold satisfying the pseudosymmetric type condition $R\cdot R = L_1 Q(g,R) + L_2 Q(S,R)$ and evaluate its characterization theorem. As special cases of $L_1$ and $L_2$ we find out the necessary and sufficient condition for a warped product manifold to satisfy various pseudosymmetric type, such as pseudosymmetry, Ricci generalized pseudosymmetry, semisymmetry due to projective curvature tensor ($P\cdot R = 0$), pseudosymmetry due to projective curvature tensor ($P\cdot R = L Q(g,R)$) etc. Finally we present some suitable examples of warped product manifolds satisfying such pseudosymmetric type conditions.


title: Polychromatic colorings of complete graphs with respect to 1-,2-factors and Hamiltonian cycles

abstract: If G is a graph and H is a set of subgraphs of G, then an edge-coloring of G is called H-polychromatic if every graph from H gets all colors present in G on its edges. The H-polychromatic number of G, denoted poly_H(G), is the largest number of colors in an H-polychromatic coloring. In this paper, poly_H(G) is determined exactly when G is a complete graph and H is the family of all 1-factors. In addition poly_H(G) is found up to an additive constant term when G is a complete graph and H is the family of all 2-factors, or the family of all Hamiltonian cycles.


title: Ideal Magnetohydrodynamics with Radiative Terms: Energy Conditions

abstract: Nowadays, the magnetic and radiation fields are very important to understand the matter accretion into compact objects, the dynamics of binary systems, the equilibrium configurations of neutron stars, the photon diffusion, etc. The energy and the momentum associated to these fields, along with the matter one, need to satisfy some conditions that guarantee an appropriate physical behavior of the source and its gravitational field. Based on this fact, we present the energy conditions for a perfect fluid with magnetic and radiation field, in which the radiation part of the energy-momentum tensor is assumed to be approximately isotropic, in accordance with the optically thick regime. In order to find these conditions, the stress tensor of the system is written in an orthonormal basis in which it becomes diagonal, and the energy conditions are computed through contractions of the energy-momentum tensor with the four velocity vector of an arbitrary observer. Finally, the conditions for a magnetized fluid are presented as a particular case in which the radiation contribution is zero.


title: A randomized FEAST algorithm for generalized eigenvalue problems

abstract: The FEAST algorithm, due to Polizzi, is a typical contour-integral based eigensolver for computing the eigenvalues, along with their eigenvectors, inside a given region in the complex plane. It was formulated under the circumstance that the considered eigenproblem is Hermitian. The FEAST algorithm is stable and accurate, and has attracted much attention in recent years. However, it was observed that the FEAST algorithm may fail to find the target eigenpairs when applying it to the non-Hermitian problems. Efforts have been made to adapt the FEAST algorithm to non-Hermitian cases. In this work, we develop a new non-Hermitian scheme for the FEAST algorithm. The mathematical framework will be established, and the convergence analysis of our new method will be studied. Numerical experiments are reported to demonstrate the effectiveness of our method and to validate the convergence properties.


title: Gradient Coding

abstract: We propose a novel coding theoretic framework for mitigating stragglers in distributed learning. We show how carefully replicating data blocks and coding across gradients can provide tolerance to failures and stragglers for Synchronous Gradient Descent. We implement our schemes in python (using MPI) to run on Amazon EC2, and show how we compare against baseline approaches in running time and generalization error.

categories: stat.ML cs.DC cs.IT cs.LG math.IT stat.CO


title: An Extension of Generalized Linear Models to Finite Mixture Outcome Distributions

abstract: Finite mixture distributions arise in sampling a heterogeneous population. Data drawn from such a population will exhibit extra variability relative to any single subpopulation. Statistical models based on finite mixtures can assist in the analysis of categorical and count outcomes when standard generalized linear models (GLMs) cannot adequately account for variability observed in the data. We propose an extension of GLM where the response is assumed to follow a finite mixture distribution, while the regression of interest is linked to the mixture's mean. This approach may be preferred over a finite mixture of regressions when the population mean is the quantity of interest; here, only a single regression function must be specified and interpreted in the analysis. A technical challenge is that the mean of a finite mixture is a composite parameter which does not appear explicitly in the density. The proposed model is completely likelihood-based and maintains the link to the regression through a certain random effects structure. We consider typical GLM cases where means are either real-valued, constrained to be positive, or constrained to be on the unit interval. The resulting model is applied to two example datasets through a Bayesian analysis: one with success/failure outcomes and one with count outcomes. Supporting the extra variation is seen to improve residual plots and to appropriately widen prediction intervals.


title: Persistence of translational symmetry in the BCS model with radial pair interaction

abstract: We consider the two-dimensional BCS functional with a radial pair interaction. We show that the translational symmetry is not broken in a certain temperature interval below the critical temperature. In the case of vanishing angular momentum our results carry over to the three-dimensional case.


title: On the generalized porous medium equation in Fourier-Besov spaces

abstract: We study a kind of generalized porous medium equation with fractional Laplacian and abstract pressure term. For a large class of equations corresponding to the form: $u_t+\nu \Lambda^{\beta}u=\nabla\cdot(u\nabla Pu)$, we get their local well-posedness in Fourier-Besov spaces for large initial data. If the initial data is small, then the solution becomes global. Furthermore, we prove a blowup criterion for the solutions.


title: Radiation Pattern of Radio and Optical Components of Extended Radio Sources

abstract: The relation between parameters the D/sqrt(I) and Ic/Isum and radiation patterns of the optical and radio components of an extended radio source is analyzed, where D and I are the apparent size and observed radiation intensity of the source or its components respectively. The parameters of the pattern in the optical and radio (1.4 GHz) ranges are estimated. The radiation pattern of extended radio-emitting regions is close to spherical and the radiation of the central component is concentrated in a 24 degrees wide beam. Its luminosity is a factor of 4.58 higher than that of the extended component of the radio source. The radiation pattern of the optical component of the radio source turned out to be unexpectedly non-spherical: the main lobe of the pattern is about 26 degrees wide. The g-band luminosity is 6.4-12.3 times higher than the luminosity of the spherical fraction of the "optical" radiation pattern. A list of 116 new giant radio sources is presented.


title: Arduino control of a pulsatile flow rig

abstract: This note describes the design and testing of a programmable pulsatile flow pump using an Arduino micro-controller. The goal of this work is to build a compact and affordable system that can relatively easily be programmed to generate physiological waveforms. The system described here was designed to be used in an in-vitro set-up for vascular access hemodynamics research, and hence incorporates a gear pump that can deliver up to 1.5 l/min in a test flow loop. After a number of simple identification experiments to assess the dynamic behaviour of the system, a feed-forward control routine was implemented. The resulting system was shown to be able to produce the targeted patient-specific waveform with less than 3.6% error. Finally, we outline how to further increase the accuracy of the system, and how to adapt it to specific user needs.


title: Particle velocity based hydrofracturing algorithm for a penny-shaped crack

abstract: A universal particle velocity based algorithm for simulating hydraulic fracture with leak-off, previously demonstrated for the PKN and KGD models, is extended to obtain solutions for a penny-shaped crack. The numerical scheme is capable of dealing with both the viscosity and toughness dominated regimes, with the fracture being driven by a power-law fluid. The computational approach utilizes two dependent variables; the fracture aperture and the reduced particle velocity. The latter allows for the application of a local condition of the Stefan type (the speed equation) to trace the fracture front. The obtained numerical solutions are carefully tested using various methods, and are shown to achieve a high level of accuracy. Simple, accurate, semi-analytical approximations of the solution are provided for the zero leak-off case. A comparison with other results available in the literature is undertaken.


title: GraCT: A Grammar based Compressed representation of Trajectories

abstract: We present a compressed data structure to store free trajectories of moving objects (ships over the sea, for example) allowing spatio-temporal queries. Our method, GraCT, uses a $k^2$-tree to store the absolute positions of all objects at regular time intervals (snapshots), whereas the positions between snapshots are represented as logs of relative movements compressed with Re-Pair. Our experimental evaluation shows important savings in space and time with respect to a fair baseline.


title: Negative definite functions for C*-dynamical systems

abstract: Given an action $\alpha$ of a discrete group $G$ on a unital C*-algebra $A$, we introduce a natural concept of $\alpha$-negative definiteness for functions from $G$ to $A$, and examine some of the first consequences of such a notion. In particular, we prove analogs of theorems due to Delorme-Guichardet and Schoenberg in the classical case where $A$ is trivial. We also give a characterization of the Haagerup property for the action $\alpha$ when $G$ is countable.


title: MITHRA 1.0: A full-wave simulation tool for Free Electron Lasers

abstract: Free Electron Lasers (FELs) are a solution for providing intense, coherent and bright radiation in the hard X-ray regime. Due to the low wall-plug efficiency of FEL facilities, it is crucial and additionally very useful to develop complete and accurate simulation tools for better optimizing a FEL interaction. The highly sophisticated dynamics involved in a FEL process was the main obstacle hindering the development of general simulation tools for this problem. We present a numerical algorithm based on finite difference time domain/Particle in cell (FDTD/PIC) in a Lorentz boosted coordinate system which is able to fulfil a full-wave simulation of a FEL process. The developed software offers a suitable tool for the analysis of FEL interactions without considering any of the usual approximations. A coordinate transformation to bunch rest frame makes the very different length scales of bunch size, optical wavelengths and the undulator period transform to values with the same order. Consequently, FDTD/PIC simulations in conjunction with efficient parallelization techniques make the full-wave simulation feasible using the available computational resources. Several examples of free electron lasers are analyzed using the developed software, the results are benchmarked based on standard FEL codes and discussed in detail.


title: Cyclic Cohomology Groups of Some Self-similar Sets

abstract: We define a variant of the Young integration on some kinds of self-similar sets which are called cellular self-similar sets. This variant is an analogue of the Young integration defined on the unit interval. We give the criteria of the variant on cellular self-similar sets, and also show that the variant is a cyclic 1-cocycle of the algebra of complex-valued H\"older continuous functions on the cellular self-similar sets. This suggests that the cocycle is a variant of currents.


title: Monet: A User-oriented Behavior-based Malware Variants Detection System for Android

abstract: Android, the most popular mobile OS, has around 78% of the mobile market share. Due to its popularity, it attracts many malware attacks. In fact, people have discovered around one million new malware samples per quarter, and it was reported that over 98% of these new malware samples are in fact "derivatives" (or variants) from existing malware families. In this paper, we first show that runtime behaviors of malware's core functionalities are in fact similar within a malware family. Hence, we propose a framework to combine "runtime behavior" with "static structures" to detect malware variants. We present the design and implementation of MONET, which has a client and a backend server module. The client module is a lightweight, in-device app for behavior monitoring and signature generation, and we realize this using two novel interception techniques. The backend server is responsible for large scale malware detection. We collect 3723 malware samples and top 500 benign apps to carry out extensive experiments of detecting malware variants and defending against malware transformation. Our experiments show that MONET can achieve around 99% accuracy in detecting malware variants. Furthermore, it can defend against 10 different obfuscation and transformation techniques, while only incurs around 7% performance overhead and about 3% battery overhead. More importantly, MONET will automatically alert users with intrusion details so to prevent further malicious behaviors.


title: Microphase Separation in Random Multiblock Copolymers

abstract: Microphase separation in random multiblock copolymers is studied with mean-field theory assuming that long blocks of a copolymer are strongly segregated, whereas short blocks are able to penetrate into "alien" domains and exchange between the domains and interfacial layer. A bidisperse copolymer with blocks of only two sizes (long and short) is considered as a model of multiblock copolymers with high polydispersity in the block size. Short blocks of the copolymer play an important role in microphase separation. First, their penetration into the "alien" domains leads to the formation of joint long blocks in their own domains. Second, short blocks localized at the interface considerably change the interfacial tension. The possibility of penetration of short blocks into the "alien" domains is controlled by the product chi*Nsh (chi is the Flory-Huggins interaction parameter, Nsh is the short block length). At not very large chi*Nsh, the domain size is larger than that for a regular copolymer consisting of the same long blocks as in the considered random copolymer. At a fixed mean block size, the domain size grows with an increase in the block size dispersity, the rate of the growth being dependent of the more detailed parameters of the block size distribution.


title: Differential flatness for neuroscience population dynamics -- A preliminary study

abstract: The present document is devoted to structural properties of neural population dynamics and especially their differential flatness. Several applications of differential flatness in the present context can be envisioned, among which: trajectory tracking, feedforward to feedback switching, cyclic character, positivity and boundedness.


title: SUSY $N$-supergroups and their real forms

abstract: We study SUSY $N$-supergroups, $N=1,2$, their classification and explicit realization, together with their real forms. In the end, we give the supergroup of SUSY preserving automorphism of $\mathbf{C}^{1|1}$ and we identify it with a subsupergroup of the SUSY preserving automorphisms of $\mathbf{P}^{1|1}$.


title: Label Visualization and Exploration in IR

abstract: There is a renaissance in visual analytics systems for data analysis and sharing, in particular, in the current wave of big data applications. We introduce RAVE, a prototype that automates the generation of an interface that uses facets and visualization techniques for exploring and analyzing relevance assessments data sets collected via crowdsourcing. We present a technical description of the main components and demonstrate its use.


title: Deformation in amorphous-crystalline nanolaminates -- an effective-temperature theory and interaction between defects

abstract: Experiments and atomic-scale simulations suggest that the transmission of plasticity carriers in deforming amorphous-crystalline nanolaminates is mediated by the biphase interface between the amorphous and crystalline layers. In this paper, we present a micromechanics model for these biphase nanolaminates that describes defect interactions through the amorphous-crystalline interface (ACI). The model is based on an effective-temperature framework to achieve a unified description of the slow, configurational atomic rearrangements in both phases when driven out of equilibrium. We show how the second law of thermodynamics constrains the density of defects and the rate of configurational rearrangements, and apply this framework to dislocations in crystalline solids and shear transformation zones (STZs) in amorphous materials. The effective-temperature formulation enables us to interpret the observed movement of dislocations to the ACI and the production of STZs at the interface as a "diffusion" of configurational disorder across the material. We demonstrate favorable agreement with experimental findings reported in (Kim et al., Adv. Funct. Mater., 2011), and demonstrate how the ACI acts as a sink of dislocations and a source of STZs.


title: Limits in Categories of Vietoris Coalgebras

abstract: Motivated by the need to reason about hybrid systems, we study limits in categories of coalgebras whose underlying functor is a Vietoris polynomial one - intuitively, the topological analogue of a Kripke polynomial functor. Among other results, we prove that every Vietoris polynomial functor admits a final coalgebra if it respects certain conditions concerning separation axioms and compactness. When the functor is restricted to some of the categories induced by these conditions the resulting categories of coalgebras are even complete. As a practical application, we use these developments in the specification and analysis of non-deterministic hybrid systems, in particular to obtain suitable notions of stability, and behaviour.


title: Anytime Monte Carlo

abstract: Monte Carlo algorithms simulate some prescribed number of samples, taking some random real time to complete the computations necessary. This work considers the converse: to impose a real-time budget on the computation, which results in the number of samples simulated being random. To complicate matters, the real time taken for each simulation may depend on the sample produced, so that the samples themselves are not independent of their number, and a length bias with respect to compute time is apparent. This is especially problematic when a Markov chain Monte Carlo (MCMC) algorithm is used and the final state of the Markov chain -- rather than an average over all states -- is required, which is the case in parallel tempering implementations of MCMC. The length bias does not diminish with the compute budget in this case. It also occurs in sequential Monte Carlo (SMC) algorithms, which is the focus of this paper. We propose an anytime framework to address the concern, using a continuous-time Markov jump process to study the progress of the computation in real time. We first show that for any MCMC algorithm, the length bias of the final state's distribution due to the imposed real-time computing budget can be eliminated by using a multiple chain construction. The utility of this construction is then demonstrated on a large-scale SMC^2 implementation, using four billion particles distributed across a cluster of 128 graphics processing units on the Amazon EC2 service. The anytime framework imposes a real-time budget on the MCMC move steps within the SMC$^{2}$ algorithm, ensuring that all processors are simultaneously ready for the resampling step, demonstrably reducing idleness to due waiting times and providing substantial control over the total compute budget.


title: NuSTAR and XMM-Newton observations of the Arches cluster in 2015: fading hard X-ray emission from the molecular cloud

abstract: We present results of long NuSTAR (200 ks) and XMM-Newton (100 ks) observations of the Arches stellar cluster, a source of bright thermal (kT~2 keV) X-rays with prominent Fe XXV K_alpha 6.7 keV line emission and a nearby molecular cloud, characterized by an extended non-thermal hard X-ray continuum and fluorescent Fe K_alpha 6.4 keV line of a neutral or low ionization state material around the cluster. Our analysis demonstrates that the non-thermal emission of the Arches cloud underwent a dramatic change, with its homogeneous morphology, traced by fluorescent Fe K_alpha line emission, vanishing after 2012, revealing three bright clumps. The declining trend of the cloud emission, if linearly fitted, is consistent with half-life decay time of ~8 years. Such strong variations have been observed in several other molecular clouds in the Galactic Centre, including the giant molecular cloud Sgr B2, and point toward a similar propagation of illuminating fronts, presumably induced by the past flaring activity of Sgr A*.


title: Rings beyond the giant planets

abstract: Until 2013, only the giant planets were known to host ring systems. In June 2013, a stellar occulation revealed the presence of narrow and dense rings around Chariklo, a small Centaur object that orbits between Saturn and Uranus. Meanwhile, the Cassini spacecraft revealed evidence for the possible past presence of rings around the Saturnian satellites Rhea and Iapetus. Mars and Pluto are expected to have tenuous dusty rings, though they have so far evaded detection. More remotely, transit events observed around a star in 2007 may have revealed for the first time exoplanetary rings around a giant planet orbiting that star. So, evidence is building to show that rings are more common features in the universe than previously thought. Several interesting issues arising from the discovery (or suspicion) of new ring systems are described in this chapter.


title: A new approach to ultrasensitive gravitational wave detection

abstract: We suggest here a method to detect gravitational waves (GW) different from the interferometric approach. It is based on two critical steps: conversion of the GW action into rotational motion and subsequent conversion into electric current. The ability to detect extremely tiny currents empowers this suggested approach in which the gravitational wave moves charges around closed loops. This new method may allow detection of gravitational waves with strain sensitivities beyond the reach of the interferometric approaches.


title: Sur l'annulation de la valeur centrale de la fonction L de Hecke

abstract: In this note, following results from Henri Cohen and Winfried Kohnen, we show that for all integer k greater than 12 and divisible by 4, there exists a cuspidal eigenform of weight k for the full modular group SL2(Z) such that its Hecke L-function does not vanish on k/2.


title: Fisher information of two superconducting charge qubits under dephasing noisy

abstract: The dynamics of a charged two-qubit system prepared initially in a maximum entangled state is discussed, where each qubit interacts independently with a dephasing channel. The Fisher information is used to estimate the channel and the energy parameters. Moreover, the contribution from different parts of Fisher information is discussed. We show that, the degree of estimation of any parameter depends on the initial value of this parameter. It turns out that the upper bounds of Fisher information with respect to the dephasing parameter increase and that corresponding to coupling parameter decreases as the initial energy of each qubit increases.


title: Ribbon-move-unknotting-number-two 2-knots, pass-move-unknotting-number-two 1-knots, and high dimensional analogue

abstract: The (ordinary) unknotting-number of 1-dimensional knots, which is defined by using the crossing-change, is a very basic and important invariant. It is very natural to consider the `unknotting-number' associated with other local-moves on n-dimensional knots, where n is a natural number. In this paper we prove the following facts. For the ribbon-move on 2-knots, which is a kind of local-move on knots, we have the following: There is a ribbon-move-unknotting-number-two 2-knot. The ribbon-move-unknotting-number of 2-knots is unbounded. For the pass-move on 1-knots, which is a kind of local-move on knots, we have the following: There is a pass-move-unknotting-number-two 1-knot whose (ordinary) unknotting-number is 4. For any natural number n, there is a 1-knot whose pass-move-unknotting-number is>n and whose (ordinary) unknotting-number is 4n. For the high-dimensional-pass-move on high-dimensional knots, which is a kind of local-move on knots, we have the following: There is a (2k+1, 2k+2)-pass-move-unknotting-number-two (4k+2)-knot. The (2k+1, 2k+2)-pass-move-unknotting-number of (4k+2)-knot is unbounded. There is a (2k+1, 2k+1)-pass-move-unknotting-number-two (4k+1)-knot. The (2k+1, 2k+1)-pass-move-unknotting-number of (4k+1)-knot is unbounded. There is a (4k+1)-knot whose twist-move-unknotting-number is n for any natural number n.


title: Richard Dedekind: Style and Influence

abstract: This text is based on an invited talk at the Dedekind Symposium at Braunschweig in October 2016. It summarizes views from my recent commented edition of Dedekinds two books on the foundations of mathematics.


title: Yosida Duality

abstract: In this note we prove Yosida duality --- that is: the category of compact Hausdorff spaces with continuous maps is dually equivalent to the category of uniformly complete Archimedean Riesz spaces with distinguished units and unit-preserving Riesz homomorphisms between them.


title: Knowledge Elicitation via Sequential Probabilistic Inference for High-Dimensional Prediction

abstract: Prediction in a small-sized sample with a large number of covariates, the "small n, large p" problem, is challenging. This setting is encountered in multiple applications, such as precision medicine, where obtaining additional samples can be extremely costly or even impossible, and extensive research effort has recently been dedicated to finding principled solutions for accurate prediction. However, a valuable source of additional information, domain experts, has not yet been efficiently exploited. We formulate knowledge elicitation generally as a probabilistic inference process, where expert knowledge is sequentially queried to improve predictions. In the specific case of sparse linear regression, where we assume the expert has knowledge about the values of the regression coefficients or about the relevance of the features, we propose an algorithm and computational approximation for fast and efficient interaction, which sequentially identifies the most informative features on which to query expert knowledge. Evaluations of our method in experiments with simulated and real users show improved prediction accuracy already with a small effort from the expert.


title: Uniqueness of the $1$D Compressible to Incompressible Limit

abstract: Consider two compressible immiscible fluids in 1D in the isentropic approximation. The first fluid is surrounded and in contact with the second one. As the Mach number of the first fluid vanishes, the coupled dynamics of the two fluids results as the compressible to incompressible limit and is known to satisfy an ODE--PDE system. Below, a characterization of this limit is provided, ensuring its uniqueness.


title: Devil's Staircase Continuum in the Chiral Clock Spin Glass with Competing Ferromagnetic-Antiferromagnetic and Left-Right Chiral Interactions

abstract: The chiral clock spin-glass model with q=5 states, with both competing ferromagnetic-antiferromagnetic and left-right chiral frustrations, is studied in d=3 spatial dimensions by renormalization-group theory. The global phase diagram is calculated in temperature, antiferromagnetic bond concentration p, random chirality strength, and right-chirality concentration c. The system has a ferromagnetic phase, a multitude of different chiral phases, a chiral spin-glass phase, and a critical (algebraically) ordered phase. The ferromagnetic and chiral phases accumulate at the disordered phase boundary and form a spectrum of devil's staircases, where different ordered phases characteristically intercede at all scales of phase-diagram space. Shallow and deep reentrances of the disordered phase, bordered by fragments of regular and temperature-inverted devil's staircases, are seen. The extremely rich phase diagrams are presented as continuously and qualitatively changing videos.


title: Hadronic light-by-light contribution to the muon g-2

abstract: We have computed the hadronic light-by-light (LbL) contribution to the muon anomalous magnetic moment $a_\mu$ in the frame of Chiral Perturbation Theory with the inclusion of the lightest resonance multiplets as dynamical fields (R$\chi$T). It is essential to give a more accurate prediction of this hadronic contribution due to the future projects of J-Parc and FNAL on reducing the uncertainty in this observable. We, therefore, computed the pseudoscalar transition form factor and proposed the measurement of the $e^+e^-\to\mu^+\mu^-\pi^0$ cross section and dimuon invariant mass spectrum to determine more accurately its parameters. Then, we evaluated the pion exchange contribution to $a_\mu$, obtaining $(6.66\pm0.21)\cdot10^{-10}$. By comparing the pion exchange contribution and the pion-pole approximation to the corresponding transition form factor ($\pi$TFF) we recalled that the latter underestimates the complete $\pi$TFF by (15-20)\%. Then, we obtained the $\eta^({}'^)$ TFF, obtaining a total contribution of the lightest pseudoscalar exchanges of $(10.47\pm0.54)\cdot10^{-10}$, in agreement with previous results and with smaller error.


title: A counterexample to a strong variant of the Polynomial Freiman-Ruzsa conjecture in Euclidean space

abstract: The Polynomial Freiman-Ruzsa conjecture is one of the central open problems in additive combinatorics. If true, it would give tight quantitative bounds relating combinatorial and algebraic notions of approximate subgroups. In this note, we restrict our attention to subsets of Euclidean space. In this regime, the original conjecture considers approximate algebraic subgroups as the set of lattice points in a convex body. Green asked in 2007 whether this can be simplified to a generalized arithmetic progression, while not losing more than a polynomial factor in the underlying parameters. We give a negative answer to this question, based on a recent reverse Minkowski theorem combined with estimates for random lattices.


title: An Extended Cubic B-spline Finite Element Method for Solving Generalized Burgers-sher Equation

abstract: In this study collocation method based on the extended B-spline functions for the numerical solutions of the Generalized Burhers Fisher equation is set up. The approximate solution of the equation is constructed with the combination of the extended B-splines. Some initial boundary value problems are solved by the proposed method. The accuracy and validity and of the method is demonstrated by measuring the error between the numerical solutions and the analytical solutions, if exist.


title: Comment on "Motional Averaging of Nuclear Resonance in a Field Gradient"

abstract: In the Letter by Nanette N. Jarenwattananon and Louis-S. Bouchard [PRL 114, 197601 (2015)] an NMR experiment on gases in the presence of a magnetic-field gradient has been considered. As distinct from the traditional description of molecular self-diffusion, the authors calculate the decoherence of the signal taking into account the histories of molecular displacements. For this purpose the generalized Langevin equation (GLE) is applied. We show that the use of this equation is inappropriate. The calculations performed in the Letter are not correct and do not lead to the reported revised expression for line broadening that takes into account the autocorrelation effects in the diffusion process. The surprising temperature behavior of the observed NMR signal is thus not explained. In particular, the linewidth does not follow the power law f ~ T^(-1/2) at high temperatures. We give also Remarks on the Jarenwattananon and Bouchard Reply [Phys. Rev. Lett. 117, 249702 (2016)].


title: On the Capacity and Depth of Compact Surfaces

abstract: K. Borsuk in 1979, in the Topological Conference in Moscow, introduced the concept of capacity and depth of a compactum. In this paper, we compute the capacity and depth of compact surfaces. We show that the capacity and depth of every compact orientable surface of genus $g\geq 0$ is equal to $g+2$. Also, we prove that the capacity and depth of a compact non-orientable surface of genus $g>0$ is $[\frac{g}{2}]+2$.


title: Train Timetabling on double track and multiple station capacity railway with useful upper and lower bounds

abstract: Train scheduling is one of the significant issues in the railway industry in recent years since it has an important role in efficacy of railway infrastructure. In this paper, the timetabling problem of a multiple tracked railway network is discussed. More specifically, a general model is presented here in which a set of operational and safety requirements is considered. The model handles the trains overtaking in stations and considers the stations capacity. The objective function is to minimize the total travel time. Unfortunately, the problem is NP-hard and real size problems cannot be solved in an acceptable amount of time. In order to reduce the processing time, we presented some heuristic rules, which reduce the number of binary variables. These rules are based on problem's parameters such as travel time, dwell time and safety time of stations and try to remove the impracticable areas of the solution space. Furthermore, a Lagrangian Relaxation algorithm model is presented in order to find a lower-bound. Finally, comprehensive numerical experiments on the Tehran Metro case are reported. Results show the efficiency of the heuristic rules and also the Lagrangian Relaxation method in a way that for all analyzed problems the optimum value are obtained.


title: Effect of meson cloud on the jet nuclear modification factor in $pA$ collisions

abstract: We study the effect of the nucleon meson cloud on centrality dependence of the jet nuclear modification factor $R_{pA}$. We find that the meson-baryon Fock components may lead to a noticeable deviation of $R_{pA}$ from unity. Our results for $R_{pA}$ show the same tendency as that observed by ATLAS in $p+Pb$ collisions at $\sqrt{s}=5.02$ TeV. The meson cloud suppresses the central jet events and enhances the peripheral jet events. But quantitatively the effect is somewhat smaller than in the data.


title: Kondo Destruction in RKKY-Coupled Kondo Lattice and Multi-Impurity Systems

abstract: In a Kondo lattice, the spin exchange coupling between a local spin and the conduction electrons acquires nonlocal contributions due to conduction electron scattering from surrounding local spins and the subsequent RKKY interaction. It leads to a hitherto unrecognized interference of Kondo screening and the RKKY interaction beyond the Doniach scenario. We develop a renormalization group theory for the RKKY-modified Kondo vertex. The Kondo temperature, $T_K(y)$, is suppressed in a universal way, controlled by the antiferromagnetic RKKY coupling parameter $y$. Complete spin screening ceases to exist beyond a critical RKKY strength $y_c$ even in the absence of magnetic ordering. At this breakdown point, $T_K(y)$ remains nonzero and is not defined for larger RKKY couplings, $y>y_c$. The results are in quantitative agreement with STM spectroscopy experiments on tunable two-impurity Kondo systems. The possible implications for quantum critical scenarios in heavy-fermion systems are discussed.


title: A Primal-Dual Approximation Algorithm for Min-Sum Single-Machine Scheduling Problems

abstract: We consider the following single-machine scheduling problem, which is often denoted $1||\sum f_{j}$: we are given $n$ jobs to be scheduled on a single machine, where each job $j$ has an integral processing time $p_j$, and there is a nondecreasing, nonnegative cost function $f_j(C_{j})$ that specifies the cost of finishing $j$ at time $C_{j}$; the objective is to minimize $\sum_{j=1}^n f_j(C_j)$. Bansal \& Pruhs recently gave the first constant approximation algorithm with a performance guarantee of 16. We improve on this result by giving a primal-dual pseudo-polynomial-time algorithm based on the recently introduced knapsack-cover inequalities. The algorithm finds a schedule of cost at most four times the constructed dual solution. Although we show that this bound is tight for our algorithm, we leave open the question of whether the integrality gap of the LP is less than 4. Finally, we show how the technique can be adapted to yield, for any $\epsilon >0$, a $(4+\epsilon )$-approximation algorithm for this problem.


title: Three pearls of Bernoulli numbers

abstract: The Bernoulli numbers are fascinating and ubiquitous numbers, they occur in several domains of Mathematics like Number theory (FLT), Group theory, Calculus and even in Physics. Since Bernoulli's work, they are yet studied to understand their deep nature and particularly to find relationships between them. In this paper, we give, firstly, a short response to a problem stated, in 1971, by Carlitz and studied by many authors like Prodinger \cite{PRO}, the second pearl is an answer to a question raised, in 2008, by Tom Apostol .The third pearl is another proof of a relationship already given in 2011, by the authors.


title: Estimating covariance functions of multivariate skew-Gaussian random fields on the sphere

abstract: This paper considers a multivariate spatial random field, with each component having univariate marginal distributions of the skew-Gaussian type. We assume that the field is defined spatially on the unit sphere embedded in $\mathbb{R}^3$, allowing for modeling data available over large portions of planet Earth. This model admits explicit expressions for the marginal and cross covariances. However, the $n$-dimensional distributions of the field are difficult to evaluate, because it requires the sum of $2^n$ terms involving the cumulative and probability density functions of a $n$-dimensional Gaussian distribution. Since in this case inference based on the full likelihood is computationally unfeasible, we propose a composite likelihood approach based on pairs of spatial observations. This last being possible thanks to the fact that we have a closed form expression for the bivariate distribution. We illustrate the effectiveness of the method through simulation experiments and the analysis of a real data set of minimum and maximum temperatures.


title: Hydrodynamic-type systems describing 2-dimensional polynomially integrable geodesic flows

abstract: Starting from a homogeneous polynomial in momenta of arbitrary order we extract multi-component hydrodynamic-type systems which describe 2-dimensional geodesic flows admitting the initial polynomial as integral. All these hydrodynamic-type systems are semi-Hamiltonian, thus implying that they are integrable according to the generalized hodograph method. Moreover, they are integrable in a constructive sense as polynomial first integrals allow to construct generating equations of conservation laws. According to the multiplicity of the roots of the polynomial integral, we separate integrable particular cases.


title: Oblivious Sorting and Queues

abstract: We present a deterministic oblivious LIFO (Stack), FIFO, double-ended and double-ended priority queue as well as an oblivious mergesort and quicksort algorithm. Our techniques and ideas include concatenating queues end-to-end, size balancing of multiple arrays, several multi-level partitionings of an array. Our queues are the first to enable executions of pop and push operations without any change of the data structure (controlled by a parameter). This enables interesting applications in computing on encrypted data such as hiding confidential expressions. Mergesort becomes practical using our LIFO queue, ie. it improves prior work (STOC '14) by a factor of (more than) 1000 in terms of comparisons for all practically relevant queue sizes. We are the first to present double-ended (priority) and LIFO queues as well as oblivious quicksort which is asymptotically optimal. Aside from theortical analysis, we also provide an empirical evaluation of all queues.


title: Existence of martingale solutions and the incompressible limit for stochastic compressible flows on the whole space

abstract: We give an existence and asymptotic result for the so-called finite energy weak martingale solution of the compressible isentropic Navier--Stokes system driven by some random force in the whole spatial region. In particular, given a general nonlinear multiplicative noise, we establish the convergence to the incompressible system as the Mach number, representing the ratio between the average flow velocity and the speed of sound, approaches zero.


title: Obfuscation using Encryption

abstract: Protecting source code against reverse engineering and theft is an important problem. The goal is to carry out computations using confidential algorithms on an untrusted party while ensuring confidentiality of algorithms. This problem has been addressed for Boolean circuits known as `circuit privacy'. Circuits corresponding to real-world programs are impractical. Well-known obfuscation techniques are highly practicable, but provide only limited security, e.g., no piracy protection. In this work, we modify source code yielding programs with adjustable performance and security guarantees ranging from indistinguishability obfuscators to (non-secure) ordinary obfuscation. The idea is to artificially generate `misleading' statements. Their results are combined with the outcome of a confidential statement using encrypted \emph{selector variables}. Thus, an attacker must `guess' the encrypted selector variables to disguise the confidential source code. We evaluated our method using more than ten programmers as well as pattern mining across open source code repositories to gain insights of (micro-)coding patterns that are relevant for generating misleading statements. The evaluation reveals that our approach is effective in that it successfully preserves source code confidentiality.


title: Location, identification, and representability of monotone operators in locally convex spaces

abstract: In this paper we study, in the relaxed context of locally convex spaces, intrinsic properties of monotone operators needed for the sum conjecture for maximal monotone operators to hold under classical interiority-type domain constraints.


title: Dual Moments and Risk Attitudes

abstract: In decision under risk, the primal moments of mean and variance play a central role to define the local index of absolute risk aversion. In this paper, we show that in canonical non-EU models dual moments have to be used instead of, or on par with, their primal counterparts to obtain an equivalent index of absolute risk aversion.


title: The evolution of the radial gradient of Oxygen abundance in spiral galaxies

abstract: The aim of this work is to present our new series of chemical evolution models computed for spiral and low mass galaxies of different total masses and star formation efficiencies. We analyze the results of models, in particular the evolution of the radial gradient of oxygen abundance. Furthermore, we study the role of the infall rate and of the star formation history on the variations of this radial gradient. The relations between the O/H radial gradient and other spiral galaxies characteristics as the size or the stellar mass are also shown. We find that the radial gradient is mainly a scale effect which basically does not change with the redshift (or time) if it is measured within the optical radius. Moreover, when it is measured as a function of a normalized radius, show a similar value for all galaxies masses, showing a correlation with a dispersion around an average value which is due to the differences star formation efficiencies, in agreement with the idea of an universal O/H radial gradient


title: An Empirical Study of ADMM for Nonconvex Problems

abstract: The alternating direction method of multipliers (ADMM) is a common optimization tool for solving constrained and non-differentiable problems. We provide an empirical study of the practical performance of ADMM on several nonconvex applications, including l0 regularized linear regression, l0 regularized image denoising, phase retrieval, and eigenvector computation. Our experiments suggest that ADMM performs well on a broad class of non-convex problems. Moreover, recently proposed adaptive ADMM methods, which automatically tune penalty parameters as the method runs, can improve algorithm efficiency and solution quality compared to ADMM with a non-tuned penalty.


title: Non-negative Factorization of the Occurrence Tensor from Financial Contracts

abstract: We propose an algorithm for the non-negative factorization of an occurrence tensor built from heterogeneous networks. We use l0 norm to model sparse errors over discrete values (occurrences), and use decomposed factors to model the embedded groups of nodes. An efficient splitting method is developed to optimize the nonconvex and nonsmooth objective. We study both synthetic problems and a new dataset built from financial documents, resMBS.


title: Harmonic Maass forms associated to real quadratic fields

abstract: In this paper, we explicitly construct harmonic Maass forms that map to the weight one theta series associated by Hecke to odd ray class group characters of real quadratic fields. From this construction, we give precise arithmetic information contained in the Fourier coefficients of the holomorphic part of the harmonic Maass form, establishing the main part of a conjecture of the second author.


title: Fractional constitutive equation (FACE) for non-Newtonian fluid flow: Theoretical description

abstract: Non-Newtonian fluid flow might be driven by spatially nonlocal velocity, the dynamics of which can be described by promising fractional derivative models. This short communication proposes a space FrActional-order Constitutive Equation (FACE) that links viscous shear stress with the velocity gradient, and then interprets physical properties of non-Newtonian fluids for steady pipe flow. Results show that the generalized FACE model contains previous non-Newtonian fluid flow models as end-members by simply adjusting the order of the fractional index, and a preliminary test shows that the FACE model conveniently captures the observed growth of shear stress for various velocity gradients. Further analysis of the velocity profile, frictional head loss, and Reynolds number using the FACE model also leads to analytical tools and criterion that can significantly extend standard models in quantifying the complex dynamics of non-Newtonian fluid flow with a wide range of spatially nonlocal velocities.


title: FOCA: A Methodology for Ontology Evaluation

abstract: Modeling an ontology is a hard and time-consuming task. Although methodologies are useful for ontologists to create good ontologies, they do not help with the task of evaluating the quality of the ontology to be reused. For these reasons, it is imperative to evaluate the quality of the ontology after constructing it or before reusing it. Few studies usually present only a set of criteria and questions, but no guidelines to evaluate the ontology. The effort to evaluate an ontology is very high as there is a huge dependence on the evaluator's expertise to understand the criteria and questions in depth. Moreover, the evaluation is still very subjective. This study presents a novel methodology for ontology evaluation, taking into account three fundamental principles: i) it is based on the Goal, Question, Metric approach for empirical evaluation; ii) the goals of the methodologies are based on the roles of knowledge representations combined with specific evaluation criteria; iii) each ontology is evaluated according to the type of ontology. The methodology was empirically evaluated using different ontologists and ontologies of the same domain. The main contributions of this study are: i) defining a step-by-step approach to evaluate the quality of an ontology; ii) proposing an evaluation based on the roles of knowledge representations; iii) the explicit difference of the evaluation according to the type of the ontology iii) a questionnaire to evaluate the ontologies; iv) a statistical model that automatically calculates the quality of the ontologies.


title: Resonant Raman imaging of MoS2-substrate interaction

abstract: We report a study of long-range MoS2-substrate interaction using resonant Raman imaging. We observed a strong thickness-dependent peak shift of a Raman-forbidden mode that can be used as a new method of determining the thickness of multilayered MoS2 flakes. In addition, dependence of the Raman scattering intensity on thickness is explained by the interference enhancement theory. Finally, the resonant Raman spectra on different substrates are analysed.


title: Stability of periodic orbits in no-slip billiards

abstract: Rigid bodies collision maps in dimension two, under a natural set of physical requirements, can be classified into two types: the standard specular reflection map and a second which we call, after Broomhead and Gutkin, no-slip. This leads to the study of no-slip billiards--planar billiard systems in which the moving particle is a disc (with rotationally symmetric mass distribution) whose translational and rotational velocities can both change at each collision with the boundary of the billiard domain. In this paper we greatly extend previous results on boundedness of orbits (Broomhead and Gutkin) and linear stability of periodic orbits for a Sinai-type billiard (Wojtkowski) for no-slip billiards. We show among other facts that: (i) for billiard domains in the plane having piecewise smooth boundary and at least one corner of inner angle less than $\pi$, no-slip billiard dynamics will always contain elliptic period-$2$ orbits; (ii) polygonal no-slip billiards always admit small invariant open sets and thus cannot be ergodic with respect to the canonical invariant billiard measure; (iii) the no-slip version of a Sinai billiard must contain linearly stable periodic orbits of period $2$ and, more generally, we provide a curvature threshold at which a commonly occurring period-$2$ orbit shifts from being hyperbolic to being elliptic; (iv) finally, we make a number of observations concerning periodic orbits in a class of polygonal billiards.


title: The Majorana STM as a perfect detector of odd-frequency superconductivity

abstract: We propose a novel scanning tunneling microscope (STM) device in which the tunneling tip is formed by a Majorana bound state (MBS). This peculiar bound state exists at the boundary of a one-dimensional topological superconductor. Since the MBS has to be effectively spinless and local, we argue that it is the smallest unit that shows itself odd-frequency superconducting pairing. Odd-frequency superconductivity is characterized by an anomalous Green function which is an odd function of the time arguments of the two electrons forming the Cooper pair. Interestingly, our Majorana STM can be used as the perfect detector of odd-frequency superconductivity. The reason is that a supercurrent between the Majorana STM and any other superconductor can only flow if the latter system exhibits itself odd-frequency pairing. To illustrate our general idea, we consider the tunneling problem of the Majorana STM coupled to a quantum dot in vicinity to a conventional superconductor. In such a (superconducting) quantum dot, the effective pairing can be tuned from even- to odd-frequency behavior if an external magnetic field is applied to it.


title: Limited Feedback in MISO Systems with Finite-Bit ADCs

abstract: We analyze limited feedback in systems where a multiple-antenna transmitter sends signals to single-antenna receivers with finite-bit ADCs. If channel state information (CSI) is not available with high resolution at the transmitter and the precoding is not well designed, the inter-user interference is a big decoding challenge for receivers with low-resolution quantization. In this paper, we derive achievable rates with finite-bit ADCs and finite-bit CSI feedback. The performance loss compared to the case with perfect CSI is then analyzed. The results show that the number of bits per feedback should increase linearly with the ADC resolution to restrict the loss.


title: A large arboreal Galois representation for a cubic postcritically finite polynomial

abstract: We give a complete description of the arboreal Galois representation of a certain postcritically finite cubic polynomial over a large class of number fields and for a large class of basepoints. This is the first such example that is not conjugate to a power map, Chebyshev polynomial, or Latt\`es map. The associated Galois action on an infinite ternary rooted tree has Hausdorff dimension bounded strictly between that of the infinite wreath product of cyclic groups and that of the infinite wreath product of symmetric groups. We deduce a zero-density result for prime divisors in an orbit under this polynomial. We also obtain a zero-density result for the set of places of convergence of Newton's method for a certain cubic polynomial, thus resolving the first nontrivial case of a conjecture of Faber and Voloch.


title: Accurate representation of formation energies of crystalline alloys with many components

abstract: In this paper I propose a new model for representing the formation energies of multicomponent crystalline alloys as a function of atom types. In the cases when displacements of atoms from their equilibrium positions are not large, the proposed method has a similar accuracy as the state-of-the-art cluster expansion method, and a better accuracy when the fitting dataset size is small. The proposed model has only two tunable parameters - one for the interaction range and one for the interaction complexity.


title: Information Theory of Molecular Communication: Directions and Challenges

abstract: Molecular Communication (MC) is a communication strategy that uses molecules as carriers of information, and is widely used by biological cells. As an interdisciplinary topic, it has been studied by biologists, communication theorists and a growing number of information theorists. This paper aims to specifically bring MC to the attention of information theorists. To do this, we first highlight the unique mathematical challenges of studying the capacity of molecular channels. Addressing these problems require use of known, or development of new mathematical tools. Toward this goal, we review a subjective selection of the existing literature on information theoretic aspect of molecular communication. The emphasis here is on the mathematical techniques used, rather than on the setup or modeling of a specific paper. Finally, as an example, we propose a concrete information theoretic problem that was motivated by our study of molecular communication.


title: An Energy-Efficient VCO-Based Matrix Multiplier Block to Support On-Chip Image Analysis

abstract: Images typically are represented as uniformly sampled data in the form of matrix of pixels/voxels. Therefore, matrix multiply-and-accumulate (MAC) forms the core of most state-of-the-art image analysis algorithms. While digital implementation of MAC has generally been the preferred approach, high power consumption is an impediment to adopting it for medical image analysis. In this work, we present a time-domain signal processing architecture which performs MAC operations with 7bit accuracy while consuming 400X lower energy than digital implementation. The proposed architecture performs analog computation using mostly digital circuits and is suitable for scaled CMOS technologies. The proposed time-domain MAC architecture is expected to play a central role in empowering the advancement of various on-chip image analysis operations.


title: Community Detection Algorithm Evaluation using Size and Hashtags

abstract: Understanding community structure in social media is critical due to its broad applications such as friend recommendations, link predictions and collaborative filtering. However, there is no widely accepted definition of community in literature. Existing work use structure related metrics such as modularity and function related metrics such as ground truth to measure the performance of community detection algorithms, while ignoring an important metric, size of the community. [1] suggests that the size of community with strong ties in social media should be limited to 150. As we discovered in this paper, the majority of the communities obtained by many popular community detection algorithms are either very small or very large. Too small communities don't have practical value and too large communities contain weak connections therefore not stable. In this paper, we compare various community detection algorithms considering the following metrics: size of the communities, coverage of the communities, extended modularity, triangle participation ratio, and user interest in the same community. We also propose a simple clique based algorithm for community detection as a baseline for the comparison. Experimental results show that both our proposed algorithm and the well-accepted disjoint algorithm InfoMap perform well in all the metrics.


title: Quantum dynamical entropy, chaotic unitaries and complex Hadamard matrices

abstract: We introduce two information-theoretical invariants for the projective unitary group acting on a finite-dimensional complex Hilbert space: PVM- and POVM-dynamical (quantum) entropies. They quantify the randomness of the successive quantum measurement results in the case where the evolution of the system between each two consecutive measurements is described by a given unitary operator. We study the class of chaotic unitaries, i.e., the ones of maximal entropy or, equivalently, such that they can be represented by suitably rescaled complex Hadamard matrices in some orthonormal bases. We provide necessary conditions for a unitary operator to be chaotic, which become also sufficient for qubits and qutrits. These conditions are expressed in terms of the relation between the trace and the determinant of the operator. We also compute the volume of the set of chaotic unitaries in dimensions two and three, and the average PVM-dynamical entropy over the unitary group in dimension two. We prove that this mean value behaves as the logarithm of the dimension of the Hilbert space, which implies that the probability that the dynamical entropy of a unitary is almost as large as possible approaches unity as the dimension tends to infinity.


title: Technical Report: A Generalized Matching Pursuit Approach for Graph-Structured Sparsity

abstract: Sparsity-constrained optimization is an important and challenging problem that has wide applicability in data mining, machine learning, and statistics. In this paper, we focus on sparsity-constrained optimization in cases where the cost function is a general nonlinear function and, in particular, the sparsity constraint is defined by a graph-structured sparsity model. Existing methods explore this problem in the context of sparse estimation in linear models. To the best of our knowledge, this is the first work to present an efficient approximation algorithm, namely, Graph-structured Matching Pursuit (Graph-Mp), to optimize a general nonlinear function subject to graph-structured constraints. We prove that our algorithm enjoys the strong guarantees analogous to those designed for linear models in terms of convergence rate and approximation accuracy. As a case study, we specialize Graph-Mp to optimize a number of well-known graph scan statistic models for the connected subgraph detection task, and empirical evidence demonstrates that our general algorithm performs superior over state-of-the-art methods that are designed specifically for the task of connected subgraph detection.


title: Multiple Instance Learning: A Survey of Problem Characteristics and Applications

abstract: Multiple instance learning (MIL) is a form of weakly supervised learning where training instances are arranged in sets, called bags, and a label is provided for the entire bag. This formulation is gaining interest because it naturally fits various problems and allows to leverage weakly labeled data. Consequently, it has been used in diverse application fields such as computer vision and document classification. However, learning from bags raises important challenges that are unique to MIL. This paper provides a comprehensive survey of the characteristics which define and differentiate the types of MIL problems. Until now, these problem characteristics have not been formally identified and described. As a result, the variations in performance of MIL algorithms from one data set to another are difficult to explain. In this paper, MIL problem characteristics are grouped into four broad categories: the composition of the bags, the types of data distribution, the ambiguity of instance labels, and the task to be performed. Methods specialized to address each category are reviewed. Then, the extent to which these characteristics manifest themselves in key MIL application areas are described. Finally, experiments are conducted to compare the performance of 16 state-of-the-art MIL methods on selected problem characteristics. This paper provides insight on how the problem characteristics affect MIL algorithms, recommendations for future benchmarking and promising avenues for research.


title: Palomar Transient Factory and RR Lyrae: Metallicity-Light Curve Relation Based on ab-Type RR Lyrae in Kepler Field

abstract: The wide-field synoptic sky surveys, known as the Palomar Transient Factory (PTF) and the intermediate Palomar Transient Factory (iPTF), will accumulate a large number of known and new RR Lyrae. These RR Lyrae are good tracers to study the substructure of the Galactic halo if their distance, metallicity, and galactocentric velocity can be measured. Candidates of halo RR Lyrae can be identified from their distance and metallicity before requesting spectroscopic observations for confirmation. This is because both quantities can be obtained via their photometric light curves, because the absolute V-band magnitude for RR Lyrae is correlated with metallicity, and the metallicity can be estimated using a metallicity-light curve relation. To fully utilize the PTF and iPTF light-curve data in related future work, it is necessary to derive the metallicity-light curve relation in the native PTF/iPTF R-band photometric system. In this work, we derived such a relation using the known ab-type RR Lyrae located in the Kepler field, and it is found to be $[Fe/H]_{PTF} = -4.089 - 7.346 P + 1.280 \phi_{31}$ (where $P$ is pulsational period and $\phi_{31}$ is one of the Fourier parameters describing the shape of the light curve), with a dispersion of 0.118 dex. We tested our metallicity-light curve relation with new spectroscopic observations of a few RR Lyrae in the Kepler field, as well as several data sets available in the literature. Our tests demonstrated that the derived metallicity-light curve relation could be used to estimate metallicities for the majority of the RR Lyrae, which are in agreement with the published values.


title: Nilpotent orbit theorem in $p$-adic Hodge theory

abstract: We state and prove three orbit theorems on the period domains for the $p$-adic Hodge structure analogous to the complex case. We shall consider the variation of de Rham (resp. \'etale) cohomology in a family of projective varieties $f:\mathfrak{X} \to S$ defined over a p-adic field. First, we show that any nilpotent orbit in the period domain of p-adic Hodge structures converges to a semistable point (filtration) in the period domain of the p-adic Hodge structure. Furthermore, the nilpotent orbits of the limit point are asymptotic to the twisted period map [Theorem \ref{thm:nilpotent-orbit}]. The orbit theorems come with some estimates of the distance between the nilpotent orbit and the twisted period map. The distance estimate in the p-adic nilpotent orbit theorem is given concerning the non-archimedean metric and is based on the p-adic Fourier analysis of Amice-Schneider. The result is analogous to the orbit theorems of W. Schmid [\cite{Sch}-1973] on complex Hodge structures. Our proof is based on a \textit{Geometric Invariant Theory} (GIT) criterion for semi-stability (Kempf-Ness theorem) and estimates from the (Amice-Schneider) p-adic Fourier theory. We also state the $SL_2$-orbit theorem in the p-adic case, [Theorem \ref{th:homomorphism}]. Finally, we explain how the nilpotent orbit theorem should be modified and stated for a variation of the mixed Hodge structure [Theorem \ref{thm:mixed-orbit}].}


title: On the question of genericity of hyperbolic knots

abstract: A well-known conjecture in knot theory says that the percentage of hyperbolic knots amongst all of the prime knots of $n$ or fewer crossings approaches $100$ as $n$ approaches infinity. In this paper, it is proved that this conjecture contradicts several other plausible conjectures, including the 120-year-old conjecture on additivity of the crossing number of knots under connected sum and the conjecture that the crossing number of a satellite knot is not less than that of its companion.


title: Flexible and scalable particle-in-cell methods for massively parallel computations

abstract: Particle-in-cell methods couple mesh-based methods for the solution of continuum mechanics problems, with the ability to advect and evolve particles. They have a long history and many applications in scientific computing. However, they have most often only been implemented for either sequential codes, or parallel codes with static meshes that are statically partitioned. In contrast, many mesh-based codes today use adaptively changing, dynamically partitioned meshes, and can scale to thousands or tens of thousands of processors. Consequently, there is a need to revisit the data structures and algorithms necessary to use particle methods with modern, mesh-based methods. Here we review commonly encountered requirements of particle-in-cell methods, and describe efficient ways to implement them in the context of large-scale parallel finite-element codes that use dynamically changing meshes. We also provide practical experience for how to address bottlenecks that impede the efficient implementation of these algorithms and demonstrate with numerical tests both that our algorithms can be implemented with optimal complexity and that they are suitable for very large-scale, practical applications. We provide a reference implementation in ASPECT, an open source code for geodynamic mantle-convection simulations built on the deal.II library.


title: The Lackadaisical Quantum Walker is NOT Lazy at all

abstract: In this paper, we study the properties of lackadaisical quantum walks on a line. This model is first proposed in~\cite{wong2015grover} as a quantum analogue of lazy random walks where each vertex is attached $\tau$ self-loops. We derive an analytic expression for the localization probability of the walker at the origin after infinite steps, and obtain the peak velocities of the walker. We also calculate rigorously the wave function of the walker starting from the origin and obtain a long time approximation for the entire probability density function. As an application of the density function, we prove that lackadaisical quantum walks spread ballistically for arbitrary $\tau$, and give an analytic solution for the variance of the walker's probability distribution.


title: Rangzen: Anonymously Getting the Word Out in a Blackout

abstract: In recent years governments have shown themselves willing to impose blackouts to shut off key communication infrastructure during times of civil strife, and to surveil citizen communications whenever possible. However, it is exactly during such strife that citizens need reliable and anonymous communications the most. In this paper, we present Rangzen, a system for anonymous broadcast messaging during network blackouts. Rangzen is distinctive in both aim and design. Our aim is to provide an anonymous, one-to-many messaging layer that requires only users' smartphones and can withstand network-level attacks. Our design is a delay-tolerant mesh network which deprioritizes adversarial messages by means of a social graph while preserving user anonymity. We built a complete implementation that runs on Android smartphones, present benchmarks of its performance and battery usage, and present simulation results suggesting Rangzen's efficacy at scale.


title: On Jacobian group and complexity of the generalized Petersen graph GP(n,k) through Chebyshev polynomials

abstract: In the present paper we find a simple algorithm for counting Jacobian group of the generalized Petersen graph GP(n,k). Also, we obtain a closed formula for the number of spanning trees of this graph in terms of Chebyshev polynomials.


title: A probabilistic graphical model approach in 30 m land cover mapping with multiple data sources

abstract: There is a trend to acquire high accuracy land-cover maps using multi-source classification methods, most of which are based on data fusion, especially pixel- or feature-level fusions. A probabilistic graphical model (PGM) approach is proposed in this research for 30 m resolution land-cover mapping with multi-temporal Landsat and MODerate Resolution Imaging Spectroradiometer (MODIS) data. Independent classifiers were applied to two single-date Landsat 8 scenes and the MODIS time-series data, respectively, for probability estimation. A PGM was created for each pixel in Landsat 8 data. Conditional probability distributions were computed based on data quality and reliability by using information selectively. Using the administrative territory of Beijing City (Area-1) and a coastal region of Shandong province, China (Area-2) as study areas, multiple land-cover maps were generated for comparison. Quantitative results show the effectiveness of the proposed method. Overall accuracies promoted from 74.0% (maps acquired from single-temporal Landsat images) to 81.8% (output of the PGM) for Area-1. Improvements can also be seen when using MODIS data and only a single-temporal Landsat image as input (overall accuracy: 78.4% versus 74.0% for Area-1, and 86.8% versus 83.0% for Area-2). Information from MODIS data did not help much when the PGM was applied to cloud free regions of. One of the advantages of the proposed method is that it can be applied where multi-temporal data cannot be simply stacked as a multi-layered image.


title: Beta-plane turbulence above monoscale topography

abstract: Using a one-layer QG model, we study the effect of random monoscale topography on forced beta-plane turbulence. The forcing is a uniform steady wind stress that produces both a uniform large-scale zonal flow $U(t)$ and smaller-scale macroturbulence (both standing and transient eddies). The flow $U(t)$ is retarded by Ekman drag and by the domain-averaged topographic form stress produced by the eddies. The topographic form stress typically balances most of the applied wind stress, while the Ekman drag provides all of the energy dissipation required to balance the wind work. A collection of statistically equilibrated solutions delineates the main flow regimes and the dependence of the time-mean $U$ on the problem parameters and the statistical properties of the topography. If $\beta$ is smaller than the topographic PV gradient then the flow consists of stagnant pools attached to pockets of closed geostrophic contours. The stagnant dead zones are bordered by jets and the flow through the domain is concentrated into a narrow channel of open geostrophic contours. If $\beta$ is comparable to, or larger than, the topographic PV gradient then all geostrophic contours are open and the flow is uniformly distributed throughout the domain. In this case there is an "eddy saturation" regime in which $U$ is insensitive to changes in the wind stress. We show that eddy saturation requires strong transient eddies that act as PV diffusion. This PV diffusion does not alter the energy of the standing eddies, but it does increase the topographic form stress by enhancing the correlation between topographic slope and the standing-eddy pressure field. Last, using bounds based on the energy and enstrophy we show that as the wind stress increases the flow transitions from a regime in which form stress balances the wind stress to a regime in which the form stress is very small and large transport ensues.


title: Sample complexity of the distinct elements problem

abstract: We consider the distinct elements problem, where the goal is to estimate the number of distinct colors in an urn containing $ k $ balls based on $n$ samples drawn with replacements. Based on discrete polynomial approximation and interpolation, we propose an estimator with additive error guarantee that achieves the optimal sample complexity within $O(\log\log k)$ factors, and in fact within constant factors for most cases. The estimator can be computed in $O(n)$ time for an accurate estimation. The result also applies to sampling without replacement provided the sample size is a vanishing fraction of the urn size. One of the key auxiliary results is a sharp bound on the minimum singular values of a real rectangular Vandermonde matrix, which might be of independent interest.


title: Revisiting Jackson-Hunt calculations: Unified theoretical analysis for generic multi-phase growth in a multi-component system

abstract: A straight-forward extension of the Jackson-Hunt theory for directionally solidifying multi-phase growth where the number of components exceeds the number of solid phases becomes difficult on account of the absence of the required number of equations to determine the boundary layer compositions ahead of the interface. In this paper, we therefore revisit the Jackson-Hunt(JH) type calculations for any given situation of multi-phase growth in a multi-component system and self-consistently derive the variations of the compositions of the solid phases as well as their volume fractions, which grow such that the composite solid-liquid interface is isothermal. This allows us to unify the (JH) calculation schemes for both in-variant as well as multi-variant eutectic reactions. The derived analytical expressions are then utilized to study the effect of dissimilar solute diffusivities and interfacial energies on the undercoolings and the solidified fractions. We also perform phase field simulations to confirm our theoretical predictions and find a good agreement between our analytical calculations and model predictions for model symmetric alloys as well as for a particular Ni-Al-Zr alloy.


title: A chain transitive accessible partially hyperbolic diffeomorphism

abstract: In this paper, we construct a partially hyperbolic skew-product diffeomorphism $f$ on $\mathbb{T}^3$, such that $f$ is accessible and chain transitive, but not transitive.


title: Collider signals of $W'$ and $Z'$ bosons in the gauge-Higgs unification

abstract: In the $SO(5)\times U(1)$ gauge-Higgs unification (GHU), Kaluza-Klein (KK) excited states of charged and neutral vector bosons, $W^{(1)}$, $W_R^{(1)}$, $Z^{(1)}$, $\gamma^{(1)}$ and $Z_R^{(1)}$, can be observed as $W'$ and $Z'$ signals in collider experiments. In this paper we evaluate the decay rates of the $W'$ and $Z'$, and $s$-channel cross sections mediated by $W'$ and $Z'$ bosons with final states involving the standard model (SM) fermion pair ($\ell \nu$, $\ell\bar{\ell}$, $q\bar{q}'$), $WH$, $ZH$, $WW$ and $WZ$. $W'$ and $Z'$ resonances appear around 6.0 TeV (8.5 TeV) for $\theta_H=0.115$ (0.0737) where $\theta_H$ is the Aharonov-Bohm phase in the fifth dimension in GHU. For decay rates we find $\Gamma(W'\to WH) \simeq \Gamma(W'\to WZ)$ ($W'=W^{(1)}$, $W_R^{(1)}$), $\Gamma(W^{(1)}\to WH,WZ) \sim\Gamma(W_R^{(1)}\to WH,WZ)$, $\Gamma(Z^{(1)}\to ZH) \simeq\sum_{Z'=Z^{(1)},\gamma^{(1)}}\Gamma(Z'\to WW)$, and $\Gamma(Z_R^{(1)}\to ZH) \simeq \Gamma(Z_R^{(1)}\to WW)$. $W'$ and $Z'$ signals of GHU can be best found at the LHC experiment in the processes $pp\to W' (Z')+X$ followed by $W'\to t\bar{b}$, $WH$, and $Z'\to e^+ e^-$, $\mu^+\mu^-, ZH$. For the lighter $Z'$ ($\theta_H = 0.115$) case, with forthcoming 30 fb$^{-1}$ data of the 13 TeV LHC experiment we expect about ten $\mu^+\mu^-$ events for the invariant mass range 3 to 7 TeV, though the number of the events becomes much smaller when $\theta_H=0.0737$. In the process with $WZ$ in the final state, it is confirmed that the unitarity is preserved, provided that all KK excited intermediate states are taken into account. Deviation of the $WWZ$ coupling from the SM is very tiny. Exotic partners $t_T^{(1)}$ and $b_Y^{(1)}$ of the top and bottom quarks with electric charge $+5/3$ and $-4/3$ have mass $M_{t_T^{(1)},b_Y^{(1)}}\simeq 4-5$ TeV, becoming the lightest non-SM particles in GHU which can be singly produced in collider experiments.


title: Thermodynamics of a time dependent and dissipative oval billiard: a heat transfer and billiard approach

abstract: We study some statistical properties for the behavior of the average squared velocity -- hence the temperature -- for an ensemble of classical particles moving in a billiard whose boundary is time dependent. We assume the collisions of the particles with the boundary of the billiard are inelastic leading the average squared velocity to reach a steady state dynamics for large enough time. The description of the stationary state is made by using two different approaches: (i) heat transfer motivated by the Fourier law and, (ii) billiard dynamics using either numerical simulations and theoretical description.


title: Curvature Induced Topological Defects of $p$-wave Superfluid on a Sphere

abstract: We study the ground state of spinless fermions living on a sphere across $p$-wave Feschbach resonances. By construsting a microscopic model of fermions on a general curved surface, we show that the Guassian curvature induces an emergent magnetic field coupled to the $p\pm ip$ order parameters. In the case of a sphere, the magnetic field corresponds to a Dirac monopole field, which causes topological defects in the superfluid ground state. Using the BCS mean field theory, we calculate its many-body ground state self consistently and give the phase diagram. The ground state may exhibit two types of topological defects, two voritces on the south and north pole or a domain wall which separates $p_\theta+ ip_\phi$ and $p_\theta-ip_\phi$ superfluids.


title: Detecting planet pairs in mean motion resonances via astrometry method

abstract: GAIA leads us to step into a new era with a high astrometry precision of 10 uas. Under such a precision, astrometry will play important roles in detecting and characterizing exoplanets. Specially, we can identify planet pairs in mean motion resonances(MMRs) via astrometry, which constrains the formation and evolution of planetary systems. In accordance with observations, we consider two Jupiters or two super-Earths systems in 1:2, 2:3 and 3:4 MMRs. Our simulations show the false alarm probabilities(FAPs) of a third planet are extremely small while the real two planets can be good fitted with signal-to-noise ratio(SNR)> 3. The probability of reconstructing a resonant system is related with the eccentricities and resonance intensity. Generally, when SNR >= 10, if eccentricities of both planets are larger than 0.01 and the resonance is quite strong, the probabilities to reconstruct the planet pair in MMRs >= 80%. Jupiter pairs in MMRs are reconstructed more easily than super-Earth pairs with similar SNR when we consider the dynamical stability. FAPs are also calculated when we detect planet pairs in or near MMRs. FAPs for 1:2 MMR are largest, i.e., FAPs > 15% when SNR <= 10. Extrapolating from the Kepler planet pairs near MMRs and assuming SNR to be 3, we will discover and reconstruct a few tens of Jupiter pairs and hundreds of super-Earth pairs in 2:3 and 1:2 MMRs within 30 pc. We also compare the differences between even and uneven data cadence and find that planets are better measured with more uniform phase coverage.


title: A Novel Motion Detection Method Resistant to Severe Illumination Changes

abstract: Recently, there has been a considerable attention given to the motion detection problem due to the explosive growth of its applications in video analysis and surveillance systems. While the previous approaches can produce good results, an accurate detection of motion remains a challenging task due to the difficulties raised by illumination variations, occlusion, camouflage, burst physical motion, dynamic texture, and environmental changes such as those on climate changes, sunlight changes during a day, etc. In this paper, we propose a novel per-pixel motion descriptor for both motion detection and dynamic texture segmentation which outperforms the current methods in the literature particularly in severe scenarios. The proposed descriptor is based on two complementary three-dimensional-discrete wavelet transform (3D-DWT) and three-dimensional wavelet leader. In this approach, a feature vector is extracted for each pixel by applying a novel three dimensional wavelet-based motion descriptor. Then, the extracted features are clustered by a clustering method such as well-known k-means algorithm or Gaussian Mixture Model (GMM). The experimental results demonstrate the effectiveness of our proposed method compared to the other motion detection approaches from the literature. The application of the proposed method and additional experimental results for the different datasets are available at (http://dspl.ce.sharif.edu/motiondetector.html).


title: A survey of configurations

abstract: To any finite ordered subset and any finite partition of a group a set of tuples of positive integers, named as configurations, is associated that describes the group's behavior. The present paper provides an exposition of this notion and related topics, including amenability, paradoxical decomposition and the group theoretical properties.


title: Locations of Roberge-Weiss transition endpoints in lattice QCD with $N_f=2$ improved KS quarks

abstract: Result on the locations of the tricritical points of $N_f=2$ lattice QCD with imaginary chemical potential is presented. Simulations are carried out with Symanzik improved gauge action and Asqtad fermion action. With imaginary chemical potential $i\mu_I=i\pi T$, previous studies show that the Roberge-Weiss (RW) transition endpoints are triple points at both large and small quark masses, and second order transition points at intermediate quark masses. The triple and second order endpoints are separated by two tricritical ones. Our simulations are carried out at 7 values of quark mass $am$ ranging from 0.024 to 0.070 on lattice volume $12^3\times 4, 16^3\times 4,\, 20^3\times4$. The susceptibility and Binder cumulant of the imaginary part of Polyakov loop are employed to determine the nature of RW transition endpoints. The simulations suggest that the two tricritical points are within the range $0.024-0.026$ and $0.040-0.050$, respectively.


title: Spatial Charge Inhomogeneity and Defect States in Topological Dirac Semimetal Thin Films

abstract: The close approach of the Fermi energy EF of a Dirac semimetal to the Dirac point ED uncovers new physics such as velocity renormalization,1,2,3 and the Dirac plasma 4,5 at |EF -ED| < kBT, where kBT is the thermal energy. In graphene, substrate disorder drives fluctuations in EF. Three-dimensional topological Dirac semimetals (TDS)6,7 obviate the substrate, and should show reduced EF fluctuations due to better metallic screening and higher dielectric constants. Here we map the potential fluctuations in TDS Na3Bi using a scanning tunneling microscope. The rms potential fluctuations are significantly smaller than room temperature ({\Delta}EF,rms = 4-6 meV = 40-70 K) and comparable to the highest quality graphene on h-BN;8 far smaller than graphene on SiO2,9,10 or the Dirac surface state of a topological insulator.11 Surface Na vacancies produce a novel resonance close to the Dirac point with surprisingly large spatial extent and provides a unique way to tune the surface density of states in a TDS thin-film material.


title: Superconvergence analysis of DG-FEM based on the polynomial preserving recovery for Helmholtz equation with high wave number

abstract: We study superconvergence property of the linear discontinuous Galerkin finite element method with the polynomial preserving recovery (PPR) and Richardson extrapolation for the two dimensional Helmholtz equation. The error estimate with explicit dependence on the wave number $k$, the penalty parameter $\mu$ and the mesh condition parameter $\alpha$ is derived. First, we prove that under the assumption $k(kh)^2\leq C_0$ ($h$ is the mesh size) and certain mesh condition, the estimate between the finite element solution and the linear interpolation of the exact solution is superconvergent under the $\norme{\cdot}$-seminorm. Second, we prove a superconvergence result for the recovered gradient by PPR. Furthermore, we estimate the error between the finite element gradient and recovered gradient, which motivate us to define the a posteriori error estimator. Finally, Some numerical examples are provided to confirm the theoretical results of superconvergence analysis. All theoretical findings are verified by numerical tests.


title: Metallicities in long GRB host galaxies at z$<$0.5 calculated by the detailed modelling of optical and infrared line ratios

abstract: We revisited the line spectra emitted from long GRB (LGRB) host galaxies at z<0.5 in order to calculate by the detailed modelling of the line ratios the physical conditions and relative abundances in LGRB hosts in this redshift range. We have found lower metallicities than in LGRB hosts at higher z. New results about metallicities and physical conditions in the different regions throughout the LGRB 980425 host at z=0.0085 are presented. In particular, we have found that the effective starburst temperature in the supernova (SN) region is the highest throughout the host galaxy. The low ionization parameter reveals that the radiation source is far or somehow prevented from reaching the emitting gas in the SN region. The models constrained by a few oxygen, nitrogen and sulphur line ratios to Hb in LGRB 980425 host satisfactorily reproduce the HeII/Hb and [ArIII]/Hb line ratios. The modelling of the observed [SIV]10.51\mu/[SIII]18.71\mum and [NeIII]10.6\mum/[NeII]12.81\mum line ratios from LGRB 031203 host galaxy at z=0.105 shows that the mid-IR lines are emitted from geometrically thin shock dominated filaments which are not reached by the photoionizing flux, while the optical lines are emitted from the radiation dominated outflowing clouds.


title: Surface magnetism of cool stars

abstract: Magnetic fields are essential ingredients of many physical processes in the interiors and envelopes of cool stars. Yet their direct detection and characterisation is notoriously difficult, requiring high-quality observations and advanced analysis techniques. Significant progress has been recently achieved by several types of direct magnetic field studies on the surfaces of cool active stars. In particular, complementary techniques of the field topology mapping with polarisation data and total magnetic flux measurements from intensity spectra have been systematically applied to different classes of active stars leading to interesting and occasionally controversial results. In this paper we summarise the current status of direct magnetic field studies of cool stars, and investigations of surface inhomogeneities caused by the field, based on the material presented at the Cool Stars 19 splinter session.


title: Central limit theorems for supercritical superprocesses with immigration

abstract: In this paper, we establish a central limit theorem for a large class of general supercritical superprocesses with immigration with spatially dependent branching mechanisms satisfying a second moment condition. This central limit theorem extends and generalizes the results obtained by Ren, Song and Zhang[2015]. We first give law of large numbers for supercritical superprocesses with immigration since there is few convergence result on immigration superprocesses, then based on these results, we establish the central limit theorem.


title: On groups of H\"older diffeomorphisms and their regularity

abstract: We study the set $\mathcal D^{n,\beta}(\mathbb R^d)$ of orientation preserving diffeomorphisms of $\mathbb R^d$ which differ from the identity by a H\"older $C^{n,\beta}_0$-mapping, where $n \in \mathbb N_{\ge 1}$ and $\beta \in (0,1]$. We show that $\mathcal D^{n,\beta}(\mathbb R^d)$ forms a group, but left translations in $\mathcal D^{n,\beta}(\mathbb R^d)$ are in general discontinuous. The groups $\mathcal D^{n,\beta-}(\mathbb R^d) := \bigcap_{\alpha < \beta} \mathcal D^{n,\alpha}(\mathbb R^d)$ (with its natural Fr\'echet topology) and $\mathcal D^{n,\beta+}(\mathbb R^d) := \bigcup_{\alpha > \beta} \mathcal D^{n,\alpha}(\mathbb R^d)$ (with its natural inductive locally convex topology) however are $C^{0,\omega}$ Lie groups for any slowly vanishing modulus of continuity $\omega$. In particular, $\mathcal D^{n,\beta-}(\mathbb R^d)$ is a topological group and a so-called half-Lie group (with smooth right translations). We prove that the H\"older spaces $C^{n,\beta}_0$ are ODE closed, in the sense that pointwise time-dependent $C^{n,\beta}_0$-vector fields $u$ have unique flows $\Phi$ in $\mathcal D^{n,\beta}(\mathbb R^d)$. This includes, in particular, all Bochner integrable functions $u \in L^1([0,1],C^{n,\beta}_0(\mathbb R^d,\mathbb R^d))$. For the latter and $n\ge 2$, we show that the flow map $L^1([0,1],C^{n,\beta}_0(\mathbb R^d,\mathbb R^d)) \to C([0,1],\mathcal D^{n,\alpha}(\mathbb R^d))$, $u \mapsto \Phi$, is continuous (even $C^{0,\beta-\alpha}$), for every $\alpha < \beta$. As an application we prove that the corresponding Trouv\'e group $\mathcal G_{n,\beta}(\mathbb R^d)$ from image analysis coincides with the connected component of the identity of $\mathcal D^{n,\beta}(\mathbb R^d)$.


title: Massive Neutral Gas Outflow in Reddened Quasar SDSS J072910.34+333634.3

abstract: SDSS J072910.34+333634.3 is a reddened quasar at z=0.96. The archivel Keck/ESI spectrum and our new P200/TripleSpec spectrum reveal an absorption line system in He I*, Ca II H\&K and Na ID. The absorption line system has a width of $\sim$600 km/s and a blueshift velocity of $\sim$800 km/s relative to the core of narrow emission lines, indicating an outflow. Using the Ca II doublet, we determined that the outflowing gas covers $\sim$70\% of the continuum. On the other hand, the HST/ACS image which taken in rest-frame 4130 \AA\ show that the fraction of the quasar emission in ESI aperture was $<$40\%. We thus conclude that the absorbing gas covers a significant fraction of extended starlight emission, and the best-estimated fraction of $\sim$50\% yields a lower limit of the crosssectional area of the outflowing gas to be $>$8 kpc$^2$. The strong Na\&Ca absorption suggests that the absorbing gas is thick and mostly neutral, which is also supported by dust extinction $A_V\sim3$. Using the best-estimated hydrogen column density $N_H\sim3\times10^{22}$ cm$^{-2}$, the total mass of the outflowing gas is $>3\times10^9 M_\odot$. The outflow is likely to be driven by AGN because of the $\sim$800 km/s blueshift velocity, suggesting SDSS J072910.34+333634.3 is undergoing one of the most violent AGN feedback we have seen. In the future, one can find more such massive neutral gas outflows in other reddened quasars using similar method, and this can shed new light on the study of AGN feedback.


title: Analog curved spacetimes in the reversed dissipation regime of cavity optomechanics

abstract: In this paper, we theoretically propose an optomechanical scheme to explore the possibility of simulating the propagation of the collective excitations of the photon fluid in a curved spacetime. For this purpose, we introduce two theoretical models for two-dimensional photon gas in a planar optomechanical microcavity and a two-dimensional array of coupled optomechanical systems. In the reversed dissipation regime (RDR) of cavity optomechanics where the mechanical oscillator reaches equilibrium with its thermal reservoir much faster than the cavity modes, the mechanical degrees of freedom can adiabatically be eliminated. The adiabatic elimination of the mechanical mode provides an effective nonlinear Kerr-type photon-photon interaction. Using the nonlinear Schr\"{o}dinger equation (NLSE), we show that the phase fluctuations in the two-dimensional photon fluid obey the Klein-Gordon equation for a massless scalar field propagating in a curved spacetime. The results reveal that the photon fluid as well as the corresponding metric can be controlled by manipulating the system parameters.


title: Exact recovery low-rank matrix via transformed affine matrix rank minimization

abstract: The goal of affine matrix rank minimization problem is to reconstruct a low-rank or approximately low-rank matrix under linear constraints. In general, this problem is combinatorial and NP-hard. In this paper, a nonconvex fraction function is studied to approximate the rank of a matrix and translate this NP-hard problem into a transformed affine matrix rank minimization problem. The equivalence between these two problems is established, and we proved that the uniqueness of the global minimizer of transformed affine matrix rank minimization problem also solves affine matrix rank minimization problem if some conditions are satisfied. Moreover, we also proved that the optimal solution to the transformed affine matrix rank minimization problem can be approximately obtained by solving its regularization problem for some proper smaller $\lambda>0$. Lastly, the DC algorithm is utilized to solve the regularization transformed affine matrix rank minimization problem and the numerical experiments on image inpainting problems show that our method performs effectively in recovering low-rank images compared with some state-of-art algorithms.


title: Quantum field theory and classical optics: determining the fine structure constant

abstract: The properties of the vacuum are described by quantum physics including the response to external fields such as electromagnetic radiation. Of the two parameters that govern the details of the electromagnetic field dynamics in vacuum, one is fixed by the requirement of Lorentz invariance $c= 1/\sqrt{\varepsilon_{0} \mu_{0}}$. The other one, $Z_{0}= \sqrt{\mu_{0}/\varepsilon_{0}} = 1/(c\varepsilon_{0})$ and its relation to the quantum vacuum, is discussed in this contribution. Deriving $\varepsilon_{0}$ from the properties of the quantum vacuum implies the derivation of the fine structure constant.