title: What's In A Patch, I: Tensors, Differential Geometry and Statistical Shading Analysis

abstract: We develop a linear algebraic framework for the shape-from-shading problem, because tensors arise when scalar (e.g. image) and vector (e.g. surface normal) fields are differentiated multiple times. The work is in two parts. In this first part we investigate when image derivatives exhibit invariance to changing illumination by calculating the statistics of image derivatives under general distributions on the light source. We computationally validate the hypothesis that image orientations (derivatives) provide increased invariance to illumination by showing (for a Lambertian model) that a shape-from-shading algorithm matching gradients instead of intensities provides more accurate reconstructions when illumination is incorrectly estimated under a flatness prior.


title: Mn-induced magnetic symmetry breaking and its correlation with the metal-insulator transition in bilayered Sr3(Ru1-xMnx)2O7

abstract: Bilayered Sr3Ru2O7 is an unusual metamagnetic metal with inherently antiferromagnetic (AFM) and ferromagnetic (FM) fluctuations. Partial substitution of Ru by Mn results in the establishment of metal-insulator transition (MIT) at TMIT and AFM ordering at TM in Sr3(Ru1-xMnx)2O7. Using elastic neutron scattering we determined the effect of Mn doping on the magnetic structure and in-plane magnetic correlation lengths in Sr3(Ru1-xMnx)2O7 (x = 0.06 and 0.12). With increasing Mn doping (x) from 0.06 to 0.12 or decreasing temperatures for x=0.12, an evolution from an in-plane short-range to long-range double-stripe AFM ground state occurs. For both compounds, the onset of magnetic correlation with an anisotropic behavior coincides with the sharp rise of the electrical resistivity and the specific heat. Since it does not induce measurable lattice distortion, the double-stripe magnetic order with anisotropic spin texture breaks the symmetry from C4v crystal lattice to C2v magnetic sublattice. These observations shed new light on an age-old question of Slater versus Mott-type MIT.


title: Coupling MOS Quantum Dot and Phosphorus Donor Qubit Systems

abstract: Si-MOS based QD qubits are attractive due to their similarity to the current semiconductor industry. We introduce a highly tunable MOS foundry compatible qubit design that couples an electrostatic quantum dot (QD) with an implanted donor. We show for the first time coherent two-axis control of a two-electron spin qubit that evolves under the QD-donor exchange interaction and the hyperfine interaction with the donor nucleus. The two interactions are tuned electrically with surface gate voltages to provide control of both qubit axes. Qubit decoherence is influenced by charge noise, which is of similar strength as epitaxial systems like GaAs and Si/SiGe.


title: The population of white dwarf-main sequence binaries in the SDSS DR 12

abstract: We present a Monte Carlo population synthesis study of white dwarf-main sequence (WD+MS) binaries in the Galactic disk aimed at reproducing the ensemble properties of the entire population observed by the Sloan Digital Sky Survey (SDSS) Data Release 12. Our simulations take into account all known observational biases and use the most up-to-date stellar evolutionary models. This allows us to perform a sound comparison between the simulations and the observational data. We find that the properties of the simulated and observed parameter distributions agree best when assuming low values of the common envelope efficiency (0.2-0.3), a result that is in agreement with previous findings obtained by observational and population synthesis studies of close SDSS WD+MS binaries. We also show that all synthetic populations that result from adopting an initial mass ratio distribution with a positive slope are excluded by observations. Finally, we confirm that the properties of the simulated WD+MS binary populations are nearly independent of the age adopted for the thin disk, on the contribution of WD+MS binaries from the thick disk (0-17 per cent of the total population) and on the assumed fraction of the internal energy that is used to eject the envelope during the common envelope phase (0.1-0.5).


title: First simultaneous extraction of spin-dependent parton distributions and fragmentation functions from a global QCD analysis

abstract: We perform the first global QCD analysis of polarized inclusive and semi-inclusive deep-inelastic scattering and single-inclusive $e^+e^-$ annihilation data, fitting simultaneously the parton distribution and fragmentation functions using the iterative Monte Carlo method. Without imposing SU(3) symmetry relations, we find the strange polarization to be very small, consistent with zero for both inclusive and semi-inclusive data, which provides a resolution to the strange quark polarization puzzle. The combined analysis also allows the direct extraction from data of the isovector and octet axial charges, and is consistent with a small SU(2) flavor asymmetry of the polarized sea.


title: The nuts and bolts of the BMS Bootstrap

abstract: In this paper, we elaborate on aspects of the recently introduced BMS bootstrap programme. We consider two-dimensional (2d) field theories with BMS3 symmetry and extensively use highest weight representations to uncover the BMS version of crossing symmetry in 4-point functions that are constrained by symmetry. The BMS bootstrap equation is formulated and then analytic expressions for BMS blocks are constructed by looking at the limit of large central charges. These results are also applicable to 2d Galilean Conformal Field Theories through the isomorphism between the BMS3 and 2d Galilean Conformal Algebras. We recover our previously obtained results in the non-relativistic limit of the corresponding ones in 2d relativistic CFTs. This provides a comprehensive check of our previous analysis. We also explore the chiral limit of BMS3 where the BMS algebra reduces to a single copy of the Virasoro algebra and show that our analysis is consistent with earlier work in this direction.


title: $D \to \rho \,\ell^+\ell^-$ Decays in the QCD Factorization Approach

abstract: We consider rare semileptonic decays of a heavy $D$-meson into a light vector meson in the framework of QCD factorization. In contrast to the corresponding $B$-meson decays, the naive factorization hypothesis does not even serve as a first approximation. Rather, the decay amplitudes appear to be dominated by non-factorizable dynamics, e.g. through annihilation topologies, which are particularly sensitive to long-distance hadronic contributions. We therefore pay particular attention to intermediate vector-meson resonances appearing in quark-loop and annihilation topologies. Compared to the analogous $B$-meson decays, we identify a number of effects that result in very large theoretical uncertainties for differential decay rates. Some of these effects are found to cancel in the ratio of partially integrated decay rates for transversely and longitudinally polarized $\rho$ mesons. On the phenomenological side this implies a very limited potential to constrain physics beyond the Standard Model by means of these decays.


title: Warm-hot Gas in X-ray Bright Galaxy Clusters and the H I-deficient Circumgalactic Medium in Dense Environments

abstract: We analyze the intracluster medium (ICM) and circumgalactic medium (CGM) in 7 X-ray detected galaxy clusters using spectra of background QSOs (HST-COS/STIS), optical spectroscopy of the cluster galaxies (MMT/Hectospec and SDSS), and X-ray imaging/spectroscopy (XMM-Newton and Chandra). First, we report a very low covering fraction of H I absorption in the CGM of these cluster galaxies, f_c = 25% +25%/-15%, to stringent detection limits (log N(H I) < 13 cm^-2). As field galaxies have an H I covering fraction of ~100% at similar radii, the dearth of CGM H I in our data indicates that the cluster environment has effectively stripped or overionized the gaseous halos of these cluster galaxies. Second, we assess the contribution of warm-hot (10^5 - 10^6 K) gas to the ICM as traced by O VI and broad Ly-alpha (BLA) absorption. Despite the high signal-to-noise of our data, we do not detect O VI in any cluster, and we only detect BLA features in the QSO spectrum probing one cluster. We estimate that the total column density of warm-hot gas along this line of sight totals to ~3% of that contained in the hot T > 10^7 K X-ray emitting phase. Residing at high relative velocities, these features may trace pre-shocked material outside the cluster. Comparing gaseous galaxy halos from the low-density 'field' to galaxy groups and high-density clusters, we find that the CGM is progressively depleted of H I with increasing environmental density, and the CGM is most severely transformed in galaxy clusters. This CGM transformation may play a key role in environmental galaxy quenching.


title: Computed Axial Lithography (CAL): Toward Single Step 3D Printing of Arbitrary Geometries

abstract: Most additive manufacturing processes today operate by printing voxels (3D pixels) serially point-by-point to build up a 3D part. In some more recently-developed techniques, for example optical printing methods such as projection stereolithography [Zheng et al. 2012], [Tumbleston et al. 2015], parts are printed layer-by-layer by curing full 2d (very thin in one dimension) layers of the 3d part in each print step. There does not yet exist a technique which is able to print arbitrarily-defined 3D geometries in a single print step. If such a technique existed, it could be used to expand the range of printable geometries in additive manufacturing and relax constraints on factors such as overhangs in topology optimization. It could also vastly increase print speed for 3D parts. In this work, we develop the principles for an approach for single exposure 3D printing of arbitrarily defined geometries. The approach, termed Computed Axial Lithgography (CAL), is based on tomographic reconstruction, with mathematical optimization to generate a set of projections to optically define an arbitrary dose distribution within a target volume. We demonstrate the potential ability of the technique to print 3D parts using a prototype CAL system based on sequential illumination from many angles. We also propose new hardware designs which will help us to realize true single-shot arbitrary-geometry 3D CAL.


title: Probing spinon nodal structures in three-dimensional Kitaev spin liquids

abstract: We propose that resonant inelastic X-ray scattering (RIXS) is an effective probe of the fractionalized excitations in three-dimensional (3D) Kitaev spin liquids. While the non-spin-conserving RIXS responses are dominated by the gauge-flux excitations and reproduce the inelastic-neutron-scattering response, the spin-conserving (SC) RIXS response picks up the Majorana-fermion excitations and detects whether they are gapless at Weyl points, nodal lines, or Fermi surfaces. As a signature of symmetry fractionalization, the SC RIXS response is suppressed around the $\Gamma$ point. On a technical level, we calculate the exact SC RIXS responses of the Kitaev models on the hyperhoneycomb, stripyhoneycomb, hyperhexagon, and hyperoctagon lattices, arguing that our main results also apply to generic 3D Kitaev spin liquids beyond these exactly solvable models.


title: Highly connected 7-manifolds and non-negative curvature

abstract: In this article, a six-parameter family of highly connected 7-manifolds which admit an SO(3)-invariant metric of non-negative sectional curvature is constructed and the Eells-Kuiper invariant of each is computed. In particular, it follows that all exotic spheres in dimension 7 admit an SO(3)-invariant metric of non-negative curvature.


title: Completeness Theorems for Pomset Languages and Concurrent Kleene Algebras

abstract: Pomsets constitute one of the most basic models of concurrency. A pomset is a generalisation of a word over an alphabet in that letters may be partially ordered. A term $t$ using the bi-Kleene operations $0,1, +, \cdot\, ,^*, \parallel, ^{(*)}$ defines a language $ \mathopen{[\![ } t \mathclose{]\!] } $ of pomsets in a natural way. We prove that every valid universal equality over pomset languages using these operations is a consequence of the equational theory of regular languages (in which parallel multiplication and iteration are undefined) plus that of the commutative-regular languages (in which sequential multiplication and iteration are undefined). We also show that the class of $\textit{rational}$ pomset languages (that is, those languages generated from singleton pomsets using the bi-Kleene operations) is closed under all Boolean operations. An $ \textit{ideal}$ of a pomset $p$ is a pomset using the letters of $p$, but having an ordering at least as strict as $p$. A bi-Kleene term $t$ thus defines the set $ \textbf{Id} (\mathopen{[\![ } t \mathclose{]\!] }) $ of ideals of pomsets in $ \mathopen{[\![ } t \mathclose{]\!] } $. We prove that if $t$ does not contain commutative iteration $^{(*)}$ (in our terminology, $t$ is bw-rational) then $\textbf{Id} (\mathopen{[\![ } t \mathclose{]\!] }) \cap \textbf{Pom}_{sp}$, where $ \textbf{Pom}_{sp}$ is the set of pomsets generated from singleton pomsets using sequential and parallel multiplication ($ \cdot$ and $ \parallel$) is defined by a bw-rational term, and if two such terms $t,t'$ define the same ideal language, then $t'=t$ is provable from the Kleene axioms for $0,1, +, \cdot\, ,^*$ plus the commutative idempotent semiring axioms for $0,1, +, \parallel$ plus the exchange law $ (u \parallel v)\cdot ( x \parallel y) \le v \cdot y \parallel u \cdot x $.


title: Automated Cryptographic Analysis of the Pedersen Commitment Scheme

abstract: Aiming for strong security assurance, recently there has been an increasing interest in formal verification of cryptographic constructions. This paper presents a mechanised formal verification of the popular Pedersen commitment protocol, proving its security properties of correctness, perfect hiding, and computational binding. To formally verify the protocol, we extended the theory of EasyCrypt, a framework which allows for reasoning in the computational model, to support the discrete logarithm and an abstraction of commitment protocols. Commitments are building blocks of many cryptographic constructions, for example, verifiable secret sharing, zero-knowledge proofs, and e-voting. Our work paves the way for the verification of those more complex constructions.


title: Analysis of the non-Markovianity for electron transfer reactions in an oligothiophene-fullerene heterojunction

abstract: The non-Markovianity of the electron transfer in an oligothiophene-fullerene heterojunction described by a spin-boson model is analyzed using the time dependent decoherence canonical rates and the volume of accessible states in the Bloch sphere. The dynamical map of the reduced electronic system is computed by the hierarchical equations of motion methodology (HEOM) providing an exact dynamics. Transitory witness of non-Markovianity is linked to the bath dynamics analyzed from the HEOM auxiliary matrices. The signature of the collective bath mode detected from HEOM in each electronic state is compared with predictions of the effective mode extracted from the spectral density. We show that including this main reaction coordinate in a one-dimensional vibronic system coupled to a residual bath satisfactorily describes the electron transfer by a simple Markovian Redfield equation. Non-Markovianity is computed for three inter fragment distances and compared with a priori criterion based on the system and bath characteristic timescales.


title: Scalability analysis of large-scale LoRaWAN networks in ns-3

abstract: As LoRaWAN networks are actively being deployed in the field, it is important to comprehend the limitations of this Low Power Wide Area Network technology. Previous work has raised questions in terms of the scalability and capacity of LoRaWAN networks as the number of end devices grows to hundreds or thousands per gateway. Some works have modeled LoRaWAN networks as pure ALOHA networks, which fails to capture important characteristics such as the capture effect and the effects of interference. Other works provide a more comprehensive model by relying on empirical and stochastic techniques. This work uses a different approach where a LoRa error model is constructed from extensive complex baseband bit error rate simulations and used as an interference model. The error model is combined with the LoRaWAN MAC protocol in an ns-3 module that enables to study multi channel, multi spreading factor, multi gateway, bi-directional LoRaWAN networks with thousands of end devices. Using the lorawan ns-3 module, a scalability analysis of LoRaWAN shows the detrimental impact of downstream traffic on the delivery ratio of confirmed upstream traffic. The analysis shows that increasing gateway density can ameliorate but not eliminate this effect, as stringent duty cycle requirements for gateways continue to limit downstream opportunities.


title: Lie algebras of vector fields on smooth affine varieties

abstract: We reprove the results of Jordan [18] and Siebert [31] and show that the Lie algebra of polynomial vector fields on an irreducible affine variety X is simple if and only if X is a smooth variety. Given proof is self-contained and does not depend on papers mentioned above. Besides, the structure of the module of polynomial functions on an irreducible smooth affine variety over the Lie algebra of vector fields is studied. Examples of Lie algebras of polynomial vector fields on an N-dimensional sphere, non-singular hyperelliptic curves and linear algebraic groups are considered.


title: The Alexander Polynomial of a Rational Link

abstract: We relate some terms on the boundary of the Newton polygon of the Alexander polynomial $\Delta(x,y)$ of a rational link to the number and length of monochromatic twist sites in a particular diagram that we call the standard form. Normalize $\Delta(x,y)$ so that no $x^{-1}$ or $y^{-1}$ terms appear, but $x^{-1}\Delta(x,y)$ and $y^{-1}\Delta(x,y)$ have negative exponents, and so that terms of even total degree are positive and terms with odd total degree are negative. If the rational link has a reduced alternating diagram with no self crossings, then $\Delta(-1, 0) = 1$. If the standard form of the rational link has $m$ monochromatic twist sites, and the $j^{\textrm{th}}$ monochromatic twist site has $\hat{q}_j$ crossings, then $\Delta(-1, 0) = \prod_{j=1}^{m}(\hat{q}_j+1)$. Our proof employs Kauffman's clock moves and a lattice for the terms of $\Delta(x,y)$ in which the $y$-power cannot decrease.


title: Tensors, Differential Geometry and Statistical Shading Analysis

abstract: We develop a linear algebraic framework for the shape-from-shading problem, because tensors arise when scalar (e.g. image) and vector (e.g. surface normal) fields are differentiated multiple times. Using this framework, we first investigate when image derivatives exhibit invariance to changing illumination by calculating the statistics of image derivatives under general distributions on the light source. Second, we apply that framework to develop Taylor-like expansions, and build a boot-strapping algorithm to find the polynomial surface solutions (under any light source) consistent with a given patch to arbitrary order. A generic constraint on the light source restricts these solutions to a 2-D subspace, plus an unknown rotation matrix. It is this unknown matrix that encapsulates the ambiguity in the problem. Finally, we use the framework to computationally validate the hypothesis that image orientations (derivatives) provide increased invariance to illumination by showing (for a Lambertian model) that a shape-from-shading algorithm matching gradients instead of intensities provides more accurate reconstructions when illumination is incorrectly estimated under a flatness prior.


title: A Tube Dynamics Perspective Governing Stability Transitions: An Example Based on Snap-through Buckling

abstract: The equilibrium configuration of an engineering structure, able to withstand a certain loading condition, is usually associated with a local minimum of the underlying potential energy. However, in the nonlinear context, there may be other equilibria present, and this brings with it the possibility of a transition to an alternative (remote) minimum. That is, given a sufficient disturbance, the structure might buckle, perhaps suddenly, to another shape. This paper considers the dynamic mechanisms under which such transitions (typically via saddle points) occur. A two-mode Hamiltonian is developed for a shallow arch/buckled beam. The resulting form of the potential energy---two stable wells connected by rank-1 saddle points---shows an analogy with resonance transitions in celestial mechanics or molecular reconfigurations in chemistry, whereas here the transition corresponds to switching between two stable structural configurations. Then, from Hamilton's equations, the analytical equilibria are determined and linearization of the equations of motion about the saddle is obtained. After computing the eigenvalues and eigenvectors of the coefficient matrix associated with the linearization, a symplectic transformation is given which puts the Hamiltonian into normal form and simplifies the equations, allowing us to use the conceptual framework known as tube dynamics. The flow in the equilibrium region of phase space as well as the invariant manifold tubes in position space are discussed. Also, we account for the addition of damping in the tube dynamics framework, which leads to a richer set of behaviors in transition dynamics than previously explored.


title: Motion-Compensated Autonomous Scanning for Tumour Localisation using Intraoperative Ultrasound

abstract: Intraoperative ultrasound facilitates localisation of tumour boundaries during minimally invasive procedures. Autonomous ultrasound scanning systems have been recently proposed to improve scanning accuracy and reduce surgeons' cognitive load. However, current methods mainly consider static scanning environments typically with the probe pressing against the tissue surface. In this work, a motion-compensated autonomous ultrasound scanning system using the da Vinci Research Kit (dVRK) is proposed. An optimal scanning trajectory is generated considering both the tissue surface shape and the ultrasound transducer dimensions. A robust vision-based approach is proposed to learn the underlying tissue motion characteristics. The learned motion model is then incorporated into the visual servoing framework. The proposed system has been validated with both phantom and ex vivo experiments using the ground truth motion data for comparison.


title: Derivation of the Navier - Stokes - Poisson system with radiation for an accretion disk

abstract: We study the 3-D compressible barotropic radiation fluid dynamics system describing the motion of the compressible rotating viscous fluid with gravitation and radiation confined to a straight layer. We show that weak solutions in the 3-D domain converge to the strong solution of the rotating 2-D Navier-Stokes-Poisson system with radiation for all times less than the maximal life time of the strong solution of the 2-D system when the Froude number is small or to the strong solution of the rotating pure 2-D Navier- Stokes system with radiation.


title: On wormholes with long throats and the stability problem

abstract: We construct explicit examples of globally regular static, spherically symmetric solutions of general relativity with a phantom scalar field as the source of gravity, describing traversable wormholes with flat asymptotic regions on both sides of the throat as well as regular black holes, in particular, those called black universes. To explain why such phantom fields are not observed under usual conditions, we invoke the concept of "invisible ghosts," which means that the phantom field decays quickly enough at infinity and is there too weak to be observed. This approach leads to wormhole models in which the spherical radius is almost constant in some range of the radial coordinate near the throat, forming a "long throat". We discuss the peculiar features and difficulties of the stability problem for such configurations. It is shown that the limiting case of a "maximally long throat" has the form of an unstable model with the Nariai metric. This allows us to conjecture that a long throat does not stabilize wormholes with a scalar source.


title: Machine Learning Molecular Dynamics for the Simulation of Infrared Spectra

abstract: Machine learning has emerged as an invaluable tool in many research areas. In the present work, we harness this power to predict highly accurate molecular infrared spectra with unprecedented computational efficiency. To account for vibrational anharmonic and dynamical effects -- typically neglected by conventional quantum chemistry approaches -- we base our machine learning strategy on ab initio molecular dynamics simulations. While these simulations are usually extremely time consuming even for small molecules, we overcome these limitations by leveraging the power of a variety of machine learning techniques, not only accelerating simulations by several orders of magnitude, but also greatly extending the size of systems that can be treated. To this end, we develop a molecular dipole moment model based on environment dependent neural network charges and combine it with the neural network potentials of Behler and Parrinello. Contrary to the prevalent big data philosophy, we are able to obtain very accurate machine learning models for the prediction of infrared spectra based on only a few hundreds of electronic structure reference points. This is made possible through the introduction of a fully automated sampling scheme and the use of molecular forces during neural network potential training. We demonstrate the power of our machine learning approach by applying it to model the infrared spectra of a methanol molecule, n-alkanes containing up to 200 atoms and the protonated alanine tripeptide, which at the same time represents the first application of machine learning techniques to simulate the dynamics of a peptide. In all these case studies we find excellent agreement between the infrared spectra predicted via machine learning models and the respective theoretical and experimental spectra.


title: The Effect of Collective Attention on Controversial Debates on Social Media

abstract: We study the evolution of long-lived controversial debates as manifested on Twitter from 2011 to 2016. Specifically, we explore how the structure of interactions and content of discussion varies with the level of collective attention, as evidenced by the number of users discussing a topic. Spikes in the volume of users typically correspond to external events that increase the public attention on the topic -- as, for instance, discussions about `gun control' often erupt after a mass shooting. This work is the first to study the dynamic evolution of polarized online debates at such scale. By employing a wide array of network and content analysis measures, we find consistent evidence that increased collective attention is associated with increased network polarization and network concentration within each side of the debate; and overall more uniform lexicon usage across all users.


title: Castelnuovo-Mumford regularity of Koszul cycles and Koszul homologies

abstract: We extend to one dimensional quotients the result of A. Conca and S. Murai on the convexity of the regularity of Koszul cycles. By providing a relation between the regularity of Koszul cycles and Koszul homologies we prove a sharp regularity bound for the Koszul homologies of a homogeneous ideal in a polynomial ring under the same conditions.


title: Graphene channels interfaced with an array of individual quantum dots

abstract: Surface graphene guides were interfaced with an array of individual semiconductor quantum dots, whose position was commensurate with the optical guide modes. The surface guide served as a channel for a Field Effect Transistor (FET) while the dots were placed within the capacitor formed between the graphene channel and the gate electrode. We report on negative differential photo-related conductance under light and a diminishing fluorescence effect as a function of bias. We suggest that the quenched fluorescence may be hindered, to some degree, by incorporating the QD in a resonator, which is tuned to the emission wavelength.


title: Four NP-complete problems about generalizations of perfect graphs

abstract: We show that the following problems are NP-complete. 1. Can the vertex set of a graph be partitioned into two sets such that each set induces a perfect graph? 2. Is the difference between the chromatic number and clique number at most $1$ for every induced subgraph of a graph? 3. Can the vertex set of every induced subgraph of a graph be partitioned into two sets such that the first set induces a perfect graph, and the clique number of the graph induced by the second set is smaller than that of the original induced subgraph? 4. Does a graph contain a stable set whose deletion results in a perfect graph? The proofs of the NP-completeness of the four problems follow the same pattern: Showing that all the four problems are NP-complete when restricted to triangle-free graphs by using results of Maffray and Preissmann on $3$-colorability and $4$-colorability of triangle-free graphs


title: Coupling conditions for globally stable and robust synchrony of chaotic systems

abstract: We propose a set of general coupling conditions to select a coupling profile (a set of coupling matrices) from the linear flow matrix (LFM) of dynamical systems for realizing global stability of complete synchronization (CS) in identical systems and robustness to parameter perturbation. The coupling matrices define the coupling links between any two oscillators in a network that consists of a conventional diffusive coupling link (self-coupling link) as well as a cross-coupling link. The addition of a selective cross-coupling link in particular plays constructive roles that ensure the global stability of synchrony and furthermore enables robustness of synchrony against small to non-small parameter perturbation. We elaborate the general conditions for the selection of coupling profiles for two coupled systems, three- and four-node network motifs analytically as well as numerically using benchmark models, the Lorenz system, the Hindmarsh-Rose neuron model, the Shimizu-Morioka laser model, the R\"ossler system and a Sprott system. The role of the cross-coupling link is, particularly, exemplified with an example of a larger network where it saves the network from a breakdown of synchrony against large parameter perturbation in any node. The perturbed node in the network transits from CS to generalized synchronization (GS) when all the other nodes remain in CS. The GS is manifested by an amplified response of the perturbed node in a coherent state.


title: Unveiling the role of competing fluctuations at an unconventional quantum critical point

abstract: Quantum critical points (QCPs) are widely accepted as a source of a diverse set of collective quantum phases of matter. A central question is how the order parameters of phases near a QCP interact and determine the fundamental character of the critical dynamics which drive the quantum critical behavior. One of the most interesting proposals for the quantum critical behavior that occurs in correlated electron systems is that the behavior may arise from local, as opposed to long wavelength, critical fluctuations of the order parameter. The local criticality is believed to give rise to energy over temperature ($E/T$) scaling of the dynamic susceptibility with a fractional exponent near the quantum critical point (QCP). Here we show that $E/T$ scaling is indeed observed for CeCu$_{6-x}$Ag$_x$ but on closer inspection, the fluctuations can be separated into two components, implying that multiple order parameters play an important role in the unconventional critical behavior. Additionally, when the fluctuations corresponding to the magnetically ordered side of the phase diagram are separated, they are found to be three dimensional and to obey the scaling behavior expected for long wavelength fluctuations near an itinerant antiferromagnetic QCP.


title: The oscillatory behavior of light in the composite Goos-Haenchen shift

abstract: For incidence in the critical region, the propagation of gaussian lasers through triangular dielectric blocks is characterized by the joint action of angular deviations and lateral displacements. This mixed effect, known as composite Goos-Haenchen shift, produces a lateral displacement dependent on the axial coordinate, recently confirmed by a weak measurement experiment. We discuss under which conditions this axial lateral displacement, which only exists for the composite Goos-Haenchen shift, presents an oscillatory behavior. This oscillation phenomenon shows a peculiar behavior of light for critical incidence and, if experimentally tested, could stimulate further theoretical studies and lead to interesting optical applications.


title: Lifted Polymatroid Inequalities for Mean-Risk Optimization with Indicator Variables

abstract: We investigate a mixed 0-1 conic quadratic optimization problem with indicator variables arising in mean-risk optimization. The indicator variables are often used to model non-convexities such as fixed charges or cardinality constraints. Observing that the problem reduces to a submodular function minimization for its binary restriction, we derive three classes of strong convex valid inequalities by lifting the polymatroid inequalities on the binary variables. Computational experiments demonstrate the effectiveness of the inequalities in strengthening the convex relaxations and, thereby, improving the solution times for mean-risk problems with fixed charges and cardinality constraints significantly.


title: Network Design with Probabilistic Capacities

abstract: We consider a network design problem with random arc capacities and give a formulation with a probabilistic capacity constraint on each cut of the network. To handle the exponentially-many probabilistic constraints a separation procedure that solves a nonlinear minimum cut problem is introduced. For the case with independent arc capacities, we exploit the supermodularity of the set function defining the constraints and generate cutting planes based on the supermodular covering knapsack polytope. For the general correlated case, we give a reformulation of the constraints that allows to uncover and utilize the submodularity of a related function. The computational results indicate that exploiting the underlying submodularity and supermodularity arising with the probabilistic constraints provides significant advantages over the classical approaches.


title: A discussion about LNG Experiment: Irreversible or Reversible Generation of the OR Logic Gate?

abstract: In a recent paper M. Lopez-Suarez, I. Neri, and L. Gammaitoni (LNG) present a concrete realization of the Boolean OR irreversible gate, but contrary to the standard Landauer principle, with an arbitrary small dissipation of energy. A Popperian good falsification! In this paper we discuss a theoretical description of the LNG device which is indeed a 3in/3out self--reversible realization of the involved OR gate, satisfying in this way the Landauer principle of no dispersion of energy, contrary to the LNG conclusions. The different point of view is due to a different interpretation of the two outputs corresponding to the inputs 10 and 01, which are considered by LNG indistinguishable so producing a non reversible realization of the standard 2in/1out gate. On the contrary, always considering these two outputs indistinguishable, by a suitable normalization function of the cantilever angles, the experimental results obtained by the LNG device coincide with the OR connective obtained from the third output of the self-reversible 3in/3out CL gate by the Inputs-Ancilla->Garbage-Output procedure. Thus, by the self-reversibility this realization is without dissipation of energy according to the Landauer principle. Furthermore, using the self-reversible Toffoli gate it is possible to obtain from the LNG device the realization of the connective AND adopting another normalization function on the cantilever angles. Finally, by other suitable normalization procedures on cantilever angles it is possible to obtain also the other logic NOR and NAND connectives, and in a more sophisticated way the XOR and NXOR connectives in a self-reversible way. All this leads to introduce a universal logic machine consisting of the LNG device plus a memory containing all the necessary angle normalization functions to produce in a self-reversible way, by choosing one of these latter, the logic connectives now listed.


title: Submodularity in conic quadratic mixed 0-1 optimization

abstract: We describe strong convex valid inequalities for conic quadratic mixed 0-1 optimization. These inequalities can be utilized for solving numerous practical nonlinear discrete optimization problems from value-at-risk minimization to queueing system design, from robust interdiction to assortment optimization through appropriate conic quadratic mixed 0-1 relaxations. The inequalities exploit the submodularity of the binary restrictions and are based on the polymatroid inequalities over binaries for the diagonal case. We prove that the convex inequalities completely describe the convex hull of a single conic quadratic constraint as well as the rotated cone constraint over binary variables and unbounded continuous variables. We then generalize and strengthen the inequalities by incorporating additional constraints of the optimization problem. Computational experiments on mean-risk optimization with correlations, assortment optimization, and robust conic quadratic optimization indicate that the new inequalities strengthen the convex relaxations substantially and lead to significant performance improvements.


title: DeepGO: Predicting protein functions from sequence and interactions using a deep ontology-aware classifier

abstract: A large number of protein sequences are becoming available through the application of novel high-throughput sequencing technologies. Experimental functional characterization of these proteins is time-consuming and expensive, and is often only done rigorously for few selected model organisms. Computational function prediction approaches have been suggested to fill this gap. The functions of proteins are classified using the Gene Ontology (GO), which contains over 40,000 classes. Additionally, proteins have multiple functions, making function prediction a large-scale, multi-class, multi-label problem. We have developed a novel method to predict protein function from sequence. We use deep learning to learn features from protein sequences as well as a cross-species protein-protein interaction network. Our approach specifically outputs information in the structure of the GO and utilizes the dependencies between GO classes as background information to construct a deep learning model. We evaluate our method using the standards established by the Computational Assessment of Function Annotation (CAFA) and demonstrate a significant improvement over baseline methods such as BLAST, with significant improvement for predicting cellular locations.


title: Path Cover and Path Pack Inequalities for the Capacitated Fixed-Charge Network Flow Problem

abstract: Capacitated fixed-charge network flows are used to model a variety of problems in telecommunication, facility location, production planning and supply chain management. In this paper, we investigate capacitated path substructures and derive strong and easy-to-compute \emph{path cover and path pack inequalities}. These inequalities are based on an explicit characterization of the submodular inequalities through a fast computation of parametric minimum cuts on a path, and they generalize the well-known flow cover and flow pack inequalities for the single-node relaxations of fixed-charge flow models. We provide necessary and sufficient facet conditions. Computational results demonstrate the effectiveness of the inequalities when used as cuts in a branch-and-cut algorithm.


title: Genetic load makes cancer cells more sensitive to common drugs: evidence from Cancer Cell Line Encyclopedia

abstract: Genetic alterations initiate tumors and enable the evolution of drug resistance. The pro-cancer view of mutations is however incomplete, and several studies show that mutational load can reduce tumor fitness. Given its negative effect, genetic load should make tumors more sensitive to anticancer drugs. Here, we test this hypothesis across all major types of cancer from the Cancer Cell Line Encyclopedia, which provides genetic and expression data of 496 cell lines together with their response to 24 common anticancer drugs. We found that the efficacy of 9 out of 24 drugs showed significant association with genetic load in a pan-cancer analysis. The associations for some tissue-drug combinations were remarkably strong, with genetic load explaining up to 83% of the variance in the drug response. Overall, the role of genetic load depended on both the drug and the tissue type with 10 tissues being particularly vulnerable to genetic load. We also identified changes in gene expression associated with increased genetic load, which included cell-cycle checkpoints, DNA damage and apoptosis. Our results show that genetic load is an important component of tumor fitness and can predict drug sensitivity. Beyond being a biomarker, genetic load might be a new, unexplored vulnerability of cancer.


title: LCDet: Low-Complexity Fully-Convolutional Neural Networks for Object Detection in Embedded Systems

abstract: Deep convolutional Neural Networks (CNN) are the state-of-the-art performers for object detection task. It is well known that object detection requires more computation and memory than image classification. Thus the consolidation of a CNN-based object detection for an embedded system is more challenging. In this work, we propose LCDet, a fully-convolutional neural network for generic object detection that aims to work in embedded systems. We design and develop an end-to-end TensorFlow(TF)-based model. Additionally, we employ 8-bit quantization on the learned weights. We use face detection as a use case. Our TF-Slim based network can predict different faces of different shapes and sizes in a single forward pass. Our experimental results show that the proposed method achieves comparative accuracy comparing with state-of-the-art CNN-based face detection methods, while reducing the model size by 3x and memory-BW by ~4x comparing with one of the best real-time CNN-based object detector such as YOLO. TF 8-bit quantized model provides additional 4x memory reduction while keeping the accuracy as good as the floating point model. The proposed model thus becomes amenable for embedded implementations.


title: Convolutions for Berezin quantization and Berezin-Lieb inequalities

abstract: Notions and results from quantum harmonic analysis, such as the convolution between functions and operators or between two operators, is identified as the appropriate setting for Berezin quantization and Berezin-Lieb inequalities. Based on this insight we provide a rigorous approach to generalized phase-space representation introduced by Klauder-Skagerstam and their variants of Berezin-Lieb inequalities in this setting. Hence our presentation of the results of Klauder-Skagerstam gives a more conceptual framework, which yields as a byproduct an interesting perspective on the connection between Berezin quantization and Weyl quantization.


title: A Comprehensive Introduction to the Theory of Word-Representable Graphs

abstract: Letters $x$ and $y$ alternate in a word $w$ if after deleting in $w$ all letters but the copies of $x$ and $y$ we either obtain a word $xyxy\cdots$ (of even or odd length) or a word $yxyx\cdots$ (of even or odd length). A graph $G=(V,E)$ is word-representable if and only if there exists a word $w$ over the alphabet $V$ such that letters $x$ and $y$ alternate in $w$ if and only if $xy\in E$. Word-representable graphs generalize several important classes of graphs such as circle graphs, $3$-colorable graphs and comparability graphs. This paper offers a comprehensive introduction to the theory of word-representable graphs including the most recent developments in the area.


title: The Large N Limit with Vanishing Leading Order Condensate for Zero Pion Mass

abstract: It is conventionally assumed that the negative mass squared term in the linear sigma model version of the pion Lagrangian is $M^2 \sim \Lambda_{\rm QCD}^2$ in powers of $N_c$. We consider the case where $M^2 \sim \Lambda^2_{\rm QCD}/N_c$ so that to leading order in $N_c$ this symmetry breaking term vanishes. We present some arguments why this might be plausible. One might think that such a radical assumption would contradict lattice Monte Carlo data on QCD as function of $N_c$. We show that the linear sigma model gives a fair description of the data of DeGrand and Liu both for $N_c = 3$, and for variable $N_c$. The values of quark masses considered by DeGrand and Liu, and by Bali et. al. turn out to be too large to resolve the case we consider from that of the conventional large $N_c$ limit. We argue that for quark masses $m_{q} \ll \Lambda_{\rm QCD}/N_c^{3/2}$, both the baryon mass and nucleon size scale as $\sqrt{N_c}$. For $m_{q} \gg \Lambda_{\rm QCD}/N_c^{3/2}$ the conventional large-$N_c$ counting holds. The physical values of quark masses for QCD correspond to the small quark mass limit. We find pion-nucleon coupling strengths are reduced to order ${\cal O}(1)$ rather than ${\cal O}(N_c)$. Under the assumption that in the large $N_c$ limit the sigma meson mass is larger than that of the omega, and that the omega-nucleon coupling constant is larger than that of the sigma, we argue that the nucleon-nucleon large range potential is weakly attractive and admits an interaction energy of order $\Lambda_{\rm QCD}/N_c^{5/2} \sim 10$ MeV. With these assumptions on coupling and masses, there is no strong long range attractive channel for nucleon-nucleon interactions, so that nuclear matter at densities much smaller than that where nucleons strongly interact is a weakly interacting configuration of nucleons with strongly interacting localized cores.


title: Low energy constituent quark and pion effective couplings in a weak external magnetic field

abstract: An effective model with pions and constituent quarks in the presence of a weak external background electromagnetic field is derived by starting from a dressed one gluon exchange quark-quark interaction.By applying the auxiliary field and background field methods, the structureless pion limit is considered to extract effective pion and constituent quark couplings in the presence of a weak magnetic field. The leading terms of a large quark and gluon masses expansion are obtained by resolving effective coupling constants which turn out to depend on a weak magnetic field. Two pion field definitions are considered for that. Several relations between the effective coupling constants and parameters can be derived exactly or in the limit of very large quark mass at zero and weak constant magnetic field. Among these ratios, the Gell Mann Oakes Renner and the quark level Goldberger Treiman relations are obtained. In addition to that, in the pion sector, the leading terms of Chiral Perturbation Theory coupled to the electromagnetic field is recovered. Some numerical estimates are provided for the effective coupling constants and parameters.


title: Lagrangian Reachabililty

abstract: We introduce LRT, a new Lagrangian-based ReachTube computation algorithm that conservatively approximates the set of reachable states of a nonlinear dynamical system. LRT makes use of the Cauchy-Green stretching factor (SF), which is derived from an over-approximation of the gradient of the solution flows. The SF measures the discrepancy between two states propagated by the system solution from two initial states lying in a well-defined region, thereby allowing LRT to compute a reachtube with a ball-overestimate in a metric where the computed enclosure is as tight as possible. To evaluate its performance, we implemented a prototype of LRT in C++/Matlab, and ran it on a set of well-established benchmarks. Our results show that LRT compares very favorably with respect to the CAPD and Flow* tools.


title: Are eikonal quasinormal modes linked to the unstable circular null geodesics?

abstract: In Phys. Rev. D 79, 064016 (2009) it was claimed that quasinormal modes which any stationary, spherically symmetric and asymptotically flat black hole emits in the eikonal regime are determined by the parameters of the circular null geodesic: the real and imaginary parts of the quasinormal mode are multiples of the frequency and instability timescale of the circular null geodesics respectively. We shall consider asymptotically flat black hole in the Einstein-Lovelock theory, find analytical expressions for gravitational quasinormal modes in the eikonal regime and analyze the null geodesics. Comparison of the both phenomena shows that the expected link between the null geodesics and quasinormal modes is violated in the Einstein-Lovelock theory. Nevertheless, the correspondence exists for a number of other cases and here we formulate its actual limits.


title: The Integer Cuboid Table

abstract: Integer cuboids are rectangular Diophantine parallelepipeds It has been discovered that these cuboids come in 3 varieties: Euler or body type, edge type, and face type. In all three cases, one edge or diagonal is irrational, all six others are rational. We discuss an exhaustive computer search procedure which uses the Pythagorean group Py(n) to locate all possible cuboids with a given edge n. Over the range of 44 to 200,000,000,027 for the smallest edge, 167,043 cuboids were discovered. They are listed in the Integer Cuboid Table.


title: Univalent functions in model spaces: revisited

abstract: Motivated by a problem in approximation theory, we find a necessary and sufficient condition for a model (backward shift invariant) subspace $K_\varTheta = H^2\ominus \varTheta H^2$ of the Hardy space $H^2$ to contain a bounded univalent function.


title: Universal short time quantum critical dynamics of finite size systems

abstract: We investigate the short time quantum critical dynamics in the imaginary time relaxation processes of finite size systems. Universal scaling behaviors exist in the imaginary time evolution and in particular, the system undergoes a critical initial slip stage characterized by an exponent $\theta$, in which an initial power-law increase emerges in the imaginary time correlation function when the initial state has zero order parameter and vanishing correlation length. Under different initial conditions, the quantum critical point and critical exponents can be determined from the universal scaling behaviors. We apply the method to the one- and two-dimensional transverse field Ising models using quantum Monte Carlo simulations. In the one-dimensional case, we locate the quantum critical point at $(h/J)_{c}=1.00003(8)$ \thirdrevise{in the thermodynamic limit}, and estimate the critical initial slip exponent $\theta=0.3734(2)$, static exponent $\beta/\nu=0.1251(2)$ \thirdrevise{by analyzing data on chains of length $L=32\sim 256$ and $L=48\sim 256$, respectively}. For the two-dimensional square-lattice system, the critical coupling ratio is given by $3.04451(7)$ \thirdrevise{in the thermodynamic limit} while the critical exponents are \thirdrevise{$\theta=0.209(4)$ and $\beta/\nu=0.518(1)$ estimated by data on systems of size $L=24\sim 64$ and $L=32\sim 64$, correspondingly.} Remarkably, the critical initial slip exponents obtained in both models are notably distinct from their classical counterparts, owing to the essential differences between classical and quantum dynamics. The short time critical dynamics and the imaginary time relaxation QMC approach can be readily adapted to various models.


title: Free fermions and the classical compact groups

abstract: There is a close connection between the ground state of non-interacting fermions in a box with classical (absorbing, reflecting, and periodic) boundary conditions and the eigenvalue statistics of the classical compact groups. The associated determinantal point processes can be extended in two natural directions: i) we consider the full family of admissible quantum boundary conditions (i.e., self-adjoint extensions) for the Laplacian on a bounded interval, and the corresponding projection correlation kernels; ii) we construct the grand canonical extensions at finite temperature of the projection kernels, interpolating from Poisson to random matrix eigenvalue statistics. The scaling limits in the bulk and at the edges are studied in a unified framework, and the question of universality is addressed. Whether the finite temperature determinantal processes correspond to the eigenvalue statistics of some matrix models is, a priori, not obvious. We complete the picture by constructing a finite temperature extension of the Haar measure on the classical compact groups. The eigenvalue statistics of the resulting grand canonical matrix models (of random size) corresponds exactly to the grand canonical measure of non-interacting free fermions with classical boundary conditions.

categories: math-ph cond-mat.stat-mech math.MP math.PR


title: Sub-sampled Cubic Regularization for Non-convex Optimization

abstract: We consider the minimization of non-convex functions that typically arise in machine learning. Specifically, we focus our attention on a variant of trust region methods known as cubic regularization. This approach is particularly attractive because it escapes strict saddle points and it provides stronger convergence guarantees than first- and second-order as well as classical trust region methods. However, it suffers from a high computational complexity that makes it impractical for large-scale learning. Here, we propose a novel method that uses sub-sampling to lower this computational cost. By the use of concentration inequalities we provide a sampling scheme that gives sufficiently accurate gradient and Hessian approximations to retain the strong global and local convergence guarantees of cubically regularized methods. To the best of our knowledge this is the first work that gives global convergence guarantees for a sub-sampled variant of cubic regularization on non-convex functions. Furthermore, we provide experimental results supporting our theory.


title: Analytic techniques for option pricing under a hyperexponential L\'{e}vy model

abstract: We develop series expansions in powers of $q^{-1}$ and $q^{-1/2}$ of solutions of the equation $\psi(z) = q$, where $\psi(z)$ is the Laplace exponent of a hyperexponential L\'{e}vy process. As a direct consequence we derive analytic expressions for the prices of European call and put options and their Greeks (Theta, Delta, and Gamma) and a full asymptotic expansion of the short-time Black-Scholes at-the-money implied volatility. Further we demonstrate how the speed of numerical algorithms for pricing exotic options, which are based on the Laplace transform, may be increased.


title: Rise of the humanbot

abstract: The accelerated path of technological development, particularly at the interface between hardware and biology has been suggested as evidence for future major technological breakthroughs associated to our potential to overcome biological constraints. This includes the potential of becoming immortal, having expanded cognitive capacities thanks to hardware implants or the creation of intelligent machines. Here I argue that several relevant evolutionary and structural constraints might prevent achieving most (if not all) these innovations. Instead, the coming future will bring novelties that will challenge many other aspects of our life and that can be seen as other feasible singularities. One particularly important one has to do with the evolving interactions between humans and non-intelligent robots capable of learning and communication. Here I argue that a long term interaction can lead to a new class of "agent" (the humanbot). The way shared memories get tangled over time will inevitably have important consequences for both sides of the pair, whose identity as separated entities might become blurred and ultimately vanish. Understanding such hybrid systems requires a second-order neuroscience approach while posing serious conceptual challenges, including the definition of consciousness.


title: Steady Prandtl Layers over a Moving Boundary: Non-Shear Euler flows

abstract: In this article we establish the validity of Prandtl layer expansions around Euler flows which are not shear. The presence of non-shear flows at the leading order creates a singularity of $o(\frac{1}{\sqrt{\epsilon}})$. A new $y$-weighted positivity estimate is developed to control this leading-order growth at the far field.


title: Engineering Record And Replay For Deployability: Extended Technical Report

abstract: The ability to record and replay program executions with low overhead enables many applications, such as reverse-execution debugging, debugging of hard-to-reproduce test failures, and "black box" forensic analysis of failures in deployed systems. Existing record-and-replay approaches limit deployability by recording an entire virtual machine (heavyweight), modifying the OS kernel (adding deployment and maintenance costs), requiring pervasive code instrumentation (imposing significant performance and complexity overhead), or modifying compilers and runtime systems (limiting generality). We investigated whether it is possible to build a practical record-and-replay system avoiding all these issues. The answer turns out to be yes - if the CPU and operating system meet certain non-obvious constraints. Fortunately modern Intel CPUs, Linux kernels and user-space frameworks do meet these constraints, although this has only become true recently. With some novel optimizations, our system 'rr' records and replays real-world low-parallelism workloads with low overhead, with an entirely user-space implementation, using stock hardware, compilers, runtimes and operating systems. "rr" forms the basis of an open-source reverse-execution debugger seeing significant use in practice. We present the design and implementation of 'rr', describe its performance on a variety of workloads, and identify constraints on hardware and operating system design required to support our approach.


title: Unusual structures inherent in point pattern data predict colon cancer patient survival

abstract: Cancer patient diagnosis and prognosis is informed by assessment of morphological properties observed in patient tissue. Pathologists normally carry out this assessment, yet advances in computational image analysis provide opportunities for quantitative assessment of tissue. A key aspect of that quantitative assessment is the development of algorithms able to link image data to patient survival. Here, we develop a point process methodology able to describe patterns in cell distribution within cancerous tissue samples. In particular, we consider the Palm intensities of two Neyman Scott point processes, and a void process developed herein to reflect the spatial patterning of the cells. An approximate-likelihood technique is taken in order to fit point process models to patient data and the predictive performance of each model is determined. We demonstrate that based solely on the spatial arrangement of cells we are able to predict patient survival.


title: Pressure tuning of light-induced superconductivity in K3C60

abstract: Optical excitation at terahertz frequencies has emerged as an effective means to manipulate complex solids dynamically. In the molecular solid K3C60, coherent excitation of intramolecular vibrations was shown to transform the high temperature metal into a non-equilibrium state with the optical conductivity of a superconductor. Here we tune this effect with hydrostatic pressure, and we find it to disappear around 0.3 GPa. Reduction with pressure underscores the similarity with the equilibrium superconducting phase of K3C60, in which a larger electronic bandwidth is detrimental for pairing. Crucially, our observation excludes alternative interpretations based on a high-mobility metallic phase. The pressure dependence also suggests that transient, incipient superconductivity occurs far above the 150 K hypothesised previously, and rather extends all the way to room temperature.


title: Subregular Complexity and Deep Learning

abstract: This paper argues that the judicial use of formal language theory and grammatical inference are invaluable tools in understanding how deep neural networks can and cannot represent and learn long-term dependencies in temporal sequences. Learning experiments were conducted with two types of Recurrent Neural Networks (RNNs) on six formal languages drawn from the Strictly Local (SL) and Strictly Piecewise (SP) classes. The networks were Simple RNNs (s-RNNs) and Long Short-Term Memory RNNs (LSTMs) of varying sizes. The SL and SP classes are among the simplest in a mathematically well-understood hierarchy of subregular classes. They encode local and long-term dependencies, respectively. The grammatical inference algorithm Regular Positive and Negative Inference (RPNI) provided a baseline. According to earlier research, the LSTM architecture should be capable of learning long-term dependencies and should outperform s-RNNs. The results of these experiments challenge this narrative. First, the LSTMs' performance was generally worse in the SP experiments than in the SL ones. Second, the s-RNNs out-performed the LSTMs on the most complex SP experiment and performed comparably to them on the others.


title: Modelling wave-induced sea ice breakup in the marginal ice zone

abstract: A model of ice floe breakup under ocean wave forcing in the marginal ice zone (MIZ) is proposed to investigate how floe size distribution (FSD) evolves under repeated wave breakup events. A three-dimensional linear model of ocean wave scattering by a finite array of compliant circular ice floes is coupled to a flexural failure model, which breaks a floe into two floes provided the two-dimensional stress field satisfies a breakup criterion. A closed-feedback loop algorithm is devised, which (i)~solves wave scattering problem for a given FSD under time-harmonic plane wave forcing, (ii)~computes the stress field in all the floes, (iii)~fractures the floes satisfying the breakup criterion and (iv)~generates an updated FSD, initialising the geometry for the next iteration of the loop.The FSD after 50 breakup events is uni-modal and near normal, or bi-modal. Multiple scattering is found to enhance breakup for long waves and thin ice, but to reduce breakup for short waves and thick ice. A breakup front marches forward in the latter regime, as wave-induced fracture weakens the ice cover allowing waves to travel deeper into the MIZ.


title: Inverse Multipath Fingerprinting for Millimeter Wave V2I Beam Alignment

abstract: Efficient beam alignment is a crucial component in millimeter wave systems with analog beamforming, especially in fast-changing vehicular settings. This paper proposes a position-aided approach where the vehicle's position (e.g., available via GPS) is used to query the multipath fingerprint database, which provides prior knowledge of potential pointing directions for reliable beam alignment. The approach is the inverse of fingerprinting localization, where the measured multipath signature is compared to the fingerprint database to retrieve the most likely position. The power loss probability is introduced as a metric to quantify misalignment accuracy and is used for optimizing candidate beam selection. Two candidate beam selection methods are developed, where one is a heuristic while the other minimizes the misalignment probability. The proposed beam alignment is evaluated using realistic channels generated from a commercial ray-tracing simulator. Using the generated channels, an extensive investigation is provided, which includes the required measurement sample size to build an effective fingerprint, the impact of measurement noise, the sensitivity to changes in traffic density, and beam alignment overhead comparison with IEEE 802.11ad as the baseline. Using the concept of beam coherence time, which is the duration between two consecutive beam alignments, and parameters of IEEE 802.11ad, the overhead is compared in the mobility context. The results show that while the proposed approach provides increasing rates with larger antenna arrays, IEEE 802.11ad has decreasing rates due to the larger beam training overhead that eats up a large portion of the beam coherence time, which becomes shorter with increasing mobility.


title: Banks as Tanks: A Continuous-Time Model of Financial Clearing

abstract: We present a simple continuous-time model of clearing in financial networks. Financial firms are represented as "tanks" filled with fluid (money), flowing in and out. Once "pipes" connecting "tanks" are open, the system reaches the clearing payment vector in finite time. This approach provides a simple recursive solution to a classical static model of financial clearing in bankruptcy, and suggests a practical payment mechanism. With sufficient resources, a system of mutual obligations can be restructured into an equivalent system that has a cascade structure: there is a group of banks that paid off their debts, another group that owes money only to banks in the first group, and so on. Technically, we use the machinery of Markov chains to analyze evolution of a deterministic dynamical system.


title: Electromagnetic properties of neutrinos: three new phenomena in neutrino spin oscillations

abstract: In studies of neutrino electromagnetic properties we discuss three very interesting aspects related to neutrino spin oscillations. First we consider neutrino mixing and oscillations in the mass and flavour bases under the influence of a constant magnetic field with nonzero transversal and longitudinal components. Then we discuss the effect of neutrino spin oscillations induced by electroweak interactions of neutrino with moving matter in case there is matter transversal current or polarization. In the final part of the paper we discuss recently developed approach to description of neutrino spin and spin-flavour oscillations in a constant magnetic field that is based on the use of the exact neutrino stationary states in the magnetic field.


title: Feasibility of the optical fiber clock

abstract: We explore the feasibility of a compact high-precision Hg atomic clock based on a hollow core optical fiber. We evaluate the sensitivity of the $^1S_0$-$^3P_0$ clock transition in Hg and other divalent atoms to the fiber inner core surface at non-zero temperatures. The Casimir-Polder interaction induced $^1S_0$-$^3P_0$ transition frequency shift is calculated for the atom inside the hollow capillary as a function of atomic position, capillary material, and geometric parameters. For $^{199}\mathrm{Hg}$ atoms on the axis of a silica capillary with inner radius $\geq 15 \,\mu \mathrm{m}$ and optimally chosen thickness $d\sim 1 \,\mu \mathrm{m}$, the atom-surface interaction induced $^1S_0$-$^3P_0$ clock transition frequency shift can be kept on the level $\delta\nu/\nu_{\mathrm{Hg}} \sim10^{-19}$. We also estimate the atom loss and heating due to the collisions with the buffer gas, lattice intensity noise induced heating, spontaneous photon scattering, and residual birefringence induced frequency shifts.


title: Optimal segmentation of directed graph and the minimum number of feedback arcs

abstract: The minimum feedback arc set problem asks to delete a minimum number of arcs (directed edges) from a digraph (directed graph) to make it free of any directed cycles. In this work we approach this fundamental cycle-constrained optimization problem by considering a generalized task of dividing the digraph into D layers of equal size. We solve the D-segmentation problem by the replica-symmetric mean field theory and belief-propagation heuristic algorithms. The minimum feedback arc density of a given random digraph ensemble is then obtained by extrapolating the theoretical results to the limit of large D. A divide-and-conquer algorithm (nested-BPR) is devised to solve the minimum feedback arc set problem with very good performance and high efficiency.


title: Optimal Resource Allocation for Power-Efficient MC-NOMA with Imperfect Channel State Information

abstract: In this paper, we study power-efficient resource allocation for multicarrier non-orthogonal multiple access (MC-NOMA) systems. The resource allocation algorithm design is formulated as a non-convex optimization problem which jointly designs the power allocation, rate allocation, user scheduling, and successive interference cancellation (SIC) decoding policy for minimizing the total transmit power. The proposed framework takes into account the imperfection of channel state information at transmitter (CSIT) and quality of service (QoS) requirements of users. To facilitate the design of optimal SIC decoding policy on each subcarrier, we define a channel-to-noise ratio outage threshold. Subsequently, the considered non-convex optimization problem is recast as a generalized linear multiplicative programming problem, for which a globally optimal solution is obtained via employing the branch-and-bound approach. The optimal resource allocation policy serves as a system performance benchmark due to its high computational complexity. To strike a balance between system performance and computational complexity, we propose a suboptimal iterative resource allocation algorithm based on difference of convex programming. Simulation results demonstrate that the suboptimal scheme achieves a close-to-optimal performance. Also, both proposed schemes provide significant transmit power savings than that of conventional orthogonal multiple access (OMA) schemes.


title: On angled bounce-off impact of a drop impinging on a flowing soap film

abstract: Small drops impinging angularly on thin flowing soap films frequently demonstrate the rare emergence of bulk elastic effects working in-tandem with the more common-place hydrodynamic interactions. Three collision regimes are observable: (a) drop piercing through the film, (b) it coalescing with the flow, and (c) it bouncing off the film surface. During impact, the drop deforms along with a bulk elastic deformation of the film. For impacts that are close-to-tangential, the bounce-off regime predominates. We outline a reduced order analytical framework assuming a deformable drop and a deformable three-dimensional film, and the idealization invokes a phase-based parametric study. Angular inclination of the film and the ratio of post and pre impact drop sizes entail the phase parameters. We also perform experiments with vertically descending droplets impacting against an inclined soap film, flowing under constant pressure head. Model predicted phase domain for bounce-off compares well to our experimental findings. Additionally, the experiments exhibit momentum transfer to the film in the form of shed vortex dipole, along with propagation of free surface waves. On consulting prior published work, we note that for locomotion of water-walking insects using an impulsive action, the momentum distribution to the shed vortices and waves are both significant, taking up respectively 2/3-rd and 1/3-rd of the imparted streamwise momentum. In view of the potentially similar impulse actions, this theory is applied to the bounce-off examples in our experiments, and the resultant shed vortex dipole momenta are compared to the momenta computed from particle imaging velocimetry data. The magnitudes reveal identical order ($10^{-7}$ N$\cdot$s), suggesting that the bounce-off regime can be tapped as a simple analogue for interfacial bio-locomotion relying on impulse reactions.


title: Invariance of topological indices under Hilbert space truncation

abstract: We show that the topological index of a wavefunction, computed in the space of twisted boundary phases, is preserved under Hilbert space truncation, provided the truncated state remains normalizable. If truncation affects the boundary condition of the resulting state, the invariant index may acquire a different physical interpretation. If the index is symmetry protected, the truncation should preserve the protecting symmetry. We discuss implications of this invariance using paradigmatic integer and fractional Chern insulators, $Z_2$ topological insulators, and Spin-$1$ AKLT and Heisenberg chains, as well as its relation with the notion of bulk entanglement. As a possible application, we propose a partial quantum tomography scheme from which the topological index of a generic multi-component wavefunction can be extracted by measuring only a small subset of wavefunction components, equivalent to the measurement of a bulk entanglement topological index.


title: Kernel clustering: density biases and solutions

abstract: Kernel methods are popular in clustering due to their generality and discriminating power. However, we show that many kernel clustering criteria have density biases theoretically explaining some practically significant artifacts empirically observed in the past. For example, we provide conditions and formally prove the density mode isolation bias in kernel K-means for a common class of kernels. We call it Breiman's bias due to its similarity to the histogram mode isolation previously discovered by Breiman in decision tree learning with Gini impurity. We also extend our analysis to other popular kernel clustering methods, e.g. average/normalized cut or dominant sets, where density biases can take different forms. For example, splitting isolated points by cut-based criteria is essentially the sparsest subset bias, which is the opposite of the density mode bias. Our findings suggest that a principled solution for density biases in kernel clustering should directly address data inhomogeneity. We show that density equalization can be implicitly achieved using either locally adaptive weights or locally adaptive kernels. Moreover, density equalization makes many popular kernel clustering objectives equivalent. Our synthetic and real data experiments illustrate density biases and proposed solutions. We anticipate that theoretical understanding of kernel clustering limitations and their principled solutions will be important for a broad spectrum of data analysis applications across the disciplines.


title: Optimal Ballistic Transport and Hopf-Lax Formulae on Wasserstein Space

abstract: We investigate the optimal mass transport problem associated to the following "ballistic" cost functional on phase space $M\times M^*$, $$ b_T(v, x):=\inf\{\langle v, \gamma (0)\rangle +\int_0^TL(\gamma (t), {\dot \gamma}(t))\, dt, \gamma \in C^1([0, T), M), \gamma(T)=x\}, $$ where $M=\mathbb{R}^d$, $T>0$, and $L:M\times M \to \mathbb{R}$ is a Lagrangian that is jointly convex in both variables. Under suitable conditions on the initial and final probability measures, we use convex duality \`a la Bolza and Monge-Kantorovich theory to lift classical Hopf-Lax formulae from state space to Wasserstein space. This allows us to relate optimal transport maps for the ballistic cost to those associated with the fixed-end cost defined on $M\times M$ by $$ c_T(x,y):=\inf\{\int_0^TL(\gamma(t), {\dot \gamma}(t))\, dt, \gamma\in C^1([0, T), M), \gamma(0)=x, \gamma(T)=y\}. $$ We also point to links with the theory of mean field games.


title: A Novel Neural Network Model for Joint POS Tagging and Graph-based Dependency Parsing

abstract: We present a novel neural network model that learns POS tagging and graph-based dependency parsing jointly. Our model uses bidirectional LSTMs to learn feature representations shared for both POS tagging and dependency parsing tasks, thus handling the feature-engineering problem. Our extensive experiments, on 19 languages from the Universal Dependencies project, show that our model outperforms the state-of-the-art neural network-based Stack-propagation model for joint POS tagging and transition-based dependency parsing, resulting in a new state of the art. Our code is open-source and available together with pre-trained models at: https://github.com/datquocnguyen/jPTDP


title: LoRa Backscatter: Enabling The Vision of Ubiquitous Connectivity

abstract: The vision of embedding connectivity into billions of everyday objects runs into the reality of existing communication technologies --- there is no existing wireless technology that can provide reliable and long-range communication at tens of microwatts of power as well as cost less than a dime. While backscatter is low-power and low-cost, it is known to be limited to short ranges. This paper overturns this conventional wisdom about backscatter and presents the first wide-area backscatter system. Our design can successfully backscatter from any location between an RF source and receiver, separated by 475 m, while being compatible with commodity LoRa hardware. Further, when our backscatter device is co-located with the RF source, the receiver can be as far as 2.8 km away. We deploy our system in a 4,800 $ft^{2}$ (446 $m^{2}$) house spread across three floors, a 13,024 $ft^{2}$ (1210 $m^{2}$) office area covering 41 rooms, as well as a one-acre (4046 $m^{2}$) vegetable farm and show that we can achieve reliable coverage, using only a single RF source and receiver. We also build a contact lens prototype as well as a flexible epidermal patch device attached to the human skin. We show that these devices can reliably backscatter data across a 3,328 $ft^{2}$ (309 $m^{2}$) room. Finally, we present a design sketch of a LoRa backscatter IC that shows that it costs less than a dime at scale and consumes only 9.25 $\mu$W of power, which is more than 1000x lower power than LoRa radio chipsets.


title: Convergence Results on Pulse Coupled Oscillator Protocols in Locally Connected Networks

abstract: This work provides new insights on the convergence of a locally connected network of pulse coupled oscillator (PCOs) (i.e., a bio-inspired model for communication networks) to synchronous and desynchronous states, and their implication in terms of the decentralized synchronization and scheduling in communication networks. Bio-inspired techniques have been advocated by many as fault-tolerant and scalable alternatives to produce self-organization in communication networks. The PCO dynamics in particular have been the source of inspiration for many network synchronization and scheduling protocols. However, their convergence properties, especially in locally connected networks, have not been fully understood, prohibiting the migration into mainstream standards. This work provides further results on the convergence of PCOs in locally connected networks and the achievable convergence accuracy under propagation delays. For synchronization, almost sure convergence is proved for $3$ nodes and accuracy results are obtained for general locally connected networks whereas, for scheduling (or desynchronization), results are derived for locally connected networks with mild conditions on the overlapping set of maximal cliques. These issues have not been fully addressed before in the literature.


title: ExoMol Line List XXI: Nitric Oxide (NO)

abstract: Line lists for the ground electronic ground state for six major isotopologues of nitric oxide are presented. The line lists are constructed using empirical energy levels (and line positions) and high-level {\it ab inito} intensities. The energy levels were obtained using a combination of two approaches, from an effective Hamiltonian and from solving the rovibronic Schr\"{o}dinger equation variationally. The effective hamiltonian model was obtained through a fit to the experimental line positions of NO available in the literature for all six isotopologues using the programs SPFIT and SPCAT. The variational model was built through a least squares fit of the \textit{ab inito} potential and spin-orbit curves to the experimentally derived energies and experimental line positions of the main isotopologue only using the Duo program. The \textit{ab inito} potential energy, spin-orbit and dipole moment curves (PEC, SOC and DMC) are computed using high-level {\it ab inito} methods and the MARVEL method is used to obtain energies of NO from experimental transition frequencies. Each line list covers a wavenumber range from 0 - 40,000 \cm with approximately 22,000 rovibronic states and 2.3-2.6 million transitions extending to $J_{max} = 184.5$ and $v_{max} = 51$. Partition functions are also calculated up to a temperature of 5000 K. The calculated absorption line intensities at 296 K using these line lists show excellent agreement with those included in the HITRAN and HITEMP databases. The computed NO line lists are the most comprehensive to date, covering a wider wavenumber and temperature range compared to both the HITRAN and HITEMP databases. These line lists are also more accurate than those used in HITEMP. The full line lists are available from the CDS and ExoMol databases; data will also be available from CDMS.


title: Reduced Wigner coefficients for Lie superalgebra gl(m|n) corresponding to unitary representations and beyond

abstract: In this paper fundamental Wigner coefficients are determined algebraically by considering the eigenvalues of certain generalized Casimir invariants. Here this method is applied in the context of both type 1 and type 2 unitary representations of the Lie superalgebra gl(mjn). Extensions to the non-unitary case are investigated. A symmetry relation between two classes of Wigner coefficients is given in terms of a ratio of dimensions.


title: Differentially Private Release of Public Transport Data: The Opal Use Case

abstract: This document describes the application of a differentially private algorithm to release public transport usage data from Transport for New South Wales (TfNSW), Australia. The data consists of two separate weeks of "tap-on/tap-off" data of individuals who used any of the four different modes of public transport from TfNSW: buses, light rail, train and ferries. These taps are recorded through the smart ticketing system, known as Opal, available in the state of New South Wales, Australia.


title: Cartan Subalgebras for Quantum Symmetric Pair Coideals

abstract: There is renewed interest in the coideal subalgebras used to form quantum symmetric pairs because of recent discoveries showing that they play a fundamental role in the representation theory of quantized enveloping algebras. However, there is still no general theory of finite-dimensional modules for these coideals. In this paper, we establish an important step in this direction: we show that every quantum symmetric pair coideal subalgebra admits a quantum Cartan subalgebra which is a polynomial ring that specializes to its classical counterpart. The construction builds on Kostant and Sugiura's classification of Cartan subalgebras for real semisimple Lie algebras via strongly orthogonal systems of positive roots. We show that these quantum Cartan subalgebras act semisimply on finite-dimensional unitary modules and identify particularly nice generators of the quantum Cartan subalgebra for a family of examples.


title: Constraint Energy Minimizing Generalized Multiscale Finite Element Method in the Mixed Formulation

abstract: This paper presents a novel mass-conservative mixed multiscale method for solving flow equations in heterogeneous porous media. The media properties (the permeability) contain multiple scales and high contrast. The proposed method solves the flow equation in a mixed formulation on a coarse grid by constructing multiscale basis functions. The resulting velocity field is mass conservative on the fine grid. Our main goal is to obtain first-order convergence in terms of the mesh size which is independent of local contrast. This is achieved, first, by constructing some auxiliary spaces, which contain global information that can not be localized, in general. This is built on our previous work on the Generalized Multiscale Finite Element Method (GMsFEM). In the auxiliary space, multiscale basis functions corresponding to small (contrast-dependent) eigenvalues are selected. These basis functions represent the high-conductivity channels (which connect the boundaries of a coarse block). Next, we solve local problems to construct multiscale basis functions for the velocity field. These local problems are formulated in the oversampled domain taking into account some constraints with respect to auxiliary spaces. The latter allows fast spatial decay of local solutions and, thus, allows taking smaller oversampled regions. The number of basis functions depends on small eigenvalues of the local spectral problems. Moreover, multiscale pressure basis functions are needed in constructing the velocity space. Our multiscale spaces have a minimal dimension, which is needed to avoid contrast-dependence in the convergence. The method's convergence requires an oversampling of several layers. We present an analysis of our approach. Our numerical results confirm that the convergence rate is first order with respect to the mesh size and independent of the contrast.


title: Utility Maximizing Sequential Sensing Over a Finite Horizon

abstract: We consider the problem of optimally utilizing $N$ resources, each in an unknown binary state. The state of each resource can be inferred from state-dependent noisy measurements. Depending on its state, utilizing a resource results in either a reward or a penalty per unit time. The objective is a sequential strategy governing the decision of sensing and exploitation at each time to maximize the expected utility (i.e., total reward minus total penalty and sensing cost) over a finite horizon $L$. We formulate the problem as a Partially Observable Markov Decision Process (POMDP) and show that the optimal strategy is based on two time-varying thresholds for each resource and an optimal selection rule for which resource to sense. Since a full characterization of the optimal strategy is generally intractable, we develop a low-complexity policy that is shown by simulations to offer near optimal performance. This problem finds applications in opportunistic spectrum access, marketing strategies and other sequential resource allocation problems.


title: Classical causal models for Bell and Kochen-Specker inequality violations require fine-tuning

abstract: Nonlocality and contextuality are at the root of conceptual puzzles in quantum mechanics, and are key resources for quantum advantage in information-processing tasks. Bell nonlocality is best understood as the incompatibility between quantum correlations and the classical theory of causality, applied to relativistic causal structure. Contextuality, on the other hand, is on a more controversial foundation. In this work, I provide a common conceptual ground between nonlocality and contextuality as violations of classical causality. First, I show that Bell inequalities can be derived solely from the assumptions of no-signalling and no-fine-tuning of the causal model. This removes two extra assumptions from a recent result from Wood and Spekkens, and remarkably, does not require any assumption related to independence of measurement settings -- unlike all other derivations of Bell inequalities. I then introduce a formalism to represent contextuality scenarios within causal models and show that all classical causal models for violations of a Kochen-Specker inequality require fine-tuning. Thus the quantum violation of classical causality goes beyond the case of space-like separated systems, and manifests already in scenarios involving single systems.


title: Calibration of the NDHA model to describe N2O dynamics from respirometric assays

abstract: The NDHA model comprehensively describes nitrous oxide (N2O) producing pathways by both autotrophic ammonium oxidizing and heterotrophic bacteria. The model was calibrated via a set of targeted extant respirometric assays using enriched nitrifying biomass from a lab-scale reactor. Biomass response to ammonium, hydroxylamine, nitrite and N2O additions under aerobic and anaerobic conditions were tracked with continuous measurement of dissolved oxygen (DO) and N2O. The sequential addition of substrate pulses allowed the isolation of oxygen-consuming processes. The parameters to be estimated were determined by the information content of the datasets using identifiability analysis. Dynamic DO profiles were used to calibrate five parameters corresponding to endogenous, nitrite oxidation and ammonium oxidation processes. The subsequent N2O calibration was not significantly affected by the uncertainty propagated from the DO calibration because of the high accuracy of the estimates. Five parameters describing the individual contribution of three biological N2O pathways were estimated accurately (variance/mean < 10% for all estimated parameters). The NDHA model response was evaluated with statistical metrics (F-test, autocorrelation function). The 95% confidence intervals of DO and N2O predictions based on the uncertainty obtained during calibration are studied for the first time. The measured data fall within the 95% confidence interval of the predictions, indicating a good model description. Overall, accurate parameter estimation and identifiability analysis of ammonium removal significantly decreases the uncertainty propagated to N2O production, which is expected to benefit N2O model discrimination studies and reliable full scale applications.


title: Sharp bounds for the Randic index of graphs with given minimum and maximum degree

abstract: The Randi{\' c} index of a graph $G$, written $R(G)$, is the sum of $\frac 1{\sqrt{d(u)d(v)}}$ over all edges $uv$ in $E(G)$. %let $R(G)=\sum_{uv \in E(G)} \frac 1{\sqrt{d(u)d(v)}}$, which is called the Randi{\' c} index of it. Let $d$ and $D$ be positive integers $d < D$. In this paper, we prove that if $G$ is a graph with minimum degree $d$ and maximum degree $D$, then $R(G) \ge \frac{\sqrt{dD}}{d+D}n$; equality holds only when $G$ is an $n$-vertex $(d,D)$-biregular. Furthermore, we show that if $G$ is an $n$-vertex connected graph with minimum degree $d$ and maximum degree $D$, then $R(G) \le \frac n2- \sum_{i=d}^{D-1}\frac 12 \left( \frac 1{\sqrt{i}} - \frac 1{\sqrt{i+1}}\right)^2$; it is sharp for infinitely many $n$, and we characterize when equality holds in the bound.


title: Regularizing nonlinear Schroedinger equations through partial off-axis variations

abstract: We study a class of focusing nonlinear Schroedinger-type equations derived recently by Dumas, Lannes and Szeftel within the mathematical description of high intensity laser beams [7]. These equations incorporate the possibility of a (partial) off-axis variation of the group velocity of such laser beams through a second order partial differential operator acting in some, but not necessarily all, spatial directions. We study the well-posedness theory for such models and obtain a regularizing effect, even in the case of only partial off-axis dependence. This provides an answer to an open problem posed in [7].


title: Measurement of the directional sensitivity of DMTPC detectors

abstract: The Dark Matter Time Projection Chamber (DMTPC) is a direction-sensitive detector designed to measure the direction of recoiling $^{19}$F and $^{12}$C nuclei in low-pressure CF$_4$ gas using optical and charge readout systems. In this paper, we employ measurements from two DMTPC detectors, with operating pressures of 30-60 torr, to develop and validate a model of the directional response and performance of such detectors as a function of recoil energy. Using our model as a benchmark, we formulate the necessary specifications for a scalable directional detector with sensitivity comparable to that of current-generation counting (non-directional) experiments, which measure only recoil energy. Assuming the performance of existing DMTPC detectors, as well as current limits on the spin-dependent WIMP-nucleus cross section, we find that a 10-20 kg scale direction-sensitive detector is capable of correlating the measured direction of nuclear recoils with the predicted direction of incident dark matter particles and providing decisive (3$\sigma$) confirmation that a candidate signal from a non-directional experiment was indeed induced by elastic scattering of dark matter particles off of target nuclei.


title: Search for $CP$ Violation and Measurement of the Branching Fraction in the Decay $D^{0} \to K^0_S K^0_S$

abstract: We report a study of the decay $D^0 \to K^0_S K^0_S$ using 921~fb$^{-1}$ of data collected at or near the $\Upsilon(4S)$ and $\Upsilon(5S)$ resonances with the Belle detector at the KEKB asymmetric energy $e^+e^-$ collider. The measured time-integrated $CP$ asymmetry is $ A_{CP}(D^0 \to K^0_S K^0_S) = (-0.02 \pm 1.53 \pm 0.02 \pm 0.17) \%$, and the branching fraction is $\mathcal{B} (D^{0}\rightarrow K_{S}^{0}K_{S}^{0})$ = (1.321 $\pm$ 0.023 $\pm$ 0.036 $\pm$ 0.044) $\times$ 10$^{-4}$, where the first uncertainty is statistical, the second is systematic, and the third is due to the normalization mode ($D^0 \to K_S^0 \pi^0$). These results are significantly more precise than previous measurements available for this mode. The $A_{CP}$ measurement is consistent with the standard model expectation.


title: Community Detection for Multilayer Heterogeneous Network

abstract: Many real world networks consist of multiple types of nodes with edges that are heterogeneous in nature. However, most of the existing work for community detection only focused on homogeneous network consisting of a single layer. In this paper, we propose a modified Degree-Corrected Stochastic Model (DCBM) for modeling multilayer heterogeneous network. We develop a spectral clustering method that can unify the information contained in each sub-network, and demonstrate its efficiency to detect communities on simulated data and on Authorship/Citation network data. As a by-product, we present a novel algorithm called BiScore for clustering bipartite network under DCBM, and show that under mild conditions BiScore is guaranteed to yield consistent results.


title: Suppression of the Landau-Zener transition probability by a weak classical noise

abstract: When the drive which causes the level crossing in a qubit is slow, the probability, P_{LZ}, of the Landau-Zener transition is close to 1. We show that in this regime, which is most promising for applications, the noise due to the coupling to the environment, reduces the average P_{LZ}. At the same time, the survival probability, 1-P_{LZ}, which is exponentially small for a slow drive, can be completely dominated by noise-induced correction. Our main message is that the effect of a weak classical noise can be captured analytically by treating it as a perturbation in the Schroedinger equation. This allows us to study the dependence of the noise-induced correction to P_{LZ} on the correlation time of the noise. As this correlation time exceeds the bare Landau-Zener transition time, the effect of noise becomes negligible. We consider two conventional realizations of noise: Gaussian noise and telegraph noise.


title: An invitation to 2D TQFT and quantization of Hitchin spectral curves

abstract: This article consists of two parts. In Part 1, we present a formulation of two-dimensional topological quantum field theories in terms of a functor from a category of Ribbon graphs to the endofuntor category of a monoidal category. The key point is that the category of ribbon graphs produces all Frobenius objects. Necessary backgrounds from Frobenius algebras, topological quantum field theories, and cohomological field theories are reviewed. A result on Frobenius algebra twisted topological recursion is included at the end of Part 1. In Part 2, we explain a geometric theory of quantum curves. The focus is placed on the process of quantization as a passage from families of Hitchin spectral curves to families of opers. To make the presentation simpler, we unfold the story using SL_2(\mathbb{C})-opers and rank 2 Higgs bundles defined on a compact Riemann surface $C$ of genus greater than $1$. In this case, quantum curves, opers, and projective structures in $C$ all become the same notion. Background materials on projective coordinate systems, Higgs bundles, opers, and non-Abelian Hodge correspondence are explained.


title: Changes in Ugandan Climate Rainfall at the Village and Forest Level

abstract: In 2013, the US National Oceanographic and Atmospheric Agency refined the historical rainfall estimates over the African Continent and produced the African Rainfall Climate version 2.0 (ARC2) estimator. ARC2 offers a nearly complete record of daily rainfall estimates since 1983 at 0.1 x 0.1 degree resolution. Despite short-term anomalies in twice-yearly rainy season intensities in Uganda, we identify an overall decrease in average rainfall of about 12% during the past 34 years. Spatiotemporally, these decreases are greatest in agricultural regions of central and western Uganda, but are also reflective of rainfall decreases in the gorilla habitat within the Bwindi Forest in Southwest Uganda. The findings carry significant implications for agriculture, food security, and wildlife habitat.


title: Ultra-low-loss Polaritons in Isotopically Pure Materials: A New Approach

abstract: Conventional optical components are limited to size-scales much larger than the wavelength of light, as changes in the amplitude, phase and polarization of the electromagnetic fields are accrued gradually along an optical path. However, advances in nanophotonics have produced ultra-thin, co-called "flat" optical components that beget abrupt changes in these properties over distances significantly shorter than the free space wavelength. While high optical losses still plague many approaches, phonon polariton (PhP) materials have demonstrated long lifetimes for sub-diffractional modes in comparison to plasmon-polariton-based nanophotonics. We experimentally observe a three-fold improvement in polariton lifetime through isotopic enrichment of hexagonal boron nitride (hBN). Commensurate increases in the polariton propagation length are demonstrated via direct imaging of polaritonic standing waves by means of infrared nano-optics. Our results provide the foundation for a materials-growth-directed approach towards realizing the loss control necessary for the development of PhP-based nanophotonic devices.