title: A General Definition of Network Communities and the Corresponding Detection Algorithm

abstract: Network structures, consisting of nodes and edges, have applications in almost all subjects. A set of nodes is called a community if the nodes have strong interrelations. Industries (including cell phone carriers and online social media companies) need community structures to allocate network resources and provide proper and accurate services. However, all the current detection algorithms are motivated by the practical problems, whose applicabilities in other fields are open to question. Thence, for a new community problem, researchers need to derive algorithms ad hoc, which is arduous and even unnecessary. In this paper, we represent a general procedure to find community structures in practice. We mainly focus on two typical types of networks: transmission networks and similarity networks. We reduce them to a unified graph model, based on which we propose a general method to define and detect communities. Readers can specialize our general algorithm to accommodate their problems. In the end, we also give a demonstration to show how the algorithm works.


title: Protecting Pegged Currency Markets from Speculative Investors

abstract: We consider a stochastic game between a trader and a central bank in a target zone market with a lower currency peg. This currency peg is maintained by the central bank through the generation of permanent price impact, thereby aggregating an ever increasing risky position in foreign reserves. We describe this situation mathematically by means of two coupled singular control problems, where the common singularity arises from a local time along a random curve. Our first result identifies a certain local time as that central bank strategy for which this risk position is minimized. We then consider the worst-case situation the central bank may face by identifying that strategy of the strategic investor that maximizes the expected inventory of the central bank under a cost criterion, thus establishing a Stackelberg equilibrium in our model.


title: Strong Sure Screening of Ultra-high Dimensional Data with Interaction Effects

abstract: Ultrahigh dimensional data sets are becoming increasingly prevalent in areas such as bioinformatics, medical imaging, and social network analysis. Sure independent screening of such data is commonly used to analyze such data. Nevertheless, few methods exist for screening for interactions among predictors. Moreover, extant interaction screening methods prove to be highly inaccurate when applied to data sets exhibiting strong interactive effects, but weak marginal effects, on the response. We propose a new interaction screening procedure based on joint cumulants which is not inhibited by such limitations. Under a collection of sensible conditions, we demonstrate that our interaction screening procedure has the strong sure screening property. Four simulations are used to investigate the performance of our method relative to two other interaction screening methods. We also apply a two-stage analysis to a real data example by first employing our proposed method, and then further examining a subset of selected covariates using multifactor dimensionality reduction.


title: An undergraduate laboratory study of the polarisation of annihilation photons using Compton scattering

abstract: An experiment for the advanced undergraduate laboratory which allows students to study the effect of photon polarisation in Compton scattering and to explore q\ uantum entanglement is described. The quantum entangled photons are produced through electron-positron annihilation in the $S$-state, and their polarisations a\ re analysed using the Compton scattering cross-section dependence on the photon polarisation. The equipment necessary for this experiment is available at a typ\ ical undergraduate physics laboratory. Finite geometry effects are discussed and investigated with the use of a Geant4-based simulation.


title: Composition inversion in mixtures of binary colloids and polymer

abstract: Understanding the phase behaviour of mixtures continues to pose challenges, even for systems that might be considered "simple". Here we consider a very simple mixture of two colloidal and one non-adsorbing polymer species which can be simplified even further to a size-asymmetrical binary mixture, in which the effective colloid-colloid interactions depend on the polymer concentration. We show that this basic system exhibits surprisingly rich phase behaviour. In particular, we enquire whether such a system features only a liquid-vapor phase separation (as in one-component colloid-polymer mixtures) or whether, additionally, liquid-liquid demixing of two colloidal phases can occur. Particle-resolved experiments show demixing-like behaviour, but when combined with bespoke Monte Carlo simulations, this proves illusory, and we reveal that only a single liquid-vapor transition occurs. Progressive migration of the small particles to the liquid phase as the polymer concentration increases gives rise to composition inversion - a maximum in the large particle concentration in the liquid phase. Near criticality the density fluctuations are found to be dominated by the larger colloids.


title: Equilibrium electric current of massive electrons with anomalous magnetic moments induced by a magnetic field and the electroweak interaction with matter

abstract: We study the possibility of the existence of the electric current, formed by massive electrons and positrons, flowing along an external magnetic field. The charged fermions are supposed to have nonzero anomalous magnetic moments and electroweakly interact with background matter. The expression for the current is obtained on the basis of the exact solution of the Dirac equation in the corresponding external fields. We demonstrate that, in the state of equilibrium, such a current is vanishing for any characteristics of the electron-positron plasma as well as the external fields. Our results are compared with the recent findings of other authors.


title: Data is the Fuel of Organizations: Opportunities and Challenges in Afghanistan

abstract: In this paper, the author at first briefly outlines the value of data in organizations and the opportunities and challenges in Afghanistan. Then the author takes the Kankor (National University Entrance Exam) data, particularly names of participants, locations, high schools and higher education institutions into account and explains how these data, that organizations in Afghanistan do not use for anything, can be useful in several cases and areas. The application of these data is shown through cases such as Auto filling missing values, identifying names of people, locations, and institutions from unstructured text, generating fake data to benchmark the database and web application performance and appearance, comparing and matching high school data with Kankor data, producing the top-n male and female names very common in Afghanistan or province-wise, and the data mining application in education and higher education institutions.


title: Multivariate analysis of short time series in terms of ensembles of correlation matrices

abstract: When dealing with non-stationary systems, for which many time series are available, it is common to divide time in epochs, i.e. smaller time intervals and deal with short time series in the hope to have some form of approximate stationarity on that time scale. We can then study time evolution by looking at properties as a function of the epochs. This leads to singular correlation matrices and thus poor statistics. In the present paper, we propose an ensemble technique to deal with a large set of short time series without any consideration of non-stationarity. Given a singular data matrix, we randomly select subsets of time series and thus create an ensemble of non-singular correlation matrices. As the selection possibilities are binomially large, we will obtain good statistics for eigenvalues of correlation matrices, which are typically not independent. Once we defined the ensemble, we analyze its behavior for constant and block-diagonal correlations and compare numerics with analytic results for the corresponding correlated Wishart ensembles. We discuss differences resulting from spurious correlations due to repetitive use of time-series. The usefulness of this technique should extend beyond the stationary case if, on the time scale of the epochs, we have quasi-stationarity at least for most epochs.


title: PointCNN: Convolution On $\mathcal{X}$-Transformed Points

abstract: We present a simple and general framework for feature learning from point clouds. The key to the success of CNNs is the convolution operator that is capable of leveraging spatially-local correlation in data represented densely in grids (e.g. images). However, point clouds are irregular and unordered, thus directly convolving kernels against features associated with the points, will result in desertion of shape information and variance to point ordering. To address these problems, we propose to learn an $\mathcal{X}$-transformation from the input points, to simultaneously promote two causes. The first is the weighting of the input features associated with the points, and the second is the permutation of the points into a latent and potentially canonical order. Element-wise product and sum operations of the typical convolution operator are subsequently applied on the $\mathcal{X}$-transformed features. The proposed method is a generalization of typical CNNs to feature learning from point clouds, thus we call it PointCNN. Experiments show that PointCNN achieves on par or better performance than state-of-the-art methods on multiple challenging benchmark datasets and tasks.


title: Contact Localization through Spatially Overlapping Piezoresistive Signals

abstract: Achieving high spatial resolution in contact sensing for robotic manipulation often comes at the price of increased complexity in fabrication and integration. One traditional approach is to fabricate a large number of taxels, each delivering an individual, isolated response to a stimulus. In contrast, we propose a method where the sensor simply consists of a continuous volume of piezoresistive elastomer with a number of electrodes embedded inside. We measure piezoresistive effects between all pairs of electrodes in the set, and count on this rich signal set containing the information needed to pinpoint contact location with high accuracy using regression algorithms. In our validation experiments, we demonstrate submillimeter median accuracy in locating contact on a 10mm by 16mm sensor using only four electrodes (creating six unique pairs). In addition to extracting more information from fewer wires, this approach lends itself to simple fabrication methods and makes no assumptions about the underlying geometry, simplifying future integration on robot fingers.


title: A New Correlation Coefficient for Aggregating Non-strict and Incomplete Rankings

abstract: We introduce a correlation coefficient that is designed to deal with a variety of ranking formats including those containing non-strict (i.e., with-ties) and incomplete (i.e., unknown) preferences. The correlation coefficient is designed to enforce a neutral treatment of incompleteness whereby no assumptions are made about individual preferences involving unranked objects. The new measure, which can be regarded as a generalization of the seminal Kendall tau correlation coefficient, is proven to satisfy a set of metric-like axioms and to be equivalent to a recently developed ranking distance function associated with Kemeny aggregation. In an effort to further unify and enhance both robust ranking methodologies, this work proves the equivalence of an additional distance and correlation-coefficient pairing in the space of non-strict incomplete rankings. These connections induce new exact optimization methodologies: a specialized branch and bound algorithm and an exact integer programming formulation. Moreover, the bridging of these complementary theories reinforces the singular suitability of the featured correlation coefficient to solve the general consensus ranking problem. The latter premise is bolstered by an accompanying set of experiments on random instances, which are generated via a herein developed sampling technique connected with the classic Mallows distribution of ranking data. Associated experiments with the branch and bound algorithm demonstrate that, as data becomes noisier, the featured correlation coefficient yields relatively fewer alternative optimal solutions and that the aggregate rankings tend to be closer to an underlying ground truth shared by a majority.


title: Simultaneous broadband observations and high-resolution X-ray spectroscopy of the transitional millisecond pulsar PSR J1023+0038

abstract: We report on the first simultaneous XMM-Newton, NuSTAR and Swift observations of the transitional millisecond pulsar PSR J1023+0038 in the X-ray active state. Our multi-wavelength campaign allowed us to investigate with unprecedented detail possible spectral variability over a broad energy range in the X-rays, as well as correlations and lags among emissions in different bands. The soft and hard X-ray emissions are significantly correlated, with no lags between the two bands. On the other hand, the X-ray emission does not correlate with the UV emission. We refine our model for the observed mode switching in terms of rapid transitions between a weak propeller regime and a rotation-powered radio pulsar state, and report on a detailed high-resolution X-ray spectroscopy using all XMM-Newton Reflection Grating Spectrometer data acquired since 2013. We discuss our results in the context of the recent discoveries on the system and of the state of the art simulations on transitional millisecond pulsars, and show how the properties of the narrow emission lines in the soft X-ray spectrum are consistent with an origin within the accretion disc.


title: A Generalized Muon Trajectory Estimation Algorithm with Energy Loss for Application to Muon Tomography

abstract: This work presents a generalized muon trajectory estimation (GMTE) algorithm to estimate the path of a muon in either uniform or nonuniform media. The use of cosmic ray muons in nuclear nonproliferation and safeguards verification applications has recently gained attention due to the nonintrusive and passive nature of the inspection, penetrating capabilities, as well as recent advances in detectors that measure position and direction of the individual muons before and after traversing the imaged object. However, muon image reconstruction techniques are limited in resolution due to low muon flux and the effects of multiple Coulomb scattering (MCS). Current reconstruction algorithms rely on overly simple assumptions for muon path estimation through the imaged object. For robust muon tomography, efficient and flexible physics based algorithms are needed to model the MCS process and accurately estimate the most probable trajectory of a muon as it traverses an object. In the present work, the use of a Bayesian framework and a Gaussian approximation of MCS are explored for estimation of the most likely path of a cosmic ray muon traversing uniform or nonuniform media and undergoing MCS. The algorithm's precision is compared to Monte Carlo simulated muon trajectories. It was found that the algorithm is expected to be able to predict muon tracks to less than 1.5 mm RMS for 0.5 GeV muons and 0.25 mm RMS for 3 GeV muons, a 50 percent improvement compared to SLP and 15 percent improvement when compared to PoCA. Further, a 30 percent increase in useful muon flux was observed relative to PoCA. Muon track prediction improved for higher muon energies or smaller penetration depth where energy loss is not significant. The effect of energy loss due to ionization is investigated, and a linear energy loss relation that is easy to use is proposed.


title: Linearized Filtering of Affine Processes Using Stochastic Riccati Equations

abstract: We consider an affine process $X$ which is only observed up to an additive white noise, and we ask for its law, for some time $t > 0 $, conditional on all observations up to this time $ t $. This is a general, possibly high dimensional filtering problem which is not even locally approximately Gaussian, whence essentially only particle filtering methods remain as solution techniques. In this work we present an efficient numerical solution by introducing an approximate filter for which conditional characteristic functions can be calculated by solving a system of generalized Riccati differential equations depending on the observation and the process characteristics of the signal $X$. The quality of the approximation can be controlled by easily observable quantities in terms of a macro location of the signal in state space. Asymptotic techniques as well as maximization techniques can be directly applied to the solutions of the Riccati equations leading to novel very tractable filtering formulas. The efficiency of the method is illustrated with numerical experiments for Cox-Ingersoll-Ross and Wishart processes, for which Gaussian approximations usually fail.


title: Band dependent inter-layer $f$-electron hybridization in CeRhIn$_5$

abstract: A key issue in heavy fermion research is how subtle changes in the hybridization between the 4$f$ (5$f$) and conduction electrons can result in fundamentally different ground states. CeRhIn$_5$ stands out as a particularly notable example: replacing Rh by either Co or Ir, located above or below Rh in the periodic table, antiferromagnetism gives way to superconductivity. In this photoemission study of CeRhIn$_5$, we demonstrate that the use of resonant ARPES with polarized light allows to extract detailed information on the 4$f$ crystal field states and details on the 4$f$ and conduction electron hybridization which together determine the ground state. We directly observe weakly dispersive Kondo resonances of $f$-electrons and identify two of the three Ce $4f_{5/2}^{1}$ crystal-electric-field levels and band-dependent hybridization, which signals that the hybridization occurs primarily between the Ce $4f$ states in the CeIn$_3$ layer and two more three-dimensional bands composed of the Rh $4d$ and In $5p$ orbitals in the RhIn$_2$ layer. Our results allow to connect the properties observed at elevated temperatures with the unusual low-temperature properties of this enigmatic heavy fermion compound.


abstract: We report observations of Ramsey interferences and spin echoes from electron spins inside a levitating macroscopic particle. The experiment is realized using nitrogen-vacancy (NV) centers hosted in a micron-sized diamond stored in a Paul trap both under atmospheric conditions and under vacuum. Spin echoes are used to show that the Paul trap preserves the coherence time of the embedded electron spins for more than microseconds. Conversely, the NV spin is employed to demonstrate high angular stability of the diamond even under vacuum. These results are significant steps towards strong coupling of NV spins to the rotational mode of levitating diamonds.


title: Koopman spectra in reproducing kernel Hilbert spaces

abstract: Every invertible, measure-preserving dynamical system induces a Koopman operator, which is a linear, unitary evolution operator acting on the $L^2$ space of observables associated with the invariant measure. Koopman eigenfunctions represent the quasiperiodic, or non-mixing, component of the dynamics. The extraction of these eigenfunctions and their associated eigenfrequencies from a given time series is a non-trivial problem when the underlying system has a dense point spectrum, or a continuous spectrum behaving similarly to noise. This paper describes methods for identifying Koopman eigenfrequencies and eigenfunctions from a discretely sampled time series generated by such a system with unknown dynamics. Our main result gives necessary and sufficient conditions for a Fourier function, defined on $N$ states sampled along an orbit of the dynamics, to be extensible to a Koopman eigenfunction on the whole state space, lying in a reproducing kernel Hilbert space (RKHS). In particular, we show that such an extension exists if and only if the RKHS norm of the Fourier function does not diverge as $ N \to \infty $, in which case the corresponding Fourier frequency is also a Koopman eigenfrequency. If such an RKHS extension does not exist, we can still construct an $L^2$ approximation of the eigenfunction. Numerical experiments on mixed-spectrum systems with weak periodic components demonstrate that this approach has significantly higher skill in identifying Koopman eigenfrequencies compared to conventional spectral estimation techniques based on the discrete Fourier transform.


title: qrypt0 - encrypted short messages exchanged between offline computers

abstract: A system is described for exchanging encrypted short messages between computers which remain permanently isolated from any network accessible to the attacker. The main advantage is effective protection of these computers from malware which could circumvent the encryption. For transmission, the ciphertext is passed between isolated and connected computers in the form of a QR code, which is displayed on and scanned from a screen. The security of qrypt0 therefore rests on the cryptography and the computer's physical isolation rather than on the computer security of the encrypting device.


title: Parametrized Equation of State for QCD from 3D Ising Model

abstract: The only first principle knowledge of the QCD equation of state at finite baryonic density is given from Lattice QCD as a Taylor expansion around $\mu_B = 0$. The coefficients of such an expansion are currently available up to order ${\cal O}(\mu_B^6)$. The expected critical behavior of QCD is in the same static universality class as the 3D Ising model. By means of a suitable parametrization for the scaling equation of state of 3D Ising and a parametrized map to connect to QCD, we present an equation of state matching first principle Lattice QCD calculations, which spans the values of baryonic densities explored in the BES-II program, and includes the correct scaling behavior in the proximity of the critical point. This EoS can serve as an important ingredient for the fluid dynamical simulations of heavy ion collisions at BES energies needed as a basis for the calculation of observables. Future comparisons between such calculations and BES-II data can constrain the parameters in the EoS -- including the parameters that describe the location of the critical point. This contribution reports on work done within the Fluctuations/Equation of State working group of the BEST Collaboration.


title: Phase transitions on C*-algebras arising from number fields and the generalized Furstenberg conjecture

abstract: In recent work, Cuntz, Deninger and Laca have studied the Toeplitz type C*-algebra associated to the affine monoid of algebraic integers in a number field, under a time evolution determined by the absolute norm. The KMS equilibrium states of their system are parametrized by traces on the C*-algebras of the semidirect products $J \rtimes O^*$ resulting from the multiplicative action of the units $O^*$ on integral ideals $J$ representing each ideal class. At each fixed inverse temperature $\beta > 2$, the extremal equilibrium states correspond to extremal traces of $C^*(J\rtimes O^*)$. Here we undertake the study of these traces using the transposed action of $O^*$ on the duals $\hat J$ of the ideals and the recent characterization of traces on transformation group C*-algebras due to Neshveyev. We show that the extremal traces of $C^*(J\rtimes O^*)$ are parametrized by pairs consisting of an ergodic invariant measure for the action of $O^*$ on $\hat{J}$ together with a character of the isotropy subgroup associated to the support of this measure. For every ideal the dual group $\hat {J}$ is a d-torus on which $O^*$ acts by linear toral automorphisms. Hence, the problem of classifying all extremal traces is a generalized version of Furstenberg's celebrated $\times 2$ $\times 3$ conjecture. We classify the results for various number fields in terms of ideal class group, degree, and unit rank, and we point along the way the trivial, the intractable, and the conjecturally classifiable cases. At the topological level, it is possible to characterize the number fields for which infinite $O^*$-invariant sets are dense in $\hat{J} $, thanks to a theorem of Berend; as an application we give a description of the primitive ideal space of $C^*(J\rtimes O^*)$ for those number fields.


title: Role of Internal Motions and Molecular Geometry on the NMR Relaxation of Hydrocarbons

abstract: The role of internal motions and molecular geometry on $^1$H NMR relaxation times $T_{1,2}$ in hydrocarbons is investigated using MD (molecular dynamics) simulations of the autocorrelation functions for in{\it tra}molecular $G_R(t)$ and in{\it ter}molecular $G_T(t)$ $^1$H-$^1$H dipole-dipole interactions arising from rotational ($R$) and translational ($T$) diffusion, respectively. We show that molecules with increased molecular symmetry such as neopentane, benzene, and isooctane show better agreement with traditional hard-sphere models than their corresponding straight-chain $n$-alkane, and furthermore that spherically-symmetric neopentane agrees well with the Stokes-Einstein theory. The influence of internal motions on the dynamics and $T_{1,2}$ relaxation of $n$-alkanes are investigated by simulating rigid $n$-alkanes and comparing with flexible (i.e. non-rigid) $n$-alkanes. Internal motions cause the rotational and translational correlation-times $\tau_{R,T}$ to get significantly shorter and the relaxation times $T_{1,2}$ to get significantly longer, especially for longer-chain $n$-alkanes. Site-by-site simulations of $^1$H's along the chains indicate significant variations in $\tau_{R,T}$ and $T_{1,2}$ across the chain, especially for longer-chain $n$-alkanes. The extent of the stretched (i.e. multi-exponential) decay in the autocorrelation functions $G_{R,T}(t)$ are quantified using inverse Laplace transforms, for both rigid and flexible molecules, and on a site-by-site bases. Comparison of $T_{1,2}$ measurements with the site-by-site simulations indicate that cross-relaxation (partially) averages-out the variations in $\tau_{R,T}$ and $T_{1,2}$ across the chain of long-chain $n$-alkanes. This work also has implications on the role of nano-pore confinement on the NMR relaxation of fluids in the organic-matter pores of kerogen and bitumen.


title: Vietnamese Open Information Extraction

abstract: Open information extraction (OIE) is the process to extract relations and their arguments automatically from textual documents without the need to restrict the search to predefined relations. In recent years, several OIE systems for the English language have been created but there is not any system for the Vietnamese language. In this paper, we propose a method of OIE for Vietnamese using a clause-based approach. Accordingly, we exploit Vietnamese dependency parsing using grammar clauses that strives to consider all possible relations in a sentence. The corresponding clause types are identified by their propositions as extractable relations based on their grammatical functions of constituents. As a result, our system is the first OIE system named vnOIE for the Vietnamese language that can generate open relations and their arguments from Vietnamese text with highly scalable extraction while being domain independent. Experimental results show that our OIE system achieves promising results with a precision of 83.71%.


title: The dense cores and filamentary structure of the molecular cloud in Corona Australis. Herschel SPIRE and PACS observations from the Herschel Gould Belt Survey

abstract: We present a catalogue of prestellar and starless cores within the Corona Australis molecular cloud using photometric data from the Herschel Space Observatory. At a distance of d~130 pc, Corona Australis is one of the closest star-forming regions. Herschel has taken multi-wavelength data of Corona Australis with both the SPIRE and PACS photometric cameras in a parallel mode with wavelengths in the range 70 {\mu}m to 500 {\mu}m. A complete sample of starless and prestellar cores and embedded protostars is identified. Other results from the Herschel Gould Belt Survey have shown spatial correlation between the distribution of dense cores and the filamentary structure within the molecular clouds. We go further and show correlations between the properties of these cores and their spatial distribution within the clouds, with a particular focus on the mass distribution of the dense cores with respect to their filamentary proximity. We find that only lower-mass starless cores form away from filaments, while all of the higher-mass prestellar cores form in close proximity to, or directly on the filamentary structure. This result supports the paradigm that prestellar cores mostly form on filaments. We analyse the mass distribution across the molecular cloud, finding evidence that the region around the Coronet appears to be at a more dynamically advanced evolutionary stage to the rest of the clumps within the cloud.


title: Estimating the Physical State of a Laboratory Slow Slipping Fault from Seismic Signals

abstract: Over the last two decades, strain and GPS measurements have shown that slow slip on earthquake faults is a widespread phenomenon. Slow slip is also inferred from correlated small amplitude seismic signals known as nonvolcanic tremor and low frequency earthquakes (LFEs). Slow slip has been reproduced in laboratory and simulation studies, however the fundamental physics of these phenomena and their relationship to dynamic earthquake rupture remains poorly understood. Here we show that, in a laboratory setting, continuous seismic waves are imprinted with fundamental signatures of the fault's physical state. Using machine learning on continuous seismic waves, we can infer several bulk characteristics of the fault (friction, shear displacement, gouge thickness), at any time during the slow slip cycle. This analysis also allows us to infer many properties of the future behavior of the fault, including the time remaining before the next slow slip event. Our work suggests that by applying machine learning approaches to continuous seismic data, new insight into the physics of slow slip could be obtained in Earth.


title: Algorithmic Bio-surveillance For Precise Spatio-temporal Prediction of Zoonotic Emergence

abstract: Viral zoonoses have emerged as the key drivers of recent pandemics. Human infection by zoonotic viruses are either spillover events -- isolated infections that fail to cause a widespread contagion -- or species jumps, where successful adaptation to the new host leads to a pandemic. Despite expensive bio-surveillance efforts, historically emergence response has been reactive, and post-hoc. Here we use machine inference to demonstrate a high accuracy predictive bio-surveillance capability, designed to pro-actively localize an impending species jump via automated interrogation of massive sequence databases of viral proteins. Our results suggest that a jump might not purely be the result of an isolated unfortunate cross-infection localized in space and time; there are subtle yet detectable patterns of genotypic changes accumulating in the global viral population leading up to emergence. Using tens of thousands of protein sequences simultaneously, we train models that track maximum achievable accuracy for disambiguating host tropism from the primary structure of surface proteins, and show that the inverse classification accuracy is a quantitative indicator of jump risk. We validate our claim in the context of the 2009 swine flu outbreak, and the 2004 emergence of H5N1 subspecies of Influenza A from avian reservoirs; illustrating that interrogation of the global viral population can unambiguously track a near monotonic risk elevation over several preceding years leading to eventual emergence.


title: Pebbling on Graph Products and other Binary Graph Constructions

abstract: Pebbling on graphs is a two-player game which involves repeatedly moving a pebble from one vertex to another by removing another pebble from the first vertex. The pebbling number $\pi(G)$ is the least number of pebbles required so that, regardless of the initial configuration of pebbles, a pebble can reach any vertex. Graham conjectured that the pebbling number for the cartesian product, $G \hspace{1mm}\square\hspace{1mm} H$, is bounded above by $\pi(G) \pi(H)$. We show that $\pi(G\hspace{1mm}\square\hspace{1mm} H) \le 2\pi(G) \pi(H)$ and, more sharply, that $\pi(G \hspace{1mm}\square\hspace{1mm} H) \le (\pi(G)+|G|) \pi(H)$. Furthermore, we provide similar results for other graph products and graph operations.


title: Statistical Learning For DC Optimal Power Flow

abstract: The optimal power flow problem plays an important role in the market clearing and operation of electric power systems. However, with increasing uncertainty from renewable energy operation, the optimal operating point of the system changes more significantly in real-time. In this paper, we aim at developing control policies that are able to track the optimal set-point with high probability. The approach is based on the observation that the OPF solution corresponding to a certain uncertainty realization is a basic feasible solution, which provides an affine control policy. The optimality of this basis policy is restricted to uncertainty realizations that share the same set of active constraints. We propose an ensemble control policy that combines several basis policies to improve performance. Although the number of possible bases is exponential in the size of the system, we show that only a few of them are relevant to system operation. We adopt a statistical learning approach to learn these important bases, and provide theoretical results that validate our observations. For most systems, we observe that efficient ensemble policies constructed using as few as ten bases, are able to obtain optimal solutions with high probability.


title: Life Beyond the Solar System: Observation and Modeling of Exoplanet Environments

abstract: The search for life on planets outside our solar system has largely been the province of the astrophysics community until recently. A major development since the NASA Astrobiology Strategy 2015 document (AS15) has been the integration of other NASA science disciplines (planetary science, heliophysics, Earth science) with ongoing exoplanet research in astrophysics. The NASA Nexus for Exoplanet System Science (NExSS) provides a forum for scientists to collaborate across disciplines to accelerate progress in the search for life elsewhere. Here we describe recent developments in these other disciplines, with a focus on exoplanet properties and environments, and the prospects for future progress that will be achieved by integrating emerging knowledge from astrophysics with insights from these fields.


title: Technology Needs for Detecting Life Beyond the Solar System: A White Paper in Support of the Astrobiology Science Strategy

abstract: In support of the Astrobiology Science Strategy, this whitepaper outlines some key technology challenges pertaining to the remote search for life in exoplanetary systems. Finding evidence for life on rocky planets outside of our solar system requires new technical capabilities for the key measurements of spectral signatures of biosignature gases, and of planetary mass measurement. Spectra of Earth-like planets can be directly measured in reflected stellar light in the visible band or near-infrared using a factor 1e-10 starlight suppression with occulters, either internal (coronagraph) or external (starshade). Absorption and emission (reflected and thermal) spectra can be obtained in the mid-infrared of rocky planets transiting M-dwarfs via spectroscopy of the transit and secondary eclipse, respectively. Mass can be measured from the star's reflex motion, the reflex motion of a star, via either precision radial velocity methods or astrometry. Several technology gaps must be closed to provide astronomers the necessary capabilities to obtain these key measurements for small planets orbiting within the predicted temperate zones around nearby stars. These involved performance improvements, in some cases, 1-2 orders of magnitude from state-of-the-art or involve performances never demonstrated. The technologies advancing to close these gaps have been identified through the NASA Exoplanet Exploration Program's annual Technology Selection and Prioritization Process in collaboration with the larger exoplanet science and technology community


title: A Wideband Self-Consistent Disk-Averaged Spectrum of Jupiter Near 30 GHz and Its Implications for NH$_{3}$ Saturation in the Upper Troposphere

abstract: We present a new set of measurements obtained with the Combined Array for Research in Millimeter-wave Astronomy (CARMA) of Jupiter's microwave thermal emission near the 1.3 cm ammonia (NH$_{3}$) absorption band. We use these observations to investigate the ammonia mole fraction in the upper troposphere, near $0.3 < P < 2$ bar, based on radiative transfer modeling. We find that the ammonia mole fraction must be $\sim2.4\times 10^{-4}$ below the NH$_{3}$ ice cloud, i.e., at $0.8 < P < 8$ bar, in agreement with results by de Pater et al. (2001, 2016a). We find the NH$_{3}$ cloud-forming region between $0.3 < P < 0.8$ bar to be globally sub-saturated by 55% on average, in accordance with the result in Gibson et al. (2005). Although our data are not very sensitive to the region above the cloud layer, we are able to set an upper limit of $2.4\times 10^{-7}$ on the mole fraction here, a factor of $\sim$10 above the saturated vapor curve.


title: Matter-neutrino resonance in a multi-angle neutrino bulb model

abstract: Simulations of neutrino flavor evolution in compact merger environments have shown that neutrino flavor, and hence nucleosynthesis, can be strongly affected by the presence of matter-neutrino resonances (MNRs), where there is a cancelation between the matter and the neutrino potential. Simulations performed thus far follow flavor evolution along a single neutrino trajectory, but self-consistency requires all trajectories to be treated simultaneously, and it has not been known whether MNR phenomena would still occur in multi-angle models. In this paper, we present the first fully multi-angle calculations of MNR. We find that familiar MNR phenomena, where neutrinos transform to a greater extent than anti-neutrinos and a feedback mechanism maintains the cancellation between the matter and neutrino potential, still occurs for a subset of angular bins, although the flavor transformation is not as efficient as in the single-angle case. In addition, we find other types of flavor transformation that are not seen in single-angle simulations. These flavor transformation phenomena appear to be robust and are present for a wide range of model parameters, as long as an MNR is present. Although computational constraints currently limit us to models with spherical symmetry, our results suggest that the presence of an MNR generally leads to large-scale neutrino flavor evolution in multi-angle systems.


title: Blockchain moderated by empty blocks to reduce the energetic impact of crypto-moneys

abstract: While cryptocurrencies and blockchain applications continue to gain popularity, their energy cost is evidently becoming unsustainable. In most instances, the main cost comes from the required amount of energy for the Proof-of-Work, and this cost is inherent to the design. In addition, useless costs from discarded work (e.g., the so-called Forks) and lack of scalability (in number of users and in rapid transactions) limit their practical effectiveness. In this paper, we present an innovative scheme which eliminates the nonce and thus the burden of the Proof-of-Work which is the main cause of the energy waste in cryptocurrencies such as Bitcoin. We prove that our scheme guarantees a tunable and bounded average number of simultaneous mining whatever the size of the population in competition, thus by making the use of nonce-based techniques unnecessary, achieves scalability without the cost of consuming a large volume of energy. The technique used in the proof of our scheme is based on the analogy of the analysis of a green leader election. The additional difference with Proof-of-Work schemes (beyond the suppression of the nonce field that is triggering most of the waste), is the introduction of (what we denote as) "empty blocks" which aim are to call regular blocks following a staircase set of values. Our scheme reduces the risk of Forks and provides tunable scalability for the number of users and the speed of block generation. We also prove using game theoretical analysis that our scheme is resilient to unfair competitive investments (e.g., "51 percent" attack) and block nursing.


title: Multivariate approximations in Wasserstein distance by Stein's method and Bismut's formula

abstract: Stein's method has been widely used for probability approximations. However, in the multi-dimensional setting, most of the results are for multivariate normal approximation or for test functions with bounded second- or higher-order derivatives. For a class of multivariate limiting distributions, we use Bismut's formula in Malliavin calculus to control the derivatives of the Stein equation solutions by the first derivative of the test function. Combined with Stein's exchangeable pair approach, we obtain a general theorem for multivariate approximations with near optimal error bounds on the Wasserstein distance.We apply the theorem to the unadjusted Langevin algorithm.


title: The Grassmannian complex and Goncharov's motivic complex in weight 4

abstract: For a field $F$ and a given integer $n>1$, Goncharov has given a complex $\Gamma_F(n)$ which he calls motivic and which he expects to rationally compute the weight $n$ motivic cohomology of $\text{Spec }F$, and hence its algebraic $K$-groups in Adams weight $n$, and he was also led to---conjecturally quasiisomorphic---`thickened' complexes thereof. These complexes involve tensor products of higher Bloch groups, the latter having been linked to the geometry of certain configurations in Goncharov's proof of Zagier's Polylogarithm Conjecture for weight 3, and an analogous picture has long been envisioned by Goncharov for higher weight as well. We provide a partial morphism in weight 4 by giving three out of four maps for configurations in general position.


title: Generalised Lyapunov Functions and Functionally Generated Trading Strategies

abstract: This paper investigates the dependence of functional portfolio generation, introduced by Fernholz (1999), on an extra finite variation process. The framework of Karatzas and Ruf (2017) is used to formulate conditions on trading strategies to be strong arbitrage relative to the market over sufficiently large time horizons. A mollification argument and Komlos theorem yield a general class of potential arbitrage strategies. These theoretical results are complemented by several empirical examples using data from the S&P 500 stocks.


title: Verification of real-time WSA-ENLIL+Cone simulations of CME arrival-time at the CCMC from 2010-2016

abstract: The Wang-Sheeley-Arge (WSA)-ENLIL+Cone model is used extensively in space weather operations world-wide to model CME propagation. As such, it is important to assess its performance. We present validation results of the WSA-ENLIL+Cone model installed at the Community Coordinated Modeling Center (CCMC) and executed in real-time by the CCMC space weather team. CCMC uses the WSA-ENLIL+Cone model to predict CME arrivals at NASA missions throughout the inner heliosphere. In this work we compare model predicted CME arrival-times to in-situ ICME leading edge measurements at STEREO-A, STEREO-B, and Earth (Wind and ACE) for simulations completed between March 2010-December 2016 (over 1,800 CMEs). We report hit, miss, false alarm, and correct rejection statistics for all three locations. For all predicted CME arrivals, the hit rate is 0.5, and the false alarm rate is 0.1. For the 273 events where the CME was predicted to arrive at Earth, STEREO-A, or STEREO-B, and was actually observed (hit event), the mean absolute arrival-time prediction error was 10.4 +/- 0.9 hours, with a tendency to early prediction error of -4.0 hours. We show the dependence of the arrival-time error on CME input parameters. We also explore the impact of the multi-spacecraft observations used to initialize the model CME inputs by comparing model verification results before and after the STEREO-B communication loss (since September 2014) and STEREO-A sidelobe operations (August 2014-December 2015). There is an increase of 1.7 hours in the CME arrival time error during single, or limited two-viewpoint periods, compared to the three-spacecraft viewpoint period. This trend would apply to a future space weather mission at L5 or L4 as another coronagraph viewpoint to reduce CME arrival time errors compared to a single L1 viewpoint.


title: The Zilber-Pink Conjecture and the Generalized Cosmetic Surgery Conjecture

abstract: In this paper, we generalize the Cosmetic Surgery Conjecture to an $n$-cusped hyperbolic $3$-manifold and prove it under the assumption of another well-known conjecture in number theory, so called the Zilber-Pink Conjecture. For $n=1$ and $2$, we show them without the assumption.


title: Evidence for Declination Dependence of the Ultrahigh Energy Cosmic Ray Spectrum in the Northern Hemisphere

abstract: Telescope Array (TA) is the largest experiment in the Northern Hemisphere studying ultrahigh energy cosmic rays. TA measurements of the cosmic ray spectrum using the surface detector have the best statistical power in the experiment, and observe the ankle of the spectrum and the high energy cutoff. When the data are divided into two declination bands, above and below 24.8 degrees, the cutoff appears at $10^{19.64 \pm 0.04}$ ($10^{19.84 \pm 0.02}$) eV in the lower (higher) band, an energy difference of 58\%. The global significance of the difference is 4.3 standard deviations. The lack of an instrumental cause of this difference implies it is astrophysical in nature.


title: A counterexample to Matsumoto's conjecture regarding absolute length vs. relative length in Finsler manifolds

abstract: Matsumoto conjectured that for any Finsler manifold $(M, F)$ for which the restriction of the fundamental tensor to the indicatrix of $F$ is positive definite, the absolute length $F(X)$ of any tangent vector $X \in T_xM$ is the global minimum for the relative length $|X|_y$ as $y$ varies along the indicatrix $I_x \subset T_xM$ of $F$. In this note, we disprove this conjecture by presenting a counterexample.


title: Partially Specified Spatial Autoregressive Model with Artificial Neural Network

abstract: Spatial autoregressive model, introduced by Clif and Ord in 1970s has been widely applied in many areas of science and econometrics such as regional economics, public finance, political sciences, agricultural economics, environmental studies and transportation analyses. As information technology grows rapidly, observations are seldom independent from others so a space autoregressive models can take this dependence into account and draw more reliable conclusions between covariates and the target variable itself. Based on the classical spatial model, Su and Jin proposed a semi-parametric model named as partially specified spatial autoregressive model (PSAR) to allow for more flexibility in modeling. And to estimate this nonparametric component, we use the neural network model which adds more flexibility to the classical model and allows for variations in the choice of activation functions according to different types of data. This paper extends an artificial neural network model to a partially specified space autoregressive model and proposes maximum likelihood estimators instead of quasi-maximum likelihood estimates. We establish the consistency and asymptotic normality of the estimators in this model. These results are obtained under some standard conditions in spatial models as well as neural network models. To illustrate, we investigate the quality of the normal approximation for finite samples by means of numerical simulation studies with three common choices of error distributions (standard normal, student-t distribution and the Laplace distribution). We apply our proposed model to a soil-water tension problem and a criminal study in Chicago. The results showed that our model can capture the spatial dependence between units as well as the unknown correlation structure between the target variable and covariates.


title: Multigrid for Staggered Lattice Fermions

abstract: Critical slowing down in Krylov methods for the Dirac operator presents a major obstacle to further advances in lattice field theory as it approaches the continuum solution. Here we formulate a multi-grid algorithm for the Kogut-Susskind (or staggered) fermion discretization which has proven difficult relative to Wilson multigrid due to its first-order anti-Hermitian structure. The solution is to introduce a novel spectral transformation by the K\"ahler-Dirac spin structure prior to the Galerkin projection. We present numerical results for the two-dimensional, two-flavor Schwinger model, however, the general formalism is agnostic to dimension and is directly applicable to four-dimensional lattice QCD.


title: A singularity-free analytic solution of rise dynamics of a liquid in a vertical cylindrical capillary

abstract: Capillary driven flow is a famous problem in fluid dynamics which dates back to Leonardo da Vinci. In this paper, we apply an analytic approximation method for highly nonlinear problem, namely the homotopy analysis method (HAM), to a model of the meniscus movement in a uniform vertical circular tube. Convergent explicit series solution is successfully obtained. Our results agree well with the numerical results given by the symbolic computing software Mathematica using six-order Runge-Kutta methods. More importantly, our analytic solution is valid in the whole region of physical parameters, and therefore can predict whether the path of liquid is monotonic or oscillatory. This kind of solution, to the best knowledge of the authors, has never been reported in the past, which might greatly deepen our understandings about capillarity.


title: On Small Beams with Large Topological Charge II: Photons, Electrons and Gravitational Waves

abstract: Beams of light with a large topological charge significantly change their spatial structure when they are focused strongly. Physically, it can be explained by an emerging electromagnetic field component in the direction of propagation, which is neglected in the simplified scalar wave picture in optics. Here we ask: Is this a specific photonic behavior, or can similar phenomena also be predicted for other species of particles? We show that the same modification of the spatial structure exists for relativistic electrons as well as for focused gravitational waves. However, this is for different physical reasons: For electrons, which are described by the Dirac equation, the spatial structure changes due to a Spin-Orbit coupling in the relativistic regime. In gravitational waves described with linearized general relativity, the curvature of space-time between the transverse and propagation direction leads to the modification of the spatial structure. Thus, this universal phenomenon exists for both massive and massless elementary particles with Spin 1/2, 1 and 2. It would be very interesting whether other types of particles such as composite systems (neutrons or C$_{60}$) or neutrinos show a similar behaviour and how this phenomenon can be explained in a unified physical way.


title: Estimating Heterogeneous Consumer Preferences for Restaurants and Travel Time Using Mobile Location Data

abstract: This paper analyzes consumer choices over lunchtime restaurants using data from a sample of several thousand anonymous mobile phone users in the San Francisco Bay Area. The data is used to identify users' approximate typical morning location, as well as their choices of lunchtime restaurants. We build a model where restaurants have latent characteristics (whose distribution may depend on restaurant observables, such as star ratings, food category, and price range), each user has preferences for these latent characteristics, and these preferences are heterogeneous across users. Similarly, each item has latent characteristics that describe users' willingness to travel to the restaurant, and each user has individual-specific preferences for those latent characteristics. Thus, both users' willingness to travel and their base utility for each restaurant vary across user-restaurant pairs. We use a Bayesian approach to estimation. To make the estimation computationally feasible, we rely on variational inference to approximate the posterior distribution, as well as stochastic gradient descent as a computational approach. Our model performs better than more standard competing models such as multinomial logit and nested logit models, in part due to the personalization of the estimates. We analyze how consumers re-allocate their demand after a restaurant closes to nearby restaurants versus more distant restaurants with similar characteristics, and we compare our predictions to actual outcomes. Finally, we show how the model can be used to analyze counterfactual questions such as what type of restaurant would attract the most consumers in a given location.


title: Semi-Supervised Convolutional Neural Networks for Human Activity Recognition

abstract: Labeled data used for training activity recognition classifiers are usually limited in terms of size and diversity. Thus, the learned model may not generalize well when used in real-world use cases. Semi-supervised learning augments labeled examples with unlabeled examples, often resulting in improved performance. However, the semi-supervised methods studied in the activity recognition literatures assume that feature engineering is already done. In this paper, we lift this assumption and present two semi-supervised methods based on convolutional neural networks (CNNs) to learn discriminative hidden features. Our semi-supervised CNNs learn from both labeled and unlabeled data while also performing feature learning on raw sensor data. In experiments on three real world datasets, we show that our CNNs outperform supervised methods and traditional semi-supervised learning methods by up to 18% in mean F1-score (Fm).


title: Using the Geometric Phase to Optimise Planar Somersaults

abstract: We derive the equations of motion for the planar somersault, which consist of two additive terms. The first is the dynamic phase that is proportional to the angular momentum, and the second is the geometric phase that is independent of angular momentum and depends solely on the details of the shape change. Next, we import digitised footage of an elite athlete performing 3.5 forward somersaults off the 3m springboard, and use the data to validate our model. We show that reversing and reordering certain sections of the digitised dive can maximise the geometric phase without affecting the dynamic phase, thereby increasing the overall rotation achieved. Finally, we propose a theoretical planar somersault consisting of four shape changing states, where the optimisation lies in finding the shape change strategy that maximises the overall rotation of the dive. This is achieved by balancing the rotational contributions from the dynamic and geometric phases, in which we show the geometric phase plays a small but important role in the optimisation process.


title: Dynamic Graph CNN for Learning on Point Clouds

abstract: Point clouds provide a flexible geometric representation suitable for countless applications in computer graphics; they also comprise the raw output of most 3D data acquisition devices. While hand-designed features on point clouds have long been proposed in graphics and vision, however, the recent overwhelming success of convolutional neural networks (CNNs) for image analysis suggests the value of adapting insight from CNN to the point cloud world. Point clouds inherently lack topological information so designing a model to recover topology can enrich the representation power of point clouds. To this end, we propose a new neural network module dubbed EdgeConv suitable for CNN-based high-level tasks on point clouds including classification and segmentation. EdgeConv acts on graphs dynamically computed in each layer of the network. It is differentiable and can be plugged into existing architectures. Compared to existing modules operating in extrinsic space or treating each point independently, EdgeConv has several appealing properties: It incorporates local neighborhood information; it can be stacked applied to learn global shape properties; and in multi-layer systems affinity in feature space captures semantic characteristics over potentially long distances in the original embedding. We show the performance of our model on standard benchmarks including ModelNet40, ShapeNetPart, and S3DIS.


title: Singular integrals of stable subordinator

abstract: It is well known that $\int_{0}^{1} t^{-\theta} d t<\infty$ for $\theta \in (0,1)$ and $\int_{0}^{1} t^{-\theta} d t=\infty$ for $\theta \in [1,\infty)$. Since $t$ can be taken as an $\alpha$-stable subordinator with $\alpha=1$, it is natural to ask whether $\int_{0}^{1} t^{-\theta} d S_{t}$ has a similar property when $S_{t}$ is an $\alpha$-stable subordinator with $\alpha \in (0,1)$. We show that $\theta=\frac 1\alpha$ is the border line such that $\int_{0}^{1} t^{-\theta} d S_{t}$ is finite a.s. for $\theta \in (0, \frac 1\alpha)$ and blows up a.s. for $\theta \in [\frac1\alpha,\infty)$. When $\alpha=1$, our result recovers that of $\int_{0}^{1} t^{-\theta} d t$. Moreover, we give a $p$-th moment estimate for the integral when $\theta \in (0,\frac 1\alpha)$.


title: New Constraints on Gliese 876 -- Exemplar of Mean-Motion Resonance

abstract: Gliese 876 harbors one of the most dynamically rich and well-studied exoplanetary systems. The nearby M4V dwarf hosts four known planets, the outer three of which are trapped in a Laplace mean-motion resonance. A thorough characterization of the complex resonant perturbations exhibited by the orbiting planets, and the chaotic dynamics therein, is key to a complete picture of the system's formation and evolutionary history. Here we present a reanalysis of the system using six years of new radial velocity (RV) data from four instruments. This new data augments and more than doubles the size of the decades-long collection of existing velocity measurements. We provide updated estimates of the system parameters by employing a computationally efficient Wisdom-Holman N-body symplectic integrator, coupled with a Gaussian Process (GP) regression model to account for correlated stellar noise. Experiments with synthetic RV data show that the dynamical characterization of the system can differ depending on whether a white noise or correlated noise model is adopted. Despite there being a region of stability for an additional planet in the resonant chain, we find no evidence for one. Our new parameter estimates place the system even deeper into resonance than previously thought and suggest that the system might be in a low energy, quasi-regular double apsidal corotation resonance. This result and others will be used in a subsequent study on the primordial migration processes responsible for the formation of the resonant chain.


title: Sensitivity of codispersion to noise and error in ecological and environmental data

abstract: Codispersion analysis is a new statistical method developed to assess spatial covariation between two spatial processes that may not be isotropic or stationary. Its application to anisotropic ecological datasets have provided new insights into mechanisms underlying observed patterns of species distributions and the relationship between individual species and underlying environmental gradients. However, the performance of the codispersion coefficient when there is noise or measurement error ("contamination") in the data has been addressed only theoretically. Here, we use Monte Carlo simulations and real datasets to investigate the sensitivity of codispersion to four types of contamination commonly seen in many real-world environmental and ecological studies. Three of these involved examining codispersion of a spatial dataset with a contaminated version of itself. The fourth examined differences in codisperson between plants and soil conditions, where the estimates of soil characteristics were based on complete or thinned datasets. In all cases, we found that estimates of codispersion were robust when contamination, such as data thinning, was relatively low (<15\%), but were sensitive to larger percentages of contamination. We also present a useful method for imputing missing spatial data and discuss several aspects of the codispersion coefficient when applied to noisy data to gain more insight about the performance of codispersion in practice.


title: Effects of controlling parameter on symbolic nonlinear complexity detection

abstract: Symbolic transformation, a coarse-graining process, is a crucial prerequisite for and has evidential influence to the symbolic time series analysis. We employ Shannon entropy for a parameter-dependent symbolization, KW (Kurths-Wessel) symbolic method, to test the effects of controlling parameter on its symbolic nonlinear complexity detection. Two chaotic models, logistic and Henon series, and heartbeats of CHF (Congestive Heart Failure) patients, healthy young and elderly subjects from PhysioNet are applied to test the KW symbolic entropy. The complexity-loss theory about aging and diseases in heart rates is validated and reasons that may account for some paradoxes in nonlinear analysis are discussed. Tests results suggest that due to different structural or dynamical properties of different nonlinear systems, controlling parameter of the KW symbolization should be adjusted accordingly to have reliable symbolic nonlinear analysis.


title: SkyMapper Southern Survey: First Data Release (DR1)

abstract: We present the first data release (DR1) of the SkyMapper Southern Survey, a hemispheric survey carried out with the SkyMapper Telescope at Siding Spring Observatory in Australia. Here, we present the survey strategy, data processing, catalogue construction and database schema. The DR1 dataset includes over 66,000 images from the Shallow Survey component, covering an area of 17,200 deg$^2$ in all six SkyMapper passbands $uvgriz$, while the full area covered by any passband exceeds 20,000 deg$^2$. The catalogues contain over 285 million unique astrophysical objects, complete to roughly 18 mag in all bands. We compare our $griz$ point-source photometry with PanSTARRS1 DR1 and note an RMS scatter of 2%. The internal reproducibility of SkyMapper photometry is on the order of 1%. Astrometric precision is better than 0.2 arcsec based on comparison with Gaia DR1. We describe the end-user database, through which data are presented to the world community, and provide some illustrative science queries.


title: QED Plasma at High Temperature

abstract: We demonstrate that the early universe behaved as a relativistic QED (Quantum Electrodynamics) plasma around the nucleosynthesis time while the temperature of the universe was below the neutrino decoupling temperature in the early universe. QED coupling constant becomes a temperature dependent parameter due to the radiative corrections to vacuum polarization in the early universe at nucleosynthesis temperature. Renormalization scheme of QED is used to calculate the effective parameters of relativistic plasma in the early universe. Renormalization constants of QED serve as effective parameters of the theory and are used to determine the behavior of matter. We explicitly compute the parameters of QED plasma such as Debye length and the plasma frequency as a function of temperature. Light is slowed down and trapped due to bending in such a medium and the frequency of electromagnetic radiation becomes a function of temperature as well.


title: Large spectral gaps for Steklov eigenvalues under volume constraints and under localized conformal deformations

abstract: In this paper we construct compact manifolds with fixed boundary geometry which admit Riemannian metrics of unit volume with arbitrarily large Steklov spectral gap. We also study the effect of localized conformal deformations that fix the boundary geometry. For instance, we prove that it is possible to make the spectral gap arbitrarily large using conformal deformations which are localized on domains of small measure, as long as the support of the deformations contains and connects each component of the boundary.


title: On the Fejes T\'oth Problem about the Sum of Angles Between Lines

abstract: In 1959 Fejes T\'oth posed a conjecture that the sum of pairwise non-obtuse angles between $N$ unit vectors in $\mathbb S^d$ is maximized by periodically repeated elements of the standard orthonormal basis. We obtain new improved upper bounds for this sum, as well as for the corresponding energy integral. We also provide several new approaches to the only settled case of the conjecture: $d=1$.


title: Morphism spaces in stable categories of Frobenius algebras

abstract: We present an explicit formula computing morphism spaces in the stable category of a Frobenius algebra.


title: Improved list-decodability of random linear binary codes

abstract: There has been a great deal of work establishing that random linear codes are as list-decodable as uniformly random codes, in the sense that a random linear binary code of rate $1 - H(p) - \epsilon$ is $(p,O(1/\epsilon))$-list-decodable with high probability. In this work, we show that such codes are $(p, H(p)/\epsilon + 2)$-list-decodable with high probability, for any $p \in (0, 1/2)$ and $\epsilon > 0$. In addition to improving the constant in known list-size bounds, our argument, which is quite simple, works simultaneously for all values of $p$, while previous works obtaining $L = O(1/\epsilon)$ patched together different arguments to cover different parameter regimes. Our approach is to strengthen an existential argument of (Guruswami, H{\aa}stad, Sudan and Zuckerman, IEEE Trans. IT, 2002) to hold with high probability. To complement our upper bound for random linear codes, we also improve an argument of (Guruswami, Narayanan, IEEE Trans. IT, 2014) to obtain an essentially tight lower bound of $1/\epsilon$ on the list size of uniformly random codes; this implies that random linear codes are in fact more list-decodable than uniformly random codes, in the sense that the list sizes are strictly smaller. To demonstrate the applicability of these techniques, we use them to (a) obtain more information about the distribution of list sizes of random linear codes and (b) to prove a similar result for random linear rank-metric codes.


title: Discovering intermediate-mass black hole lenses through gravitational wave lensing

abstract: Intermediate-mass black holes are the missing link that connects stellar-mass to supermassive black holes and are key to understanding galaxy evolution. Gravitational waves, like photons, can be lensed, leading to discernable effects such as diffraction or repeated signals. We investigate the detectability of intermediate-mass black hole deflectors in the LIGO-Virgo detector network. In particular, we simulate gravitational waves with variable source distributions lensed by an astrophysical population of intermediate-mass black holes, and use standard LIGO tools to infer the properties of these lenses. We find detections of intermediate-mass black holes at $98\%$ confidence level over a wide range of binary and lens parameters. Therefore, we conclude that intermediate-mass black holes could be detected through lensing of gravitational waves in the LIGO-Virgo detector network.


title: Why does hydronium diffuse faster than hydroxide in liquid water?

abstract: Proton transfer via hydronium and hydroxide ions in water is ubiquitous. It underlies acid-base chemistry, certain enzyme reactions, and even infection by the flu. Despite two-centuries of investigation, the mechanism underlying why hydronium diffuses faster than hydroxide in water is still not well understood. Herein, we employ state of the art Density Functional Theory based molecular dynamics, with corrections for nonlocal van der Waals interactions, and self-interaction in the electronic ground state, to model water and the hydrated water ions. At this level of theory, structural diffusion of hydronium preserves the previously recognized concerted behavior. However, by contrast, proton transfer via hydroxide is dominated by stepwise events, arising from a stabilized hyper-coordination solvation structure that discourages proton transfer. Specifically, the latter exhibits non-planar geometry, which agrees with neutron scattering results. Asymmetry in the temporal correlation of proton transfer enables hydronium to diffuse faster than hydroxide.

categories: cond-mat.soft cond-mat.mtrl-sci physics.chem-ph physics.comp-ph


title: Voltage-induced precessional switching at zero-bias magnetic field in a conically magnetized free layer

abstract: Voltage-induced magnetization dynamics in a conically magnetized free layer with an elliptic cylinder shape is theoretically studied on the basis of the macrospin model. It is found that an application of voltage pulse can induce the precessional switching of magnetization even at zero-bias magnetic field, which is of substantial importance for device applications such as voltage-controlled nonvolatile memory. Analytical expressions of the conditions for precessional switching are derived.


title: Symmetry in the extra-dimensional Yang-Mills theory and its Kaluza-Klein effective description

abstract: We construct an effective Lagrangian for Yang-Mills theories with $n$ extra dimensions. We start from a field theory governed by the extra-dimensional Poincar\'e group ISO$(1,3+n)$ and by the extended gauge group $SU(N,{\cal M}^{4+n})$, characterized by an unknown energy scale $\Lambda$ and assumed to be valid at energies far below this scale. Assuming that the size of the extra dimensions is much larger than the distance scale at which this theory is valid, we construct an effective theory with symmetry groups ISO$(1,3)$, $SU(N,{\cal M}^{4})$. Such theories are connected by a canonical transformation that hides the extended symmetries ISO$(1,3+n)$, $SU(N,{\cal M}^{4+n})$ into the standard ISO(1,3), $SU(N,{\cal M}^{4})$, thus generating KK gauge masses. Using a set of orthogonal functions $\{f^{(\underline{0})}, f^{(\underline{m})}(\bar x)\}$, defined by the Casimir invariant $\bar{P}^2$ of the translations subgroup $T(n)\subset$ISO$(n)$, we expand the degrees of freedom of ISO(1,3+n), $SU(N,{\cal M}^{4+n})$ in general Fourier series, whose coefficients are the degrees of freedom of ISO(1,3), $SU(N,{\cal M}^{4})$. These functions, which correspond to the projection on the basis $\{|\bar{x} \big >\}$ of the discrete basis $\{|0\big >,|p^{(\underline{m})}\big >\}$ generated by $\bar {P}^2$, are central to define the effective theory. Components along the ground state $f^{(\underline{0})}=\big <\bar x|0\big>$ do not receive mass at the compactification scale, so they are identified with the standard Yang-Mills fields; components along excited states $f^{(\underline{m})}=\big<\bar x|p^{(\underline{m})}\big>$ receive mass at this scale, so they correspond to KK excitations. Associated with any direction $|p^{(\underline{m})}\neq0\big>$ there are a massive gauge field and a pseudo-Goldstone boson. We stress resemblances of this mass-generating mechanism with the Englert-Higgs mechanism.


title: Server-Aided Revocable Predicate Encryption: Formalization and Lattice-Based Instantiation

abstract: Efficient user revocation is a necessary but challenging problem in many multi-user cryptosystems. Among known approaches, server-aided revocation yields a promising solution, because it allows to outsource the major workloads of system users to a computationally powerful third party, called the server, whose only requirement is to carry out the computations correctly. Such a revocation mechanism was considered in the settings of identity-based encryption and attribute-based encryption by Qin et al. (ESORICS 2015) and Cui et al. (ESORICS 2016), respectively. In this work, we consider the server-aided revocation mechanism in the more elaborate setting of predicate encryption (PE). The latter, introduced by Katz, Sahai, and Waters (EUROCRYPT 2008), provides fine-grained and role-based access to encrypted data and can be viewed as a generalization of identity-based and attribute-based encryption. Our contribution is two-fold. First, we formalize the model of server-aided revocable predicate encryption (SR-PE), with rigorous definitions and security notions. Our model can be seen as a non-trivial adaptation of Cui et al.'s work into the PE context. Second, we put forward a lattice-based instantiation of SR-PE. The scheme employs the PE scheme of Agrawal, Freeman and Vaikuntanathan (ASIACRYPT 2011) and the complete subtree method of Naor, Naor, and Lotspiech (CRYPTO 2001) as the two main ingredients, which work smoothly together thanks to a few additional techniques. Our scheme is proven secure in the standard model (in a selective manner), based on the hardness of the Learning With Errors (LWE) problem.


title: Floating potential of emitting surfaces in plasmas with respect to the space potential

abstract: The potential difference between a floating emitting surface and the plasma surrounding it has been described by several sheath models, including the space-charge-limited sheath, the electron sheath with high emission current, and the inverse sheath produced by charge-exchange ion trapping. Our measurements reveal that each of these models has its own regime of validity. We determine the potential of an emissive filament relative to the plasma potential, emphasizing variations in emitted current density and neutral particle density. The potential of a filament in a diffuse plasma is first shown to vanish, consistent with the electron sheath model and increasing electron emission. In a denser plasma with ample neutral pressure, the floating filament potential is positive, as predicted by a derived ion trapping condition. Lastly, the filament floated negatively in a third plasma, where flowing ions and electrons and nonnegligible electric fields may have disrupted ion trapping. Depending on the regime chosen, emitting surfaces can float positively or negatively with respect to the plasma potential.


title: Four-Qubit Monogamy and Four-Way Entanglement

abstract: We examine the various properties of the three four-qubit monogamy relations, all of which introduce the power factors in the three-way entanglement to reduce the tripartite contributions. On the analytic ground as much as possible we try to find the minimal power factors, which make the monogamy relations hold if the power factors are larger than the minimal powers. Motivated to the three-qubit monogamy inequality we also examine whether those four-qubit monogamy relations provide the SLOCC-invariant four-way entanglement measures or not. Our analysis indicate that this is impossible provided that the monogamy inequalities are derived merely by introducing weighting power factors.


title: Detectability of Small-Scale Dark Matter Clumps with Pulsar Timing Arrays

abstract: We examine the capability of pulsar timing arrays (PTAs) to detect very small-scale clumps of dark matter (DM), which are a natural outcome of the standard cold dark matter (CDM) paradigm. A clump streaming near the Earth or a pulsar induces an impulsive acceleration to encode residuals on pulsar timing data. We show that, assuming the standard abundance of DM clumps predicted by the CDM model, small-scale DM clumps with masses from $\sim 10^{-11} M_\odot$ to $\sim 10^{-8} \ M_\odot$ can be detectable by a PTA observation for a few decades with ${\cal O}(100)$ of pulsars with a timing noise of ${\cal O}(10)$ ns located at $\gtrsim 3$ kpc away from the Galactic center, as long as these mass scales are larger than the cutoff scale of the halo mass function that is determined by the particle nature of DM. Our result suggests that PTAs can provide a unique opportunity for testing one of the most fundamental predictions of the CDM paradigm. In addition, the detections and non-detections can constrain the cutoff mass scale inherent to the DM model.


title: Feeding Hand-Crafted Features for Enhancing the Performance of Convolutional Neural Networks

abstract: Since the convolutional neural network (CNN) is be- lieved to find right features for a given problem, the study of hand-crafted features is somewhat neglected these days. In this paper, we show that finding an appropriate feature for the given problem may be still important as they can en- hance the performance of CNN-based algorithms. Specif- ically, we show that feeding an appropriate feature to the CNN enhances its performance in some face related works such as age/gender estimation, face detection and emotion recognition. We use Gabor filter bank responses for these tasks, feeding them to the CNN along with the input image. The stack of image and Gabor responses can be fed to the CNN as a tensor input, or as a fused image which is a weighted sum of image and Gabor responses. The Gabor filter parameters can also be tuned depending on the given problem, for increasing the performance. From the extensive experiments, it is shown that the proposed methods provide better performance than the conventional CNN-based methods that use only the input images.


title: Probing vectorial near field of light: imaging theory and design principles of nanoprobes

abstract: Near-field microscopy is widely used for characterizing electromagnetic fields at nanoscale, where nanoprobes afford the opportunity to extract subwavelength optical quantities, including the amplitude, phase, polarization and chirality. However, owing to the complexity of various nanoprobes, a general and intuitive theory is highly needed to assess the vectorial field response of the nanoprobes and interpret the mechanism of the probe-field interaction. Here, we develop a general imaging theory based on the reciprocity of electromagnetism and multipole expansion analysis. The proposed theory closely resembles the multipolar Hamiltonian for light-matter interaction energy, revealing the coupling mechanism of the probe-field interaction. Based on this theory, we introduce a new paradigm for the design of functional nanoprobes by analyzing the reciprocal dipole moments, and establish effective design principles for the imaging of vectorial near fields. Moreover, we numerically analyze the responses of two typical probes, which can quantitatively reproduce and well explain the experimental results of previously reported measurements of optical magnetism and transverse spin angular momentum. Our work provides a powerful tool for the design and analysis of new functional probes that may enable the probing of various physical quantities of the vectorial near field.


title: Mid-infrared frequency comb generation via cascaded quadratic nonlinearities in quasi-phase-matched waveguides

abstract: We experimentally demonstrate a simple configuration for mid-infrared (MIR) frequency comb generation in quasi-phase-matched lithium niobate waveguides using the cascaded-$\chi^{(2)}$ nonlinearity. With nanojoule-scale pulses from an Er:fiber laser, we observe octave-spanning supercontinuum in the near-infrared with dispersive-wave generation in the 2.5--3 $\text{\mu}$m region and intra-pulse difference-frequency generation in the 4--5 $\text{\mu}$m region. By engineering the quasi-phase-matched grating profiles, tunable, narrow-band MIR and broadband MIR spectra are both observed in this geometry. Finally, we perform numerical modeling using a nonlinear envelope equation, which shows good quantitative agreement with the experiment---and can be used to inform waveguide designs to tailor the MIR frequency combs. Our results identify a path to a simple single-branch approach to mid-infrared frequency comb generation in a compact platform using commercial Er:fiber technology.


title: Zero-Delay Gaussian Joint Source-Channel Coding for the Interference Channel

abstract: This paper studies zero-delay joint source channel coding (JSCC) for transmission of correlated Gaussian sources over a Gaussian interference channel (GIC). We propose to adopt delay-free hybrid digital and analog (HDA) scheme, which is, transmitting the superposition of scaled source and its quantized version after applying scalar quantization to the source at each transmitter. At the corresponding receiver, two kinds of estimators are presented. It is shown that both the schemes, when optimized, beat the uncoded transmission if the channel signal-to-noise ratio (CSNR) is higher than a threshold value for different correlation coefficients and interference values.


title: Relaxed Schroedinger bridges and robust network routing

abstract: We consider network routing under random link failures with a desired final distribution. We provide a mathematical formulation of a relaxed transport problem where the final distribution only needs to be close to the desired one. The problem is a maximum entropy problem for path distributions with an extra terminal cost. We show that the unique solution may be obtained solving a generalized Schroedinger system. An iterative algorithm to compute the solution is provided. It contracts the Hilbert metric with contraction ratio less than 1/2 leading to extremely fast convergence.


title: Structured Triplet Learning with POS-tag Guided Attention for Visual Question Answering

abstract: Visual question answering (VQA) is of significant interest due to its potential to be a strong test of image understanding systems and to probe the connection between language and vision. Despite much recent progress, general VQA is far from a solved problem. In this paper, we focus on the VQA multiple-choice task, and provide some good practices for designing an effective VQA model that can capture language-vision interactions and perform joint reasoning. We explore mechanisms of incorporating part-of-speech (POS) tag guided attention, convolutional n-grams, triplet attention interactions between the image, question and candidate answer, and structured learning for triplets based on image-question pairs. We evaluate our models on two popular datasets: Visual7W and VQA Real Multiple Choice. Our final model achieves the state-of-the-art performance of 68.2% on Visual7W, and a very competitive performance of 69.6% on the test-standard split of VQA Real Multiple Choice.


title: Chorded pancyclicity in $k$-partite graphs

abstract: We prove that for any integers $p\geq k\geq 3$ and any $k$-tuple of positive integers $(n_1,\ldots ,n_k)$ such that $p=\sum _{i=1}^k{n_i}$ and $n_1\geq n_2\geq \ldots \geq n_k$, the condition $n_1\leq {p\over 2}$ is necessary and sufficient for every subgraph of the complete $k$-partite graph $K(n_1,\ldots ,n_k)$ with at least \[{{4 -2p+2n_1+\sum _{i=1}^{k} n_i(p-n_i)}\over 2}\] edges to be chorded pancyclic. Removing all but one edge incident with any vertex of minimum degree in $K(n_1,\ldots ,n_k)$ shows that this result is best possible. Our result implies that for any integers, $k\geq 3$ and $n\geq 1$, a balanced $k$-partite graph of order $kn$ with has at least ${{(k^2-k)n^2-2n(k-1)+4}\over 2}$ edges is chorded pancyclic. In the case $k=3$, this result strengthens a previous one by Adamus, who in 2009 showed that a balanced tripartite graph of order $3n$, $n \geq 2$, with at least $3n^2 - 2n + 2$ edges is pancyclic.


title: Fabrication and Characterization of Superconducting Quantum Interference Device using (Bi_{1-x}Sb_x)_2Se_3 Topological Insulator Nanoribbons

abstract: We report on the fabrication and electrical transport properties of superconducting quantum interference devices (SQUIDs) made from a (Bi_{1-x}Sb_x)_2Se_3 topological insulator (TI) nanoribbon (NR) connected with Pb0.5In0.5 superconducting electrodes. Below the transition temperature of the superconducting Pb0.5In0.5 electrodes, periodic oscillations of the critical current are observed in the TI NR SQUID under a magnetic field applied perpendicular to the plane owing to flux quantization. Also the output voltage modulates as a function of the external magnetic field. Moreover, the SQUID the SQUID shows a voltage modulation envelope, which is considered to represent the Fraunhofer-like patterns of each single junction. These properties of the TI NR SQUID would provide a useful method to explore Majorana fermions.


title: Random matrix theory for an adiabatically-varying oceanic acoustic waveguide

abstract: Problem of sound propagation in the ocean is considered. A novel approach of K. Hegewisch and S. Tomsovic for statistical modelling of acoustic wavefields in the random ocean is examined. The approach is based on construction of a wavefield propagator by means of random matrix theory. It is shown that this approach can be generalized onto acoustic waveguides with adiabatic longitudinal variations. Efficient generalization is obtained by means of stepwise approximation of the propagator. Accuracy of the generalized approach is confirmed numerically for a model of an underwater sound channel crossing a cold synoptic eddy. It is found that the eddy leads to substantial suppression of sound scattering.

categories: physics.ao-ph nlin.CD physics.comp-ph physics.flu-dyn


title: An accurate spectral method for Maxwell equations in Cole-Cole dispersive media

abstract: In this paper, we propose an accurate numerical means built upon a spectral-Galerkin method in spatial discretization and an enriched multi-step spectral-collocation approach in temporal direction, for Maxwell equations in Cole-Cole dispersive media in two-dimensional setting. Our starting point is to derive a new model involving only one unknown field from the original model with three unknown fields: electric, magnetic fields and the induced electric polarisation (described by a global temporal convolution of the electric field). This results in a second-order integral-differential equation with a weakly singular integral kernel expressed by the Mittag-Lefler (ML) function. The most interesting but challenging issue resides in how to efficiently deal with the singularity in time induced by the ML function which is an infinite series of singular power functions with different nature. With this in mind, we introduce a spectral-Galerkin method using Fourier-like basis functions for spatial discretization, leading to a sequence of decoupled temporal integral-differential equations (IDE) with the same weakly singular kernel involving the ML function as the original two-dimensional problem. With a careful study of the regularity of IDE, we incorporate several leading singular terms into the numerical scheme and approximate much regular part of the solution. Then we solve to IDE by a multi-step well-conditioned collocation scheme together with mapping technique to increase the accuracy and enhance the resolution. We show such an enriched collocation method is convergent and accurate. % analysis of the scheme is carried out.


title: Equidistribution of toral eigenfunctions along hypersurfaces

abstract: We prove a new quantum variance estimate for toral eigenfunctions. As an application, we show that, given any orthonormal basis of toral eigenfunctions and any smooth embedded hypersurface with nonvanishing principal curvatures, there exists a density one subsequence of eigenfunctions that equidistribute along the hypersurface. This is an analogue of the Quantum Ergodic Restriction theorems in the case of the flat torus, which in particular verifies the Bourgain-Rudnick's conjecture on $L^2$-restriction estimates for a density one subsequence of eigenfunctions in any dimension. Using our quantum variance estimates, we also obtain equidistribution of eigenfunctions against measures whose supports have Fourier dimension larger than $d-2$. In the end, we also describe a few quantitative results specific to dimension $2$.


title: The Travelling Salesman Problem and Adiabatic Quantum Computation: An Algorithm

abstract: An explicit algorithm for the travelling salesman problem is constructed in the framework of adiabatic quantum computation, AQC. The initial Hamiltonian for the AQC process admits canonical coherent states as the ground state, and the target Hamiltonian has the shortest tour as the desirable ground state. Some estimates/bounds are also given for the computational complexity of the algorithm with particular emphasis on the required energy resources, besides the space and time complexity, for the physical process of (quantum) computation in general.


title: Scalable and accurate deep learning for electronic health records

abstract: Predictive modeling with electronic health record (EHR) data is anticipated to drive personalized medicine and improve healthcare quality. Constructing predictive statistical models typically requires extraction of curated predictor variables from normalized EHR data, a labor-intensive process that discards the vast majority of information in each patient's record. We propose a representation of patients' entire, raw EHR records based on the Fast Healthcare Interoperability Resources (FHIR) format. We demonstrate that deep learning methods using this representation are capable of accurately predicting multiple medical events from multiple centers without site-specific data harmonization. We validated our approach using de-identified EHR data from two U.S. academic medical centers with 216,221 adult patients hospitalized for at least 24 hours. In the sequential format we propose, this volume of EHR data unrolled into a total of 46,864,534,945 data points, including clinical notes. Deep learning models achieved high accuracy for tasks such as predicting in-hospital mortality (AUROC across sites 0.93-0.94), 30-day unplanned readmission (AUROC 0.75-0.76), prolonged length of stay (AUROC 0.85-0.86), and all of a patient's final discharge diagnoses (frequency-weighted AUROC 0.90). These models outperformed state-of-the-art traditional predictive models in all cases. We also present a case-study of a neural-network attribution system, which illustrates how clinicians can gain some transparency into the predictions. We believe that this approach can be used to create accurate and scalable predictions for a variety of clinical scenarios, complete with explanations that directly highlight evidence in the patient's chart.


title: Improving Review Representations with User Attention and Product Attention for Sentiment Classification

abstract: Neural network methods have achieved great success in reviews sentiment classification. Recently, some works achieved improvement by incorporating user and product information to generate a review representation. However, in reviews, we observe that some words or sentences show strong user's preference, and some others tend to indicate product's characteristic. The two kinds of information play different roles in determining the sentiment label of a review. Therefore, it is not reasonable to encode user and product information together into one representation. In this paper, we propose a novel framework to encode user and product information. Firstly, we apply two individual hierarchical neural networks to generate two representations, with user attention or with product attention. Then, we design a combined strategy to make full use of the two representations for training and final prediction. The experimental results show that our model obviously outperforms other state-of-the-art methods on IMDB and Yelp datasets. Through the visualization of attention over words related to user or product, we validate our observation mentioned above.


title: Downlink Power Control in Massive MIMO Networks with Distributed Antenna Arrays

abstract: In this paper, we investigate downlink power control in massive multiple-input multiple-output (MIMO) networks with distributed antenna arrays. The base station (BS) in each cell consists of multiple antenna arrays, which are deployed in arbitrary locations within the cell. Due to the spatial separation between antenna arrays, the large-scale propagation effect is different from a user to different antenna arrays in a cell, which makes power control a challenging problem as compared to conventional massive MIMO. We assume that the BS in each cell obtains the channel estimates via uplink pilots. Based on the channel estimates, the BSs perform maximum ratio transmission for the downlink. We then derive a closed-form spectral efficiency (SE) expression, where the channels are subject to correlated fading. Utilizing the derived expression, we propose a max-min power control algorithm to ensure that each user in the network receives a uniform quality of service. Numerical results demonstrate that, for the network considered in this work, optimizing for max-min SE through the max-min power control improves the sum SE of the network as compared to equal power allocation.


title: Opinion Dynamics with Varying Susceptibility to Persuasion

abstract: A long line of work in social psychology has studied variations in people's susceptibility to persuasion -- the extent to which they are willing to modify their opinions on a topic. This body of literature suggests an interesting perspective on theoretical models of opinion formation by interacting parties in a network: in addition to considering interventions that directly modify people's intrinsic opinions, it is also natural to consider interventions that modify people's susceptibility to persuasion. In this work, we adopt a popular model for social opinion dynamics, and we formalize the opinion maximization and minimization problems where interventions happen at the level of susceptibility. We show that modeling interventions at the level of susceptibility lead to an interesting family of new questions in network opinion dynamics. We find that the questions are quite different depending on whether there is an overall budget constraining the number of agents we can target or not. We give a polynomial-time algorithm for finding the optimal target-set to optimize the sum of opinions when there are no budget constraints on the size of the target-set. We show that this problem is NP-hard when there is a budget, and that the objective function is neither submodular nor supermodular. Finally, we propose a heuristic for the budgeted opinion optimization and show its efficacy at finding target-sets that optimize the sum of opinions compared on real world networks, including a Twitter network with real opinion estimates.


title: Behavior Trees as a Representation for Medical Procedures

abstract: Behavior trees (BTs) emerged from video game development as a graphical language for modeling intelligent agent behavior. BTs have several properties which are attractive for modeling medical procedures including human-readability, authoring tools, and composability. This paper will illustrate construction of BTs for exemplary medical procedures. We are pleased to acknowledge support from National Science Foundation grant #IIS-1637444 and collaborations on that project with Johns Hopkins University and Worcester Polytechnic Institute.


title: Further Progress on the GM-MDS Conjecture for Reed-Solomon Codes

abstract: Designing good error correcting codes whose generator matrix has a support constraint, i.e., one for which only certain entries of the generator matrix are allowed to be non-zero, has found many recent applications, including in distributed coding and storage, multiple access networks, and weakly secure data exchange. The dual problem, where the parity check matrix has a support constraint, comes up in the design of locally repairable codes. The central problem here is to design codes with the largest possible minimum distance, subject to the given support constraint on the generator matrix. An upper bound on the minimum distance can be obtained through a set of singleton bounds, which can be alternatively thought of as a cut-set bound. Furthermore, it is well known that, if the field size is large enough, any random generator matrix obeying the support constraint will achieve the maximum minimum distance with high probability. Since random codes are not easy to decode, structured codes with efficient decoders, e.g., Reed-Solomon codes, are much more desirable. The GM-MDS conjecture of Dau et al states that the maximum minimum distance over all codes satisfying the generator matrix support constraint can be obtained by a Reed Solomon code. If true, this would have significant consequences. The conjecture has been proven for several special case: when the dimension of the code k is less than or equal to five, when the number of distinct support sets on the rows of the generator matrix m, say, is less than or equal to three, or when the generator matrix is sparsest and balanced. In this paper, we report on further progress on the GM-MDS conjecture. In particular, we show that the conjecture is true for all m less than equal to six. This generalizes all previous known results (except for the sparsest and balanced case, which is a very special support constraint).


title: Directed flow of open charm in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV using a quark coalescence model

abstract: The directed flow ($v_{1}$) of open charm meson ($D^{0}$) is studied in Au+Au collisions at $\sqrt{s_{NN}}$ = 200~GeV using A Multi-Phase Transport (AMPT) model framework with partonic interactions (string melting version). Within this framework, it is found that although the initial spatial eccentricity ($\epsilon_{1}$) of charm quark is smaller than light quarks, the charm quark $v_{1}$ magnitude is found to be approximately 7 times larger than that of the light u quark at large rapidity. This indicates that the charm quarks can retain more information from initial condition than the light quarks. We have studied the directed flow of $D^{0}$ as a function of rapidity and transverse momentum using quark coalescence as the mechanism for hadron production. Like charm quark, the $D^{0}$ $v_{1}$ magnitude is found to be about 7 times larger than that of the light ($\pi$) hadrons at large rapidity.