title: Laser induced distortion of band structure in solids: an analytic model

abstract: We consider a spatially periodic (cosine) potential as a model for a crystalline solid that interacts with a harmonically oscillating external electric field. This problem is periodic both in space and time and can be solved analytically using the Kramers-Henneberger co-moving frame. By analyzing the stability of the closely related Mathieu-type differential equation, the electronic band structure can be obtained. We demonstrate that by changing the field intensity, the width of the zero-field band gaps can be drastically modified, including the special case when the external field causes the band gaps to disappear


title: Contrastive Training for Models of Information Cascades

abstract: This paper proposes a model of information cascades as directed spanning trees (DSTs) over observed documents. In addition, we propose a contrastive training procedure that exploits partial temporal ordering of node infections in lieu of labeled training links. This combination of model and unsupervised training makes it possible to improve on models that use infection times alone and to exploit arbitrary features of the nodes and of the text content of messages in information cascades. With only basic node and time lag features similar to previous models, the DST model achieves performance with unsupervised training comparable to strong baselines on a blog network inference task. Unsupervised training with additional content features achieves significantly better results, reaching half the accuracy of a fully supervised model.


title: Geometry of epithelial cells provides a robust method for image based inference of stress within tissues

abstract: Cellular mechanics plays an important role in epithelial morphogenesis, a process wherein cells reshape and rearrange to produce tissue-scale deformations. However, the study of tissue-scale mechanics is impaired by the difficulty of direct measurement of stress in-vivo. Alternative, image-based inference schemes aim to estimate stress from snapshots of cellular geometry but are challenged by sensitivity to fluctuations and measurement noise as well as the dependence on boundary conditions. Here we overcome these difficulties by introducing a new variational approach - the Geometrical Variation Method (GVM) - which exploits the fundamental duality between stress and cellular geometry that exists in the state of mechanical equilibrium of discrete mechanical networks that approximate cellular tissues. In the Geometrical Variation Method, the two dimensional apical geometry of an epithelial tissue is approximated by a 2D tiling with Circular Arc Polygons (CAP) in which the arcs represent intercellular interfaces defined by the balance of local line tension and pressure differentials between adjacent cells. We take advantage of local constraints that mechanical equilibrium imposes on CAP geometry to define a variational procedure that extracts the best fitting equilibrium configuration from images of epithelial monolayers. The GVM-based stress inference algorithm has been validated by the comparison of the predicted cellular and mesoscopic scale stress and measured myosin II patterns in the epithelial tissue during Drosophila embryogenesis. GVM prediction of mesoscopic stress tensor correlates at the 80% level with the measured myosin distribution and reveals that most of the myosin II activity is involved in a static internal force balance within the epithelial layer. Lastly, this study provides a practical method for non-destructive estimation of stress in live epithelial tissues.


title: Testing the standard model with $D_{(s)} \to K_1 (\to K\pi\pi) \gamma$ decays

abstract: The photon polarization in $D_{(s)} \to K_1 (\to K\pi\pi) \gamma$ decays can be extracted from an up-down asymmetry in the $K \pi \pi$ system, along the lines of the method known to $B \to K_1 (\to K\pi\pi) \gamma$ decays. Charm physics is advantageous as partner decays exist: $D^+ \to K_1^+ (\to K\pi\pi) \gamma$, which is standard model-like, and $D_s \to K_1^+ (\to K\pi\pi) \gamma$, which is sensitive to physics beyond the standard model in $|\Delta c| =|\Delta u|=1$ transitions. The standard model predicts their photon polarizations to be equal up to U-spin breaking corrections, while new physics in the dipole operators can split them apart at order one level. We estimate the proportionality factor in the asymmetry multiplying the polarization parameter from axial vectors $K_1(1270)$ and $K_1(1400)$ to be sizable, up to the few ${\cal{O}}(10)\%$ range. The actual value of the hadronic factor matters for the experimental sensitivity, but is not needed as an input to perform the null test.


title: A Score Based Test for Functional Linear Concurrent Regression

abstract: We propose a novel method for testing the null hypothesis of no effect of a covariate on the response in the context of functional linear concurrent regression. We establish an equivalent random effects formulation of our functional regression model under which our testing problem reduces to testing for zero variance component for random effects. For this purpose, we use a one-sided score test approach, which is an extension of the classical score test. We provide theoretical justification as to why our testing procedure has the right levels (asymptotically) under null using standard assumptions. Using numerical simulations, we show that our testing method has the desired type I error rate and gives higher power compared to a bootstrapped F test currently existing in the literature. Our model and testing procedure are shown to give good performances even when the data is sparsely observed, and the covariate is contaminated with noise. Applications of the proposed testing method are demonstrated on gait study and a dietary calcium absorption data.


title: Rapid adiabatic gating for capacitively coupled quantum dot hybrid qubits without barrier control

abstract: We theoretically examine the capacitive coupling between two quantum dot hybrid qubits, each consisting of three electrons in a double quantum dot, as a function of the energy detuning of the double dot potentials. We show that a shaped detuning pulse can produce a two-qubit maximally entangling operation in $\sim$50ns without having to simultaneously change tunnel couplings. Simulations of the entangling operation in the presence of experimentally realistic charge noise yield two-qubit fidelities over 90%.


title: Femural Head Autosegmentation for 3D Radiotherapy Planning: Preliminary Results

abstract: Contouring of organs at risk is an important but time consuming part of radiotherapy treatment planning. Several authors proposed methods for automatic delineation but the clinical experts eye remains the gold standard method. In this paper, we present a totally visual software for automated delineation of the femural head. The software was successfully characterized in pelvic CT Scan of prostate patients (n=11). The automatic delineation was compared with manual and approved delineation through blind test evaluated by a panel of seniors radiation oncologists (n=9). Clinical experts evaluated that no any contouring correction were need in 77.8% and 67.8% of manual and automatic delineation respectively. Our results show that the software is robust, the automated delineation was reproducible in all patient, and its performance was similar to manually delineation.


title: Recovering the QNEC from the ANEC

abstract: We study relative entropy in QFT, comparing the vacuum state to a special family of purifications determined by an input state and constructed using relative modular flow. We use this to prove a conjecture by Wall that relates the shape derivative of relative entropy to a variational expression over the averaged null energy of possible purifications. This variational expression can be used to easily prove the quantum null energy condition. We formulate Wall's conjecture as a theorem pertaining to operator algebras satisfying the properties of a half-sided modular inclusion, with the additional assumption that the input state has finite averaged null energy. We also give a new derivation of the strong superadditivity property of relative entropy in this context. We speculate about possible connections to the recent methods used to strengthen monotonicity of relative entropy with recovery maps.


title: A new, clean catalogue of extragalactic non-nuclear X-ray sources in nearby galaxies

abstract: We have created a new, clean catalogue of extragalactic non-nuclear X-ray sources by correlating the 3XMM-DR4 data release of the XMM-Newton Serendipitous Source Catalogue with the Third Reference Catalogue of Bright Galaxies and the Catalogue of Neighbouring Galaxies, using an improved version of the method presented in Walton et al. (2011). Our catalogue contains 1,314 sources, of which 384 are candidate ultraluminous X-ray sources (ULXs). The resulting catalogue improves upon previous catalogues in its handling of spurious detections by taking into account XMM-Newton quality flags. We estimate the contamination of ULXs by background sources to be 24 per cent. We define a 'complete' subsample as those ULXs in galaxies for which the sensitivity limit is below $10^{39}$ erg/s and use it to examine the hardness ratio properties between ULX and non-ULX sources, and ULXs in different classes of host galaxy. We find that ULXs have a similar hardness ratio distribution to lower-luminosity sources, consistent with previous studies. We also find that ULXs in spiral and elliptical host galaxies have similar distributions to each other independent of host galaxy morphology, however our results do support previous indications that the population of ULXs is more luminous in star-forming host galaxies than in non-star-forming galaxies. Our catalogue contains further interesting subpopulations for future study, including Eddington Threshold sources and highly variable ULXs. We also examine the highest-luminosity (L$_X$ > $5 \times 10^{40}$ erg/s) ULXs in our catalogue in search of intermediate-mass black hole candidates, and find nine new possible candidates.


title: IceCube Neutrinos from Hadronically Powered Gamma-Ray Galaxies

abstract: In this work we use a multi-messenger approach to determine if the high energy diffuse neutrino flux observed by the IceCube Observatory can originate from $\gamma$-ray sources powered by Cosmic Rays interactions with gas. Typical representatives of such sources are Starburst and Ultra-Luminous Infrared Galaxies. Using the three most recent calculations of the non-blazar contribution to the extragalactic $\gamma$-ray background measured by the Fermi-LAT collaboration, we find that a hard power-law spectrum with spectral index $\alpha \leq 2.12$ is compatible with all the estimations for the allowed contribution from non-blazar sources, within 1$\sigma$. Using such a spectrum we are able to interpret the IceCube results, showing that various classes of hadronically powered $\gamma$-ray galaxies can provide the dominant contribution to the astrophysical signal above 100 TeV and about half of the contribution to the energy flux between 10-100 TeV. With the addition of neutrinos from the Galactic plane, it is possible to saturate the IceCube signal at high energy. Our result shows that these sources are still well motivated candidates.


title: Microwave impedance readout of a hafnium microbridge detector

abstract: We present proof-of-operation for a new method of electron thermometry using microwave impedance of a hafnium micro-absorber. The new method leads to an ultimate THz-range detector suitable for microwave readout and frequency division multiplexing. The sensing part of the device is a hot-electron-gas absorber responding to the incident radiation by variation of its impedance measured at probing frequency about 1.5 GHz. The absorber is a microbridge made from hafnium (Tc = 375 mK, RN = 30 Ohm) sized 2.5 um by 2.5 um by 50 nm and integrated with a planar 600-700 GHz antenna placed near the open end of a quarter-wave CPW resonator (Q-factor about 10^4). All elements of the circuit, except the microbridge, are made from 100-nm thick Nb, including the resonator, which is weakly coupled to a throughput line. The device was tested at 50-350 mK smoothly responding with its transmission coefficient S21 to applied microwave power at the resonance frequency. We have found that the power absorbed by the bridge fits to the model of hot electron gas, P=k(Te^n-Tph^n) (n = 5...6). The idle NEP down to about 10^-18 W/Hz^(-1/2) and the corresponding cross-over temperature for photon background about 5 K are estimated from the measured data. The saturation power of about 1 pW and possibility of moderate gain are anticipated for a practicable device operating at temperature 200 mK. Since the optimum readout frequency is found exactly at the resonance, the detector is insensitive to most phase instabilities at the probing frequency.


title: Enhanced mixing in magnetized fingering convection, and implications for RGB stars

abstract: Double-diffusive convection has been well studied in geophysical contexts, but detailed investigations of the regimes characteristic of stellar or planetary interiors have only recently become feasible. Since most astrophysical fluids are electrically conducting, it is possible that magnetic fields play a role in either enhancing or suppressing double-diffusive convection, but to date there have been no numerical investigations of such possibilities. Here we study the effects of a vertical background magnetic field (aligned with the gravitational axis) on the linear stability and nonlinear saturation of fingering (thermohaline) convection, through a combination of theoretical work and direct numerical simulations (DNSs). We find that a vertical magnetic field rigidifies the fingers along the vertical direction which has the remarkable effect of enhancing vertical mixing. We propose a simple analytical model for mixing by magnetized fingering convection, and argue that magnetic effects may help explain discrepancies between theoretical and observed mixing rates in low-mass red giant branch (RGB) stars. Other implications of our findings are also discussed.


title: The integrated properties of the molecular clouds from the JCMT CO(3-2) High Resolution Survey

abstract: We define the molecular cloud properties of the Milky Way first quadrant using data from the JCMT CO(3-2) High Resolution Survey. We apply the Spectral Clustering for Interstellar Molecular Emission Segmentation (SCIMES) algorithm to extract objects from the full-resolution dataset, creating the first catalog of molecular clouds with a large dynamic range in spatial scale. We identify $>85\,000$ clouds with two clear sub-samples: $\sim35\,500$ well-resolved objects and $\sim540$ clouds with well-defined distance estimations. Only 35% of the cataloged clouds (as well as the total flux encompassed by them) appear enclosed within the Milky Way spiral arms. The scaling relationships between clouds with known distances are comparable to the characteristics of the clouds identified in previous surveys. However, these relations between integrated properties, especially from the full catalog, show a large intrinsic scatter ($\sim0.5$ dex), comparable to other cloud catalogs of the Milky Way and nearby galaxies. The mass distribution of molecular clouds follows a truncated-power law relationship over three orders of magnitude in mass with a form $dN/dM \propto M^{-1.7}$ with a clearly defined truncation at an upper mass of $M_0 \sim 3 \times 10^6~M_\odot$, consistent with theoretical models of cloud formation controlled by stellar feedback and shear. Similarly, the cloud population shows a power-law distribution of size with $dN/dR \propto R^{-2.8}$ with a truncation at $R_0 = 70$ pc.


title: Foliations and Conjugacy II: The Mendes Conjecture for Time-One Maps of Flows

abstract: A diffeomorphism $f:\mathbb{R}^2\to\mathbb{R}^2$ in the plane is Anosov if it has a hyperbolic splitting at every point of the plane. The two known topological conjugacy classes of such diffeomorphisms are linear hyperbolic automorphisms and translations (the existence of Anosov structures for plane translations was originally shown by W. White). P. Mendes conjectured that these are the only topological conjugacy classes for Anosov diffeomorphisms in the plane. We prove that this claim holds when the Anosov diffeomorphism is the time-one map of a flow, via a theorem about foliations invariant under a time one map.


title: Learning representations of molecules and materials with atomistic neural networks

abstract: Deep Learning has been shown to learn efficient representations for structured data such as image, text or audio. In this chapter, we present neural network architectures that are able to learn efficient representations of molecules and materials. In particular, the continuous-filter convolutional network SchNet accurately predicts chemical properties across compositional and configurational space on a variety of datasets. Beyond that, we analyze the obtained representations to find evidence that their spatial and chemical properties agree with chemical intuition.


title: Improved stabilization technique for frictional contact problems solved with hp-BEM

abstract: We improve the residual based stabilization technique for Signorini contact problems with Tresca friction in linear elasticity solved with $hp$-mixed BEM which has been recently analyzed by Banz et al.~in Numer.~Math.~135 (2017) pp.~217--263. The stabilization allows us to circumvent the discrete inf-sup condition and thus the primal and dual sets can be discretized independently. Compared to the above mentioned paper we are able to remove the dependency of the scaling parameter on the unknown Sobolev regularity of the exact solution and can thus also improve the convergence rate in the a priori error estimate. The second improvement is a rigorous a priori and a posteriori error analysis when the boundary integral operators in the stabilization term are approximated. The latter is of fundamental importance to keep the computational time small. We present numerical results in two and three dimensions to underline our theoretical findings, show the superiority of the $hp$-adaptive stabilized mixed scheme and the effect induced by approximating the stabilization term. Moreover, we show the applicability of the proposed method to the Coulomb frictional case for which we extend the a posteriori error analysis


title: A Survey on High-Capacity OFDM-based Passive Optical Networks

abstract: The exponential growth of demand for high-speed internet and high bandwidth applications such as interactive entertainment in access networks mandates the requirement for higher optical signal speeds. On the other hand, since cost and energy efficiency should be concurrently preserved in access networks, PONs have emerged as a breakthrough solution. This survey shows that future-proof PONs should be supported by an adaptively modulated WDM-OFDM architecture to maximize signal capacity and transmission-reach, while the employment of reflective semiconductor optical amplification and re-modulation can potentially prevent the employment of an additional light source.


title: ECG Arrhythmia Classification Using Transfer Learning from 2-Dimensional Deep CNN Features

abstract: Due to the recent advances in the area of deep learning, it has been demonstrated that a deep neural network, trained on a huge amount of data, can recognize cardiac arrhythmias better than cardiologists. Moreover, traditionally feature extraction was considered an integral part of ECG pattern recognition; however, recent findings have shown that deep neural networks can carry out the task of feature extraction directly from the data itself. In order to use deep neural networks for their accuracy and feature extraction, high volume of training data is required, which in the case of independent studies is not pragmatic. To arise to this challenge, in this work, the identification and classification of four ECG patterns are studied from a transfer learning perspective, transferring knowledge learned from the image classification domain to the ECG signal classification domain. It is demonstrated that feature maps learned in a deep neural network trained on great amounts of generic input images can be used as general descriptors for the ECG signal spectrograms and result in features that enable classification of arrhythmias. Overall, an accuracy of 97.23 percent is achieved in classifying near 7000 instances by ten-fold cross validation.


title: Rotation measurements using a a resonant fiber optic gyroscope based on Kagome fiber

abstract: We build a resonant fiber optic gyro based on Kagome hollow-core fiber. A semi-bulk cavity architecture based on a 18-m-long Kagome fiber permits to achieve a cavity finesse of 23 with a resonance linewidth of 700 kHz. An optimized Pound-Drever-Hall servo-locking scheme is used to probe the cavity in reflection. Closed-loop operation of the gyroscope permits to reach an angular random walk as small as 0.004$^\circ/\sqrt{\mathrm{h}}$ and a bias stability of 0.45$^\circ$/h over 0.5 s of integration time.


title: Clebsch-Lagrange variational principle and geometric constraint analysis of relativistic field theories

abstract: Inspired by the Clebsch optimal control problem, we introduce a new variational principle that is suitable for capturing the geometry of relativistic field theories with constraints related to a gauge symmetry. Its special feature is that the Lagrange multipliers couple to the configuration variables via the symmetry group action. The resulting constraints are formulated as a condition on the momentum map of the gauge group action on the phase space of the system. We discuss the Hamiltonian picture and the reduction of the gauge symmetry by stages in this geometric setting. We show that the Yang-Mills-Higgs action and the Einstein--Hilbert action fit into this new framework after a $ (1+3) $-splitting. Moreover, we recover the Gau{\ss} constraint of Yang-Mills-Higgs theory and the diffeomorphism constraint of general relativity as momentum map constraints.


title: Polymorphism of Bulk Boron Nitride

abstract: Boron nitride (BN) is a material with outstanding technological promise because of its exceptional thermochemical stability, structural, electronic and thermal conductivity properties, and extreme hardness. Yet, the relative thermodynamic stability of its most common polymorphs (diamond-like cubic and graphite-like hexagonal) has not been resolved satisfactorily because of the crucial role played by kinetic factors in the formation of BN phases at high temperatures and pressures (experiments), and by competing bonding, electrostatic and many-body dispersion forces in BN cohesion (theory). This lack of understanding hampers the development of potential technological applications, and challenges the boundaries of fundamental science. Here, we use high-level first-principles theories that correctly reproduce all important electronic interactions (the adiabatic-connection fluctuation-dissipation theorem in the random phase approximation) to estimate with unprecedented accuracy the energy differences between BN polymorphs, and thus overcome the accuracy hurdle that hindered previous theoretical studies. We show that the ground-state phase of BN is cubic and that the frequently observed hexagonal polymorph becomes entropically stabilized over the cubic at temperatures slightly above ambient conditions ($T_{\rm c \to h} = 335 \pm 30$ K). We also reveal a low-symmetry monoclinic phase that is extremely competitive with the other low-energy polymorphs and which could explain the origins of the experimentally observed "compressed h-BN" phase. Our theoretical findings therefore should stimulate new experimental efforts in bulk BN as well as promote the use of high-level theories in modelling of technologically relevant van der Waals materials.


title: Anomaly Generation using Generative Adversarial Networks in Host Based Intrusion Detection

abstract: Generative adversarial networks have been able to generate striking results in various domains. This generation capability can be general while the networks gain deep understanding regarding the data distribution. In many domains, this data distribution consists of anomalies and normal data, with the anomalies commonly occurring relatively less, creating datasets that are imbalanced. The capabilities that generative adversarial networks offer can be leveraged to examine these anomalies and help alleviate the challenge that imbalanced datasets propose via creating synthetic anomalies. This anomaly generation can be specifically beneficial in domains that have costly data creation processes as well as inherently imbalanced datasets. One of the domains that fits this description is the host-based intrusion detection domain. In this work, ADFA-LD dataset is chosen as the dataset of interest containing system calls of small foot-print next generation attacks. The data is first converted into images, and then a Cycle-GAN is used to create images of anomalous data from images of normal data. The generated data is combined with the original dataset and is used to train a model to detect anomalies. By doing so, it is shown that the classification results are improved, with the AUC rising from 0.55 to 0.71, and the anomaly detection rate rising from 17.07% to 80.49%. The results are also compared to SMOTE, showing the potential presented by generative adversarial networks in anomaly generation.


title: Slice theorem and orbit type stratification in infinite dimensions

abstract: We establish a general slice theorem for the action of a locally convex Lie group on a locally convex manifold, which generalizes the classical slice theorem of Palais to infinite dimensions. We discuss two important settings under which the assumptions of this theorem are fulfilled. First, using Gl\"ockner's inverse function theorem, we show that the linear action of a compact Lie group on a Fr\'echet space admits a slice. Second, using the Nash--Moser theorem, we establish a slice theorem for the tame action of a tame Fr\'echet Lie group on a tame Fr\'echet manifold. For this purpose, we develop the concept of a graded Riemannian metric, which allows the construction of a path-length metric compatible with the manifold topology and of a local addition. Finally, generalizing a classical result in finite dimensions, we prove that the existence of a slice implies that the decomposition of the manifold into orbit types of the group action is a stratification.


title: Stability boundaries of a Mathieu equation having PT symmetry

abstract: I have applied multiple-scale perturbation theory to a generalized complex $PT$-symmetric Mathieu equation in order to find the stability boundaries between bounded and unbounded solutions. The analysis suggests that the non-Hermitian parameter present in the equation can be used to control the shape and curvature of these boundaries. Although this was suggested earlier by several authors, analytic formulas for the boundary curves were not given. This paper is a first attempt to fill this gap in the theory.


title: Predictive Learning on Hidden Tree-Structured Ising Models

abstract: We provide high-probability sample complexity guarantees for exact structure recovery and accurate predictive learning using noise-corrupted samples from an acyclic (tree-shaped) graphical model. The hidden variables follow a tree-structured Ising model distribution, whereas the observable variables are generated by a binary symmetric channel taking the hidden variables as its input (flipping each bit independently with some constant probability $q\in [0,1/2)$). In the absence of noise, predictive learning on Ising models was recently studied by Bresler and Karzand (2020); this paper quantifies how noise in the hidden model impacts the tasks of structure recovery and marginal distribution estimation by proving upper and lower bounds on the sample complexity. Our results generalize state-of-the-art bounds reported in prior work, and they exactly recover the noiseless case ($q=0$). In fact, for any tree with $p$ vertices and probability of incorrect recovery $\delta>0$, the sufficient number of samples remains logarithmic as in the noiseless case, i.e., $\mathcal{O}(\log(p/\delta))$, while the dependence on $q$ is $\mathcal{O}\big( 1/(1-2q)^{4} \big)$, for both aforementioned tasks. We also present a new equivalent of Isserlis' Theorem for sign-valued tree-structured distributions, yielding a new low-complexity algorithm for higher-order moment estimation.

categories: stat.ML cs.IT cs.LG math.IT math.ST stat.TH


title: Numerical Methods for Fast Nonlinear Fourier Transformation, Part I: Exponential Runge-Kutta and Linear Multistep Methods

abstract: The main objective of this series of papers is to explore the entire landscape of numerical methods for fast nonlinear Fourier transformation (NFT) within the class of integrators known as the exponential integrators. In this paper, we explore the theoretical aspects of exponential Runge-Kutta (RK) and linear multistep (LM) methods, in particular, the stability and convergence of these methods via the transfer matrix formulation. The analysis carried out in the paper shows that while the exponential LM methods are naturally amenable to FFT-based fast polynomial arithmetic, the RK methods require equispaced nodes to achieve that. Therefore, each these family of methods is capable of yielding a family of fast NFT algorithms such that the scattering coefficients can be computed with a complexity of $\mathscr{O}(N\log^2N)$ and a rate of convergence given by $\mathscr{O}(N^{-p})$ where $N$ is the number of samples of the signal and $p$ is order of the underlying discretization scheme. Further, while RK methods can accommodate vanishing as well as periodic boundary conditions, the LM methods can only handle the former type of boundary conditions without requiring a starting procedure. The ideas presented in this paper extend naturally to the family of integrators known as general linear methods which will be explored in a forthcoming paper.


title: Size effect of soft phonon dispersion law in nanosized ferroics

abstract: Using Landau-Ginsburg-Devonshire theory, we derive and analyze analytical expressions for the frequency dispersion of soft phonon modes in nanosized ferroics and perform numerical calculations for a thin SrTiO3 film. We revealed the pronounced "true" size effect in the dependence of soft phonon spatial dispersion on the film thickness and predict that it can lead to the "apparent" or "false" size effect of dynamic flexoelectric coupling constants


title: Thermal resistance by transition between collective and non-collective phonon flows in graphitic materials

abstract: Phonons in graphitic materials exhibit strong normal scattering (N-scattering) compared to umklapp scattering (U-scattering). The strong N-scattering cause collective phonon flow, unlike the relatively common cases where U-scattering is dominant. If graphitic materials have finite size and contact with hot and cold reservoirs emitting phonons with non-collective distribution, N-scattering change the non-collective phonon flow to the collective phonon flow near the interface between graphitic material and a heat reservoir. We study the thermal resistance by N-scattering during the transition between non-collective and collective phonon flows. Our Monte Carlo solution of Peierls-Boltzmann transport equation shows that the N-scattering in graphitic materials reduce heat flux from the ballistic case by around 15%, 30%, and 40% at 100, 200, and 300 K, respectively. This is significantly larger than ~ 5% reduction of Debye crystal with similar Debye temperature (~ 2300 K). We associate the large reduction of heat flux by N-scattering with the non-linear dispersion and multiple phonon branches with different group velocities of graphitic materials.


title: Schur Multipliers and Second Quandle Homology

abstract: We define a map from second quandle homology to the Schur multiplier and examine its properties. Furthermore, we express the second homology of Alexander quandles in terms of exterior algebras. Additionally, we present a self-contained proof of its structure and provide some computational examples.


title: A classification of monotone ribbons with full Schur support with application to the classification of full equivalence classes

abstract: We consider ribbon shapes, not necessarily connected, whose rows, with at least two boxes in each, are in monotone length order. These ribbons are uniquely defined by a pair of partitions: the row partition consisting of the row lengths in decreasing order, and the overlapping partition whose entries count the total number of columns with two boxes in the successive ribbon shapes obtained by sequentially subtracting the longest row. The support of such ribbon Schur functions, considered as a subposet of the dominance order lattice on partitions, has the row partition as bottom element, and, as top element, the partition whose two parts consist of the total number of columns, and the total number of columns of length two respectively. We give a complete system of linear inequalities in terms of the partition pair defining the aforesaid ribbon shape under which the ribbon Schur function attains all the Schur interval when expanded in the basis of Schur functions. We then conclude that the Gaetz-Hardt-Sridhar necessary condition for a connected ribbon to have full equivalence class is equivalent to the condition for a monotone connected ribbon to have full Schur support. That is, the set of partitions with full equivalence class is a subset of those monotone connected ribbons with full Schur support. M. Gaetz, W. Hardt and S. Sridhar conjectured that the necessary condition is also sufficient which translates now to every monotone connected ribbon with full Schur support has full equivalence class. The main tool of our analysis is the structure of the companion tableau of a ribbon Littlewood-Richardson (LR) tableau detected by the descent set defined by the composition whose parts are the ribbon row lengths.


title: Rotation Invariant Descriptors for Galaxy Morphological Classification

abstract: The detection of objects that are multi-oriented is a difficult pattern recognition problem. In this paper, we propose to evaluate the performance of different families of descriptors for the classification of galaxy morphologies. We investigate the performance of the Hu moments, Flusser moments, Zernike moments, Fourier-Mellin moments, and ring projection techniques based on 1D moment and the Fourier transform. We consider two main datasets for the performance evaluation. The first dataset is an artificial dataset based on representative templates from 11 types of galaxies, which are evaluated with different transformations (noise, smoothing), alone or combined. The evaluation is based on image retrieval performance to estimate the robustness of the rotation invariant descriptors with this type of images. The second dataset is composed of real images extracted from the Galaxy Zoo 2 project. The binary classification of elliptical and spiral galaxies is achieved with pre-processing steps including morphological filtering and a Laplacian pyramid. For the binary classification, we compare the different set of features with Support Vector Machines (SVM), Extreme Learning Machine, and different types of linear discriminant analysis techniques. The results support the conclusion that the proposed framework for the binary classification of elliptical and spiral galaxies provides an area under the ROC curve reaching 99.54%, proving the robustness of the approach for helping astronomers to study galaxies.


title: Singular symplectic cotangent bundle reduction of gauge field theory

abstract: We prove a theorem on singular symplectic cotangent bundle reduction in the Fr\'echet setting and apply it to Yang-Mills-Higgs theory with special emphasis on the Higgs sector of the Glashow-Weinberg-Salam model. For the latter model we give a detailed description of the reduced phase space and show that the singular structure is encoded in a finite-dimensional Lie group action.


title: Non-local Operational Anisotropic Diffusion Filter

abstract: High-frequency noise is present in several modalities of medical images. It originates from the acquisition process and may be related to the scanner configurations, the scanned body, or to other external factors. This way, prospective filters are an important tool to improve the image quality. In this paper, we propose a non-local weighted operational anisotropic diffusion filter and evaluate its effect on magnetic resonance images and on kV/CBCT radiotherapy images. We also provide a detailed analysis of non-local parameter settings. Results show that the new filter enhances previous local implementations and has potential application in radiotherapy treatments.


title: Inter-Core Crosstalk in Multicore Fibers: Impact on 56-Gbaud/{\lambda}/Core PAM-4 Transmission

abstract: We experimentally demonstrate the impact of inter-core crosstalk in multicore fibers on 56Gbaud PAM-4 signal quality after 2.5-km transmission over a weakly-coupled and uncoupled sevencore fibers, revealing the crosstalk dependence on carrier central wavelength in range of 1540-1560 nm.


title: Capillary transport in low saturated sands: superfast non-linear diffusion model versus direct experimental observations

abstract: We have established in a pilot study, that the spreading of liquids in sandy porous materials at low levels of saturation, typically less than ten percent of the available void space, has very distinctive features in comparison to that at higher saturation levels. In particular, it has been shown, on theoretical grounds, that the spreading is controlled by a special type of diffusional process, that its physics can be captured by an equation of the super-fast diffusion class, and that these findings were supported by first-of-a-kind experiments. In this paper, we take these findings to the next level including deeper examination and exposition of the theory, an expanded set of experiments to address scaling properties, and systematic evaluations of the predictive performance against these experimental data, keeping in mind also potential practical applications.


title: Joint Computation Offloading and Resource Allocation in Cloud Based Wireless HetNets

abstract: In this paper, we study the joint computation offloading and resource allocation problem in the two-tier wireless heterogeneous network (HetNet). Our design aims to optimize the computation offloading to the cloud jointly with the subchannel allocation to minimize the maximum (min-max) weighted energy consumption subject to practical constraints on bandwidth, computing resource and allowable latency for the multi-user multitask computation system. To tackle this non-convex mixed integer non-linear problem (MINLP), we employ the bisection search method to solve it where we propose a novel approach to transform and verify the feasibility of the underlying problem in each iteration. In addition, we propose a low-complexity algorithm, which can decrease the number of binary optimization variables and enable more scalable computation offloading optimization in the practical wireless HetNets. Numerical studies confirm that the proposed design achieves the energy saving gains about 55% in comparison with the local computation scheme under the strict required latency of 0.1s.


title: Patient-Centric Cellular Networks Optimization using Big Data Analytics

abstract: Big data analytics is one of the state-of-the-art tools to optimize networks and transform them from merely being a blind tube that conveys data, into a cognitive, conscious, and self-optimizing entity that can intelligently adapt according to the needs of its users. This, in fact, can be regarded as one of the highest forthcoming priorities of future networks. In this paper, we propose a system for Out-Patient (OP) centric Long Term Evolution-Advanced (LTE-A) network optimization. Big data harvested from the OPs' medical records, along with current readings from their body sensors are processed and analyzed to predict the likelihood of a life-threatening medical condition, for instance, an imminent stroke. This prediction is used to ensure that the OP is assigned an optimal LTE-A Physical Resource Blocks (PRBs) to transmit their critical data to their healthcare provider with minimal delay. To the best of our knowledge, this is the first time big data analytics are utilized to optimize a cellular network in an OP-conscious manner. The PRBs assignment is optimized using Mixed Integer Linear Programming (MILP) and a real-time heuristic. Two approaches are proposed, the Weighted Sum Rate Maximization (WSRMax) approach and the Proportional Fairness (PF) approach. The approaches increased the OPs' average SINR by 26.6% and 40.5%, respectively. The WSRMax approach increased the system's total SINR to a level higher than that of the PF approach, however, the PF approach reported higher SINRs for the OPs, better fairness and a lower margin of error.


title: Charged Higgs Bosons in Naturally Aligned Two Higgs Doublet Models at the LHC

abstract: Measurements of a Higgs boson at the Large Hadron Collider (LHC) have become increasingly consistent with the predictions of the Standard Model (SM). This fact puts severe constraints on many potential low-energy extensions of the Higgs sector of the SM. In the well-known Two Higgs Doublet Model (2HDM), an `alignment limit' of parameters readily furnishes one SM-like scalar, and can be achieved naturally through an underlying symmetry. Among the other physical states of the 2HDM, a charged scalar $H^\pm$ would provide striking evidence of new physics if observed. We propose a novel technique for the observation of the process $pp\to t b H^\pm \to t\bar{t}b\bar{b}$ in the dileptonic decay channel at the LHC. The reconstruction of events in this channel is complicated by multiple $b$-jets and unobserved neutrinos in the final state. To determine the neutrino momenta, we implement a neutrino weighting procedure to study, for the first time, the $t\bar{t}b\bar{b}$ signature. We further train a pair of boosted decision trees to reconstruct and classify signal events. We determine the resulting reach within the context of naturally aligned 2HDMs, such as the Maximally Symmetric Two Higgs Doublet Model (MS-2HDM). By testing at the integrated luminosity of 150 fb$^{-1}$ achieved in Run 2 of the LHC, we find that this channel may restrict the parameter space of a Type-II MS-2HDM with charged Higgs masses as high as 680 GeV.


title: Telecom Compatibility Validation of Quantum Key Distribution Co-existing with 112 Gbps/{\lambda}/core Data Transmission in Non-Trench and Trench-Assistant Multicore Fibers

abstract: We experimentally characterize photon leakage from 112Gbps data channels in both non-trench and trench-assistant 7-core fibers, demonstrating telecom compatibility for QKD co-existing with high-speed data transmission when a proper core/wavelength allocation is carried out.


title: Deep Neural Networks Meet CSI-Based Authentication

abstract: The first step of a secure communication is authenticating legible users and detecting the malicious ones. In the last recent years, some promising schemes proposed using wireless medium network's features, in particular, channel state information (CSI) as a means for authentication. These schemes mainly compare user's previous CSI with the new received CSI to determine if the user is in fact what it is claiming to be. Despite high accuracy, these approaches lack the stability in authentication when the users rotate in their positions. This is due to a significant change in CSI when a user rotates which mislead the authenticator when it compares the new CSI with the previous ones. Our approach presents a way of extracting features from raw CSI measurements which are stable towards rotation. We extract these features by the means of a deep neural network. We also present a scenario in which users can be {efficiently} authenticated while they are at certain locations in an environment (even if they rotate); and, they will be rejected if they change their location. Also, experimental results are presented to show the performance of the proposed scheme.


title: Comments on One-Form Global Symmetries and Their Gauging in 3d and 4d

abstract: We study 3d and 4d systems with a one-form global symmetry, explore their consequences, and analyze their gauging. For simplicity, we focus on $\mathbb{Z}_N$ one-form symmetries. A 3d topological quantum field theory (TQFT) $\mathcal{T}$ with such a symmetry has $N$ special lines that generate it. The braiding of these lines and their spins are characterized by a single integer $p$ modulo $2N$. Surprisingly, if $\gcd(N,p)=1$ the TQFT factorizes $\mathcal{T}=\mathcal{T}'\otimes \mathcal{A}^{N,p}$. Here $\mathcal{T}'$ is a decoupled TQFT, whose lines are neutral under the global symmetry and $\mathcal{A}^{N,p}$ is a minimal TQFT with the $\mathbb{Z}_N$ one-form symmetry of label $p$. The parameter $p$ labels the obstruction to gauging the $\mathbb{Z}_N$ one-form symmetry; i.e.\ it characterizes the 't Hooft anomaly of the global symmetry. When $p=0$ mod $2N$, the symmetry can be gauged. Otherwise, it cannot be gauged unless we couple the system to a 4d bulk with gauge fields extended to the bulk. This understanding allows us to consider $SU(N)$ and $PSU(N)$ 4d gauge theories. Their dynamics is gapped and it is associated with confinement and oblique confinement -- probe quarks are confined. In the $PSU(N)$ theory the low-energy theory can include a discrete gauge theory. We will study the behavior of the theory with a space-dependent $\theta$-parameter, which leads to interfaces. Typically, the theory on the interface is not confining. Furthermore, the liberated probe quarks are anyons on the interface. The $PSU(N)$ theory is obtained by gauging the $\mathbb{Z}_N$ one-form symmetry of the $SU(N)$ theory. Our understanding of the symmetries in 3d TQFTs allows us to describe the interface in the $PSU(N)$ theory.


title: Pible: Battery-Free Mote for Perpetual Indoor BLE Applications

abstract: Smart building applications require a large-scale deployment of sensors distributed across the environment. Recent innovations in smart environments are driven by wireless networked sensors as they are easy to deploy. However, replacing these batteries at scale is a non-trivial, labor-intensive task. Energy harvesting has emerged as a potential solution to avoid battery replacement but requires compromises such as application specific design, simplified communication protocol or reduced quality of service. We explore the design space of battery-free sensor nodes using commercial off the shelf components, and present Pible: a Perpetual Indoor BLE sensor node that leverages ambient light and can support numerous smart building applications. We analyze node-lifetime, quality of service and light availability trade-offs and present a predictive algorithm that adapts to changing lighting conditions to maximize node lifetime and application quality of service. Using a 20 node, 15-day deployment in a real building under varying lighting conditions, we show feasible applications that can be implemented using Pible and the boundary conditions under which they can fail.


title: Text Data Augmentation Made Simple By Leveraging NLP Cloud APIs

abstract: In practice, it is common to find oneself with far too little text data to train a deep neural network. This "Big Data Wall" represents a challenge for minority language communities on the Internet, organizations, laboratories and companies that compete the GAFAM (Google, Amazon, Facebook, Apple, Microsoft). While most of the research effort in text data augmentation aims on the long-term goal of finding end-to-end learning solutions, which is equivalent to "using neural networks to feed neural networks", this engineering work focuses on the use of practical, robust, scalable and easy-to-implement data augmentation pre-processing techniques similar to those that are successful in computer vision. Several text augmentation techniques have been experimented. Some existing ones have been tested for comparison purposes such as noise injection or the use of regular expressions. Others are modified or improved techniques like lexical replacement. Finally more innovative ones, such as the generation of paraphrases using back-translation or by the transformation of syntactic trees, are based on robust, scalable, and easy-to-use NLP Cloud APIs. All the text augmentation techniques studied, with an amplification factor of only 5, increased the accuracy of the results in a range of 4.3% to 21.6%, with significant statistical fluctuations, on a standardized task of text polarity prediction. Some standard deep neural network architectures were tested: the multilayer perceptron (MLP), the long short-term memory recurrent network (LSTM) and the bidirectional LSTM (biLSTM). Classical XGBoost algorithm has been tested with up to 2.5% improvements.


title: Parallel contact-aware simulations of deformable particles in 3D Stokes flow

abstract: We present a parallel-scalable method for simulating non-dilute suspensions of deformable particles immersed in Stokesian fluid in three dimensions. A critical component in these simulations is robust and accurate collision handling. This work complements our previous work [L. Lu, A. Rahimian, and D. Zorin. Contact-aware simulations of particulate Stokesian suspensions. Journal of Computational Physics 347C: 160-182] by extending it to 3D and by introducing new parallel algorithms for collision detection and handling. We use a well-established boundary integral formulation with spectral Galerkin method to solve the fluid flow. The key idea is to ensure an interference-free particle configuration by introducing explicit contact constraints into the system. While such constraints are typically unnecessary in the formulation they make it possible to eliminate catastrophic loss of accuracy in the discretized problem by preventing contact explicitly. The incorporation of contact constraints results in a significant increase in stable time-step size for locally-implicit time-stepping and a reduction in the necessary number of discretization points for stability. Our method maintains the accuracy of previous methods at a significantly lower cost for dense suspensions and the time step size is independent from the volume fraction. Our method permits simulations with high volume fractions; we report results with up to 60% volume fraction. We demonstrated the parallel scaling of the algorithms on up to 16K CPU cores.


title: Stability of the centers of the symplectic groups rings $\matbb{Z}[Sp_{2n}(q)]$

abstract: We investigate the structure constants of the center $\mathcal{H}_n$ of the group algebra $Sp_{n}(q)$ over a finite field. The reflection length on the group $GL_{2n}(q)$ induces a filtration on the algebras $\mathcal{H}_n$. We prove that the structure constants of the associated filtered algebra $\mathcal{S}_n$ are independent of $n$. As a technical tool in the proof, we determine the growth of the centralizers under the embedding $Sp_m(q)\subset Sp_{m+l}(q)$ and we show that the index of the centralizer of $g\in Sp_m(q)$ in the centralizer of $g\in Sp_{m+k}$ is equal to $q^{2ld}|Sp_{r+l}(q)||Sp_{r}(q)|^{-1}$ for some $d$ and $r$ which are uniquely determined by the conjugacy class of $g$ in $GL_{2n}(q).$


title: On Universality of Classical Probability with Contextually Labeled Random Variables: Response to A. Khrennikov

abstract: In his constructive and well-informed commentary, Andrei Khrennikov acknowledges a privileged status of classical probability theory with respect to statistical analysis. He also sees advantages offered by the Contextuality-by-Default theory, notably, that it `demystifies quantum mechanics by highlighting the role of contextuality,' and that it can detect and measure contextuality in inconsistently connected systems. He argues, however, that classical probability theory may have difficulties in describing empirical phenomena if they are described entirely in terms of observable events. We disagree: contexts in which random variables are recorded are as observable as the variables' values. Khrennikov also argues that the Contextuality-by-Default theory suffers the problem of non-uniqueness of couplings. We disagree that this is a problem: couplings are all possible ways of imposing counterfactual joint distributions on random variables that de facto are not jointly distributed. The uniqueness of modeling experiments by means of quantum formalisms brought up by Khrennikov is achieved for the price of additional, substantive assumptions. This is consistent with our view of quantum theory as a special-purpose generator of classical probabilities. Khrennikov raises the issue of `mental signaling,' by which he means inconsistent connectedness in behavioral systems. Our position is that it is as inherent to behavioral systems as their stochasticity.


title: Context is Key: New Approaches to Neural Coherence Modeling

abstract: We formulate coherence modeling as a regression task and propose two novel methods to combine techniques from our setup with pairwise approaches. The first of our methods is a model that we call "first-next," which operates similarly to selection sorting but conditions decision-making on information about already-sorted sentences. The second consists of a technique for adding context to regression-based models by concatenating sentence-level representations with an encoding of its corresponding out-of-order paragraph. This latter model achieves Kendall-tau distance and positional accuracy scores that match or exceed the current state-of-the-art on these metrics. Our results suggest that many of the gains that come from more complex, machine-translation inspired approaches can be achieved with simpler, more efficient models.


title: Impact of Intervals on the Emotional Effect in Western Music

abstract: Every art form ultimately aims to invoke an emotional response over the audience, and music is no different. While the precise perception of music is a highly subjective topic, there is an agreement in the "feeling" of a piece of music in broad terms. Based on this observation, in this study, we aimed to determine the emotional feeling associated with short passages of music; specifically by analyzing the melodic aspects. We have used the dataset put together by Eerola et. al. which is comprised of labeled short passages of film music. Our initial survey of the dataset indicated that other than "happy" and "sad" labels do not possess a melodic structure. We transcribed the main melody of the happy and sad tracks and used the intervals between the notes to classify them. Our experiments have shown that treating a melody as a bag-of-intervals do not possess any predictive power whatsoever, whereas counting intervals with respect to the key of the melody yielded a classifier with 85% accuracy.


title: Noise investigation of a dual-frequency VECSEL for application to Cesium clocks

abstract: We theoretically and experimentally study the noise of a class-A dual-frequency vertical external cavity surface emitting laser operating at Cesium clock wavelength. The intensity noises of the two orthogonally polarized modes and the phase noise of their beatnote are investigated. The intensity noises of the two modes and their correlations are well predicted by a theory based on coupled rate equations. The phase noise of the beatnote is well described by considering both thermal effects and the effect of phase-amplitude coupling. The good agreement between theory and experiment indicates possible ways to further decrease the laser noises.


title: Refraction efficiency of Huygens' and bianisotropic terahertz metasurfaces

abstract: Metasurfaces are an enabling technology for complex wave manipulation functions, including in the terahertz frequency range, where they are expected to advance security, imaging, sensing, and communications technology. For operation in transmission, Huygens' metasurfaces are commonly used, since their good impedance match to the surrounding media minimizes reflections and maximizes transmission. Recent theoretical work has shown that Huygens' metasurfaces are non-optimal, particularly for large angles of refraction, and that to eliminate reflections and spurious diffracted beams it is necessary to use a bianisotropic metasurface. However, it remains to be demonstrated how significant the efficiency improvement is when using bianisotropic metasurfaces, considering all the non-ideal features that arise when implementing the metasurface design with real meta-atoms. Here we compare concrete terahertz metasurface designs based on the Huygens' and Omega-type bianisotropic approaches, demonstrating anomalous refraction angles for 55 degrees, and 70 degrees. We show that for the lower angle of 55 degrees, there is no significant improvement when using the bianisotropic design, whereas for refraction at 70 degrees the bianisotropic design shows much higher efficiency and fidelity of refraction into the designed direction. We also demonstrate the strong perturbations caused by near-field interaction, both between and within cells, which we compensate using numerical optimization.


title: Polarisabilities from Compton Scattering on 3He

abstract: This executive summary of recent theory progress in Compton scattering off 3He focuses on determining neutron polarisabilities; see ref. [2] and references therein for details and a better bibliography. Prepared for the Proceedings of the 22nd International Conference on Few-Body Problems in Physics, Caen 9-13 July 2018.


title: First-order synchronization transition in a large population of strongly coupled relaxation oscillators

abstract: Onset and loss of synchronization in coupled oscillators are of fundamental importance in understanding emergent behavior in natural and man-made systems, which range from neural networks to power grids. We report on experiments with hundreds of strongly coupled photochemical relaxation oscillators that exhibit a discontinuous synchronization transition with hysteresis, as opposed to the paradigmatic continuous transition expected from the widely used weak coupling theory. The resulting first-order transition is robust with respect to changes in network connectivity and natural frequency distribution. This allows us to identify the relaxation character of the oscillators as the essential parameter that determines the nature of the synchronization transition. We further support this hypothesis by revealing the mechanism of the transition, which cannot be accounted for by standard phase reduction techniques.


title: Guided Exploration of Human Intentions for Human-Robot Interaction

abstract: Robot understanding of human intentions is essential for fluid human-robot interaction. Intentions, however, cannot be directly observed and must be inferred from behaviors. We learn a model of adaptive human behavior conditioned on the intention as a latent variable. We then embed the human behavior model into a principled probabilistic decision model, which enables the robot to (i) explore actively in order to infer human intentions and (ii) choose actions that maximize its performance. Furthermore, the robot learns from the demonstrated actions of human experts to further improve exploration. Preliminary experiments in simulation indicate that our approach, when applied to autonomous driving, improves the efficiency and safety of driving in common interactive driving scenarios.


title: Hyperuniformity on spherical surfaces

abstract: In this work we present a study on the characterization of ordered and disordered hyperuniform point distributions on spherical surfaces. In spite of the extensive literature on disordered hyperuniform systems in Euclidean geometries, to date few works have dealt with the problem of hyperuniformity in curved spaces. As a matter of fact, some systems that display disordered hyperuniformity, like the space distribution of photoreceptors in avian retina, actually occur on curved surfaces. Here we will focus on the local particle number variance and its dependence on the size of the sampling window (which we take to be a spherical cap) for regular and uniform point distributions, as well as for equilibrium configurations of fluid particles interacting through Lennard-Jones, dipole-dipole and charge-charge potentials. We will show how the scaling of the local number variance enables the characterization of hyperuniform point patterns also on spherical surfaces.


title: Dynamical system induced by quantum walk

abstract: We consider the Grover walk model on a connected finite graph with two infinite length tails and we set an $\ell^\infty$-infinite external source from one of the tails as the initial state. We show that for any connected internal graph, a stationary state exists, moreover a perfect transmission to the opposite tail always occurs in the long time limit. We also show that the lower bound of the norm of the stationary measure restricted to the internal graph is proportion to the number of edges of this graph. Furthermore when we add more tails (e.g., $r$-tails) to the internal graph, then we find that from the temporal and spatial global view point, the scattering to each tail in the long time limit coincides with the local one-step scattering manner of the Grover walk at a vertex whose degree is $(r+1)$.


title: Compatibility Between Wind Turbines and the Radio Astronomy Service

abstract: Modern radio astronomical facilities are able to detect extremely weak electromagnetic signals not only from the universe but also from man-made radio frequency interference of various origins. These range from wanted signals to unwanted out-of-band emission of radio services and applications to electromagnetic interference produced by all kinds of electronic and electric devices. Energy harvesting wind turbines are not only equipped with electric power conversion hardware but also copious amounts of electronics to control and monitor the turbines. A wind turbine in the vicinity of a radio telescope could therefore lead to harmful interference, corrupting the measured astronomical data. Many observatories seek to coordinate placement of new wind farms with wind turbine manufacturers and operators, as well as with the local planning authorities, to avoid such a situation. In our study, we provide examples as well as guidelines for the determination of the separation distances between wind turbines and radio observatories, to enable a benign co-existence for both. The proposed calculations entail three basic steps. At first, the anticipated maximum emitted power level based on the European EN 550011 (CISPR-11) standard, which applies to industrial devices, is determined. Then secondly, the propagation loss along the path to the radio receiver is computed via a model provided by the international telecommunication union. Finally, the received power is compared to the permitted power limit that pertains in the protected radio astronomical observing band under consideration. This procedure may be carried out for each location around a telescope site, in order to obtain a map of potentially problematic wind turbine positions.


title: First determination of the $\rho $ parameter at $\sqrt{s} = 13$ TeV -- probing the existence of a colourless three-gluon bound state

abstract: The TOTEM experiment at the LHC has performed the first measurement at $\sqrt{s} = 13$ TeV of the $\rho$ parameter, the real to imaginary ratio of the nuclear elastic scattering amplitude at $t=0$, obtaining the following results: $\rho = 0.09 \pm 0.01$ and $\rho = 0.10 \pm 0.01$, depending on different physics assumptions and mathematical modelling. The unprecedented precision of the $\rho$ measurement, combined with the TOTEM total cross-section measurements in an energy range larger than 10 TeV (from 2.76 to 13 TeV), has implied the exclusion of all the models classified and published by COMPETE. The $\rho$ results obtained by TOTEM are compatible with the predictions, from alternative theoretical models both in the Regge-like framework and in the QCD framework, of a colourless 3-gluon bound state exchange in the $t$-channel of the proton-proton elastic scattering. On the contrary, if shown that the 3-gluon bound state $t$-channel exchange is not of importance for the description of elastic scattering, the $\rho$ value determined by TOTEM would represent a first evidence of a slowing down of the total cross-section growth at higher energies. The very low-$|t|$ reach allowed also to determine the absolute normalisation using the Coulomb amplitude for the first time at the LHC and obtain a new total proton-proton cross-section measurement $\sigma_{tot} = 110.3 \pm 3.5$ mb, completely independent from the previous TOTEM determination. Combining the two TOTEM results yields $\sigma_{tot} = 110.5 \pm 2.4$ mb.


title: Reducible and irreducible approximation of complex symmetric operators

abstract: This paper aims to study reducible and irreducible approximation in the set $\textsl{CSO}$ of all complex symmetric operators on a separable, complex Hilbert space $\mathcal H$. When ${\rm dim} \mathcal H=\infty$, it is proved that both those reducible ones and those irreducible ones are norm dense in $\textsl{CSO}$. When ${\rm dim} \mathcal H<\infty$, irreducible complex symmetric operators constitute an open, dense subset of $\textsl{CSO}$.


title: Impact of the distribution of recovery rates on disease spreading in complex networks

abstract: We study a general epidemic model with arbitrary recovery rate distributions. This simple deviation from the standard setup is sufficient to prove that heterogeneity in the dynamical parameters can be as important as the more studied structural heterogeneity. Our analytical solution is able to predict the shift in the critical properties induced by heterogeneous recovery rates. Additionally, we show that the critical value of infectivity tends to be smaller than the one predicted by quenched mean-field approaches in the homogeneous case and that it can be linked to the variance of the recovery rates. We then illustrate the role of dynamical--structural correlations, which allow for a complete change in the critical behavior. We show that it is possible for a power-law network topology to behave similarly to a homogeneous structure by an appropriate tuning of its recovery rates, and vice versa. Finally, we show how heterogeneity in recovery rates affects the network localization properties of the spreading process.


title: Quantum Algorithms, Architecture, and Error Correction

abstract: Quantum algorithms have the potential to provide exponential speedups over some of the best known classical algorithms. These speedups may enable quantum devices to solve currently intractable problems such as those in the fields of optimization, material science, chemistry, and biology. Thus, the realization of large-scale, reliable quantum-computers will likely have a significant impact on the world. For this reason, the focus of this dissertation is on the development of quantum-computing applications and robust, scalable quantum-architectures. I begin by presenting an overview of the language of quantum computation. I then, in joint work with Ojas Parekh, analyze the performance of the quantum approximate optimization algorithm (QAOA) on a graph problem called Max Cut. Next, I present a new stabilizer simulation algorithm that gives improved runtime performance for topological stabilizer codes. After that, in joint work with Andrew Landahl, I present a new set of procedures for performing logical operations called "color-code lattice-surgery." Finally, I describe a software package I developed for studying, developing, and evaluating quantum error-correcting codes under realistic noise.


title: A multi-dimensional, moment-accelerated deterministic particle method for time-dependent, multi-frequency thermal radiative transfer problems

abstract: Thermal Radiative Transfer (TRT) is the dominant energy transfer mechanism in high-energy density physics with applications in inertial confinement fusion and astrophysics. The stiff interactions between the material and radiation fields make TRT problems challenging to model. In this study, we propose a multi-dimensional extension of the deterministic particle (DP) method. The DP method combines aspects from both particle and deterministic methods. If the emission source is known \apriori, and no physical scattering is present, the intensity of a particle can be integrated analytically. This introduces no statistical noise compared to Monte-Carlo methods, while maintaining the flexibility of particle methods. The method is closely related to the popular method of long characteristics. The combination of the DP-method with a discretely-consistent, nonlinear, gray low-order system enables an efficient solution algorithm for multi-frequency TRT problems. We demonstrate with numerical examples that the use of a linear-source approximation based on spatial moments improves the behavior of our method in the thick diffusion limit significantly.


title: Spectral Analysis for Non-Hermitian Matrices and Directed Graphs

abstract: We generalize classical results in spectral graph theory and linear algebra more broadly, from the case where the underlying matrix is Hermitian to the case where it is non-Hermitian. New admissibility conditions are introduced to replace the Hermiticity condition. We prove new variational estimates of the Rayleigh quotient for non-Hermitian matrices. As an application, a new Delsarte-Hoffman-type bound on the size of the largest independent set in a directed graph is developed. Our techniques consist in quantifying the impact of breaking the Hermitian symmetry of a matrix and are broadly applicable.


title: Guessing probability under unlimited known-plaintext attack on secret keys for Y00 quantum stream cipher by quantum multiple hypotheses testing

abstract: Although quantum key distribution is regarded as promising secure communication, security of Y00 protocol proposed by Yuen in 2000 for the affinity to conventional optical communication is not well-understood yet; its security has been evaluated only by the eavesdropper's error probabilities of detecting individual signals or masking size, the number of hidden signal levels under quantum and classical noise. Our study is the first challenge of evaluating the guessing probabilities on shared secret keys for pseudorandom number generators in a simplified Y00 communication system based on quantum multiple hypotheses testing theory. The result is that even unlimitedly long known-plaintext attack only lets the eavesdropper guess the shared secret keys of limited lengths with a probability strictly < 1. This study will give some insights for detailed future works on this quantum communication protocol.


title: Large angular scale multipoles at redshift ~0.8

abstract: We prepare the full sky radio galaxy map ($|b|>10^{\circ}$) using the north NVSS and south SUMSS galaxy catalogs and study the large scale multipoles anomalies. These galaxies are roughly at redshift $z \sim 0.8$ and therefore tracing the matter distribution at very large scales. The quadruple and octopole from radio galaxy catalog are consistent with $\Lambda$CDM for a reasonable value of galaxy bias and we do not find dipole--quadruple--octopole alignment as seen in CMB temperature maps. The quadrupole direction is roughly $46^{\circ}$ away from dipole, and octopole direction is approximately $33^{\circ}$ from dipole. The angle between quadrupole and octopole is around $70^\circ$ degree. We have large errors in multipole directions due to shot noise, even so with this data we are able to rule out dipole--quadruple and quadruple--octopole alignment. The magnitude of all multipoles, except dipole, are roughly consistent with $\Lambda$CDM for reasonable galaxy bias. The dipole magnitude remains inconsistent with CMB as reported in previous studies. The results may impose stringent constraints on cosmological models with large scale anisotropy features.


title: Spectral Theory in a Twisted Groupoid Setting: Spectral Decompositions, Localization and Fredholmness

abstract: We study bounded operators defined in terms of the regular representations of the $C^*$-algebra of an amenable, Hausdorff, second countable locally compact groupoid endowed with a continuous $2$-cocycle. We concentrate on spectral quantities associated to natural quotients of this twisted algebra, such as the essential spectrum, the essential numerical range, and Fredholm properties. We obtain decompositions for the regular representations associated to units of the groupoid belonging to a free locally closed orbit, in terms of spectral quantities attached to points (or orbits) in the boundary of this main orbit. As examples, we discuss various classes of magnetic pseudo-differential operators on nilpotent groups. We also prove localization and non-propagation properties associated to suitable parts of the essential spectrum. These are applied to twisted groupoids having a totally intransitive groupoid restriction at the boundary.


title: The Jiminy Advisor: Moral Agreements Among Stakeholders Based on Norms and Argumentation

abstract: An autonomous system is constructed by a manufacturer, operates in a society subject to norms and laws, and interacts with end users. All of these actors are stakeholders affected by the behavior of the autonomous system. We address the challenge of how the ethical views of such stakeholders can be integrated in the behavior of an autonomous system. We propose an ethical recommendation component called Jiminy which uses techniques from normative systems and formal argumentation to reach moral agreements among stakeholders. A Jiminy represents the ethical views of each stakeholder by using normative systems, and has three ways of resolving moral dilemmas that involve the opinions of the stakeholders. First, the Jiminy considers how the arguments of the stakeholders relate to one another, which may already resolve the dilemma. Secondly, the Jiminy combines the normative systems of the stakeholders such that the combined expertise of the stakeholders may resolve the dilemma. Thirdly, and only if these two other methods have failed, the Jiminy uses context-sensitive rules to decide which of the stakeholders take preference over the others. At the abstract level, these three methods are characterized by adding arguments, adding attacks between arguments, and revising attacks between arguments. We show how a Jiminy can be used not only for ethical reasoning and collaborative decision-making, but also to provide explanations about ethical behavior.


title: Radial basis function-generated finite differences with Bessel weights for the 2D Helmholtz equation

abstract: In this paper we obtain approximated numerical solutions for the 2D Helmholtz equation using a radial basis function-generated finite difference scheme (RBF-FD), where weights are calculated by applying an oscillatory radial basis function given in terms of Bessel functions of the first kind. The problem of obtaining weights by local interpolation is ill-conditioned; we overcome this difficulty by means of regularization of the interpolation matrix by perturbing its diagonal. The condition number of this perturbed matrix is controlled according to a prescribed value of a regularization parameter. Different numerical tests are performed in order to study convergence and algorithmic complexity. As a result, we verify that dispersion and pollution effects are mitigated.


title: Vanishing OPE Coefficients in 4d N=2 SCFTs

abstract: We compute the superconformal characters of various short multiplets in 4d N=2 superconformal algebra, from which selection rules for operator products are obtained. Combining with the superconformal index, we show that a particular short multiplet appearing in the n-fold product of stress-tensor multiplet is absent in the $(A_1, A_{2n})$ Argyres-Douglas (AD) theory. This implies that the operator product expansion (OPE) coefficients involving this multiplet vanish whenever the central charge $c$ is identical to that of the AD theory. Similarly, by considering the n-th power of the current multiplet, we show that a particular short multiplet and OPE coefficients vanish for a class of AD theories with ADE flavor symmetry. We also consider the generalized AD theory of type $(A_{k-1}, A_{n-1})$ for coprime k, n and compute its Macdonald index using the associated W-algebra under a mild assumption. This allows us to show that a number of short multiplets and OPE coefficients vanish in this theory. We also provide a Mathematica file along with this paper, where we implement the algorithm by Cordova-Dumitrescu-Intriligator to compute the spectrum of 4d N=2 superconformal multiplets as well as their superconformal character.


title: Generative Adversarial Networks for Recovering Missing Spectral Information

abstract: Ultra-wideband (UWB) radar systems nowadays typical operate in the low frequency spectrum to achieve penetration capability. However, this spectrum is also shared by many others communication systems, which causes missing information in the frequency bands. To recover this missing spectral information, we propose a generative adversarial network, called SARGAN, that learns the relationship between original and missing band signals by observing these training pairs in a clever way. Initial results shows that this approach is promising in tackling this challenging missing band problem.


title: A new investigation of half-lives for the decay modes of $^{50}$V

abstract: A new search for the decay modes of the 4-fold forbidden non-unique decay of $^{50}$V has been performed at the Gran Sasso Underground Laboratory (LNGS). In total an exposure of 197 kg $\times$ d has been accumulated. The half-life for the electron capture into the first excited state of $^{50}$Ti has been measured with the highest precision to date as $2.67_{-0.18}^{+0.16} \times 10^{17}$ yr (68% C.I.) in which systematics uncertainties dominate. The search for the $\beta$-decay into the first excited state of $^{50}$Cr resulted in a lower limit of ${1.9} \times 10^{19}$ yr (90% C.I.), which is an improvement of almost one order of magnitude compared to existing results. The sensitivity of the new measurement is now in the region of theoretical predictions.


title: Transverse momentum $p_{\mathrm{T}}$ spectra of strange particles production in different collisions at $\sqrt{s_{NN}}= 2.76, 5.02,$ and $7$ TeV

abstract: We analyse the transverse momentum $p_{\mathrm{T}}$ spectra of strange particles \Kslxi produced in $\textsf{Pb+Pb}$ collision at $\sqrt{s_{NN}}= 2.76$ TeV, $\textsf{p+Pb}$ collision at $\sqrt{s_{NN}}= 5.02$ TeV, and $\textsf{p+p}$ collision at $\sqrt{s_{NN}}= 7$ TeV in different multiplicity events measured by the \textit{CMS} experiment at the Large Hadron Collider. The $p_{\mathrm{T}}$ spectra of strange particles are fitted by Tsallis statistics and Boltzmann statistics, respectively. The fitting parameters are studied as a function of the multiplicity events for all systems. The Tsallis temperature ($T_{\mathtt{Ts}}$), Boltzmann temperature ($T_{\mathtt{Boltz}}$), and radius of the system ($R$) increase with both the mass and strangeness number of the particle and also increase with the multiplicity events. The non-extensive parameter ($q$) decreases with the increase in the mass of the particle and also decrease with the increase in the multiplicity events which means that the system tends to thermodynamic stabilization. The extracted temperatures from the two statistics for the strange particles are exhibited a linear correlation.


title: The order of the non-abelian tensor product of groups

abstract: Let $G$ and $H$ be groups that act compatibly on each other. We denote by $[G,H]$ the derivative subgroup of $G$ under $H$. We prove that if the set $\{g^{-1}g^h \mid g \in G, h \in H\}$ has $m$ elements, then the derivative $[G,H]$ is finite with $m$-bounded order. Moreover, we show that if the set of all tensors $T_{\otimes}(G,H) = \{g\otimes h \mid g \in G, h\in H\}$ has $m$ elements, then the non-abelian tensor product $G \otimes H$ is finite with $m$-bounded order. We also examine some finiteness conditions for the non-abelian tensor square of groups.


title: An efficient supervised dictionary learning method for audio signal recognition

abstract: Machine hearing or listening represents an emerging area. Conventional approaches rely on the design of handcrafted features specialized to a specific audio task and that can hardly generalized to other audio fields. For example, Mel-Frequency Cepstral Coefficients (MFCCs) and its variants were successfully applied to computational auditory scene recognition while Chroma vectors are good at music chord recognition. Unfortunately, these predefined features may be of variable discrimination power while extended to other tasks or even within the same task due to different nature of clips. Motivated by this need of a principled framework across domain applications for machine listening, we propose a generic and data-driven representation learning approach. For this sake, a novel and efficient supervised dictionary learning method is presented. The method learns dissimilar dictionaries, one per each class, in order to extract heterogeneous information for classification. In other words, we are seeking to minimize the intra-class homogeneity and maximize class separability. This is made possible by promoting pairwise orthogonality between class specific dictionaries and controlling the sparsity structure of the audio clip's decomposition over these dictionaries. The resulting optimization problem is non-convex and solved using a proximal gradient descent method. Experiments are performed on both computational auditory scene (East Anglia and Rouen) and synthetic music chord recognition datasets. Obtained results show that our method is capable to reach state-of-the-art hand-crafted features for both applications.


title: Product-free sets in the free semigroup

abstract: In this paper, we study product-free subsets of the free semigroup over a finite alphabet $A$. We prove that the maximum density of a product-free subset of the free semigroup over $A$, with respect to the natural measure that assigns a weight of $|A|^{-n}$ to each word of length $n$, is precisely $1/2$.


title: An Exchange Mechanism to Coordinate Flexibility in Residential Energy Cooperatives

abstract: Energy cooperatives (ECs) such as residential and industrial microgrids have the potential to mitigate increasing fluctuations in renewable electricity generation, but only if their joint response is coordinated. However, the coordination and control of independently operated flexible resources (e.g., storage, demand response) imposes critical challenges arising from the heterogeneity of the resources, conflict of interests, and impact on the grid. Correspondingly, overcoming these challenges with a general and fair yet efficient exchange mechanism that coordinates these distributed resources will accommodate renewable fluctuations on a local level, thereby supporting the energy transition. In this paper, we introduce such an exchange mechanism. It incorporates a payment structure that encourages prosumers to participate in the exchange by increasing their utility above baseline alternatives. The allocation from the proposed mechanism increases the system efficiency (utilitarian social welfare) and distributes profits more fairly (measured by Nash social welfare) than individual flexibility activation. A case study analyzing the mechanism performance and resulting payments in numerical experiments over real demand and generation profiles of the Pecan Street dataset elucidates the efficacy to promote cooperation between co-located flexibilities in residential cooperatives through local exchange.


title: Considering Race a Problem of Transfer Learning

abstract: As biometric applications are fielded to serve large population groups, issues of performance differences between individual sub-groups are becoming increasingly important. In this paper we examine cases where we believe race is one such factor. We look in particular at two forms of problem; facial classification and image synthesis. We take the novel approach of considering race as a boundary for transfer learning in both the task (facial classification) and the domain (synthesis over distinct datasets). We demonstrate a series of techniques to improve transfer learning of facial classification; outperforming similar models trained in the target's own domain. We conduct a study to evaluate the performance drop of Generative Adversarial Networks trained to conduct image synthesis, in this process, we produce a new annotation for the Celeb-A dataset by race. These networks are trained solely on one race and tested on another - demonstrating the subsets of the CelebA to be distinct domains for this task.


title: The Color and Binarity of (486958) 2014 MU69 and Other Long-Range New Horizons Kuiper Belt Targets

abstract: The Hubble Space Telescope (HST) measured the colors of eight Kuiper Belt Objects (KBOs) that will be observed by the New Horizons spacecraft including its 2019 close fly-by target the Cold Classical KBO (486958) 2014 MU69. We find that the photometric colors of all eight objects are red, typical of the Cold Classical dynamical population within which most reside. Because 2014 MU69 has a similar color to that of other KBOs in the Cold Classical region of the Kuiper Belt, it may be possible to use the upcoming high-resolution New Horizons observations of 2014 MU69 to draw conclusions about the greater Cold Classical population. Additionally, HST found none of these KBOs to be binary within separations of ~0.06 arcsec (~2000 km at 44 AU range) and {\Delta}m less than or equal to 0.5. This conclusion is consistent with the lower fraction of binaries found at relatively wide separations. A few objects appear to have significant photometric variability, but our observations are not of sufficient signal-to-noise or time duration for further interpretation.


title: Constant-mean-curvature Slicing of the Swiss-cheese Universe

abstract: A sequence of Constant-Mean-Curvature(CMC) slices in the Swiss-Cheese(SC) Universe is investigated. We focus on the CMC slices which smoothly connect to the homogeneous time slices in the Einstein-de Sitter region in the SC universe. It is shown that the slices do not pass through the black hole region but white hole region.


title: Gradient Descent Happens in a Tiny Subspace

abstract: We show that in a variety of large-scale deep learning scenarios the gradient dynamically converges to a very small subspace after a short period of training. The subspace is spanned by a few top eigenvectors of the Hessian (equal to the number of classes in the dataset), and is mostly preserved over long periods of training. A simple argument then suggests that gradient descent may happen mostly in this subspace. We give an example of this effect in a solvable model of classification, and we comment on possible implications for optimization and learning.


title: Aerial navigation in obstructed environments with embedded nonlinear model predictive control

abstract: We propose a methodology for autonomous aerial navigation and obstacle avoidance of micro aerial vehicles (MAV) using nonlinear model predictive control (NMPC) and we demonstrate its effectiveness with laboratory experiments. The proposed methodology can accommodate obstacles of arbitrary, potentially non-convex, geometry. The NMPC problem is solved using PANOC: a fast numerical optimization method which is completely matrix-free, is not sensitive to ill conditioning, involves only simple algebraic operations and is suitable for embedded NMPC. A C89 implementation of PANOC solves the NMPC problem at a rate of 20Hz on board a lab-scale MAV. The MAV performs smooth maneuvers moving around an obstacle. For increased autonomy, we propose a simple method to compensate for the reduction of thrust over time, which comes from the depletion of the MAV's battery, by estimating the thrust constant.


title: Lattice-Boltzmann simulations of electrowetting phenomena

abstract: We present a lattice-Boltzmann method that can simulate the coupled hydrodynamics and electrostatics equations of motion of a two-phase fluid as a means to model electrowetting phenomena. Our method has the advantage of modelling the electrostatic fields within the lattice-Boltzmann algorithm itself, eliminating the need for a hybrid method. We validate our method by reproducing the static equilibrium configuration of a droplet subject to an applied voltage and show that the apparent contact angle of the drop depends on the voltage following the Young-Lippmann equation up to contact angles of $\approx 50^\circ$. At higher voltages, we observe a saturation of the contact angle caused by the competition between electric and capillary stresses. We also study the stability of a dielectric film trapped between a conducting fluid and a solid electrode and find a good agreement with analytical predictions based on lubrication theory. Finally, we investigate the film dynamics at long times and report observations of film breakup and entrapment similar to previously reported experimental results.


title: Quasiperiodicity and valence fluctuation in the spin-1/2 Falicov-Kimball model

abstract: We study the spin-1/2 Falicov-Kimball model with conduction and localized $f$ electrons on the Penrose lattice using the real-space dynamical mean-field theory. By changing the $f$ electron level, the $f$ electron density at each site changes continuously, in contrast to periodic systems with first-order valence transitions. In the intermediate valence regime, the local $f$ electron number strongly depends on a wider range of the Penrose structure surrounding its lattice site, in spite of the local interaction between the conduction and $f$ electrons. The temperature dependence of the magnetic response is also discussed.


title: The HST Lightcurve of (486958) 2014 MU69

abstract: We report HST lightcurve observations of the New Horizons (NH) spacecraft encounter KBO (486958) 2014 MU69 acquired near opposition in July 2017. In order to plan the optimum flyby sequence the NH mission planners needed to learn as much as possible about the target in advance of the encounter. Specifically, from lightcurve data, encounter timing could be adjusted to accommodate a highly elongated, binary, or rapidly rotating target. HST astrometric (Porter et al. 2018) and stellar occultation (Buie et al. 2018) observations constrained MU69's orbit and diameter (21-41 km for an albedo of 0.15-0.04), respectively. Photometry from the astrometric dataset suggested a variability of $\ge$0.3 mags, but they did not determine the period or provide shape information. We strategically spaced 24 HST orbits over 9 days to investigate rotation periods from approximately 2-100 hours and to better constrain the lightcurve amplitude. Until NH detected MU69 in its optical navigation images beginning in August 2018, this HST campaign provided the most accurate photometry to date. The mean variation in our data is 0.15 mags which suggests that MU69 is either nearly spherical (a:b axis ratio of 1:1.15), or its pole vector is pointed near the line of sight to Earth; this interpretation does not preclude a near-contact binary or bi-lobed object. However, image stacks do conclude that MU69 does not have a binary companion $\ge$2000 km with a sensitivity to 29th magnitude. We report with confidence that MU69 is not both rapidly rotating and highly elongated. We note that our results are consistent with the fly-by imagery and orientation of MU69 (Stern et al. 2019). The combined dataset also suggests that within the KBO lightcurve literature there are likely other objects which share a geometric configuration like MU69 resulting in an underestimate of the contact binary fraction for the CC Kuiper Belt.


title: Chaotic few-body vortex dynamics in rotating Bose--Einstein condensates

abstract: We investigate a small vortex-lattice system of four co-rotating vortices in an atomic Bose--Einstein condensate and find that the vortex dynamics display chaotic behaviour after a system quench introduced by reversing the direction of circulation of a single vortex through a phase-imprinting process. By tracking the vortex trajectories and Lyapunov exponent, we show the onset of chaotic dynamics is not immediate, but occurs at later times and is accelerated by the close-approach and separation of all vortices in a scattering event. The techniques we develop could potentially be applied to create locally induced chaotic dynamics in larger lattice systems as a stepping stone to study the role of chaotic events in turbulent vortex dynamics.


title: $\ell^2$ Decoupling in $\mathbb R^2$ for curves with Vanishing Curvature

abstract: We expand the class of curves $(\varphi_1(t),\varphi_2(t)),\ t\in[0,1]$ for which the $\ell^2$ decoupling conjecture holds for $2\leq p\leq 6$. Our class of curves includes all real-analytic regular curves with isolated points of vanishing curvature and all curves of the form $(t,t^{1+\nu})$ for $\nu\in (0,\infty)$.


title: A rigidity theorem for ideal surfaces with flat boundary

abstract: We consider surfaces with boundary satisfying a sixth order nonlinear elliptic partial differential equation corresponding to extremising the $L^2$-norm of the gradient of the mean curvature. We show that such surfaces with small $L^2$-norm of the second fundamental form and satisfying so-called `flat boundary conditions' are necessarily planar.


title: Regularization Properties of the Krylov Iterative Solvers CGME and LSMR For Linear Discrete Ill-Posed Problems with an Application to Truncated Randomized SVDs

abstract: For the large-scale linear discrete ill-posed problem $\min\|Ax-b\|$ or $Ax=b$ with $b$ contaminated by Gaussian white noise, there are four commonly used Krylov solvers: LSQR and its mathematically equivalent CGLS, the Conjugate Gradient (CG) method applied to $A^TAx=A^Tb$, CGME, the CG method applied to $\min\|AA^Ty-b\|$ or $AA^Ty=b$ with $x=A^Ty$, and LSMR, the minimal residual (MINRES) method applied to $A^TAx=A^Tb$. These methods have intrinsic regularizing effects, where the number $k$ of iterations plays the role of the regularization parameter. In this paper, we establish a number of regularization properties of CGME and LSMR, including the filtered SVD expansion of CGME iterates, and prove that the 2-norm filtering best regularized solutions by CGME and LSMR are less accurate than and at least as accurate as those by LSQR, respectively. We also prove that the semi-convergence of CGME and LSMR always occurs no later and sooner than that of LSQR, respectively. As a byproduct, using the analysis approach for CGME, we improve a fundamental result on the accuracy of the truncated rank $k$ approximate SVD of $A$ generated by randomized algorithms, and reveal how the truncation step damages the accuracy. Numerical experiments justify our results on CGME and LSMR.


title: Granular Graphene: direct observation of zigzag and armchair edge waves

abstract: We propose a mechanical granular graphene obtained by replacing the carbon atoms with macroscopic spherical stainless steel beads in contact. The experimental measured dispersion relation is presented, in conjunction with evidence of the Dirac points. In addition, wave propagation along the zigzag and a robust turning effect of edge waves from the zigzag to the armchair boundary is experimentally revealed, even in the absence of a full band gap for bulk modes. Our work shows that mechanical granular graphene can serve as an excellent experimental platform to study novel Dirac, topological and nonlinear wave phenomena.


title: DEAP-3600 Recent Dark Matter Results

abstract: The DEAP-3600 experiment is searching for WIMP dark matter with a 3.3 tonne single phase liquid argon (LAr) target, located at SNOLAB. The construction and filling of DEAP-3600 was completed in 2016, and the experiment is currently taking physics data. First results were recently published, which demonstrated stable detector operations and the power of pulse shape discrimination to reject electron-recoil backgrounds in LAr. In addition, the most sensitive WIMP exclusion with a LAr target was achieved at the time of publication.


title: Strong lensing constraints on modified gravity models

abstract: We impose the first strong-lensing constraints on a wide class of modified gravity models where an extra field that modifies gravity also couples to photons (either directly or indirectly through a coupling with baryons) and thus modifies lensing. We use the nonsingular isothermal ellipsoid (NIE) profile as an effective potential, which produces flat galactic rotation curves. If a concrete modified gravity model gives a flat rotation curve, then the parameter $\Gamma$ that characterizes the lensing effect must take some definite value. We find that $\Gamma = 1.24\pm0.65$ at $1\sigma$, consistent with general relativity ($\Gamma = 1$). This constrains the parameter space in some recently proposed models.