title: PassBio: Privacy-Preserving User-Centric Biometric Authentication

abstract: The proliferation of online biometric authentication has necessitated security requirements of biometric templates. The existing secure biometric authentication schemes feature a server-centric model, where a service provider maintains a biometric database and is fully responsible for the security of the templates. The end-users have to fully trust the server in storing, processing and managing their private templates. As a result, the end-users' templates could be compromised by outside attackers or even the service provider itself. In this paper, we propose a user-centric biometric authentication scheme (PassBio) that enables end-users to encrypt their own templates with our proposed light-weighted encryption scheme. During authentication, all the templates remain encrypted such that the server will never see them directly. However, the server is able to determine whether the distance of two encrypted templates is within a pre-defined threshold. Our security analysis shows that no critical information of the templates can be revealed under both passive and active attacks. PassBio follows a "compute-then-compare" computational model over encrypted data. More specifically, our proposed Threshold Predicate Encryption (TPE) scheme can encrypt two vectors x and y in such a manner that the inner product of x and y can be evaluated and compared to a pre-defined threshold. TPE guarantees that only the comparison result is revealed and no key information about x and y can be learned. Furthermore, we show that TPE can be utilized as a flexible building block to evaluate different distance metrics such as Hamming distance and Euclidean distance over encrypted data. Such a compute-then-compare computational model, enabled by TPE, can be widely applied in many interesting applications such as searching over encrypted data while ensuring data security and privacy.


title: Robust Multilingual Part-of-Speech Tagging via Adversarial Training

abstract: Adversarial training (AT) is a powerful regularization method for neural networks, aiming to achieve robustness to input perturbations. Yet, the specific effects of the robustness obtained from AT are still unclear in the context of natural language processing. In this paper, we propose and analyze a neural POS tagging model that exploits AT. In our experiments on the Penn Treebank WSJ corpus and the Universal Dependencies (UD) dataset (27 languages), we find that AT not only improves the overall tagging accuracy, but also 1) prevents over-fitting well in low resource languages and 2) boosts tagging accuracy for rare / unseen words. We also demonstrate that 3) the improved tagging performance by AT contributes to the downstream task of dependency parsing, and that 4) AT helps the model to learn cleaner word representations. 5) The proposed AT model is generally effective in different sequence labeling tasks. These positive results motivate further use of AT for natural language tasks.


title: Strongly graded groupoids and strongly graded Steinberg algebras

abstract: We study strongly graded groupoids, which are topological groupoids $\mathcal G$ equipped with a continuous, surjective functor $\kappa: \mathcal G \to \Gamma$, to a discrete group $\Gamma$, such that $\kappa^{-1}(\gamma)\kappa^{-1}(\delta) = \kappa^{-1}(\gamma \delta)$, for all $\gamma, \delta \in \Gamma$. We introduce the category of graded $\mathcal G$-sheaves, and prove an analogue of Dade's Theorem: $\mathcal G$ is strongly graded if and only if every graded $\mathcal G$-sheaf is induced by a $\mathcal G_{\epsilon}$-sheaf. The Steinberg algebra of a graded ample groupoid is graded, and we prove that the algebra is strongly graded if and only if the groupoid is. Applying this result, we obtain a complete graphical characterisation of strongly graded Leavitt path and Kumjian-Pask algebras.


title: On a small-scale EUV wave: the driving mechanism and the associated oscillating filament

abstract: We present observations of a small-scale Extreme-ultraviolet (EUV) wave that was associated with a mini-filament eruption and a GOES B1.9 micro-fare in the quiet Sun region. The initiation of the event was due to the photospheric magnetic emergence and cancellation in the eruption source region, which first caused the ejection of a small plasma ejecta, then the ejecta impacted on a nearby mini-filament and thereby led to the filament's eruption and the associated fare. During the filament eruption, an EUV wave at a speed of 182 { 317 km/s was formed ahead of an expanding coronal loop, which propagated faster than the expanding loop and showed obvious deceleration and refection during the propagation. In addition, the EUV wave further resulted in the transverse oscillation of a remote filament whose period and damping time are 15 and 60 minutes, respectively. Based on the observational results, we propose that the small-scale EUV wave should be a fast-mode magnetosonic wave that was driven by the the expanding coronal loop. Moreover, with the application of filament seismology, it is estimated that the radial magnetic field strength is about 7 Gauss. The observations also suggest that small-scale EUV waves associated with miniature solar eruptions share similar driving mechanism and observational characteristics with their large-scale counterparts.


title: The Maximal Rank Conjecture

abstract: Let C be a general curve of genus g, embedded in P^r via a general linear series of degree d. In this paper, we prove the Maximal Rank Conjecture, which determines the Hilbert function of C.


title: A Family of Constrained Adaptive filtering Algorithms Based on Logarithmic Cost

abstract: This paper introduces a novel constraint adaptive filtering algorithm based on a relative logarithmic cost function which is termed as Constrained Least Mean Logarithmic Square (CLMLS). The proposed CLMLS algorithm elegantly adjusts the cost function based on the amount of error thereby achieves better performance compared to the conventional Constrained LMS (CLMS) algorithm. With no assumption on input, the mean square stability analysis of the proposed CLMLS algorithm is presented using the energy conservation approach. The analytical expressions for the transient and steady state MSD are derived and these analytical results are validated through extensive simulations. The proposed CLMLS algorithm is extended to the sparse case by incorporating the $\ell_1$-norm penalty into the CLMLS cost function. detailed Simulations confirms the superiority of the sparse CLMLS over the state-of-the-art.


title: Theory and evidence of global Rossby waves in upper main-sequence stars: r-mode oscillations in many Kepler stars

abstract: Asteroseismic inference from pressure modes (p modes) and buoyancy, or gravity, modes (g modes) is ubiquitous for stars across the Hertzsprung--Russell diagram. Until now, however, discussion of r modes (global Rossby waves) has been rare. Here we derive the expected frequency ranges of r modes in the observational frame by considering the visibility of these modes. We find that the frequencies of r modes of azimuthal order $m$ appear as groups at slightly lower frequency than $m$ times the rotation frequency. Comparing the visibility curves for r modes with Fourier amplitude spectra of Kepler light curves of upper main-sequence B, A and F stars, we find that r modes are present in many $\gamma$ Dor stars (as first discovered by Van Reeth et al. 2016), spotted stars, and so-called Heartbeat stars, which are highly eccentric binary stars. We also find a signature of r modes in a frequently bursting Be star observed by Kepler. In the amplitude spectra of moderately to rapidly rotating $\gamma$ Dor stars, r-mode frequency groups appear at lower frequency than prograde g-mode frequency groups, while in the amplitude spectra of spotted early A to B stars, groups of symmetric (with respect to the equator) r-mode frequencies appear just below the frequency of a structured peak that we suggest represents an approximate stellar rotation rate. In many Heartbeat stars, a group of frequencies can be fitted with symmetric $m=1$ r modes, which can be used to obtain rotation frequencies of these stars.


title: An Optimal Convergence Rate for the Gaussian Regularized Shannon Sampling Series

abstract: We consider the reconstruction of a bandlimited function from its finite localized sample data. Truncating the classical Shannon sampling series results in an unsatisfactory convergence rate due to the slow decay of the sinc function. To overcome this drawback, a simple and highly effective method, called the Gaussian regularization of the Shannon series, was proposed in the engineering and has received remarkable attention. It works by multiplying the sinc function in the Shannon series with a regularized Gaussian function. Recently, it was proved that the upper error bound of this method can achieve a convergence rate of the order $O(\frac{1}{\sqrt{n}}\exp(-\frac{\pi-\delta}{2}n))$, where $0<\delta<\pi$ is the bandwidth and $n$ the number of sample data. The convergence rate is by far the best convergence rate among all regularized methods for the Shannon sampling series. The main objective of this article is to present the theoretical justification and numerical verification that the convergence rate is optimal when $0<\delta<\pi/2$ by estimating the lower error bound of the truncated Gaussian regularized Shannon sampling series.


title: Skeleton-stabilized IsoGeometric Analysis: High-regularity Interior-Penalty methods for incompressible viscous flow problems

abstract: A Skeleton-stabilized IsoGeometric Analysis (SIGA) technique is proposed for incompressible viscous flow problems with moderate Reynolds number. The proposed method allows utilizing identical finite dimensional spaces (with arbitrary B-splines/NURBS order and regularity) for the approximation of the pressure and velocity components. The key idea is to stabilize the jumps of high-order derivatives of variables over the skeleton of the mesh. For B-splines/NURBS basis functions of degree $k$ with $C^{\alpha}$-regularity ($0 \leq \alpha < k$), only the derivative of order $\alpha +1$ has to be controlled. This stabilization technique thus can be viewed as a high-regularity generalization of the (Continuous) Interior-Penalty Finite Element Method. Numerical experiments are performed for the Stokes and Navier-Stokes equations in two and three dimensions. Oscillation-free solutions and optimal convergence rates are obtained. In terms of the sparsity pattern of the algebraic system, we demonstrate that the block matrix associated with the stabilization term has a considerably smaller bandwidth when using B-splines than when using Lagrange basis functions, even in the case of $C^0$-continuity. This important property makes the proposed isogeometric framework practical from a computational effort point of view.


title: Magnification by Gravity

abstract: Gravitational lensing time-delays using standard candles can provide tremendous leverage to cosmology. A recently discovered gravitationally lensed type Ia supernova, iPTF16geu, may yield a precise and accurate measurement of the Hubble constant.


title: Revisiting Higgs inflation in the context of collapse theories

abstract: In this work we consider the Higgs inflation scenario, but in contrast with past works, the present analysis is done in the context of a spontaneous collapse theory for the quantum state of the inflaton field. In particular, we will rely on a previously studied adaptation of the Continuous Spontaneous Localization model for the treatment of inflationary cosmology. We will show that with the introduction of the dynamical collapse hypothesis, some of the most serious problems of the Higgs inflation proposal can be resolved in a natural way.


title: pyLEMMINGS: Large Margin Multiple Instance Classification and Ranking for Bioinformatics Applications

abstract: Motivation: A major challenge in the development of machine learning based methods in computational biology is that data may not be accurately labeled due to the time and resources required for experimentally annotating properties of proteins and DNA sequences. Standard supervised learning algorithms assume accurate instance-level labeling of training data. Multiple instance learning is a paradigm for handling such labeling ambiguities. However, the widely used large-margin classification methods for multiple instance learning are heuristic in nature with high computational requirements. In this paper, we present stochastic sub-gradient optimization large margin algorithms for multiple instance classification and ranking, and provide them in a software suite called pyLEMMINGS. Results: We have tested pyLEMMINGS on a number of bioinformatics problems as well as benchmark datasets. pyLEMMINGS has successfully been able to identify functionally important segments of proteins: binding sites in Calmodulin binding proteins, prion forming regions, and amyloid cores. pyLEMMINGS achieves state-of-the-art performance in all these tasks, demonstrating the value of multiple instance learning. Furthermore, our method has shown more than 100-fold improvement in terms of running time as compared to heuristic solutions with improved accuracy over benchmark datasets. Availability and Implementation: pyLEMMINGS python package is available for download at: http://faculty.pieas.edu.pk/fayyaz/software.html#pylemmings.


title: Finite volume effects on the chiral phase transition from Dyson-Schwinger equations of QCD

abstract: Within the framework of Dyson-Schwinger equations of QCD, we study the finite volume effects on the chiral phase transition, especially the influence on the position of the possible pseudo-critical end point (pCEP). The results show that in the chiral limit case (the current quark mass $m=0$), the absolute value of quark condensate decreases for smaller volumes, and more interestingly, so does the pseudo-critical temperature $T_c(\mu=0)$, which is in agreement with the Polyakov Nambu--Jona-Lasinio model result and opposite to the Polyakov linear sigma model prediction. These conclusions hold for $m>0$ case in our calculations. Moreover, the results of pCEP as a function of different volumes show that $T$ of pCEP also decreases for smaller volumes, but $\mu$ of pCEP will increase, which are qualitatively more close to Polyakov linear sigma model results. For our model setup, results for systems with a size larger than (5 fm)$^3$ closely approximate those from infinite volume, but if the volume is smaller, the corrections are non-negligible, even significantly affect signatures of the results from an infinite system. There also exists some possibility that, if the system size is too small, the whole phase transition would be crossover, which means no pCEP exists at all. It is no doubt that, finite volume effects deserve further researches.


title: Adversarial Symmetric Variational Autoencoder

abstract: A new form of variational autoencoder (VAE) is developed, in which the joint distribution of data and codes is considered in two (symmetric) forms: ($i$) from observed data fed through the encoder to yield codes, and ($ii$) from latent codes drawn from a simple prior and propagated through the decoder to manifest data. Lower bounds are learned for marginal log-likelihood fits observed data and latent codes. When learning with the variational bound, one seeks to minimize the symmetric Kullback-Leibler divergence of joint density functions from ($i$) and ($ii$), while simultaneously seeking to maximize the two marginal log-likelihoods. To facilitate learning, a new form of adversarial training is developed. An extensive set of experiments is performed, in which we demonstrate state-of-the-art data reconstruction and generation on several image benchmark datasets.


title: Generative Steganography with Kerckhoffs' Principle

abstract: The distortion in steganography that usually comes from the modification or recoding on the cover image during the embedding process leaves the steganalyzer with possibility of discriminating. Faced with such a risk, we propose generative steganography with Kerckhoffs' principle (GSK) in this letter. In GSK, the secret messages are generated by a cover image using a generator rather than embedded into the cover, thus resulting in no modifications in the cover. To ensure the security, the generators are trained to meet Kerckhoffs' principle based on generative adversarial networks (GAN). Everything about the GSK system, except the extraction key, is public knowledge for the receivers. The secret messages can be outputted by the generator if and only if the extraction key and the cover image are both inputted. In the generator training procedures, there are two GANs, Message- GAN and Cover-GAN, designed to work jointly making the generated results under the control of the extraction key and the cover image. We provide experimental results on the training process and give an example of the working process by adopting a generator trained on MNIST, which demonstrate that GSK can use a cover image without any modification to generate messages, and without the extraction key or the cover image, only meaningless results would be obtained.


title: Rydberg-atom-based scheme of nonadiabatic geometric quantum computation

abstract: Nonadiabatic geometric quantum computation provides a means to perform fast and robust quantum gates. It has been implemented in various physical systems, such as trapped ions, nuclear magnetic resonance and superconducting circuits. Another system being adequate for implementation of nonadiabatic geometric quantum computation may be Rydberg atoms, since their internal states have very long coherence time and the Rydberg-mediated interaction facilitates the implementation of a two-qubit gate. Here, we propose a scheme of nonadiabatic geometric quantum computation based on Rydberg atoms, which combines the robustness of nonadiabatic geometric gates with the merits of Rydberg atoms.


title: Hardy Space Decompositions of $L^p(\mathbb{R}^n)$ for $0<p<1$ with Rational Approximation

abstract: This paper aims to obtain decompositions of higher dimensional $L^p(\mathbb{R}^n)$ functions into sums of non-tangential boundary limits of the corresponding Hardy space functions on tubes for the index range $0<p<1$. In the one-dimensional case, Deng and Qian \cite{DQ} recently obtained such Hardy space decomposition result: for any function $f\in L^p(\mathbb{R}),\ 0<p<1$, there exist functions $f_1$ and $f_2$ such that $f=f_1+f_2$, where $f_1$ and $f_2$ are, respectively, the non-tangential boundary limits of some Hardy space functions in the upper-half and lower-half planes. In the present paper, we generalize the one-dimensional Hardy space decomposition result to the higher dimensions, and discuss the uniqueness issue of such decomposition.


title: Strong Exciton-Photon Coupling and lasing behavior in All-Inorganic CsPbBr3 Micro/nanowire Fabry-Perot cavity

abstract: Macroscopic spontaneous coherence of exciton-polariton in semiconductor cavity is one important research field in condensed matter physics. All-inorganic micro/nanowire Fabry-P\'erot cavity with excellent optical performance makes it possible to realize strong coupling and lasing of exciton-photon at room temperature. In this work, we demonstrated strong coupling of exciton-photon and polariton lasing in CsPbBr3 micro/nanowires synthesized by CVD method. By exploring spatial resolved PL spectra of CsPbBr3 cavity, we observed mode volume dependent coupling strength , as well as significant increase in group index. Moreover, low threshold polariton lasing was achieved at room temperature within strong coupling regime; the polariton characteristic is confirmed by comparing lasing spectra with waveguided output spectra and the dramatically reduced lasing threshold. Our present results provide new avenues to achieve high coupling strengths potentially enabling application of exciting phenomena such as Bose-Einstein condensation of polaritons, efficient light-emitting diodes and lasers.


title: Size-scaling limits of impulsive elastic energy release from a resilin-like elastomer

abstract: Elastically-driven motion has been used as a strategy to achieve high speeds in small organisms and engineered micro-robotic devices. We examine the size-scaling relations determining the limit of elastic energy release from elastomer bands with mechanical properties similar to the biological protein resilin. The maximum center-of-mass velocity of the elastomer bands was found to be size-scale independent, while smaller bands demonstrated larger accelerations and shorter durations of elastic energy release. Scaling relationships determined from these measurements are consistent with the performance of small organisms which utilize elastic elements to power motion. Engineered devices found in the literature do not follow the same size-scaling relationships, which suggests an opportunity for improved design of engineered devices.


title: Lattice Surgery on the Raussendorf Lattice

abstract: Lattice surgery is a method to perform quantum computation fault-tolerantly by using operations on boundary qubits between different patches of the planar code. This technique allows for universal planar-code computation without eliminating the intrinsic two-dimensional nearest-neighbor properties of the surface code that eases physical hardware implementations. Lattice-surgery approaches to algorithmic compilation and optimization have been demonstrated to be more resource efficient for resource-intensive components of a fault-tolerant algorithm, and consequently may be preferable over braid-based logic. Lattice surgery can be extended to the Raussendorf lattice, providing a measurement-based approach to the surface code. In this paper we describe how lattice surgery can be performed on the Raussendorf lattice and therefore give a viable alternative to computation using braiding in measurement based implementations of topological codes.


title: Localized massive halo properties in Bahamas and Macsis simulations: scalings, log-normality, and covariance

abstract: Using tens of thousands of halos realized in the BAHAMAS and MACSIS simulations produced with a consistent astrophysics treatment that includes AGN feedback, we validate a multi-property statistical model for the stellar and hot gas mass behavior in halos hosting groups and clusters of galaxies. The large sample size allows us to extract fine-scale mass--property relations (MPRs) by performing local linear regression (LLR) on individual halo stellar mass (${\rm M}_{\rm star}$) and hot gas mass (${\rm M}_{\rm gas}$) as a function of total halo mass (${\rm M}_{\rm halo}$). We find that: 1) both the local slope and variance of the MPRs run with mass (primarily) and redshift (secondarily); 2) the conditional likelihood, $p({\rm M}_{\rm star},\ {\rm M}_{\rm gas} | \ {\rm M}_{\rm halo}, z)$ is accurately described by a multivariate, log-normal distribution, and; 3) the covariance of ${\rm M}_{\rm star}$ and ${\rm M}_{\rm gas}$ at fixed ${\rm M}_{\rm halo}$ is generally negative, reflecting a partially closed baryon box model for high mass halos. We validate the analytical population model of Evrard et al. (2014), finding sub-percent accuracy in the log-mean halo mass selected at fixed property, $\langle \ln {\rm M}_{\rm halo} | {\rm M}_{\rm gas} \rangle$ or $\langle \ln {\rm M}_{\rm halo} | {\rm M}_{\rm star} \rangle$, when scale-dependent MPR parameters are employed. This work highlights the potential importance of allowing for running in the slope and scatter of MPRs when modeling cluster counts for cosmological studies. We tabulate LLR fit parameters as a function of halo mass at $z=0$, $0.5$ and 1 for two popular mass conventions.


title: Higgs mode in the d-wave superconductor Bi2Sr2CaCu2O8+x driven by an intense terahertz pulse

abstract: We investigated the terahertz (THz)-pulse driven nonlinear response in the d-wave cuprate superconductor Bi2Sr2CaCu2O8+x (Bi2212) using a THz pump near-infrared probe scheme in the time domain. We have observed an oscillatory behavior of the optical reflectivity that follows the THz electric field squared and is strongly enhanced below Tc. The corresponding third-order nonlinear effect exhibits both A1g and B1g symmetry components, which are decomposed from polarization-resolved measurements. Comparison with a BCS calculation of the nonlinear susceptibility indicates that the A1g component is associated with the Higgs mode of the d-wave order parameter.


title: On Fermat Diophantine functional equations and little Picard theorem

abstract: We discuss equivalence conditions on the non-existence of non-trivial meromorphic solution to the Fermat Diophantine equations $f^m(z)+g^n(z)=1$ with integers $m,n\geq2$, from which other approaches to prove little Picard theorem are provided.


title: Global Weak Solutions for the Compressible Active Liquid Crystal System

abstract: We study the hydrodynamics of compressible flows of active liquid crystals in the Beris-Edwards hydrodynamics framework, using the Landau-de Gennes $Q$-tensor order parameter to describe liquid crystalline ordering. We prove the existence of global weak solutions for this active system in three space dimensions by the three-level approximations and weak convergence argument. New techniques and estimates are developed to overcome the difficulties caused by the active terms.


title: The dominance of the $\nu(0d_{5/2})^2$ configuration in the $N=8$ shell in $^{12}$Be from the breakup reaction on a proton target at intermediate energy

abstract: The momentum distribution of $^{11}$Be fragments produced by the breakup of $^{12}$Be interacting with a proton target at 700.5 MeV/$u$ energy has been measured at GSI Darmstadt. To obtain the structure information on the anomaly of the $N=8$ neutron shell, the momentum distribution of $^{11}$Be fragments from the one-neutron knockout $^{12}$Be(p,pn) reaction, measured in inverse kinematics, has been analysed in the distorted wave impulse approximation (DWIA) based on a quasi-free scattering scenario. The DWIA analysis shows a surprisingly strong contribution of the neutron $0d_{5/2}$ orbital in $^{12}$Be to the transverse momentum distribution of the $^{11}$Be fragments. The single-neutron $0d_{5/2}$ spectroscopic factor deduced from the present knock-out data is 1.39(10), which is significantly larger than that deduced recently from data of $^{12}$Be breakup on a carbon target. This result provides a strong experimental evidence for the dominance of the neutron $\nu(0d_{5/2})^2$ configuration in the ground state of $^{12}$Be.


title: APO Time Resolved Color Photometry of Highly-Elongated Interstellar Object 1I/'Oumuamua

abstract: We report on $g$, $r$ and $i$ band observations of the Interstellar Object 'Oumuamua (1I) taken on 2017 October 29 from 04:28 to 08:40 UTC by the Apache Point Observatory (APO) 3.5m telescope's ARCTIC camera. We find that 1I's colors are $g-r=0.41\pm0.24$ and $r-i=0.23\pm0.25$, consistent with the visible spectra of Masiero (2017), Ye et al. (2017) and Fitzsimmons et al. (2017), and most comparable to the population of Solar System C/D asteroids, Trojans, or comets. We find no evidence of any cometary activity at a heliocentric distance of 1.46 au, approximately 1.5 months after 1I's closest approach distance to the Sun. Significant brightness variability was seen in the $r$ observations, with the object becoming notably brighter towards the end of the run. By combining our APO photometric time series data with the Discovery Channel Telescope (DCT) data of Knight et al. (2017), taken 20 h later on 2017 October 30, we construct an almost complete light curve with a most probable lightcurve period of $P \simeq 4~{\rm h}$. Our results imply a double peaked rotation period of 8.1 $\pm$ 0.02 h, with a peak-to-peak amplitude of 1.5 - 2.1 mags. Assuming that 1I's shape can be approximated by an ellipsoid, the amplitude constraint implies that 1I has an axial ratio of 3.5 to 10.3, which is strikingly elongated. Assuming that 1I is rotating above its critical break up limit, our results are compatible with 1I having having modest cohesive strength and may have obtained its elongated shape during a tidal disruption event before being ejected from its home system. Astrometry useful for constraining 1I's orbit was also obtained and published in Weaver et al. (2017).


title: Tracer Transport within an Unstructured Grid Ocean Model using Characteristic Discontinuous Galerkin Advection

abstract: In a previous article a characteristic discontinuous Galerkin (CDG) advection scheme was presented for tracer transport. The scheme is conservative, unconditionally stable with respect to time step and scales sub-linearly with the number of tracers being advected. Here we present the implementation of the CDG advection scheme for tracer transport within MPAS-Ocean, a Boussinesque unstructured grid ocean model with an arbitrary Lagrangian Eulerian vertical coordinate. The scheme is implemented in both the vertical and horizontal dimensions, and special care is taken to ensure that the scheme remains conservative in the context of moving vertical layers. Consistency is ensured with respect to the dynamics by a renormalization of the fluxes with respect to the volume fluxes derived from the continuity equation. For spherical implementations, the intersection of the flux swept regions and the Eulerian grid are determined for great circle arcs, and the fluxes and element assembly are performed on the plane via a length preserving projection. Solutions are presented for a suite of test cases and comparisons made to the existing flux corrected transport scheme in MPAS-Ocean.


title: Anomalous vibrational properties in the continuum limit of glasses

abstract: The low-temperature thermal properties of glasses are anomalous with respect to those of crystals. These thermal anomalies indicate that the low-frequency vibrational properties of glasses differ from those of crystals. Recent studies revealed that, in the simplest model of glasses, i.e., the harmonic potential system, phonon modes coexist with soft localized modes in the low-frequency (continuum) limit. However, the nature of low-frequency vibrational modes of more realistic models is still controversial. In the present work, we study the Lennard-Jones (LJ) system using large-scale molecular-dynamics (MD) simulation and establish that the vibrational property of the LJ glass converges to coexistence of the phonon modes and the soft localized modes in the continuum limit as in the case of the harmonic potential system. Importantly, we find that the low-frequency vibrations are rather sensitive to the numerical scheme of potential truncation which is usually implemented in the MD simulation, and this is the reason why contradictory arguments have been reported by previous works. We also discuss the physical origin of this sensitiveness by means of a linear stability analysis.


title: Type III Solar Radio Burst Source Region Splitting Due to a Quasi-Separatrix Layer

abstract: We present low-frequency (80-240 MHz) radio imaging of type III solar radio bursts observed by the Murchison Widefield Array (MWA) on 2015/09/21. The source region for each burst splits from one dominant component at higher frequencies into two increasingly-separated components at lower frequencies. For channels below ~132 MHz, the two components repetitively diverge at high speeds (0.1-0.4 c) along directions tangent to the limb, with each episode lasting just ~2 s. We argue that both effects result from the strong magnetic field connectivity gradient that the burst-driving electron beams move into. Persistence mapping of extreme ultraviolet (EUV) jets observed by the Solar Dynamics Observatory reveals quasi-separatrix layers (QSLs) associated with coronal null points, including separatrix dome, spine, and curtain structures. Electrons are accelerated at the flare site toward an open QSL, where the beams follow diverging field lines to produce the source splitting, with larger separations at larger heights (lower frequencies). The splitting motion within individual frequency bands is interpreted as a projected time-of-flight effect, whereby electrons traveling along the outer field lines take slightly longer to excite emission at adjacent positions. Given this interpretation, we estimate an average beam speed of 0.2 c. We also qualitatively describe the quiescent corona, noting in particular that a disk-center coronal hole transitions from being dark at higher frequencies to bright at lower frequencies, turning over around 120 MHz. These observations are compared to synthetic images based on the Magnetohydrodynamic Algorithm outside a Sphere (MAS) model, which we use to flux-calibrate the burst data.


title: Magnon Splitting Induced by Charge Transfer in the Three-Orbital Hubbard Model

abstract: Understanding spin excitations and their connection to unconventional superconductivity have remained a central issue since the discovery of the cuprates. Direct measurement of the dynamical spin structure factor in the parent compounds can provide key information on important interactions relevant in the doped regime, and variations in the magnon dispersion have been linked closely to differences in crystal structure between families of cuprate compounds. Here, we elucidate the relationship between spin excitations and various controlling factors thought to be significant in high-$T_c$ materials by systematically evaluating the dynamical spin structure factor for the three-orbital Hubbard model, revealing differences in the spin dispersion along the Brillouin zone axis and the diagonal. Generally, we find that the absolute energy scale and momentum dependence of the excitations primarily are sensitive to the effective charge-transfer energy, while changes in the on-site Coulomb interactions have little effect on the details of the dispersion. In particular, our result highlights the splitting between spin excitations along the axial and diagonal directions in the Brillouin zone. This splitting decreases with increasing charge-transfer energy and correlates with changes in the apical oxygen position, and general structural variations, for different cuprate families.


title: Poisson Statistics in the Non-Homogeneous Hierarchical Anderson Model

abstract: In this article we study the problem of localization of eigenvalues for the non-homogeneous hierarchical Anderson model. More specifically, given the hierarchical Anderson model with spectral dimension $0<d<1$ with a random potential acting on the diagonal of non i.i.d. random variables, sufficient conditions on the disorder are provided in order to obtain the two main results: the weak convergence of the counting measure for almost all realization of the random potential and the weak convergence of the re-scaled eigenvalue counting measure to a Poisson point process. The technical part improves the already existing arguments of Kritchevski , who studied the hierarchical model with a disorder acting on the diagonal, with independent and identically distributed random variables, by using the argument of Minami . At the end of this article, we study an application example that allows us to understand some relations between the spectral dimension of the hierarchical Laplacian and the magnitude of the disorder.


title: A phase space description of the FLRW quantum cosmology in Ho$\check{\rm r}$ava-Lifshitz type gravity

abstract: Quantum cosmology of the Friedmann-Lema\^{i}tre-Robertson-Walker model with cosmological constant in the Ho$\check{\rm r}$ava-Lifshitz type gravity is studied in the phase space by means of the Wigner function. The modification of the usual general relativity description by the Ho$\check{\rm r}$ava-Lifshitz type gravity induces a new scenario for the origin of the Universe with an embryonic era where the Universe can exist classically before the tunneling process takes place and which gives rise to the current evolution of the Universe. The Wigner functions corresponding to the Hartle-Hawking, Vilenkin and Linde boundary conditions are obtained by means of numerical calculations. In particular three cases were studied for the potential of the Wheeler-DeWitt equation: tunneling barrier with and without embryonic era and when the potential barrier is not present. The quantum behavior of these three cases are analyzed using the Wigner function for the three boundary conditions considered.


title: Statistically Optimal and Computationally Efficient Low Rank Tensor Completion from Noisy Entries

abstract: In this article, we develop methods for estimating a low rank tensor from noisy observations on a subset of its entries to achieve both statistical and computational efficiencies. There have been a lot of recent interests in this problem of noisy tensor completion. Much of the attention has been focused on the fundamental computational challenges often associated with problems involving higher order tensors, yet very little is known about their statistical performance. To fill in this void, in this article, we characterize the fundamental statistical limits of noisy tensor completion by establishing minimax optimal rates of convergence for estimating a $k$th order low rank tensor under the general $\ell_p$ ($1\le p\le 2$) norm which suggest significant room for improvement over the existing approaches. Furthermore, we propose a polynomial-time computable estimating procedure based upon power iteration and a second-order spectral initialization that achieves the optimal rates of convergence. Our method is fairly easy to implement and numerical experiments are presented to further demonstrate the practical merits of our estimator.

categories: stat.ML cs.IT math.IT math.ST stat.ME stat.TH


title: Tree-Based Unrooted Nonbinary Phylogenetic Networks

abstract: Phylogenetic networks are a generalisation of phylogenetic trees that allow for more complex evolutionary histories that include hybridisation-like processes. It is of considerable interest whether a network can be considered `tree-like' or not, which lead to the introduction of \textit{tree-based} networks in the rooted, binary context. Tree-based networks are those networks which can be constructed by adding additional edges into a given phylogenetic tree, called the \textit{base tree}. Previous extensions have considered extending to the binary, unrooted case and the nonbinary, rooted case. We extend tree-based networks to the context of unrooted, nonbinary networks in three ways, depending on the types of additional edges that are permitted. A phylogenetic network in which every embedded tree is a base tree is termed a \textit{fully tree-based} network. We also extend this concept to unrooted, nonbinary phylogenetic networks and classify the resulting networks. We also derive some results on the colourability of tree-based networks, which can be useful to determine whether a network is tree-based.


title: A Two-Stage Approach for Routing Multiple Unmanned Aerial Vehicles with Stochastic Fuel Consumption

abstract: The past decade has seen a substantial increase in the use of small unmanned aerial vehicles (UAVs) in both civil and military applications. This article addresses an important aspect of refueling in the context of routing multiple small UAVs to complete a surveillance or data collection mission. Specifically, this article formulates a multiple-UAV routing problem with the refueling constraint of minimizing the overall fuel consumption for all of the vehicles as a two-stage stochastic optimization problem with uncertainty associated with the fuel consumption of each vehicle. The two-stage model allows for the application of sample average approximation (SAA). Although the SAA solution asymptotically converges to the optimal solution for the two-stage model, the SAA run time can be prohibitive for medium- and large-scale test instances. Hence, we develop a tabu-search-based heuristic that exploits the model structure while considering the uncertainty in fuel consumption. Extensive computational experiments corroborate the benefits of the two-stage model compared to a deterministic model and the effectiveness of the heuristic for obtaining high-quality solutions.


title: Operational effects of the UNOT gate on classical and quantum correlations

abstract: The NOT gate that flips a classical bit is ubiquitous in classical information processing. However its quantum analogue, the universal NOT (UNOT) gate that flips a quantum spin in any alignment into its antipodal counterpart is strictly forbidden. Here we explore the connection between this discrepancy and how UNOT gates affect classical and quantum correlations. We show that while a UNOT gate always preserves classical correlations between two spins, it can non-locally increase or decrease their shared discord in ways that allow violation of the data processing inequality. We experimentally illustrate this using a multi-level trapped \Yb ion that allows simulation of anti-unitary operations.


title: On a theorem of Magidor

abstract: Assuming $\kappa$ is a supercompact cardinal and $\lambda$ is an inaccessible cardinal above it, we present an idea due to Magidor, to find a generic extension in which $\kappa=\aleph_\omega$ and $\lambda=\aleph_{\omega+1}.$


title: Spontaneous lateral atomic recoil force close to a photonic topological material

abstract: We investigate the quantum recoil force acting on an excited atom close to the surface of a nonreciprocal photonic topological insulator (PTI). The main atomic emission channel is the unidirectional surface-plasmon propagating at the PTI-vacuum interface, and we show that it enables a spontaneous lateral recoil force that scales at short distance as $1/d^4$, where $d$ is the atom-PTI separation. Remarkably, the sign of the recoil force is polarization and orientation-independent, and it occurs in a translation-invariant homogeneous system in thermal equilibrium. Surprisingly, the recoil force persists for very small values of the gyration pseudovector, which, for a biased plasma, corresponds to very low cyclotron frequencies. The ultra-strong recoil force is rooted on the quasi-hyperbolic dispersion of the surface-plasmons. We consider both an initially excited atom and a continuous pump scenario, the latter giving rise to a continuous lateral force whose direction can be changed at will by simply varying the orientation of the biasing magnetic field. Our predictions may be tested in experiments with cold Rydberg atoms and superconducting qubits.


title: New relativistic effective interaction for finite nuclei, infinite nuclear matter and neutron stars

abstract: We carry out the study for finite nuclei, infinite nuclear matter and neutron star properties with the newly developed relativistic force named as the Institute Of Physics Bhubaneswar-I(IOPB-I). Using this force, we calculate the binding energies, charge radii and neutron skin thickness for some selected nuclei. From the ground state properties of superheavy element i.e. Z=120, it is noticed that considerable shell gaps appear at neutron numbers N=172, 184 and 198, showing the magicity of these numbers. The low density behavior of the equation of state for pure neutron matter is compatible with other microscopic models. Along with the nuclear symmetry energy, its slope and curvature parameters at the saturation density are consistent with those extracted from various experimental data. We calculate the neutron star properties with the equation of state composed of nucleons and leptons in $\it beta-equilibrium$ which are in good agreement with the X-ray observations by Steiner and N\"{a}ttil\"{a}. We find that the maximum mass of the neutron star to be 2.15$M_{\odot}$ and stellar radius 11.936 km . Moreover, the radius and tidal deformability of a {\it canonical} neutron star mass 1.4$M_\odot$ come out to be 13.242 km and 3.910$\times$10$^{36}$ g cm$^2$ s$^2$ respectively within this parameter set.


title: Fluctuation-induced forces on an atom near a photonic topological material

abstract: We theoretically study the Casimir-Polder force on an atom in a arbitrary initial state in a rather general electromagnetic environment wherein the materials may have a nonreciprocal bianisotropic dispersive response. It is shown that under the Markov approximation the force has resonant and nonresonant contributions. We obtain explicit expressions for the optical force both in terms of the system Green function and of the electromagnetic modes. We apply the theory to the particular case wherein a two-level system interacts with a topological gyrotropic material, showing that the nonreciprocity enables exotic light-matter interactions and the opportunity to sculpt and tune the Casimir-Polder forces on the nanoscale. With a quasi-static approximation, we obtain a simple analytical expression for the optical force and unveil the crucial role of surface plasmons in fluctuation induced forces. Finally, we derive the Green function for a gyrotropic material half-space in terms of a Sommerfeld integral.


title: Corks with large shadow-complexity and exotic 4-manifolds

abstract: We construct an infinite family $\{ C_{n,k}\}_{k=1}^{\infty}$ of corks of Mazur type satisfying $2n\leq \mathrm{sc}^{\mathrm{sp}}(C_{n,k})\leq O(n^{3/2})$ for any positive integer $n$. Furthermore, using these corks, we construct an infinite family $\{(W_{n,k},W'_{n,k})\}_{k=1}^{\infty}$ of exotic pairs of $4$-manifolds with boundary whose special shadow-complexities satisfy the above inequalities. We also discuss exotic pairs with small shadow-complexity.


title: W$^*$-superrigidity for coinduced actions

abstract: We prove W$^*$-superrigidity for a large class of coinduced actions. We prove that if $\Sigma$ is an amenable almost-malnormal subgroup of an infinite conjugagy class (icc) property (T) countable group $\Gamma$, the coinduced action $\Gamma\curvearrowright X$ from an arbitrary probability measure preserving action $\Sigma\curvearrowright X_0$ is W$^*$-superrigid. We also prove a similar statement if $\Gamma$ is an icc non-amenable group which is measure equivalent to a product of two infinite groups. In particular, we obtain that any Bernoulli action of such a group $\Gamma$ is W$^*$-superrigid.


title: Excitation of a global plasma mode by an intense electron beam in a dc discharge

abstract: Interaction of an intense electron beam with a finite-length, inhomogeneous plasma is investigated numerically. The plasma density profile is maximal in the middle and decays towards the plasma edges. Two regimes of the two-stream instability are observed. In one regime, the frequency of the instability is the plasma frequency at the density maximum and plasma waves are excited in the middle of the plasma. In the other regime, the frequency of the instability matches the local plasma frequency near the edges of the plasma and the intense plasma oscillations occur near plasma boundaries. The latter regime appears sporadically and only for strong electron beam currents. This instability generates copious amount of suprathermal electrons. The energy transfer to suprathermal electrons is the saturation mechanism of the instability.


title: Capturing Localized Image Artifacts through a CNN-based Hyper-image Representation

abstract: Training deep CNNs to capture localized image artifacts on a relatively small dataset is a challenging task. With enough images at hand, one can hope that a deep CNN characterizes localized artifacts over the entire data and their effect on the output. However, on smaller datasets, such deep CNNs may overfit and shallow ones find it hard to capture local artifacts. Thus some image-based small-data applications first train their framework on a collection of patches (instead of the entire image) to better learn the representation of localized artifacts. Then the output is obtained by averaging the patch-level results. Such an approach ignores the spatial correlation among patches and how various patch locations affect the output. It also fails in cases where few patches mainly contribute to the image label. To combat these scenarios, we develop the notion of hyper-image representations. Our CNN has two stages. The first stage is trained on patches. The second stage utilizes the last layer representation developed in the first stage to form a hyper-image, which is used to train the second stage. We show that this approach is able to develop a better mapping between the image and its output. We analyze additional properties of our approach and show its effectiveness on one synthetic and two real-world vision tasks - no-reference image quality estimation and image tampering detection - by its performance improvement over existing strong baselines.


title: Floquet states of kicked particle in a singular potential: Exponential and power-law profiles

abstract: It is well known that, in the chaotic regime, all the Floquet states of kicked rotor system display an exponential profile resulting from dynamical localization. If the kicked rotor is placed in an additional stationary infinite potential well, its Floquet states display power-law profile. It has also been suggested in general that the Floquet states of periodically kicked systems with singularities in the potential would have power-law profile. In this work, we study the Floquet states of a kicked particle in finite potential barrier. By varying the height of finite potential barrier, the nature of transition in the Floquet state from exponential to power-law decay profile is studied. We map this system to a tight binding model and show that the nature of decay profile depends on energy band spanned by the Floquet states (in unperturbed basis) relative to the potential height. This property can also be inferred from the statistics of Floquet eigenvalues and eigenvectors. This leads to an unusual scenario in which the level spacing distribution, as a window in to the spectral correlations, is not a unique characteristic for the entire system.


title: A rapid cosmic-ray increase in BC 3372-3371 from ancient buried tree rings in China

abstract: Cosmic rays interact with the Earth's atmosphere to produce $^{14}$C, which can be absorbed by trees. Therefore, rapid increases of $^{14}$C in tree rings can be used to probe previous cosmic-ray events. By this method, three $^{14}$C rapidly increasing events have been found. Plausible causes of these events include large solar proton events, supernovae or short gamma-ray bursts. However, due to the lack of measurements of $^{14}$C by year, the occurrence frequency of such $^{14}$C rapidly increasing events is poorly known. In addition, rapid increases may be hidden in the IntCal13 data with five-year resolution. Here we report the result of $^{14}$C measurements using an ancient buried tree during the period between BC 3388 and 3358. We find a rapid increase of about 9\textperthousand~ in the $^{14}$C content from BC 3372 to BC 3371. We suggest that this event could originate from a large solar proton event.


title: Confinement Effect Emission from Infiltrated ZnO in PS-b-PMMA Nanostructures

abstract: We have characterized the growth of ZnO using sequential infiltration synthesis (SiS) on PS-b-PMMA block copolymers (BCP) of spherical and cylindrical sub-20nm morphologies and studied how the photoluminescence of these nanostructures varies per its seed layer. Investigation of these structures was done using atomic force microscopy (AFM), spectrofluorometry, Raman spectroscopy, and scanning electron microscopy (SEM). We report blue-shifted photoemission at 335 nm (3.70 eV), suggesting quantum confinement effects. This UV-photoluminescence can be translated into emitter sizes of roughly 1.5nm in radius. Furthermore, samples of ZnO prepared with an alumina seed layer showed additional defect state photoemission at 470 nm and 520 nm for spherical and cylindrical BCP morphologies, respectively. Defect photoemission was not observed in samples prepared without a seed layer. Raman and EXAFS data suggest lack of long range order between the ZnO nanostructures during early stages of infiltrated ZnO growth and therefore supports the blue shift emission is due to confinement. Our work demonstrates that ZnO nanostructures grown on PS-b-PMMA via infiltration are advantageous in uniformity and size, and exhibit unique fluorescence properties. These observations suggest that infiltrated ZnO in PS-b-PMMA nanostructures lends itself to a new regime of applications in photonics and quantum materials.


title: Poisson-Boltzmann theory with non-linear ion correlations

abstract: The Poisson-Boltzmann (PB) theory is widely used to depict ionic systems. As a mean-field theory, the PB theory neglects the correlation effect in the ionic atmosphere and leads to deviations from experimental results as the concentration or charge valance increases. A modified PB theory including ion correlation effect while retaining its simplicity is critical for many important applications in which ion correlation effect can be significant. In this paper, we present a new model to incorporate ion correlations into the original PB equation by utilizing the Green's function with a non-linear form of the self-energy, which is different from the linear self-energy equation obtained by the Field-Theoretic (FT) approach. Both equations are solved numerically and compared with our molecular dynamics (MD) simulation. The co-ion distribution calculated by the FT approach deviates significantly from the MD simulation, while our results for both counter-ion and co-ion distributions are justified by the MD simulation.


title: Modelling stochastic resonance in humans: the influence of lapse rate

abstract: Adding noise to a sensory signal generally decreases human performance. However noise can improve performance too, due to a process called stochastic resonance (SR). This paradoxical effect may be exploited in psychophysical experiments, to provide additional insights into how the sensory system deals with noise. Here, I develop a model for stochastic resonance to study the influence of noise on human perception, in which the biological parameter of `lapse rate' was included. I show that the inclusion of lapse rate allows for the occurrence of stochastic resonance in terms of the performance metric d'. At the same time, I show that high levels of lapse rate cause stochastic resonance to disappear. It is also shown that noise generated in the brain (i.e., internal noise) may obscure any effect of stochastic resonance in experimental settings. I further relate the model to a standard equivalent noise model, the linear amplifier model, and show that the lapse rate can function to scale the threshold versus noise (TvN) curve, similar to the efficiency parameter in equivalent noise (EN) models. Therefore, lapse rate provides a psychophysical explanation for reduced efficiency in EN paradigms. Furthermore, I note that ignoring lapse rate may lead to an overestimation of internal noise in equivalent noise paradigms. Overall, describing stochastic resonance in terms of signal detection theory, with the inclusion of lapse rate, may provide valuable new insights into how human performance depends on internal and external noise.


title: From Word Segmentation to POS Tagging for Vietnamese

abstract: This paper presents an empirical comparison of two strategies for Vietnamese Part-of-Speech (POS) tagging from unsegmented text: (i) a pipeline strategy where we consider the output of a word segmenter as the input of a POS tagger, and (ii) a joint strategy where we predict a combined segmentation and POS tag for each syllable. We also make a comparison between state-of-the-art (SOTA) feature-based and neural network-based models. On the benchmark Vietnamese treebank (Nguyen et al., 2009), experimental results show that the pipeline strategy produces better scores of POS tagging from unsegmented text than the joint strategy, and the highest accuracy is obtained by using a feature-based model.


title: Sparse High-Dimensional Linear Regression. Algorithmic Barriers and a Local Search Algorithm

abstract: We consider a sparse high dimensional regression model where the goal is to recover a $k$-sparse unknown vector $\beta^*$ from $n$ noisy linear observations of the form $Y=X\beta^*+W \in \mathbb{R}^n$ where $X \in \mathbb{R}^{n \times p}$ has iid $N(0,1)$ entries and $W \in \mathbb{R}^n$ has iid $N(0,\sigma^2)$ entries. Under certain assumptions on the parameters, an intriguing assymptotic gap appears between the minimum value of $n$, call it $n^*$, for which the recovery is information theoretically possible, and the minimum value of $n$, call it $n_{\mathrm{alg}}$, for which an efficient algorithm is known to provably recover $\beta^*$. In \cite{gamarnikzadik} it was conjectured that the gap is not artificial, in the sense that for sample sizes $n \in [n^*,n_{\mathrm{alg}}]$ the problem is algorithmically hard. We support this conjecture in two ways. Firstly, we show that the optimal solution of the LASSO provably fails to $\ell_2$-stably recover the unknown vector $\beta^*$ when $n \in [n^*,c n_{\mathrm{alg}}]$, for some sufficiently small constant $c>0$. Secondly, we establish that $n_{\mathrm{alg}}$, up to a multiplicative constant factor, is a phase transition point for the appearance of a certain Overlap Gap Property (OGP) over the space of $k$-sparse vectors. The presence of such an Overlap Gap Property phase transition, which originates in statistical physics, is known to provide evidence of an algorithmic hardness. Finally we show that if $n>C n_{\mathrm{alg}}$ for some large enough constant $C>0$, a very simple algorithm based on a local search improvement rule is able both to $\ell_2$-stably recover the unknown vector $\beta^*$ and to infer correctly its support, adding it to the list of provably successful algorithms for the high dimensional linear regression problem.


title: Chandra Observation of PWN G16.73+0.08 in SNR G16.7+0.1

abstract: We present X-ray observations of PWN G16.73+0.08/SNR G16.7+0.1 using archival data of {\it Chandra} ACIS. The X-ray emission peak location of this pulsar wind nebula is found to be offset by 24 arcsec from the centre of the 1.4-GHz emission of this nebula. The X-ray nebula is elongated in the direction from the X-ray peak to the 1.4-GHz emission centre. This offset suggests that G16.73+0.08 is an evolved pulsar wind nebula interacting with the supernova remnant reverse shock. We identify a point source, CXO J182058.16-142001.5, near the location of the X-ray peak. The spectrum of the X-ray nebula can be described by an absorbed power law of photon index $0.98^{+0.79}_{-0.71}$ and hydrogen column density $N_{\rm H}=4.99^{+2.75}_{-2.28}\times 10^{22}$ cm$^{-2}$. CXO J182058.16-142001.5 is likely a pulsar. We estimate its spin-down power to be about $2.6\times 10 ^{36}$ erg s$^{-1}$. Assuming its age at 3000 and 10,000 years, its dipole magnetic field strength at the polar surface is estimated to be about $4.2 \times 10^{13}$ G and $1.3 \times 10^{13}$ G, respectively.


title: Impartial Triangular Chocolate Bar Games

abstract: Chocolate bar games are variants of the game of Nim in which the goal is to leave your opponent with the single bitter part of the chocolate bar. The rectangular chocolate bar game is a thinly disguised form of classical multi-heap Nim. In this work, we investigate the mathematical structure of triangular chocolate bar games in which the triangular chocolate bar can be cut in three directions. In the triangular chocolate bar game, a position is a $\mathcal{P}$-position if and only if $x \oplus y \oplus z = 0$, where the numbers $x,y,z$ stand for the maximum number of times that the chocolate bar can be cut in each direction. Moreover, the Grundy number of a position $(x,y,z)$ is not always equal to $x \oplus y \oplus z $, and a generic formula for Grundy numbers in not known. Therefore, the mathematical structure of triangular chocolate bar game is different from that of classical Nim.


title: Scalable Peaceman-Rachford Splitting Method with Proximal Terms

abstract: Along with developing of Peaceman-Rachford Splittling Method (PRSM), many batch algorithms based on it have been studied very deeply. But almost no algorithm focused on the performance of stochastic version of PRSM. In this paper, we propose a new stochastic algorithm based on PRSM, prove its convergence rate in ergodic sense, and test its performance on both artificial and real data. We show that our proposed algorithm, Stochastic Scalable PRSM (SS-PRSM), enjoys the $O(1/K)$ convergence rate, which is the same as those newest stochastic algorithms that based on ADMM but faster than general Stochastic ADMM (which is $O(1/\sqrt{K})$). Our algorithm also owns wide flexibility, outperforms many state-of-the-art stochastic algorithms coming from ADMM, and has low memory cost in large-scale splitting optimization problems.


title: Classical Structured Prediction Losses for Sequence to Sequence Learning

abstract: There has been much recent work on training neural attention models at the sequence-level using either reinforcement learning-style methods or by optimizing the beam. In this paper, we survey a range of classical objective functions that have been widely used to train linear models for structured prediction and apply them to neural sequence to sequence models. Our experiments show that these losses can perform surprisingly well by slightly outperforming beam search optimization in a like for like setup. We also report new state of the art results on both IWSLT'14 German-English translation as well as Gigaword abstractive summarization. On the larger WMT'14 English-French translation task, sequence-level training achieves 41.5 BLEU which is on par with the state of the art.


title: New inequalities for operator concave functions involving positive linear maps

abstract: The purpose of this paper is to present some general inequalities for operator concave functions which include some known inequalities as a particular case. Among other things, we prove that if $A\in \mathcal{B}\left( \mathcal{H} \right)$ is a positive operator such that $mI\le A\le MI$ for some scalars $0<m<M$ and $\Phi $ is a normalized positive linear map on $\mathcal{B}\left( \mathcal{H} \right)$, then \[\begin{aligned} {{\left( \frac{M+m}{2\sqrt{Mm}} \right)}^{r}}&\ge {{\left( \frac{\frac{1}{\sqrt{Mm}}\Phi \left( A \right)+\sqrt{Mm}\Phi \left( {{A}^{-1}} \right)}{2} \right)}^{r}} & \ge \frac{\frac{1}{{{\left( Mm \right)}^{\frac{r}{2}}}}\Phi {{\left( A \right)}^{r}}+{{\left( Mm \right)}^{\frac{r}{2}}}\Phi {{\left( {{A}^{-1}} \right)}^{r}}}{2} & \ge \Phi {{\left( A \right)}^{r}}\sharp\Phi {{\left( {{A}^{-1}} \right)}^{r}}, \end{aligned}\] where $0\le r\le 1$, which nicely extend the operator Kantorovich inequality.


title: Flat maps to and from Noetherian rings

abstract: We investigate flat maps where the source or target is a Noetherian ring, giving necessary and/or sufficient conditions on a ring for such maps to exist. Along the way, we develop some general facts about flat ring maps, and exhibit many examples, including a new class of zero-dimensional local rings.


title: Compactified Cosmological Simulations of the Infinite Universe

abstract: We present a novel $N$-body simulation method that compactifies the infinite spatial extent of the Universe into a finite sphere with isotropic boundary conditions to follow the evolution of the large-scale structure. Our approach eliminates the need for periodic boundary conditions, a mere numerical convenience which is not supported by observation and which modifies the law of force on large scales in an unrealistic fashion. We demonstrate that our method outclasses standard simulations executed on workstation-scale hardware in dynamic range, it is balanced in following a comparable number of high and low $k$ modes and, its fundamental geometry and topology match observations. Our approach is also capable of simulating an expanding, infinite universe in static coordinates with Newtonian dynamics. The price of these achievements is that most of the simulated volume has smoothly varying mass and spatial resolution, an approximation that carries different systematics than periodic simulations. Our initial implementation of the method is called StePS which stands for Stereographically Projected Cosmological Simulations. It uses stereographic projection for space compactification and naive $\mathcal{O}(N^2)$ force calculation which is nevertheless faster to arrive at a correlation function of the same quality than any standard (tree or P$^3$M) algorithm with similar spatial and mass resolution. The $N^2$ force calculation is easy to adapt to modern graphics cards, hence our code can function as a high-speed prediction tool for modern large-scale surveys. To learn about the limits of the respective methods, we compare StePS with GADGET-2 \citep{Gadget2_2005MNRAS.364.1105S} running matching initial conditions.


title: Wythoff's Game with a Pass

abstract: The author introduces the Wythoff game with a pass, a variant of the classical Wythoff game. The latter game is played with two piles of stones, and two players take turns removing stones from one or both piles; when removing stones from both piles, the numbers of stones removed from each pile must be equal. The player who removes the last stone or stones is the winner. An equivalent description of the game is that a single chess queen is placed somewhere on a large grid of squares, and each player can move the queen towards the upper left corner of the grid: leftward, upward, or to the upper left, any number of steps. The winner is the player who moves the queen into the upper left corner. In Wythoff's game with a pass, we modify the standard rules of the game so as to allow for a one-time pass, i.e., a pass move which may be used at most once in a game, and not from a terminal position. Once the pass has been used by either player, it is no longer available. The author discovered that, when a pass move is available, the set of $\mathcal{P}$-positions (previsous player's position) of the Wythoff's game with a pass is equal to $ T_0 \cup T_1\cup B -A$, where $T_0 = \{(x,y): \mathcal{G}(x,y)=0 \}$, $T_0 = \{(x,y): \mathcal{G}(x,y)=1 \}$, $B$ is a finite subset of $\{(x,y): \mathcal{G}(x,y)=4 \}$, $A$ is a finite subset of $T_0$ and $\mathcal{G}$ is the Grundy number of the classical Wythoff game. Here the size of $B$ is $6$, and the size of $A$ is $7$. According to the research of U.Blass and A.S.Fraenkel, for a position $(x,y) \in T_1 $ there is always a position $(x^{\prime},y^{\prime}) \in T_0$ such that $\mid x- x^{\prime} \mid \leq 2$ and $\mid y - y^{\prime} \mid \leq 4$, and hence the graph of the set of $\mathcal{P}$-positions of the Wythoff's game with a pass is very similar to the graph of the set of $\mathcal{P}$-positions of the classical Wythoff's game.


title: The Descartes circles theorem and division by zero calculus

abstract: From the viewpoint of the division by zero $(0/0=1/0=z/0=0)$ and the division by zero calculus, we will show that in the very beautiful theorem by Descartes on three touching circles is valid for lines and points for circles except for one case. However, for the exceptional case, we can obtain interesting results from the division by zero calculus.


title: Emergence of a Dimer-Dimer Interaction in the Low-Energy Effective Quantum-Dimer Model of a Diamond-Like-Decorated Square-Lattice Heisenberg Antiferromagnets with Further Neighbor Couplings

abstract: We study spin-1/2 Heisenberg antiferromagnets on a diamond-like-decorated square lattice perturbed by two kinds of further neighbor couplings. In our previous study [J. Phys. Soc. Jpn. 85, 094002 (2016)], the second-order effective Hamiltonian for the Heisenberg model perturbed by a further neighbor coupling was found to be a square-lattice quantum-dimer model with a finite hopping amplitude, t>0, and no dimer-dimer interaction, v=0. In this study, we introduce another kind of further neighbor coupling and show that it leads to an attractive interaction between dimers, which suggests the stabilization of the columnar phase of the square-lattice quantum-dimer model. The calculated v/t is presented as a function of the ratio of the two exchange parameters in the Heisenberg model.


title: Electronic orders near the type-II van Hove singularity in BC$_3$

abstract: Using the functional renormalization group, we investigate the electron instability in the single-sheet BC$_3$ when the electron filling is near a type-II van Hove singularity. For a finite Hubbard interaction, the ferromagnetic-like spin density wave order dominates in the immediate vicinity of the singularity. Elsewhere near the singularity the p-wave superconductivity prevails. We also find that a small nearest-neighbor Coulomb repulsion can enhance the superconductivity. Our results show that BC$_3$ would be a promising candidate to realize topological $p+ip'$ superconductivity, but the transition temperature is practically sizable only if the local interaction is moderately strong.


title: Dynamic Fusion Networks for Machine Reading Comprehension

abstract: This paper presents a novel neural model - Dynamic Fusion Network (DFN), for machine reading comprehension (MRC). DFNs differ from most state-of-the-art models in their use of a dynamic multi-strategy attention process, in which passages, questions and answer candidates are jointly fused into attention vectors, along with a dynamic multi-step reasoning module for generating answers. With the use of reinforcement learning, for each input sample that consists of a question, a passage and a list of candidate answers, an instance of DFN with a sample-specific network architecture can be dynamically constructed by determining what attention strategy to apply and how many reasoning steps to take. Experiments show that DFNs achieve the best result reported on RACE, a challenging MRC dataset that contains real human reading questions in a wide variety of types. A detailed empirical analysis also demonstrates that DFNs can produce attention vectors that summarize information from questions, passages and answer candidates more effectively than other popular MRC models.


title: Near-optimal sample complexity for convex tensor completion

abstract: We analyze low rank tensor completion (TC) using noisy measurements of a subset of the tensor. Assuming a rank-$r$, order-$d$, $N \times N \times \cdots \times N$ tensor where $r=O(1)$, the best sampling complexity that was achieved is $O(N^{\frac{d}{2}})$, which is obtained by solving a tensor nuclear-norm minimization problem. However, this bound is significantly larger than the number of free variables in a low rank tensor which is $O(dN)$. In this paper, we show that by using an atomic-norm whose atoms are rank-$1$ sign tensors, one can obtain a sample complexity of $O(dN)$. Moreover, we generalize the matrix max-norm definition to tensors, which results in a max-quasi-norm (max-qnorm) whose unit ball has small Rademacher complexity. We prove that solving a constrained least squares estimation using either the convex atomic-norm or the nonconvex max-qnorm results in optimal sample complexity for the problem of low-rank tensor completion. Furthermore, we show that these bounds are nearly minimax rate-optimal. We also provide promising numerical results for max-qnorm constrained tensor completion, showing improved recovery results compared to matricization and alternating least squares.


title: Coherent wave transmission in quasi-one-dimensional systems with L\'evy disorder

abstract: We study the random fluctuations of the transmission in disordered quasi-one-dimensional systems such as disordered waveguides and/or quantum wires whose random configurations of disorder are characterized by density distributions with a long tail known as L\'evy distributions. The presence of L\'evy disorder leads to large fluctuations of the transmission and anomalous localization, in relation to the standard exponential localization (Anderson localization). We calculate the complete distribution of the transmission fluctuations for different number of transmission channels in the presence and absence of time-reversal symmetry. Significant differences in the transmission statistics between disordered systems with Anderson and anomalous localizations are revealed. The theoretical predictions are independently confirmed by tight binding numerical simulations.


title: A proof of the Jacobian conjecture

abstract: Based on the results people have obtained, we try to prove the Jacobian conjecture, but there is a gap in the proof.


title: Experimental realization of Feynman's ratchet

abstract: Feynman's ratchet is a microscopic machine in contact with two heat reservoirs, at temperatures $T_A$ and $T_B$, that was proposed by Richard Feynman to illustrate the second law of thermodynamics. In equilibrium ($T_A=T_B$), thermal fluctuations prevent the ratchet from generating directed motion. When the ratchet is maintained away from equilibrium by a temperature difference ($T_A \ne T_B$), it can operate as a heat engine, rectifying thermal fluctuations to perform work. While it has attracted much interest, the operation of Feynman's ratchet as a heat engine has not been realized experimentally, due to technical challenges. In this work, we realize Feynman's ratchet with a colloidal particle in a one dimensional optical trap in contact with two heat reservoirs: one is the surrounding water, while the effect of the other reservoir is generated by a novel feedback mechanism, using the Metropolis algorithm to impose detailed balance. We verify that the system does not produce work when $T_A=T_B$, and that it becomes a microscopic heat engine when $T_A \ne T_B$. We analyze work, heat and entropy production as functions of the temperature difference and external load. Our experimental realization of Feynman's ratchet and the Metropolis algorithm can also be used to study the thermodynamics of feedback control and information processing, the working mechanism of molecular motors, and controllable particle transportation.


title: Straggler Mitigation in Distributed Optimization Through Data Encoding

abstract: Slow running or straggler tasks can significantly reduce computation speed in distributed computation. Recently, coding-theory-inspired approaches have been applied to mitigate the effect of straggling, through embedding redundancy in certain linear computational steps of the optimization algorithm, thus completing the computation without waiting for the stragglers. In this paper, we propose an alternate approach where we embed the redundancy directly in the data itself, and allow the computation to proceed completely oblivious to encoding. We propose several encoding schemes, and demonstrate that popular batch algorithms, such as gradient descent and L-BFGS, applied in a coding-oblivious manner, deterministically achieve sample path linear convergence to an approximate solution of the original problem, using an arbitrarily varying subset of the nodes at each iteration. Moreover, this approximation can be controlled by the amount of redundancy and the number of nodes used in each iteration. We provide experimental results demonstrating the advantage of the approach over uncoded and data replication strategies.


title: Vortex dipole precession in a Fermi superfluid with Rashba spin-orbit coupling

abstract: We present a closed-form expression for the time dynamics of a balanced Fermi superfluid with s-wave pairing and Rashba spin-orbit (SO) coupling. Solving the associated self-consistency conditions for an initial number density given by a Gaussian and an initial gap containing a vortex dipole, we show that Rashba SO coupling results in precession of the vortex dipole. The integrated self-consistent gap decays as a function of time, characterising a "transient Fermi superfluid". Our analytic solution forms a starting point for studies of vortex dynamics -- in particular braiding of Majorana fermions -- in Fermi superfluids in the presence of spin population imbalance and two-dimensional SO coupling, which are both required to realize a topological atomic Fermi gas in experiments.


title: DataVizard: Recommending Visual Presentations for Structured Data

abstract: Selecting the appropriate visual presentation of the data such that it preserves the semantics of the underlying data and at the same time provides an intuitive summary of the data is an important, often the final step of data analytics. Unfortunately, this is also a step involving significant human effort starting from selection of groups of columns in the structured results from analytics stages, to the selection of right visualization by experimenting with various alternatives. In this paper, we describe our \emph{DataVizard} system aimed at reducing this overhead by automatically recommending the most appropriate visual presentation for the structured result. Specifically, we consider the following two scenarios: first, when one needs to visualize the results of a structured query such as SQL; and the second, when one has acquired a data table with an associated short description (e.g., tables from the Web). Using a corpus of real-world database queries (and their results) and a number of statistical tables crawled from the Web, we show that DataVizard is capable of recommending visual presentations with high accuracy. We also present the results of a user survey that we conducted in order to assess user views of the suitability of the presented charts vis-a-vis the plain text captions of the data.


title: Superconductivity in Sr$_2$RuO$_4$ thin film under biaxial strain

abstract: Motivated by the success of experimental manipulation of the band structure through biaxial strain in Sr$_2$RuO$_4$ thin film grown on a mismatched substrate, we investigate theoretically the effects of biaxial strain on the electronic instabilities, such as superconductivity (SC) and spin density wave (SDW), by functional renormalization group. According to the experiment, the positive strain (from lattice expansion) causes charge transfer to the $\gamma$-band and consequently Lifshitz reconstruction of the Fermi surface. Our theoretical calculations show that within a limited range of positive strain a p-wave superconducting order is realized. However, as the strain is increased further the system develops into the SDW state well before the Lifshitz transition is reached. We also consider the effect of negative strains (from lattice constriction). As the strain increases, there is a transition from p-wave SC state to nodal s-wave SC state. The theoretical results are discussed in comparison to experiment and can be checked by further experiments.


title: A Robust Variable Step Size Fractional Least Mean Square (RVSS-FLMS) Algorithm

abstract: In this paper, we propose an adaptive framework for the variable step size of the fractional least mean square (FLMS) algorithm. The proposed algorithm named the robust variable step size-FLMS (RVSS-FLMS), dynamically updates the step size of the FLMS to achieve high convergence rate with low steady state error. For the evaluation purpose, the problem of system identification is considered. The experiments clearly show that the proposed approach achieves better convergence rate compared to the FLMS and adaptive step-size modified FLMS (AMFLMS).


title: Modeling and Performance Comparison of Privacy Approaches for Location Based Services

abstract: In pervasive computing environment, Location Based Services (LBSs) are getting popularity among users because of their usefulness in day-to-day life. LBSs are information services that use geospatial data of mobile device and smart phone users to provide information, entertainment and security in real time. A key concern in such pervasive computing environment is the need to reveal the user's exact location which may allow an adversary to infer private information about the user. To address the privacy concerns of LBS users, a large number of security approaches have been proposed based on the concept of k-anonymity. The central idea in location k-anonymity is to find a set of k-1 users confined in a given geographical area of the actual user, such that the location of these k users are indistinguishable from one another, thus protecting the identity of the user. Although a number of performance parameters like success rate, amount of privacy achieved are used to measure the performance of the k-anonymity approaches, they make the implicit, unrealistic assumption that the k-1 users are readily available. As such these approaches ignore the turnaround time to process a user request, which is crucial for a real-time application like LBS. In this work, we model the k-anonymity approaches using queuing theory to compute the average sojourn time of users and queue length of the system. To demonstrate that queuing theory can be used to model all k-anonymity approaches, we consider graph-based k-anonymity approaches. The proposed analytical model is further validated with experimental results.


title: Study on a Spinorial Representation of Linear Canonical Transformations

abstract: This work is a continuation of our previous works concerning linear canonical transformations and phase space representation of quantum theory. It is mainly focused on the description of an approach which allows to establish spinorial representation of linear canonical transformations. This description is started with the presentation of a suitable parameterization of linear canonical transformations which permits to represent them with special pseudo-orthogonal transformations in an operator space. Then the establishment of the spinorial representation is deduced using the well-known relation existing between special pseudo-orthogonal and spin groups. The cases of one dimension and general multidimensional theory are both studied.


title: Molecular Configuration Around Single Vacancy in Solid CO2 at T = 0 K Studied by Monte Carlo Simulation Technique

abstract: Energetically minimum configurations of the first- and second-nearest neighbor (NN) CO2 molecules surrounding a vacancy, created by removing a single molecule from the Pa3 structure, were calculated by use of the Monte Carlo simulation technique at zero temperature. It was found that among the 1st NN molecules, only six NN molecules, closest to the oxygen atoms of the removed CO2 molecule, significantly changed the positions and orientations, while the other six molecules were not influenced by the creation of the vacancy. The configurations of the molecules no longer had the three-fold symmetry that the Pa3 structure possessed. The influence of the vacancy on the 2nd NN molecules was relatively small.


title: Approximation by generalized bivariate Kantorovich sampling type series

abstract: The purpose of this paper is to construct a bivariate generalization of new family of Kantorovich type sampling operators $(K_w^{\varphi}f)_{w>0}.$ First, we give the pointwise convergence theorem and a Voronovskaja type theorem for these Kantorovich generalized sampling series. Further, we obtain the degree of approximation by means of modulus of continuity and quantitative version of Voronovskaja type theorem for the family $(K_w^{\varphi}f)_{w>0}.$ Finally, we give some examples of kernels such as box spline kernels and Bochner-Riesz kernel to which the theory can be applied.


title: Robust Online Speed Scaling With Deadline Uncertainty

abstract: A speed scaling problem is considered, where time is divided into slots, and jobs with payoff $v$ arrive at the beginning of the slot with associated deadlines $d$. Each job takes one slot to be processed, and multiple jobs can be processed by the server in each slot with energy cost $g(k)$ for processing $k$ jobs in one slot. The payoff is accrued by the algorithm only if the job is processed by its deadline. We consider a robust version of this speed scaling problem, where a job on its arrival reveals its payoff $v$, however, the deadline is hidden to the online algorithm, which could potentially be chosen adversarially and known to the optimal offline algorithm. The objective is to derive a robust (to deadlines) and optimal online algorithm that achieves the best competitive ratio. We propose an algorithm (called min-LCR) and show that it is an optimal online algorithm for any convex energy cost function $g(.)$. We do so without actually evaluating the optimal competitive ratio, and give a general proof that works for any convex $g$, which is rather novel. For the popular choice of energy cost function $g(k) = k^\alpha, \alpha \ge 2$, we give concrete bounds on the competitive ratio of the algorithm, which ranges between $2.618$ and $3$ depending on the value of $\alpha$. The best known online algorithm for the same problem, but where deadlines are revealed to the online algorithm has competitive ratio of $2$ and a lower bound of $\sqrt{2}$. Thus, importantly, lack of deadline knowledge does not make the problem degenerate, and the effect of deadline information on the optimal competitive ratio is limited.


title: Quantum transport senses community structure in networks

abstract: Quantum time evolution exhibits rich physics, attributable to the interplay between the density and phase of a wave function. However, unlike classical heat diffusion, the wave nature of quantum mechanics has not yet been extensively explored in modern data analysis. We propose that the Laplace transform of quantum transport (QT) can be used to construct an ensemble of maps from a given complex network to a circle $S^1$, such that closely-related nodes on the network are grouped into sharply concentrated clusters on $S^1$. The resulting QT clustering (QTC) algorithm is as powerful as the state-of-the-art spectral clustering in discerning complex geometric patterns and more robust when clusters show strong density variations or heterogeneity in size. The observed phenomenon of QTC can be interpreted as a collective behavior of the microscopic nodes that evolve as macroscopic cluster orbitals in an effective tight-binding model recapitulating the network. Python source code implementing the algorithm and examples are available at https://github.com/jssong-lab/QTC.

categories: quant-ph cond-mat.other cs.DS q-bio.QM stat.ML


title: Hadronic vacuum polarization contribution to the anomalous magnetic moments of leptons from first principles

abstract: We compute the leading, strong-interaction contribution to the anomalous magnetic moment of the electron, muon and tau using lattice quantum chromodynamics (QCD) simulations. Calculations include the effects of $u$, $d$, $s$ and $c$ quarks and are performed directly at the physical values of the quark masses and in volumes of linear extent larger than $6\,\mathrm{fm}$. All connected and disconnected Wick contractions are calculated. Continuum limits are carried out using six lattice spacings. We obtain $a_e^\mathrm{LO-HVP}=189.3(2.6)(5.6)\times 10^{-14}$, $a_\mu^\mathrm{LO-HVP}=711.1(7.5)(17.4)\times 10^{-10}$ and $a_\tau^\mathrm{LO-HVP}=341.0(0.8)(3.2)\times 10^{-8}$, where the first error is statistical and the second is systematic.


title: SkipFlow: Incorporating Neural Coherence Features for End-to-End Automatic Text Scoring

abstract: Deep learning has demonstrated tremendous potential for Automatic Text Scoring (ATS) tasks. In this paper, we describe a new neural architecture that enhances vanilla neural network models with auxiliary neural coherence features. Our new method proposes a new \textsc{SkipFlow} mechanism that models relationships between snapshots of the hidden representations of a long short-term memory (LSTM) network as it reads. Subsequently, the semantic relationships between multiple snapshots are used as auxiliary features for prediction. This has two main benefits. Firstly, essays are typically long sequences and therefore the memorization capability of the LSTM network may be insufficient. Implicit access to multiple snapshots can alleviate this problem by acting as a protection against vanishing gradients. The parameters of the \textsc{SkipFlow} mechanism also acts as an auxiliary memory. Secondly, modeling relationships between multiple positions allows our model to learn features that represent and approximate textual coherence. In our model, we call this \textit{neural coherence} features. Overall, we present a unified deep learning architecture that generates neural coherence features as it reads in an end-to-end fashion. Our approach demonstrates state-of-the-art performance on the benchmark ASAP dataset, outperforming not only feature engineering baselines but also other deep learning models.


title: Anomalous dispersion relations in the staggered flux state

abstract: We study the quasiparticle properties in the d-wave superconductivity, antiferromagnetism, and the staggered flux state within a renormalized mean-field theory based on the two-dimensional t-J model. In particular, we focus on the anomalous quasiparticle dispersion relations around the antinodal region of the Brillouin zone argued in Hashimoto et al., Nature Phys. 6, 414 (2010). We obtain the qualitatively consistent results with their observations when we take account of the SF state. The present analysis shows that the SF order can be a possible candidate of symmetry-breaking pseudogap states coexisting with the dSC.


title: Multiplication structure of the cohomology ring of real toric spaces

abstract: A real toric space is a topological space which admits a well-behaved $\mathbb{Z}_2^k$-action. Real moment-angle complexes and real toric varieties are typical examples of real toric spaces. A real toric space is determined by a pair of a simplicial complex $K$ and a characteristic matrix $\Lambda$. In this paper, we provide an explicit $R$-cohomology ring formula of a real toric space in terms of $K$ and $\Lambda$, where $R$ is a commutative ring with unity in which $2$ is a unit. Interestingly, it has a natural $(\mathbb{Z} \oplus \operatorname*{row} \Lambda)$-grading. As corollaries, we compute the cohomology rings of (generalized) real Bott manifolds in terms of binary matroids, and we also provide a criterion for real toric spaces to be cohomology symplectic.


title: Solid solution decomposition and Guinier-Preston zone formation in Al-Cu alloys: A kinetic theory with anisotropic interactions

abstract: Using methods of statistical kinetic theory parametrized with first-principles interatomic interactions that include chemical and strain contributions, we investigated the kinetics of decomposition and microstructure formation in Al-Cu alloys as a function of temperature and alloy concentration. We show that the decomposition of the solid solution forming platelets of copper, known as Guinier-Preston (GP) zones, includes several stages and that the transition from GP1 to GP2 zones is determined mainly by kinetic factors. With increasing temperature, the model predicts a gradual transition from platelet-like precipitates to equiaxial ones and at intermediate temperatures both precipitate morphologies may coexist.


title: Equidistribution of closed geodesics along random walk trajectories with respect to the harmonic invariant measure

abstract: We prove that for suitable random walks on isometry groups of $CAT(-1)$ spaces, typical sample paths eventually land on loxodromic elements which equidistribute with respect to a flow invariant measure on the unit tangent bundle of the quotient space.


title: A generalised Euler-Poincar\'e formula for associahedra

abstract: We derive a formula for the number of flip-equivalence classes of tilings of an $n$-gon by collections of tiles of shape dictated by an integer partition $\lambda$. The proof uses the Euler-Poincar\'e formula; and the formula itself generalises the Euler-Poincar\'e formula for associahedra.


title: Unified Pragmatic Models for Generating and Following Instructions

abstract: We show that explicit pragmatic inference aids in correctly generating and following natural language instructions for complex, sequential tasks. Our pragmatics-enabled models reason about why speakers produce certain instructions, and about how listeners will react upon hearing them. Like previous pragmatic models, we use learned base listener and speaker models to build a pragmatic speaker that uses the base listener to simulate the interpretation of candidate descriptions, and a pragmatic listener that reasons counterfactually about alternative descriptions. We extend these models to tasks with sequential structure. Evaluation of language generation and interpretation shows that pragmatic inference improves state-of-the-art listener models (at correctly interpreting human instructions) and speaker models (at producing instructions correctly interpreted by humans) in diverse settings.


title: Concurrent Pump Scheduling and Storage Level Optimization Using Meta-Models and Evolutionary Algorithms

abstract: In spite of the growing computational power offered by the commodity hardware, fast pump scheduling of complex water distribution systems is still a challenge. In this paper, the Artificial Neural Network (ANN) meta-modeling technique has been employed with a Genetic Algorithm (GA) for simultaneously optimizing the pump operation and the tank levels at the ends of the cycle. The generalized GA+ANN algorithm has been tested on a real system in the UK. Comparing to the existing operation, the daily cost is reduced by about 10-15%, while the number of pump switches are kept below 4 switches-per-day. In addition, tank levels are optimized ensure a periodic behavior, which results in a predictable and stable performance over repeated cycles.


title: Morphology of Critically-Sized Crystalline Nuclei at Shear-Induced Crystal Nucleation in Amorphous Solid

abstract: In this work we study morphological characteristics of the critically-sized crystalline nuclei at initial stage of the shear-induced crystallization of a model single-component amorphous (glassy) system. These characteristics are estimated quantitatively through statistical treatment of the non-equilibrium molecular dynamics simulation results for the system under steady shear at various (fixed) values of the shear rate $\dot\gamma$ and at different temperatures. It is found that the sheared glassy system is crystallized through nucleation mechanism. From analysis of a time-dependent trajectories of the largest crystalline nuclei, the critical size $n_{c}$ and the nucleation time $\tau_{c}$ were defined. It is shown that the critically-sized nuclei in the system are oriented within the shear-gradient $xy$-plane at moderate and high shear rates; and a tilt angle of the oriented nuclei depends on the shear rate. At extremely high shear rates and at shear deformation of the system more than $60$\%, the tilt angle of the nuclei tends to take the value $\simeq45^\circ$ respective to the shear direction. We found that this feature depends weakly on the temperature. Asphericity of the nucleus shape increases with increasing shear rate, that is verified by increasing value of the asphericity parameter and by the contour of the pair distribution function calculated for the particles of the critically-sized nuclei. The critical size increases with increasing shear rate according to the power-law, $n_{c}\propto(\dot\gamma\tau_{c})^{1/3}$, whereas the shape of the critically-sized nucleus changes from spherical to the elongated ellipsoidal. We found that the $n_{c}$-dependencies of the nuclei deformation parameter evaluated for the system at different temperatures and shear rates are collapsed into unified master-curve.

categories: cond-mat.mtrl-sci cond-mat.dis-nn cond-mat.mes-hall physics.app-ph physics.atm-clus


title: Strong consistency and optimality for generalized estimating equations with stochastic covariates

abstract: In this article we study the existence and strong consistency of GEE estimators, when the generalized estimating functions are martingales with random coefficients. Furthermore, we characterize estimating functions which are asymptotically optimal.


title: Anomaly free cosmological perturbations with generalised holonomy correction in loop quantum cosmology

abstract: In the spatially flat case of loop quantum cosmology, the connection $\bar{k}$ is usually replaced by the holonomy $\frac{\sin(\bar{\mu}k)}{\bar{\mu}}$ in the effective theory. In this paper, instead of the $\bar{\mu}$ scheme, we use a generalised, undertermined function $g(\bar{k},\bar{p})$ to represent the holonomy and by using the approach of anomaly free constraint algebra we fix all the counter terms in the constraints and find the restriction on the form of $g(\bar{k},\bar{p})$, then we derive the gauge invariant equations of motion of the scalar, tensor and vector perturbations and study the inflationary power spectra with generalised holonomy correction.


title: Feature importance scores and lossless feature pruning using Banzhaf power indices

abstract: Understanding the influence of features in machine learning is crucial to interpreting models and selecting the best features for classification. In this work we propose the use of principles from coalitional game theory to reason about importance of features. In particular, we propose the use of the Banzhaf power index as a measure of influence of features on the outcome of a classifier. We show that features having Banzhaf power index of zero can be losslessly pruned without damage to classifier accuracy. Computing the power indices does not require having access to data samples. However, if samples are available, the indices can be empirically estimated. We compute Banzhaf power indices for a neural network classifier on real-life data, and compare the results with gradient-based feature saliency, and coefficients of a logistic regression model with $L_1$ regularization.