title: Identifying Generalization Properties in Neural Networks

abstract: While it has not yet been proven, empirical evidence suggests that model generalization is related to local properties of the optima which can be described via the Hessian. We connect model generalization with the local property of a solution under the PAC-Bayes paradigm. In particular, we prove that model generalization ability is related to the Hessian, the higher-order "smoothness" terms characterized by the Lipschitz constant of the Hessian, and the scales of the parameters. Guided by the proof, we propose a metric to score the generalization capability of the model, as well as an algorithm that optimizes the perturbed model accordingly.


title: TESS in the Solar System

abstract: The Transiting Exoplanet Survey Satellite (TESS), launched successfully on 18th of April, 2018, will observe nearly the full sky and will provide time-series imaging data in ~27-day-long campaigns. TESS is equipped with 4 cameras; each has a field-of-view of 24x24 degrees. During the first two years of the primary mission, one of these cameras, Camera #1, is going to observe fields centered at an ecliptic latitude of 18 degrees. While the ecliptic plane itself is not covered during the primary mission, the characteristic scale height of the main asteroid belt and Kuiper belt implies that a significant amount of small solar system bodies will cross the field-of-view of this camera. Based on the comparison of the expected amount of information of TESS and Kepler/K2, we can compute the cumulative etendues of the two optical setups. This comparison results in roughly comparable optical etendues, however the net etendue is significantly larger in the case of TESS since all of the imaging data provided by the 30-minute cadence frames are downlinked rather than the pre-selected stamps of Kepler/K2. In addition, many principles of the data acquisition and optical setup are clearly different, including the level of confusing background sources, full-frame integration and cadence, the field-of-view centroid with respect to the apparent position of the Sun, as well as the differences in the duration of the campaigns. As one would expect, TESS will yield time-series photometry and hence rotational properties for only brighter objects, but in terms of spatial and phase space coverage, this sample will be more homogeneous and more complete. Here we review the main analogues and differences between the Kepler/K2 mission and the TESS mission, focusing on scientific implications and possible yields related to our Solar System.


title: Examples of badly approximable vectors over number fields

abstract: We consider approximation of vectors $\mathbf{z}\in F\otimes\mathbb{R}\cong\mathbb{R}^r\times\mathbb{C}^s$ by elements of a number field $F$ and construct examples of badly approximable vectors. These examples come from compact subspaces of $SL_2(\mathcal{O}_F)\backslash SL_2(F\otimes\mathbb{R})$ naturally associated to (totally indefinite, anisotropic) $F$-rational binary quadratic and Hermitian forms, a generalization of the well-known fact that quadratic irrationals are badly approximable over $\mathbb{Q}$.


title: Distances for WiFi Based Topological Indoor Mapping

abstract: For localization and mapping of indoor environments through WiFi signals, locations are often represented as likelihoods of the received signal strength indicator. In this work we compare various measures of distance between such likelihoods in combination with different methods for estimation and representation. In particular, we show that among the considered distance measures the Earth Mover's Distance seems the most beneficial for the localization task. Combined with kernel density estimation we were able to retain the topological structure of rooms in a real-world office scenario.


title: Exploiting Tournament Selection for Efficient Parallel Genetic Programming

abstract: Genetic Programming (GP) is a computationally intensive technique which is naturally parallel in nature. Consequently, many attempts have been made to improve its run-time from exploiting highly parallel hardware such as GPUs. However, a second methodology of improving the speed of GP is through efficiency techniques such as subtree caching. However achieving parallel performance and efficiency is a difficult task. This paper will demonstrate an efficiency saving for GP compatible with the harnessing of parallel CPU hardware by exploiting tournament selection. Significant efficiency savings are demonstrated whilst retaining the capability of a high performance parallel implementation of GP. Indeed, a 74% improvement in the speed of GP is achieved with a peak rate of 96 billion GPop/s for classification type problems.


title: Unfolding the complexity of the global value chain: Strengths and entropy in the single-layer, multiplex, and multi-layer international trade networks

abstract: The worldwide trade network has been widely studied through different data sets and network representations with a view to better understanding interactions among countries and products. Here we investigate international trade through the lenses of the single-layer, multiplex, and multi-layer networks. We discuss differences among the three network frameworks in terms of their relative advantages in capturing salient topological features of trade. We draw on the World Input-Output Database to build the three networks. We then uncover sources of heterogeneity in the way strength is allocated among countries and transactions by computing the strength distribution and entropy in each network. Additionally, we trace how entropy evolved, and show how the observed peaks can be associated with the onset of the global economic downturn. Findings suggest how more complex representations of trade, such as the multi-layer network, enable us to disambiguate the distinct roles of intra- and cross-industry transactions in driving the evolution of entropy at a more aggregate level. We discuss our results and the implications of our comparative analysis of networks for research on international trade and other empirical domains across the natural and social sciences.


title: Egocentric Vision-based Future Vehicle Localization for Intelligent Driving Assistance Systems

abstract: Predicting the future location of vehicles is essential for safety-critical applications such as advanced driver assistance systems (ADAS) and autonomous driving. This paper introduces a novel approach to simultaneously predict both the location and scale of target vehicles in the first-person (egocentric) view of an ego-vehicle. We present a multi-stream recurrent neural network (RNN) encoder-decoder model that separately captures both object location and scale and pixel-level observations for future vehicle localization. We show that incorporating dense optical flow improves prediction results significantly since it captures information about motion as well as appearance change. We also find that explicitly modeling future motion of the ego-vehicle improves the prediction accuracy, which could be especially beneficial in intelligent and automated vehicles that have motion planning capability. To evaluate the performance of our approach, we present a new dataset of first-person videos collected from a variety of scenarios at road intersections, which are particularly challenging moments for prediction because vehicle trajectories are diverse and dynamic.


title: New Scalar Field Quartessence

abstract: We propose a cosmological scenario involving a scalar field, $\varphi$, that is a source of Dark Matter as well as of Dark Energy. Besides $\varphi$, the Lagrangian of the field theory envisaged in our scenario contains a second field $\chi$, for simplicity assumed to be a scalar, too. For fixed values of $\chi$, the potential term decays exponentially at large positive values of $\varphi$. While $\varphi$ is not coupled to Standard Model fields, $\chi$ is assumed to be coupled to them, and the Green functions of $\chi$ depend on the cosmological redshift in the expanding universe. We assume that the term in the Lagrangian coupling $\chi$ to $\varphi$ is such that, at redshifts $z$ larger than some critical redshift $z_c$, $\varphi$ is trapped near $\varphi=0$, and oscillations of $\varphi$ about $\varphi=0$ describing massive scalar particles give rise to Dark Matter. At redshifts below $z_c$, the field $\varphi$ is no longer trapped near the origin and starts to `roll' towards large field values. A homogenous component of $\varphi$ emerges that acts as Dark Energy. Within over-dense regions, such as galaxies and galaxy clusters, the redshifting of $\chi$ stops, and $\varphi$ therefore remains trapped near $\varphi=0$ as long as $z_c$ is smaller than the redshift when structures on galactic scales decouple from the Hubble flow. Thus, at the present time, $\varphi$ describes both Dark Energy and Dark Matter.


title: Faster super-resolution imaging with auto-correlation two-step deconvolution

abstract: Despite super-resolution fluorescence blinking microscopes break the diffraction limit, the intense phototoxic illumination and long-term image sequences thus far still pose to major challenges in visualizing live-organisms. Here, we proposed a super-resolution method based on auto-correlation two-step deconvolution (SACD) to enhance the temporal resolution at lower signal intensity levels. Unlike conventional techniques, such as super-resolution optical fluctuation imaging (SOFI) or stochastic optical reconstruction microscopy (STORM), our model allows 16 frames to generate super-resolution images, without noticeable degradation in recording quality. We demonstrate SACD both in simulated predictions and experimental validations, with the resulting spatial resolution of 64 nm and 2-10 fold speed improvements. The use of low signal-to-noise ratio acquisition of image sequences, our SACD enables fast, background-free, super-resolution imaging that maybe advance implemented as a suitable tool for rapid live-cells imaging.


title: Operating Nanobeams in a Quantum Fluid

abstract: Microelectromechanical (MEMS) and nanoelectromechanical systems (NEMS) are ideal candidates for exploring quantum fluids since they can be manufactured reproducibly, cover the frequency range from hundreds of kilohertz up to gigahertz and usually have very low power dissipation. Their small size offers the possibility of probing the superfluid on scales comparable to, and below, the coherence length. That said, there have been hitherto no successful measurements of NEMS resonators in the liquid phases of helium. Here we report the operation of doubly-clamped aluminum nanobeams in superfluid $^4$He at temperatures spanning the superfluid transition. The devices are shown to be very sensitive detectors of the superfluid density and the normal fluid damping. However, a further and very important outcome of this work is the knowledge that now we have demonstrated that these devices can be successfully operated in superfluid $^4$He, it is straightforward to apply them in superfluid $^3$He which can be routinely cooled to below 100\,$\mu$K. This brings us into the regime where nanomechanical devices operating at a few MHz frequencies may enter their mechanical quantum ground state.


title: Learning, Planning, and Control in a Monolithic Neural Event Inference Architecture

abstract: We introduce REPRISE, a REtrospective and PRospective Inference SchEme, which learns temporal event-predictive models of dynamical systems. REPRISE infers the unobservable contextual event state and accompanying temporal predictive models that best explain the recently encountered sensorimotor experiences retrospectively. Meanwhile, it optimizes upcoming motor activities prospectively in a goal-directed manner. Here, REPRISE is implemented by a recurrent neural network (RNN), which learns temporal forward models of the sensorimotor contingencies generated by different simulated dynamic vehicles. The RNN is augmented with contextual neurons, which enable the encoding of distinct, but related, sensorimotor dynamics as compact event codes. We show that REPRISE concurrently learns to separate and approximate the encountered sensorimotor dynamics: it analyzes sensorimotor error signals adapting both internal contextual neural activities and connection weight values. Moreover, we show that REPRISE can exploit the learned model to induce goal-directed, model-predictive control, that is, approximate active inference: Given a goal state, the system imagines a motor command sequence optimizing it with the prospective objective to minimize the distance to the goal. The RNN activities thus continuously imagine the upcoming future and reflect on the recent past, optimizing the predictive model, the hidden neural state activities, and the upcoming motor activities. As a result, event-predictive neural encodings develop, which allow the invocation of highly effective and adaptive goal-directed sensorimotor control.


title: The Role of Mesons in Muon g-2

abstract: The muon anomaly $a_\mu=(g_\mu-2)/2$ showing a persisting 3 to 4 $\sigma$ deviation between the SM prediction and the experiment is one of the most promising signals for physics beyond the SM. As is well known, the hadronic uncertainties are limiting the accuracy of the Standard Model prediction. Therefore a big effort is going on to improve the evaluations of hadronic effects in order to keep up with the 4-fold improved precision expected from the new Fermilab measurement in the near future. A novel complementary type experiment planned at J-PARC in Japan, operating with ultra cold muons, is expected to be able to achieve the same accuracy but with completely different systematics. So exciting times in searching for New Physics are under way. I discuss the role of meson physics in calculations of the hadronic part of the muon g-2. The improvement is expected to substantiate the present deviation $\Delta a_\mu^{\rm New \ Physics}=\Delta a_\mu^{\rm Experiment}- \Delta a_\mu^{\rm Standard \ Model}$ to a 6 to 10 standard deviation effect, provided hadronic uncertainties can be reduce by a factor two. This concerns the hadronic vacuum polarization as well as the hadronic light-by-light scattering contributions, both to a large extent determined by the low lying meson spectrum. Better meson production data and progress in modeling meson form factors could greatly help to improve the precision and reliability of the SM prediction of $a_\mu$ and thereby provide more information on what is missing in the SM.


title: Dimension rigidity for cookie-cutter Cantor sets

abstract: We show that two cookie-cutter Cantor sets with the same symbolic coding are differentiably equivalent if and only if their Hausdorff dimensions are equal.


title: Bias corrected minimum distance estimator for short and long memory processes

abstract: This work proposes a new minimum distance estimator (MDE) for the parameters of short and long memory models. This bias corrected minimum distance estimator (BCMDE) considers a correction in the usual MDE to account for the bias of the sample autocorrelation function when the mean is unknown. We prove the weak consistency of the BCMDE for the general fractional autoregressive moving average (ARFIMA(p, d, q)) model and derive its asymptotic distribution for some particular cases. Simulation studies show that the BCMDE presents a good performance compared to other procedures frequently used in the literature, such as the maximum likelihood estimator, the Whittle estimator and the MDE. The results also show that the BCMDE presents, in general, the smallest mean squared error and is less biased than the MDE when the mean is a non-trivial function of time.


title: Evidence for an $\eta_c(1S) \pi^-$ resonance in $B^0 \to \eta_c(1S) K^+\pi^-$ decays

abstract: A Dalitz plot analysis of $B^0 \to \eta_c(1S) K^+\pi^-$ decays is performed using data samples of $pp$ collisions collected with the LHCb detector at centre-of-mass energies of $\sqrt{s}=7,~8$ and $13$ TeV, corresponding to a total integrated luminosity of $4.7~\text{fb}^{-1}$. A satisfactory description of the data is obtained when including a contribution representing an exotic $\eta_c(1S) \pi^-$ resonant state. The significance of this exotic resonance is more than three standard deviations, while its mass and width are $4096 \pm 20~^{+18}_{-22}$ MeV and $152 \pm 58~^{+60}_{-35}$ MeV, respectively. The spin-parity assignments $J^P=0^+$ and $J^{P}=1^-$ are both consistent with the data. In addition, the first measurement of the $B^0 \to \eta_c(1S) K^+\pi^-$ branching fraction is performed and gives $\displaystyle \mathcal{B}(B^0 \to \eta_c(1S) K^+\pi^-) = (5.73 \pm 0.24 \pm 0.13 \pm 0.66) \times 10^{-4}$, where the first uncertainty is statistical, the second systematic, and the third is due to limited knowledge of external branching fractions.


title: Deep Part Induction from Articulated Object Pairs

abstract: Object functionality is often expressed through part articulation -- as when the two rigid parts of a scissor pivot against each other to perform the cutting function. Such articulations are often similar across objects within the same functional category. In this paper, we explore how the observation of different articulation states provides evidence for part structure and motion of 3D objects. Our method takes as input a pair of unsegmented shapes representing two different articulation states of two functionally related objects, and induces their common parts along with their underlying rigid motion. This is a challenging setting, as we assume no prior shape structure, no prior shape category information, no consistent shape orientation, the articulation states may belong to objects of different geometry, plus we allow inputs to be noisy and partial scans, or point clouds lifted from RGB images. Our method learns a neural network architecture with three modules that respectively propose correspondences, estimate 3D deformation flows, and perform segmentation. To achieve optimal performance, our architecture alternates between correspondence, deformation flow, and segmentation prediction iteratively in an ICP-like fashion. Our results demonstrate that our method significantly outperforms state-of-the-art techniques in the task of discovering articulated parts of objects. In addition, our part induction is object-class agnostic and successfully generalizes to new and unseen objects.


title: Loss asymmetries in quantum travelling wave parametric amplifiers

abstract: We study theoretically how loss impacts the amplification and squeezing performance of a generic quantum travelling wave parametric amplifier. Unlike previous studies, we analyze how having different levels of loss at signal and idler frequencies can dramatically alter properties compared to the case of frequency-independent loss. We find that loss asymmetries increase the amplifier's added noise in comparison to the symmetric loss case. More surprisingly, even small levels of loss asymmetry can completely destroy any quantum squeezing of symmetric collective output quadratures, while nonetheless leaving the output state strongly entangled.


title: Randomization Tests for Weak Null Hypotheses in Randomized Experiments

abstract: The Fisher randomization test (FRT) is appropriate for any test statistic, under a sharp null hypothesis that can recover all missing potential outcomes. However, it is often sought after to test a weak null hypothesis that the treatment does not affect the units on average. To use the FRT for a weak null hypothesis, we must address two issues. First, we need to impute the missing potential outcomes although the weak null hypothesis cannot determine all of them. Second, we need to choose a proper test statistic. For a general weak null hypothesis, we propose an approach to imputing missing potential outcomes under a compatible sharp null hypothesis. Building on this imputation scheme, we advocate a studentized statistic. The resulting FRT has multiple desirable features. First, it is model-free. Second, it is finite-sample exact under the sharp null hypothesis that we use to impute the potential outcomes. Third, it conservatively controls large-sample type I errors under the weak null hypothesis of interest. Therefore, our FRT is agnostic to the treatment effect heterogeneity. We establish a unified theory for general factorial experiments. We also extend it to stratified and clustered experiments.


title: Time dependent fracture under unloading in a fiber bundle model

abstract: We investigate the fracture of heterogeneous materials occurring under unloading from an initial load. Based on a fiber bundle model of time dependent fracture, we show that depending on the unloading rate the system has two phases: for rapid unloading the system suffers only partial failure and it has an infinite lifetime, while at slow unloading macroscopic failure occurs in a finite time. The transition between the two phases proved to be analogous to continuous phase transitions. Computer simulations revealed that during unloading the fracture proceeds in bursts of local breakings triggered by slowly accumulating damage. In both phases the time evolution starts with a relaxation of the bursting activity characterized by a universal power law decay of the burst rate. In the phase of finite lifetime the initial slowdown is followed by an acceleration towards macroscopic failure where the increasing rate of bursts obeys the (inverse) Omori law of earthquakes. We pointed out a strong correlation between the time where the event rate reaches a minimum value and of the lifetime of the system which allows for forecasting of the imminent catastrophic failure.


title: Supersingular Elliptic Curves and Moonshine

abstract: We generalize a theorem of Ogg on supersingular $j$-invariants to supersingular elliptic curves with level. Ogg observed that the level one case yields a characterization of the primes dividing the order of the monster. We show that the corresponding analyses for higher levels give analogous characterizations of the primes dividing the orders of other sporadic simple groups (e.g., baby monster, Fischer's largest group). This situates Ogg's theorem in a broader setting. More generally, we characterize, in terms of supersingular elliptic curves with level, the primes arising as orders of Fricke elements in centralizer subgroups of the monster. We also present a connection between supersingular elliptic curves and umbral moonshine. Finally, we present a procedure for explicitly computing invariants of supersingular elliptic curves with level structure.


title: From Kondo effect to weak-link regime in quantum spin-1/2 spin chains

abstract: We analyze the crossover from Kondo to weak-link regime by means of a model of tunable bond impurities in the middle of a spin-1/2 XXZ Heisenberg chain. We study the Kondo screening cloud and estimate the Kondo length by combining perturbative renormalization group approach with the exact numerical calculation of the integrated real-space spin-spin correlation functions. We show that, when the spin impurity is symmetrically coupled to the two parts of the chain with realistic values of the Kondo coupling strengths and spin-parity symmetry is preserved, the Kondo length takes values within the reach of nowadays experimental technology in ultracold-atom setups. In the case of non-symmetric Kondo couplings and/or spin parity broken by a nonzero magnetic field applied to the impurity, we discuss how Kondo screening redistributes among the chain as a function of the asymmetry in the couplings and map out the shrinking of the Kondo length when the magnetic field induces a crossover from Kondo impurity to weak-link physics.


title: Metacommutation in Central Simple Algebras

abstract: In a quaternion order of class number one, an element can be factored in multiple ways depending on the order of the factorization of its reduced norm. The fact that multiplication is not commutative causes an element to induce a permutation on the set of primes of a given reduced norm. We discuss this permutation and previously known results about the cycle structure, sign, and number of fixed points for quaternion orders. We generalize these results to other orders in central simple algebras over global fields.


title: Towards Accountable AI: Hybrid Human-Machine Analyses for Characterizing System Failure

abstract: As machine learning systems move from computer-science laboratories into the open world, their accountability becomes a high priority problem. Accountability requires deep understanding of system behavior and its failures. Current evaluation methods such as single-score error metrics and confusion matrices provide aggregate views of system performance that hide important shortcomings. Understanding details about failures is important for identifying pathways for refinement, communicating the reliability of systems in different settings, and for specifying appropriate human oversight and engagement. Characterization of failures and shortcomings is particularly complex for systems composed of multiple machine learned components. For such systems, existing evaluation methods have limited expressiveness in describing and explaining the relationship among input content, the internal states of system components, and final output quality. We present Pandora, a set of hybrid human-machine methods and tools for describing and explaining system failures. Pandora leverages both human and system-generated observations to summarize conditions of system malfunction with respect to the input content and system architecture. We share results of a case study with a machine learning pipeline for image captioning that show how detailed performance views can be beneficial for analysis and debugging.


title: Mean Estimation with Sub-Gaussian Rates in Polynomial Time

abstract: We study polynomial time algorithms for estimating the mean of a heavy-tailed multivariate random vector. We assume only that the random vector $X$ has finite mean and covariance. In this setting, the radius of confidence intervals achieved by the empirical mean are large compared to the case that $X$ is Gaussian or sub-Gaussian. We offer the first polynomial time algorithm to estimate the mean with sub-Gaussian-size confidence intervals under such mild assumptions. Our algorithm is based on a new semidefinite programming relaxation of a high-dimensional median. Previous estimators which assumed only existence of finitely-many moments of $X$ either sacrifice sub-Gaussian performance or are only known to be computable via brute-force search procedures requiring time exponential in the dimension.


title: Personalized Top-N Sequential Recommendation via Convolutional Sequence Embedding

abstract: Top-$N$ sequential recommendation models each user as a sequence of items interacted in the past and aims to predict top-$N$ ranked items that a user will likely interact in a `near future'. The order of interaction implies that sequential patterns play an important role where more recent items in a sequence have a larger impact on the next item. In this paper, we propose a Convolutional Sequence Embedding Recommendation Model (\emph{Caser}) as a solution to address this requirement. The idea is to embed a sequence of recent items into an `image' in the time and latent spaces and learn sequential patterns as local features of the image using convolutional filters. This approach provides a unified and flexible network structure for capturing both general preferences and sequential patterns. The experiments on public datasets demonstrated that Caser consistently outperforms state-of-the-art sequential recommendation methods on a variety of common evaluation metrics.


title: Congruences on infinite partition and partial Brauer monoids

abstract: We give a complete description of the congruences on the partition monoid $P_X$ and the partial Brauer monoid $PB_X$, where $X$ is an arbitrary infinite set, and also of the lattices formed by all such congruences. Our results complement those from a recent article of East, Mitchell, Ruskuc and Torpey, which deals with the finite case. As a consequence of our classification result, we show that the congruence lattices of $P_X$ and $PB_X$ are isomorphic to each other, and are distributive and well quasi-ordered. We also calculate the smallest number of pairs of partitions required to generate any congruence; when this number is infinite, it depends on the cofinality of certain limit cardinals.


title: Ranking Distillation: Learning Compact Ranking Models With High Performance for Recommender System

abstract: We propose a novel way to train ranking models, such as recommender systems, that are both effective and efficient. Knowledge distillation (KD) was shown to be successful in image recognition to achieve both effectiveness and efficiency. We propose a KD technique for learning to rank problems, called \emph{ranking distillation (RD)}. Specifically, we train a smaller student model to learn to rank documents/items from both the training data and the supervision of a larger teacher model. The student model achieves a similar ranking performance to that of the large teacher model, but its smaller model size makes the online inference more efficient. RD is flexible because it is orthogonal to the choices of ranking models for the teacher and student. We address the challenges of RD for ranking problems. The experiments on public data sets and state-of-the-art recommendation models showed that RD achieves its design purposes: the student model learnt with RD has a model size less than half of the teacher model while achieving a ranking performance similar to the teacher model and much better than the student model learnt without RD.


title: On truncations of the Chalker-Coddington model

abstract: The supersymmetric reformulation of physical observables in the Chalker-Coddington model (CC) for the plateau transition in the integer quantum Hall effect leads to a reformulation of its critical properties in terms of a 2D non-compact loop model or a 1D non-compact $gl(2|2)$ spin chain. Following a proposal by Ikhlef, Fendley and Cardy, we define and study a series of truncations of these loop models and spin chains, involving a finite and growing number of degrees of freedom per site. The case of the first truncation is solved analytically using the Bethe-ansatz. It is shown to exhibit many of the qualitative features expected for the untruncated theory, including a quadratic spectrum of exponents with a continuous component, and a normalizable ground state below that continuum. Quantitative properties are however at odds with the results of simulations on the CC model. Higher truncations are studied only numerically. While their properties are found to get closer to those of the CC model, it is not clear whether this is a genuine effect, or the result of strong finite-size corrections.


title: CRN++: Molecular Programming Language

abstract: Synthetic biology is a rapidly emerging research area, with expected wide-ranging impact in biology, nanofabrication, and medicine. A key technical challenge lies in embedding computation in molecular contexts where electronic micro-controllers cannot be inserted. This necessitates effective representation of computation using molecular components. While previous work established the Turing-completeness of chemical reactions, defining representations that are faithful, efficient, and practical remains challenging. This paper introduces CRN++, a new language for programming deterministic (mass-action) chemical kinetics to perform computation. We present its syntax and semantics, and build a compiler translating CRN++ programs into chemical reactions, thereby laying the foundation of a comprehensive framework for molecular programming. Our language addresses the key challenge of embedding familiar imperative constructs into a set of chemical reactions happening simultaneously and manipulating real-valued concentrations. Although some deviation from ideal output value cannot be avoided, we develop methods to minimize the error, and implement error analysis tools. We demonstrate the feasibility of using CRN++ on a suite of well-known algorithms for discrete and real-valued computation. CRN++ can be easily extended to support new commands or chemical reaction implementations, and thus provides a foundation for developing more robust and practical molecular programs.


title: Dwarf Galaxies as Cosmological Probes

abstract: The Lambda Cold Dark Matter (LCDM) paradigm makes specific predictions for the abundance, structure, substructure and clustering of dark matter halos, the sites of galaxy formation. These predictions can be directly tested, in the low-mass halo regime, by dark matter-dominated dwarf galaxies. A number of potential challenges to LCDM have been identified when confronting the expected properties of dwarfs with observation. I review our understanding of a few of these issues, including the `missing satellites' and the `too-big-to-fail' problems, and argue that neither poses an insurmountable challenge to LCDM. Solving these problems requires that most dwarf galaxies inhabit halos of similar mass, and that there is a relatively sharp minimum halo mass threshold to form luminous galaxies. These predictions are eminently falsifiable. In particular, LCDM predicts a large number of `dark' low-mass halos, some of which should have retained enough primordial gas to be detectable in deep 21 cm or H$_\alpha$ surveys. Detecting this predicted population of `mini-halos' would be a major discovery and a resounding success for LCDM on small scales.


title: Optimal transport with discrete long range mean field interactions

abstract: We study an optimal transport problem where, at some intermediate time, the mass is accelerated by either an external force field, or self-interacting. We obtain regularity of the velocity potential, intermediate density, and optimal transport map, under conditions on the interaction potential that are related to the so-called Ma-Trudinger-Wang condition from optimal transport.


title: On the Cheng-Yau gradient estimate for Carnot groups and sub-Riemannian manifolds

abstract: In this note we show how results in \cite{BaudoinBonnefont2016, BaudoinGarofalo2013, CoulhonJiangKoskelaSikora2017} yield the Cheng-Yau estimate on two classes of sub-Riemannian manifolds: Carnot groups and sub-Riemannian manifolds satisfying a generalized curvature-dimension inequality.


title: Hurwitz Theory of Elliptic Orbifolds, II

abstract: An elliptic orbifold is the quotient of an elliptic curve by a finite group. In 2001, Eskin and Okounkov proved that generating functions for the number of branched covers of an elliptic curve with specified ramification are quasimodular forms for $SL_2(\mathbb{Z}).$ In 2006, they generalized this theorem to the enumeration of branched covers of the quotient of an elliptic curve by $\pm 1$, proving quasi-modularity for $\Gamma_1(2)$. In 2017, the author generalized their work to the quotient of an elliptic curve by $\langle \zeta_N\rangle$ for $N=3, 4, 6$, proving quasimodularity for $\Gamma_1(N)$. In these works, both Eskin-Okounkov and the author had to assume that there was at least one orbifold point of order $N$ over which there was no ramification. Here we remove that assumption, with the caveat that the generating functions are only quasimodular for $\Gamma(N)$. We deduce the following corollary: Let $h_6(\vec{\kappa},q)$ be the generating function whose $q^n$ coefficient is the number of surface triangulations with $2n$ triangles, such that the set of non-zero curvatures is $\kappa_i$. Here the curvature of a vertex is six minus its valence. Then under the substitution $q=e^{2\pi i \tau}$, the function $h_6(\vec{\kappa},q)$ is a quasimodular form for $\Gamma_1(6)$ with weight bounded in terms of $\vec{\kappa}$. This statement in turn implies that the Masur-Veech volume of any stratum of sextic differentials is polynomial in $\pi$.


title: Predicting Periodicity with Temporal Difference Learning

abstract: Temporal difference (TD) learning is an important approach in reinforcement learning, as it combines ideas from dynamic programming and Monte Carlo methods in a way that allows for online and incremental model-free learning. A key idea of TD learning is that it is learning predictive knowledge about the environment in the form of value functions, from which it can derive its behavior to address long-term sequential decision making problems. The agent's horizon of interest, that is, how immediate or long-term a TD learning agent predicts into the future, is adjusted through a discount rate parameter. In this paper, we introduce an alternative view on the discount rate, with insight from digital signal processing, to include complex-valued discounting. Our results show that setting the discount rate to appropriately chosen complex numbers allows for online and incremental estimation of the Discrete Fourier Transform (DFT) of a signal of interest with TD learning. We thereby extend the types of knowledge representable by value functions, which we show are particularly useful for identifying periodic effects in the reward sequence.


title: Deep Generative Classifiers for Thoracic Disease Diagnosis with Chest X-ray Images

abstract: Thoracic diseases are very serious health problems that plague a large number of people. Chest X-ray is currently one of the most popular methods to diagnose thoracic diseases, playing an important role in the healthcare workflow. However, reading the chest X-ray images and giving an accurate diagnosis remain challenging tasks for expert radiologists. With the success of deep learning in computer vision, a growing number of deep neural network architectures were applied to chest X-ray image classification. However, most of the previous deep neural network classifiers were based on deterministic architectures which are usually very noise-sensitive and are likely to aggravate the overfitting issue. In this paper, to make a deep architecture more robust to noise and to reduce overfitting, we propose using deep generative classifiers to automatically diagnose thorax diseases from the chest X-ray images. Unlike the traditional deterministic classifier, a deep generative classifier has a distribution middle layer in the deep neural network. A sampling layer then draws a random sample from the distribution layer and input it to the following layer for classification. The classifier is generative because the class label is generated from samples of a related distribution. Through training the model with a certain amount of randomness, the deep generative classifiers are expected to be robust to noise and can reduce overfitting and then achieve good performances. We implemented our deep generative classifiers based on a number of well-known deterministic neural network architectures, and tested our models on the chest X-ray14 dataset. The results demonstrated the superiority of deep generative classifiers compared with the corresponding deep deterministic classifiers.


title: Rapid change of superconductivity and electron-phonon coupling through 19% doping in Bi2212

abstract: Electron-boson coupling plays a key role in superconductivity for many systems. However, in copper-based high-temperature ($T_c$) superconductors, its relation to superconductivity remains controversial despite strong spectroscopic fingerprints. Here we use angle-resolved photoemission spectroscopy to find a striking correlation between the superconducting gap and the bosonic coupling strength near the Brillouin zone boundary in Bi$_2$Sr$_2$CaCu$_2$O$_{8+\delta}$. The bosonic coupling strength rapidly increases from the overdoped Fermi-liquid regime to the optimally doped strange metal, concomitant with the quadrupled superconducting gap and the doubled gap-to-Tc ratio across the pseudogap boundary. This synchronized lattice and electronic response suggests that the effects of electronic interaction and the electron-phonon coupling become intimately entangled upon entering the strange metal regime, which may in turn drive a stronger superconductivity.


title: Morphological, nanostructural, and compositional evolution during phase separation of a model Ni-Al-Mo superalloy: Atom-probe tomographic experiments and lattice-kinetic Monte Carlo simulations

abstract: The details of phase separation of a Ni-6.5Al-9.9Mo aged at 978 K for aging times ranging from 0.125 to 1024 h are investigated by atom-probe tomography and lattice-kinetic Monte Carlo (LKMC) simulations. On the basis of the temporal evolution of the nanostructure, three experimental regimes are identified: (1) concomitant precipitate nucleation and growth (t less than 0.25 h); (2) concurrent coagulation and coalescence (t 0.25 to 16 h); and (3) quasi-stationary coarsening of gamma prime(L12)- precipitates (t 16 to 1024 h). The temporal dependencies of the mean precipitate radius and precipitate number density, Nv(t), are determined experimentally, 0.344 (0.012) and -0.95(0.02), respectively, following the predictions of quasi-stationary coarsening models. In this alloy aged at 978 K, Al partitions strongly to the gamma prime(L12)-phase with a partitioning coefficient 4.06(0.04), whereas Mo and Ni partition to the gamma(f.c.c.)-matrix with values of 0.61(0.01) and 0.90(0.01), respectively. In the quasi-stationary regime(t larger than 16h), the temporal exponents of the Al, Mo, and Ni supersaturations in both the gamma(f.c.c.)-matrix and gamma prime(L12)-precipitates are in reasonable agreement with a multi-component coarsening model's prediction of -0.33. Quantitative analyses of the edge-to-edge inter-precipitate distances demonstrate that coagulation and coalescence are consequences of the overlap of the diffusion fields surrounding the gamma prime(L12)-precipitates. Both 3-D APT and LKMC results demonstrate that the interfacial compositional width decrease with increasing values. And the interfacial compositional width at infinite ageing time are estimated to be 1.89(0.22) nm, 2.09(0.12) nm and 2.64(0.03) nm for Ni, Al and Mo, respectively.


title: Unraveling materials Berry curvature and Chern-Simons numbers from real-time evolution of Bloch states

abstract: Materials can be classified by the topological character of their electronic structure and, in this perspective, global attributes immune to local deformations have been discussed in terms of Berry curvature and Chern numbers. Except for instructional simple models, linear response theories have been ubiquitously employed in calculations of topological properties of real materials. Here we propose a completely different and versatile approach to get the topological characteristics of materials by calculating physical observables from the real-time evolving Bloch states: the cell-averaged current density reveals the anomalous velocities whose integration leads to the conductivity quantum. Results for prototypical cases are shown, including a spin-frozen valley-Hall and a quantum anomalous Hall insulator. The advantage of this method is best illustrated by the example of a quantum spin Hall insulator: the quantized spin Hall conductivity is straightforwardly obtained irrespective of the non-Abelian nature in its Berry curvature. Moreover, the method can be extended to the description of real observables in non-equilibrium states of topological materials.


title: Notes on quasi-Polish spaces

abstract: Quasi-Polish spaces were introduced by de Brecht as a possibly non-Hausdorff generalization of Polish spaces sharing many of their descriptive set-theoretic properties. We give a self-contained exposition of the basic theory of quasi-Polish spaces, based on their "logical" characterization as $\mathbf{\Pi}^0_2$ subspaces of countable powers of Sierpinski space, with several new proofs emphasizing this point of view as well as making more extensive use of Baire category techniques.


title: Distribution-Free Prediction Sets for Two-Layer Hierarchical Models

abstract: We consider the problem of constructing distribution-free prediction sets for data from two-layer hierarchical distributions. For iid data, prediction sets can be constructed using the method of conformal prediction. The validity of conformal prediction hinges on the exchangeability of the data, which does not hold when groups of observations come from distinct distributions, such as multiple observations on each patient in a medical database. We extend conformal methods to this hierarchical setting. We develop CDF pooling, single subsampling, and repeated subsampling approaches to construct prediction sets in unsupervised and supervised settings. We compare these approaches in terms of coverage and average set size. If asymptotic coverage is acceptable, we recommend CDF pooling for its balance between empirical coverage and average set size. If we desire coverage guarantees, then we recommend the repeated subsampling approach.


title: Quantum Computational Supremacy

abstract: The field of quantum algorithms aims to find ways to speed up the solution of computational problems by using a quantum computer. A key milestone in this field will be when a universal quantum computer performs a computational task that is beyond the capability of any classical computer, an event known as quantum supremacy. This would be easier to achieve experimentally than full-scale quantum computing, but involves new theoretical challenges. Here we present the leading proposals to achieve quantum supremacy, and discuss how we can reliably compare the power of a classical computer to the power of a quantum computer.


title: Derivations on almost complex manifolds

abstract: In this short note, we propose an unified method to derive formulas for derivations conjugated by exponential functions on an almost complex manifold. In v3, we corrected some mistakes in previous versions.


title: DynaSOAr: Accelerating Single-Method Multiple-Objects Applications on GPUs

abstract: Object-oriented programming (OOP) has long been regarded as too inefficient for SIMD high-performance computing, despite the fact that many important HPC applications have an inherent object structure. We discovered a broad subset of OOP that can be implemented efficiently on massively parallel SIMD accelerators. We call it Single-Method Multiple-Objects (SMMO), because parallelism is expressed by running a method on all objects of a type. To make fast GPU programming available to domain experts who are less experienced in GPU programming, we developed DynaSOAr, a CUDA framework for SMMO applications. DynaSOAr improves the usage of allocated memory with an SOA data layout and achieves low memory fragmentation through efficient management of free and allocated memory blocks with lock-free, hierarchical bitmaps.


title: DP-3-coloring of planar graphs without $4,9$-cycles and two cycles from $\{5,6,7,8\}$

abstract: A generalization of list-coloring, now known as DP-coloring, was recently introduced by Dvo\v{r}\'{a}k and Postle. Essentially, DP-coloring assigns an arbitrary matching between lists of colors at adjacent vertices, as opposed to only matching identical colors as is done for list-coloring. Several results on list-coloring of planar graphs have since been extended to the setting of DP-coloring. We note that list-coloring results do not always extend to DP-coloring results. Our main result in this paper is to prove that every planar graph without cycles of length $\{4, a, b, 9\}$ for $a, b \in \{6, 7, 8\}$ is DP-$3$-colorable, extending three existing results on $3$-choosability of planar graphs.


title: The Broad Absorption Line Tidal Disruption Event iPTF15af: Optical and Ultraviolet Evolution

abstract: We present multi-wavelength observations of the tidal disruption event (TDE) iPTF15af, discovered by the intermediate Palomar Transient Factory (iPTF) survey at redshift $z=0.07897$. The optical and ultraviolet (UV) light curves of the transient show a slow decay over five months, in agreement with previous optically discovered TDEs. It also has a comparable black-body peak luminosity of $L_{\rm{peak}} \approx 1.5 \times 10^{44}$ erg/s. The inferred temperature from the optical and UV data shows a value of (3$-$5) $\times 10^4$ K. The transient is not detected in X-rays up to $L_X < 3 \times 10^{42}$erg/s within the first five months after discovery. The optical spectra exhibit two distinct broad emission lines in the He II region, and at later times also H$\alpha$ emission. Additionally, emission from [N III] and [O III] is detected, likely produced by the Bowen fluorescence effect. UV spectra reveal broad emission and absorption lines associated with high-ionization states of N V, C IV, Si IV, and possibly P V. These features, analogous to those of broad absorption line quasars (BAL QSOs), require an absorber with column densities $N_{\rm{H}} > 10^{23}$ cm$^{-2}$. This optically thick gas would also explain the non-detection in soft X-rays. The profile of the absorption lines with the highest column density material at the largest velocity is opposite that of BAL QSOs. We suggest that radiation pressure generated by the TDE flare at early times could have provided the initial acceleration mechanism for this gas. Spectral UV line monitoring of future TDEs could test this proposal.


title: A Coupled Evolutionary Network for Age Estimation

abstract: Age estimation of unknown persons is a challenging pattern analysis task due to the lacking of training data and various aging mechanisms for different people. Label distribution learning-based methods usually make distribution assumptions to simplify age estimation. However, age label distributions are often complex and difficult to be modeled in a parameter way. Inspired by the biological evolutionary mechanism, we propose a Coupled Evolutionary Network (CEN) with two concurrent evolutionary processes: evolutionary label distribution learning and evolutionary slack regression. Evolutionary network learns and refines age label distributions in an iteratively learning way. Evolutionary label distribution learning adaptively learns and constantly refines the age label distributions without making strong assumptions on the distribution patterns. To further utilize the ordered and continuous information of age labels, we accordingly propose an evolutionary slack regression to convert the discrete age label regression into the continuous age interval regression. Experimental results on Morph, ChaLearn15 and MegaAge-Asian datasets show the superiority of our method.


title: Determination of Bending Angle of Light Deflection Subject to Possible Weak and Strong Quantum Gravity Effects

abstract: Explicit expressions for the bending angle of light deflection arising from phenomenologically deformed black-hole metrics, subject to possible weak and strong quantum gravity effects, respectively, are obtained, by a highly effective method. The accuracy and effectiveness of these expressions are then illustrated by numerically solving the differential equation governing the deflection angle directly in the weak quantum effect situation.


title: Small eigenvalues of closed Riemann surfaces for large genus

abstract: In this article we study the asymptotic behavior of small eigenvalues of Riemann surfaces for large genus. We show that for any positive integer $k$, as the genus $g$ goes to infinity, the smallest $k$-th eigenvalue of Riemann surfaces in any thick part of moduli space of Riemann surfaces of genus $g$ is uniformly comparable to $\frac{1}{g^2}$ in $g$. In the proof of the upper bound, for any constant $\epsilon>0$, we will construct a closed Riemann surface of genus $g$ in any $\epsilon$-thick part of moduli space such that it admits a pants decomposition whose boundary curves all have length equal to $\epsilon$, and the number of separating systole curves in this surface is uniformly comparable to $g$.


title: Tight Continuous-Time Reachtubes for Lagrangian Reachability

abstract: We introduce continuous Lagrangian reachability (CLRT), a new algorithm for the computation of a tight and continuous-time reachtube for the solution flows of a nonlinear, time-variant dynamical system. CLRT employs finite strain theory to determine the deformation of the solution set from time $t_i$ to time $t_{i+1}$. We have developed simple explicit analytic formulas for the optimal metric for this deformation; this is superior to prior work, which used semi-definite programming. CLRT also uses infinitesimal strain theory to derive an optimal time increment $h_i$ between $t_i$ and $t_{i+1}$, nonlinear optimization to minimally bloat (i.e., using a minimal radius) the state set at time $t_i$ such that it includes all the states of the solution flow in the interval $[t_i,t_{i+1}]$. We use $\delta$-satisfiability to ensure the correctness of the bloating. Our results on a series of benchmarks show that CLRT performs favorably compared to state-of-the-art tools such as CAPD in terms of the continuous reachtube volumes they compute.


title: Blue Large-Amplitude Pulsators (BLAPs): possible origin, evolutionary status, and nature of their pulsations

abstract: The Blue Large-Amplitude Pulsators (BLAPs) constitute a new class of pulsating stars. They are hot stars with effective temperatures of $T_{\rm eff}\sim 30\,000$ K and surface gravities of $\log g \sim 4.9$, that pulsate with periods in the range $\sim 20-40$ min. In Romero et al. (2018), we proposed that BLAPs are hot low-mass He-core pre-white dwarf (WD) stars that pulsate either in high-order non-radial $g$(gravity) modes or low-order radial modes, including the fundamental radial mode. The theoretical modes with periods in the observed range are unstable due to the $\kappa$ mechanism associated with the $Z$ bump in the opacity at $\log T \sim 5.25$. In this work, we extend the study of Romero et al. (2018) by assessing the rate of period changes of nonradial $g$ modes and radial modes and comparing them with the values measured for BLAPs, in an attempt to validate the proposed evolutionary scenario, and to discern whether the observed modes are high-order $g$ modes or radial modes.


title: Intrinsic flow structure and multifractality in two-dimensional bacterial turbulence

abstract: The active interaction between the bacteria and fluid generates turbulent structures even at zero Reynolds number. Velocity of such a flow obtained experimentally has been quantitatively investigated based on streamline segment analysis. There is a clear transition at about $16$ times of the organism body length separating two different scale regimes, which may be attributed to the different influence of the viscous effect. Surprisingly the scaling extracted from the streamline segment indicates the existence of scale similarity even at the zero Reynolds number limit. Moreover the multifractal feature can be quantitatively described via a lognormal formula with the Hurst number $H=0.76$ and the intermittency parameter $\mu=0.20$, which is coincidentally in agreement with the three-dimensional hydrodynamic turbulence result. The direction of cascade is measured via the filter-space technique. An inverse energy cascade is confirmed. For the enstrophy, a forward cascade is observed when $r/R\le 3$, and an inverse one is observed when $r/R>3$. Additionally, the lognormal statistics is verified for the coarse-grained energy dissipation and enstrophy, which supports the lognormal formula to fit the measured scaling exponent.


title: Uplink Resource Allocation for Multiple Access Computational Offloading (Extended Version)

abstract: The mobile edge computing framework offers the opportunity to reduce the energy that devices must expend to complete computational tasks. The extent of that energy reduction depends on the nature of the tasks, and on the choice of the multiple access scheme. In this paper, we first address the uplink communication resource allocation for offloading systems that exploit the full capabilities of the multiple access channel (FullMA). For indivisible tasks we provide a closed-form optimal solution of the energy minimization problem when a given set of users with different latency constraints are offloading, and a tailored greedy search algorithm for finding a good set of offloading users. For divisible tasks we develop a low-complexity algorithm to find a stationary solution. To highlight the impact of the choice of multiple access scheme, we also consider the TDMA scheme, which, in general, cannot exploit the full capabilities of the channel, and we develop low-complexity optimal resource allocation algorithms for indivisible and divisible tasks under that scheme. The energy reduction facilitated by FullMA is illustrated in our numerical experiments. Further, those results show that the proposed algorithms outperform existing algorithms in terms of energy consumption and computational cost.


title: Conv-TasNet: Surpassing Ideal Time-Frequency Magnitude Masking for Speech Separation

abstract: Single-channel, speaker-independent speech separation methods have recently seen great progress. However, the accuracy, latency, and computational cost of such methods remain insufficient. The majority of the previous methods have formulated the separation problem through the time-frequency representation of the mixed signal, which has several drawbacks, including the decoupling of the phase and magnitude of the signal, the suboptimality of time-frequency representation for speech separation, and the long latency in calculating the spectrograms. To address these shortcomings, we propose a fully-convolutional time-domain audio separation network (Conv-TasNet), a deep learning framework for end-to-end time-domain speech separation. Conv-TasNet uses a linear encoder to generate a representation of the speech waveform optimized for separating individual speakers. Speaker separation is achieved by applying a set of weighting functions (masks) to the encoder output. The modified encoder representations are then inverted back to the waveforms using a linear decoder. The masks are found using a temporal convolutional network (TCN) consisting of stacked 1-D dilated convolutional blocks, which allows the network to model the long-term dependencies of the speech signal while maintaining a small model size. The proposed Conv-TasNet system significantly outperforms previous time-frequency masking methods in separating two- and three-speaker mixtures. Additionally, Conv-TasNet surpasses several ideal time-frequency magnitude masks in two-speaker speech separation as evaluated by both objective distortion measures and subjective quality assessment by human listeners. Finally, Conv-TasNet has a significantly smaller model size and a shorter minimum latency, making it a suitable solution for both offline and real-time speech separation applications.


title: Surface Majorana Flat Bands in $j=\frac{3}{2}$ Superconductors with Singlet-Quintet Mixing

abstract: Recent experiments have revealed the evidence of nodal-line superconductivity in half-Heusler superconductors, e.g. YPtBi. Theories have suggested the topological nature of such nodal-line superconductivity and proposed the existence of surface Majorana flat bands on the (111) surface of half-Heusler superconductors. Due to the divergent density of states of the surface Majorana flat bands, the surface order parameter and the surface impurity play essential roles in determining the surface properties. In this work, we studied the effect of the surface order parameter and the surface impurity on the surface Majorana flat bands of half-Heusler superconductors based on the Luttinger model. To be specific, we consider the topological nodal-line superconducting phase induced by the singlet-quintet pairing mixing, classify all the possible translationally invariant order parameters for the surface states according to irreducible representations of $C_{3v}$ point group, and demonstrate that any energetically favorable order parameter needs to break time-reversal symmetry. We further discuss the energy splitting in the energy spectrum of surface Majorana flat bands induced by different order parameters and non-magnetic or magnetic impurities. We proposed that the splitting in the energy spectrum can serve as the fingerprint of the pairing symmetry and mean-field order parameters. Our theoretical prediction can be examined in the future scanning tunneling microscopy experiments.


title: Towards Discrete Solution: A Sparse Preserving Method for Correspondence Problem

abstract: Many problems of interest in computer vision can be formulated as a problem of finding consistent correspondences between two feature sets. Feature correspondence (matching) problem with one-to-one mapping constraint is usually formulated as an Integral Quadratic Programming (IQP) problem with permutation (or orthogonal) constraint. Since it is NP-hard, relaxation models are required. One main challenge for optimizing IQP matching problem is how to incorporate the discrete one-to-one mapping (permutation) constraint in its quadratic objective optimization. In this paper, we present a new relaxation model, called Sparse Constraint Preserving Matching (SPM), for IQP matching problem. SPM is motivated by our observation that the discrete permutation constraint can be well encoded via a sparse constraint. Comparing with traditional relaxation models, SPM can incorporate the discrete one-to-one mapping constraint straightly via a sparse constraint and thus provides a tighter relaxation for original IQP matching problem. A simple yet effective update algorithm has been derived to solve the proposed SPM model. Experimental results on several feature matching tasks demonstrate the effectiveness and efficiency of SPM method.


title: Taylor dispersion in two-dimensional bacterial turbulence

abstract: In this work, single particle dispersion was analyzed for a bacterial turbulence by retrieving the virtual Lagrangian trajectory via numerical integration of the Lagrangian equation. High-order displacement functions were calculated for cases with and without mean velocity effect. Two-regime power-law behavior for short and long time evolutions were identified experimentally, which were separated by the Lagrangian integral time. For the case with the mean velocity effect, the experimental Hurst numbers were determined to be $0.94$ and $0.97$ for short and long times evolutions, respectively. For the case without the mean velocity effect, the values were $0.88$ and $0.58$. Moreover, very weak intermittency correction was detected. All measured Hurst number were above $1/2$, the value of the normal diffusion, which verifies the super-diffusion behavior of living fluid. This behavior increases the efficiency of bacteria to obtain food.


title: One-particle Green's function of interacting two electrons using analytic solutions for a three-body problem: comparison with exact Kohn--Sham system

abstract: For a three-electron system with finite-strength interactions confined to a one-dimensional harmonic trap, we solve the Schroedinger equation analytically to obtain the exact solutions, from which we construct explicitly the simultaneous eigenstates of the energy and total spin for the first time. The solutions for the three-electron system allow us to derive analytic expressions for the exact one-particle Green's function (GF) for the corresponding two-electron system. We calculate the GF in frequency domain to examine systematically its behavior depending on the electronic interactions. We also compare the pole structure of non-interacting GF using the exact Kohn--Sham (KS) potential with that of the exact GF to find that the discrepancy of the energy gap between the KS system and the original system is larger for a stronger interaction. We perform numerical examination on the behavior of GFs in real space to demonstrate that the exact and KS GFs can have shapes quite different from each other. Our simple model will help to understand generic characteristics of interacting GFs.


title: Parity of the Partition Function p(n,k)

abstract: Let p(n, k) denote the number of partitions of n into parts less than or equal to k. We show several properties of this function modulo 2. First, we prove that for fixed positive integers k and m, p(n,k) is periodic modulo m. Using this, we are able to find lower and upper bounds for the number of odd values of the function for a fixed k.


title: Absolute moments in terms of characteristic functions

abstract: The absolute moments of probability distributions are much more complicated than conventional ones. By using a direct and simpler approach, we retreat P. L. Hsu's (1951, J. Chinese Math. Soc., Vol. 1, pp. 257-280) formulas in terms of the characteristic function (which have been ignored in the literature) and provide some new results as well. The case of nonnegative random variables is also investigated through both characteristic function and Laplace-Stieltjes transform. Besides, we prove that the distribution of a nonnegative random variable with a finite fractional moment can be completely determined by a proper subset of the translated fractional moments. This improves significantly P. Hall's (1983, Z. W., Vol. 62, 355-359) result for distributions on the right-half line.


title: A note on Castelnuovo-Mumford regularity and Hilbert coefficients

abstract: New upper and lower bounds on the Castelnuovo-Mumford regularity are given in terms of the Hilbert coefficients. Examples are provided to show that these bounds are in some sense nearly sharp.


title: Trojan Quantum Walks

abstract: We investigate the transport properties and entanglement between spin and position of one-dimensional quantum walks starting from a qubit over position states following a delta-like (local state) and Gaussian (delocalized state) distributions. We find out that if the initial state is delocalized enough and a NOT gate reflects this state backwards, then the interference pattern extinguishes the position dispersion without preventing the propagation of the state. This effect allows the creation of a Trojan wave packet, a non-spreading and non-stationary double-peak quantum state.


title: Data Shuffling in Wireless Distributed Computing via Low-Rank Optimization

abstract: Intelligent mobile platforms such as smart vehicles and drones have recently become the focus of attention for onboard deployment of machine learning mechanisms to enable low latency decisions with low risk of privacy breach. However, most such machine learning algorithms are both computation-and-memory intensive, which makes it highly difficult to implement the requisite computations on a single device of limited computation, memory, and energy resources. Wireless distributed computing presents new opportunities by pooling the computation and storage resources among devices. For low-latency applications, the key bottleneck lies in the exchange of intermediate results among mobile devices for data shuffling. To improve communication efficiency, we propose a co-channel communication model and design transceivers by exploiting the locally computed intermediate values as side information. A low-rank optimization model is proposed to maximize the achieved degrees-of-freedom (DoF) by establishing the interference alignment condition for data shuffling. Unfortunately, existing approaches to approximate the rank function fail to yield satisfactory performance due to the poor structure in the formulated low-rank optimization problem. In this paper, we develop an efficient DC algorithm to solve the presented low-rank optimization problem by proposing a novel DC representation for the rank function. Numerical experiments demonstrate that the proposed DC approach can significantly improve the communication efficiency whereas the achievable DoF almost remains unchanged when the number of mobile devices grows.


title: Powers of Monomial Ideals and Combinatorics

abstract: This is an exposition of some new results on associated primes and the depth of different kinds of powers of monomial ideals in order to show a deep connection between commutative algebra and some objects in combinatorics such as simplicial complexes, integral points in polytopes and graphs.


title: Joint Distributions of Permutation Statistics and the Parabolic Cylinder Functions

abstract: In this paper, we introduce a context-free grammar $G\colon x \rightarrow xy,\, y \rightarrow zu,\, z \rightarrow zw,\, w \rightarrow xv,\, u \rightarrow xyz^{-1}v,\, v \rightarrow x^{-1}zwu$ over the variable set $V=\{x,y,z,w,u,v\}$. We use this grammar to study joint distributions of several permutation statistics related to descents, rises, peaks and valleys. By considering the pattern of an exterior peak, we introduce the exterior peaks of pattern 132 and of pattern 231. Similarly, peaks can also be classified according to their patterns. Let $D$ be the formal derivative operator with respect to the grammar $G$. By using a grammatical labeling, we show that $D^n(z)$ is the generating function of the number of permutations on $[n]=\{1,2,\ldots,n\}$ with given numbers of exterior peaks of pattern 132 and of pattern 231, and proper double descents. By solving a cylinder differential equation, we obtain an explicit formula of the generating function of $D^n(z)$, which can be viewed as a unification of the results of Elizalde-Noy, Barry, Basset, Fu and Gessel. Specializations lead to the joint distributions of certain consecutive patterns in permutations, as studied by Elizalde-Noy and Kitaev. By a different labeling with respect to the same grammar $G$, we derive the joint distribution of peaks of pattern 132 and of pattern 231, double descents and double rises, with the generating function also expressed by the parabolic cylinder functions. This formula serves as a refinement of the work of Carlitz-Scoville. Furthermore, we obtain the joint distribution of exterior peaks of pattern 132 and of pattern 231 over alternating permutations.


title: The sharp Remez-type inequality for even trigonometric polynomials on the period

abstract: We prove that $$\max_{t \in [-\pi,\pi]}{|Q(t)|} \leq T_{2n}(\sec(s/4)) = \frac 12 ((\sec(s/4) + \tan(s/4))^{2n} + (\sec(s/4) - \tan(s/4))^{2n})$$ for every even trigonometric polynomial $Q$ of degree at most $n$ with complex coefficients satisfying $$m(\{t \in [-\pi,\pi]: |Q(t)| \leq 1\}) \geq 2\pi-s\,, \qquad s \in (0,2\pi)\,,$$ where $m(A)$ denotes the Lebesgue measure of a measurable set $A \subset {\Bbb R}$ and $T_{2n}$ is the Chebysev polynomial of degree $2n$ on $[-1,1]$ defined by $T_{2n}(\cos t) = \cos(2nt)$ for $t \in {\Bbb R}$. This inequality is sharp. We also prove that $$\max_{t \in [-\pi,\pi]}{|Q(t)|} \leq T_{2n}(\sec(s/2)) = \frac 12 ((\sec(s/2) + \tan(s/2))^{2n} + (\sec(s/2) - \tan(s/2))^{2n})$$ for every trigonometric polynomial $Q$ of degree at most $n$ with complex coefficients satisfying $$m(\{t \in [-\pi,\pi]: |Q(t)| \leq 1\}) \geq 2\pi-s\,, \qquad s \in (0,\pi)\,.$$


title: Morita Equivalent Blocks of Symmetric Groups

abstract: A well-known result of Scopes states that there are only finitely many Morita equivalence classes of $p$-blocks of symmetric groups with a given weight (or defect). In this note we investigate a lower bound on the number of those Morita equivalence classes.


title: Compounding of Wealth in Proof-of-Stake Cryptocurrencies

abstract: Proof-of-stake (PoS) is a promising approach for designing efficient blockchains, where block proposers are randomly chosen with probability proportional to their stake. A primary concern with PoS systems is the "rich getting richer" phenomenon, whereby wealthier nodes are more likely to get elected, and hence reap the block reward, making them even wealthier. In this paper, we introduce the notion of equitability, which quantifies how much a proposer can amplify her stake compared to her initial investment. Even with everyone following protocol (i.e., honest behavior), we show that existing methods of allocating block rewards lead to poor equitability, as does initializing systems with small stake pools and/or large rewards relative to the stake pool. We identify a \emph{geometric} reward function, which we prove is maximally equitable over all choices of reward functions under honest behavior and bound the deviation for strategic actions; the proofs involve the study of optimization problems and stochastic dominances of Polya urn processes, and are of independent mathematical interest. These results allow us to provide a systematic framework to choose the parameters of a practical incentive system for PoS cryptocurrencies.


title: Signatures of charge-order correlations in transport properties of electron-doped cuprate superconductors

abstract: The high-temperature superconductivity in copper oxides emerges under strong influence of spin correlations in doped Mott insulators. Recent discoveries of charge-order (CO) correlations in Y-based hole-doped cuprates as well as in electron-doped cuprates suggest that charge correlations should also play an important role on the electronic states of cuprates, although those correlations have been so far detected mainly by x-ray scattering measurements. Here we show signatures of CO correlations in transport properties of electron-doped cuprates as anomalous enhancement of the metal-to-insulator crossover temperature (Tmin) appears in the limited doping range near the onset of superconductivity, while it decreases exactly when superconductivity sets in. We explain this non-monotonous peak-like behavior of Tmin in terms of the evolution of the electronic states through development of CO correlations and appearance of the hole pockets in the folded Fermi surface, which impact on transport properties consecutively at different locations in the momentum space.


title: Joint Routing and Resource Allocation for Millimeter Wave Picocellular Backhaul

abstract: Picocellular architectures are essential for providing the spatial reuse required to satisfy the ever-increasing demand for mobile data. A key deployment challenge is to provide backhaul connections with sufficiently high data rate. Providing wired support (e.g., using optical fiber) to pico base stations deployed opportunistically on lampposts and rooftops is impractical, hence wireless backhaul becomes an attractive approach. A multihop mesh network comprised of directional millimeter wave links is considered here for this purpose. The backhaul design problem is formulated as one of joint routing and resource allocation, accounting for mutual interference across simultaneously active links. A computationally tractable formulation is developed by leveraging the localized nature of interference and the provable existence of a sparse optimal allocation. Numerical results are provided for millimeter (mm) wave mesh networks, which are well suited for scaling backhaul data rates due to abundance of spectrum, and the ability to form highly directional, electronically steerable beams.


title: Local Density Estimation in High Dimensions

abstract: An important question that arises in the study of high dimensional vector representations learned from data is: given a set $\mathcal{D}$ of vectors and a query $q$, estimate the number of points within a specified distance threshold of $q$. We develop two estimators, LSH Count and Multi-Probe Count that use locality sensitive hashing to preprocess the data to accurately and efficiently estimate the answers to such questions via importance sampling. A key innovation is the ability to maintain a small number of hash tables via preprocessing data structures and algorithms that sample from multiple buckets in each hash table. We give bounds on the space requirements and sample complexity of our schemes, and demonstrate their effectiveness in experiments on a standard word embedding dataset.


title: Thermally and field-driven mobility of emergent magnetic charges in square artificial spin ice

abstract: Designing and constructing model systems that embody the statistical mechanics of frustration is now possible using nanotechnology. We have arranged nanomagnets on a two-dimensional square lattice to form an artificial spin ice, and studied its fractional excitations, emergent magnetic monopoles, and how they respond to a driving field using X-ray magnetic microscopy. We observe a regime in which the monopole drift velocity is linear in field above a critical field for the onset of motion. The temperature dependence of the critical field can be described by introducing an interaction term into the Bean-Livingston model of field-assisted barrier hopping. By analogy with electrical charge drift motion, we define and measure a monopole mobility that is larger both for higher temperatures and stronger interactions between nanomagnets. The mobility in this linear regime is described by a creep model of zero-dimensional charges moving within a network of quasi-one-dimensional objects.


title: PP-DBLP: Modeling and Generating Attributed Public-Private Networks with DBLP

abstract: In many online social networks (e.g., Facebook, Google+, Twitter, and Instagram), users prefer to hide her/his partial or all relationships, which makes such private relationships not visible to public users or even friends. This leads to a new graph model called public-private networks, where each user has her/his own perspective of the network including the private connections. Recently, public-private network analysis has attracted significant research interest in the literature. A great deal of important graph computing problems (e.g., shortest paths, centrality, PageRank, and reachability tree) has been studied. However, due to the limited data sources and privacy concerns, proposed approaches are not tested on real-world datasets, but on synthetic datasets by randomly selecting vertices as private ones. Therefore, real-world datasets of public-private networks are essential and urgently needed for such algorithms in the evaluation of efficiency and effectiveness. In this paper, we generate four public-private networks from real-world DBLP records, called PPDBLP. We take published articles as public information and regard ongoing collaborations as the hidden information, which is only known by the authors. Our released datasets of PPDBLP offer the prospects for verifying various kinds of efficient public-private analysis algorithms in a fair way. In addition, motivated by widely existing attributed graphs, we propose an advanced model of attributed public-private graphs where vertices have not only private edges but also private attributes. We also discuss open problems on attributed public-private graphs. Preliminary experimental results on our generated real-world datasets verify the effectiveness and efficiency of public-private models and state-of-the-art algorithms.


title: Multipartite Dark Matter with Scalars, Fermions and signatures at LHC

abstract: Basic idea of this analysis is to achieve a two-component dark matter (DM) framework composed of a scalar and a fermion, with non-negligible DM-DM interaction contributing to thermal freeze out (hence relic density), but hiding them from direct detection bounds. We therefore augment the Standard Model (SM) with a scalar singlet ($S$) and three vectorlike fermions: two singlets ($\chi_1,\chi_2$) and a doublet ($N$). Stability of the two DM components is achieved by a discrete $\mathcal{Z}_2 \times {\mathcal{Z}^\prime}_2$ symmetry, under which the additional fields transform suitably. Fermion fields having same $\mathcal{Z}_2 \times {\mathcal{Z}^\prime}_2$ charge ($N,\chi_1$ in the model) mix after electroweak symmetry breaking (EWSB) and the lightest component becomes one of the DM candidates, while scalar singlet $S$ is the other DM component connected to visible sector by Higgs portal coupling. The heavy fermion ($\chi_2$) plays the role of mediator to connect the two DM candidates through Yukawa interaction. This opens up a large parameter space for the heavier DM component through DM-DM conversion. Hadronically quiet dilepton signature, arising from the fermion dark sector, can be observed at Large Hadron Collider (LHC) aided by the presence of a lighter scalar DM component, satisfying relic density and direct search bounds through DM-DM conversion.


title: Compact Low-Profile Wearable Antennas For Breast Cancer Detection

abstract: Many lives can be saved if tumors are detected in early stages, which can result in a bigger chance for recovery. Many patients find it irritating to get regular checkups due to the fact that the majority of the monitoring systems are complicated, not available everywhere and not mobile. Furthermore, for medical field applications, micro-strip antennas are efficient and have flexible properties that are utilized in imaging, diagnosis and treatment. It is known that breast cancer is the most common type of cancer in the world, and the earlier its been detected the better. In the early stages of breast cancer, getting rid of from the tumors is much easier and more guaranteed. Nowadays, the main method that is used in the hospitals for breast cancer detection is the Ultra-Wideband method (UWB). However, Many patients find it irritating to get regular check ups due to the fact that the majority of the monitoring systems are complicated and not mobile.


title: Direct Measurement of Quasar Outflow Wind Acceleration

abstract: We search for velocity changes (i.e., acceleration/deceleration) of narrow absorption lines (NALs) that are intrinsic to the quasars, using spectra of 6 bright quasars that have been observed more than once with 8-10m class telescopes. While variations in line strength and profile are frequently reported (especially in broader absorption lines), definitive evidence for velocity shifts has not been found with only a few exceptions. Direct velocity shift measurements are valuable constraints on the acceleration mechanisms. In this study, we determine velocity shifts by comparing the absorption profiles of NALs at two epochs separated by more than 10 years in the observed frame, using the cross-correlation function method and we estimate the uncertainties using Monte Carlo simulations. We do not detect any significant shifts but we obtain 3$\sigma$ upper limits on the acceleration of intrinsic NALs (compared to intervening NALs in same quasars) of $\sim$0.7 km s$^{-1}$ yr$^{-1}$ ($\sim$0.002 cm s$^{-2}$). We discuss possible scenarios for non-detection of NAL acceleration/deceleration and examine resulting constraints on the physical conditions in accretion disk wind.


title: Taming the War in Memory: A Resilient Mitigation Strategy Against Memory Safety Attacks in CPS

abstract: Memory-safety attacks have been one of the most critical threats against computing systems. Although a wide-range of defense techniques have been developed against these attacks, the existing mitigation strategies have several limitations. In particular, most of the existing mitigation approaches are based on aborting or restarting the victim program when a memory-safety attack is detected, thus making the system unavailable. This might not be acceptable in systems with stringent timing constraints, such as cyber-physical systems (CPS), since the system unavailability leaves the control system in an unsafe state. To address this problem, we propose CIMA -- a resilient and light-weight mitigation technique that prevents invalid memory accesses at runtime. CIMA manipulates the compiler-generated control flow graph to automatically detect and bypass unsafe memory accesses at runtime, thereby mitigating memory-safety attacks along the process. An appealing feature of CIMA is that it also ensures system availability and resilience of the CPS even under the presence of memory-safety attacks. To this end, we design our experimental setup based on a realistic Secure Water Treatment (SWaT) and Secure Urban Transportation System (SecUTS) testbeds and evaluate the effectiveness and the efficiency of our approach. The experimental results reveal that CIMA handles memory-safety attacks effectively with low overhead. Moreover, it meets the real-time constraints and physical-state resiliency of the CPS under test.


title: Biquadratic fields having a non-principal euclidean ideal class

abstract: H. W. Lenstra \cite{lenstra} introduced the notion of an Euclidean ideal class, which is a generalization of norm-Euclidean ideals in number fields. Later, families of number fields of small degree were obtained with an Euclidean ideal class (for instance, in \cite{hester1} and \cite{cathy}). In this paper, we construct certain new families of biquadratic number fields having a non-principal Euclidean ideal class and this extends the previously known families given by H. Graves \cite{hester1} and C. Hsu \cite{cathy}.


title: Analysis on the Integrability of the Rayleigh-Plesset Equation with Painlev\'e Test and Lie Symmetry Groups

abstract: The Rayleigh-Plesset equation (RPE) is a nonlinear ordinary differential equation (ODE) of second order that governs the dynamics of a spherical bubble and plays an essential role in interpreting many real-world phenomena involving the presence of bubbles in engineering and medical fields. In the present paper, we present a relatively comprehensive analysis of the analytical solutions of the RPE. The integrability of the Rayleigh-Plesset equation is investigated and discussed using the Painlev\'e test. Lie symmetry groups are employed subsequentially to obtain several exact solutions to the simplified Rayleigh-Plesset equations.


title: Sim-to-Real Transfer of Robot Learning with Variable Length Inputs

abstract: Current end-to-end deep Reinforcement Learning (RL) approaches require jointly learning perception, decision-making and low-level control from very sparse reward signals and high-dimensional inputs, with little capability of incorporating prior knowledge. This results in prohibitively long training times for use on real-world robotic tasks. Existing algorithms capable of extracting task-level representations from high-dimensional inputs, e.g. object detection, often produce outputs of varying lengths, restricting their use in RL methods due to the need for neural networks to have fixed length inputs. In this work, we propose a framework that combines deep sets encoding, which allows for variable-length abstract representations, with modular RL that utilizes these representations, decoupling high-level decision making from low-level control. We successfully demonstrate our approach on the robot manipulation task of object sorting, showing that this method can learn effective policies within mere minutes of highly simplified simulation. The learned policies can be directly deployed on a robot without further training, and generalize to variations of the task unseen during training.


title: $L_1$ Shortest Path Queries in Simple Polygons

abstract: Let $P$ be a simple polygon of $n$ vertices. We consider two-point $L_1$ shortest path queries in $P$. We build a data structure of $O(n)$ size in $O(n)$ time such that given any two query points $s$ and $t$, the length of an $L_1$ shortest path from $s$ to $t$ in $P$ can be computed in $O(\log n)$ time, or in $O(1)$ time if both $s$ and $t$ are vertices of $P$, and an actual shortest path can be output in additional linear time in the number of edges of the path. To achieve the result, we propose a mountain decomposition of simple polygons, which may be interesting in its own right. Most importantly, our approach is much simpler than the previous work on this problem.


title: Optimal Guaranteed Cost Control of Discrete-Time Linear Systems subject to Structured Uncertainties

abstract: In this paper, we propose state- and static output-feedback generalized guaranteed cost control (GCC) approaches for discrete-time linear systems subject to norm-bounded structured parametric uncertainties. This method enables the convex synthesis for a more general class of systems, where uncertainties are uncorrelated block diagonal, and no feed-through uncertainty is multiplicative with control input ones. It also provides necessary and sufficient conditions for state-feedback and sufficient conditions for static output-feedback. We also present a comparative study among the proposed controllers, standard Linear Quadratic Regulator, and Guaranteed Cost Controller found in the literature.


title: Building Context-aware Clause Representations for Situation Entity Type Classification

abstract: Capabilities to categorize a clause based on the type of situation entity (e.g., events, states and generic statements) the clause introduces to the discourse can benefit many NLP applications. Observing that the situation entity type of a clause depends on discourse functions the clause plays in a paragraph and the interpretation of discourse functions depends heavily on paragraph-wide contexts, we propose to build context-aware clause representations for predicting situation entity types of clauses. Specifically, we propose a hierarchical recurrent neural network model to read a whole paragraph at a time and jointly learn representations for all the clauses in the paragraph by extensively modeling context influences and inter-dependencies of clauses. Experimental results show that our model achieves the state-of-the-art performance for clause-level situation entity classification on the genre-rich MASC+Wiki corpus, which approaches human-level performance.


title: Materials knowledge system for nonlinear composites

abstract: In this contribution, we present a new Materials Knowledge System framework for microstructure-sensitive predictions of effective stress--strain responses in composite materials. The model is developed for composites with a wide range of combinations of strain hardening laws and topologies of the constituents. The theoretical foundation of the model is inspired by statistical continuum theories, leveraged by mean-field approximation of self-consistent models, and calibrated to data obtained from micromechanical finite element simulations. The model also relies on newly formulated data-driven linkages between micromechanical responses (phase-average strain rates and effective strength) and microstructure as well as strength contrast of the constituents. The paper describes in detail the theoretical development of the model, its implementation into an efficient computational plasticity framework, calibration of the linkages, and demonstration of the model predictions on two-phase composites with isotropic constituents exhibiting linear and power-law strain hardening laws. It is shown that the model reproduces finite element results reasonably well with significant savings of the computational cost.


title: A Quantitative Evaluation of Natural Language Question Interpretation for Question Answering Systems

abstract: Systematic benchmark evaluation plays an important role in the process of improving technologies for Question Answering (QA) systems. While currently there are a number of existing evaluation methods for natural language (NL) QA systems, most of them consider only the final answers, limiting their utility within a black box style evaluation. Herein, we propose a subdivided evaluation approach to enable finer-grained evaluation of QA systems, and present an evaluation tool which targets the NL question (NLQ) interpretation step, an initial step of a QA pipeline. The results of experiments using two public benchmark datasets suggest that we can get a deeper insight about the performance of a QA system using the proposed approach, which should provide a better guidance for improving the systems, than using black box style approaches.


title: A geometrical 1% distance to the short-period binary Cepheid V1334 Cygni

abstract: Cepheid stars play a considerable role as extragalactic distances indicators, thanks to the simple empirical relation between their pulsation period and their luminosity. They overlap with that of secondary distance indicators, such as Type Ia supernovae, whose distance scale is tied to Cepheid luminosities. However, the Period-Luminosity (P-L) relation still lacks a calibration to better than 5%. Using an original combination of interferometric astrometry with optical and ultraviolet spectroscopy, we measured the geometrical distance d = 720.35+/-7.84 pc of the 3.33 d period Cepheid V1334 Cyg with an unprecedented accuracy of +/-1 %, providing the most accurate distance for a Cepheid. Placing this star in the P-L diagram provides an independent test of existing period-luminosity relations. We show that the secondary star has a significant impact on the integrated magnitude, particularly at visible wavelengths. Binarity in future high precision calibrations of the P-L relations is not negligible, at least in the short-period regime. Subtracting the companion flux leaves V1334 Cyg in marginal agreement with existing photometric-based P-L relations, indicating either an overall calibration bias or a significant intrinsic dispersion at a few percent level. Our work also enabled us to determine the dynamical masses of both components, M1 = 4.288 +/- 0.133 Msun (Cepheid) and M2 = 4.040 +/- 0.048 Msun (companion), providing the most accurate masses for a Galactic binary Cepheid system.


title: Apollonian sets in taxicab geometry

abstract: Fix two points $p$ and $q$ in the plane and a positive number $k \neq 1$. A result credited to Apollonius of Perga states that the set of points $x$ that satisfy $d(x, p)/d(x, q) = k$ forms a circle. In this paper we study the analogous set in taxicab geometry. We find that while Apollonian sets are not taxicab circles, more complicated Apollonian sets can be characterized in terms of simpler ones.


title: A Novel Probabilistic Framework to Study the Impact of PV-battery Systems on Low-Voltage Distribution Networks

abstract: Battery storage, particularly residential battery storage coupled with rooftop PV, is emerging as an essential component of the smart grid technology mix. However, including battery storage and other flexible resources like electric vehicles and loads with thermal inertia into a probabilistic analysis based on Monte Carlo (MC) simulation is challenging, because their operational profiles are determined by computationally intensive optimization. Additionally, MC analysis requires a large pool of statistically-representative demand profiles to sample from. As a result, the analysis of the network impact of PV-battery systems has attracted little attention in the existing literature. To fill these knowledge gaps, this paper proposes a novel probabilistic framework to study the impact of PV-battery systems on low-voltage distribution networks. Specifically, the framework incorporates home energy management(HEM) operational decisions within the MC time series power flow analysis. First, using available smart meter data, we use a Bayesian nonparametric model to generate statistically-representative synthetic demand and PV profiles. Second, a policy function approximation that emulates battery scheduling decisions is used to make the simulation of optimization-based HEM feasible within the MC framework. The efficacy of our method is demonstrated on three representative low-voltage feeders, where the computation time to execute our MC framework is 5% of that when using explicit optimization methods in each MC sample. The assessment results show that uncoordinated battery scheduling has a limited beneficial impact, which is against the conjecture that batteries will serendipitously mitigate the technical problems induced by PV generation.


title: New narrow nucleon resonances $N^{\ast}(1685)$ and $N^{\ast}(1726)$ within the chiral quark-soliton model

abstract: We investigate the strong and radiative decay widths of the narrow nucleon resonances $N^*(1685)$ and $N^{\ast}(1726)$ within the framework of the SU(3) chiral quark-soliton model. All the relevant parameters are taken from those used to describe the properties of the baryon octet and decuplet in previous works. The masses of the antidecuplet nucleon and the eikosiheptaplet (27-plet) nucleon with spin 3/2 are determined respectively to be $(1690.2\pm 10.5)\, \mathrm{MeV}$ and $(1719.6\pm7.4)\,\mathrm{MeV}$. The decay width for $N^{\ast}(1685)\to \eta + N$ is found to be approximately three times larger than that for $N^{\ast}(1685)\to \pi +N$. The width of the decay $N^{\ast}\left(1726\right)3/2^+\to \eta +N$ is even about 31 times larger than that of $N^{\ast}\left(1726\right)3/2^+\to \pi + N$. The ratio of the radiative decays for $N^*(1685)$ is obtained to be $\Gamma_{nn^*(1685)}/\Gamma_{pp^*(1685)}=8.62\pm3.45$ which explains very well the neutron anomaly. In contrast, we find $\Gamma_{pp^*(1726)}/\Gamma_{nn^*(1726)}=3.72\pm0.64$, which indicates that the production of $N^*(1726)$ is more likely to be observed in the proton channel. We also examined the decay modes of these narrow nucleon resonances with the strangeness hadrons involved.


title: Percolation on Homology Generators in Codimension One

abstract: This paper introduces a new percolation model motivated from polymer materials. The mathematical model is defined over a random cubical set in the $d$-dimensional space $\mathbb{R}^d$ and focuses on generations and percolations of $(d-1)$-dimensional holes as higher dimensional topological objects. Here, the random cubical set is constructed by the union of unit faces in dimension $d-1$ which appear randomly and independently with probability $p$, and holes are formulated by the homology generators. Under this model, the upper and lower estimates of the critical probability $p_c^{\rm hole}$ of the hole percolation are shown in this paper, implying the existence of the phase transition. The uniqueness of infinite hole cluster is also proven. This result shows that, when $p > p_c^{\rm hole}$, the probability $P_p(x^*\overset{\rm hole}{\longleftrightarrow} y^*)$ that two points in the dual lattice $(\mathbb{Z}^d)^*$ belong to the same hole cluster is uniformly greater than 0.