title: Structure Learning of Deep Networks via DNA Computing Algorithm

abstract: Convolutional Neural Network (CNN) has gained state-of-the-art results in many pattern recognition and computer vision tasks. However, most of the CNN structures are manually designed by experienced researchers. Therefore, auto- matically building high performance networks becomes an important problem. In this paper, we introduce the idea of using DNA computing algorithm to automatically learn high-performance architectures. In DNA computing algorithm, we use short DNA strands to represent layers and long DNA strands to represent overall networks. We found that most of the learned models perform similarly, and only those performing worse during the first runs of training will perform worse finally than others. The indicates that: 1) Using DNA computing algorithm to learn deep architectures is feasible; 2) Local minima should not be a problem of deep networks; 3) We can use early stop to kill the models with the bad performance just after several runs of training. In our experiments, an accuracy 99.73% was obtained on the MNIST data set and an accuracy 95.10% was obtained on the CIFAR-10 data set.


title: Boundary H\"older gradient estimates for the Monge-Amp\`ere equation

abstract: We investigate global H\"older gradient estimates for solutions to the Monge-Amp\`ere equation $$\mathrm{det}\;D^2 u=f\quad\mathrm{in}\;\Omega,$$ where the right-hand side $f$ is bounded away from $0$ and $\infty$. We consider two main situations when a) the domain $\Omega$ is uniformly convex and b) $\Omega$ is flat.


title: On Arveson's Boundary Theorem

abstract: This short note aims to give an insight to Arveson's boundary theorem by means of non-commutative Poisson boundaries and its applications.


title: SpiderBoost and Momentum: Faster Stochastic Variance Reduction Algorithms

abstract: SARAH and SPIDER are two recently developed stochastic variance-reduced algorithms, and SPIDER has been shown to achieve a near-optimal first-order oracle complexity in smooth nonconvex optimization. However, SPIDER uses an accuracy-dependent stepsize that slows down the convergence in practice, and cannot handle objective functions that involve nonsmooth regularizers. In this paper, we propose SpiderBoost as an improved scheme, which allows to use a much larger constant-level stepsize while maintaining the same near-optimal oracle complexity, and can be extended with proximal mapping to handle composite optimization (which is nonsmooth and nonconvex) with provable convergence guarantee. In particular, we show that proximal SpiderBoost achieves an oracle complexity of $\mathcal{O}(\min\{n^{1/2}\epsilon^{-2},\epsilon^{-3}\})$ in composite nonconvex optimization, improving the state-of-the-art result by a factor of $\mathcal{O}(\min\{n^{1/6},\epsilon^{-1/3}\})$. We further develop a novel momentum scheme to accelerate SpiderBoost for composite optimization, which achieves the near-optimal oracle complexity in theory and substantial improvement in experiments.


title: On the dynamics of ferrofluids: Global weak solutions to the Rosensweig system and rigorous convergence to equilibrium

abstract: This article establishes the global existence of weak solutions to a model proposed by Rosensweig (Rosensweig, Ferrohydrodynamics (1985)) for the dynamics of ferrofluids. The system is expressed by the conservation of linear momentum, the incompressibility condition, the conservation of angular momentum, and the evolution of the magnetization. The existence proof is inspired by the DiPerna-Lions theory of renormalized solutions. In addition, the rigorous relaxation limit of the equations of ferrohydrodynamics towards the quasi-equilibrium is investigated. The proof relies on the relative entropy method, which involves constructing a suitable functional, analyzing its time evolution and obtaining convergence results for the sequence of approximating solutions.


title: A New Class of Symmetric Distributions Including the Elliptically Symmetric Logistic

abstract: We introduce a new broad and exible class of multivariate elliptically symmetric distributions in- cluding the elliptically symmetric logistic and multivariate normal. Various probabilistic properties of the new distribution are studied, including the distribution of linear transformations, marginal distributions, conditional distributions, moments, stochastic representations and characteristic function.


title: Search for Multi-messenger Sources of Gravitational Waves and High-energy Neutrinos with Advanced LIGO during its first Observing Run, ANTARES and IceCube

abstract: Astrophysical sources of gravitational waves, such as binary neutron star and black hole mergers or core-collapse supernovae, can drive relativistic outflows, giving rise to non-thermal high-energy emission. High-energy neutrinos are signatures of such outflows. The detection of gravitational waves and high-energy neutrinos from common sources could help establish the connection between the dynamics of the progenitor and the properties of the outflow. We searched for associated emission of gravitational waves and high-energy neutrinos from astrophysical transients with minimal assumptions using data from Advanced LIGO from its first observing run O1, and data from the ANTARES and IceCube neutrino observatories from the same time period. We focused on candidate events whose astrophysical origin could not be determined from a single messenger. We found no significant coincident candidate, which we used to constrain the rate density of astrophysical sources dependent on their gravitational wave and neutrino emission processes.


title: Double-slit experiment at fermi scale: coherent photoproduction in heavy-ion collisions

abstract: The double-slit experiment has become a classic thought experiment, for its clarity in expressing the central puzzle of quantum mechanics -- wave-particle complementarity. Such wave-particle duality continues to be challenged and investigated in a broad range of entities with electrons, neutrons, helium atoms, C$_{60}$ fullerenes, Bose-Einstein condensates and biological molecules. All existing experiments are performed with the scale larger than angstrom. In this article, we present a double-slit scenario at fermi scale with new entities -- coherent photon products in heavy-ion collisions. Virtual photons from the electromagnetic fields of relativistic heavy ions can fluctuate to quark-antiquark pairs, scatter off a target nucleus and emerge as vector mesons. The two colliding nuclei can take turns to act as targets, forming a double-slit interference pattern. Furthermore, the `which-way' information can be partially solved by the violent strong interactions in the proposed scenario, which demonstrates a key concept of quantum mechanics -- complementary principle.


title: Spectral Embedding Norm: Looking Deep into the Spectrum of the Graph Laplacian

abstract: The extraction of clusters from a dataset which includes multiple clusters and a significant background component is a non-trivial task of practical importance. In image analysis this manifests for example in anomaly detection and target detection. The traditional spectral clustering algorithm, which relies on the leading $K$ eigenvectors to detect $K$ clusters, fails in such cases. In this paper we propose the {\it spectral embedding norm} which sums the squared values of the first $I$ normalized eigenvectors, where $I$ can be significantly larger than $K$. We prove that this quantity can be used to separate clusters from the background in unbalanced settings, including extreme cases such as outlier detection. The performance of the algorithm is not sensitive to the choice of $I$, and we demonstrate its application on synthetic and real-world remote sensing and neuroimaging datasets.


title: Electroweak monopoles and electroweak baryogenesis

abstract: We describe electroweak monopoles within the Born-Infeld extension of $SU(2)_L\times U(1)_Y$ electroweak theory. We argue for topological stability of these monopoles and computed their mass in terms of the Born-Infeld mass parameters. We then propose a new mechanism for electroweak baryogenesis which takes advantage of the following salient features of the electroweak monopoles: (i) monopoles support extra CP violation in the topological sector of the electroweak theory; (ii) they mediate unsuppressed baryon number violating interactions; (iii) non-thermal production of monopoles during the electroweak phase transitions generates departure from thermal equilibrium. We demonstrate that the observed baryon asymmetry of the universe can be explained in our theory in the presence of electroweak monopoles of mass $M\sim 10^{4}$ TeV.


title: COUSTIC: Combinatorial Double auction for Task Assignment in Device-to-Device Clouds

abstract: With the emerging technologies of Internet of Things (IOTs), the capabilities of mobile devices have increased tremendously. However, in the big data era, to complete tasks on one device is still challenging. As an emerging technology, crowdsourcing utilizing crowds of devices to facilitate large scale sensing tasks has gaining more and more research attention. Most of existing works either assume devices are willing to cooperate utilizing centralized mechanisms or design incentive algorithms using double auctions. Which is not practical to deal with the case when there is a lack of centralized controller for the former, and not suitable to the case when the seller device is also resource constrained for the later. In this paper, we propose a truthful incentive mechanism with combinatorial double auction for crowd sensing task assignment in device-to-device (D2D) clouds, where a single mobile device with intensive sensing task can hire a group of idle neighboring devices. With this new mechanism, time critical sensing tasks can be handled in time with a distributed nature. We prove that the proposed mechanism is truthful, individual rational, budget balance and computational efficient. Our simulation results demonstrate that combinatorial double auction mechanism gets a 26.3% and 15.8% gains in comparison to existing double auction scheme and the centralized maximum matching based algorithm respectively.


title: Antimagic orientations of disconnected even regular graphs

abstract: A $labeling$ of a digraph $D$ with $m$ arcs is a bijection from the set of arcs of $D$ to $\{1,2,\ldots,m\}$. A labeling of $D$ is $antimagic$ if no two vertices in $D$ have the same vertex-sum, where the vertex-sum of a vertex $u \in V(D)$ for a labeling is the sum of labels of all arcs entering $u$ minus the sum of labels of all arcs leaving $u$. An antimagic orientation $D$ of a graph $G$ is $antimagic$ if $D$ has an antimagic labeling. Hefetz, M$\ddot{u}$tze and Schwartz in [J. Graph Theory 64(2010)219-232] raised the question: Does every graph admits an antimagic orientation? It had been proved that for any integer $d$, every 2$d$-regular graph with at most two odd components has an antimagic orientation. In this paper, we consider the 2$d$-regular graph with many odd components. We show that every 2$d$-regular graph with any odd components has an antimagic orientation provide each odd component with enough order.


title: Complex Axis and de Medeiros' Campo Vetorial

abstract: The Complex Axis theorem states that any endomorphism of a finite-dimensional complex vector space affords an eigen-vector (or "invariant axis"). A geometric proof of this geometric result was given by A. de Medeiros, transforming the endomorphism into a topological self-map with Lefschetz number not equal to zero. We give a dual version of this proof, which may be more uniform, and does not rely on the need to do any calculation of an Euler characteristic or Lefschetz number. A vector field on Projective space is read off directly from the coordinates ("entries") of the given endomorphism (complex square matrix). A bordism is defined between such vector fields by means of Stokes' Theorem applied to a real manifold-with-boundary. This is the principle behind Hopf's lemma relating the Gauss map and the index of a vector field. All vector fields of the de Medeiros type are co-bordant to the Milnor-Hopf vector field. This latter comes from a non-derogatory, real diagonal endomorphism, so clearly possesses an eigen-vector. Therefore so has the given arbitrary endomorphism. The main theorem on complex polynomials naturally follows, using the companion matrix, secular polynomial reciprocity. The geometric Complex Axis derivation is meant to avoid determinants or "general position" arguments.


title: A Novel Mobile Edge Network Architecture with Joint Caching-Delivering and Horizontal Cooperation

abstract: Mobile edge caching/computing has been emerging as a promising paradigm to provide new services (e.g., ultra-high rate, ultra-reliable, and/or low-latency communications) in future wireless networks. In this paper, we introduce a novel mobile edge caching network architecture that leverages the optimal joint caching-delivering with horizontal cooperation among mobile edge nodes (MENs), namely JOCAD. Under this architecture, MENs cooperate with each other in both caching and delivering contents, aiming to simultaneously minimize the total average delay for the mobile users and mitigate the network traffic on the backhaul link. Extensive simulations demonstrate that the proposed solutions can reduce the total average delay for the whole network up to 40% compared with the most frequency-of-access policy, and up to 25% compared with locally optimal caching policy (i.e., without collaboration). Furthermore, the proposed solutions also increase the cache hit rate for the network by four times, thereby dramatically reducing the traffic load on the backhaul network.


title: Raman Spectroscopy of Diesel and Gasoline Engine-Out Soot Using Different Laser Power

abstract: We studied engine-out soot samples collected from a heavy-duty direct-injection diesel engine and a port-fuel injection gasoline spark-ignition engine. The two types of soot samples were characterized using Raman spectroscopy with different laser power. A Matlab program using least-square-method with trust-region-reflective algorithm was developed for curve fitting. We used a DOE (design of experiments) method to avoid local convergence. This method was used for two-band fitting and three-band fitting. The fitting results were used to determine the intensity ratio of D and G Raman bands. We find that high laser power may cause oxidation of soot samples, which gives higher D/G intensity ratio. Diesel soot has consistently higher amorphous/graphitic carbon ratio and thus higher oxidation reactivity, in comparison to gasoline soot, which is revealed by the higher D/G intensity ratio in Raman spectra measured under the same laser power.


title: Subgradient Descent Learns Orthogonal Dictionaries

abstract: This paper concerns dictionary learning, i.e., sparse coding, a fundamental representation learning problem. We show that a subgradient descent algorithm, with random initialization, can provably recover orthogonal dictionaries on a natural nonsmooth, nonconvex $\ell_1$ minimization formulation of the problem, under mild statistical assumptions on the data. This is in contrast to previous provable methods that require either expensive computation or delicate initialization schemes. Our analysis develops several tools for characterizing landscapes of nonsmooth functions, which might be of independent interest for provable training of deep networks with nonsmooth activations (e.g., ReLU), among numerous other applications. Preliminary experiments corroborate our analysis and show that our algorithm works well empirically in recovering orthogonal dictionaries.


title: K for the Price of 1: Parameter-efficient Multi-task and Transfer Learning

abstract: We introduce a novel method that enables parameter-efficient transfer and multi-task learning with deep neural networks. The basic approach is to learn a model patch - a small set of parameters - that will specialize to each task, instead of fine-tuning the last layer or the entire network. For instance, we show that learning a set of scales and biases is sufficient to convert a pretrained network to perform well on qualitatively different problems (e.g. converting a Single Shot MultiBox Detection (SSD) model into a 1000-class image classification model while reusing 98% of parameters of the SSD feature extractor). Similarly, we show that re-learning existing low-parameter layers (such as depth-wise convolutions) while keeping the rest of the network frozen also improves transfer-learning accuracy significantly. Our approach allows both simultaneous (multi-task) as well as sequential transfer learning. In several multi-task learning problems, despite using much fewer parameters than traditional logits-only fine-tuning, we match single-task performance.


title: Giant negative mobility of inertial particles caused by the periodic potential in steady laminar flows

abstract: Transport of an inertial particle advected by a two-dimensional steady laminar flow is numerically investigated in the presences of a constant force and a periodic potential. Within particular parameter regimes this system exhibits absolute negative mobility, which means that the particle can travel in a direction opposite to the constant force. It is found that the profile of the periodic potential plays an important role in the nonlinear response regime. Absolute negative mobility can be drastically enhanced by applying appropriate periodic potential, the parameter regime for this phenomenon becomes larger and the amplitude of negative mobility grows exceedingly large (giant negative mobility). In addition, giant positive mobility is also observed in the presence of appropriate periodic potential.


title: High-Dimensional Smoothing Splines and Application in Alzheimer's Disease Prediction Using Magnetic Resonance Imaging

abstract: Recent evidence has shown that structural magnetic resonance imaging (MRI) is an effective tool for Alzheimer's disease (AD) prediction and diagnosis. While traditional MRI-based diagnosis uses images acquired at a single time point, a longitudinal study is more sensitive and accurate in detecting early pathological changes of the AD. Two main difficulties arise in longitudinal MRI-based diagnosis: (1) the inconsistent longitudinal scans among subjects (i.e., different scanning time and different total number of scans); (2) the heterogeneous progressions of high-dimensional regions of interest (ROIs) in MRI. In this work, we propose a novel feature selection and estimation method which can be applied to extract features from the heterogeneous longitudinal MRI. A key ingredient of our method is the combination of smoothing splines and the $l_1$-penalty. We perform experiments on the Alzheimer's Disease Neuroimaging Initiative (ADNI) database. The results corroborate the advantages of the proposed method for AD prediction in longitudinal studies.


title: Integral and Absolute Hodge Classes

abstract: Despite the failure of the integral Hodge conjecture, we show that the rational Hodge conjecture implies an integral version (modulo torsion) of the absolute Hodge conjecture.


title: Injectivity in higher order complex domains

abstract: This is an expository survey of the Jacobian problem for the class of Pluriharmonic functions.


title: Learning with Interpretable Structure from Gated RNN

abstract: The interpretability of deep learning models has raised extended attention these years. It will be beneficial if we can learn an interpretable structure from deep learning models. In this paper, we focus on Recurrent Neural Networks~(RNNs) especially gated RNNs whose inner mechanism is still not clearly understood. We find that Finite State Automaton~(FSA) that processes sequential data has more interpretable inner mechanism according to the definition of interpretability and can be learned from RNNs as the interpretable structure. We propose two methods to learn FSA from RNN based on two different clustering methods. With the learned FSA and via experiments on artificial and real datasets, we find that FSA is more trustable than the RNN from which it learned, which gives FSA a chance to substitute RNNs in applications involving humans' lives or dangerous facilities. Besides, we analyze how the number of gates affects the performance of RNN. Our result suggests that gate in RNN is important but the less the better, which could be a guidance to design other RNNs. Finally, we observe that the FSA learned from RNN gives semantic aggregated states and its transition graph shows us a very interesting vision of how RNNs intrinsically handle text classification tasks.


title: Calculation of Critical Nucleation Rates by the Persistent Embryo Method: Application to Quasi Hard Sphere Models

abstract: We study crystal nucleation of the Weeks-Chandler-Andersen (WCA) model, using the recently introduced Persistent Embryo Method (PEM). The method provides detailed characterization of pre-critical, critical and post-critical nuclei, as well as nucleation rates that compare favorably with those obtained using other methods (umbrella sampling, forward flux sampling or seeding). We further map our results to a hard sphere model allowing to compare with other existing predictions. Implications for experiments are also discussed.


title: Quantum data compression by principal component analysis

abstract: Data compression can be achieved by reducing the dimensionality of high-dimensional but approximately low-rank datasets, which may in fact be described by the variation of a much smaller number of parameters. It often serves as a preprocessing step to surmount the curse of dimensionality and to gain efficiency, and thus it plays an important role in machine learning and data mining. In this paper, we present a quantum algorithm that compresses an exponentially large high-dimensional but approximately low-rank dataset in quantum parallel, by dimensionality reduction (DR) based on principal component analysis (PCA), the most popular classical DR algorithm. We show that the proposed algorithm achieves exponential speedup over the classical PCA algorithm when the original dataset are projected onto a polylogarithmically low-dimensional space. The compressed dataset can then be further processed to implement other tasks of interest, with significantly less quantum resources. As examples, we apply this algorithm to reduce data dimensionality for two important quantum machine learning algorithms, quantum support vector machine and quantum linear regression for prediction. This work demonstrates that quantum machine learning can be released from the curse of dimensionality to solve problems of practical importance.


title: First-principles investigation of Sc-III/IV under High Pressure

abstract: Using ab initio evolutionary structure prediction method in conjunction with density functional theory, we performed a systematic investigation on the structural transition of elemental scandium under pressure up to 250 GPa. Our prediction successfully reproduced several allotropes which have been reported in the previous literature, including the Sc-I, Sc-II and Sc-V. Moreover, we observed a series of energetically degenerate and geometrically similar structures at 110-195 GPa, which are likely to explain the unsolved phases III and IV reported by Akahama [Phys. Rev. Lett.,94, 19, 195503, (2005)]. A detailed comparison on powder X-ray diffraction pattern (PXRD) suggested that the Ccca-20 phase may account for the observed Sc-III, while Sc-IV is likely to be explained by a mixture of multiple energetically competing structures. We also used the candidate Sc-III structure as the model system to explore its superconducting behavior under pressures between 80-130 GPa. The predicted superconducting transition temperature Tc values are in satisfactory agreement with previous experimental results.


title: The effect of baseline layouts on the EoR foreground wedge: a semi-analytical approach

abstract: The 2D power spectrum is a cornerstone of the modern toolkit for analysis of the low-frequency radio interferometric observations of the 21 cm signal arising from the early Universe. Its familiar form disentangles a great deal of systematic information concerning both the sky and telescope, and displays as a foreground-dominated `brick' and `wedge' on large line-of-sight scales, and a complementary `window' on smaller scales. This paper builds on many previous works in the literature which seek to elucidate the varied instrumental and foreground factors which contribute to these familiar structures in the 2D power spectrum. In particular, we consider the effects of $uv$-sampling on the emergence of the wedge. Our results verify the expectation that arbitrarily dense instrument layouts in principal restore the missing information that leads to mode-mixing, and can therefore mitigate the wedge feature. We derive rule-of-thumb estimates for the required baseline density for complete wedge mitigation, showing that these will be unachievable in practice. We also discuss the optimal shape of the layout, showing that logarithmic regularity in the radial separation of baselines is favourable. While complete suppression of foreground leakage into the wedge is practically unachievable, we find that designing layouts which promote radial density and regularity is able to reduce the amplitude of foreground power by 1-3 orders of magnitude.


title: A Remarkably Loud Quasi-Periodicity after a Star is Disrupted by a Massive Black Hole

abstract: The immense tidal forces of massive black holes can rip apart stars that come too close to them. As the resulting stellar debris spirals inwards, it heats up and emits x-rays when near the black hole. Here, we report the discovery of an exceptionally stable 131-second x-ray quasi-periodicity from a black hole after it disrupted a star. Using a black hole mass indicated from host galaxy scaling relations implies that, (1) this periodicity originates from very close to the black hole's event horizon, and (2) the black hole is rapidly spinning. Our findings suggest that other disruption events with similar highly sensitive observations likely also exhibit quasi-periodicities that encode information about the fundamental properties of their black holes.


title: Preliminary results of using k-Nearest Neighbours Regression to estimate the redshift of radio selected datasets

abstract: In the near future, all-sky radio surveys are set to produce catalogues of tens of millions of sources with limited multi-wavelength photometry. Spectroscopic redshifts will only be possible for a small fraction of these new-found sources. In this paper, we provide the first in-depth investigation into the use of k-Nearest Neighbours Regression for the estimation of redshift of these sources. We use the Australia Telescope Large Area Survey radio data, combined with the Spitzer Wide-Area Infrared Extragalactic Survey infra-red, the Dark Energy Survey optical and the Australian Dark Energy Survey spectroscopic survey data. We then reduce the depth of photometry to match what is expected from upcoming Evolutionary Map of the Universe survey, testing against both data sets. To examine the generalisation of our methods, we test one of the sub-fields of Australia Telescope Large Area Survey against the other. We achieve an outlier rate of ~10% across all tests, showing that the k-Nearest Neighbours regression algorithm is an acceptable method of estimating redshift, and would perform better given a sample training set with uniform redshift coverage.


title: Coherent optical and acoustic phonons generated at lattice-matched GaP/Si(001) heterointerfaces

abstract: Thin GaP films can be grown on an exact Si(001) substrate with nearly perfect lattice match. We perform all-optical pump-probe measurements to investigate the ultrafast electron-phonon coupling at the buried interface of GaP/Si. Above-bandgap excitation with a femtosecond laser pulse can induce coherent longitudinal optical (LO) phonons both in the GaP opverlayer and in the Si substrate. The coupling of the GaP LO phonons with photoexcited plasma is reduced significantly with decreasing the GaP layer thickness from 56 to 16 nm due to the quasi-two-dimensional confinement of the plasma. The same laser pulse can also generate coherent longitudinal acoustic (LA) phonons in the form of a strain pulse. The strain induces not only a periodic modulation in the optical reflectivity as they propagate in the semiconductors, but also a sharp spike when it arrives at the GaP layer boundaries. The acoustic pulse induced at the GaP/Si interface is remarkably stronger than that at the Si surface, suggesting a possible application of the GaP/Si heterostructure as an opto-acoustic transducer. The amplitude and the phase of the reflectivity modulation varies with the GaP layer thickness, which can be understood in terms of the interference caused by the multiple acoustic pulses generated at the top surface and at the buried interface.


title: On the Zeroes of Half-Integral Weight Eisenstein Series on $\Gamma_0(4)$

abstract: We locate all but $O(\sqrt{k\log{k}})$ zeroes of the half integral weight Eisenstein series $E_\infty(z,k)$ of $\Gamma_0(4)$ for $k$ sufficiently large. To do this, we relate $E_\infty(z,k)$ to $\Gamma_0(4)$'s other Eisenstein series, $E_0(z,k)$ and $E_\frac{1}{2}(z,k)$, which are easier to study in the regions of which zeroes occur. We will use trigonometric approximations of $E_0(z,k)$ and $E_\frac{1}{2}(z,k)$ in order to locate the zeroes.


title: Positive one-point commuting difference operators

abstract: In this paper we study commuting difference operators containing a shift operator with only positive degrees. We construct examples of such operators in the case of hyperelliptic spectral curves.


title: Beamforming Optimization for Intelligent Reflecting Surface with Discrete Phase Shifts

abstract: Intelligent reflecting surface (IRS) is a promising technology for achieving high spectrum efficiency in future wireless networks by leveraging massive low-cost reflecting elements with each reflecting the incident signal with a proper phase shift. However, prior works on IRS are mainly based on the optimization of infinite-resolution phase shifters which are practically infeasible due to hardware imperfections. In contrast, we study in this paper an IRS-aided wireless network, where an IRS with only finite-resolution phase shifter available at each element is deployed to assist in the communication from a multi-antenna access point (AP) to a single-antenna user. We aim to minimize the transmit power at the AP by jointly optimizing the transmit beamforming at the AP and reflect beamforming at the IRS, subject to the signal-to-noise ratio (SNR) constraint and practical discrete phase shift constraints. We first propose a suboptimal but low-complexity algorithm by exploiting the alternating optimization technique. Then, we reveal that as in the case with continuous phase shifts, the IRS with discrete phase shifts also achieves the squared power gain for asymptotically large number of reflecting elements, despite suffering a performance loss that depends only on the resolution of phase shifters.


title: The outer halo globular cluster system of M31 - III. Relationship to the stellar halo

abstract: We utilise the final catalogue from the Pan-Andromeda Archaeological Survey to investigate the links between the globular cluster system and field halo in M31 at projected radii $R_p=25-150$ kpc. In this region the cluster radial density profile exhibits a power-law decline with index $\Gamma=-2.37\pm0.17$, matching that for the stellar halo component with [Fe/H] $<-1.1$. Spatial density maps reveal a striking correspondence between the most luminous substructures in the metal-poor field halo and the positions of many globular clusters. By comparing the density of metal-poor halo stars local to each cluster with the azimuthal distribution at commensurate radius, we reject the possibility of no correlation between clusters and field overdensities with high confidence. We use our stellar density measurements and previous kinematic data to demonstrate that $\approx35-60\%$ of clusters exhibit properties consistent with having been accreted into the outskirts of M31 at late times with their parent dwarfs. Conversely, at least $\sim40\%$ of remote clusters show no evidence for a link with halo substructure. The radial density profile for this subgroup is featureless and closely mirrors that observed for the apparently smooth component of the metal-poor stellar halo. We speculate that these clusters are associated with the smooth halo; if so, their properties appear consistent with a scenario where the smooth halo was built up at early times via the destruction of primitive satellites. In this picture the entire M31 globular cluster system outside $R_p=25$ kpc comprises objects accumulated from external galaxies over a Hubble time of growth.


title: An embedded X-ray source shines through the aspherical AT2018cow: revealing the inner workings of the most luminous fast-evolving optical transients

abstract: We present the first extensive radio to gamma-ray observations of a fast-rising blue optical transient (FBOT), AT2018cow, over its first ~100 days. AT2018cow rose over a few days to a peak luminosity $L_{pk}\sim4\times 10^{44}$ erg/s exceeding those of superluminous supernovae (SNe), before declining as $\propto t^{-2}$. Initial spectra at $\lesssim 15$ days were mostly featureless and indicated large expansion velocities v~0.1c and temperatures reaching 30000 K. Later spectra revealed a persistent optically-thick photosphere and the emergence of H and He emission features with v~sim 4000 km/s with no evidence for ejecta cooling. Our broad-band monitoring revealed a hard X-ray spectral component at $E\ge 10$ keV, in addition to luminous and highly variable soft X-rays, with properties unprecedented among astronomical transients. An abrupt change in the X-ray decay rate and variability appears to accompany the change in optical spectral properties. AT2018cow showed bright radio emission consistent with the interaction of a blastwave with $v_{sh}$~0.1c with a dense environment ($\dot M\sim10^{-3}-10^{-4}\,M_{\odot}yr^{-1}$ for $v_w=1000$ km\s). While these properties exclude Ni-powered transients, our multi-wavelength analysis instead indicates that AT2018cow harbored a "central engine", either a compact object (magnetar or black hole) or an embedded internal shock produced by interaction with a compact, dense circumstellar medium. The engine released $\sim10^{50}-10^{51.5}$ erg over $\sim10^3-10^5$ s and resides within low-mass fast-moving material with equatorial-polar density asymmetry ($M_{ej,fast}\lesssim0.3\,\rm{M_{\odot}}$). Successful SNe from low-mass H-rich stars (like electron-capture SNe) or failed explosions from blue supergiants satisfy these constraints. Intermediate-mass black-holes are disfavored by the large environmental density probed by the radio observations.


title: Searching for GeV gamma-ray emission from the bulge of M31

abstract: The three major large-scale, diffuse $\gamma$-ray structures of the Milky Way are the Galactic disk, a bulge-like GeV excess towards the Galactic center, and the Fermi bubble. Whether such structures can also be present in other normal galaxies remains an open question. M31, as the nearest massive normal galaxy, holds promise for spatially-resolving the $\gamma$-ray emission. Based on more than 8 years of Fermi-LAT observations, we use (1) disk, (2) bulge, and (3) disk-plus-bulge templates to model the spatial distribution of the $\gamma$-ray emission from M31. Among these, the disk-plus-bulge template delivers the best-fit, in which the bulge component has a TS value 25.7 and a photon-index of $2.57\pm0.17$, providing strong evidence for a centrally-concentrated gamma-ray emission from M31, that is analogous to the Galactic center excess. The total 0.2-300 GeV gamma-ray luminosity from this bulge component is $(1.16\pm0.14)\times10^{38} erg s^{-1}$, which would require $\sim1.5\times10^5$ millisecond pulsars, if they were the dominant source. We also search for a Fermi bubble-like structure in M31 using the full dataset (pass8), but no significant evidence is found. In addition, a likelihood analysis using only photons with the most accurate reconstructed direction (i.e., PSF3-only data) reveals a 4.8 $\sigma$ point-like source located at $\sim$ 10kpc to the northwest of the M31 disk, with a luminosity of $(0.97\pm0.27)\times10^{38} erg s^{-1}$ and a photon-index of $2.31\pm0.18$. Lacking of a counterpart on the southeast side of the disk, the relation between this point-like source and a bubble-like structure remains elusive.


title: Joint modelling of progression-free and overall survival and computation of correlation measures

abstract: In this paper, we derive the joint distribution of progression-free and overall survival as a function of transition probabilities in a multistate model. No assumptions on copulae or latent event times are needed and the model is allowed to be non-Markov. From the joint distribution, statistics of interest can then be computed. As an example, we provide closed formulas and statistical inference for Pearson's correlation coefficient between progression-free and overall survival in a parametric framework. The example is inspired by recent approaches to quantify the dependence between progression-free survival, a common primary outcome in phase III trials in oncology, and overall survival. We complement these approaches by providing methods of statistical inference while at the same time working within a much more parsimonious modelling framework. Our approach is completely general and can be applied to other measures of dependence. We also discuss extensions to nonparametric inference. Our analytical results are illustrated using a large randomized clinical trial in breast cancer.


title: Urban Healthcare Big Data System Based on Crowdsourced and Cloud-Based Air Quality Indicators

abstract: The ever-accelerating process of globalization enables more than half the population to live in cities. Thus, the air quality in cities exerts critical influence on the health status of more and more urban residents. In this article, based on urban air quality data collected through meteorological sites, mobile crowdsourcing, and IoT sensing, along with users' body signals, we propose an urban healthcare big data system named UH-BigDataSys. In this article, we first introduce a method of integrating multi-source air quality data for the data preparation of artificial-intelligence-based smart urban services. Then a testbed of UH-BigDataSys is set up with the deployment of air-quality-aware healthcare applications. Finally, we provide health guidance for urban residents in aspects of respiratory diseases, outdoor travel, sleep quality, and so on. The ultimate goal of UH-BigDataSys is for urban residents to lead healthier lives.


title: Generalized Beamspace Modulation Using Multiplexing: A Breakthrough in mmWave MIMO

abstract: Spatial multiplexing (SMX) multiple-input multiple-output (MIMO) over the best beamspace was considered as the best solution for millimeter wave (mmWave) communications regarding spectral efficiency (SE), referred as the best beamspace selection (BBS) solution. The equivalent MIMO water-filling (WF-MIMO) channel capacity was treated as an unsurpassed SE upper bound. Recently, researchers have proposed various schemes trying to approach the benchmark and the performance bound. But, are they the real limit of mmWave MIMO systems with reduced radio-frequency (RF) chains? In this paper, we challenge the benchmark and the corresponding bound by proposing a better transmission scheme that achieves higher SE, namely the Generalized Beamspace Modulation using Multiplexing (GBMM). Inspired by the concept of spatial modulation, besides the selected beamspace, the selection operation is used to carry information. We prove that GBMM is superior to BBS in terms of SE and can break through the well known `upper bound'. That is, GBMM renews the upper bound of the SE. We investigate SE-oriented precoder activation probability optimization, fully-digital precoder design, optimal power allocation and hybrid precoder design for GBMM. A gradient ascent algorithm is developed to find the optimal solution, which is applicable in all signal-to-noise-ratio (SNR) regimes. The best solution is derived in the high SNR regime. Additionally, we investigate the hybrid receiver design and deduce the minimum number of receive RF chains configured to gain from GBMM in achievable SE. We propose a coding approach to realize the optimized precoder activation. An extension to mmWave broadband communications is also discussed. Comparisons with the benchmark (i.e., WF-MIMO channel capacity) are made under different system configurations to show the superiority of GBMM.


title: Convolutional Deblurring for Natural Imaging

abstract: In this paper, we propose a novel design of image deblurring in the form of one-shot convolution filtering that can directly convolve with naturally blurred images for restoration. The problem of optical blurring is a common disadvantage to many imaging applications that suffer from optical imperfections. Despite numerous deconvolution methods that blindly estimate blurring in either inclusive or exclusive forms, they are practically challenging due to high computational cost and low image reconstruction quality. Both conditions of high accuracy and high speed are prerequisites for high-throughput imaging platforms in digital archiving. In such platforms, deblurring is required after image acquisition before being stored, previewed, or processed for high-level interpretation. Therefore, on-the-fly correction of such images is important to avoid possible time delays, mitigate computational expenses, and increase image perception quality. We bridge this gap by synthesizing a deconvolution kernel as a linear combination of Finite Impulse Response (FIR) even-derivative filters that can be directly convolved with blurry input images to boost the frequency fall-off of the Point Spread Function (PSF) associated with the optical blur. We employ a Gaussian low-pass filter to decouple the image denoising problem for image edge deblurring. Furthermore, we propose a blind approach to estimate the PSF statistics for two Gaussian and Laplacian models that are common in many imaging pipelines. Thorough experiments are designed to test and validate the efficiency of the proposed method using 2054 naturally blurred images across six imaging applications and seven state-of-the-art deconvolution methods.


title: How Not To Do Mean-Variance Analysis

abstract: We use the 2014 market history of two high-returning biotechnology exchange-traded funds to illustrate how ex post mean-variance analysis should not be done. Unfortunately, the way it should not be done is the way it generally is done -- to our knowledge.


title: Speaker Selective Beamformer with Keyword Mask Estimation

abstract: This paper addresses the problem of automatic speech recognition (ASR) of a target speaker in background speech. The novelty of our approach is that we focus on a wakeup keyword, which is usually used for activating ASR systems like smart speakers. The proposed method firstly utilizes a DNN-based mask estimator to separate the mixture signal into the keyword signal uttered by the target speaker and the remaining background speech. Then the separated signals are used for calculating a beamforming filter to enhance the subsequent utterances from the target speaker. Experimental evaluations show that the trained DNN-based mask can selectively separate the keyword and background speech from the mixture signal. The effectiveness of the proposed method is also verified with Japanese ASR experiments, and we confirm that the character error rates are significantly improved by the proposed method for both simulated and real recorded test sets.


title: Demographic stochasticity and extinction in populations with Allee effect

abstract: We study simple stochastic scenarios, based on birth-and-death Markovian processes, that describe populations with Allee effect, to account for the role of demographic stochasticity. In the mean-field deterministic limit we recover well-known deterministic evolution equations widely employed in population ecology. The mean-time to extinction is in general obtained by the Wentzel-Kramers-Brillouin (WKB) approximation for populations with strong and weak Allee effects. An exact solution for the mean time to extinction can be found via a recursive equation for special cases of the stochastic dynamics. We study the conditions for the validity of the WKB solution and analyze the boundary between the weak and strong Allee effect by comparing exact solutions with numerical simulations.


title: Non-Hermitian quantum systems and their geometric phases

abstract: We discuss the basic theoretical framework for non-Hermitian quantum systems with particular emphasis on the diagonalizability of non-Hermitian Hamiltonians and their $GL(1,\mathbb{C})$ gauge freedom, which are relevant to the adiabatic evolution of non-Hermitian quantum systems. We find that the adiabatic evolution is possible only when the eigen-energies are real. The accompanying geometric phase is found to be generally complex and associated with not only the phase of a wavefunction but also its amplitude. The condition for the real geometric phase is laid out. Our results are illustrated with two examples of non-Hermitian $\mathcal{PT}$ symmetric systems, the two-dimensional non-Hermitian Dirac fermion model and bosonic Bogoliubov quasi-particles.


title: Online adaptive basis enrichment for mixed CEM-GMsFEM

abstract: In this research, an online basis enrichment strategy for the constraint energy minimizing generalized multiscale finite element method in mixed formulation is proposed. The online approach is based on the technique of oversampling. One makes use of the information of residual and the data in the partial differential equation such as the source function. The analysis presented shows that the proposed online enrichment leads to a fast convergence from multiscale approximation to the fine-scale solution. The error reduction can be made sufficiently large by suitably selecting oversampling regions and the number of oversampling layers. Also, the convergence rate of the enrichment can be tuned by a user-defined parameter. Numerical results are provided to illustrate the efficiency of the proposed method.


title: Law and Adversarial Machine Learning

abstract: When machine learning systems fail because of adversarial manipulation, how should society expect the law to respond? Through scenarios grounded in adversarial ML literature, we explore how some aspects of computer crime, copyright, and tort law interface with perturbation, poisoning, model stealing and model inversion attacks to show how some attacks are more likely to result in liability than others. We end with a call for action to ML researchers to invest in transparent benchmarks of attacks and defenses; architect ML systems with forensics in mind and finally, think more about adversarial machine learning in the context of civil liberties. The paper is targeted towards ML researchers who have no legal background.


title: Detection of Coronal Magnetic Activity in Nearby Active Supermassive Black Holes

abstract: Central supermassive black holes of active galactic nuclei host hot plasma with a temperature of $10^9$ K, namely coronae. Like the Sun, black hole coronae are theoretically believed to be heated by their magnetic activity, which have never been observed yet. Here we report the detection of coronal radio synchrotron emission from two nearby Seyfert galaxies using the Atacama Large Millimeter/submillimeter Array, the Karl G. Jansky Very Large Array, and Australia Telescope Compact Array. The coronal magnetic field of both systems is estimated to be $\sim10$ Gauss on scales of $\sim40$ Schwarzschild radii from the central black holes. This magnetic field strength is weaker than the prediction from the magnetically heated corona scenario. We also find that coronae of Seyferts are composed of both thermal and non-thermal electrons. This may imply a possible contribution of Seyferts to the cosmic MeV gamma-ray background radiation.


title: Cicero: Multi-Turn, Contextual Argumentation for Accurate Crowdsourcing

abstract: Traditional approaches for ensuring high quality crowdwork have failed to achieve high-accuracy on difficult problems. Aggregating redundant answers often fails on the hardest problems when the majority is confused. Argumentation has been shown to be effective in mitigating these drawbacks. However, existing argumentation systems only support limited interactions and show workers general justifications, not context-specific arguments targeted to their reasoning. This paper presents Cicero, a new workflow that improves crowd accuracy on difficult tasks by engaging workers in multi-turn, contextual discussions through real-time, synchronous argumentation. Our experiments show that compared to previous argumentation systems which only improve the average individual worker accuracy by 6.8 percentage points on the Relation Extraction domain, our workflow achieves 16.7 percentage point improvement. Furthermore, previous argumentation approaches don't apply to tasks with many possible answers; in contrast, Cicero works well in these cases, raising accuracy from 66.7% to 98.8% on the Codenames domain.


title: OFDM based Sparse Time Dispersive Channel Estimation with Additional Spectral Knowledge

abstract: A new model for sparse time dispersive channels in pilot aided OFDM systems is developed by considering prior knowledge on channel time dispersions. Weighted atomic norm minimization is implemented in the model which enables a more accurate channel estimation. The channel response is identified by solving a Least Squares problem. In this work, we assume that time dispersions' associated frequencies can take any value with a minimum distance on the normalized interval $[0,1)$. The performance of the new model is compared with conventional approaches. With respect to pilot number and SNR, the simulation results reveal that the new model performs superior to that of conventional methods. It is shown that both a lower energy and pilot number are required to achieve the same symbol error rate (SER) reported in previous works.


title: Numerical approximation of optimal convex shapes

abstract: This article investigates the numerical approximation of shape optimization problems with PDE constraint on classes of convex domains. The convexity constraint provides a compactness property which implies well posedness of the problem. Moreover, we prove the convergence of discretizations in two-dimensional situations. A numerical algorithm is devised that iteratively solves the discrete formulation. Numerical experiments show that optimal convex shapes are generally non-smooth and that three-dimensional problems require an appropriate relaxation of the convexity condition.


title: Path-shortening realizations of nonadiabatic holonomic gates

abstract: Nonadiabatic holonomic quantum computation uses non-Abelian geometric phases to implement a universal set of quantum gates that are robust against control imperfections and decoherence. Until now, a number of three-level-based schemes of nonadiabatic holonomic computation have been put forward, and several of them have been experimentally realized. However, all these works are based on the same class of nonadiabatic paths, which originates from the first nonadiabatic holonomic proposal. Here, we propose a universal set of nonadiabatic holonomic gates based on an extended class of nonadiabatic paths. We find that nonadiabatic holonomic gates can be realized with paths shorter than the known ones, which provides the possibility of realizing nonadiabatic holonomic gates with less exposure to decoherence. Furthermore, inspired by the form of this new type of paths, we find a way to eliminate decoherence from nonadiabatic holonomic gates without resorting to redundancies.


title: Modified averaged vector field methods preserving multiple invariants for conservative stochastic differential equations

abstract: A novel class of conservative numerical methods for general conservative Stratonovich stochastic differential equations with multiple invariants is proposed and analyzed. These methods, which are called modified averaged vector field methods, are constructed by modifying the averaged vector field methods to preserve multiple invariants simultaneously. Based on the prior estimate for high order moments of the modification coefficient, the mean square convergence order $1$ of proposed methods is proved in the case of commutative noises. In addition, the effect of quadrature formula on the mean square convergence order and the preservation of invariants for the modified averaged vector field methods is considered. Numerical experiments are performed to verify the theoretical analyses and to show the superiority of the proposed methods in long time simulation.


title: Batch-Parallel Euler Tour Trees

abstract: The dynamic trees problem is to maintain a forest undergoing edge insertions and deletions while supporting queries for information such as connectivity. There are many existing data structures for this problem, but few of them are capable of exploiting parallelism in the batch-setting, in which large batches of edges are inserted or deleted from the forest at once. In this paper, we demonstrate that the Euler tour tree, an existing sequential dynamic trees data structure, can be parallelized in the batch setting. For a batch of $k$ updates over a forest of $n$ vertices, our parallel Euler tour trees perform $O(k \log (1 + n/k))$ expected work with $O(\log n)$ depth with high probability. Our work bound is asymptotically optimal, and we improve on the depth bound achieved by Acar et al. for the batch-parallel dynamic trees problem. The main building block for parallelizing Euler tour trees is a batch-parallel skip list data structure, which we believe may be of independent interest. Euler tour trees require a sequence data structure capable of joins and splits. Sequentially, balanced binary trees are used, but they are difficult to join or split in parallel. We show that skip lists, on the other hand, support batches of joins or splits of size $k$ over $n$ elements with $O(k \log (1 + n/k))$ work in expectation and $O(\log n)$ depth with high probability. We also achieve the same efficiency bounds for augmented skip lists, which allows us to augment our Euler tour trees to support subtree queries. Our data structures achieve between 67--96x self-relative speedup on 72 cores with hyper-threading on large batch sizes. Our data structures also outperform the fastest existing sequential dynamic trees data structures empirically.


title: On the squeezing function and Fridman invariants

abstract: For a domain $D \subset \mathbb C^n$, the relationship between the squeezing function and the Fridman invariants is clarified. Furthermore, localization properties of these functions are obtained. As applications, some known results concerning their boundary behavior are extended.


title: Solidification dynamic of an impacted drop

abstract: This paper is dedicated to the solidification of a water drop impacting a cold solid surface. In a first part, we establish a 1D solidification model, derived from the Stefan problem, that aims at predicting the freezing dynamic of a liquid on a cold substrate, taking into account the thermal properties of this substrate. This model is then experimentally validated through a 1D solidification setup, using different liquids and substrates. In a second part, we show that during the actual drop spreading, a thin layer of ice develops between the water and the substrate, and pins the contact line at its edge when the drop reaches its maximal diameter. The liquid film then remains still on its ice and keeps freezing. This configuration lasts until the contact line eventually depins and the liquid film retracts on the ice. We measure and interpret this crucial time of freezing during which the main ice layer is built. Finally, we compare our 1D model prediction to the thickness of this ice pancake and we find a very good agreement. This allows us to provide a general expression for the frozen drop main thickness, using the drop impact and liquid parameters.


title: All-optical storage of phase-sensitive quantum states of light

abstract: We experimentally demonstrate storage and on-demand release of phase-sensitive, photon-number superposition states of the form $\alpha |0\rangle + \beta e^{i\theta} |1\rangle$ for an optical quantized oscillator mode. For this purpose, we introduce a phase-probing mechanism to a storage system composed of two concatenated optical cavities, which was previously employed for storage of phase-insensitive single-photon states [Phys. Rev. X 3, 041028 (2013)]. This is the first demonstration of all-optically storing highly nonclassical and phase-sensitive quantum states of light. The strong nonclassicality of the states after storage becomes manifest as a negative region in the corresponding Wigner function shifted away from the origin in phase space. This negativity is otherwise, without the phase information of the memory system, unobtainable. While our scheme includes the possibility of optical storage, on-demand release and synchronization of arbitrary single-rail qubit states, it is not limited to such states. In fact, our technique is extendible to more general phase-sensitive states such as multiphoton superposition or entangled states, and thus it represents a significant step toward advanced optical quantum information processing, where highly non-classical states are utilized as resources.


title: Anomalous time-scaling of extreme events in infinite systems and Birkhoff sums of infinite observables

abstract: We establish quantitative results for the statistical be\-ha\-vi\-our of \emph{infinite systems}. We consider two kinds of infinite system: i) a conservative dynamical system $(f,X,\mu)$ preserving a $\sigma$-finite measure $\mu$ such that $\mu(X)=\infty$; ii) the case where $\mu$ is a probability measure but we consider the statistical behaviour of an observable $\phi\colon X\to[0,\infty)$ which is non-integrable: $\int \phi \, d\mu=\infty$. In the first part of this work we study the behaviour of Birkhoff sums of systems of the kind ii). For certain weakly chaotic systems, we show that these sums can be strongly oscillating. However, if the system has superpolynomial decay of correlations or has a Markov structure, then we show this oscillation cannot happen. In this case we prove asymptotic relations between the behaviour of $\phi $, the local dimension of $\mu$, and on the growth of Birkhoff sums (as time tends to infinity). We then establish several important consequences which apply to infinite systems of the kind i). This includes showing anomalous scalings in extreme event limit laws, or entrance time statistics. We apply our findings to non-uniformly hyperbolic systems preserving an infinite measure, establishing anomalous scalings in the case of logarithm laws of entrance times, dynamical Borel--Cantelli lemmas, almost sure growth rates of extremes, and dynamical run length functions.


title: Wearable Affective Robot

abstract: With the development of the artificial intelligence (AI), the AI applications have influenced and changed people's daily life greatly. Here, a wearable affective robot that integrates the affective robot, social robot, brain wearable, and wearable 2.0 is proposed for the first time. The proposed wearable affective robot is intended for a wide population, and we believe that it can improve the human health on the spirit level, meeting the fashion requirements at the same time. In this paper, the architecture and design of an innovative wearable affective robot, which is dubbed as Fitbot, are introduced in terms of hardware and algorithm's perspectives. In addition, the important functional component of the robot-brain wearable device is introduced from the aspect of the hardware design, EEG data acquisition and analysis, user behavior perception, and algorithm deployment, etc. Then, the EEG based cognition of user's behavior is realized. Through the continuous acquisition of the in-depth, in-breadth data, the Fitbot we present can gradually enrich user's life modeling and enable the wearable robot to recognize user's intention and further understand the behavioral motivation behind the user's emotion. The learning algorithm for the life modeling embedded in Fitbot can achieve better user's experience of affective social interaction. Finally, the application service scenarios and some challenging issues of a wearable affective robot are discussed.


title: Simultaneous imaging of voltage and current density of flowing electrons in two dimensions

abstract: Electron transport in nanoscale devices can often result in nontrivial spatial patterns of voltage and current that reflect a variety of physical phenomena, particularly in nonlocal transport regimes. While numerous techniques have been devised to image electron flows, the need remains for a nanoscale probe capable of simultaneously imaging current and voltage distributions with high sensitivity and minimal invasiveness, in magnetic field, across a broad range of temperatures, and beneath an insulating surface. Here we present such a technique for spatially mapping electron flows based on a nanotube single-electron transistor, which achieves high sensitivity for both voltage and current imaging. In a series of experiments using high-mobility graphene devices, we demonstrate the ability of our technique to visualize local aspects of intrinsically nonlocal transport, as in ballistic flows, which are not easily resolvable via existing methods. This technique should both aid in understanding the physics of two-dimensional electronic devices, as well as enable new classes of experiments that image electron flow through buried nanostructures in the quantum and interaction-dominated regimes.


title: Gaussian Message Passing for Overloaded Massive MIMO-NOMA

abstract: This paper considers a low-complexity Gaussian Message Passing (GMP) scheme for a coded massive Multiple-Input Multiple-Output (MIMO) systems with Non-Orthogonal Multiple Access (massive MIMO-NOMA), in which a base station with $N_s$ antennas serves $N_u$ sources simultaneously in the same frequency. Both $N_u$ and $N_s$ are large numbers, and we consider the overloaded cases with $N_u>N_s$. The GMP for MIMO-NOMA is a message passing algorithm operating on a fully-connected loopy factor graph, which is well understood to fail to converge due to the correlation problem. In this paper, we utilize the large-scale property of the system to simplify the convergence analysis of the GMP under the overloaded condition. First, we prove that the \emph{variances} of the GMP definitely converge to the mean square error (MSE) of Linear Minimum Mean Square Error (LMMSE) multi-user detection. Secondly, the \emph{means} of the traditional GMP will fail to converge when $ N_u/N_s< (\sqrt{2}-1)^{-2}\approx5.83$. Therefore, we propose and derive a new convergent GMP called scale-and-add GMP (SA-GMP), which always converges to the LMMSE multi-user detection performance for any $N_u/N_s>1$, and show that it has a faster convergence speed than the traditional GMP with the same complexity. Finally, numerical results are provided to verify the validity and accuracy of the theoretical results presented.


title: Achieving Efficient and Secure Data Acquisition for Cloud-supported Internet of Things in Smart Grid

abstract: Cloud-supported Internet of Things (Cloud-IoT) has been broadly deployed in smart grid systems. The IoT front-ends are responsible for data acquisition and status supervision, while the substantial amount of data is stored and managed in the cloud server. Achieving data security and system efficiency in the data acquisition and transmission process are of great significance and challenging, because the power grid-related data is sensitive and in huge amount. In this paper, we present an efficient and secure data acquisition scheme based on CP-ABE (Ciphertext Policy Attribute Based Encryption). Data acquired from the terminals will be partitioned into blocks and encrypted with its corresponding access sub-tree in sequence, thereby the data encryption and data transmission can be processed in parallel. Furthermore, we protect the information about the access tree with threshold secret sharing method, which can preserve the data privacy and integrity from users with the unauthorized sets of attributes. The formal analysis demonstrates that the proposed scheme can fulfill the security requirements of the Cloud-supported IoT in smart grid. The numerical analysis and experimental results indicate that our scheme can effectively reduce the time cost compared with other popular approaches.


title: Ultrafast hot carrier dynamics of ZrTe$_5$ from time-resolved optical reflectivity

abstract: We investigate the hot carrier dynamics of ZrTe$_5$ by ultrafast time-resolved optical reflectivity. Our results reveal a phonon-mediated across-gap recombination, consistent with its temperature-dependent gap nature as observed previously by photoemission. In addition, two distinct relaxations with a kink feature right after initial photoexcitation are well resolved, suggesting the complexity of electron thermalization process. Our findings indicate that correlated many-body effects play important role for the transient dynamics of ZrTe$_5$.


title: Towards Delay-Tolerant Flexible Data Access Control for Smart Grid with Renewable Energy Resources

abstract: In the Smart Grid with Renewable Energy Resources (RERs), the Residential Units (RUs) with Distributed Energy Resources (DERs) are considered to be both power consumers and suppliers. Specifically, RUs with excessive renewable generations can trade with the utility in deficit of power supplies for mutual benefits. It causes two challenging issues. First, the trading data of RUs is quite sensitive, which should be only accessed by authorized users with fine-grained policies. Second, the behaviors of the RUs to generate trading data are spontaneous and unpredictable, then the problem is how to guarantee system efficiency and delay tolerance simultaneously. In this paper, we propose a delay-tolerant flexible data access control scheme based on Key Policy Attribute Based Encryption (KP-ABE) for Smart Grid with Renewable Energy Resources (RERs). We adopt the secret sharing scheme (SSS) to realize a flexible access control with encryption delay tolerance. Furthermore, there is no central trusted server to perform the encryption/decryption. We reduce the computation cost on RUs and operators via a semi-trusted model. The analysis shows that the proposed scheme can meet the data security requirement of the Smart Grid with RERs, and it also has less cost compared with other popular models.


title: The surface diffusion flow with elasticity in three dimensions

abstract: We establish short-time existence of a smooth solution to the surface diffusion equation with an elastic term and without an additional curvature regularization in three space dimensions. We also prove the asymptotic stability of strictly stable stationary sets.


title: An Efficient Encryption Scheme with Verifiable Outsourced Decryption in Mobile Cloud Computing

abstract: With the increasing number of mobile applications and the popularity of cloud computing, the combination of these two techniques that named mobile cloud computing (MCC) attracts great attention in recent years. A promising public key encryption scheme, Attribute-Based Encryption (ABE), especially the Ciphertext Policy Attribute-Based Encryption (CP-ABE), has been used for realizing fine-grained access control on encrypted data stored in MCC. However, the computational overhead of encryption and decryption grow with the complexity of the access policy. Thus, maintaining data security as well as efficiency of data processing in MCC are important and challenging issues. In this paper, we propose an efficient encryption method based on CP-ABE, which can lower the overhead on data owners. To further reduce the decryption overhead on data receivers, we additionally propose a verifiable outsourced decryption scheme. By security analysis and performance evaluation, the proposed scheme is proved to be secure as well as efficient.


title: Attack Graph Convolutional Networks by Adding Fake Nodes

abstract: In this paper, we study the robustness of graph convolutional networks (GCNs). Previous work have shown that GCNs are vulnerable to adversarial perturbation on adjacency or feature matrices of existing nodes; however, such attacks are usually unrealistic in real applications. For instance, in social network applications, the attacker will need to hack into either the client or server to change existing links or features. In this paper, we propose a new type of "fake node attacks" to attack GCNs by adding malicious fake nodes. This is much more realistic than previous attacks; in social network applications, the attacker only needs to register a set of fake accounts and link to existing ones. To conduct fake node attacks, a greedy algorithm is proposed to generate edges of malicious nodes and their corresponding features aiming to minimize the classification accuracy on the target nodes. In addition, we introduce a discriminator to classify malicious nodes from real nodes, and propose a Greedy-GAN attack to simultaneously update the discriminator and the attacker, to make malicious nodes indistinguishable from the real ones. Our non-targeted attack decreases the accuracy of GCN down to 0.03, and our targeted attack reaches a success rate of 78% on a group of 100 nodes, and 90% on average for attacking a single target node.


title: Word Embedding based Edit Distance

abstract: Text similarity calculation is a fundamental problem in natural language processing and related fields. In recent years, deep neural networks have been developed to perform the task and high performances have been achieved. The neural networks are usually trained with labeled data in supervised learning, and creation of labeled data is usually very costly. In this short paper, we address unsupervised learning for text similarity calculation. We propose a new method called Word Embedding based Edit Distance (WED), which incorporates word embedding into edit distance. Experiments on three benchmark datasets show WED outperforms state-of-the-art unsupervised methods including edit distance, TF-IDF based cosine, word embedding based cosine, Jaccard index, etc.


title: Neutrino phenomenology in the flavored NMSSM without domain wall problems

abstract: We propose a next-to-minimal supersymmetric Standard Model (NMSSM) extended by an $\mathbb{A}_{4}\times \boldsymbol{Z}_{3}$ flavor symmetry and three right-handed neutrinos providing a detailed study of the neutrino sector and a solution to the domain wall problem. In this proposal, neutrino masses are generated through Type I seesaw mechanism while the mixing angles are described by the trimaximal mixing realized using the NMSSM singlet S and only two flavon fields. The phenomenology of neutrino parameters is studied for normal and inverted mass hierarchies. In particular, we numerically evaluated the observables related to neutrino masses and mixing, namely, $\sum m_{i}$, $m_{ee}$, $m_{\nu_{e}}$, and $\delta_{CP}$ where we find that the ranges of $m_{ee}$ and $m_{\nu_{e}}$ are accessible by current and future experiments while the obtained ranges of $\sum m_{i}$ and $\delta_{CP}$ lie within the current experimental data. Another attractive feature we discussed in this paper is the circumvention of the domain wall problem induced by the spontaneous breaking of the $\mathbb{A}_{4}\times \boldsymbol{Z}_{3}$ discrete symmetry. We first showed that the domain walls in the charged lepton sector occur at high energy scale leading to unproblematic domain walls, while in the neutrino sector they are inevitable. Then, to solve this problem, we reconsidered the well-known approach that relies on the explicit breaking of the discrete symmetry through the insertion of Planck-suppressed operators induced by supergravity. \keywords{Neutrino physics, Discrete flavor symmetry,Trimaximal mixing, Domain walls}


title: Hubble drift in Palatini $f(\mathcal{R})$-theories

abstract: In a Palatini $f(\mathcal{R})$-model, we define chonodynamical effects due to the choice of atomic clocks as standard reference clocks and we develop a formalism able to quantitatively separate them from the usual effective dark sources one has in extended theories. We apply the formalism to Hubble drift and briefly discuss the issue about the physical frame. In particular, we argue that there is no physical frame in the sense one does different things in different frames and that, in a sense, is the physical characteristic of extended gravity. As an example, we discuss how Jordan frame may be well suited to discuss cosmology, though it fails within the solar system.


title: Fault Diagnosis and Bad Data Detection of Power Transmission Network - A Time Domain Approach

abstract: Fault analysis and bad data are often processed in separate manners. In this paper it is proved that fault as well as bad current measurement data can be modeled as control failure for the power transmission network and any fault on the transmission line can be treated as multiple bad data. Subsequently a linear observer theory is designed in order to identify the fault type and bad data simultaneously. The state space model based observer theory allows a particular failure mode manifest itself as residual which remains in a fixed direction. Moreover coordinate transformation is performed to allow the residual for each failure mode to generate specific geometry characteristic in separate output dimensions. The design approach based on the observer theory is presented in this paper. The design allows 1) bad data detection for current measurement, and 2) fault location, and fault resistance estimation (as a byproduct) where the fault location accuracy is not affected by fault resistance. However it loses freedom in designing the eigenvalues in the excessive subspace. While the theoretical framework is general, the analysis and design are dedicated to transmission lines.


title: Asymptotic models for free boundary flow in porous media

abstract: We provide rigorous asymptotic models for the free boundary Darcy and Forchheimer problem under the assumption of weak nonlinear interaction, in a regime in which the steepness parameter of the interface is considered to be very small. The models we derive capture the nonlinear interaction of the original free boundary Darcy and Forchheimer problem up to quadratic terms. Furthermore, we provide models that consider both the two-dimensional and three-dimensional cases, with and without bottom topography.


title: A novel bidding method for combined heat and power units in district heating systems

abstract: We propose a bidding method for the participation of combined heat and power (CHP) units in the day-ahead electricity market. More specifically, we consider a district heating system where heat can be produced by CHP units or heat-only units, e.g., gas or wood chip boilers. We use a mixed-integer linear program to determine the optimal operation of the portfolio of production units and storages on a daily basis. Based on the optimal production of subsets of units, we can derive the bidding prices and amounts of electricity offered by the CHP units for the day-ahead market. The novelty about our approach is that the prices are derived by iteratively replacing the production of heat-only units through CHP production. This results in an algorithm with a robust bidding strategy that does not increase the system costs even if the bids are not won. We analyze our method on a small realistic test case to illustrate our method and compare it with other bidding strategies from literature, which consider CHP units individually. The analysis shows that considering a portfolio of units in a district heating system and determining bids based on replacement of heat production of other units leads to better results.


title: Catching the Birth of a Dark Molecular Cloud for the First Time

abstract: The majority of hydrogen in the interstellar medium (ISM) is in atomic form. The transition from atoms to molecules and, in particular, the formation of the H$_2$ molecule, is a key step in cosmic structure formation en route to stars. Quantifying H$_2$ formation in space is difficult, due to the confusion in the emission of atomic hydrogen (HI) and the lack of a H$_2$ signal from the cold ISM. Here we present the discovery of a rare, isolated dark cloud currently undergoing H$_2$ formation, as evidenced by a prominent "ring" of HI self-absorption. Through a combined analysis of HI narrow self-absorption, CO emission, dust emission, and extinction, we directly measured, for the first time, the [HI]/[H$_2$] abundance varying from 2% to 0.2%, within one region. These measured HI abundances are orders of magnitude higher than usually assumed initial conditions for protoplanetary disk models. None of the fast cloud formation model could produce such low atomic hydrogen abundance. We derived a cloud formation timescale of 6$\times$10$^6$ years, consistent with the global Galactic star formation rate, and favoring the classical star formation picture over fast star formation models. Our measurements also help constrain the H$_2$ formation rate, under various ISM conditions.


title: Note on (active-)QRAM-style data access as a quantum circuit

abstract: We observe how an active (i.e., requring $2^n$ parallel control operations) QRAM-like effect $$\sum_{y=0}^{N-1} |y\rangle\langle y| \otimes U^y_{\text{result},\text{memory}_y}$$ can be realized, as a quantum circuit of depth $O(n+\sqrt m)$ (where $m$ is the size of the result register) plus the maximum over all~$z$ of the circuit depths of controlled-$U^z$ operations.


title: Quenched normal approximation for random sequences of transformations

abstract: We study random compositions of transformations having certain uniform fiberwise properties and prove bounds which in combination with other results yield a quenched central limit theorem equipped with a convergence rate, also in the multivariate case, assuming fiberwise centering. For the most part we work with non-stationary randomness and non-invariant, non-product measures. Independently, we believe our work sheds light on the mechanisms that make quenched central limit theorems work, by dissecting the problem into three separate parts.


title: Low-Power Wide-Area Networks for Sustainable IoT

abstract: Low-power wide-area (LPWA) networks are attracting extensive attention because of their abilities to offer low-cost and massive connectivity to Internet of Things (IoT) devices distributed over wide geographical areas. This article provides a brief overview on the existing LPWA technologies and useful insights to aid the large-scale deployment of LPWA networks. Particularly, we first review the currently competing candidates of LPWA networks, such as narrowband IoT (NB-IoT) and long range (LoRa), in terms of technical fundamentals and large-scale deployment potential. Then we present two implementation examples on LPWA networks. By analyzing the field-test results, we identify several challenges that prevent LPWA technologies moving from the theory to wide-spread practice.


title: Chromospheric observations and magnetic configuration of a supergranular structure

abstract: Unipolar magnetic regions are often associated with supergranular cells. The chromosphere above these regions is regulated by the magnetic field, but the field structure is poorly known. In unipolar regions, the fibrillar arrangement does not always coincide with magnetic field lines, and polarimetric observations are needed to establish the chromospheric magnetic topology. In an active region close to the limb, we observed a unipolar annular network of supergranular size. This supergranular structure harbours a radial distribution of the fibrils converging towards its centre. We observed this structure at different heights by taking data in the FeI 6301-6302 {\AA}, H-\alpha, CaII 8542 \AA\ and the CaII H&K spectral lines with the CRISP and CHROMIS instruments at the Swedish 1-m Solar Telescope. We performed Milne-Eddington inversions of the spectropolarimetric data of FeI and applied the weak field approximation to CaII 8542 \AA\ data to retrieve the magnetic field in the photosphere and chromosphere. We used magnetograms of CRISP, HINODE/SP and HMI to calculate the magnetic flux. We investigated the velocity using the line-of-sight velocities computed from the Milne-Eddington inversion and from Doppler shift of the K$_3$ feature in the CaII K spectral line. To describe the typical spectral profiles characterising the chromosphere above the supergranular structure, we performed a K-mean clustering of the spectra in CaIIK. The photospheric magnetic flux is not balanced. The magnetic field vector at chromospheric heights, retrieved by the weak field approximation, indicates that the field lines within the supegranular cell tend to point inwards, and might form a canopy above the unipolar region. In the centre of the supergranular cell hosting the unipolar region, we observe a persistent chromospheric brightening coinciding with a strong gradient in the line-of-sight velocity.


title: Cheeger estimates of Dirichlet-to-Neumann operators on infinite subgraphs of graphs

abstract: In this paper, we study the Dirichlet-to-Neumann operators on infinite subgraphs of graphs. For an infinite graph, we prove Cheeger-type estimates for the bottom spectrum of the Dirichlet-to-Neumann operator, and the higher order Cheeger estimates for higher order eigenvalues of the Dirichlet-to-Neumann operator.


title: A two-phase stochastic programming approach to biomass supply planning for combined heat and power plants

abstract: Due to the new carbon neutral policies, many district heating operators start operating their combined heat and power (CHP) plants using different types of biomass instead of fossil fuel. The contracts with the biomass suppliers are negotiated months in advance and involve many uncertainties from the energy producer's side. The demand for biomass is uncertain at that time, and heat demand and electricity prices vary drastically during the planning period. Furthermore, the optimal operation of combined heat and power plants has to consider the existing synergies between the power and heating systems. We propose a solution method using stochastic optimization to support the biomass supply planning for combined heat and power plants. Our two-phase approach determines mid-term decisions about biomass supply contracts as well as short-term decisions regarding the optimal production of the producer to ensure profitability and feasibility. We present results based on two realistic test cases.


title: Unique continuation and classification of blow-up profiles for elliptic systems with Neumann boundary coupling and applications to higher order fractional equations

abstract: In this paper we develop a monotonicity formula for elliptic systems with Neumann boundary coupling, proving unique continuation and classification of blow-up profiles. As an application, we obtain strong unique continuation for some fourth order equations and higher order fractional problems.


title: The probability distribution of Brownian motion in periodic potentials

abstract: We calculate the probability distribution function (PDF) of an overdamped Brownian particle moving in a periodic potential energy landscape $U(x)$. The PDF is found by solving the corresponding Smoluchowski diffusion equation. We derive the solution for any periodic even function $U(x)$, and demonstrate that it is asymptotically (at large time $t$) correct up to terms decaying faster than $\sim t^{-3/2}$. As part of the derivation, we also recover the Lifson-Jackson formula for the effective diffusion coefficient of the dynamics. The derived solution exhibits agreement with Langevin dynamics simulations when (i) the periodic length is much larger than the ballistic length of the dynamics, and (ii) when the potential barrier $\Delta U=\max(U(x))-\min(U(x))$ is not much larger than the thermal energy $k_BT$.


title: Quasianalytic ultradifferentiability cannot be tested in lower dimensions

abstract: We show that, in contrast to the real analytic case, quasianalytic ultradifferentiability can never be tested in lower dimensions. Our results are based on a construction due to Jaffe.


title: Realizing negative Poisson's ratio in unstressed spring networks

abstract: When randomly displacing the nodes of a crystalline and unstressed spring network, we find that the Possion's ratio decreases with the increase of structural disorder and even becomes negative. Employing our finding that longer springs tend to contribute more to the shear modulus but less to the bulk modulus, we are able to achieve negative Poisson's ratio with lower structural disorder by attributing each spring a length dependent stiffness. Even with perfect crystalline structure, the network can have negative Possion's ratio, if the stiffness of each spring is set by its virtual length after a virtual network distortion. We also reveal that the nonaffine contribution arising from the structural or spring constant disorder produced in some cooperative way by network distortion is essential to the emergence of negative Poisson's ratio.


title: Expedition: A Time-Aware Exploratory Search System Designed for Scholars

abstract: Archives are an important source of study for various scholars. Digitization and the web have made archives more accessible and led to the development of several time-aware exploratory search systems. However these systems have been designed for more general users rather than scholars. Scholars have more complex information needs in comparison to general users. They also require support for corpus creation during their exploration process. In this paper we present Expedition - a time-aware exploratory search system that addresses the requirements and information needs of scholars. Expedition possesses a suite of ad-hoc and diversity based retrieval models to address complex information needs; a newspaper-style user interface to allow for larger textual previews and comparisons; entity filters to more naturally refine a result list and an interactive annotated timeline which can be used to better identify periods of importance.


title: Geometry and clustering with metrics derived from separable Bregman divergences

abstract: Separable Bregman divergences induce Riemannian metric spaces that are isometric to the Euclidean space after monotone embeddings. We investigate fixed rate quantization and its codebook Voronoi diagrams, and report on experimental performances of partition-based, hierarchical, and soft clustering algorithms with respect to these Riemann-Bregman distances.


title: An Incremental Truth Inference Approach to Aggregate Crowdsourcing Contributions in Games with a Purpose

abstract: We introduce our approach for incremental truth inference over the contributions provided by players of Games with a Purpose: we motivate the need for such a method with the specificity of GWAP vs. traditional crowdsourcing; we explain and formalize the proposed process and we explain its positive consequences; finally, we illustrate the results of an experimental comparison with state-of-the-art approaches, performed on data collected through two different GWAPs, thus showing the properties of our proposed framework.


title: Perturbation of the scattering resonances of an open cavity by small particles. Part I: The transverse magnetic polarization case

abstract: This paper aims at providing a small-volume expansion framework for the scattering resonances of an open cavity perturbed by small particles. The induced shift of the scattering frequencies by the small particles is derived without neglecting the radiation effect. The formula holds for arbitrary-shaped particles. It shows a strong enhancement in the frequency shift in the case of plasmonic particles. The formula is used to image small particles located near the boundary of an open resonator which admits whispering-gallery modes. Numerical examples of interest for applications are presented.


title: A fast radio burst with a low dispersion measure

abstract: Fast radio bursts (FRBs) are millisecond pulses of radio emission of seemingly extragalactic origin. More than 50 FRBs have now been detected, with only one seen to repeat. Here we present a new FRB discovery, FRB 110214, which was detected in the high latitude portion of the High Time Resolution Universe South survey at the Parkes telescope. FRB 110214 has one of the lowest dispersion measures of any known FRB (DM = 168.9$\pm$0.5 pc cm$^{-3}$), and was detected in two beams of the Parkes multi-beam receiver. A triangulation of the burst origin on the sky identified three possible regions in the beam pattern where it may have originated, all in sidelobes of the primary detection beam. Depending on the true location of the burst the intrinsic fluence is estimated to fall in the range of 50 -- 2000 Jy ms, making FRB 110214 one of the highest-fluence FRBs detected with the Parkes telescope. No repeating pulses were seen in almost 100 hours of follow-up observations with the Parkes telescope down to a limiting fluence of 0.3 Jy ms for a 2-ms pulse. Similar low-DM, ultra-bright FRBs may be detected in telescope sidelobes in the future, making careful modeling of multi-beam instrument beam patterns of utmost importance for upcoming FRB surveys.


title: Operational planning and bidding for district heating systems with uncertain renewable energy production

abstract: In countries with an extended use of district heating (DH), the integrated operation of DH and power systems can increase the flexibility of the power system achieving a higher integration of renewable energy sources (RES). DH operators can not only provide flexibility to the power system by acting on the electricity market, but also profit from the situation to lower the overall system cost. However, the operational planning and bidding includes several uncertain components at the time of planning: electricity prices as well as heat and power production from RES. In this publication, we propose a planning method that supports DH operators by scheduling the production and creating bids for the day-ahead and balancing electricity markets. The method is based on stochastic programming and extends bidding strategies for virtual power plants to the DH application. The uncertain factors are considered explicitly through scenario generation. We apply our solution approach to a real case study in Denmark and perform an extensive analysis of the production and trading behaviour of the DH system. The analysis provides insights on how DH system can provide regulating power as well as the impact of uncertainties and renewable sources on the planning. Furthermore, the case study shows the benefit in terms of cost reductions from considering a portfolio of units and both markets to adapt to RES production and market states.


title: Adversarially Robust Optimization with Gaussian Processes

abstract: In this paper, we consider the problem of Gaussian process (GP) optimization with an added robustness requirement: The returned point may be perturbed by an adversary, and we require the function value to remain as high as possible even after this perturbation. This problem is motivated by settings in which the underlying functions during optimization and implementation stages are different, or when one is interested in finding an entire region of good inputs rather than only a single point. We show that standard GP optimization algorithms do not exhibit the desired robustness properties, and provide a novel confidence-bound based algorithm StableOpt for this purpose. We rigorously establish the required number of samples for StableOpt to find a near-optimal point, and we complement this guarantee with an algorithm-independent lower bound. We experimentally demonstrate several potential applications of interest using real-world data sets, and we show that StableOpt consistently succeeds in finding a stable maximizer where several baseline methods fail.


title: Effective models of WIMP Direct Detection in DAMA/LIBRA-phase2

abstract: Since the DAMA/LIBRA collaboration released updated results from their search for the annual modulation signal expected from Dark Matter (DM) scattering in their NaI detectors, we have fitted the updated DAMA result for the modulation amplitudes for a Weakly Interacting Massive Particle (WIMP) signal, parameterizing the interaction with nuclei in terms of the most general effective Lagrangian for a WIMP particle spin 1/2, systematically assuming dominance of one of the 14 possible interaction terms and a standard Maxwellian for the WIMP velocity distribution. We find that most of the couplings of the non-relativistic effective Hamiltonian can provide a better fit compared to the standard Spin Independent interaction case, and with a reduced fine-tuning of the three parameters (WIMP mass, WIMP-nucleon effective cross-section and ratio between the WIMP-neutron and the WIMP-proton couplings). In addition, effective models for which the cross section depends explicitly on the WIMP incoming velocity can provide a better fit of the DAMA data at large values of $m_{\chi}$ compared to the standard velocity-independent cross-section due to a different phase of the modulation amplitudes. All the best fit solutions are in tension with exclusion plots of both XENON1T and PICO60.


title: Efficient Learning of Restricted Boltzmann Machines Using Covariance Estimates

abstract: Learning RBMs using standard algorithms such as CD(k) involves gradient descent on the negative log-likelihood. One of the terms in the gradient, which involves expectation w.r.t. the model distribution, is intractable and is obtained through an MCMC estimate. In this work we show that the Hessian of the log-likelihood can be written in terms of covariances of hidden and visible units and hence, all elements of the Hessian can also be estimated using the same MCMC samples with small extra computational costs. Since inverting the Hessian may be computationally expensive, we propose an algorithm that uses inverse of the diagonal approximation of the Hessian, instead. This essentially results in parameter-specific adaptive learning rates for the gradient descent process and improves the efficiency of learning RBMs compared to the standard methods. Specifically we show that using the inverse of diagonal approximation of Hessian in the stochastic DC (difference of convex functions) program approach results in very efficient learning of RBMs.


title: Infinite family of universal profiles for heat current statistics in Fourier's law

abstract: Using tools from large deviation theory, we study fluctuations of the heat current in a model of $d$-dimensional incompressible fluid driven out of equilibrium by a temperature gradient. We find that the most probable temperature fields sustaining atypical values of the global current can be naturally classified in an infinite set of curves, allowing us to exhaustively analyze their topological properties and to define universal profiles onto which all optimal fields collapse. We also compute the statistics of empirical heat current, where we find remarkable logarithmic tails for large current fluctuations orthogonal to the thermal gradient. Finally, we determine explicitly a number of cumulants of the current distribution, finding remarkable relations between them.


title: Discrete variants of the phi4 model: exceptional discretizations, conservation laws and related topics

abstract: Exceptional dicretizations of the phi4 model are reviewed, corresponding conservation laws are reported, and the properties of static and moving discrete kinks are discussed. Different approaches to producing such discretizations are given and unifying perspectives thereof are brought forth. It is also demonstrated that the high kink mobility in the exceptional dicretizations makes it possible to analyze kink-antikink collisions in the regime of high discreteness.


title: Non-standard scaling laws of WIMP-nucleus interactions in WIMP direct detection

abstract: Guided by the non-relativistic effective field theory of interactions between Weakly Interacting Massive Particles (WIMPs) of spin 1/2 and nuclei we study direct detection exclusion plots for an example of non-standard spin-dependent interaction and compare it to the standard one. We analyze an extensive list of 15 existing experiments including the effects of momentum dependence and isospin violation. In our analysis, we fixed the dark matter velocity distribution to a Maxwellian.


title: Limits on stable iron in Type$\,$Ia supernovae from NIR spectroscopy

abstract: We obtained optical and near-infrared spectra of Type$\,$Ia supernovae (SNe$\,$Ia) at epochs ranging from 224 to 496 days after the explosion. The spectra show emission lines from forbidden transitions of singly ionised iron and cobalt atoms. We used non-local thermodynamic equilibrium (NLTE) modelling of the first and second ionisation stages of iron, nickel, and cobalt to fit the spectra using a sampling algorithm allowing us to probe a broad parameter space. We derive velocity shifts, line widths, and abundance ratios for iron and cobalt. The measured line widths and velocity shifts of the singly ionised ions suggest a shared emitting region. Our data are fully compatible with radioactive $^{56}$Ni decay as the origin for cobalt and iron. We compare the measured abundance ratios of iron and cobalt to theoretical predictions of various SN$\,$Ia explosion models. These models include, in addition to $^{56}$Ni, different amounts of $^{57}$Ni and stable $^{54,56}$Fe. We can exclude models that produced only $^{54,56}$Fe or only $^{57}$Ni in addition to $^{56}$Ni. If we consider a model that has $^{56}$Ni, $^{57}$Ni, and $^{54,56}$Fe then our data imply that these ratios are $^{54,56}$Fe / $^{56}$Ni $=0.272\pm0.086$ and $^{57}$Ni / $^{56}$Ni $=0.032\pm0.011$.