title: Star cluster disruption by a massive black hole binary

abstract: Massive black hole binaries (BHBs) are expected to form as the result of galaxy mergers; they shrink via dynamical friction and stellar scatterings, until gravitational waves (GWs) bring them to the final coalescence. It has been argued that BHBs may stall at a parsec scale and never enter the GW stage if stars are not continuously supplied to the BHB loss cone. Here we perform several N-body experiments to study the effect of an 80,000 solar masses stellar cluster (SC) infalling on a parsec-scale BHB. We explore different orbital elements for the SC and we perform runs both with and without accounting for the influence of a rigid stellar cusp (modelled as a rigid Dehnen potential). We find that the semi-major axis of the BHB shrinks by more than 10 per cent if the SC is on a nearly radial orbit; the shrinking is more efficient when a Dehnen potential is included and the orbital plane of the SC coincides with that of the BHB. In contrast, if the SC orbit has non-zero angular momentum, only a few stars enter the BHB loss cone and the resulting BHB shrinking is negligible. Our results indicate that SC disruption might significantly contribute to the shrinking of a parsec-scale BHB only if the SC approaches the BHB on a nearly radial orbit.


title: Reference Class Forecasting for Hong Kong's Major Roadworks Projects

abstract: Reference class forecasting is a method to remove optimism bias and strategic misrepresentation in infrastructure projects and programmes. In 2012 the Hong Kong government's Development Bureau commissioned a feasibility study on reference class forecasting in Hong Kong - a first for the Asia-Pacific region. This study involved 25 roadwork projects, for which forecast costs and durations were compared with actual outcomes. The analysis established and verified the statistical distribution of the forecast accuracy at various stages of project development, and benchmarked the projects against a sample of 863 similar projects. The study contributed to the understanding of how to improve forecasts by de-biasing early estimates, explicitly considering the risk appetite of decision makers, and safeguarding public funding allocation by balancing exceedance and under-use of project budgets.


title: SOPE: A Spatial Order Preserving Encryption Model for Multi-dimensional Data

abstract: Due to the increasing demand for cloud services and the threat of privacy invasion, the user is suggested to encrypt the data before it is outsourced to the remote server. The safe storage and efficient retrieval of d-dimensional data on an untrusted server has therefore crucial importance. The paper proposes a new encryption model which offers spatial order-preservation for d-dimensional data (SOPE model). The paper studies the operations for the construction of the encrypted database and suggests algorithms that exploit unique properties that this new model offers for the efficient execution of a whole range of well-known queries over the encrypted d-dimensional data. The new model utilizes well-known database indices, such as the B+-tree and the R-tree, as backbone structures in their traditional form, as it suggests no modifications to them for loading the data and for the efficient execution of the supporting query algorithms. An extensive experimental study that is also presented in the paper indicates the effectiveness and practicability of the proposed encryption model for real-life d-dimensional data applications.


title: Electrokinetic Lattice Boltzmann solver coupled to Molecular Dynamics: application to polymer translocation

abstract: We develop a theoretical and computational approach to deal with systems that involve a disparate range of spatio-temporal scales, such as those comprised of colloidal particles or polymers moving in a fluidic molecular environment. Our approach is based on a multiscale modeling that combines the slow dynamics of the large particles with the fast dynamics of the solvent into a unique framework. The former is numerically solved via Molecular Dynamics and the latter via a multi-component Lattice Boltzmann. The two techniques are coupled together to allow for a seamless exchange of information between the descriptions. Being based on a kinetic multi-component description of the fluid species, the scheme is flexible in modeling charge flow within complex geometries and ranging from large to vanishing salt concentration. The details of the scheme are presented and the method is applied to the problem of translocation of a charged polymer through a nanopores. In the end, we discuss the advantages and complexities of the approach.


title: Setting the threshold for high throughput detectors: A mathematical approach for ensembles of dynamic, heterogeneous, probabilistic anomaly detectors

abstract: Anomaly detection (AD) has garnered ample attention in security research, as such algorithms complement existing signature-based methods but promise detection of never-before-seen attacks. Cyber operations manage a high volume of heterogeneous log data; hence, AD in such operations involves multiple (e.g., per IP, per data type) ensembles of detectors modeling heterogeneous characteristics (e.g., rate, size, type) often with adaptive online models producing alerts in near real time. Because of high data volume, setting the threshold for each detector in such a system is an essential yet underdeveloped configuration issue that, if slightly mistuned, can leave the system useless, either producing a myriad of alerts and flooding downstream systems, or giving none. In this work, we build on the foundations of Ferragut et al. to provide a set of rigorous results for understanding the relationship between threshold values and alert quantities, and we propose an algorithm for setting the threshold in practice. Specifically, we give an algorithm for setting the threshold of multiple, heterogeneous, possibly dynamic detectors completely a priori, in principle. Indeed, if the underlying distribution of the incoming data is known (closely estimated), the algorithm provides provably manageable thresholds. If the distribution is unknown (e.g., has changed over time) our analysis reveals how the model distribution differs from the actual distribution, indicating a period of model refitting is necessary. We provide empirical experiments showing the efficacy of the capability by regulating the alert rate of a system with $\approx$2,500 adaptive detectors scoring over 1.5M events in 5 hours. Further, we demonstrate on the real network data and detection framework of Harshaw et al. the alternative case, showing how the inability to regulate alerts indicates the detection model is a bad fit to the data.


title: Uniform Circle Formation by Transparent Fat Robots

abstract: This paper addresses the problem of Uniform Circle Formation by n > 1 transparent disc robots (fat robots). The robots execute repetitive cycles of the states look-compute-move in semi-synchronous manner where a set of robots execute the cycle simultaneously. They do not communicate by any explicit message passing. However, they can sense or observes the positions of other robots around themselves through sensors or camera. The robots are unable to recover the past actions and observations. They have no unique identity. The robots do not have any global coordinate system. They agree upon only y-axis (South and North direction). However, they do not have chirality or common orientation of Y axis with respect to X axis. Being transparent the robots do not cause any visual obstructions for other robots. But, they act as physical obstacles for other robots. This paper proposes a collision free movement strategy for the robots to form a uniform circle (in other words convex regular polygon) executing finite number of cycles. To the best of our knowledge this is the first reported results on uniform circle formation for fat robots under the considered model.


title: Radiation Forces on Dust Envelopes

abstract: We address in detail the radiation forces on spherical dust envelopes around luminous stars, and numerical solutions for these forces, as a first step toward more general dust geometries. Two physical quantities, a normalized force and a force-averaged radius, suffice to capture the overall effects of radiation forces; these combine to yield the radiation term in the virial theorem. In addition to the optically thin and thick regimes, the wavelength dependence of dust opacity allows for an intermediate case in which starlight is easily trapped but infrared radiation readily escapes. Scattering adds a non-negligible force in this intermediate regime. We address all three regimes analytically and provide approximate formulae for the force parameters, for arbitrary optical depth and inner dust temperature. Turning to numerical codes, we examine the convergence properties of the Monte Carlo code Hyperion run in Cartesian geometry. We calibrate both Hyperion and our analytical estimate using the DUSTY code, run in spherical geometry. We find that Monte Carlo codes tend to underestimate the radiation force when the mean free path of starlight is not well resolved, as this causes the inner dust temperature, and therefore the inner Rosseland opacity, to be too low. We briefly discuss implications for more complicated radiation transfer problems.


title: Decoding the phase structure of QCD via particle production at high energy

abstract: Recent studies based on non-perturbative lattice Monte-Carlo solutions of Quantum Chromodynamics, the theory of strong interactions, demonstrated that at high temperature there is a phase change from confined hadronic matter to a deconfined quark-gluon plasma where quarks and gluons can travel distances largely exceeding the size of hadrons. The phase structure of such strongly interacting matter can be decoded via analysis of particle abundances in high energy nuclear collisions within the framework of the statistical hadronization approach. The results imply quark-hadron duality at and experimental delineation of the location of the phase boundary of strongly interacting matter.


title: Global continuity and BMO estimates for non-Newtonian fluids with perfect slip boundary conditions

abstract: We study the generalized stationary Stokes system in a bounded domain in the plane equipped with perfect slip boundary conditions. We show natural stability results in oscillatory spaces, i.e. H\"older spaces and Campanato spaces including the border-line spaces of bounded mean oscillations (BMO) and vanishing mean oscillations (VMO). In particular, we show that under appropriate assumptions gradients of solutions are globally continuous. Since the stress tensor is assumed to be governed by a general Orlicz function, our theory includes various cases of (possibly degenerate) shear thickening and shear thinning fluids; including the model case of power law fluids. The global estimates seem to be new even in case of the linear Stokes system.


title: On the nonintegrability of equations for long- and short-wave interactions

abstract: We examine the integrability of two models used for the interaction of long and short waves in dispersive media. One is more classical but arguably cannot be derived from the underlying water wave equations, while the other one was recently derived. We use the method of Zakharov and Schulman to attempt to construct conserved quantities for these systems at different orders in the magnitude of the solutions. The coupled KdV-NLS model is shown to be nonintegrable, due to the presence of fourth-order resonances. A coupled real KdV - complex KdV system is shown to suffer the same fate, except for three special choices of the coefficients, where higher-order calculations or a different approach are necessary to conclude integrability or the absence thereof.


title: 12C within the Semimicroscopic Algebraic Cluster Model

abstract: The Semimicroscopic Algebraic Cluster Model (SACM) is applied to 12C as a system of three alpha- clusters. The microscopic model space, which observes the Pauli-Exclusion-Principle (PEP), is constructed. It is shown that the 12C nucleus can effectively be treated as a two-cluster system 8Be+alpha. The experimental spectrum is well reproduced. The geometrical mapping is discussed and it is shown that the ground state must correspond to a triangular structure, which is in agreement with other microscopic calculations. The non-zero B(E2; 0_2+ --> 2_1+) transition requires a mixing of SU(3) irreducible representations (irreps) whose consequences are discussed. The Hoyle state turns out to contain large shell excitations. The results are compared to another phenomenological model, which assumes a triangular structure and, using simple symmetry arguments, can reproduce the states observed at low energy. This model does not observe the PEP and one objective of our contribution is to verify the extend of importance of the PEP.


title: DPCA: Dimensionality Reduction for Discriminative Analytics of Multiple Large-Scale Datasets

abstract: Principal component analysis (PCA) has well-documented merits for data extraction and dimensionality reduction. PCA deals with a single dataset at a time, and it is challenged when it comes to analyzing multiple datasets. Yet in certain setups, one wishes to extract the most significant information of one dataset relative to other datasets. Specifically, the interest may be on identifying, namely extracting features that are specific to a single target dataset but not the others. This paper develops a novel approach for such so-termed discriminative data analysis, and establishes its optimality in the least-squares (LS) sense under suitable data modeling assumptions. The criterion reveals linear combinations of variables by maximizing the ratio of the variance of the target data to that of the remainders. The novel approach solves a generalized eigenvalue problem by performing SVD just once. Numerical tests using synthetic and real datasets showcase the merits of the proposed approach relative to its competing alternatives.


title: A Markov Chain Theory Approach to Characterizing the Minimax Optimality of Stochastic Gradient Descent (for Least Squares)

abstract: This work provides a simplified proof of the statistical minimax optimality of (iterate averaged) stochastic gradient descent (SGD), for the special case of least squares. This result is obtained by analyzing SGD as a stochastic process and by sharply characterizing the stationary covariance matrix of this process. The finite rate optimality characterization captures the constant factors and addresses model mis-specification.


title: On the Long-Term Memory of Deep Recurrent Networks

abstract: A key attribute that drives the unprecedented success of modern Recurrent Neural Networks (RNNs) on learning tasks which involve sequential data, is their ability to model intricate long-term temporal dependencies. However, a well established measure of RNNs long-term memory capacity is lacking, and thus formal understanding of the effect of depth on their ability to correlate data throughout time is limited. Specifically, existing depth efficiency results on convolutional networks do not suffice in order to account for the success of deep RNNs on data of varying lengths. In order to address this, we introduce a measure of the network's ability to support information flow across time, referred to as the Start-End separation rank, which reflects the distance of the function realized by the recurrent network from modeling no dependency between the beginning and end of the input sequence. We prove that deep recurrent networks support Start-End separation ranks which are combinatorially higher than those supported by their shallow counterparts. Thus, we establish that depth brings forth an overwhelming advantage in the ability of recurrent networks to model long-term dependencies, and provide an exemplar of quantifying this key attribute which may be readily extended to other RNN architectures of interest, e.g. variants of LSTM networks. We obtain our results by considering a class of recurrent networks referred to as Recurrent Arithmetic Circuits, which merge the hidden state with the input via the Multiplicative Integration operation, and empirically demonstrate the discussed phenomena on common RNNs. Finally, we employ the tool of quantum Tensor Networks to gain additional graphic insight regarding the complexity brought forth by depth in recurrent networks.


title: Gravitational Wave Signatures of Highly Compact Boson Star Binaries

abstract: Solitonic boson stars are stable objects made of a complex scalar field with a compactness that can reach values comparable to that of neutron stars. A recent study of the collision of identical boson stars produced only non-rotating boson stars or black holes, suggesting that rotating boson stars may not form from binary mergers. Here we extend this study to include an analysis of the gravitational waves radiated during the coalescence of such a binary, which is crucial to distinguish these events from other binaries with LIGO and Virgo observations. Our studies reveal that the remnant's gravitational wave signature is mainly governed by its fundamental frequency as it settles down to a non-rotating boson star, emitting significant gravitational radiation during this post-merger state. We calculate how the waveforms and their post-merger frequencies depend on the compactness of the initial boson stars and estimate analytically the amount of energy radiated after the merger.


title: Radiative cooling of swept up gas in AGN-driven galactic winds and its implications for molecular outflows

abstract: We recently used hydro-chemical simulations to demonstrate that molecular outflows observed in luminous quasars can be explained by molecule formation within the AGN wind. However, these simulations cover a limited parameter space, due to their computational cost. We have therefore developed an analytic model to follow cooling in the shocked ISM layer of an AGN wind. We explore different ambient densities ($1-10^{4} \, \rm{cm}^{-3}$), density profile slopes ($0-1.5$), AGN luminosities ($10^{44}-10^{47} \, \rm{erg} \, \rm{s}^{-1}$), and metallicities ($0.1-3 \rm{Z}_{\odot}$). The swept up gas mostly cools within ~1 Myr. Based on our previous simulations, we predict that this gas would produce observable molecular outflows. The instantaneous momentum boost initially increases as the outflow decelerates. However, it reaches a maximum of $\approx$20, due to work done against the gravitational potential. The predicted time-averaged observational estimate of the molecular outflow momentum boost reaches a maximum of $\approx1-2$, partly due to our assumed molecular fraction, 0.2, but also because the instantaneous and observational, time-averaged definitions are not equivalent. Thus recent observational estimates of order unity momentum boosts do not necessarily rule out energy-driven outflows. Finally, we find that dust grains are likely to re-form by accretion of metals after the shocked ISM layer has cooled, assuming that a small fraction of dust grains swept up after this layer has cooled are able to mix into the cool phase, and assuming that grain growth remains efficient in the presence of the strong AGN radiation field. This would enable rapid molecule formation, as assumed in our models.


title: Chromatic numbers of stable Kneser hypergraphs via topological Tverberg-type theorems

abstract: Kneser's 1955 conjecture -- proven by Lov\'asz in 1978 -- asserts that in any partition of the $k$-subsets of $\{1, 2, \dots, n\}$ into $n-2k-3$ parts, one part contains two disjoint sets. Schrijver showed that one can restrict to significantly fewer $k$-sets and still observe the same intersection pattern. Alon, Frankl, and Lov\'asz proved a different generalization of Kneser's conjecture for $r$ pairwise disjoint sets. Dolnikov generalized Lov\'asz' result to arbitrary set systems, while K\v{r}\'{i}\v{z} did the same for the $r$-fold extension of Kneser's conjecture. Here we prove a common generalization of all of these results. Moreover, we prove additional strengthenings by determining the chromatic number of certain sparse stable Kneser hypergraphs, and further develop a general approach to establishing lower bounds for chromatic numbers of hypergraphs using a combination of methods from equivariant topology and intersection results for convex hulls of points in Euclidean space.


title: Malware Detection by Eating a Whole EXE

abstract: In this work we introduce malware detection from raw byte sequences as a fruitful research area to the larger machine learning community. Building a neural network for such a problem presents a number of interesting challenges that have not occurred in tasks such as image processing or NLP. In particular, we note that detection from raw bytes presents a sequence problem with over two million time steps and a problem where batch normalization appear to hinder the learning process. We present our initial work in building a solution to tackle this problem, which has linear complexity dependence on the sequence length, and allows for interpretable sub-regions of the binary to be identified. In doing so we will discuss the many challenges in building a neural network to process data at this scale, and the methods we used to work around them.


title: The Geometry of Sagittarius Stream from Pan-STARRS1 3$\pi$ RR Lyrae

abstract: We present a comprehensive and precise description of the Sagittarius (Sgr) stellar stream's 3D geometry as traced by its old stellar population. This analysis draws on the sample of ${\sim}44,000$ RR Lyrae (RRab) stars from the Pan-STARRS1 (PS1) 3$\pi$ survey (Hernitschek et al. 2016,Sesar et al. 2017b), which is ${\sim}80\%$ complete and ${\sim}90\%$ pure within 80~kpc, and extends to ${\gtrsim} 120$~kpc with a distance precision of ${\sim} 3\%$. A projection of RR Lyrae stars within $|\tilde{B}|_{\odot}<9^\circ$ of the Sgr stream's orbital plane reveals the morphology of both the leading and the trailing arms at very high contrast, across much of the sky. In particular, the map traces the stream near-contiguously through the distant apocenters. We fit a simple model for the mean distance and line-of-sight depth of the Sgr stream as a function of the orbital plane angle $\tilde{\Lambda}_{\odot}$, along with a power-law background-model for the field stars. This modeling results in estimates of the mean stream distance precise to ${\sim}1\%$ and it resolves the stream's line-of-sight depth. These improved geometric constraints can serve as new constraints for dynamical stream models.


title: Casper the Friendly Finality Gadget

abstract: We introduce Casper, a proof of stake-based finality system which overlays an existing proof of work blockchain. Casper is a partial consensus mechanism combining proof of stake algorithm research and Byzantine fault tolerant consensus theory. We introduce our system, prove some desirable features, and show defenses against long range revisions and catastrophic crashes. The Casper overlay provides almost any proof of work chain with additional protections against block reversions.


title: Comprehensive factorisation systems

abstract: We establish a correspondence between consistent comprehension schemes and complete orthogonal factorisation systems. The comprehensive factorisation of a functor between small categories arises in this way. Similar factorisation systems exist for the categories of topological spaces, simplicial sets, small multicategories and Feynman categories. In each case comprehensive factorisation induces a natural notion of universal covering, leading to a Galois-type definition of fundamental group for based objects of the category.


title: Two-Level System Damping in a Quasi-One-Dimensional Optomechanical Resonator

abstract: Nanomechanical resonators have demonstrated great potential for use as versatile tools in a number of emerging quantum technologies. For such applications, the performance of these systems is restricted by the decoherence of their fragile quantum states, necessitating a thorough understanding of their dissipative coupling to the surrounding environment. In bulk amorphous solids, these dissipation channels are dominated at low temperatures by parasitic coupling to intrinsic two-level system (TLS) defects, however, there remains a disconnect between theory and experiment on how this damping manifests in dimensionally-reduced nanomechanical resonators. Here, we present an optomechanically-mediated thermal ringdown technique, which we use to perform simultaneous measurements of the dissipation in four mechanical modes of a cryogenically-cooled silicon nanoresonator, with resonant frequencies ranging from 3 - 19 MHz. Analyzing the device's mechanical damping rate at fridge temperatures between 10 mK - 10 K, we demonstrate quantitative agreement with the standard tunneling model for TLS ensembles confined to one dimension. From these fits, we extract the defect density of states ($P_0 \sim$ 1 - 4 $\times$ 10$^{44}$ J$^{-1}$ m$^{-3}$) and deformation potentials ($\gamma \sim$ 1 - 2 eV), showing that each mechanical mode couples on average to less than a single thermally-active defect at 10 mK.

categories: cond-mat.mes-hall cond-mat.mtrl-sci physics.optics quant-ph


title: Small-mass atomic defects enhance vibrational thermal transport at disordered interfaces with ultrahigh thermal boundary conductance

abstract: The role of interfacial nonidealities and disorder on thermal transport across interfaces is traditionally assumed to add resistance to heat transfer, decreasing the thermal boundary conductance (TBC).$^1$ However, recent computational works have suggested that interfacial defects can enhance this thermal boundary conductance through emergence of unique vibrations that are intrinsic to the material interface and defect atoms,$^{2-6}$ a finding that contradicts traditional theory and conventional understanding. By manipulating the local heat flux of atomic vibrations that comprise these interfacial modes, in principle, the TBC can be increased. In this work, we provide evidence that interfacial defects can enhance the TBC across interfaces through the emergence of unique high frequency vibrational modes that arise from atomic mass defects at the interface with relatively small masses. We demonstrate ultrahigh TBC at amorphous SiOC:H/SiC:H interfaces, approaching 1 GW m$^{-2}$ K$^{-1}$, that is further increased through the introduction of nitrogen defects. The fact that disordered interfaces can exhibit such high conductances, which can be further increased with additional defects offers a unique direction in controlling interfacial thermal transport that becomes important in manipulating heat transfer across materials with high densities of interfaces.


title: High Five: Improving Gesture Recognition by Embracing Uncertainty

abstract: Sensors on mobile devices---accelerometers, gyroscopes, pressure meters, and GPS---invite new applications in gesture recognition, gaming, and fitness tracking. However, programming them remains challenging because human gestures captured by sensors are noisy. This paper illustrates that noisy gestures degrade training and classification accuracy for gesture recognition in state-of-the-art deterministic Hidden Markov Models (HMM). We introduce a new statistical quantization approach that mitigates these problems by (1) during training, producing gesture-specific codebooks, HMMs, and error models for gesture sequences; and (2) during classification, exploiting the error model to explore multiple feasible HMM state sequences. We implement classification in Uncertain<t>, a probabilistic programming system that encapsulates HMMs and error models and then automates sampling and inference in the runtime. Uncertain<T> developers directly express a choice of application-specific trade-off between recall and precision at gesture recognition time, rather than at training time. We demonstrate benefits in configurability, precision, recall, and recognition on two data sets with 25 gestures from 28 people and 4200 total gestures. Incorporating gesture error more accurately in modeling improves the average recognition rate of 20 gestures from 34\% in prior work to 62\%. Incorporating the error model during classification further improves the average gesture recognition rate to 71\%. As far as we are aware, no prior work shows how to generate an HMM error model during training and use it to improve classification rates.


title: Semileptonic $B$-meson decays to light pseudoscalar mesons on the HISQ ensembles

abstract: We report the status of an ongoing lattice-QCD calculation of form factors for exclusive semileptonic decays of $B$ mesons with both charged currents ($B\to\pi\ell\nu$, $B_s\to K\ell\nu$) and neutral currents ($B\to\pi\ell^+\ell^-$, $B\to K\ell^+\ell^-$). The results are important for constraining or revealing physics beyond the Standard Model. This work uses MILC's (2+1+1)-flavor ensembles with the HISQ action for the sea and light valence quarks and the clover action in the Fermilab interpretation for the $b$ quark. Simulations are carried out at three lattice spacings down to $0.088$ fm, with both physical and unphysical sea-quark masses. We present preliminary results for correlation-function fits.


title: Bayesian Inference over the Stiefel Manifold via the Givens Representation

abstract: We introduce an approach based on the Givens representation for posterior inference in statistical models with orthogonal matrix parameters, such as factor models and probabilistic principal component analysis (PPCA). We show how the Givens representation can be used to develop practical methods for transforming densities over the Stiefel manifold into densities over subsets of Euclidean space. We show how to deal with issues arising from the topology of the Stiefel manifold and how to inexpensively compute the change-of-measure terms. We introduce an auxiliary parameter approach that limits the impact of topological issues. We provide both analysis of our methods and numerical examples demonstrating the effectiveness of the approach. We also discuss how our Givens representation can be used to define general classes of distributions over the space of orthogonal matrices. We then give demonstrations on several examples showing how the Givens approach performs in practice in comparison with other methods.


title: Fluctuation relation for heat exchange in Markovian open quantum systems

abstract: A fluctuation relation for the heat exchange of an open quantum system under a thermalizing Markovian dynamics is derived. We show that the probability of that the system absorbs an amount of heat from its bath, at a given time interval, divided by the probability of the reverse process (releasing the same amount of heat to the bath) is given by an exponential factor which depends on the amount of heat and the difference between the temperatures of the system and the bath. We also argue that the probability of the violation of the second law of thermodynamics (here in the form of net heat transfer from a cold system to its hot bath) drops exponentially with both the amount of heat and the temperature differences.


title: Spin Filtering with Poly-T Wrapped Single Wall Carbon Nanotubes

abstract: Spin filtering is an essential operation in spintronics that allows creation and detection of spin polarized carriers. Transition metal ferromagnets are used as spin filters in most cases, though their spin filtering efficiency is only around ~50%, thereby limiting the efficiency of spintronic devices. Recently, chiral systems such as DNA have been shown to exhibit efficient spin filtering, a phenomenon often dubbed as "chirality induced spin selectivity" (CISS). In this work, we consider single wall carbon nanotubes helically wrapped with single stranded poly-T DNA. By magnetoresistance measurements we show that this system exhibits significant spin polarization of ~80%, which could be attributed to the Rashba spin-orbit interaction induced by the inversion asymmetric helical potential of the DNA. Observed spin polarization is larger than that reported before for d(GT)15 strands. Such systems allow tailoring spin polarization by chemical means and also allow extremely localized creation and detection of spin polarization without any magnetic element and could lead to extreme miniaturization and compact integration of spintronic devices and circuits.


title: Scalar-tensor gravity in the Palatini approach

abstract: The main objective of this thesis is to discuss scalar-tensor theories in the Palatini approach. Both scalar-tensor theories and Palatini formalism are means of alternating classical theory of gravity, general relativity, in order to account for phenomena being seemingly unexplainable on the ground of the Einstein theory or to serve as toy models used to test limitations of the theory in question. In the literature, both Palatini approach and scalar-tensor theories have been widely discussed, but there are very few - if none - authors writing about a merge of these two ideas. The present paper is a result of an insufficient attention given to the topic of scalar-tensor theories in Palatini formalism. In the course of the thesis action functional for scalar-tensor theories of gravity will be introduced. This action functional differs significantly from the action defined in case of scalar-tensor theories in metric approach. We aim at analysing the theory using the language of invariants, allowing us to write down all equations in a frame-independent manner. We discover that invariants defined for the metric case do not always have their counterparts in Palatini formalism. Also, two frames most frequently used in the literature are discussed: Einstein and Jordan frame. Possible applications of the theory developed in the first part of the thesis are presented. We show the equivalence between $f(R)$ and scalar-tensor theories of gravity and exploit this fact by analysing the former using methods developed for scalar-tensor theories. We conclude the thesis with calculating the Friedmann equations for an empty universe of vanishing spatial curvature and preparing set-up for analysing inflationary behaviour.


title: Stochastic Non-convex Optimization with Strong High Probability Second-order Convergence

abstract: In this paper, we study stochastic non-convex optimization with non-convex random functions. Recent studies on non-convex optimization revolve around establishing second-order convergence, i.e., converging to a nearly second-order optimal stationary points. However, existing results on stochastic non-convex optimization are limited, especially with a high probability second-order convergence. We propose a novel updating step (named NCG-S) by leveraging a stochastic gradient and a noisy negative curvature of a stochastic Hessian, where the stochastic gradient and Hessian are based on a proper mini-batch of random functions. Building on this step, we develop two algorithms and establish their high probability second-order convergence. To the best of our knowledge, the proposed stochastic algorithms are the first with a second-order convergence in {\it high probability} and a time complexity that is {\it almost linear} in the problem's dimensionality.


title: Disordered hyperuniformity in two-component non-additive hard disk plasmas

abstract: We study the behavior of a two-component plasma made up of non-additive hard disks with a logarithmic Coulomb interaction. Due to the Coulomb repulsion, long-wavelength total density fluctuations are suppressed and the system is globally hyperuniform. Short-range volume effects lead to phase separation or to hetero-coordination for positive or negative non-additivities, respectively. These effects compete with the hidden long-range order imposed by hyperuniformity. As a result, the critical behavior of the mixture is modified, with long-wavelength concentration fluctuations partially damped when the system is charged. It is also shown that the decrease of configurational entropy due to hyperuniformity originates from contributions beyond the two-particle level. Finally, despite global hyperuniformity, we show that in our system, the spatial configuration associated with each component separately is not hyperuniform, i.e., the system is not "multihyperuniform."


title: Millimeter Wave MIMO Prototype: Measurements and Experimental Results

abstract: Millimeter-wave multi-input multi-output (mm-Wave MIMO) systems are one of the candidate schemes for 5G wireless standardization efforts. In this context, the main contributions of this article are three-fold. 1) We describe parallel sets of measurements at identical transmit-receive location pairs with 2.9, 29 and 61 GHz carrier frequencies in indoor office, shopping mall, and outdoor settings. These measurements provide insights on propagation, blockage and material penetration losses, and the key elements necessary in system design to make mm-Wave systems viable in practice. 2) One of these elements is hybrid beamforming necessary for better link margins by reaping the array gain with large antenna dimensions. From the class of fully-flexible hybrid beamformers, we describe a robust class of directional beamformers towards meeting the high data-rate requirements of mm-Wave systems. 3) Leveraging these design insights, we then describe an experimental prototype system at 28 GHz that realizes high data-rates on both the downlink and uplink and robustly maintains these rates in outdoor and indoor mobility scenarios. In addition to maintaining large signal constellation sizes in spite of radio frequency challenges, this prototype leverages the directional nature of the mm-Wave channel to perform seamless beam switching and handover across mm-Wave base-stations thereby overcoming the path losses in non-line-of-sight links and blockages encountered at mm-Wave frequencies.


title: Does the Compact Radio Jet in PG 1700+518 Drive a Molecular Outflow?

abstract: Radio jets play an important role in quasar feedback, but direct observations showing how the jets interact with the multi-phase interstellar medium of galaxy disks are few and far between. In this work, we provide new millimeter interferometric observations of PG 1700+518 in order to investigate the effect of its radio jet on the surrounding molecular gas. PG 1700 is a radio-quiet, low-ionization broad absorption line quasar whose host galaxy has a nearby interacting companion. On sub-kiloparsec scales, the ionized gas is driven to high velocities by a compact radio jet that is identified by radio interferometry. We present observations from the NOrthern Extended Millimeter Array (NOEMA) interferometer with a 3.8 arcsec (16 kpc) synthesized beam where we detect the CO (1-0) emission line at $30\sigma$ significance with a total flux of $3.12\pm0.02$ Jy km s$^{-1}$ and a typical velocity dispersion of $125\pm5$ km s$^{-1}$. Despite the outflow in ionized gas, we find no concrete evidence that the CO gas is being affected by the radio jet on size scales of a kiloparsec or more. However, a $\sim\!1$ arcsec drift in the spatial centroid of the CO emission as a function of velocity across the emission line and the compact nature of the jet hint that higher spatial resolution observations may reveal a signal of interaction between the jet and molecular gas.


title: Spatial Field estimation from Samples taken at Unknown Locations generated by an Unknown Autoregressive Process

abstract: Sampling of physical fields with mobile sensors is an upcoming field of interest. This offers greater advantages in terms of cost as often just a single sensor can be used for the purpose and this can be employed almost everywhere without sensing stations and has nominal operational costs. In a mobile sensing setup, the accurate knowledge of sampling locations may lead to a manifold increase in the costs. Moreover, the inertia of the moving vehicle constrains the independence between the intersample distances, making them correlated. This work, thus, aims at estimating spatially bandlimited fields from samples, corrupted with measurement noise, collected on sampling locations obtained from an autoregressive model on the intersample distances. The autoregressive model is used to capture the correlation between the intersample distances. In this setup of sampling at unknown sampling locations obtained from an autoregressive model on intersample distances, the mean squared error between the field and its estimated version has been shown to decrease as O($1/n$), where $n$ is the average number of samples taken by the moving sensor.


title: Efficient analysis of stochastic gene dynamics in the non-adiabatic regime using piecewise deterministic Markov processes

abstract: Single-cell experiments show that gene expression is stochastic and bursty, a feature that can emerge from slow switching between promoter states with different activities. One source of long-lived promoter states is the slow binding and unbinding kinetics of transcription factors to promoters, i.e. the non-adiabatic binding regime. Here, we introduce a simple analytical framework, known as a piecewise deterministic Markov process (PDMP), that accurately describes the stochastic dynamics of gene expression in the non-adiabatic regime. We illustrate the utility of the PDMP on a non-trivial dynamical system by analyzing the properties of a titration-based oscillator in the non-adiabatic limit. We first show how to transform the underlying Chemical Master Equation into a PDMP where the slow transitions between promoter states are stochastic, but whose rates depend upon the faster deterministic dynamics of the transcription factors regulated by these promoters. We show that the PDMP accurately describes the observed periods of stochastic cycles in activator and repressor-based titration oscillators. We then generalize our PDMP analysis to more complicated versions of titration-based oscillators to explain how multiple binding sites lengthen the period and improve coherence. Last, we show how noise-induced oscillation previously observed in a titration-based oscillator arises from non-adiabatic and discrete binding events at the promoter site.


title: On the interpolation space $(L^p(\Omega), W^{1,p}(\Omega))_{s,p}$ in non-smooth domains

abstract: We show that, for certain non-smooth bounded domains $\Omega\subset\mathbb{R}^n$, the real interpolation space $(L^p(\Omega), W^{1,p}(\Omega))_{s,p}$ is the subspace $\widetilde W^{s,p}(\Omega) \subset L^p(\Omega)$ induced by the restricted fractional seminorm $$ |f|_{\widetilde W^{s,p}(\Omega)} = \Big( \int_\Omega \int_{|x-y|<\frac{d(x)}2} \frac{|f(x)-f(y)|^p}{|x-y|^{n+sp}} \, dy \,dx \Big)^\frac{1}{p}. $$ In particular, the above result includes simply connected uniform domains in the plane, for which a characterization of the interpolation space was previously unknown.


title: On Bandlimited Spatiotemporal Field Sampling with Location and Time Unaware Mobile Sensors

abstract: Sampling of a spatiotemporal field for environmental sensing is of interest. Traditionally, a few fixed stations or sampling locations aid in the reconstruction of the spatial field. Recently, there has been an interest in mobile sensing and location-unaware sensing of spatial fields. In this work, the general class of fields governed by a constant coefficient linear partial differential equations is considered. This class models many physical fields including temperature, pollution, and diffusion fields. The analysis is presented in one dimension for a first exposition. It is assumed that a mobile sensing unit is available to sample the spatiotemporal field of interest at unknown locations and unknown times -- both generated by independent and unknown renewal processes. Based on measurement-noise affected samples, a spatial field reconstruction algorithm is proposed. It is shown that if the initial condition on the field is bandlimited, then the mean squared error between the original and the estimated field decreases as $O(1/n)$, where $n$ is the average sampling density of the mobile sensor. The field reconstruction algorithm is universal, that is, it does not require the knowledge of unknown sampling locations' or unknown time instants' statistical distributions.


title: Ninth Planet or Wandering Star ?

abstract: We study the gravitational effects of two celestrial bodies on a typical object of the Kuyper Belt. The first body is a kuyperian object itself with fairly large eccentricity and perihelion but with a large mass, about 16 times the mass of the Earth. The second body is a star whose mass is 30% - 50% of the mass of the sun that passes by our solar system at a speed between 25 km/sec and 100 km/sec and at a distance of closest approach between 0.05 and 0.5 light year. As a measure of the perturbations caused by these bodies on the light kuyperian object, we analyse its eccentricity. We find that the effects due to the passage of the wandering star are permanent in the sense that the eccentricity of the kuyperian object remains anomalous long after the passage of the star. The same is true of the heavy kuyperian object: it can perturb greatly the orbit of the lighter kuyperian object, which leads to a permanent, anomalous eccentricity.


title: Straightforward measurement of thermal properties of anisotropic materials: the case of Bi2Se3 single crystal

abstract: We show that by a simple two-steps measurement protocol one can extract the thermal properties of anisotropic insulator materials. We focus on Bi2Se3 single crystal, a layered crystal with individual sheets held together by Van der Waals force, presenting strong heat diffusion anisotropy. Since the shape of the sample does not allow experiments in the perpendicular direction the challenge is to measure the in plane and out of plane thermal properties by in plane measurements. The optical properties of Bi2Se3 being known (opaque sample for the pump and probe beams), thermoreflectance microscopy experiments can be done on the bare substrate. In this case, the lateral heat diffusion is governed by the in plane thermal diffusivity. A 100 nm gold layer is then deposited on the sample. In this second experiment it is straightforward to demonstrate that the perpendicular thermal properties of the layer rule the lateral diffusivity. The good agreement of theoretical curves obtained by thermal models with experimental data permit us to evaluate the thermal conductivity coefficients of bulk Bi2Se3.


title: The connection between the peaks in velocity dispersion and star-forming clumps of turbulent galaxies

abstract: We present Keck/OSIRIS adaptive optics observations with 150-400 pc spatial sampling of 7 turbulent, clumpy disc galaxies from the DYNAMO sample ($0.07<z<0.2$). DYNAMO galaxies have previously been shown to be well matched in properties to main sequence galaxies at $z\sim1.5$. Integral field spectroscopy observations using adaptive optics are subject to a number of systematics including a variable PSF and spatial sampling, which we account for in our analysis. We present gas velocity dispersion maps corrected for these effects, and confirm that DYNAMO galaxies do have high gas velocity dispersion ($\sigma=40-80$\kms), even at high spatial sampling. We find statistically significant structure in 6 out of 7 galaxies. The most common distance between the peaks in velocity dispersion and emission line peaks is $\sim0.5$~kpc, we note this is very similar to the average size of a clump measured with HST H$\alpha$ maps. This could suggest that the peaks in velocity dispersion in clumpy galaxies likely arise due to some interaction between the clump and the surrounding ISM of the galaxy, though our observations cannot distinguish between outflows, inflows or velocity shear. Observations covering a wider area of the galaxies will be needed to confirm this result.


title: Positrons from Primordial Black Hole Microquasars and Gamma-ray Bursts

abstract: We propose several novel scenarios how capture of small sublunar-mass primordial black holes (PBHs) by compact stars, white dwarfs or neutron stars, can lead to distinct short gamma-ray bursts (sGRBs) as well as microquasars (MQs). In addition to providing new signatures, relativistic jets from these systems will accelerate positrons to high energies. We find that if PBHs constitute a sizable fraction of DM, they can significantly contribute to the excess observed in the positron flux by the Pamela, the AMS-02 and the Fermi-LAT experiments. Our proposal combines the beneficial features of astrophysical sources and dark matter.


title: Chirality provides a direct fitness advantage and facilitates intermixing in cellular aggregates

abstract: Chirality in shape and motility can evolve rapidly in microbes and cancer cells. To determine how chirality affects cell fitness, we developed a model of chiral growth in compact aggregates such as microbial colonies and solid tumors. Our model recapitulates previous experimental findings and shows that mutant cells can invade by increasing their chirality or switching their handedness. The invasion results either in a takeover or stable coexistence between the mutant and the ancestor depending on their relative chirality. For large chiralities, the coexistence is accompanied by strong intermixing between the cells, while spatial segregation occurs otherwise. We show that the competition within the aggregate is mediated by bulges in regions where the cells with different chiralities meet. The two-way coupling between aggregate shape and natural selection is described by the chiral Kardar-Parisi-Zhang equation coupled to the Burgers' equation with multiplicative noise. We solve for the key features of this theory to explain the origin of selection on chirality. Overall, our work suggests that chirality could be an important ecological trait that mediates competition, invasion, and spatial structure in cellular populations.

categories: q-bio.PE cond-mat.soft cond-mat.stat-mech nlin.PS physics.bio-ph


title: Phase-field model of vapor-liquid-solid nanowire growth

abstract: We present a multi-phase-field model to describe quantitatively nanowire growth by the vapor-liquid-solid (VLS) process. The free-energy functional of this model depends on three non-conserved order parameters that distinguish the vapor, liquid, and solid phases. The evolution equations for those order parameters describe basic kinetic processes including the rapid (quasi-instantaneous) equilibration of the liquid catalyst to a droplet shape with constant mean curvature, the slow incorporation of growth atoms at the droplet surface, and crystallization within the droplet. The standard constraint that the sum of the phase fields equals unity and the conservation of the number of catalyst atoms, which relates the catalyst volume to the concentration of growth atoms inside the droplet, are handled via separate Lagrange multipliers. An analysis of the model is presented that rigorously maps the phase-field equations to a desired set of sharp-interface equations for the evolution of the phase boundaries under the constraint of force-balance at three-phase junctions given by the Young-Herring relation that includes torque term related to the anisotropy of the solid-liquid and solid-vapor interface excess free energies. Numerical examples of growth in two dimensions are presented. The simulations reproduce many of the salient features of nanowire growth observed experimentally, including growth normal to the substrate with tapering of the side walls, transitions between different growth orientations, and crawling growth along the substrate.


title: On Comparison Of Experts

abstract: A policy maker faces a sequence of unknown outcomes. At each stage two (self-proclaimed) experts provide probabilistic forecasts on the outcome in the next stage. A comparison test is a protocol for the policy maker to (eventually) decide which of the two experts is better informed. The protocol takes as input the sequence of pairs of forecasts and actual realizations and (weakly) ranks the two experts. We propose two natural properties that such a comparison test must adhere to and show that these essentially uniquely determine the comparison test. This test is a function of the derivative of the induced pair of measures at the realization.


title: Mixed-strategy Nash equilibrium for a discontinuous symmetric $N$-player game

abstract: We consider a game in which each player must find a compromise between more daring strategies that carry a high risk for him to be eliminated, and more cautious ones that, however, reduce his final score. For two symmetric players this game was originally formulated in 1961 by Dresher, who modeled a duel between two opponents. The game has also been of interest in the description of athletic competitions. We extend here the two-player game to an arbitrary number $N$ of symmetric players. We show that there is a mixed-strategy Nash equilibrium and find its exact analytic expression, which we analyze in particular in the limit of large $N$, where mean-field behavior occurs. The original game with $N=2$ arises as a singular limit of the general case.


title: Theoretical Investigation of Li and Na Oxides Adsorption on TiC(111) Surface for Metal-Air Rechargeable Batteries

abstract: We analyze, with Density Functional Theory (DFT) calculations, the adsorption energies of Li$_2$O$_2$, Na$_2$O$_2$ and NaO$_2$ on clean and oxygen passivated TiC (111) surfaces. We show, that after deposition of two molecular layers of alkali metal oxides, the initial state of the TiC surface becomes unimportant for the adsorption energy and that all adsorption energies approach their native crystal values. The structure of the adsorbed molecular layers is analyzed and compared to their native oxide crystal structure. Finally, we discuss the possible implications for electrode optimization for Li-air and Na-air batteries.


title: Fast and Luminous Transients from the Explosions of Long Lived Massive White Dwarf Merger Remnants

abstract: We study the evolution and final outcome of long-lived (${\approx}10^5$ years) remnants from the merger of a He white dwarf (WD) with a more massive C/O or O/Ne WD. Using Modules for Experiments in Stellar Astrophysics ($\texttt{MESA}$), we show that these remnants have a red giant configuration supported by steady helium burning, adding mass to the WD core until it reaches $M_{\rm core}\approx 1.12-1.20 M_\odot$. At that point, the base of the surface convection zone extends into the burning layer, mixing the helium burning products (primarily carbon and magnesium) throughout the convective envelope. Further evolution depletes the convective envelope of helium, and dramatically slows the mass increase of the underlying WD core. The WD core mass growth re-initiates after helium depletion, as then an uncoupled carbon burning shell is ignited and proceeds to burn the fuel from the remaining metal-rich extended envelope. For large enough initial total merger masses, O/Ne WD cores would experience electron-capture triggered collapse to neutron stars (NSs) after growing to near Chandrasekhar mass ($M_{\rm Ch}$). Massive C/O WD cores could suffer the same fate after a carbon-burning flame converts them to O/Ne. The NS formation would release ${\approx}10^{50}$ ergs into the remaining extended low mass envelope. Using the STELLA radiative transfer code, we predict the resulting optical light curves from these exploded envelopes. Reaching absolute magnitudes of $M_V\approx -17$, these transients are bright for about one week, and have many features of the class of luminous, rapidly evolving transients studied by Drout and collaborators.


title: Parameter constraints from cross-correlation of CMB lensing with galaxy clustering

abstract: The lensing convergence measurable with future CMB surveys like CMB-S4 will be highly correlated with the clustering observed by deep photometric large scale structure (LSS) surveys such as the LSST, with cross-correlation coefficient as high as 95\%. This will enable use of sample variance cancellation techniques to determine cosmological parameters, and use of cross-correlation measurements to break parameter degeneracies. Assuming large sky overlap between CMB-S4 and LSST, we show that a joint analysis of CMB-S4 lensing and LSST clustering can yield very tight constraints on the matter amplitude $\sigma_8(z)$, halo bias, and $f_\mathrm{NL}$, competitive with the best stage IV experiment predictions, but using complementary methods, which may carry different and possibly lower systematics. Having no sky overlap between experiments degrades the precision of $\sigma_8(z)$ by a factor of 20, and that of $f_\mathrm{NL}$ by a factor of 1.5 to 2. Without CMB lensing, the precision always degrades by an order of magnitude or more, showing that a joint analysis is critical. Our results also suggest that CMB lensing in combination with LSS photometric surveys is a competitive probe of the evolution of structure in the redshift range $z\simeq 1-7$, probing a regime that is not well tested observationally. We explore predictions against other surveys and experiment configurations, finding that wide patches with maximal sky overlap between CMB and LSS surveys are most powerful for $\sigma_8(z)$ and $f_\mathrm{NL}$.


title: Optimal Mechanism Design with Flexible Consumers and Costly Supply

abstract: The problem of designing a profit-maximizing, Bayesian incentive compatible and individually rational mechanism with flexible consumers and costly heterogeneous supply is considered. In our setup, each consumer is associated with a flexibility set that describes the subset of goods the consumer is equally interested in. Each consumer wants to consume one good from its flexibility set. The flexibility set of a consumer and the utility it gets from consuming a good from its flexibility set are its private information. We adopt the flexibility model of [1] and focus on the case of nested flexibility sets -- each consumer's flexibility set can be one of k nested sets. Examples of settings with this inherent nested structure are provided. On the supply side, we assume that the seller has an initial stock of free supply but it can purchase more goods for each of the nested sets at fixed exogenous prices. We characterize the allocation and purchase rules for a profit-maximizing, Bayesian incentive compatible and individually rational mechanism as the solution to an integer program. The optimal payment function is pinned down by the optimal allocation rule in the form of an integral equation. We show that the nestedness of flexibility sets can be exploited to obtain a simple description of the optimal allocations, purchases and payments in terms of thresholds that can be computed through a straightforward iterative procedure.


title: Comparing Dushnik-Miller Dimension, Boolean Dimension and Local Dimension

abstract: The original notion of dimension for posets is due to Dushnik and Miller and has been studied extensively in the literature. Quite recently, there has been considerable interest in two variations of dimension known as Boolean dimension and local dimension. For a poset $P$, the Boolean dimension of $P$ and the local dimension of $P$ are both bounded from above by the dimension of $P$ and can be considerably less. Our primary goal will be to study analogies and contrasts among these three parameters. As one example, it is known that the dimension of a poset is bounded as a function of its height and the tree-width of its cover graph. The Boolean dimension of a poset is bounded in terms of the tree-width of its cover graph, independent of its height. We show that the local dimension of a poset cannot be bounded in terms of the tree-width of its cover graph, independent of height. We also prove that the local dimension of a poset is bounded in terms of the path-width of its cover graph. In several of our results, Ramsey theoretic methods will be applied.


title: Infrared Flares from M Dwarfs: a Hinderance to Future Transiting Exoplanet Studies

abstract: Many current and future exoplanet missions are pushing to infrared (IR) wavelengths where the flux contrast between the planet and star is more favorable (Deming et al. 2009), and the impact of stellar magnetic activity is decreased. Indeed, a recent analysis of starspots and faculae found these forms of stellar activity do not substantially impact the transit signatures or science potential for FGKM stars with JWST (Zellem et al. 2017). However, this is not true in the case of flares, which I demonstrate can be a hinderance to transit studies in this note.


title: Logical relations for coherence of effect subtyping

abstract: A coercion semantics of a programming language with subtyping is typically defined on typing derivations rather than on typing judgments. To avoid semantic ambiguity, such a semantics is expected to be coherent, i.e., independent of the typing derivation for a given typing judgment. In this article we present heterogeneous, biorthogonal, step-indexed logical relations for establishing the coherence of coercion semantics of programming languages with subtyping. To illustrate the effectiveness of the proof method, we develop a proof of coherence of a type-directed, selective CPS translation from a typed call-by-value lambda calculus with delimited continuations and control-effect subtyping. The article is accompanied by a Coq formalization that relies on a novel shallow embedding of a logic for reasoning about step-indexing.


title: An inexact Newton-Krylov method for stochastic eigenvalue problems

abstract: This paper aims at the efficient numerical solution of stochastic eigenvalue problems. Such problems often lead to prohibitively high dimensional systems with tensor product structure when discretized with the stochastic Galerkin method. Here, we exploit this inherent tensor product structure to develop a globalized low-rank inexact Newton method with which we tackle the stochastic eigenproblem. We illustrate the effectiveness of our solver with numerical experiments.


title: Inductive Representation Learning in Large Attributed Graphs

abstract: Graphs (networks) are ubiquitous and allow us to model entities (nodes) and the dependencies (edges) between them. Learning a useful feature representation from graph data lies at the heart and success of many machine learning tasks such as classification, anomaly detection, link prediction, among many others. Many existing techniques use random walks as a basis for learning features or estimating the parameters of a graph model for a downstream prediction task. Examples include recent node embedding methods such as DeepWalk, node2vec, as well as graph-based deep learning algorithms. However, the simple random walk used by these methods is fundamentally tied to the identity of the node. This has three main disadvantages. First, these approaches are inherently transductive and do not generalize to unseen nodes and other graphs. Second, they are not space-efficient as a feature vector is learned for each node which is impractical for large graphs. Third, most of these approaches lack support for attributed graphs. To make these methods more generally applicable, we propose a framework for inductive network representation learning based on the notion of attributed random walk that is not tied to node identity and is instead based on learning a function $\Phi : \mathrm{\rm \bf x} \rightarrow w$ that maps a node attribute vector $\mathrm{\rm \bf x}$ to a type $w$. This framework serves as a basis for generalizing existing methods such as DeepWalk, node2vec, and many other previous methods that leverage traditional random walks.


title: Lattice Boltzmann method for semiclassical fluids

abstract: We determine properties of the lattice Boltzmann method for semiclassical fluids, which is based on the Boltzmann equation and the equilibrium distribution function is given either by the Bose-Einstein or the Fermi-Dirac ones. New D-dimensional polynomials, that generalize the Hermite ones, are introduced and we find that the weight that renders the polynomials orthonormal has to be approximately equal, or equal, to the equilibrium distribution function itself for an efficient numerical implementation of the lattice Boltzmann method. In light of the new polynomials we discuss the convergence of the series expansion of the equilibrium distribution function and the obtainment of the hydrodynamic equations. A discrete quadrature is proposed and some discrete lattices in one, two and three dimensions associated to weight functions other than the Hermite weight are obtained. We derive the forcing term for the LBM, given by the Lorentz force, which dependents on the microscopic velocity, since the bosonic and fermionic particles can be charged. Motivated by the recent experimental observations of the hydrodynamic regime of electrons in graphene, we build an isothermal lattice Boltzmann method for electrons in metals in two and three dimensions. This model is validated by means of the Riemann problem and of the Poiseuille flow. As expected for electron in metals, the Ohm's law is recovered for a system analogous to a porous medium.


title: Spin-1 two-impurity Kondo problem on a lattice

abstract: We present an extensive study of the two-impurity Kondo problem for spin-1 adatoms on square lattice using an exact canonical transformation to map the problem onto an effective one-dimensional system that can be numerically solved using the density matrix renormalization group method. We provide a simple intuitive picture and identify the different regimes, depending on the distance between the two impurities, Kondo coupling $J_K$, longitudinal anisotropy $D$, and transverse anisotropy $E$. In the isotropic case, two impurities on opposite(same) sublattices have a singlet(triplet) ground state. However, the energy difference between the triplet ground state and the singlet excited state is very small and we expect an effectively four-fold degenerate ground state, i.e., two decoupled impurities. For large enough $J_K$ the impurities are practically uncorrelated forming two independent underscreened states with the conduction electrons, a clear non-perturbative effect. When the impurities are entangled in an RKKY-like state, Kondo correlations persists and the two effects coexist: the impurities are underscreened, and the dangling spin-$1/2$ degrees of freedom are responsible for the inter-impurity entanglement. We analyze the effects of magnetic anisotropy in the development of quasi-classical correlations.


title: Topological charge and cooling scales in pure SU(2) lattice gauge theory

abstract: Using Monte Carlo simulations with overrelaxation, we have equilibrated lattices up to $\beta=2.928$, size $60^4$, for pure SU(2) lattice gauge theory with the Wilson action. We calculate topological charges with the standard cooling method and find that they become more reliable with increasing $\beta$ values and lattice sizes. Continuum limit estimates of the topological susceptibility $\chi$ are obtained of which we favor $\chi^{1/4}/T_c=0.643\,(12)$, where $T_c$ is the SU(2) deconfinement temperature. Differences between cooling length scales in different topological sectors turn out to be too small to be detectable within our statistical errors.


title: Triangular fractal approximating graphs and their covering paths and cycles

abstract: We observe two kinds of fractal approximating graphs, the background structures of the generalized Sierpinski Arrowhead Curve independently of the recursive curves. Both graphs related to the generalized Sierpinski Gasket and based on a checked triangular generator pattern. In the Overall Graph we connect the corners of the up facing neighbouring dark tiles. In the Inscribed Graph we connect their centroids. We describe their cardinalities in general case with recursive and explicit formulas, the numbers of their vertices and edges, their edge covering Hamiltonian-paths and -cycles, and their tiling-paths and -cycles which cover all of their dark tiles. Some of these formations are unambiguously transformable into each other and these bijective pairs keep their basic properties after the transformation. Some of their cardinalities form new integer sequences.


title: A Mathematical Analysis of Technical Analysis

abstract: In this paper, we investigate trading strategies based on exponential moving averages (ExpMAs) of an underlying risky asset. We study both logarithmic utility maximization and long-term growth rate maximization problems and find closed-form solutions when the drift of the underlying is modeled by either an Ornstein-Uhlenbeck process or a two-state continuous-time Markov chain. For the case of an Ornstein-Uhlenbeck drift, we carry out several Monte Carlo experiments in order to investigate how the performance of optimal ExpMA strategies is affected by variations in model parameters and by transaction costs.


title: Fair division with multiple pieces

abstract: Given a set of $p$ players we consider problems concerning envy-free allocation of collections of $k$ pieces from a given set of goods or chores. We show that if $p\le n$ and each player can choose $k$ pieces out of $n$ pieces of a cake, then there exist a division of the cake and an allocation of the pieces where at least $\frac{p}{2(k^2-k+1)}$ players get their desired $k$ pieces each. We further show that if $p\le k(n-1)+1$ and each player can choose $k$ pieces, one from each of $k$ cakes that are divided into $n$ pieces each, then there exist a division of the cakes and allocation of the pieces where at least $\frac{p}{2k(k-1)}$ players get their desired $k$ pieces. Finally we prove that if $p\ge k(n-1)+1$ and each player can choose one shift in each of $k$ days that are partitioned into $n$ shifts each, then, given that the salaries of the players are fixed, there exist $n(1+\ln k)$ players covering all the shifts, and moreover, if $k=2$ then $n$ players suffice. Our proofs combine topological methods and theorems of F\"uredi, Lov\'asz and Gallai from hypergraph theory.


title: Highlights from the Pierre Auger Observatory (ICRC17)

abstract: In this contribution we summarize the highlights from the Pierre Auger Observatory presented at the 35th International Cosmic Ray Conference. We discuss the update of the measurement of the energy spectrum of cosmic rays over a wide range of energy ($10^{17.5}$ to above $10^{20}$ eV), studies of the cosmic-ray mass composition with the fluorescence and surface detector of the Observatory, the discovery of a large-scale anisotropy in the arrival direction of cosmic rays above $8\times 10^{18}$ eV and indications of anisotropy at intermediate angular scales above $4 \times 10^{19}$ eV. Moreover, we report on tests of hadronic interactions beyond LHC energies, multi-messenger analyses with neutral primaries and the progress of the upgrade of the Observatory, AugerPrime, aimed at elucidating the origin of the observed flux suppression at ultra-high energies.


title: K\"ahler Uniformization from Holographic Renormalization Group Flows of M5-branes

abstract: In this paper, we initiate the study of holographic renormalization group flows acting on the metric of four-manifolds. In particular, we derive a set of equations which govern the evolution of a generic K\"ahler four-manifold along the renormalization group flow in seven-dimensional gauged supergravity. The physical eleven-dimensional M-theory setup is given by a stack of M5-branes wrapping a calibrated K\"ahler four-cycle inside a Calabi-Yau threefold. By topologically twisting the theory in the ultraviolet, we may choose an arbitrary K\"ahler metric on the four-cycle as an asymptotic boundary condition. Along the renormalization group flow, the metric moduli are largely washed out, and at the infrared fixed point we will reach a K\"ahler-Einstein metric, which is the expected uniformization behavior.


title: Evolution of jets driven by relativistic radiation hydrodynamics as Long and Low Luminosity GRBs

abstract: We present the three-dimensional numerical simulation of jets modeled with Relativistic Radiation Hydrodynamics (RRH), that evolve across two environments: i) a stratified surrounding medium and ii) a 16TI progenitor model. We consider opacities consistent with various processes of interaction between the fluid and radiation, specifically, free-free, bound-free, bound-bound and electron scattering. We explore various initial conditions, with different radiation energy densities of the beam in hydrodynamical and radiation pressure dominated scenarios, considering only ultra-relativistic jets. In order to investigate the impact of the radiation field on the evolution of the jets, we compare our results with purely hydrodynamical jets. Comparing among jets driven by RRH, we find that radiation pressure dominated jets propagate slightly faster than gas pressure dominated ones. Finally, we construct the luminosity Light Curves (LCs) associated with all these cases. The construction of LCs uses the fluxes of the radiation field which is fully coupled to the hydrodynamics equations during the evolution. The main properties of the jets propagating on the stratified surrounding medium are that the LCs show the same order of magnitude as the gamma-ray luminosity of typical Long Gamma-Ray Bursts $10^{50}-10^{54}$erg/s and the difference between the radiation and gas temperatures is of nearly one order of magnitude. The properties of jets breaking out from the progenitor star model are that the LCs are of the order of magnitude of low-luminosity GRBs $10^{46}-10^{49}$ erg/s, and in this scenario the difference between the gas and radiation temperature is of four orders of magnitude, which is a case far from thermal equilibrium.


title: Additive Matrix Convolutions of P\'olya Ensembles and Polynomial Ensembles

abstract: Recently subclasses of polynomial ensembles for additive and multiplicative matrix convolutions were identified which were called P\'olya ensembles (or polynomial ensembles of derivative type). Those ensembles are closed under the respective convolutions and, thus, build a semi-group when adding by hand a unit element. They even have a semi-group action on the polynomial ensembles. Moreover in several works transformations of the bi-orthogonal functions and kernels of a given polynomial ensemble were derived when performing an additive or multiplicative matrix convolution with particular P\'olya ensembles. For the multiplicative matrix convolution on the complex square matrices the transformations were even done for general P\'olya ensembles. In the present work we generalize these results to the additive convolution on Hermitian matrices, on Hermitian anti-symmetric matrices, on Hermitian anti-self-dual matrices and on rectangular complex matrices. For this purpose we derive the bi-orthogonal functions and the corresponding kernel for a general P\'olya ensemble which was not done before. With the help of these results we find transformation formulas for the convolution with a fixed matrix or a random matrix drawn from a general polynomial ensemble. As an example we consider P\'olya ensembles with an associated weight which is a P\'olya frequency function of infinite order. But we also explicitly evaluate the Gaussian unitary ensemble as well as the complex Laguerre (aka Wishart, Ginibre or chiral Gaussian unitary) ensemble. All results hold for finite matrix dimension. Furthermore we derive a recursive relation between Toeplitz determinants which appears as a by-product of our results.


title: The Sunburst Arc: Direct Lyman {\alpha} escape observed in the brightest known lensed galaxy

abstract: We present rest-frame ultraviolet and optical spectroscopy of the brightest lensed galaxy yet discovered, at redshift z = 2.4. This source reveals a characteristic, triple-peaked Lyman {\alpha} profile which has been predicted by various theoretical works but to our knowledge has not been unambiguously observed previously. The feature is well fit by a superposition of two components: a double-peak profile emerging from substantial radiative transfer, and a narrow, central component resulting from directly escaping Lyman {\alpha} photons; but is poorly fit by either component alone. We demonstrate that the feature is unlikely to contain contamination from nearby sources, and that the central peak is unaffected by radiative transfer effects apart from very slight absorption. The feature is detected at signal-to-noise ratios exceeding 80 per pixel at line center, and bears strong resemblance to synthetic profiles predicted by numerical models.


title: Multimodal Probabilistic Model-Based Planning for Human-Robot Interaction

abstract: This paper presents a method for constructing human-robot interaction policies in settings where multimodality, i.e., the possibility of multiple highly distinct futures, plays a critical role in decision making. We are motivated in this work by the example of traffic weaving, e.g., at highway on-ramps/off-ramps, where entering and exiting cars must swap lanes in a short distance---a challenging negotiation even for experienced drivers due to the inherent multimodal uncertainty of who will pass whom. Our approach is to learn multimodal probability distributions over future human actions from a dataset of human-human exemplars and perform real-time robot policy construction in the resulting environment model through massively parallel sampling of human responses to candidate robot action sequences. Direct learning of these distributions is made possible by recent advances in the theory of conditional variational autoencoders (CVAEs), whereby we learn action distributions simultaneously conditioned on the present interaction history, as well as candidate future robot actions in order to take into account response dynamics. We demonstrate the efficacy of this approach with a human-in-the-loop simulation of a traffic weaving scenario.


title: Precision searches in dijets at the HL-LHC and HE-LHC

abstract: This paper explores the physics reach of the High-Luminosity Large Hadron Collider (HL-LHC) for searches of new particles decaying to two jets. We discuss inclusive searches in dijets and b-jets, as well as searches in semi-inclusive events by requiring an additional lepton that increases sensitivity to different aspects of the underlying processes. We discuss the expected exclusion limits for generic models predicting new massive particles that result in resonant structures in the dijet mass. Prospects of the Higher-Energy LHC (HE-LHC) collider are also discussed. The study is based on the Pythia8 Monte Carlo generator using representative event statistics for the HL-LHC and HE-LHC running conditions. The event samples were created using supercomputers at NERSC.


title: Signed Network Modeling Based on Structural Balance Theory

abstract: The modeling of networks, specifically generative models, have been shown to provide a plethora of information about the underlying network structures, as well as many other benefits behind their construction. Recently there has been a considerable increase in interest for the better understanding and modeling of networks, but the vast majority of this work has been for unsigned networks. However, many networks can have positive and negative links(or signed networks), especially in online social media, and they inherently have properties not found in unsigned networks due to the added complexity. Specifically, the positive to negative link ratio and the distribution of signed triangles in the networks are properties that are unique to signed networks and would need to be explicitly modeled. This is because their underlying dynamics are not random, but controlled by social theories, such as Structural Balance Theory, which loosely states that users in social networks will prefer triadic relations that involve less tension. Therefore, we propose a model based on Structural Balance Theory and the unsigned Transitive Chung-Lu model for the modeling of signed networks. Our model introduces two parameters that are able to help maintain the positive link ratio and proportion of balanced triangles. Empirical experiments on three real-world signed networks demonstrate the importance of designing models specific to signed networks based on social theories to obtain better performance in maintaining signed network properties while generating synthetic networks.


title: A Differential Evaluation Markov Chain Monte Carlo algorithm for Bayesian Model Updating

abstract: The use of the Bayesian tools in system identification and model updating paradigms has been increased in the last ten years. Usually, the Bayesian techniques can be implemented to incorporate the uncertainties associated with measurements as well as the prediction made by the finite element model (FEM) into the FEM updating procedure. In this case, the posterior distribution function describes the uncertainty in the FE model prediction and the experimental data. Due to the complexity of the modeled systems, the analytical solution for the posterior distribution function may not exist. This leads to the use of numerical methods, such as Markov Chain Monte Carlo techniques, to obtain approximate solutions for the posterior distribution function. In this paper, a Differential Evaluation Markov Chain Monte Carlo (DE-MC) method is used to approximate the posterior function and update FEMs. The main idea of the DE-MC approach is to combine the Differential Evolution, which is an effective global optimization algorithm over real parameter space, with Markov Chain Monte Carlo (MCMC) techniques to generate samples from the posterior distribution function. In this paper, the DE-MC method is discussed in detail while the performance and the accuracy of this algorithm are investigated by updating two structural examples.


title: The zeta function of stacks of $G$-zips and truncated Barsotti-Tate groups

abstract: We study stacks of truncated Barsotti-Tate groups and the $G$-zips defined by Pink, Wedhorn & Ziegler. The latter occur naturally when studying truncated Barsotti-Tate groups of height 1 with additional structure. By studying objects over finite fields and their automorphisms we determine the zeta functions of these stacks. These zeta functions can be expressed in terms of the Weyl group of the reductive group $G$ and its action on the root system. The main ingredients are the classification of $G$-zips over algebraically closed fields and their automorphism groups by Pink, Wedhorn & Ziegler, and the study of truncated Barsotti-Tate groups and their automorphism groups by Gabber & Vasiu.


title: Slow and Long-ranged Dynamical Heterogeneities in Dissipative Fluids

abstract: A two-dimensional bidisperse granular fluid is shown to exhibit pronounced long-ranged dynamical heterogeneities as dynamical arrest is approached. Here we focus on the most direct approach to study these heterogeneities: we identify clusters of slow particles and determine their size, $N_c$, and their radius of gyration, $R_G$. We show that $N_c\propto R_G^{d_f}$, providing direct evidence that the most immobile particles arrange in fractal objects with a fractal dimension, $d_f$, that is observed to increase with packing fraction $\phi$. The cluster size distribution obeys scaling, approaching an algebraic decay in the limit of structural arrest, i.e., $\phi\to\phi_c$. Alternatively, dynamical heterogeneities are analyzed via the four-point structure factor $S_4(q,t)$ and the dynamical susceptibility $\chi_4(t)$. $S_4(q,t)$ is shown to obey scaling in the full range of packing fractions, $0.6\leq\phi\leq 0.805$, and to become increasingly long-ranged as $\phi\to\phi_c$. Finite size scaling of $\chi_4(t)$ provides a consistency check for the previously analyzed divergences of $\chi_4(t)\propto (\phi-\phi_c)^{-\gamma_{\chi}}$ and the correlation length $\xi\propto (\phi-\phi_c)^{-\gamma_{\xi}}$. We check the robustness of our results with respect to our definition of mobility. The divergences and the scaling for $\phi\to\phi_c$ suggest a non-equilibrium glass transition which seems qualitatively independent of the coefficient of restitution.


title: Scalable Neural Network Decoders for Higher Dimensional Quantum Codes

abstract: Machine learning has the potential to become an important tool in quantum error correction as it allows the decoder to adapt to the error distribution of a quantum chip. An additional motivation for using neural networks is the fact that they can be evaluated by dedicated hardware which is very fast and consumes little power. Machine learning has been previously applied to decode the surface code. However, these approaches are not scalable as the training has to be redone for every system size which becomes increasingly difficult. In this work the existence of local decoders for higher dimensional codes leads us to use a low-depth convolutional neural network to locally assign a likelihood of error on each qubit. For noiseless syndrome measurements, numerical simulations show that the decoder has a threshold of around $7.1\%$ when applied to the 4D toric code. When the syndrome measurements are noisy, the decoder performs better for larger code sizes when the error probability is low. We also give theoretical and numerical analysis to show how a convolutional neural network is different from the 1-nearest neighbor algorithm, which is a baseline machine learning method.


title: Complete 3D Scene Parsing from an RGBD Image

abstract: One major goal of vision is to infer physical models of objects, surfaces, and their layout from sensors. In this paper, we aim to interpret indoor scenes from one RGBD image. Our representation encodes the layout of orthogonal walls and the extent of objects, modeled with CAD-like 3D shapes. We parse both the visible and occluded portions of the scene and all observable objects, producing a complete 3D parse. Such a scene interpretation is useful for robotics and visual reasoning, but difficult to produce due to the well-known challenge of segmentation, the high degree of occlusion, and the diversity of objects in indoor scenes. We take a data-driven approach, generating sets of potential object regions, matching to regions in training images, and transferring and aligning associated 3D models while encouraging fit to observations and spatial consistency. We use support inference to aid interpretation and propose a retrieval scheme that uses convolutional neural networks (CNNs) to classify regions and retrieve objects with similar shapes. We demonstrate the performance of our method on our newly annotated NYUd v2 dataset with detailed 3D shapes.


title: Electroforming-Free TaOx Memristors using Focused Ion Beam Irradiations

abstract: We demonstrate creation of electroforming-free TaOx memristive devices using focused ion beam irradiations to locally define conductive filaments in TaOx films. Electrical characterization shows that these irradiations directly create fully functional memristors without the need for electroforming. Ion beam forming of conductive filaments combined with state-of-the-art nano-patterning presents a CMOS compatible approach to wafer level fabrication of fully formed and operational memristors.


title: Approximation by mappings with singular Hessian minors

abstract: Let $\Omega\subset\mathbb R^n$ be a Lipschitz domain. Given $1\leq p<k\leq n$ and any $u\in W^{2,p}(\Omega)$ belonging to the little H\"older class $c^{1,\alpha}$, we construct a sequence $u_j$ in the same space with $\operatorname{rank}D^2u_j<k$ almost everywhere such that $u_j\to u$ in $C^{1,\alpha}$ and weakly in $W^{2,p}$. This result is in strong contrast with known regularity behavior of functions in $W^{2,p}$, $p\geq k$, satisfying the same rank inequality.


title: Piercing Numbers in Approval Voting

abstract: We survey a host of results from discrete geometry that have bearing on the analysis of geometric models of approval voting. Such models view the political spectrum as a geometric space, with geometric constraints on voter preferences. Results on piercing numbers then have a natural interpretation in voting theory, and we survey their implications for various classes of geometric constraints on voter approval sets.


title: Runtime Fault Detection in Programmed Molecular Systems

abstract: Watchdog timers are devices that are commonly used to monitor the health of safety-critical hardware and software systems. Their primary function is to raise an alarm if the monitored systems fail to emit periodic "heartbeats" that signal their well-being. In this paper we design and verify a molecular watchdog timer for monitoring the health of programmed molecular nanosystems. This raises new challenges because our molecular watchdog timer and the system that it monitors both operate in the probabilistic environment of chemical kinetics, where many failures are certain to occur and it is especially hard to detect the absence of a signal. Our molecular watchdog timer is the result of an incremental design process that uses goal-oriented requirements engineering, simulation, stochastic analysis, and software verification tools. We demonstrate the molecular watchdog's functionality by having it monitor a molecular oscillator. Both the molecular watchdog timer and the oscillator are implemented as chemical reaction networks, which are the current programming language of choice for many molecular programming applications.


title: A Link Quality Model for Generalised Frequency Division Multiplexing

abstract: 5G systems aim to achieve extremely high data rates, low end-to-end latency and ultra-low power consumption. Recently, there has been considerable interest in the design of 5G physical layer waveforms. One important candidate is Generalised Frequency Division Multiplexing (GFDM). In order to evaluate its performance and features, system-level studies should be undertaken in a range of scenarios. These studies, however, require highly complex computations if they are performed using bit-level simulators. In this paper, the Mutual Information (MI) based link quality model (PHY abstraction), which has been regularly used to implement system-level studies for Orthogonal Frequency Division Multiplexing (OFDM), is applied to GFDM. The performance of the GFDM waveform using this model and the bit-level simulation performance is measured using different channel types. Moreover, a system-level study for a GFDM based LTE-A system in a realistic scenario, using both a bit-level simulator and this abstraction model, has been studied and compared. The results reveal the accuracy of this model using realistic channel data. Based on these results, the PHY abstraction technique can be applied to evaluate the performance of GFDM based systems in an effective manner with low complexity. The maximum difference in the Packet Error Rate (PER) and throughput results in the abstraction case compared to bit-level simulation does not exceed 4% whilst offering a simulation time saving reduction of around 62,000 times.


title: Compressive sensing and truncated moment problems on spheres

abstract: We propose convex optimization algorithms to recover a good approximation of a point measure $\mu$ on the unit sphere $S\subseteq \mathbb{R}^n$ from its moments with respect to a set of real-valued functions $f_1,\dots, f_m$. Given a finite subset $C\subseteq S$ the algorithm produces a measure $\mu^*$ supported on $C$ and we prove that $\mu^*$ is a good approximation to $\mu$ whenever the functions $f_1,\dots, f_m$ are a sufficiently large random sample of independent Kostlan-Shub-Smale polynomials. More specifically, we give sufficient conditions for the validity of the equality $\mu=\mu^*$ when $\mu$ is supported on $C$ and prove that $\mu^*$ is close to the best approximation to $\mu$ supported on $C$ provided that all points in the support of $\mu$ are close to $C$.


title: Free deterministic equivalent Z-scores of compound Wishart models: A goodness of fit test of 2DARMA models

abstract: We introduce a new method to qualify the goodness of fit parameter estimation of compound Wishart models. Our method based on the free deterministic equivalent Z-score, which we introduce in this paper. Furthermore, an application to two dimensional autoregressive moving-average model is provided. Our proposal method is a generalization of statistical hypothesis testing to one dimensional moving average model based on fluctuations of real compound Wishart matrices, which is a recent result by Hasegawa, Sakuma and Yoshida.


title: Dynamic Social Balance and Convergent Appraisals via Homophily and Influence Mechanisms

abstract: Social balance theory describes allowable and forbidden configurations of the topologies of signed directed social appraisal networks. In this paper, we propose two discrete-time dynamical systems that explain how an appraisal network \textcolor{blue}{converges to} social balance from an initially unbalanced configuration. These two models are based on two different socio-psychological mechanisms respectively: the homophily mechanism and the influence mechanism. Our main theoretical contribution is a comprehensive analysis for both models in three steps. First, we establish the well-posedness and bounded evolution of the interpersonal appraisals. Second, we fully characterize the set of equilibrium points; for both models, each equilibrium network is composed by an arbitrary number of complete subgraphs satisfying structural balance. Third, we establish the equivalence among three distinct properties: non-vanishing appraisals, convergence to all-to-all appraisal networks, and finite-time achievement of social balance. In addition to theoretical analysis, Monte Carlo validations illustrates how the non-vanishing appraisal condition holds for generic initial conditions in both models. Moreover, numerical comparison between the two models indicate that the homophily-based model might be a more universal explanation for the formation of social balance. Finally, adopting the homophily-based model, we present numerical results on the mediation and globalization of local conflicts, the competition for allies, and the asymptotic formation of a single versus two factions.


title: Evolutionary games under incompetence

abstract: The adaptation process of a species to a new environment is a significant area of study in biology. As part of natural selection, adaptation is a mutation process which improves survival skills and reproductive functions of species. Here, we investigate this process by combining the idea of incompetence with evolutionary game theory. In the sense of evolution, incompetence and training can be interpreted as a special learning process. With focus on the social side of the problem, we analyze the influence of incompetence on behavior of species. We introduce an incompetence parameter into a learning function in a single-population game and analyze its effect on the outcome of the replicator dynamics. Incompetence can change the outcome of the game and its dynamics, indicating its significance within what are inherently imperfect natural systems.


title: $Q|SI\rangle$: A Quantum Programming Environment

abstract: This paper describes a quantum programming environment, named $Q|SI\rangle$. It is a platform embedded in the .Net language that supports quantum programming using a quantum extension of the $\mathbf{while}$-language. The framework of the platform includes a compiler of the quantum $\mathbf{while}$-language and a suite of tools for simulating quantum computation, optimizing quantum circuits, and analyzing and verifying quantum programs. Throughout the paper, using $Q|SI\rangle$ to simulate quantum behaviors on classical platforms with a combination of components is demonstrated. The scalable framework allows the user to program customized functions on the platform. The compiler works as the core of $Q|SI\rangle$ bridging the gap from quantum hardware to quantum software. The built-in decomposition algorithms enable the universal quantum computation on the present quantum hardware.


title: Effect of air confinement on thermal contact resistance in nanoscale heat transfer

abstract: We report herein the pressure dependent thermal contact resistance (Rc) between Wollaston wire thermal probe and samples which is evaluated within the framework of an analytical model. This work presents heat transfer between the Wollaston wire thermal probe and samples at nanoscale for different air pressure ranging from 1 Pa to 10e5 Pa. We make use of a scanning thermal microscopy (SThM) for the thermal analysis of two samples, fused silica (SiO2) and Titanium (Ti), with different thermal conductivities. The thermal probe's output voltage difference (deltaV) between out-off contact output voltage (Voc) and in-contact output voltage (Vic) was recorded. The result shows that the heat transfer increases with the increasing air pressure and it is found higher for higher thermal conductive material. We also propose an analytical model based on the normalized output voltage difference to extract the probe-sample thermal contact resistance. The obtained Rc values in the range of ~1.8e7 K/W to ~14.3e7 K/W validate the presented analytical model. Further analysis reveals that the thermal contact resistance between the probe and sample decreases with increasing air pressure. Such behavior is interpreted by the contribution of heat transfer through confined air on thermal contact resistance. In addition, the signature of Rc is also evidenced in the context of thermal mismatch behavior which is studied by using acoustic mismatch model (AMM) and diffuse mismatch model (DMM).